Quantitative or qualitative transcriptional diagnostic signatures? A case study for colorectal cancer

Exploring the influence of pre-analytical variables on gene expression measurements and relative expression orderings in cancer research

Gene expression profiling is an effective method for identifying predictive and prognostic biomarkers. However, measurements are prone to uncertainty and errors due to various pre-analytical variables. Systematic evaluating effects of these variables on gene expression measurements and relative expression orderings (REOs) of gene pairs, is necessary. A total of 18 datasets were collected, comprising over 800 paired samples. These paired samples were utilized to assess the impact of pre-analytical variables on gene expression measurements and REOs, including sampling methods, tumor sample heterogeneity, fixed time delays, preservation conditions, degradation levels, library preparation kits, amplification kits, RNA quantity, measuring platforms, and laboratory sites at single and multi-variable level. Low-quality samples served as the case group, while paired high-quality samples constituted the control group. In both single and multiple variable analyses, comparing each case sample to paired control sample revealed thousands of genes exhibited a twofold change in expression values. In contrast, on average, 82% and 76% of gene pairs keep consistent REO pattern between paired samples in single-variable and multi-variable analyses, respectively. Notably, the rate steadily increased after excluding gene pairs with the closest expression levels. Statistical analyses shown a higher proportion of differentially expressed genes (DEGs) than that of reversed gene pairs between case and control groups in both single-variable and multi-variable analyses. Furthermore, the proportion of reversal gene pairs among all gene pairs involving DEGs remained below 20% in the majority of comparisons. Our research demonstrates that REOs exhibit higher robustness under the influence of pre-analytical variables. These findings indicate the potential of the REOs-based approach in transcriptomics research and its applicability for biomarker studies.

Less is more: relative rank is more informative than absolute abundance for compositional NGS data

High-throughput gene expression data have been extensively generated and utilized in biological mechanism investigations, biomarker detection, disease diagnosis and prognosis. These applications encompass not only bulk transcriptome, but also single cell RNA-seq data. However, extracting reliable biological information from transcriptome data remains challenging due to the constrains of Compositional Data Analysis. Current data preprocessing methods, including dataset normalization and batch effect correction, are insufficient to address these issues and improve data quality for downstream analysis. Alternatively, qualification methods focusing on the relative order of gene expression (ROGER) are more informative than the quantification methods that rely on gene expression abundance. The Pairwise Analysis of Gene expression method is an enhancement of ROGER, designed for data integration in either sample space or feature space. In this review, we summarize the methods applied to transcriptome data analysis and discuss their potentials in predicting clinical outcomes.

Qualitative Transcriptional Signature for Predicting the Pathological Response of Colorectal Cancer to FOLFIRI Therapy

FOLFIRI (5-FU, leucovorin, irinotecan) is the first-line chemotherapy for metastatic colorectal cancer (mCRC), but response rates are under 50%. This study aimed to develop a predictive signature for FOLFIRI response in mCRC patients. Firstly, Spearman's rank correlation and Wilcoxon rank-sum test were used to select chemotherapy response genes and gene pairs, respectively. Then, an optimization procedure was used to determine the final signature. A predictive signature consisting of three gene pairs (3-GPS) was identified. In the training set, 3-GPS achieved an accuracy of 0.94. In a validation set of 60 samples, predicted responders had significantly better progression-free survival than the predicted non-responders (HR = 0.47, p = 0.01). A comparable result was observed in an additional validation set of 27 samples (HR = 0.06, p = 0.02). The co-expressed genes of the signature were enriched in pathways associated with the immunotherapy response, and they interacted extensively with FOLFIRI-related genes. Notably, the expression of signature genes significantly correlated with various immune cell types, including plasma cells and memory-resting CD4+ T cells. In conclusion, the REO-based signature effectively identifies mCRC patients likely to benefit from FOLFIRI. Furthermore, these signature genes may play a crucial role in the chemotherapy.

Randomized Phase 3 Trial of the Hypoxia Modifier Nimorazole Added to Radiation Therapy With Benefit Assessed in Hypoxic Head and Neck Cancers Determined Using a Gene Signature (NIMRAD)

PURPOSE

Tumor hypoxia is an adverse prognostic factor in head and neck squamous cell carcinoma (HNSCC). We assessed whether patients with hypoxic HNSCC benefited from the addition of nimorazole to definitive intensity modulated radiation therapy (IMRT).

METHODS AND MATERIALS

NIMRAD was a phase 3, multicenter, placebo-controlled, double-anonymized trial of patients with HNSCC unsuitable for concurrent platinum chemotherapy or cetuximab with definitive IMRT (NCT01950689). Patients were randomized 1:1 to receive IMRT (65 Gy in 30 fractions over 6 weeks) plus nimorazole (1.2 g/m2 daily, before IMRT) or placebo. The primary endpoint was freedom from locoregional progression (FFLRP) in patients with hypoxic tumors, defined as greater than or equal to the median tumor hypoxia score of the first 50 patients analyzed (≥0.079), using a validated 26-gene signature. The planned sample size was 340 patients, allowing for signature generation in 85% and an assumed hazard ratio (HR) of 0.50 for nimorazole effectiveness in the hypoxic group and requiring 66 locoregional failures to have 80% power in a 2-tail log-rank test at the 5% significance level.

RESULTS

Three hundred thirty-eight patients were randomized by 19 centers in the United Kingdom from May 2014 to May 2019, with a median follow-up of 3.1 years (95% CI, 2.9-3.4). Hypoxia scores were available for 286 (85%). The median patient age was 73 years (range, 44-88; IQR, 70-76). There were 36 (25.9%) locoregional failures in the hypoxic group, in which nimorazole + IMRT did not improve FFLRP (adjusted HR, 0.72; 95% CI, 0.36-1.44; P = .35) or overall survival (adjusted HR, 0.96; 95% CI, 0.53-1.72; P = .88) compared with placebo + IMRT. Similarly, nimorazole + IMRT did not improve FFLRP or overall survival in the whole population. In total (N = 338), 73% of patients allocated nimorazole adhered to the drug for ≥50% of IMRT fractions. Nimorazole + IMRT caused more acute nausea compared with placebo + IMRT (Common Terminology Criteria for Adverse Events version 4.0 G1+2: 56.6% vs 42.4%, G3: 10.1% vs 5.3%, respectively; P < .05).

CONCLUSIONS

Addition of the hypoxia modifier nimorazole to IMRT for locally advanced HNSCC in older and less fit patients did not improve locoregional control or survival.

Individualized detection of TMPRSS2-ERG fusion status in prostate cancer: a rank-based qualitative transcriptome signature

BACKGROUND

TMPRSS2-ERG (T2E) fusion is highly related to aggressive clinical features in prostate cancer (PC), which guides individual therapy. However, current fusion prediction tools lacked enough accuracy and biomarkers were unable to be applied to individuals across different platforms due to their quantitative nature. This study aims to identify a transcriptome signature to detect the T2E fusion status of PC at the individual level.

METHODS

Based on 272 high-throughput mRNA expression profiles from the Sboner dataset, we developed a rank-based algorithm to identify a qualitative signature to detect T2E fusion in PC. The signature was validated in 1223 samples from three external datasets (Setlur, Clarissa, and TCGA).

RESULTS

A signature, composed of five mRNAs coupled to ERG (five ERG-mRNA pairs, 5-ERG-mRPs), was developed to distinguish T2E fusion status in PC. 5-ERG-mRPs reached 84.56% accuracy in Sboner dataset, which was verified in Setlur dataset (n = 455, accuracy = 82.20%) and Clarissa dataset (n = 118, accuracy = 81.36%). Besides, for 495 samples from TCGA, two subtypes classified by 5-ERG-mRPs showed a higher level of significance in various T2E fusion features than subtypes obtained through current fusion prediction tools, such as STAR-Fusion.

CONCLUSIONS

Overall, 5-ERG-mRPs can robustly detect T2E fusion in PC at the individual level, which can be used on any gene measurement platform without specific normalization procedures. Hence, 5-ERG-mRPs may serve as an auxiliary tool for PC patient management.

Hierarchical identification of a transcriptional panel for the histological diagnosis of lung neuroendocrine tumors

Background: Lung cancer is a complex disease composed of neuroendocrine (NE) and non-NE tumors. Accurate diagnosis of lung cancer is essential in guiding therapeutic management. Several transcriptional signatures have been reported to distinguish between adenocarcinoma (ADC) and squamous cell carcinoma (SCC) belonging to non-NE tumors. This study aims to identify a transcriptional panel that could distinguish the histological subtypes of NE tumors to complement the morphology-based classification of an individual.

Methods: A public dataset with NE subtypes, including 21 small-cell lung cancer (SCLC), 56 large-cell NE carcinomas (LCNECs), and 24 carcinoids (CARCIs), and non-NE subtypes, including 85 ADC and 61 SCC, was used as a training set. In the training set, consensus clustering was first used to filter out the samples whose expression patterns disagreed with their histological subtypes. Then, a rank-based method was proposed to develop a panel of transcriptional signatures for determining the NE subtype for an individual, based on the within-sample relative gene expression orderings of gene pairs. Twenty-three public datasets with a total of 3,454 samples, which were derived from fresh-frozen, formalin-fixed paraffin-embedded, biopsies, and single cells, were used for validation. Clinical feasibility was tested in 10 SCLC biopsy specimens collected from cancer hospitals via bronchoscopy.

Results: The NEsubtype-panel was composed of three signatures that could distinguish NE from non-NE, CARCI from non-CARCI, and SCLC from LCNEC step by step and ultimately determine the histological subtype for each NE sample. The three signatures achieved high average concordance rates with 97.31%, 98.11%, and 90.63%, respectively, in the 23 public validation datasets. It is worth noting that the 10 clinic-derived SCLC samples diagnosed via immunohistochemical staining were also accurately predicted by the NEsubtype-panel. Furthermore, the subtype-specific gene expression patterns and survival analyses provided evidence for the rationality of the reclassification by the NEsubtype-panel.

Conclusion: The rank-based NEsubtype-panel could accurately distinguish lung NE from non-NE tumors and determine NE subtypes even in clinically challenging samples (such as biopsy). The panel together with our previously reported signature (KRT5-AGR2) for SCC and ADC would be an auxiliary test for the histological diagnosis of lung cancer.

lncRNA pair as candidate diagnostic signature for colorectal cancer based on the within-sample relative expression levels

Background

Early diagnosis of colorectal cancer could significantly improve the prognosis and reduce mortality. However, indeterminate diagnosis is often met in pathology diagnosis in biopsy samples. Abnormal expression of long non-coding RNA (lncRNA) is associated with the initiation and progression of colorectal cancer. It is of great value and clinical significance to explore lncRNAs as candidate diagnostic biomarkers in colorectal cancer.

Methods

Based on the within-sample relative expression levels of lncRNA pairs, we identified a group of candidate diagnostic biomarkers for colorectal cancer. In addition, we validated it in independent datasets produced by different laboratories and different platforms. We also tested it in colorectal cancer tissue samples using quantitative real-time polymerase chain reaction (RT-qPCR).

Results

A biomarker consisting of six lncRNA pairs including nine lncRNAs was identified for the diagnosis of colorectal cancer. For a total of 950 cancer samples and 247 non-cancer samples, both of the sensitivity and specificity could achieve approximately 90%. For adenoma samples, the accuracy could achieve 73%. For normal tissues from inflammatory bowel disease patients, 93% (14/15) were correctly classified as non-cancer. Furthermore, the lncRNA pair biomarker showed excellent performance in all clinical stages; even in the early stage, the accuracy could achieve 87% and 82% in stage I and II. Meanwhile, the biomarker was also robust to the microsatellite instability status. More importantly, we measured the biomarker in 35 colorectal cancer and 30 cancer-adjacent tissue samples using quantitative real-time polymerase chain reaction (RT-qPCR). The accuracy could achieve 93.3% (70/75). Specially, even in early-stage tumors (I and II), the accuracy could also achieve 90.9% (30/33). The enrichment analysis revealed that these lncRNAs were involved in highly associated cancer pathways and immune-related pathways. Immune analysis showed that these marker lncRNAs were associated with multiple immune cells, implying that they might be involved in the regulation of immune cell functions in colorectal cancer. Most of the biomarker lncRNAs were also differentially expressed between the mutant group and wild-type group of colorectal cancer driver genes.

Conclusion

We identified and validated six lncRNA pairs including nine lncRNAs as a biomarker for assisting in the diagnosis of colorectal cancer.

The prognostic and clinical significance of IFI44L aberrant downregulation in patients with oral squamous cell carcinoma

Background

It is a basic task in high-throughput gene expression profiling studies to identify differentially expressed genes (DEGs) between two phenotypes. RankComp, an algorithm, could analyze the highly stable within-sample relative expression orderings (REOs) of gene pairs in a particular type of human normal tissue that are widely reversed in the cancer condition, thereby detecting DEGs for individual disease samples measured by a particular platform.

Methods

In the present study, Gene Expression Omnibus (GEO) Series (GSE) GSE75540, GSE138206 were downloaded from GEO, by analyzing DEGs in oral squamous cell carcinoma based on online datasets using the RankComp algorithm, using the Kaplan-Meier survival analysis and Cox regression analysis to survival analysis, Gene Set Enrichment Analysis (GSEA) to explore the potential molecular mechanisms underlying.

Results

We identified 6 reverse gene pairs with stable REOs. All the 12 genes in these 6 reverse gene pairs have been reported to be associated with cancers. Notably, lower Interferon Induced Protein 44 Like (IFI44L) expression was associated with poorer overall survival (OS) and Disease-free survival (DFS) in oral squamous cell carcinoma patients, and IFI44L expression showed satisfactory predictive efficiency by receiver operating characteristic (ROC) curve. Moreover, low IFI44L expression was identified as risk factors for oral squamous cell carcinoma patients’ OS. IFI44L downregulation would lead to the activation of the FRS-mediated FGFR1, FGFR3, and downstream signaling pathways, and might play a role in the PI3K-FGFR cascades.

Conclusions

Collectively, we identified 6 reverse gene pairs with stable REOs in oral squamous cell carcinoma, which might serve as gene signatures playing a role in the diagnosis in oral squamous cell carcinoma. Moreover, high expression of IFI44L, one of the DEGs in the 6 reverse gene pairs, might be associated with favorable prognosis in oral squamous cell carcinoma patients and serve as a tumor suppressor by acting on the FRS-mediated FGFR signaling.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-09058-y.

Discovering a qualitative transcriptional signature of homologous recombination defectiveness for prostate cancer

The discovery of homologous recombination deficiency (HRD) biomarkers in prostate cancer is important for patients who will benefit from poly ADP-ribose polymerase inhibitor (PARPi). Here, we developed a transcriptional homologous recombination defectiveness (HRDness) signature, comprising 16 gene pairs (16-GPS), for prostate cancer by a relative expression ordering (REO)-based discovery procedure. Subsequently, two newly subtypes classified by 16-GPS showed a higher significance level in various clinicopathological and HRD features than subtypes obtained by other methods, such as HRDetect. HRDness subtype also displayed more aggressive features and higher genomics scores than non-HRDness in three independent datasets. HRDness prostate cancer cells were more sensitive to PARPi than non-HRDness. Moreover, the HRDness samples showed distinct multi-omics characteristics related to homologous recombination repair function loss. Overall, the newly proposed qualitative signature can robustly determine the HRD status for prostate cancer at the personalized level, and especially be an auxiliary tool for PARPi treatment strategy.

•

16 gene pairs (16-GPS) could predict HRDness for prostate cancer at individual level

•

HRDness samples classified by 16-GPS showed HRD molecular and clinical features

•

HRDness cells classified by 16-GPS tend to be sensitive to PARPi

Emergence of Cardiac Glycosides as Potential Drugs: Current and Future Scope for Cancer Therapeutics

Cardiac glycosides are natural sterols and constitute a group of secondary metabolites isolated from plants and animals. These cardiotonic agents are well recognized and accepted in the treatment of various cardiac diseases as they can increase the rate of cardiac contractions by acting on the cellular sodium potassium ATPase pump. However, a growing number of recent efforts were focused on exploring the antitumor and antiviral potential of these compounds. Several reports suggest their antitumor properties and hence, today cardiac glycosides (CG) represent the most diversified naturally derived compounds strongly recommended for the treatment of various cancers. Mutated or dysregulated transcription factors have also gained prominence as potential therapeutic targets that can be selectively targeted. Thus, we have explored the recent advances in CGs mediated cancer scope and have considered various signaling pathways, molecular aberration, transcription factors (TFs), and oncogenic genes to highlight potential therapeutic targets in cancer management.

Stratification of patients with clear cell renal cell carcinoma to facilitate drug repositioning

Clear cell renal cell carcinoma (ccRCC) is the most common histological type of kidney cancer and has high heterogeneity. Stratification of ccRCC is important since distinct subtypes differ in prognosis and treatment. Here, we applied a systems biology approach to stratify ccRCC into three molecular subtypes with different mRNA expression patterns and prognosis of patients. Further, we developed a set of biomarkers that could robustly classify the patients into each of the three subtypes and predict the prognosis of patients. Then, we reconstructed subtype-specific metabolic models and performed essential gene analysis to identify the potential drug targets. We identified four drug targets, including SOAT1, CRLS1, and ACACB, essential in all the three subtypes and GPD2, exclusively essential to subtype 1. Finally, we repositioned mitotane, an FDA-approved SOAT1 inhibitor, to treat ccRCC and showed that it decreased tumor cell viability and inhibited tumor cell growth based on in vitro experiments.

•

Three consistent molecular ccRCC subtypes were found to guide patients' prognoses

•

REOs-based biomarker was developed to robustly classify patients at individual level

•

SOAT1 is identified as a common drug target for all ccRCC subtypes

•

Mitotane was repositioned treatment of ccRCC via inhibiting SOAT1

Qualitative diagnostic signature for pancreatic ductal adenocarcinoma based on the within-sample relative expression orderings

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) accounts for about 90% of pancreatic cancer, which is one of the most aggressive malignant neoplasms with a 9.3% five-year survival rate. The pathological biopsy is the current golden standard for confirming suspicious lesions of PDAC, but it is not entirely reliable because of the insufficient sampling amount and inaccurate sampling location. Therefore, developing a robust signature to aid the accurate diagnosis of PDAC is critical.

METHODS

Based on the within-sample relative expression orderings of gene pairs, we identified a qualitative signature to discriminate both PDAC and adjacent samples from both chronic pancreatitis and normal samples in the training datasets and validated it in other independent datasets produced by different laboratories with different measuring platforms.

RESULTS

A six-gene-pair signature was identified in the training data and validated in eight independent datasets. For surgical samples, 96.63% of 356 PDAC tissues, 100% of 11 pancreatitis tissues of non-cancer patients, and 23 of 24 normal pancreatic tissues were correctly classified. Especially, 59 of 60 cancer-adjacent normal tissues of PDAC patients were correctly identified as PDAC. For biopsy samples, all of 11 PDAC biopsy tissues were correctly classified as PDAC.

CONCLUSION

The signature can distinguish both PDAC and PDAC-adjacent normal tissues from both chronic pancreatitis and normal tissues of non-cancer patients even when the sampling locations are inaccurate, which can aid the diagnosis of PDAC.

A Qualitative Transcriptional Signature for the Risk Assessment of Precancerous Colorectal Lesions

It is meaningful to assess the risk of cancer incidence among patients with precancerous colorectal lesions. Comparing the within-sample relative expression orderings (REOs) of colorectal cancer patients measured by multiple platforms with that of normal colorectal tissues, a qualitative transcriptional signature consisting of 1,840 gene pairs was identified in the training data. Within an evaluation dataset of 16 active and 18 inactive (remissive) ulcerative colitis subjects, the median incidence risk score of colorectal carcinoma was 0.6402 in active ulcerative colitis subjects, significantly higher than that in remissive subjects (0.3114). Evaluation of two other independent datasets yielded similar results. Moreover, we found that the score significantly positively correlated with the degree of dysplasia in the case of colorectal adenomas. In the merged dataset, the median incidence risk score was 0.9027 among high-grade adenoma samples, significantly higher than that among low-grade adenomas (0.8565). In summary, the developed incidence risk score could well predict the incidence risk of precancerous colorectal lesions and has value in clinical application.

Identification of the origin of brain metastases based on the relative methylation orderings of CpG sites

Accurate diagnosis of the origin of brain metastases (BMs) is crucial for tailoring an effective therapy to improve patients' prognosis. BMs of unknown origin account for approximately 2-14% of patients with BMs. Hence, the aim of this study was to identify the original cancer type of BMs based on their DNA methylation profiles. The DNA methylation profiles of glioma (GM), BM, and seven other types of primary cancers were collected. In comparison with GM, the reversal CpG site pairs were identified for each of the seven other types of primary cancers based on the within-sample relative methylation orderings (RMOs) of the CpG sites. Then, using the reversal CpG site pairs, GMs were distinguished from BMs and the seven other types of primary cancers. All 61 of the GM samples were correctly identified as GM. The cancer type was also identified for the non-GM samples. For the seven other types of primary cancers, greater than 93% of samples of each cancer type were correctly identified as their corresponding cancer type, except for breast cancer, which had an 88% accuracy. For 133 BM samples, 132 BM samples were identified as non-GM, and 95% of the 133 BM samples were correctly classified into their corresponding original cancer types. The RMO-based method can accurately identify the origin of BMs, which is important for precision treatment.

The Effects of Age, Cigarette Smoking, Sex, and Race on the Qualitative Characteristics of Lung Transcriptome

The within-sample relative expression orderings (REOs) of genes, which are stable qualitative transcriptional characteristics, can provide abundant information for a disease. Methods based on REO comparisons have been proposed for identifying differentially expressed genes (DEGs) at the individual level and for detecting disease-associated genes based on one-phenotype disease data by reusing data of normal samples from other sources. Here, we evaluated the effects of common potential confounding factors, including age, cigarette smoking, sex, and race, on the REOs of gene pairs within normal lung tissues transcriptome. Our results showed that age has little effect on REOs within lung tissues. We found that about 0.23% of the significantly stable REOs of gene pairs in nonsmokers' lung tissues are reversed in smokers' lung tissues, introduced by 344 DEGs between the two groups of samples (RankCompV2, FDR <0.05), which are enriched in metabolism of xenobiotics by cytochrome P450, glutathione metabolism, and other pathways (hypergeometric test, FDR <0.05). Comparison between the normal lung tissue samples of males and females revealed fewer reversal REOs introduced by 24 DEGs between the sex groups, among which 19 DEGs are located on sex chromosomes and 5 DEGs involving in spermatogenesis and regulation of oocyte are located on autosomes. Between the normal lung tissue samples of white and black people, we identified 22 DEGs (RankCompV2, FDR <0.05) which introduced a few reversal REOs between the two races. In summary, the REO-based study should take into account the confounding factors of cigarette smoking, sex, and race.

Early Diagnosis of Hepatocellular Carcinoma Using Machine Learning Method

Hepatocellular carcinoma (HCC) is a serious cancer which ranked the fourth in cancer-related death worldwide. Hence, more accurate diagnostic models are urgently needed to aid the early HCC diagnosis under clinical scenarios and thus improve HCC treatment and survival. Several conventional methods have been used for discriminating HCC from cirrhosis tissues in patients without HCC (CwoHCC). However, the recognition successful rates are still far from satisfactory. In this study, we applied a computational approach that based on machine learning method to a set of microarray data generated from 1091 HCC samples and 242 CwoHCC samples. The within-sample relative expression orderings (REOs) method was used to extract numerical descriptors from gene expression profiles datasets. After removing the unrelated features by using maximum redundancy minimum relevance (mRMR) with incremental feature selection, we achieved “11-gene-pair” which could produce outstanding results. We further investigated the discriminate capability of the “11-gene-pair” for HCC recognition on several independent datasets. The wonderful results were obtained, demonstrating that the selected gene pairs can be signature for HCC. The proposed computational model can discriminate HCC and adjacent non-cancerous tissues from CwoHCC even for minimum biopsy specimens and inaccurately sampled specimens, which can be practical and effective for aiding the early HCC diagnosis at individual level.

A Qualitative Transcriptional Signature for Predicting Extreme Resistance of ER-Negative Breast Cancer to Paclitaxel, Doxorubicin, and Cyclophosphamide Neoadjuvant Chemotherapy

For estrogen receptor (ER)-negative breast cancer patients, paclitaxel (P), doxorubicin (A) and cyclophosphamide (C) neoadjuvant chemotherapy (NAC) is the standard therapeutic regimen. Pathologic complete response (pCR) and residual disease (RD) are common surrogate measures of chemosensitivity. After NAC, most patients still have RD; of these, some partially respond to NAC, whereas others show extreme resistance and cannot benefit from NAC but only suffer complications resulting from drug toxicity. Here we developed a qualitative transcriptional signature, based on the within-sample relative expression ordering (REO) of gene pairs, to identify extremely resistant samples to PAC NAC. Using gene expression data for ER-negative breast cancer patients including 113 pCR samples and 137 RD samples from four datasets, we selected 61 gene pairs with reversal REO patterns between the two groups as the resistance signature, denoted as NR61. Samples with more than 37 signature gene pairs that had the same REO patterns within the extremely resistant group were defined as having extreme resistance; otherwise, they were considered responders. In the GSE25055 and GSE25065 dataset, the NR61 signature could correctly identify 44 (97.8%) of the 45 pCR samples and 22 (95.7%) of the 23 pCR samples as responder samples, respectively; it also identified 13 (16.9%) of 77 RD samples and 8 (21.1%) of 38 RD samples as extremely resistant samples, respectively. Survival analysis showed that the distant relapse-free survival (DRFS) time of the 14 extremely resistant cases was significantly shorter than that of the 108 responders (P < 0.01; HR = 3.84; 95% CI = 1.91–7.70) in GSE25055. Similar results were obtained in GSE25065. Moreover, in the integrated data of the two datasets with 94 responders and 21 extremely resistant samples identified from RD patients, the former had significantly longer DRFS than the latter (P < 0.01; HR = 2.22; 95% CI = 1.26–3.90). In summary, our signature could effectively identify patients who completely respond to PAC NAC, as well as cases of extreme resistance, which can assist decision-making on the clinical therapy for these patients.

Classification of clear cell renal cell carcinoma based on PKM alternative splicing

Clear cell renal cell carcinoma (ccRCC) accounts for 70-80% of kidney cancer diagnoses and displays high molecular and histologic heterogeneity. Hence, it is necessary to reveal the underlying molecular mechanisms involved in progression of ccRCC to better stratify the patients and design effective treatment strategies. Here, we analyzed the survival outcome of ccRCC patients as a consequence of the differential expression of four transcript isoforms of the pyruvate kinase muscle type (PKM). We first extracted a classification biomarker consisting of eight gene pairs whose within-sample relative expression orderings (REOs) could be used to robustly classify the patients into two groups with distinct molecular characteristics and survival outcomes. Next, we validated our findings in a validation cohort and an independent Japanese ccRCC cohort. We finally performed drug repositioning analysis based on transcriptomic expression profiles of drug-perturbed cancer cell lines and proposed that paracetamol, nizatidine, dimethadione and conessine can be repurposed to treat the patients in one of the subtype of ccRCC whereas chenodeoxycholic acid, fenoterol and hexylcaine can be repurposed to treat the patients in the other subtype.

Qualitative transcriptional signature for predicting pathological response of colorectal cancer to FOLFOX therapy

FOLFOX (5‐fluorouracil, leucovorin and oxaliplatin) is one of the main chemotherapy regimens for colorectal cancer (CRC), but only half of CRC patients respond to this regimen. Using gene expression profiles of 96 metastatic CRC patients treated with FOLFOX, we first selected gene pairs whose within‐sample relative expression orderings (REO) were significantly associated with the response to FOLFOX using the exact binomial test. Then, from these gene pairs, we applied an optimization procedure to obtain a subset that achieved the largest F‐score in predicting pathological response of CRC to FOLFOX. The REO‐based qualitative transcriptional signature, consisting of five gene pairs, was developed in the training dataset consisting of 96 samples with an F‐score of 0.90. In an independent test dataset consisting of 25 samples with the response information, an F‐score of 0.82 was obtained. In three other independent survival datasets, the predicted responders showed significantly better progression‐free survival than the predicted non‐responders. In addition, the signature showed a better predictive performance than two published FOLFOX signatures across different datasets and is more suitable for CRC patients treated with FOLFOX than 5‐fluorouracil‐based signatures. In conclusion, the REO‐based qualitative transcriptional signature can accurately identify metastatic CRC patients who may benefit from the FOLFOX regimen.

A qualitative transcriptional signature for the histological reclassification of lung squamous cell carcinomas and adenocarcinomas

Background

Targeted therapy for non-small cell lung cancer is histology dependent. However, histological classification by routine pathological assessment with hematoxylin-eosin staining and immunostaining for poorly differentiated tumors, particularly those from small biopsies, is still challenging. Additionally, the effectiveness of immunomarkers is limited by technical inconsistencies of immunostaining and lack of standardization for staining interpretation.

Results

Using gene expression profiles of pathologically-determined lung adenocarcinomas and squamous cell carcinomas, denoted as pADC and pSCC respectively, we developed a qualitative transcriptional signature, based on the within-sample relative gene expression orderings (REOs) of gene pairs, to distinguish ADC from SCC. The signature consists of two genes, KRT5 and AGR2, which has the stable REO pattern of KRT5 > AGR2 in pSCC and KRT5 < AGR2 in pADC. In the two test datasets with relative unambiguous NSCLC types, the apparent accuracy of the signature were 94.44 and 98.41%, respectively. In the other integrated dataset for frozen tissues, the signature reclassified 4.22% of the 805 pADC patients as SCC and 12% of the 125 pSCC patients as ADC. Similar results were observed in the clinical challenging cases, including FFPE specimens, mixed tumors, small biopsy specimens and poorly differentiated specimens. The survival analyses showed that the pADC patients reclassified as SCC had significantly shorter overall survival than the signature-confirmed pADC patients (log-rank p = 0.0123, HR = 1.89), consisting with the knowledge that SCC patients suffer poor prognoses than ADC patients. The proliferative activity, subtype-specific marker genes and consensus clustering analyses also supported the correctness of our signature.

Conclusions

The non-subjective qualitative REOs signature could effectively distinguish ADC from SCC, which would be an auxiliary test for the pathological assessment of the ambiguous cases.

A Qualitative Transcriptional Signature for Predicting Recurrence Risk for High-Grade Serous Ovarian Cancer Patients Treated With Platinum-Taxane Adjuvant Chemotherapy

Resistance to platinum and taxane adjuvant chemotherapy (ACT) is the main cause of the recurrence and poor prognosis of high-grade serous ovarian cancer (HGS-OvCa) patients receiving platinum-taxane ACT after surgery. However, currently reported quantitative transcriptional signatures, which are commonly based on risk scores summarized from gene expression, are unsuitable for clinical application because of their high sensitivity to experimental batch effects and quality uncertainties of clinical samples. Using 226 samples of HGS-OvCa patients receiving platinum-taxane ACT in TCGA, we developed a qualitative transcriptional signature, consisting of four gene pairs whose within-samples relative expression orderings could robustly predict patient recurrence-free survival (RFS). In two independent test datasets, the predicted non-responders had significantly shorter RFS than the predicted responders (log-rank p < 0.05). In a test dataset containing data for patient pathological response state, the signature reclassified 12 out of 22 pathological complete response patients as non-responders and two out of 16 pathological non-complete response patients as responders. Notably, the 12 predicted non-responders in the pathological complete response group had significantly shorter RFS than the predicted responders (log-rank p = 0.0122). This qualitative transcriptional signature, which is insensitive to experimental batch effects and quality uncertainties of clinical samples, can individually identify HGS-OvCa patients who are more likely to benefit from platinum-taxane adjuvant chemotherapy.

A qualitative transcriptional signature for predicting microsatellite instability status of right-sided Colon Cancer

Background

Microsatellite instability (MSI) accounts for about 15% of colorectal cancer and is associated with prognosis. Today, MSI is usually detected by polymerase chain reaction amplification of specific microsatellite markers. However, the instability is identified by comparing the length of microsatellite repeats in tumor and normal samples. In this work, we developed a qualitative transcriptional signature to individually predict MSI status for right-sided colon cancer (RCC) based on tumor samples.

Results

Using RCC samples, based on the relative expression orderings (REOs) of gene pairs, we extracted a signature consisting of 10 gene pairs (10-GPS) to predict MSI status for RCC through a feature selection process. A sample is predicted as MSI when the gene expression orderings of at least 7 gene pairs vote for MSI; otherwise the microsatellite stability (MSS). The classification performance reached the largest F-score in the training dataset. This signature was verified in four independent datasets of RCCs with the F-scores of 1, 0.9630, 0.9412 and 0.8798, respectively. Additionally, the hierarchical clustering analyses and molecular features also supported the correctness of the reclassifications of the MSI status by 10-GPS.

Conclusions

The qualitative transcriptional signature can be used to classify MSI status of RCC samples at the individualized level.

A qualitative transcriptional signature for determining the grade of colorectal adenocarcinoma

Histological grading (HG) is an important prognostic factor of colorectal adenocarcinoma (CRAC): the high-grade CRAC patients have poorer prognosis after tumor resection. Especially, the high-grade stage II CRAC patients are recommended to receive adjuvant chemotherapy. Due to the subjective nature of HG assessment, it is difficult to achieve consistency among pathologists, which brings patients uncertain grading outcomes and inappropriate treatments. We developed a qualitative transcriptional signature based on the within-sample relative expression orderings (REOs) of gene pairs to discriminate high-grade and low-grade CRAC. Using the stage II-III CRAC samples, we detected gene pairs with stable REOs in the high-grade samples and reversal stable REOs in the low-grade samples, and retained the gene pairs whose specific REO patterns were significantly associated with the disease-free survival of patients by univariate Cox regression model. Then, we used a forward-backward searching procedure to extract gene pairs with the highest concordance index as the final grading signature. Finally, 9 gene pairs (9-GPS) were developed to divide CRAC patients into high-grade and low-grade groups. With the signature, there were more differential expression characteristics between reclassified high-grade and low-grade groups. Significant difference of prognosis between the classified two group patients could be seen in four independent datasets. Additionally, genomic analyses showed that the classified high-grade groups were characterized by hypermutation while classified low-grade groups were characterized by frequent copy number alternations. In conclusion, the 9-GPS can provide an objective and robust grading assessment for CRAC patients, which could assist clinical treatment decision.

A qualitative transcriptional signature for the early diagnosis of colorectal cancer

Currently, using biopsy specimens for the early diagnosis of colorectal cancer (CRC) is not entirely reliable due to insufficient sampling amount and inaccurate sampling location. Thus, it is necessary to develop a signature that can accurately identify patients with CRC under these clinical scenarios. Based on the relative expression orderings of genes within individual samples, we developed a qualitative transcriptional signature to discriminate CRC tissues, including CRC adjacent normal tissues from non-CRC individuals. The signature was validated using multiple microarray and RNA sequencing data from different sources. In the training data, a signature consisting of 7 gene pairs was identified. It was well validated in both biopsy and surgical resection specimens from multiple datasets measured by different platforms. For biopsy specimens, 97.6% of 42 CRC tissues and 94.5% of 163 non-CRC (normal or inflammatory bowel disease) tissues were correctly classified. For surgically resected specimens, 99.5% of 854 CRC tissues and 96.3% of 81 CRC adjacent normal tissues were correctly identified as CRC. Notably, we additionally measured 33 CRC biopsy specimens by the Affymetrix platform and 13 CRC surgical resection specimens, with different proportions of tumor epithelial cells, ranging from 40% to 100%, by the RNA sequencing platform, and all these samples were correctly identified as CRC. The signature can be used for the early diagnosis of CRC, which is also suitable for minimum biopsy specimens and inaccurately sampled specimens, and thus has potential value for clinical application.

A Qualitative Transcriptional Signature for Predicting Recurrence Risk of Stage I-III Bladder Cancer Patients After Surgical Resection

Background: Previously reported transcriptional signatures for predicting the prognosis of stage I-III bladder cancer (BLCA) patients after surgical resection are commonly based on risk scores summarized from quantitative measurements of gene expression levels, which are highly sensitive to the measurement variation and sample quality and thus hardly applicable under clinical settings. It is necessary to develop a signature which can robustly predict recurrence risk of BLCA patients after surgical resection.

Methods: The signature is developed based on the within-sample relative expression orderings (REOs) of genes, which are qualitative transcriptional characteristics of the samples.

Results: A signature consisting of 12 gene pairs (12-GPS) was identified in training data with 158 samples. In the first validation dataset with 114 samples, the low-risk group of 54 patients had a significantly better overall survival than the high-risk group of 60 patients (HR = 3.59, 95% CI: 1.34~9.62, p = 6.61 × 10⁻⁰³). The signature was also validated in the second validation dataset with 57 samples (HR = 2.75 × 10⁰⁸, 95% CI: 0~Inf, p = 0.05). Comparison analysis showed that the transcriptional differences between the low- and high-risk groups were highly reproducible and significantly concordant with DNA methylation differences between the two groups.

Conclusions: The 12-GPS signature can robustly predict the recurrence risk of stage I-III BLCA patients after surgical resection. It can also aid the identification of reproducible transcriptional and epigenomic features characterizing BLCA metastasis.

Identification of potential key genes associated with osteosarcoma based on integrated bioinformatics analyses

Due to high rates of metastasis and poor clinical outcomes for patients, it is important to study the pathomechanisms of osteosarcoma. However, due to the fact that osteosarcoma shows significant interindividual variation and high heterogeneity, the identification of differentially expressed genes (DEGs) at the population level cannot answer many important questions related to osteosarcoma tumorigenesis. Therefore, a new strategy to identify dysregulated genes in osteosarcoma samples is required. The aim of this study was to improve our understanding of osteosarcoma pathogenesis by identifying genes with universal aberrant expression in osteosarcoma samples. Because the relative expression ordering of genes is stable in normal bone tissues but is disrupted in osteosarcoma tissues, we used the RankComp algorithm to identify DEGs in normal and osteosarcoma tissue samples. We then calculated the dysregulation frequency for each gene. Genes with deregulation frequencies above 80% were deemed to be universal DEGs. Next, coexpression, pathway enrichment, and protein-protein interaction network analyses were performed to characterize the functions of these genes. From 188 samples of osteosarcoma obtained from four datasets measured on different platforms, 51 universal DEGs were identified, including 4 universally upregulated genes and 47 universally downregulated genes. Genes that were differentially coexpressed with these universal DEGs were found to be enriched in 46 cancer-related pathways. In addition, functional and network analyses showed that genes with high dysregulation frequencies were involved in cancer-related functions. Thus, the commonly aberrant genes identified in osteosarcoma tissues may be important targets for osteosarcoma diagnosis and therapy.

A robust qualitative transcriptional signature for the correct pathological diagnosis of gastric cancer

BACKGROUND

Currently, pathological examination of gastroscopy biopsy specimens is the gold standard for gastric cancer (GC) diagnosis. However, it has a false-negative rate of 10-20% due to inaccurate sampling locations and/or insufficient sampling amount. A signature should be developed to aid the early diagnosis of GC using biopsy specimens even when they are sampled from inaccurate locations.

METHODS

We extracted a robust qualitative transcriptional signature, based on the within-sample relative expression orderings (REOs) of gene pairs, to discriminate both GC tissues and adjacent-normal tissues from non-GC gastritis, intestinal metaplasia and normal gastric tissues.

RESULTS

A signature consisting of two gene pairs for GC diagnosis was identified and validated in data of both biopsy specimens and surgical resection specimens pooled from publicly available datasets measured by different laboratories with different platforms. For gastroscopy biopsy specimens, 96.20% of 79 non-GC tissues were correctly identified as non-GC, and 96.84% of 158 GC tissues and six of seven adjacent-normal tissues were correctly identified as GC. For surgical resection specimens, 98.37% of 2560 GC tissues and 97.28% of 221 adjacent-normal tissues were correctly identified as GC. Especially, 97.67% of the 257 GC patients at stage I were exactly diagnosed as GC. We additionally measured 21 GC tissues from seven different GC patients, each with three specimens sampled from three tumor locations with different proportions of the tumor epithelial cell. All these GC tissues were correctly identified as GC, even when the proportion of the tumor epithelial cell was as low as 14%.

CONCLUSIONS

The qualitative transcriptional signature can distinguish both GC and adjacent-normal tissues from normal, gastritis and intestinal metaplasia tissues of non-GC patients even using inaccurately sampled biopsy specimens, which can be applied robustly at the individual level to aid the early GC diagnosis.

Identification and characterization of genes with absolute mRNA abundances changes in tumor cells with varied transcriptome sizes

BACKGROUND

The amount of RNA per cell, namely the transcriptome size, may vary under many biological conditions including tumor. If the transcriptome size of two cells is different, direct comparison of the expression measurements on the same amount of total RNA for two samples can only identify genes with changes in the relative mRNA abundances, i.e., cellular mRNA concentration, rather than genes with changes in the absolute mRNA abundances.

RESULTS

Our recently proposed RankCompV2 algorithm identify differentially expressed genes (DEGs) through comparing the relative expression orderings (REOs) of disease samples with that of normal samples. We reasoned that both the mRNA concentration and the absolute abundances of these DEGs must have changes in disease samples. In simulation experiments, this method showed excellent performance for identifying DEGs between normal and disease samples with different transcriptome sizes. Through analyzing data for ten cancer types, we found that a significantly higher proportion of the DEGs with absolute mRNA abundance changes overlapped or directly interacted with known cancer driver genes and anti-cancer drug targets than that of the DEGs only with mRNA concentration changes alone identified by the traditional methods. The DEGs with increased absolute mRNA abundances were enriched in DNA damage-related pathways, while DEGs with decreased absolute mRNA abundances were enriched in immune and metabolism associated pathways.

CONCLUSIONS

Both the mRNA concentration and the absolute abundances of the DEGs identified through REOs comparison change in disease samples in comparison with normal samples. In cancers these genes might play more important upstream roles in carcinogenesis.

A qualitative signature for early diagnosis of hepatocellular carcinoma based on relative expression orderings

BACKGROUND & AIMS

Currently, using biopsy specimens to confirm suspicious liver lesions of early hepatocellular carcinoma are not entirely reliable because of insufficient sampling amount and inaccurate sampling location. It is necessary to develop a signature to aid early hepatocellular carcinoma diagnosis using biopsy specimens even when the sampling location is inaccurate.

METHODS

Based on the within-sample relative expression orderings of gene pairs, we identified a simple qualitative signature to distinguish both hepatocellular carcinoma and adjacent non-tumour tissues from cirrhosis tissues of non-hepatocellular carcinoma patients.

RESULTS

A signature consisting of 19 gene pairs was identified in the training data sets and validated in 2 large collections of samples from biopsy and surgical resection specimens. For biopsy specimens, 95.7% of 141 hepatocellular carcinoma tissues and all (100%) of 108 cirrhosis tissues of non-hepatocellular carcinoma patients were correctly classified. Especially, all (100%) of 60 hepatocellular carcinoma adjacent normal tissues and 77.5% of 80 hepatocellular carcinoma adjacent cirrhosis tissues were classified to hepatocellular carcinoma. For surgical resection specimens, 99.7% of 733 hepatocellular carcinoma specimens were correctly classified to hepatocellular carcinoma, while 96.1% of 254 hepatocellular carcinoma adjacent cirrhosis tissues and 95.9% of 538 hepatocellular carcinoma adjacent normal tissues were classified to hepatocellular carcinoma. In contrast, 17.0% of 47 cirrhosis from non-hepatocellular carcinoma patients waiting for liver transplantation were classified to hepatocellular carcinoma, indicating that some patients with long-lasting cirrhosis could have already gained hepatocellular carcinoma characteristics.

CONCLUSIONS

The signature can distinguish both hepatocellular carcinoma tissues and tumour-adjacent tissues from cirrhosis tissues of non-hepatocellular carcinoma patients even using inaccurately sampled biopsy specimens, which can aid early diagnosis of hepatocellular carcinoma.