Email Alert | RSS

Journal of Peking University (Health Sciences) ›› 2021, Vol. 53 ›› Issue (3): 602-607. doi: 10.19723/j.issn.1671-167X.2021.03.028

Previous Articles     Next Articles

Exploratory screening of potential pan-cancer biomarkers based on The Cancer Genome Atlas database

ZHOU Chuan,MA Xue,XING Yun-kun,LI Lu-di,CHEN Jie,YAO Bi-yun,FU Juan-ling,ZHAO PengΔ()   

  1. Department of Toxicology, Beijing Key Laboratory of Toxicological Research and Risk Assessment for Food Safety, Peking University School of Public Health, Beijing 100191, China
  • Received:2020-11-02 Online:2021-06-18 Published:2021-06-16
  • Contact: Peng ZHAO E-mail:zhaopeng@bjmu.edu.cn
  • Supported by:
    National Natural Science Foundation of China(81370079);National Natural Science Foundation of China(81001253);Beijing Natural Science Foundation(7132122)

RICH HTML

   

PDF (PC)

Abstract:

Objective: To screen potential pan-cancer biomarkers based on The Cancer Genome Atlas (TCGA) database, and to provide help for the diagnosis and prognosis assessment of a variety of cancers. Methods: “GDC Data Transfer Tool” and “GDCRNATools” packages were used to obtain TCGA database. After data sorting, a total of 13 cancers were selected for further analysis. False disco-very rate (FDR) <0.05 and fold change (FC) >1.5 were used as the differential expression criteria to screen genes and miRNAs that were up- or down-regulated in all the 13 cancers. In the receiver operating characteristic curve (ROC curve), the area under the curve (AUC), the best cut-off value and the corresponding sensitivity and specificity were used to reflect diagnostic significance. The Kaplan-Meier method was used to calculate the survival probability and then the log-rank test was performed. Hazard ratio (HR) was calculated to reflect prognostic evaluation significance. DAVID tool were used to perform GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis for differentially expressed genes. STRING and TargetScan tools were used to analyze the regulatory network of differentially expressed genes and miRNAs. Results: A total of 48 genes and 2 miRNAs were differentially expressed in all the 13 cancers. Among them, 25 genes were up-regulated, 23 genes and 2 miRNAs were down-regulated. Most differentially expressed genes and miRNAs had good ability to distinguish between the cases and controls, with AUC, sensitivity and specificity up to 0.8-0.9. Survival analysis results show that differentially expressed genes and miRNAs were significantly associated with patient survival in a variety of cancers. Most up-regulated genes were risk factors for patient survival (HR>1), while most down-regulated genes were protective factors for patient survival (0<HR<1). The enrichment analysis of GO and KEGG showed that the differentially expressed genes were mostly enriched in biological events related to cell proliferation. In the regulatory network analysis, a total of 13 differentially expressed genes and 2 differentially expressed miRNAs had regulatory and interaction relationships. Conclusion: The 48 genes and 2 miRNAs that were differentially expressed in 13 cancers may serve as potential pan-cancer biomarkers, providing help for the diagnosis and prognosis evaluation of a variety of cancers, and providing clues for the development of broad-spectrum tumor therapeutic targets.

Key words: Pan-cancer, Biomarkers,tumor, Gene expression regulation, Genome,human

CLC Number: 

  • R730.43

Table 1

Information of the TCGA projects included in the study"

Project Disease Gene miRNAs
Case Control Case Control
TCGA-BLCA Bladder urothelial carcinoma 408 19 409 19
TCGA-BRCA Breast invasive carcinoma 1 091 113 1 078 104
TCGA-HNSC Head and neck squamous cell carcinoma 500 44 523 44
TCGA-KICH Kidney chromophobe 65 24 66 25
TCGA-KIRC Kidney renal clear cell carcinoma 530 72 516 71
TCGA-KIRP Kidney renal papillary cell carcinoma 288 32 291 34
TCGA-LIHC Liver hepatocellular carcinoma 371 50 372 50
TCGA-LUAD Lung adenocarcinoma 513 59 513 46
TCGA-LUSC Lung squamous cell carcinoma 501 49 478 45
TCGA-PRAD Prostate adenocarcinoma 495 52 494 52
TCGA-STAD Stomach adenocarcinoma 375 32 436 41
TCGA-THCA Thyroid carcinoma 502 58 506 59
TCGA-UCEC Uterine corpus endometrial carcinoma 543 35 538 33

Figure 1

Fold change of differentially expressed genes and microRNAs BLCA, bladder urothelial carcinoma; BRCA, breast invasive carcinoma; HNSC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; PRAD, prostate adenocarcinoma; STAD, stomach adenocarcinoma; THCA, thyroid carcinoma; UCEC, uterine corpus endometrial carcinoma; FC, fold change; miR, microRNA."

Figure 2

Diagnostic significance of differentially expressed genes and microRNAs AUC, area under the curve; Other abbreviations as in Figure 1."

Figure 3

Relationship with patient survival of differentially expressed genes and microRNAs HR, hazard ratio; Other abbreviations as in Figure 1."

Figure 4

GO and KEGG pathway enrichment analysis of up-regulated differentially expressed genes (▲ number of genes) GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; BP, biological process; CC, cellular component; MF, molecular function."

Figure 5

GO analysis of down-regulated differentially expressed genes (▲ number of genes) Abbreviations as in Figure 4."

Figure 6

Regulatory network of differentially expressed genes and microRNAs"

[1] Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA Cancer J Clin, 2018,68(6):394-424.
doi: 10.3322/caac.v68.6
[2] Vargas AJ, Harris CC. Biomarker development in the precision medicine era: lung cancer as a case study[J]. Nat Rev Cancer, 2016,16(8):525-537.
doi: 10.1038/nrc.2016.56
[3] 王印祥. 泛肿瘤研究和肿瘤免疫研究: 未来抗肿瘤药的发展趋势[J]. 中国药物化学杂志, 2015,25(2):149-152.
[4] Zack TI, Schumacher SE, Carter SL, et al. Pan-cancer patterns of somatic copy number alteration[J]. Nat Genet, 2013,45(10):1134-1140.
doi: 10.1038/ng.2760
[5] 陈熹, 张峻峰, 曾科, 等. 血清microRNA: 一种非侵入性的肿瘤标志物[J]. 生命科学, 2010,22(7):649-654.
[6] Li RD, Qu H, Wang SB, et al. GDCRNATools: an R/Bioconductor package for integrative analysis of lncRNA, miRNA and mRNA data in GDC[J]. Bioinformatics, 2018,34(14):2515-2517.
doi: 10.1093/bioinformatics/bty124
[7] Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources[J]. Nat Protoc, 2009,4(1):44-57.
doi: 10.1038/nprot.2008.211 pmid: 19131956
[8] Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets[J]. Nucleic Acids Res, 2019,47(D1):D607-D613.
doi: 10.1093/nar/gky1131
[9] Agarwal V, Bell GW, Nam JW, et al. Predicting effective microRNA target sites in mammalian mRNAs[J]. Elife, 2015,4:e05005.
doi: 10.7554/eLife.05005
[10] Cancer Genome Atlas Research Network, Weinstein JN, Collisson EA, et al. The Cancer Genome Atlas Pan-Cancer analysis project[J]. Nat Genet, 2013,45(10):1113-1120.
doi: 10.1038/ng.2764
[11] Penault-Llorca F, Radosevic-Robin N. Ki67 assessment in breast cancer: an update[J]. Pathology, 2017,49(2):166-171.
doi: S0031-3025(16)40573-8 pmid: 28065411
[12] Yuniati L, Scheijen B, van der Meer LT, et al. Tumor suppressors BTG1 and BTG2: Beyond growth control[J]. J Cell Physiol, 2019,234(5):5379-5389.
doi: 10.1002/jcp.27407
[13] Bhattacharjee S, Rajaraman P, Jacobs KB, et al. A subset-based approach improves power and interpretation for the combined analysis of genetic association studies of heterogeneous traits[J]. Am J Hum Genet, 2012,90(5):821-835.
doi: 10.1016/j.ajhg.2012.03.015
[14] Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation[J]. Cell, 2011,144(5):646-674.
doi: 10.1016/j.cell.2011.02.013
[15] Laine A, Westermarck J. Molecular pathways: harnessing E2F1 regulation for prosenescence therapy in p53-defective cancer cells[J]. Clin Cancer Res, 2014,20(14):3644-3650.
doi: 10.1158/1078-0432.CCR-13-1942 pmid: 24788101
[16] Xiao YS, Najeeb RM, Ma D, et al. Upregulation of CENPM promotes hepatocarcinogenesis through mutiple mechanisms[J]. J Exp Clin Cancer Res, 2019,38(1):458.
doi: 10.1186/s13046-019-1444-0
[17] Sun JB, Huang JZ, Lan J, et al. Overexpression of CENPF correlates with poor prognosis and tumor bone metastasis in breast cancer[J]. Cancer Cell Int, 2019,19(1):264.
doi: 10.1186/s12935-019-0986-8
[1] Yanhua LIU, Min LU, Xuyang ZHAO, Kuan'gen ZHANG, Rui WU, Fang MEI, Zhihao DAI, Jiangfeng YOU, Fei PEI. Effect of dephosphorylation of tumor metastasis suppressor gene LASS2 on vacuolar ATPase activity and invasiveness of prostate cancer [J]. Journal of Peking University (Health Sciences), 2025, 57(6): 1113-1123.
[2] GUO Qian, CHEN Xu-yong, SU Yin. Interleukin-2 signaling pathway regulating molecules in systemic lupus erythematosus [J]. Journal of Peking University(Health Sciences), 2016, 48(6): 1100-1104.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!
Copyright © Journal of Peking University (Health Sciences), All Rights Reserved.
Powered by Beijing Magtech Co. Ltd