Integration of miRNA expression analysis of purified leukocytes and whole blood reveals blood-borne candidate biomarkers for lung cancer.

Changes in leukocyte populations may confound the disease-associated miRNA signals in the blood of cancer patients. We aimed to develop a method to detect differentially expressed miRNAs from lung cancer whole blood samples that are not influenced by variations in leukocyte proportions. The Ref-miREO method identifies differential miRNAs unaffected by changes in leukocyte populations by comparing the within-sample relative expression orderings (REOs) of miRNAs from healthy leukocyte subtypes and those from lung cancer blood samples. Over 77% of the differential miRNAs observed between lung cancer and healthy blood samples overlapped with those between myeloid-derived and lymphoid-derived leukocytes, suggesting the potential impact of changes in leukocyte populations on miRNA profile. Ref-miREO identified 16 differential miRNAs that target 19 lung adenocarcinoma-related genes previously linked to leukocytes. These miRNAs showed enrichment in cancer-related pathways and demonstrated high potential as diagnostic biomarkers, with the LASSO regression models effectively distinguishing between healthy and lung cancer blood or serum samples (all AUC > 0.85). Additionally, 12 of these miRNAs exhibited significant prognostic correlations. The Ref-miREO method offers valuable candidates for circulating biomarker detection in cancer that are not affected by changes in leukocyte populations.

Integrated analysis of diverse cancer types reveals a breast cancer-specific serum miRNA biomarker through relative expression orderings analysis.

PURPOSE: Serum microRNA (miRNA) holds great potential as a non-invasive biomarker for diagnosing breast cancer (BrC). However, most diagnostic models rely on the absolute expression levels of miRNAs, which are susceptible to batch effects and challenging for clinical transformation. Furthermore, current studies on liquid biopsy diagnostic biomarkers for BrC mainly focus on distinguishing BrC patients from healthy controls, needing more specificity assessment. METHODS: We collected a large number of miRNA expression data involving 8465 samples from GEO, including 13 different cancer types and non-cancer controls. Based on the relative expression orderings (REOs) of miRNAs within each sample, we applied the greedy, LASSO multiple linear regression, and random forest algorithms to identify a qualitative biomarker specific to BrC by comparing BrC samples to samples of other cancers as controls. RESULTS: We developed a BrC-specific biomarker called 7-miRPairs, consisting of seven miRNA pairs. It demonstrated comparable classification performance in our analyzed machine learning algorithms while requiring fewer miRNA pairs, accurately distinguishing BrC from 12 other cancer types. The diagnostic performance of 7-miRPairs was favorable in the training set (accuracy = 98.47%, specificity = 98.14%, sensitivity = 99.25%), and similar results were obtained in the test set (accuracy = 97.22%, specificity = 96.87%, sensitivity = 98.02%). KEGG pathway enrichment analysis of the 11 miRNAs within the 7-miRPairs revealed significant enrichment of target mRNAs in pathways associated with BrC. CONCLUSION: Our study provides evidence that utilizing serum miRNA pairs can offer significant advantages for BrC-specific diagnosis in clinical practice by directly comparing serum samples with BrC to other cancer types.

Identifying individualized prognostic signature and unraveling the molecular mechanism of recurrence in early-onset colorectal cancer.

BACKGROUND: The incidence and mortality of early-onset colorectal cancer (EOCRC; < 50 years old) is increasing worldwide, with a high recurrence rate. The inherent heterogeneity of EOCRC makes its treatment challenging. Hence, to further understand the biology and reveal the molecular mechanisms of EOCRC, a recurrence risk signature is needed to guide clinical management. METHODS: Based on the relative expression orderings (REOs) of genes in each sample, a prognostic signature was developed and validated utilizing multiple independent datasets. The underlying molecular mechanisms between distinct prognostic groups were explored via integrative analysis of multi-omics data. RESULTS: The prognostic signature consisting of 6 gene pairs (6-GPS) could predict the recurrence risk for EOCRC at the individual level. High-risk EOCRC classified by 6-GPS showed a poor prognosis but a good response to adjuvant chemotherapy. Moreover, high-risk EOCRC was characterized by epithelial-mesenchymal transition (EMT) and enriched angiogenesis, and had higher mutation burden, immune cell infiltration, and PD-1/PD-L1 expression. Furthermore, we identified four genes associated with relapse-free survival in EOCRC, including SERPINE1, PECAM1, CDH1, and ANXA1. They were consistently differentially expressed at the transcriptome and proteome levels between high-risk and low-risk EOCRCs. They were also involved in regulating cancer progression and immune microenvironment in EOCRC. Notably, the expression of SERPINE1 and ANXA1 positively correlated with M2-like macrophage infiltration. CONCLUSION: Our results indicate that 6-GPS can robustly predict the recurrence risk of EOCRC, and that SERPINE1, PECAM1, CDH1, and ANXA1 may serve as potential therapeutic targets. This study provides valuable information for the precision treatment of EOCRC.

Screening prognostic markers for hepatocellular carcinoma based on pyroptosis-related lncRNA pairs.

BACKGROUND: Pyroptosis is closely related to cancer prognosis. In this study, we tried to construct an individualized prognostic risk model for hepatocellular carcinoma (HCC) based on within-sample relative expression orderings (REOs) of pyroptosis-related lncRNAs (PRlncRNAs). METHODS: RNA-seq data of 343 HCC samples derived from The Cancer Genome Atlas (TCGA) database were analyzed. PRlncRNAs were detected based on differentially expressed lncRNAs between sample groups clustered by 40 reported pyroptosis-related genes (PRGs). Univariate Cox regression was used to screen out prognosis-related PRlncRNA pairs. Then, based on REOs of prognosis-related PRlncRNA pairs, a risk model for HCC was constructed by combining LASSO and stepwise multivariate Cox regression analysis. Finally, a prognosis-related competing endogenous RNA (ceRNA) network was built based on information about lncRNA-miRNA-mRNA interactions derived from the miRNet and TargetScan databases. RESULTS: Hierarchical clustering of HCC patients according to the 40 PRGs identified two groups with a significant survival difference (Kaplan-Meier log-rank, p = 0.026). Between the two groups, 104 differentially expressed lncRNAs were identified (|log2(FC)|> 1 and FDR < 5%). Among them, 83 PRlncRNA pairs showed significant associations between their REOs within HCC samples and overall survival (Univariate Cox regression, p < 0.005). An optimal 11-PRlncRNA-pair prognostic risk model was constructed for HCC. The areas under the curves (AUCs) of time-dependent receiver operating characteristic (ROC) curves of the risk model for 1-, 3-, and 5-year survival were 0.737, 0.705, and 0.797 in the validation set, respectively. Gene Set Enrichment Analysis showed that inflammation-related interleukin signaling pathways were upregulated in the predicted high-risk group (p < 0.05). Tumor immune infiltration analysis revealed a higher abundance of regulatory T cells (Tregs) and M2 macrophages and a lower abundance of CD8 + T cells in the high-risk group, indicating that excessive pyroptosis might occur in high-risk patients. Finally, eleven lncRNA-miRNA-mRNA regulatory axes associated with pyroptosis were established. CONCLUSION: Our risk model allowed us to determine the robustness of the REO-based PRlncRNA prognostic biomarkers in the stratification of HCC patients at high and low risk. The model is also helpful for understanding the molecular mechanisms between pyroptosis and HCC prognosis. High-risk patients may have excessive pyroptosis and thus be less sensitive to immune therapy.

Construct of qualitative diagnostic biomarkers specific for glioma by pairing serum microRNAs.

BACKGROUND: Serum microRNAs (miRNAs) are promising non-invasive biomarkers for diagnosing glioma. However, most reported predictive models are constructed without a large enough sample size, and quantitative expression levels of their constituent serum miRNAs are susceptible to batch effects, decreasing their clinical applicability. METHODS: We propose a general method for detecting qualitative serum predictive biomarkers using a large cohort of miRNA-profiled serum samples (n = 15,460) based on the within-sample relative expression orderings of miRNAs. RESULTS: Two panels of miRNA pairs (miRPairs) were developed. The first was composed of five serum miRPairs (5-miRPairs), reaching 100% diagnostic accuracy in three validation sets for distinguishing glioma and non-cancer controls (n = 436: glioma = 236, non-cancers = 200). An additional validation set without glioma samples (non-cancers = 2611) showed a predictive accuracy of 95.9%. The second panel included 32 serum miRPairs (32-miRPairs), reaching 100% diagnostic performance in training set on specifically discriminating glioma from other cancer types (sensitivity = 100%, specificity = 100%, accuracy = 100%), which was reproducible in five validation datasets (n = 3387: glioma = 236, non-glioma cancers = 3151, sensitivity> 97.9%, specificity> 99.5%, accuracy> 95.7%). In other brain diseases, the 5-miRPairs classified all non-neoplastic samples as non-cancer, including stroke (n = 165), Alzheimer's disease (n = 973), and healthy samples (n = 1820), and all neoplastic samples as cancer, including meningioma (n = 16), and primary central nervous system lymphoma samples (n = 39). The 32-miRPairs predicted 82.2 and 92.3% of the two kinds of neoplastic samples as positive, respectively. Based on the Human miRNA tissue atlas database, the glioma-specific 32-miRPairs were significantly enriched in the spinal cord (p = 0.013) and brain (p = 0.015). CONCLUSIONS: The identified 5-miRPairs and 32-miRPairs provide potential population screening and cancer-specific biomarkers for glioma clinical practice.

Prognostic value of a microRNA-pair signature in laryngeal squamous cell carcinoma patients.

PURPOSE: Predicting the prognosis in laryngeal squamous cell carcinoma (LSCC) patients will improve clinical decision-making. Here, we aimed to identify a qualitative signature based on the within-sample relative expression orderings (REOs) of microRNA (miRNA) pairs to predict the overall survival (OS) of LSCC patients. METHODS: First, we constructed non-repeating miRNA pairs based on differentially expressed miRNAs (DEmiRNAs) between LSCC and normal tissues. Then, we applied a bootstrap-based feature selection method to identify a robust miRNA-pair signature. The bootstrap-based feature selection improved the stability of feature selection by an ensemble based on the data perturbation. Furthermore, a series of bioinformatics analyses were carried out to explore the potential mechanisms of the signature and potential drug targets for LSCC. RESULTS: Based on the REOs of miRNA pairs, we identified a qualitative signature that consisted of 12 miRNA pairs. The constructed signature has good performance in predicting the OS of LSCC patients. It is robust against batch effects and more suitable for individual clinical applications. Furthermore, we identified several hub genes that may be potential drug targets for LSCC. CONCLUSION: Overall, our findings provided a promising signature for predicting the OS of LSCC patients.

Identification of a Qualitative Signature for the Diagnosis of Dementia With Lewy Bodies.

Background and purpose: Diagnosis of dementia with Lewy bodies (DLB) is highly challenging, primarily due to a lack of valid and reliable diagnostic tools. To date, there is no report of qualitative signature for the diagnosis of DLB. We aimed to develop a blood-based qualitative signature for differentiating DLB patients from healthy controls. Methods: The GSE120584 dataset was downloaded from the public database Gene Expression Omnibus (GEO). We combined multiple methods to select features based on the within-sample relative expression orderings (REOs) of microRNA (miRNA) pairs. Specifically, we first quickly selected miRNA pairs related to DLB by identifying reversal stable miRNA pairs. Then, an optimal miRNA pair subset was extracted by random forest (RF) and support vector machine-recursive feature elimination (SVM-RFE) methods. Furthermore, we applied logistic regression (LR) and SVM to build several prediction models. The model performance was assessed using the receiver operating characteristic curve (ROC) analysis. Lastly, we conducted bioinformatics analyses to explore the molecular mechanisms of the discovered miRNAs. Results: A qualitative signature consisted of 17 miRNA pairs and two clinical factors was identified for discriminating DLB patients from healthy controls. The signature is robust against experimental batch effects and applicable at the individual levels. The accuracies of the-signature-based models on the test set are 82.61 and 79.35%, respectively, indicating that the signature has acceptable discrimination performance. Moreover, bioinformatics analyses revealed that predicted target genes were enriched in 11 Go terms and 2 KEGG pathways. Moreover, five potential hub genes were found for DLB, including SRF, MAPK1, YWHAE, RPS6KA3, and KDM7A. Conclusion: This study provided a blood-based qualitative signature with the potential to be used as an effective tool to improve the accuracy of DLB diagnosis.

An absolute human stemness index associated with oncogenic dedifferentiation.

The progression of cancer is accompanied by the acquisition of stemness features. Many stemness evaluation methods based on transcriptional profiles have been presented to reveal the relationship between stemness and cancer. However, instead of absolute stemness index values-the values with certain range-these methods gave the values without range, which makes them unable to intuitively evaluate the stemness. Besides, these indices were based on the absolute expression values of genes, which were found to be seriously influenced by batch effects and the composition of samples in the dataset. Recently, we have showed that the signatures based on the relative expression orderings (REOs) of gene pairs within a sample were highly robust against these factors, which makes that the REO-based signatures have been stably applied in the evaluations of the continuous scores with certain range. Here, we provided an absolute REO-based stemness index to evaluate the stemness. We found that this stemness index had higher correlation with the culture time of the differentiated stem cells than the previous stemness index. When applied to the cancer and normal tissue samples, the stemness index showed its significant difference between cancers and normal tissues and its ability to reveal the intratumor heterogeneity at stemness level. Importantly, higher stemness index was associated with poorer prognosis and greater oncogenic dedifferentiation reflected by histological grade. All results showed the capability of the REO-based stemness index to assist the assignment of tumor grade and its potential therapeutic and diagnostic implications.

A Qualitative Transcriptional Signature for Predicting CpG Island Methylator Phenotype Status of the Right-Sided Colon Cancer.

A part of colorectal cancer which is characterized by simultaneous numerous hypermethylation CpG islands sites is defined as CpG island methylator phenotype (CIMP) status. Stage II and III CIMP-positive (CIMP+) right-sided colon cancer (RCC) patients have a better prognosis than CIMP-negative (CIMP-) RCC treated with surgery alone. However, there is no gold standard available in defining CIMP status. In this work, we selected the gene pairs whose relative expression orderings (REOs) were associated with the CIMP status, to develop a qualitative transcriptional signature to individually predict CIMP status for stage II and III RCC. Based on the REOs of gene pairs, a signature composed of 19 gene pairs was developed to predict the CIMP status of RCC through a feature selection process. A sample is predicted as CIMP+ when the gene expression orderings of at least 12 gene pairs vote for CIMP+; otherwise the CIMP-. The difference of prognosis between the predicted CIMP+ and CIMP- groups was more significantly different than the original CIMP status groups. There were more differential methylation and expression characteristics between the two predicted groups. The hierarchical clustering analysis showed that the signature could perform better for predicting CIMP status of RCC than current methods. In conclusion, the qualitative transcriptional signature for classifying CIMP status at the individualized level can predict outcome and guide therapy for RCC patients.

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.

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.

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.

Individualized predictive signatures for 5-fluorouracil-based chemotherapy in right- and left-sided colon cancer.

5-Fluorouracil (5-FU)-based adjuvant chemotherapy (ACT) is widely used for the treatment of colon cancer. Colon cancers with different primary tumor locations are clinically and molecularly distinct, implied through their response to 5-FU-based ACT. In this work, using 69 and 133 samples of patients with stage II-III right-sided and left-sided colon cancer (RCC and LCC) treated with post-surgery 5-FU-based ACT, we preselected gene pairs whose relative expression orderings were significantly correlated with the disease-free survival of patients by univariate Cox proportional hazards model. Then, from the identified prognostic-related gene pairs, a forward-stepwise selection algorithm was formulated to search for an optimal subset of gene pairs that resulted in the highest concordance index, referred to as the gene pair signature (GPS). We identified prognostic signatures, 3-GPS and 5-GPS, for predicting response to 5-FU-based ACT of patients with RCC and LCC, respectively, which were validated in independent datasets of GSE14333 and GSE72970. With the aid of the signatures, the transcriptional and genomic characteristics between the predicted responders and non-responders were explored. Notably, both in RCC and LCC, the predicted responders to 5-FU-based ACT were characterized by hypermutation, whereas the predicted non-responders were characterized by frequent copy number alternations. Finally, in comparison with the established relative expression ordering-based signature, which was developed without considering the differences between RCC and LCC, the newly proposed signatures had a better predictive performance. In conclusion, 3-GPS or 5-GPS can robustly predict response to 5-FU-based ACT for patients with RCC or LCC, respectively, in an individual level.

Tamoxifen therapy benefit predictive signature coupled with prognostic signature of post-operative recurrent risk for early stage ER+ breast cancer.

Two types of prognostic signatures for predicting recurrent risk of ER+ breast cancer patients have been developed: one type for patients accepting surgery only and another type for patients receiving post-operative tamoxifen therapy. However, the first type of signature cannot distinguish high-risk patients who cannot benefit from tamoxifen therapy, while the second type of signature cannot identify patients who will be at low risk of recurrence even if they accept surgery only. In this study, we proposed to develop two coupled signatures to solve these problems based on within-sample relative expression orderings (REOs) of gene pairs. Firstly, we identified a prognostic signature of post-operative recurrent risk using 544 samples of ER+ breast cancer patients accepting surgery only. Then, applying this drug-free signature to 840 samples of patients receiving post-operative tamoxifen therapy, we recognized 553 samples of patients who would have been at high risk of recurrence if they had accepted surgery only and used these samples to develop a tamoxifen therapy benefit predictive signature. The two coupled signatures were validated in independent data. The signatures developed in this study are robust against experimental batch effects and applicable at the individual levels, which can facilitate the clinical decision of tamoxifen therapy.