An eight-long noncoding RNA expression signature for colorectal cancer patients’ prognosis

Long noncoding RNAs (lncRNAs) have recently emerged as important biomarkers of cancer progression. Here, we proposed to develop a lncRNA-based signature with a prognostic value for colorectal cancer (CRC) overall survival (OS). Through mining microarray datasets, we analyzed the lncRNA expression profiles of 122 patients with CRC from Gene Expression Omnibus. Associations between lncRNA and CRC OS were firstly evaluated through univariate Cox regression analysis. A random survival forest method was applied for further screening of the lncRNA signature, which resulted in eight lncRNAs, including PEG3-AS1, LOC100505715, MINCR, DBH-AS1, LINC00664, FAM224A, LOC642852, and LINC00662. Combination of the eight lncRNAs weighted by their multivariate Cox regression coefficients formed a prognostic signature, through which, we could divide the 122 patients with CRC into two subgroups with significantly different OS. Good robustness of the lncRNA signature's prognostic value was verified through an independent data set consisting of 55 patients with CRC. In addition, gene set enrichment analysis indicated the potential association between high prognostic value and oxygen metabolism-related processes. This result should indicate that lncRNAs could be a useful signature for CRC prognosis.     

A novel multidimensional signature predicts prognosis in hepatocellular carcinoma patients

The abnormal expression of microRNAs (miRNAs) or protein-coding genes (PCGs) have been found to be associated with the prognosis of hepatocellular carcinoma (HCC) patients. Using bioinformatics analysis methods including Cox’s proportional hazards regression analysis, the random survival forest algorithm, Kaplan–Meier, and receiver operating characteristic (ROC) curve analysis, we mined the gene expression profiles of 469 HCC patients from The Cancer Genome Atlas (n = 379) and Gene Expression Omnibus (GSE14520; n = 90) public database. We selected a signature comprising one protein-coding gene (PCG; DNA polymerase μ) and three miRNAs (hsa-miR-149-5p, hsa-miR-424-5p, hsa-miR-579-5p) with highest accurate prediction (area under the ROC curve [AUC] = 0.72; n = 189) from the training data set. The signature stratified patients into high- and low-risk groups with significantly different survival (median 27.9 vs. 55.2 months, log-rank test, p < 0.001) in the training data set, and its risk stratification ability were validated in the test data set (median 47.4 vs. 84.4 months, log-rank test, p = 0.03) and an independent data set (median 31.0 vs. 46.0 months, log-rank test, p = 0.01). Multivariable Cox regression analysis showed that the signature was an independent prognostic factor. And the signature was proved to have a better survival prediction power than tumor–node–metastasis (TNM) stage (AUC _signature = 0.72/0.64/0.62 vs. AUC _TNM = 0.65/0.61/0.61; p < 0.05). Moreover, we validated the expression of these prognostic genes from the PCG-miRNA signature in Huh-7 cell by real-time polymerase chain reaction. In conclusion, we found a signature that can predict survival of HCC patients and serve as a prognostic marker for HCC.     

A novel seven-lncRNA signature for prognosis prediction in hepatocellular carcinoma

Liver cancer is still one of the leading causes of cancer-related death worldwide. This study is dedicated to developing a multi–long noncoding RNA (lncRNA) model for risk stratification and prognosis prediction on patients with hepatocellular carcinoma (HCC). We first downloaded lncRNA expression profiles and corresponding clinical information of patients with liver cancer from The Cancer Genome Atlas database. Differentially expressed (DE) lncRNAs between HCC samples and normal samples were identified. In total, 308 patients with HCC were randomly divided into a training group (n = 154) and a testing group (n = 154). Univariate Cox regression and least absolute shrinkage and selection operator Cox regression analyses were performed to select the best survival-related candidates from these DE lncRNAs in the training set. Seven lncRNAs (AC009005.2, RP11-363N22.3, RP11-932O9.10, RP11-572O6.1, RP11-190C22.8, RP11-388C12.8, and ZFPM2-AS1) were finally identified and used to construct a seven-lncRNA signature. The signature could classify patients into high-risk and low-risk groups with significantly different overall survival. The area under the curve of receiver operating characteristic curve for the signature to predict 5-year survival reached more than 0.75. Besides, the prognostic value of the seven-lncRNA signature was independent of conventional clinical factors. The predictive performance of the signature was further validated in the testing set and the whole set. Functional enrichment analysis indicated that the seven prognostic lncRNAs may be involved in several essential biological processes and pathways. The current study demonstrated the potential clinical implications of the seven-lncRNA signature for survival prediction of patients with HCC.     

A prognostic eight-lncRNA expression signature in predicting recurrence of ER-positive breast cancer receiving endocrine therapy

Long noncoding RNAs (lncRNAs) have the main role in the tumorigenesis of breast cancer. In the present study, lncRNA expression profiling was collected to identify a lncRNA expression signature from the Gene Expression Omnibus database. An eight-lncRNA signature was established to predict the survival of patients with estrogen receptor (ER)-positive breast cancer receiving endocrine therapy. Patients were separated into a low-risk group and a high-risk group based on this signature. Patients in high-risk group have worse survival compared to those in low-risk group using Kaplan–Meier curve analysis with log-rank test. Receiver operating characteristic analysis suggested good diagnostic efficiency of the eight-lncRNA signature. When adjusting the clinical features, including age, grade, lymph node status, and tumor size, this signature was independently associated with the relapse-free survival. The prognostic value of the lncRNA prognostic model was then validated in validation sets. When validated in a cohort of patients treated with neoadjuvant chemotherapy and endocrine therapy, this signature demonstrated good performance as well. Besides, we have built a nomogram that integrated the conventional clinicopathological features and the eight-lncRNA-based signature. To sum up, our results indicated that the eight-lncRNA prognostic model was a reliable tool to group patients at high and low risk of disease relapse. This signature may have possible implication in prognostic evaluations of patients with ER-positive breast cancer receiving endocrine therapy.     

A 16-mRNA signature optimizes recurrence-free survival prediction of Stages II and III gastric cancer

High-throughput messenger RNA (mRNA) analysis has become a powerful tool for exploring tumor recurrence or metastasis mechanisms. Here, we constructed a signature to predict the recurrence risk of Stages II and III gastric cancer (GC) patients. A least absolute shrinkage and selection operator method Cox regression model was utilized to construct the signature. Using this method, a 16-mRNA signature was identified to be associated with the relapse-free survival of Stages II and III GCs in training dataset GSE62254 (n = 194). Then this signature was validated in an independent Gene Expression Omnibus cohort GSE26253 (n = 297) and a dataset of The Cancer Genome Atlas (TCGA; n = 235). This classifier could successfully screen out the high-risk Stages II and III GCs in the training cohort (hazard ratio [HR] = 40.91; 95% confidence interval [CI] = 5.58–299.7; p < .0001). Analysis in two independent validation cohorts yielded consistent results (GSE26253: HR = 1.69, 95% CI = 1.17–2.43,; p = .0045; TCGA: HR = 2.01, 95% CI = 1.13–3.56, p = .0146). Cox regression analyses revealed that the risk score derived from this signature was an independent risk factor in Stages II and III GCs. Besides, a nomogram was constructed to serve clinical practice. Through gene set variation analysis, we found several gene sets associated with chemotherapeutic drug resistance and tumor metastasis significantly enriched in high-risk patients. In summary, this 16-mRNA signature can be used as a powerful tool for prognostic evaluation and help clinicians identify high-risk patients.     

Identification of 4-lncRNA prognostic signature in head and neck squamous cell carcinoma

Deregulated long noncoding RNAs (lncRNA) have been critically implicated in tumorigenesis and serve as novel diagnostic and prognostic biomarkers. Here we sought to develop a prognostic lncRNA signature in patients with head and neck squamous cell carcinoma (HNSCC). Original RNA-seq data of 499 HNSCC samples were retrieved from The Cancer Genome Atlas database, which was randomly divided into training and testing set. Univariate Cox regression survival analysis, robust likelihood-based survival model and random sampling iterations were applied to identify prognostic lncRNA candidates in the training cohort. A prognostic risk score was developed based on the Cox coefficient of four individual lncRNA imputed as follows: (0.14546 × expression level of RP11-366H4.1) + (0.27106 × expression level of LINC01123) + (0.54316 × expression level of RP11-110I1.14) + (−0.48794 × expression level of CTD-2506J14.1). Kaplan-Meier analysis revealed that patients with high-risk score had significantly reduced overall survival as compared with those with low-risk score when patients in training, testing, and validation cohorts were stratified into high- or low-risk subgroups. Multivariate survival analysis further revealed that this 4-lncRNA signature was a novel and important prognostic factor independent of multiple clinicopathological parameters. Importantly, ROC analyses indicated that predictive accuracy and sensitivity of this 4-lncRNA signature outperformed those previously well-established prognostic factors. Noticeably, prognostic score based on quantification of these 4-lncRNA via qRT-PCR in another independent HNSCC cohort robustly stratified patients into subgroups with high or low survival. Taken together, we developed a robust 4-lncRNA prognostic signature for HNSCC that might provide a novel powerful prognostic biomarker for precision oncology.     

Development and validation of a five-lncRNA signature with prognostic value in colon cancer

Dysregulation of long noncoding RNAs (lncRNAs) has been found in a large number of human cancers, including colon cancer. Therefore, the implementation of potential lncRNAs biomarkers with prognostic prediction value are very much essential. GSE39582 data set was downloaded from database of Gene Expression Omnibus. Re-annotation analysis of lncRNA expression profiles was performed by NetAffx annotation files. Univariate and multivariate Cox proportional analyses helped select prognostic lncRNAs. Algorithm of random survival forest-variable hunting (RSF-VH) together with stepwise multivariate Cox proportional analysis were performed to establish lncRNA signature. The log-rank test was carried out to analyze and compare the Kaplan-Meier survival curves of patients’ overall survival (OS). Receiver operating characteristic (ROC) analysis was used for comparing the survival prediction regarding its specificity and sensitivity based on lncRNA risk score, followed by calculating the values of area under the curve (AUC). The single-sample GSEA (ssGSEA) analysis was used to describe biological functions associated with this signature. Finally, to determine the robustness of this model, we used the validation sets including GSE17536 and The Cancer Genome Atlas data set. After re-annotation analysis of lncRNAs, a total of 14 lncRNA probes were obtained by univariate and multivariate Cox proportional analysis. Then, the RSF-VH algorithm and stepwise multivariate Cox analysis helped to build a five-lncRNA prognostic signature for colon cancer. The patients in group with high risk showed an obviously shorter survival time compared with patients in group with low risk with AUC of 0.75. In addition, the five-lncRNA signature can be used to independently predict the survival of patients with colon cancer. The ssGSEA analysis revealed that pathways such as extracellular matrix-receptor interaction was activated with an increase in risk score. These findings determined the strong power of prognostic prediction value of this five-lncRNA signature for colon cancer.     

Prognostic significance of lncRNA DANCR expression in human cancers: a systematic review and meta-analysis

Several studies demonstrated that lncRNA differentiation antagonizing non-protein coding RNA (lncRNA DANCR) expression might have the potential capacity to predict the cancer prognosis; however, definite conclusion has not been obtained. The aim of this meta-analysis was to evaluate the prognostic value of lncRNA DANCR expression in cancers. PubMed, Web of Science, Scopus, and Embase were comprehensively searched for relevant studies. Studies meeting all inclusion standards were included into this meta-analysis. The analysis of overall survival (OS), disease-free survival (DFS), or clinicopathological features was conducted. Total 11 studies containing 1154 cancer patients were analyzed in this meta-analysis. The results showed, compared with low lncRNA DANCR expression, high lncRNA DANCR expression was significantly associated with shorter OS (hazard ratio [HR] = 1.85; 95% CI = 1.52–2.26; P<0.01) and DFS (HR = 1.82; 95% CI = 1.43–2.32; P<0.01) in cancers. Besides, high lncRNA DANCR expression predicted deeper tumor invasion (P<0.01), earlier lymph node metastasis (P<0.01), earlier distant metastasis (P<0.01), and more advanced clinical stage (P<0.01) compared with low lncRNA DANCR expression in cancer populations. High lncRNA DANCR expression was associated with worse prognosis compared with low lncRNA DANCR expression in cancers. LncRNA DANCR expression could serve as a prognostic factor of human cancers.     

Clinicopathological and prognostic value of lncRNA PANDAR expression in solid tumors: Evidence from a systematic review and meta-analysis

PANDAR (promoter of CDKN1A antisense DNA damage activated RNA) has been shown to be aberrantly expressed in many types of cancer. Considering conflicting data, the current study was aimed to assess its potential role as a prognostic marker in malignant tumors. A comprehensive literature search of PubMed, Medline, and Web of Science was performed to identify all eligible studies describing the use of PANDAR as a prognostic factor for different types of cancer. Data related to overall survival (OS) and clinicopathologic features were collected and analyzed. The pooled hazard ratio (HR) and odds radio (OR) with a 95% confidence interval (CI) were used to estimate associations. Ten original studies containing 1,231 patients were included. The results showed that in patients with cancer, high PANDAR expression is correlated with lymph node metastasis (LNM; OR = 2.57; 95% CI, 1.76–3.81; p < 0.001), tumor stage (OR = 2.90; 95% CI, 1.25–6.75; p = 0.013), and tumor size (OR = 1.79; 95% CI, 1.11–2.91; p = 0.018). However, sensitivity analysis further demonstrated a significant association between high PANDAR expression and OS, both in multivariate and univariate analysis models (pooled HR 2.01; 95% CI, 1.17–3.44 and pooled HR 2.62; 95% CI, 1.98–3.47, respectively), after omitting one study. These results suggested that PANDAR expression might be indicative of advanced disease and poor prognosis in patients with cancer. Further studies are necessary to determine the value of this risk stratification biomarker in clinical management of patients with cancer.     

A prognostic 11 long noncoding RNA expression signature for breast invasive carcinoma

Breast cancer, the most common cancer in women worldwide, is associated with high mortality. The long non-coding RNAs (lncRNAs) with a little capacity of coding proteins is playing an increasingly important role in the cancer paradigm. Accumulating evidences demonstrate that lncRNAs have crucial connections with breast cancer prognosis while the studies of lncRNAs in breast cancer are still in its primary stage. In this study, we collected 1052 clinical patient samples, a comparatively large sample size, including 13 159 lncRNA expression profiles of breast invasive carcinoma (BRCA) from The Cancer Genome Atlas database to identify prognosis-related lncRNAs. We randomly separated all of these clinical patient samples into training and testing sets. In the training set, we performed univariable Cox regression analysis for primary screening and played the model for Robust likelihood-based survival for 1000 times. Then 11 lncRNAs with a frequency more than 600 were selected for prediction of the prognosis of BRCA. Using the analysis of multivariate Cox regression, we established a signature risk-score formula for 11 lncRNA to identify the relationship between lncRNA signatures and overall survival. The 11 lncRNA signature was validated both in the testing and the complete set and could effectively classify the high-/low-risk group with different OS. We also verified our results in different stages. Moreover, we analyzed the connection between the 11 lncRNAs and the genes of ESR1, PGR, and Her2, of which protein products (ESR, PGR, and HER2) were used to classify the breast cancer subtypes widely. The results indicated correlations between 11 lncRNAs and the gene of PGR and ESR1. Thus, a prognostic model for 11 lncRNA expression was developed to classify the BRAC clinical patient samples, providing new avenues in understanding the potential therapeutic methods of breast cancer.     

Identification of autophagy-related long non-coding RNA prognostic signature for breast cancer

Autophagy-related long non-coding RNAs (lncRNAs) disorders are related to the occurrence and development of breast cancer. The purpose of this study is to explore whether autophagy-related lncRNA can predict the prognosis of breast cancer patients. The autophagy-related lncRNAs prognostic signature was constructed by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression. We identified five autophagy-related lncRNAs (MAPT-AS1, LINC01871, AL122010.1, AC090912.1, AC061992.1) associated with prognostic value, and they were used to construct an autophagy-related lncRNA prognostic signature (ALPS) model. ALPS model offered an independent prognostic value (HR = 1.664, 1.381-2.006), where this risk score of the model was significantly related to the TNM stage, ER, PR and HER2 status in breast cancer patients. Nomogram could be utilized to predict survival for patients with breast cancer. Principal component analysis and Sankey Diagram results indicated that the distribution of five lncRNAs from the ALPS model tends to be low-risk. Gene set enrichment analysis showed that the high-risk group was enriched in autophagy and cancer-related pathways, and the low-risk group was enriched in regulatory immune-related pathways. These results indicated that the ALPS model composed of five autophagy-related lncRNAs could predict the prognosis of breast cancer patients.     

RAB25 expression is epigenetically downregulated in oral and oropharyngeal squamous cell carcinoma with lymph node metastasis

Oral and oropharyngeal squamous cell carcinoma (OOSCC) have a low survival rate, mainly due to metastasis to the regional lymph nodes. For optimal treatment of these metastases, a neck dissection is required; however, inaccurate detection methods results in under- and over-treatment. New DNA prognostic methylation biomarkers might improve lymph node metastases detection. To identify epigenetically regulated genes associated with lymph node metastases, genome-wide methylation analysis was performed on 6 OOSCC with (pN+) and 6 OOSCC without (pN0) lymph node metastases and combined with a gene expression signature predictive for pN+ status in OOSCC. Selected genes were validated using an independent OOSCC cohort by immunohistochemistry and pyrosequencing, and on data retrieved from The Cancer Genome Atlas. A two-step statistical selection of differentially methylated sequences revealed 14 genes with increased methylation status and mRNA downregulation in pN+ OOSCC. RAB25, a known tumor suppressor gene, was the highest-ranking gene in the discovery set. In the validation sets, both RAB25 mRNA (P = 0.015) and protein levels (P = 0.012) were lower in pN+ OOSCC. RAB25 mRNA levels were negatively correlated with RAB25 methylation levels (P < 0.001) but RAB25 protein expression was not. Our data revealed that promoter methylation is a mechanism resulting in downregulation of RAB25 expression in pN+ OOSCC and decreased expression is associated with lymph node metastasis. Detection of RAB25 methylation might contribute to lymph node metastasis diagnosis and serve as a potential new therapeutic target in OOSCC.     

Immune-related long noncoding RNA signature for predicting survival and immune checkpoint blockade in hepatocellular carcinoma

Long noncoding RNAs (lncRNAs) show multiple functions, including immune response. Recently, the immune-related lncRNAs have been reported in some cancers. We first investigated the immune-related lncRNA signature as a potential target in hepatocellular carcinoma (HCC) survival. The training set (n = 368) and the independent external validation cohort (n = 115) were used. Immune genes and lncRNAs coexpression were constructed to identify immune-related lncRNAs. Cox regression analyses were perfumed to establish the immune-related lncRNA signature. Regulatory roles of this signature on cancer pathways and the immunologic features were investigated. The correlation between immune checkpoint inhibitors and this signature was examined. In this study, the immune-related lncRNA signature was identified in HCC, which could stratify patients into high- and low-risk groups. This immune-related lncRNA signature was correlated with disease progression and worse survival and was an independent prognostic biomarker. Our immune-related lncRNA signature was still a powerful tool in predicting survival in each stratum of age, gender, and tumor stage. This signature mediated cell cycle, glycolysis, DNA repair, mammalian target of rapamycin signaling, and immunologic characteristics (i.e., natural killer cells vs. Th1 cells down, etc). This signature was associated with immune cell infiltration (i.e., macrophages M0, Tregs, CD4 memory T cells, and macrophages M1, etc.,) and immune checkpoint blockade (ICB) immunotherapy-related molecules (i.e., PD-L1, PD-L2, and IDO1). Our findings suggested that the immune-related lncRNA signature had an important value for survival prediction and may have the potential to measure the response to ICB immunotherapy. This signature may guide the selection of the immunotherapy for HCC.     

Systematic analysis of genes involved in oral cancer metastasis to lymph nodes

Oral cancer remains a deadly disease worldwide. Lymph node metastasis and invasion is one of the causes of death from oral cancer. Elucidating the mechanism of oral cancer lymph node metastasis and identifying critical regulatory genes are important for the treatment of this disease. This study aimed to identify differentially expressed genes (gene signature) and pathways that contribute to oral cancer metastasis to lymph nodes. The GSE70604-associated study compared gene profiles in lymph nodes with metastasis of oral cancer to those of normal lymph nodes. The GSE2280-associated study compared gene profiles in primary tumor of oral cancer with lymph node metastasis to those in tumors without lymph node metastasis. There are 28 common differentially expressed genes (DEGs) showing consistent changes in both datasets in overlapping analysis. GO biological process and KEGG pathway analysis of these 28 DEGs identified the gene signature CCND1, JUN and SPP1, which are categorized as key regulatory genes involved in the focal adhesion pathway. Silencing expression of CCND1, JUN and SPP1 in the human oral cancer cell line OECM-1 confirmed that those genes play essential roles in oral cancer cell invasion. Analysis of clinical samples of oral cancer found a strong correlation of these genes with short survival, especially JUN expression associated with metastasis. Our study identified a unique gene signature – CCND1, JUN and SPP1 – which may be involved in oral cancer lymph node metastasis.     

The predictive value of a preoperative systemic immune-inflammation index and prognostic nutritional index in patients with esophageal squamous cell carcinoma

Growing evidence indicates that systemic inflammation response and malnutrition status are correlated with survival in certain types of solid tumors. The aim of this study is to evaluate the association between the systemic immune-inflammation index (SII) and prognostic nutritional index (PNI) and overall survival (OS) in patients with esophageal squamous cell carcinoma (ESCC) after esophagectomy. A consecutive series of 655 patients with resected ESCC who underwent esophagectomy were enrolled in the retrospective study. The preoperative SII was defined as platelet × neutrophil/lymphocyte counts. The PNI was calculated as albumin concentration (g/L) + 5 × total lymphocyte count (10⁹/L). The optimal cut-off values of SII, neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio (PLR), and PNI were determined by receiver operating characteristic analysis. Survival analysis was performed using the Kaplan–Meier method with a log-rank test, followed by a multivariate Cox proportional hazards model. A high SII was significantly related to tumor size, histological type, invasion depth, and TNM stage (p < 0.05). A low PNI was significantly associated with age, tumor size, invasion depth, lymph node metastasis, and TNM stage (p < 0.05). Univariate analysis revealed that age, smoking history, tumor size, invasion depth, lymph node metastasis, SII, NLR, PLR, and PNI were predictors of OS (p < 0.05). Multivariate analysis identified age (p = 0.041), tumor size (p = 0.016), invasion depth (p < 0.001), lymph node metastasis (p < 0.001), SII (p = 0.033), and PNI (p = 0.022) as independent prognostic factors correlated with OS. There was a significant inverse relationship between the SII and PNI (r = 0.309; p < 0.001). The predictive value increased when the SII and PNI were considered in combination. Our results demonstrate that the preoperative high SII and low PNI are powerful indicators of aggressive biology and poor prognosis for patients with ESCC. The combination of SII and PNI can enhance the accuracy of prognosis.     

Five key lncRNAs considered as prognostic targets for predicting pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis, and the 5‐year survival rate was only 7.7%. To improve prognosis, a screening biomarker for early diagnosis of pancreatic cancer is in urgent need. Long non‐coding RNA (lncRNA) expression profiles as potential cancer prognostic biomarkers play critical roles in development of tumorigenesis and metastasis of cancer. However, lncRNA signatures in predicting the survival of a patient with PDAC remain unknown. In the current study, we try to identify potential lncRNA biomarkers and their prognostic values in PDAC. LncRNAs expression profiles and corresponding clinical information for 182 cases with PDAC were acquired from The Cancer Genome Atlas (TCGA). A total of 14 470 lncRNA were identified in the cohort, and 175 PDAC patients had clinical variables. We obtained 108 differential expressed lncRNA via R packages. Univariate and multivariate Cox proportional hazards regression, lasso regression was performed to screen the potential prognostic lncRNA. Five lncRNAs have been recognized to significantly correlate with OS. We established a linear prognostic model of five lncRNA (C9orf139, MIR600HG, RP5‐965G21.4, RP11‐436K8.1, and CTC‐327F10.4) and divided patients into high‐ and low‐risk group according to the prognostic index. The five lncRNAs played independent prognostic biomarkers of OS of PDAC patients and the AUC of the ROC curve for the five lncRNAs signatures prediction 5‐year survival was 0.742. In addition, targeted genes of MIR600HG, C9orf139, and CTC‐327F10.4 were explored and functional enrichment was also conducted. These results suggested that this five‐lncRNAs signature could act as potential prognostic biomarkers in the prediction of PDAC patient's survival.     

Tamoxifen therapy benefit predictive signature combining with prognostic signature in surgical-only ER-positive breast cancer

Tamoxifen treatment is important assistant for estrogen-receptor-positive breast cancer (BRCA) after resection. This study aimed to identify signatures for predicting the prognosis of patients with BRCA after tamoxifen treatment. Data of gene-specific DNA methylation (DM), as well as the corresponding clinical data for the patients with BRCA, were obtained from The Cancer Genome Atlas and followed by systematic bioinformatics analyses. After mapping these DM CPG sites onto genes, we finally obtained 352 relapse-free survival (RFS) associated DM genes, with which 61,776 gene pairs were combined, including 1,614 gene pairs related to RFS. An 11 gene-pair signature was identified to cluster the 189 patients with BRCA into the surgical low-risk group (136 patients) and high-risk group (53 patients). Then, we further identified a tamoxifen-predictive signature that could classify surgical high-risk patients with significant differences on RFS. Combining surgical-only prognostic signature and tamoxifen-predictive signature, patients were clustered into surgical-only low-risk group, tamoxifen nonbenefit group, and tamoxifen benefit group. In conclusion, we identified that the gene pair PDHA2–APRT could serve as a potential prognostic biomarker for patients with BRCA after tamoxifen treatment.