An individualized gene expression signature for prediction of lung adenocarcinoma metastases

The multiple roles of interferon regulatory factor family in health and disease

Interferon Regulatory Factors (IRFs), a family of transcription factors, profoundly influence the immune system, impacting both physiological and pathological processes. This review explores the diverse functions of nine mammalian IRF members, each featuring conserved domains essential for interactions with other transcription factors and cofactors. These interactions allow IRFs to modulate a broad spectrum of physiological processes, encompassing host defense, immune response, and cell development. Conversely, their pivotal role in immune regulation implicates them in the pathophysiology of various diseases, such as infectious diseases, autoimmune disorders, metabolic diseases, and cancers. In this context, IRFs display a dichotomous nature, functioning as both tumor suppressors and promoters, contingent upon the specific disease milieu. Post-translational modifications of IRFs, including phosphorylation and ubiquitination, play a crucial role in modulating their function, stability, and activation. As prospective biomarkers and therapeutic targets, IRFs present promising opportunities for disease intervention. Further research is needed to elucidate the precise mechanisms governing IRF regulation, potentially pioneering innovative therapeutic strategies, particularly in cancer treatment, where the equilibrium of IRF activities is of paramount importance.

An individualized immune prognostic signature in nasopharyngeal carcinoma

BACKGROUND

Immune-related characteristics can serve as reliable prognostic biomarkers in various cancers. Herein, we aimed to construct an individualized immune prognostic signature in nasopharyngeal carcinoma (NPC).

METHODS

This study retrospectively included 455 NPC samples and 39 normal healthy nasopharyngeal tissue specimens. Samples from Gene Expression Omnibus (GEO) were obtained as discovery cohort to screen candidate prognostic immune-related gene pairs based on relative expression ordering of the genes. Quantitative real-time reverse transcription-PCR was used to detect the selected genes to construct an immune-related gene pair signature in training cohort, which comprised 118 clinical samples, and was then validated in validation cohort 1, comprising 92 clinical samples, and validation cohort 2, comprising 88 samples from GEO.

RESULTS

We identified 26 immune-related gene pairs as prognostic candidates in discovery cohort. A prognostic immune signature comprising 11 immune gene pairs was constructed in training cohort. In validation cohort 1, the immune signature could significantly distinguish patients with high or low risk in terms of progression-free survival (PFS) (hazard ratio [HR] 2.66, 95 % confidence interval (CI) 1.17-6.02, P=0.015) and could serve as an independent prognostic factor for PFS in multivariate analysis (HR 2.66, 95 % CI 1.17-6.02, P=0.019). Similar results were obtained using validation cohort 2, in which PFS was significantly worse in high risk group than in low risk group (HR 3.02, 95 % CI 1.12-8.18, P=0.022).

CONCLUSIONS

The constructed immune signature showed promise for estimating prognosis in NPC. It has potential for translation into clinical practice after prospective validation.

Generative pretraining from large-scale transcriptomes for single-cell deciphering

Exponential accumulation of single-cell transcriptomes poses great challenge for efficient assimilation. Here, we present an approach entitled generative pretraining from transcriptomes (tGPT) for learning feature representation of transcriptomes. tGPT is conceptually simple in that it autoregressive models the ranking of a gene in the context of its preceding neighbors. We developed tGPT with 22.3 million single-cell transcriptomes and used four single-cell datasets to evalutate its performance on single-cell analysis tasks. In addition, we examine its applications on bulk tissues. The single-cell clusters and cell lineage trajectories derived from tGPT are highly aligned with known cell labels and states. The feature patterns of tumor bulk tissues learned by tGPT are associated with a wide range of genomic alteration events, prognosis, and treatment outcome of immunotherapy. tGPT represents a new analytical paradigm for integrating and deciphering massive amounts of transcriptome data and it will facilitate the interpretation and clinical translation of single-cell transcriptomes.

Sec61γ is a vital protein in the endoplasmic reticulum membrane promoting tumor metastasis and invasion in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is one of the most common malignant tumors worldwide. Finding effective prognostic markers and therapeutic targets is of great significance for controlling metastasis and invasion clinically.

METHODS

The open copy-number aberrations and gene expression datasets were analysed, and the data of 102 LUAD patients was used for further validation. The cell proliferation, colony formation, migration, invasion assays and mice tumor models were used to detect the function of SEC61G. The epidermal growth factor receptor (EGFR) pathway was also detected to find the mechanism of Sec61γ.

RESULTS

Based on the open datasets, we found that the high level of SEC61G mRNA may drive LUAD metastasis. Furthermore, the overexpression of Sec61γ protein was significantly associated with poor prognosis and greater tumor cell proliferation and metastasis. The SEC61G knockdown could inhibit the EGFR pathway, including STAT3, AKT and PI3K, which can be reversed by Sec61γ overexpression and epithelial growth factor (EGF) supplement.

CONCLUSIONS

Sec61γ promoted the proliferation, metastasis, and invasion of LUAD through EGFR pathways. Sec61γ might be a potential target for the treatment of LUAD metastases.

Artificial intelligence-based non-small cell lung cancer transcriptome RNA-sequence analysis technology selection guide

The incidence and mortality rates of lung cancer are high worldwide, where non-small cell lung cancer (NSCLC) accounts for more than 85% of lung cancer cases. Recent non-small cell lung cancer research has been focused on analyzing patient prognosis after surgery and identifying mechanisms in connection with clinical cohort and ribonucleic acid (RNA) sequencing data, including single-cell ribonucleic acid (scRNA) sequencing data. This paper investigates statistical techniques and artificial intelligence (AI) based non-small cell lung cancer transcriptome data analysis methods divided into target and analysis technology groups. The methodologies of transcriptome data were schematically categorized so researchers can easily match analysis methods according to their goals. The most widely known and frequently utilized transcriptome analysis goal is to find essential biomarkers and classify carcinomas and cluster NSCLC subtypes. Transcriptome analysis methods are divided into three major categories: Statistical analysis, machine learning, and deep learning. Specific models and ensemble techniques typically used in NSCLC analysis are summarized in this paper, with the intent to lay a foundation for advanced research by converging and linking the various analysis methods available.

Extracellular matrix profiles determine risk and prognosis of the squamous cell carcinoma subtype of non-small cell lung carcinoma

BACKGROUND

Squamous cell carcinoma (SqCC) is a subtype of non-small cell lung cancer for which patient prognosis remains poor. The extracellular matrix (ECM) is critical in regulating cell behavior; however, its importance in tumor aggressiveness remains to be comprehensively characterized.

METHODS

Multi-omics data of SqCC human tumor specimens was combined to characterize ECM features associated with initiation and recurrence. Penalized logistic regression was used to define a matrix risk signature for SqCC tumors and its performance across a panel of tumor types and in SqCC premalignant lesions was evaluated. Consensus clustering was used to define prognostic matreotypes for SqCC tumors. Matreotype-specific tumor biology was defined by integration of bulk RNAseq with scRNAseq data, cell type deconvolution, analysis of ligand-receptor interactions and enriched biological pathways, and through cross comparison of matreotype expression profiles with aging and idiopathic pulmonary fibrosis lung profiles.

RESULTS

This analysis revealed subtype-specific ECM signatures associated with tumor initiation that were predictive of premalignant progression. We identified an ECM-enriched tumor subtype associated with the poorest prognosis. In silico analysis indicates that matrix remodeling programs differentially activate intracellular signaling in tumor and stromal cells to reinforce matrix remodeling associated with resistance and progression. The matrix subtype with the poorest prognosis resembles ECM remodeling in idiopathic pulmonary fibrosis and may represent a field of cancerization associated with elevated cancer risk.

CONCLUSIONS

Collectively, this analysis defines matrix-driven features of poor prognosis to inform precision medicine prevention and treatment strategies towards improving SqCC patient outcome.

A novel epithelial-mesenchymal transition (EMT)-related gene signature of predictive value for the survival outcomes in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a remarkably heterogeneous and aggressive disease with dismal prognosis of patients. The identification of promising prognostic biomarkers might enable effective diagnosis and treatment of LUAD. Aberrant activation of epithelial-mesenchymal transition (EMT) is required for LUAD initiation, progression and metastasis. With the purpose of identifying a robust EMT-related gene signature (E-signature) to monitor the survival outcomes of LUAD patients. In The Cancer Genome Atlas (TCGA) database, least absolute shrinkage and selection operator (LASSO) analysis and cox regression analysis were conducted to acquire prognostic and EMT-related genes. A 4 EMT-related and prognostic gene signature, comprising dickkopf-like protein 1 (DKK1), lysyl oxidase-like 2 (LOXL2), matrix Gla protein (MGP) and slit guidance ligand 3 (SLIT3), was identified. By the usage of datum derived from TCGA database and Western blotting analysis, compared with adjacent tissue samples, DKK1 and LOXL2 protein expression in LUAD tissue samples were significantly higher, whereas the trend of MGP and SLIT3 expression were opposite. Concurrent with upregulation of epithelial markers and downregulation of mesenchymal markers, knockdown of DKK1 and LOXL2 impeded the migration and invasion of LUAD cells. Simultaneously, MGP and SLIT3 silencing promoted metastasis and induce EMT of LUAD cells. In the TCGA-LUAD set, receiver operating characteristic (ROC) analysis indicated that our risk model based on the identified E-signature was superior to those reported in literatures. Additionally, the E-signature carried robust prognostic significance. The validity of prediction in the E-signature was validated by the three independent datasets obtained from Gene Expression Omnibus (GEO) database. The probabilistic nomogram including the E-signature, pathological T stage and N stage was constructed and the nomogram demonstrated satisfactory discrimination and calibration. In LUAD patients, the E-signature risk score was associated with T stage, N stage, M stage and TNM stage. GSEA (gene set enrichment analysis) analysis indicated that the E-signature might be linked to the pathways including GLYCOLYSIS, MYC TARGETS, DNA REPAIR and so on. In conclusion, our study explored an innovative EMT based prognostic signature that might serve as a potential target for personalized and precision medicine.

A novel qualitative signature based on lncRNA pairs for prognosis prediction in hepatocellular carcinoma

Background

Prognostic assessment is imperative for clinical management of patients with hepatocellular carcinoma (HCC). Most reported prognostic signatures are based on risk scores summarized from quantitative expression level of candidate genes, which are vulnerable against experimental batch effects and impractical for clinical application. We aimed to develop a robust qualitative signature to assess individual survival risk for HCC patients.

Methods

Long non-coding RNA (lncRNA) pairs correlated with overall survival (OS) were identified and an optimal combination of lncRNA pairs based on the majority voting rule was selected as a classification signature to predict the overall survival risk in the cancer genome atlas (TCGA). Then, the signature was further validated in two external datasets. Besides, biomolecular characteristics, immune infiltration status, and chemotherapeutics efficacy of different risk groups were further compared. Finally, we performed key lncRNA screening and validated it in vitro.

Results

A signature consisting of 50 lncRNA pairs (50-LPS) was identified in TCGA and successfully validated in external datasets. Patients in the high-risk group, when at least 25 of the 50-LPS voted for high risk, had significantly worse OS than the low-risk group. Multivariate Cox, receiver operating characteristic (ROC) curve and decision curve analyses (DCA) demonstrated that the 50-LPS was an independent prognostic factor and more powerful than other available clinical factors in OS prediction. Comparison analyses indicated that different risk groups had distinct biomolecular characteristics, immune infiltration status, and chemotherapeutics efficacy. TDRKH-AS1 was confirmed as a key lncRNA and associated with cell growth of HCC.

Conclusions

The 50-LPS could not only predict the prognosis of HCC patients robustly and individually, but also provide theoretical basis for therapy. Besides, TDRKH-AS1 was identified as a key lncRNA in the proliferation of HCC. The 50-LPS might guide personalized therapy for HCC patients in clinical practice.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-022-02507-z.

Nonnegative matrix factorization-based bioinformatics analysis reveals that TPX2 and SELENBP1 are two predictors of the inner sub-consensuses of lung adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is a heterogeneous disease. However the inner sub‐groups of LUAD have not been fully studied. Markers predicted the sub‐groups and prognosis of LUAD are badly needed.

Aims

To identify biomarkers associated with the sub‐groups and prognosis of LUAD.

Materials and Methods

Using nonnegative matrix factorization (NMF) clustering, LUAD patients from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) datasets and LUAD cell lines from Genomics of Drug Sensitivity in Cancer (GDSC) dataset were divided into different sub‐consensuses based on the gene expression profiling. The overall survival of LUAD patients in each sub‐consensus was determined by Kaplan‐Meier survival analysis. The common genes which were differentially expressed in each sub‐consensus of LUAD patients and LUAD cell lines were identified using TBtools. The predictive accuracy of TPX2 and SELENBP1 for theinner sub‐consensuses of LUAD was determined by Receiver operator characteristic (ROC) analysis. The Kaplan‐Meier survival analysis was also used to test the prognostic significance of TPX2 and SELENBP1 in LUAD patients.

Results

Using nonnegative matrix factorization clustering, LUAD patients in The Cancer Genome Atlas (TCGA), GSE30219, GSE42127, GSE50081, GSE68465, and GSE72094 datasets were divided into three sub‐consensuses. Sub‐consensus3 LUAD patients were with low overall survival and were with high TP53 mutations. Similarly, LUAD cell lines were also divided into three sub‐consensuses by NMF method, and sub‐consensus2 cell lines were resistant to EGFR inhibitors. Identification of the common genes which were differentially expressed in different sub‐consensuses of LUAD patients and LUAD cell lines revealed that TPX2 was highly expressed in sub‐consensus3 LUAD patients and sub‐consensus2 LUAD cell lines. On the contrary, SELENBP1 was highly expressed in sub‐consensus1 LUAD patients and sub‐consensus1 LUAD cell lines. The expression levels of TPX2 and SELENBP1 could distinguish sub‐consensus3 LUAD patients or sub‐consensus2 LUAD cell lines from other sub‐consensuses of LUAD patients or cell lines. Moreover, compared with normal lung tissues, TPX2 was highly expressed, while, SELENBP1 was lowly expressed in LUAD tissues. Furthermore, the higher expression levels of TPX2 were associated with the lower relapse‐free survival and the lower overall survival of LUAD patients. While, the higher expression levels of SELENBP1 were associated with the higher relapse‐free survival and higher overall survival. At last, we showed that TP53 mutant LUAD patients were with higher TPX2 and lower SELENBP1 expressions.

Discussion

Both iCluster and NMF method are proved to be robust LUAD classification systems. However, the LUAD patients in different iclusters had no significant clinical overall survival, while, sub‐consensus3 LUAD patients from NMF classification were with lower overall survival than other sub‐consensuses.

Conclusions

By integrated analysis of 1765 LUAD patients and 64 LUAD cell lines, we showed that NMF was a robust inner sub‐consensuses classification method of LUAD. TPX2 and SELENBP1 were differentially expressed in different LUAD sub‐ consensuses, and predicted the inner sub‐consensuses of LUAD with high accuracy. TPX2 was an unfavorable prognostic biomarker of LUAD which was up‐regulated in LUAD tissues and associated with the low overall survival of LUAD. SELENBP1 was a favorable prognostic biomarker of LUAD which was down‐regulated in LUAD tissues and associated with the prolonged overall survival of LUAD.

IRF1 regulates the progression of colorectal cancer via interferon‑induced proteins

Radiation is one of the main methods for the treatment of colorectal cancer (CRC) before or after surgery. However, radiotherapy tolerance of patients with CRC is often a major concern. Interferon regulatory factor 1 (IRF1) is a member of the IRF family and is involved in the development of multiple diseases, including tumors. The present study investigated the role of IRF1 in the development and radiation sensitivity of CRC. Immunohistochemistry was performed to examine the expression levels of IRF1 in tissue samples from patients with CRC, as well as in nude mice. MTT, 5‑ethynyl‑20‑deoxyuridine, colony formation, cell cycle alteration and apoptosis assays were performed in CRC cell lines. Western blotting and immunofluorescence were used to detect the expression levels of a series of proteins. RNA sequencing was applied to identify genes whose expression was upregulated by IRF1 overexpression. Xenograft nude mouse models and hematoxylin and eosin staining were used to validate the present findings in vivo. It was revealed that the expression levels of IRF1 were significantly lower in CRC tissues than in adjacent tissues. IRF1 upregulation inhibited cell proliferation and colony formation, caused G1 cell arrest, promoted cell apoptosis, and enhanced the sensitivity of CRC cells to X‑ray irradiation. The role of IRF1 in promoting the radiosensitivity of CRC was further demonstrated in nude mice with CRC xenografts. In addition, RNA sequencing revealed that overexpression of IRF1 in CRC cells significantly increased the expression levels of interferon‑induced protein family members interferon α inducible protein 6, interferon induced transmembrane protein 1 and interferon induced protein 35 (fold change >2.0). In summary, the present study demonstrated that the upregulation of IRF1 inhibited the progression and promoted the radiosensitivity of CRC, likely by regulating interferon‑induced proteins.

A Cancer-Specific Qualitative Method for Estimating the Proportion of Tumor-Infiltrating Immune Cells

Tumor-infiltrating immune cells are important components in the tumor microenvironment (TME) and different types of these cells exert different effects on tumor development and progression; these effects depend upon the type of cancer involved. Several methods have been developed for estimating the proportion of immune cells using bulk transcriptome data. However, there is a distinct lack of methods that are capable of predicting the immune contexture in specific types of cancer. Furthermore, the existing methods are based on absolute gene expression and are susceptible to experimental batch effects, thus resulting in incomparability across different datasets. In this study, we considered two common neoplasms as examples (colorectal cancer [CRC] and melanoma) and introduced the Tumor-infiltrating Immune Cell Proportion Estimator (TICPE), a cancer-specific qualitative method for estimating the proportion of tumor-infiltrating immune cells. The TICPE was based on the relative expression orderings (REOs) of gene pairs within a sample and is notably insensitive to batch effects. Performance evaluation using public expression data with mRNA mixtures, single-cell RNA-Seq (scRNA-Seq) data, immunohistochemistry data, and simulated bulk RNA-seq samples, indicated that the TICPE can estimate the proportion of immune cells with levels of accuracy that are clearly superior to other methods. Furthermore, we showed that the TICPE could effectively detect prognostic signals in patients with tumors and changes in the fractions of immune cells during immunotherapy in melanoma. In conclusion, our work presented a unique novel method, TICPE, to estimate the proportion of immune cells in specific cancer types and explore the effect of the infiltration of immune cells on the efficacy of immunotherapy and the prognosis of cancer. The source code for TICPE is available at https://github.com/huitingxiao/TICPE.

An individualized transcriptional signature to predict the epithelial-mesenchymal transition based on relative expression ordering

The epithelial-mesenchymal transition (EMT) process is involved in cancer cell metastasis and immune system activation. Hence, identification of gene expression signatures capable of predicting the EMT status of cancer cells is essential for development of therapeutic strategies. However, quantitative identification of EMT markers is limited by batch effects, the platform used, or normalization methods. We hypothesized that a set of EMT-related relative expression orderings are highly stable in epithelial samples yet are reversed in mesenchymal samples. To test this hypothesis, we analyzed transcriptome data for ovarian cancer cohorts from publicly available databases, to develop a qualitative 16-gene pair signature (16-GPS) that effectively distinguishes the mesenchymal from epithelial phenotype. Our method was superior to previous quantitative methods in terms of classification accuracy and applicability to individualized patients without requiring data normalization. Patients with mesenchymal-like ovarian cancer showed poorer overall survival compared to patients with epithelial-like ovarian cancer. Additionally, EMT score was positively correlated with expression of immune checkpoint genes and metastasis. We, therefore, established a robust EMT 16-GPS that is independent of detection platform, batch effects and individual variations, and which represents a qualitative signature for investigating the EMT and providing insights into immunotherapy for ovarian cancer patients.

Reference genome and annotation updates lead to contradictory prognostic predictions in gene expression signatures: a case study of resected stage I lung adenocarcinoma

RNA-sequencing enables accurate and low-cost transcriptome-wide detection. However, expression estimates vary as reference genomes and gene annotations are updated, confounding existing expression-based prognostic signatures. Herein, prognostic 9-gene pair signature (GPS) was applied to 197 patients with stage I lung adenocarcinoma derived from previous and latest data from The Cancer Genome Atlas (TCGA) processed with different reference genomes and annotations. For 9-GPS, 6.6% of patients exhibited discordant risk classifications between the two TCGA versions. Similar results were observed for other prognostic signatures, including IRGPI, 15-gene and ORACLE. We found that conflicting annotations for gene length and overlap were the major cause of their discordant risk classification. Therefore, we constructed a prognostic 40-GPS based on stable genes across GENCODE v20-v30 and validated it using public data of 471 stage I samples (log-rank P < 0.0010). Risk classification was still stable in RNA-sequencing data processed with the newest GENCODE v32 versus GENCODE v20-v30. Specifically, 40-GPS could predict survival for 30 stage I samples with formalin-fixed paraffin-embedded tissues (log-rank P = 0.0177). In conclusion, this method overcomes the vulnerability of existing prognostic signatures due to reference genome and annotation updates. 40-GPS may offer individualized clinical applications due to its prognostic accuracy and classification stability.

A qualitative transcriptional signature to reclassify histological grade of ER-positive breast cancer patients

BACKGROUND

Histological grade (HG) is commonly adopted as a prognostic factor for ER-positive breast cancer patients. However, HG evaluation methods, such as the pathological Nottingham grading system, are highly subjective with only 50-85% inter-observer agreements. Specifically, the subjectivity in the pathological assignment of the intermediate grade (HG2) breast cancers, comprising of about half of breast cancer cases, results in uncertain disease outcomes prediction. Here, we developed a qualitative transcriptional signature, based on within-sample relative expression orderings (REOs) of gene pairs, to define HG1 and HG3 and reclassify pathologically-determined HG2 (denoted as pHG2) breast cancer patients.

RESULTS

From the gene pairs with significantly stable REOs in pathologically-determined HG1 (denoted as pHG1) samples and reversely stable REOs in pathologically-determined HG3 (denoted as pHG3) samples, concordantly identified from seven datasets, we extracted a signature which could determine the HG state of samples through evaluating whether the within-sample REOs match with the patterns of the pHG1 REOs or pHG3 REOs. A sample was classified into the HG3 group if at least a half of the REOs of the 10 gene pairs signature within this sample voted for HG3; otherwise, HG1. Using four datasets including samples of early stage (I-II) ER-positive breast cancer patients who accepted surgery only, we validated that this signature was able to reclassify pHG2 patients into HG1 and HG3 groups with significantly different survival time. For the original pHG1 and pHG3 patients, the signature could also more accurately and objectively stratify them into distinct prognostic groups. And the up-regulated and down down-regulated genes in HG1 compared with HG3 involved in cell proliferation and extracellular signal transduction pathways respectively. By comparing with existing signatures, 10-GPS was with prognostic significance and was more aligned with survival of patients especially for pHG2 samples.

CONCLUSIONS

The transcriptional qualitative signature can provide an objective assessment of HG states of ER-positive breast cancer patients, especially for reclassifying patients with pHG2, to assist decision making on clinical therapy.

A qualitative transcriptional signature for predicting the biochemical recurrence risk of prostate cancer patients after radical prostatectomy

BACKGROUND

The qualitative transcriptional characteristics, the within-sample relative expression orderings (REOs) of genes, are highly robust against batch effects and sample quality variations. Hence, we develop a qualitative transcriptional signature based on REOs to predict the biochemical recurrence risk of prostate cancer (PCa) patients after radical prostatectomy.

METHODS

Gene pairs with REOs significantly correlated with the biochemical recurrence-free survival (BFS) were identified from 131 PCa samples in the training data set. From these gene pairs, we selected a qualitative transcriptional signature based on the within-sample REOs of gene pairs which could predict the recurrence risk of PCa patients after radical prostatectomy.

RESULTS

A signature consisting of 74 gene pairs, named 74-GPS, was developed for predicting the recurrence risk of PCa patients after radical prostatectomy based on the majority voting rule that a sample was assigned as high risk when at least 37 gene pairs of the 74-GPS voted for high risk; otherwise, low risk. The signature was validated in six independent datasets produced by different platforms. In each of the validation datasets, the Kaplan-Meier survival analysis showed that the average BFS of the low-risk group was significantly better than that of the high-risk group. Analyses of multiomics data of PCa samples from TCGA suggested that both the epigenomic and genomic alternations could cause the reproducible transcriptional differences between the two different prognostic groups.

CONCLUSIONS

The proposed qualitative transcriptional signature can robustly stratify PCa patients after radical prostatectomy into two groups with different recurrence risk and distinct multiomics characteristics. Hence, 74-GPS may serve as a helpful tool for guiding the management of PCa patients with radical prostatectomy at the individual level.

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.

Unique metastasis-associated lncRNA signature optimizes prediction of tumor relapse in lung adenocarcinoma

Background

Local relapses and metastases are primary causes of death in lung cancer patients. In the present study, we aimed to develop a prognostic signature based on metastasis‐associated lncRNAs in patients with lung adenocarcinoma (LUAD).

Methods

Firstly, the potential metastasis‐associated lncRNAs were identified by analyzing high‐throughput data from The Cancer Genome Atlas (TCGA), and based on which, an lncRNA signature was constructed for prediction of relapse in LUAD patients using Cox proportional hazards regression analysis. Moreover, the prognostic performance of the lncRNA signature was evaluated using Kaplan‐Meier survival analysis, time‐dependent receiver operating characteristic (ROC) curve and Cox analysis, respectively. In addition, the potential metastasis‐associated function of these six lncRNAs was confirmed by lncRNA over‐expression or depletion and in vitro transwell assays in LUAD cells.

Results

An lncRNA signature consisting of six most important prognostic factors (LINC01819, ZNF649‐AS1, HNF4A‐AS1, FAM222A‐AS1, LINC02323 and LINC00672) was developed. The signature was an independent predictor for patients' relapse‐free survival (RFS), which could provide higher tumor relapse prediction capability compared with the TNM staging system at three years and five years, respectively (P = 0.0209 and P = 0.0468). Furthermore, the combination of this lncRNA signature and TNM stage had better prognostic value than TNM stage alone at three and five years, respectively (P = 0.0006 and P = 0.0096). Additionally, all the lncRNAs of the signature had a regulatory role in the LUAD cell mobility.

Conclusions

This novel six‐lncRNA signature had considerable prognostic value for prediction of relapse in LUAD patients.

Key points

Identification of a small mutation panel of coding sequences to predict the efficacy of immunotherapy for lung adenocarcinoma

Background

Immune checkpoint inhibitors are effective in some cases of lung adenocarcinoma (LUAD). Whole-exome sequencing has revealed that the tumour mutation burden (TMB) is associated with clinical benefits among patients from immune checkpoint inhibitors. Several commercial mutation panels have been developed for estimating the TMB regardless of the cancer type. However, different cancer types have different mutational landscapes; hence, this study aimed to develop a small cancer-type-specific mutation panel for high-accuracy estimation of the TMB of LUAD patients.

Methods

We developed a small cancer-type-specific mutation panel based on coding sequences (CDSs) rather than genes, for LUAD patients. Using somatic CDSs mutation data from 486 LUAD patients in The Cancer Genome Atlas (TCGA) database, we pre-selected a set of CDSs with mutation states significantly correlated with the TMB, from which we selected a CDS mutation panel with a panel-score most significantly correlated with the TMB, using a genetic algorithm.

Results

A mutation panel containing 106 CDSs of 100 genes with only 0.34 Mb was developed, whose length was much shorter than current commercial mutation panels of 0.80–0.92 Mb. The correlation of this panel with the TMB was validated in two independent LUAD datasets with progression-free survival data for patients treated with nivolumab plus ipilimumab and pembrolizumab immunotherapies, respectively. In both test datasets, survival analyses revealed that patients with a high TMB predicted via the 106-CDS mutation panel with a cut-point of 6.20 mutations per megabase, median panel score in the training dataset, had a significantly longer progression-free survival than those with a low predicted TMB (log-rank p = 0.0018, HR = 3.35, 95% CI 1.51–7.42; log-rank p = 0.0020, HR = 5.06, 95% CI 1.63–15.69). This small panel better predicted the efficacy of immunotherapy than current commercial mutation panels.

Conclusions

The small-CDS mutation panel of only 0.34 Mb is superior to current commercial mutation panels and can better predict the efficacy of immunotherapy for LUAD patients, and its low cost and time-intensiveness make it more suitable for clinical applications.

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.

Trans-omics biomarker model improves prognostic prediction accuracy for early-stage lung adenocarcinoma

Limited studies have focused on developing prognostic models with trans-omics biomarkers for early-stage lung adenocarcinoma (LUAD). We performed integrative analysis of clinical information, DNA methylation, and gene expression data using 825 early-stage LUAD patients from 5 cohorts. Ranger algorithm was used to screen prognosis-associated biomarkers, which were confirmed with a validation phase. Clinical and biomarker information was fused using an iCluster plus algorithm, which significantly distinguished patients into high- and low-mortality risk groups (Pdiscovery = 0.01 and Pvalidation = 2.71×10-3). Further, potential functional DNA methylation-gene expression-overall survival pathways were evaluated by causal mediation analysis. The effect of DNA methylation level on LUAD survival was significantly mediated through gene expression level. By adding DNA methylation and gene expression biomarkers to a model of only clinical data, the AUCs of the trans-omics model improved by 18.3% (to 87.2%) and 16.4% (to 85.3%) in discovery and validation phases, respectively. Further, concordance index of the nomogram was 0.81 and 0.77 in discovery and validation phases, respectively. Based on systematic review of published literatures, our model was superior to all existing models for early-stage LUAD. In summary, our trans-omics model may help physicians accurately identify patients with high mortality risk.

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.

A seven-long noncoding RNA signature predicts relapse in patients with early-stage lung adenocarcinoma

BACKGROUND

Long noncoding RNA (lncRNA) has been increasingly reported to play crucial roles in cancer development. In this study, we aim to develop a lncRNA-based signature to predict the relapse of early-stage (stage I-II) lung adenocarcinoma (LUAD).

METHODS

With a lncRNA-mining strategy, lncRNA expression profiles of three LUAD cohorts were obtained from the Gene Expression Omnibus database. A risk score model was established based on the lncRNAs expression from training set (GSE31210, n = 204) and further validated in two independent testing sets (GSE50081, n = 124; and GSE30219, n = 84). The potential signaling pathways modulated by the prognostic lncRNAs were explored using bioinformatics analysis.

RESULTS

In the training set, seven lncRNAs were identified to be significantly correlated with the relapse-free survival (RFS) of early-stage LUAD, and were then aggregated to form a seven-lncRNA prognostic signature to classify patients into high-risk and low-risk groups. Individuals of training set in the high-risk group exhibited a poorer RFS than those in the low-risk group (HR: 7.574, 95% CI: 4.165-13.775; P < 0.001). The similar prognostic powers of the seven-lncRNA signature were also achieved in the two independent testing sets and in stratified analysis. Multivariate Cox regression indicated that the prognostic value of seven-lncRNA signature was independent of other clinical features. Functional enrichment analysis found that the seven-lncRNA signature may be involved in biological pathways such as cell cycle, DNA replication, and p53 signaling pathway.

CONCLUSION

Our results indicate that the seven-lncRNA signature may be an innovative biomarker to predict the relapse of early-stage LUAD.

Shared liver-like transcriptional characteristics in liver metastases and corresponding primary colorectal tumors

Background & Aims: Primary tumors of colorectal carcinoma (CRC) with liver metastasis might gain some liver-specific characteristics to adapt the liver micro-environment. This study aims to reveal potential liver-like transcriptional characteristics associated with the liver metastasis in primary colorectal carcinoma.

Methods: Among the genes up-regulated in normal liver tissues versus normal colorectal tissues, we identified “liver-specific” genes whose expression levels ranked among the bottom 10% (“unexpressed”) of all measured genes in both normal colorectal tissues and primary colorectal tumors without metastasis. These liver-specific genes were investigated for their expressions in both the primary tumors and the corresponding liver metastases of seven primary CRC patients with liver metastasis using microdissected samples.

Results: Among the 3958 genes detected to be up-regulated in normal liver tissues versus normal colorectal tissues, we identified 12 liver-specific genes and found two of them, ANGPTL3 and CFHR5, were unexpressed in microdissected primary colorectal tumors without metastasis but expressed in both microdissected liver metastases and corresponding primary colorectal tumors (Fisher's exact test, P < 0.05). Genes co-expressed with ANGPTL3 and CFHR5 were significantly enriched in metabolism pathways characterizing liver tissues, including “starch and sucrose metabolism” and “drug metabolism-cytochrome P450”.

Conclusions: For primary CRC with liver metastasis, both the liver metastases and corresponding primary colorectal tumors may express some liver-specific genes which may help the tumor cells adapt the liver micro-environment.