Cancer is one of the leading causes of death worldwide, bringing a significant burden to human health and society. Accurate cancer diagnosis and biomarkers that can be used as robust therapeutic targets are of great importance as they facilitate early and effective therapies. Shared etiology among cancers suggests the existence of pan-cancer biomarkers, performance of which could benefit from the large sample size and the heterogeneity of the studied patients. In this study, we conducted a systematic RNA-seq study of 9,213 tumors and 723 para-cancerous tissue samples of 28 solid tumors from the Cancer Genome Atlas (TCGA) database, and 7,008 normal tissue samples from the Genotype-Tissue Expression (GTEx) database. By differential gene expression analysis, we identified 214 up-regulated and 186 downregulated differentially expressed genes (DEGs) in more than 80% of the studied tumors, respectively, and obtained 20 highly linked up- and downregulated hub genes from them. These markers have rarely been reported in multiple tumors simultaneously. We further constructed pan-cancer diagnostic models to classify tumors and para-cancerous tissues using 10 up-regulated hub genes with an AUC of 0.894. Survival analysis revealed that these hub genes were significantly associated with the overall survival of cancer patients. In addition, drug sensitivity predictions for these hub genes in a variety of tumors obtained several broad-spectrum anti-cancer drugs targeting pan-cancer. Furthermore, we predicted immunotherapy sensitivity for cancers based on tumor mutational burden (TMB) and the expression of immune checkpoint genes (ICGs), providing a theoretical basis for the treatment of tumors. In summary, we identified a set of biomarkers that were differentially expressed in multiple types of cancers, and these biomarkers can be potentially used for diagnosis and used as therapeutic targets.