Computational identification of DNA damage-relevant lncRNAs for predicting therapeutic efficacy and clinical outcomes in cancer

OBJECTIVES

The role of long non-coding RNAs (lncRNAs) in cancer treatment, particularly in modulating DNA repair programs, is an emerging field that warrants systematic exploration. This study aimed to systematically identify the lncRNA regulators that potentially regulate DNA damage response (DDR).

METHODS

Using genome-wide mRNA and lncRNA expression profiles of the same tumor patients, we proposed a novel computational framework. This framework performed Gene Set Variation Analysis to calculate DDR pathway enrichment score, which relies on weighting by tumor purity to obtain DDR activity score for each patient. Then, an in-depth differential expression profiling was conducted to identify pathway activity lncRNAs between high- and low-activity groups, utilizing a bootstrap-based randomization method.

RESULTS

We comprehensively charted the landscape of DDR-relevant lncRNAs across 23 epithelial-based cancer types. Its effectiveness was validated by assessing the consistency of these lncRNAs within various datasets of the same cancer type (hypergeometric test P < 0.001), examining the expression perturbation of these lncRNAs in response to treatment and demonstrating its application in prioritizing clinically-related lncRNAs. Furthermore, leveraging 820 epithelial ovarian cancer patients from four public datasets, we applied these lncRNAs identified by DDRLnc to establish and validate a risk stratification model to evaluate the benefits of platinum-based adjuvant chemotherapy for the improvement of clinical treatment outcomes.

CONCLUSIONS

Comprehensive pan-cancer analysis illustrates the utility of computational framework in identifying potentially lncRNAs involved in DDR, thereby offering novel insights into the complex realm of cancer research, and it will become a valuable tool for identifying potential therapeutic targets for cancer.

Predictive value of DNA methylation patterns in AML patients treated with an azacytidine containing induction regimen

BACKGROUND

Acute myeloid leukemia (AML) is a heterogeneous disease with a poor prognosis. Dysregulation of the epigenetic machinery is a significant contributor to disease development. Some AML patients benefit from treatment with hypomethylating agents (HMAs), but no predictive biomarkers for therapy response exist. Here, we investigated whether unbiased genome-wide assessment of pre-treatment DNA-methylation profiles in AML bone marrow blasts can help to identify patients who will achieve a remission after an azacytidine-containing induction regimen.

RESULTS

A total of n = 155 patients with newly diagnosed AML treated in the AMLSG 12-09 trial were randomly assigned to a screening and a refinement and validation cohort. The cohorts were divided according to azacytidine-containing induction regimens and response status. Methylation status was assessed for 664,227 500-bp-regions using methyl-CpG immunoprecipitation-seq, resulting in 1755 differentially methylated regions (DMRs). Top regions were distilled and included genes such as WNT10A and GATA3. 80% of regions identified as a hit were represented on HumanMethlyation 450k Bead Chips. Quantitative methylation analysis confirmed 90% of these regions (36 of 40 DMRs). A classifier was trained using penalized logistic regression and fivefold cross validation containing 17 CpGs. Validation based on mass spectra generated by MALDI-TOF failed (AUC 0.59). However, discriminative ability was maintained by adding neighboring CpGs. A recomposed classifier with 12 CpGs resulted in an AUC of 0.77. When evaluated in the non-azacytidine containing group, the AUC was 0.76.

CONCLUSIONS

Our analysis evaluated the value of a whole genome methyl-CpG screening assay for the identification of informative methylation changes. We also compared the informative content and discriminatory power of regions and single CpGs for predicting response to therapy. The relevance of the identified DMRs is supported by their association with key regulatory processes of oncogenic transformation and support the idea of relevant DMRs being enriched at distinct loci rather than evenly distribution across the genome. Predictive response to therapy could be established but lacked specificity for treatment with azacytidine. Our results suggest that a predictive epigenotype carries its methylation information at a complex, genome-wide level, that is confined to regions, rather than to single CpGs. With increasing application of combinatorial regimens, response prediction may become even more complicated.

Lack of predictive tools for conventional and targeted cancer therapy: Barriers to biomarker development and clinical translation

Predictive tools, utilising biomarkers, aim to objectively assessthe potentialresponse toa particular clinical intervention in order to direct treatment.Conventional cancer therapy remains poorly served by predictive biomarkers, despite being the mainstay of treatment for most patients. In contrast, targeted therapy benefits from a clearly defined protein target for potential biomarker assessment. We discuss potential data sources of predictive biomarkers for conventional and targeted therapy, including patient clinical data andmulti-omicbiomarkers (genomic, transcriptomic and protein expression).Key examples, either clinically adopted or demonstrating promise for clinical translation, are highlighted. Following this, we provide an outline of potential barriers to predictive biomarker development; broadly discussing themes of approaches to translational research and study/trial design, and the impact of cellular and molecular tumor heterogeneity. Future avenues of research are also highlighted.

The tumor immune microenvironment and immune-related signature predict the chemotherapy response in patients with osteosarcoma

BACKGROUND

Genome-wide expression profiles have been shown to predict the response to chemotherapy. The purpose of this study was to develop a novel predictive signature for chemotherapy in patients with osteosarcoma.

METHODS

We analysed the relevance of immune cell infiltration and gene expression profiles of the tumor samples of good responders with those of poor responders from the TARGET and GEO databases. Immune cell infiltration was evaluated using a single-sample gene set enrichment analysis (ssGSEA) and the CIBERSORT algorithm between good and poor chemotherapy responders. Differentially expressed genes were identified based on the chemotherapy response. LASSO regression and binary logistic regression analyses were applied to select the differentially expressed immune-related genes (IRGs) and developed a predictive signature in the training cohort. A receiver operating characteristic (ROC) curve analysis was employed to assess and validate the predictive accuracy of the predictive signature in the validation cohort.

RESULTS

The analysis of immune infiltration showed a positive relationship between high-level immune infiltration and good responders, and T follicular helper cells and CD8 T cells were significantly more abundant in good responders with osteosarcoma. Two hundred eighteen differentially expressed genes were detected between good and poor responders, and a five IRGs panel comprising TNFRSF9, CD70, EGFR, PDGFD and S100A6 was determined to show predictive power for the chemotherapy response. A chemotherapy-associated predictive signature was developed based on these five IRGs. The accuracy of the predictive signature was 0.832 for the training cohort and 0.720 for the validation cohort according to ROC analysis.

CONCLUSIONS

The novel predictive signature constructed with five IRGs can be effectively utilized to predict chemotherapy responsiveness and help improve the efficacy of chemotherapy in patients with osteosarcoma.

Construction and optimization of gene expression signatures for prediction of survival in two-arm clinical trials

Background

Gene expression signatures for the prediction of differential survival of patients undergoing anti-cancer therapies are of great interest because they can be used to prospectively stratify patients entering new clinical trials, or to determine optimal treatment for patients in more routine clinical settings. Unlike prognostic signatures however, predictive signatures require training set data from clinical studies with at least two treatment arms. As two-arm studies with gene expression profiling have been rarer than similar one-arm studies, the methodology for constructing and optimizing predictive signatures has been less prominently explored than for prognostic signatures.

Results

Focusing on two “use cases” of two-arm clinical trials, one for metastatic colorectal cancer (CRC) patients treated with the anti-angiogenic molecule aflibercept, and the other for triple negative breast cancer (TNBC) patients treated with the small molecule iniparib, we present derivation steps and quantitative and graphical tools for the construction and optimization of signatures for the prediction of progression-free survival based on cross-validated multivariate Cox models. This general methodology is organized around two more specific approaches which we have called subtype correlation (subC) and mechanism-of-action (MOA) modeling, each of which leverage a priori knowledge of molecular subtypes of tumors or drug MOA for a given indication. The tools and concepts presented here include the so-called differential log-hazard ratio, the survival scatter plot, the hazard ratio receiver operating characteristic, the area between curves and the patient selection matrix. In the CRC use case for instance, the resulting signature stratifies the patient population into “sensitive” and “relatively-resistant” groups achieving a more than two-fold difference in the aflibercept-to-control hazard ratios across signature-defined patient groups. Through cross-validation and resampling the probability of generalization of the signature to similar CRC data sets is predicted to be high.

Conclusions

The tools presented here should be of general use for building and using predictive multivariate signatures in oncology and in other therapeutic areas.

The role of molecular enrichment on future therapies in hepatocellular carcinoma

Hepatocellular carcinomas (HCCs) are characterised by considerable phenotypic and molecular heterogeneity. Treating HCC and designing clinical trials are particularly challenging because co-existing liver disease, present in most patients, limits aggressive therapeutic options. Positive results in recent phase III clinical trials have confirmed the high value of anti-angiogenic therapies for HCC in both first (sorafenib and lenvatinib) and second line (regorafenib and cabozantinib) treatment modalities. However, failure of several large randomised controlled clinical trials over the last 10 years underlines the necessity for innovative treatment strategies and implementation of translational findings to overcome the unmet clinical need. Furthermore, the promising results from novel immunotherapies are likely to complement the landscape of active compounds for HCC and will require a completely different approach to patients, as well as the development of prognostic/predictive biomarkers. Given our increasing understanding of the most abundant molecular alterations in HCC, effective enrichment of patients based on clinical and molecular biomarkers, as well as adaptive clinical trials, are now feasible and should be implemented. Herein, we aim to review important aspects of precision medicine approaches in HCC that might contribute to improving the molecular subclassification of patients in a clinical trial setting and pave the way for novel therapeutic strategies.