A genomic approach to improve prognosis and predict therapeutic response in chronic lymphocytic leukemia

Multi-omics reveals clinically relevant proliferative drive associated with mTOR-MYC-OXPHOS activity in chronic lymphocytic leukemia

Chronic Lymphocytic Leukemia (CLL) has a complex pattern of driver mutations and much of its clinical diversity remains unexplained. We devised a method for simultaneous subgroup discovery across multiple data types and applied it to genomic, transcriptomic, DNA methylation and ex-vivo drug response data from 217 Chronic Lymphocytic Leukemia (CLL) cases. We uncovered a biological axis of heterogeneity strongly associated with clinical behavior and orthogonal to the known biomarkers. We validated its presence and clinical relevance in four independent cohorts (n=547 patients). We find that this axis captures the proliferative drive (PD) of CLL cells, as it associates with lymphocyte doubling rate, global hypomethylation, accumulation of driver aberrations and response to pro-proliferative stimuli. CLL-PD was linked to the activation of mTOR-MYC-oxidative phosphorylation (OXPHOS) through transcriptomic, proteomic and single cell resolution analysis. CLL-PD is a key determinant of disease outcome in CLL. Our multi-table integration approach may be applicable to other tumors whose inter-individual differences are currently unexplained.

Comparing the Predictivity of Human Placental Gene, microRNA, and CpG Methylation Signatures in Relation to Perinatal Outcomes

Molecular signatures are being increasingly integrated into predictive biology applications. However, there are limited studies comparing the overall predictivity of transcriptomic vs. epigenomic signatures in relation to perinatal outcomes. This study set out to evaluate mRNA and microRNA (miRNA) expression and cytosine-guanine dinucleotide (CpG) methylation signatures in human placental tissues and relate these to perinatal outcomes known to influence maternal/fetal health; namely, birth weight, placenta weight, placental damage, and placental inflammation. The following hypotheses were tested: (1) different molecular signatures will demonstrate varying levels of predictivity towards perinatal outcomes, and (2) these signatures will show disruptions from an example exposure (i.e., cadmium) known to elicit perinatal toxicity. Multi-omic placental profiles from 390 infants in the Extremely Low Gestational Age Newborns cohort were used to develop molecular signatures that predict each perinatal outcome. Epigenomic signatures (i.e., miRNA and CpG methylation) consistently demonstrated the highest levels of predictivity, with model performance metrics including R^2 (predicted vs. observed) values of 0.36-0.57 for continuous outcomes and balanced accuracy values of 0.49-0.77 for categorical outcomes. Top-ranking predictors included miRNAs involved in injury and inflammation. To demonstrate the utility of these predictive signatures in screening of potentially harmful exogenous insults, top-ranking miRNA predictors were analyzed in a separate pregnancy cohort and related to cadmium. Key predictive miRNAs demonstrated altered expression in association with cadmium exposure, including miR-210, known to impact placental cell growth, blood vessel development, and fetal weight. These findings inform future predictive biology applications, where additional benefit will be gained by including epigenetic markers.

HELQ and EGR3 expression correlate with IGHV mutation status and prognosis in chronic lymphocytic leukemia

Background

IGHV mutation status is a crucial prognostic biomarker for CLL. In the present study, we investigated the transcriptomic signatures associating with IGHV mutation status and CLL prognosis.

Methods

The co-expression modules and hub genes correlating with IGHV status, were identified using the GSE28654, by ‘WGCNA’ package and R software (version 4.0.2). The over-representation analysis was performed to reveal enriched cell pathways for genes of correlating modules. Then 9 external cohorts were used to validate the correlation of hub genes expression with IGHV status or clinical features (treatment response, transformation to Richter syndrome, etc.). Moreover, to elucidate the significance of hub genes on disease course and prognosis of CLL patients, the Kaplan–Meier analysis for the OS and TTFT of were performed between subgroups dichotomized by the median expression value of individual hub genes.

Results

2 co-expression modules and 9 hub genes ((FCRL1/FCRL2/HELQ/EGR3LPL/LDOC1/ZNF667/SOWAHC/SEPTIN10) correlating with IGHV status were identified by WGCNA, and validated by external datasets. The modules were found to be enriched in NF-kappaB, HIF-1 and other important pathways, involving cell proliferation and apoptosis. The expression of hub genes was revealed to be significantly different, not only between CLL and normal B cell, but also between various types of lymphoid neoplasms. HELQ expression was found to be related with response of immunochemotherapy treatment significantly (p = 0.0413), while HELQ and ZNF667 were expressed differently between stable CLL and Richter syndrome patients (p < 0.0001 and p = 0.0278, respectively). By survival analysis of subgroups, EGR3 expression was indicated to be significantly associated with TTFT by 2 independent cohorts (GSE39671, p = 0.0311; GSE22762, p = 0.0135). While the expression of HELQ and EGR3 was found to be associated with OS (p = 0.0291 and 0.0114 respectively).The Kras, Hedgehog and IL6-JAK-STAT3 pathways were found to be associating with the expression of hub genes, resulting from GSEA.

Conclusions

The expression of HELQ and EGR3 were correlated with IGHV mutation status in CLL patients. Additionally, the expression of HELQ/EGR3 were prognostic markers for CLL associating with targetable cell signaling pathways.

Prognostic models for newly-diagnosed chronic lymphocytic leukaemia in adults: a systematic review and meta-analysis

BACKGROUND

Chronic lymphocytic leukaemia (CLL) is the most common cancer of the lymphatic system in Western countries. Several clinical and biological factors for CLL have been identified. However, it remains unclear which of the available prognostic models combining those factors can be used in clinical practice to predict long-term outcome in people newly-diagnosed with CLL.

OBJECTIVES

To identify, describe and appraise all prognostic models developed to predict overall survival (OS), progression-free survival (PFS) or treatment-free survival (TFS) in newly-diagnosed (previously untreated) adults with CLL, and meta-analyse their predictive performances.

SEARCH METHODS

We searched MEDLINE (from January 1950 to June 2019 via Ovid), Embase (from 1974 to June 2019) and registries of ongoing trials (to 5 March 2020) for development and validation studies of prognostic models for untreated adults with CLL. In addition, we screened the reference lists and citation indices of included studies.

SELECTION CRITERIA

We included all prognostic models developed for CLL which predict OS, PFS, or TFS, provided they combined prognostic factors known before treatment initiation, and any studies that tested the performance of these models in individuals other than the ones included in model development (i.e. 'external model validation studies'). We included studies of adults with confirmed B-cell CLL who had not received treatment prior to the start of the study. We did not restrict the search based on study design.

DATA COLLECTION AND ANALYSIS

We developed a data extraction form to collect information based on the Checklist for Critical Appraisal and Data Extraction for Systematic Reviews of Prediction Modelling Studies (CHARMS). Independent pairs of review authors screened references, extracted data and assessed risk of bias according to the Prediction model Risk Of Bias ASsessment Tool (PROBAST). For models that were externally validated at least three times, we aimed to perform a quantitative meta-analysis of their predictive performance, notably their calibration (proportion of people predicted to experience the outcome who do so) and discrimination (ability to differentiate between people with and without the event) using a random-effects model. When a model categorised individuals into risk categories, we pooled outcome frequencies per risk group (low, intermediate, high and very high). We did not apply GRADE as guidance is not yet available for reviews of prognostic models.

MAIN RESULTS

From 52 eligible studies, we identified 12 externally validated models: six were developed for OS, one for PFS and five for TFS. In general, reporting of the studies was poor, especially predictive performance measures for calibration and discrimination; but also basic information, such as eligibility criteria and the recruitment period of participants was often missing. We rated almost all studies at high or unclear risk of bias according to PROBAST. Overall, the applicability of the models and their validation studies was low or unclear; the most common reasons were inappropriate handling of missing data and serious reporting deficiencies concerning eligibility criteria, recruitment period, observation time and prediction performance measures. We report the results for three models predicting OS, which had available data from more than three external validation studies: CLL International Prognostic Index (CLL-IPI) This score includes five prognostic factors: age, clinical stage, IgHV mutational status, B2-microglobulin and TP53 status. Calibration: for the low-, intermediate- and high-risk groups, the pooled five-year survival per risk group from validation studies corresponded to the frequencies observed in the model development study. In the very high-risk group, predicted survival from CLL-IPI was lower than observed from external validation studies. Discrimination: the pooled c-statistic of seven external validation studies (3307 participants, 917 events) was 0.72 (95% confidence interval (CI) 0.67 to 0.77). The 95% prediction interval (PI) of this model for the c-statistic, which describes the expected interval for the model's discriminative ability in a new external validation study, ranged from 0.59 to 0.83. Barcelona-Brno score Aimed at simplifying the CLL-IPI, this score includes three prognostic factors: IgHV mutational status, del(17p) and del(11q). Calibration: for the low- and intermediate-risk group, the pooled survival per risk group corresponded to the frequencies observed in the model development study, although the score seems to overestimate survival for the high-risk group. Discrimination: the pooled c-statistic of four external validation studies (1755 participants, 416 events) was 0.64 (95% CI 0.60 to 0.67); 95% PI 0.59 to 0.68. MDACC 2007 index score The authors presented two versions of this model including six prognostic factors to predict OS: age, B2-microglobulin, absolute lymphocyte count, gender, clinical stage and number of nodal groups. Only one validation study was available for the more comprehensive version of the model, a formula with a nomogram, while seven studies (5127 participants, 994 events) validated the simplified version of the model, the index score. Calibration: for the low- and intermediate-risk groups, the pooled survival per risk group corresponded to the frequencies observed in the model development study, although the score seems to overestimate survival for the high-risk group. Discrimination: the pooled c-statistic of the seven external validation studies for the index score was 0.65 (95% CI 0.60 to 0.70); 95% PI 0.51 to 0.77.

AUTHORS' CONCLUSIONS

Despite the large number of published studies of prognostic models for OS, PFS or TFS for newly-diagnosed, untreated adults with CLL, only a minority of these (N = 12) have been externally validated for their respective primary outcome. Three models have undergone sufficient external validation to enable meta-analysis of the model's ability to predict survival outcomes. Lack of reporting prevented us from summarising calibration as recommended. Of the three models, the CLL-IPI shows the best discrimination, despite overestimation. However, performance of the models may change for individuals with CLL who receive improved treatment options, as the models included in this review were tested mostly on retrospective cohorts receiving a traditional treatment regimen. In conclusion, this review shows a clear need to improve the conducting and reporting of both prognostic model development and external validation studies. For prognostic models to be used as tools in clinical practice, the development of the models (and their subsequent validation studies) should adapt to include the latest therapy options to accurately predict performance. Adaptations should be timely.

P-STAT3 and IL-17 in tumor tissues enhances the prognostic value of CEA and CA125 in patients with lung adenocarcinoma

AIM

The present study aimed to examine the capability of p- signal transducer and activator of transcription (STAT)3 and interleukin-17 (IL-17), along with two known tumor markers carcinoembryonic antigen (CEA) and carbohydrate antigen 125 (CA125), for disease prognosis. Moreover, the associations among biomarkers and clinicopathological parameters were evaluated to uncover the potential mechanisms responsible for their correlations with lung adenocarcinoma (LAD) prognosis.

METHODS

Five LAD-related parameters were used in the study: CEA, CA125, STAT3, p-STAT3, and IL-17. Spearman and chi-square correlation tests were used to explore the relationships between some clinicopathological variables and parameter expression levels and the associations among these five parameters.

RESULTS

The disease-specific survival decreased with the positive expression of CEA, CA125, p-STAT3, and IL-17, with no significant difference in the expression level of STAT3. Combinations of p-STAT3 and IL-17, CEA and p-STAT3, CEA and IL-17, CA125 and p-STAT3, and CA125 and IL-17 had higher predictive values in LAD prognosis. The correlation analyses indicated the synergic activities of STAT3, p-STAT3, and IL-17 and the coordinated expression of CEA, CA125, p-STAT3, and IL-17. The tumor-node-metastasis (TNM) stage significantly correlated with the levels of CA125 and p-STAT3.

CONCLUSIONS

Elevated levels of CEA, CA125, p-STAT3, and IL-17 alone and/or combinations of p-STAT3 and IL-17, CEA and p-STAT3, CEA and IL-17, CA125 and p-STAT3, and CA125 and IL-17 were recommended as the prognostic predictors of unfavorable clinical outcomes in patients with postoperative LAD. Also, p-STAT3 and IL-17 combined with CA125 and CEA helped in predicting the overall survival of patients with LAD and informing the TNM stage.

Reflections on and future of hematologic malignancies research in the Veterans Health Administration

Research in the Veterans Health Administration (VHA) has played an integral part in learning about cancer biology and treatment. Here we provide examples of past research performed in the VHA focusing on hematologic malignancies, and identify future opportunities for areas of research in this group of uncommon diseases that have specific importance for Veterans and the VHA. Veterans treated in the VHA and in the private sector deserve information that is focused on them, and is not an extrapolation from the larger population. Only by building upon and expanding existing research within the VHA can Veteran-specific results be collected and best practices be developed. In turn, such advances will benefit Veterans affected by these cancers with an improved quality of life and a longer lifespan.

Time to Treatment Prediction in Chronic Lymphocytic Leukemia Based on New Transcriptional Patterns

Chronic lymphocytic leukemia (CLL) is the most frequent lymphoproliferative syndrome in western countries. CLL evolution is frequently indolent, and treatment is mostly reserved for those patients with signs or symptoms of disease progression. In this work, we used RNA sequencing data from the International Cancer Genome Consortium CLL cohort to determine new gene expression patterns that correlate with clinical evolution.We determined that a 290-gene expression signature, in addition to immunoglobulin heavy chain variable region (IGHV) mutation status, stratifies patients into four groups with notably different time to first treatment. This finding was confirmed in an independent cohort. Similarly, we present a machine learning algorithm that predicts the need for treatment within the first 5 years following diagnosis using expression data from 2,198 genes. This predictor achieved 90% precision and 89% accuracy when classifying independent CLL cases. Our findings indicate that CLL progression risk largely correlates with particular transcriptomic patterns and paves the way for the identification of high-risk patients who might benefit from prompt therapy following diagnosis.

Subtype assignment of CLL based on B-cell subset associated gene signatures from normal bone marrow - A proof of concept study

Diagnostic and prognostic evaluation of chronic lymphocytic leukemia (CLL) involves blood cell counts, immunophenotyping, IgVH mutation status, and cytogenetic analyses. We generated B-cell associated gene-signatures (BAGS) based on six naturally occurring B-cell subsets within normal bone marrow. Our hypothesis is that by segregating CLL according to BAGS, we can identify subtypes with prognostic implications in support of pathogenetic value of BAGS. Microarray-based gene-expression samples from eight independent CLL cohorts (1,024 untreated patients) were BAGS-stratified into pre-BI, pre-BII, immature, naïve, memory, or plasma cell subtypes; the majority falling within the memory (24.5-45.8%) or naïve (14.5-32.3%) categories. For a subset of CLL patients (n = 296), time to treatment (TTT) was shorter amongst early differentiation subtypes (pre-BI/pre-BII/immature) compared to late subtypes (memory/plasma cell, HR: 0.53 [0.35-0.78]). Particularly, pre-BII subtype patients had the shortest TTT among all subtypes. Correlates derived for BAGS subtype and IgVH mutation (n = 405) revealed an elevated mutation frequency in late vs. early subtypes (71% vs. 45%, P < .001). Predictions for BAGS subtype resistance towards rituximab and cyclophosphamide varied for rituximab, whereas all subtypes were sensitive to cyclophosphamide. This study supports our hypothesis that BAGS-subtyping may be of tangible prognostic and pathogenetic value for CLL patients.

sPAGM: inferring subpathway activity by integrating gene and miRNA expression-robust functional signature identification for melanoma prognoses

MicroRNAs (miRNAs) regulate biological pathways by inhibiting gene expression. However, most current analytical methods fail to consider miRNAs, when inferring functional or pathway activities. In this study, we developed a model called sPAGM to infer subpathway activities by integrating gene and miRNA expressions. In this model, we reconstructed subpathway graphs by embedding miRNA components, and characterized subpathway activity (sPA) scores by simultaneously considering the expression levels of miRNAs and genes. The results showed that the sPA scores could distinguish different samples across tumor types, as well as samples between tumor and normal conditions. Moreover, the sPAGM model displayed more specificities than the entire pathway-based analyses. This model was applied to melanoma tumors to perform a prognosis analysis, which identified a robust 55-subpathway signature. By using The Cancer Genome Atlas and independently verified data sets, the subpathway-based signature significantly predicted the patients’ prognoses, which were independent of clinical variables. In the prognostic performance comparison, the sPAGM model was superior to the gene-only and miRNA-only methods. Finally, we dissected the functional roles and interactions of components within the subpathway signature. Taken together, the sPAGM model provided a framework for inferring subpathway activities and identifying functional signatures for clinical applications.

Aberrantly expressed LGR4 empowers Wnt signaling in multiple myeloma by hijacking osteoblast-derived R-spondins

The unrestrained growth of tumor cells is generally attributed to mutations in essential growth control genes, but tumor cells are also affected by, or even addicted to, signals from the microenvironment. As therapeutic targets, these extrinsic signals may be equally significant as mutated oncogenes. In multiple myeloma (MM), a plasma cell malignancy, most tumors display hallmarks of active Wnt signaling but lack activating Wnt-pathway mutations, suggesting activation by autocrine Wnt ligands and/or paracrine Wnts emanating from the bone marrow (BM) niche. Here, we report a pivotal role for the R-spondin/leucine-rich repeat-containing G protein-coupled receptor 4 (LGR4) axis in driving aberrant Wnt/β-catenin signaling in MM. We show that LGR4 is expressed by MM plasma cells, but not by normal plasma cells or B cells. This aberrant LGR4 expression is driven by IL-6/STAT3 signaling and allows MM cells to hijack R-spondins produced by (pre)osteoblasts in the BM niche, resulting in Wnt (co)receptor stabilization and a dramatically increased sensitivity to auto- and paracrine Wnts. Our study identifies aberrant R-spondin/LGR4 signaling with consequent deregulation of Wnt (co)receptor turnover as a driver of oncogenic Wnt/β-catenin signaling in MM cells. These results advocate targeting of the LGR4/R-spondin interaction as a therapeutic strategy in MM.

Colony-Stimulating Factor-1 Receptor Is Required for Nurse-like Cell Survival in Chronic Lymphocytic Leukemia

PURPOSE

Monocytes and their progeny are abundant constituents of the tumor microenvironment in lymphoproliferative disorders, including chronic lymphocytic leukemia (CLL). Monocyte-derived cells, including nurse-like cells (NLC) in CLL, promote lymphocyte proliferation and survival, confer resistance to chemotherapy, and are associated with more rapid disease progression. Colony-stimulating factor-1 receptor (CSF-1R) regulates the homeostatic survival of tissue-resident macrophages. Therefore, we sought to determine whether CSF-1R is similarly required for NLC survival.

EXPERIMENTAL DESIGN

CSF-1R expression by NLC was examined by flow cytometry and IHC. CSF-1R blocking studies were performed using an antagonistic mAb to examine its role in NLC generation and in CLL survival. A rational search strategy was performed to identify a novel tyrosine kinase inhibitor (TKI) targeting CSF-1R. The influence of TKI-mediated CSF-1R inhibition on NLC and CLL viability was examined.

RESULTS

We demonstrated that the generation and survival of NLC in CLL is dependent upon CSF-1R signaling. CSF-1R blockade is associated with significant depletion of NLC and consequently inhibits CLL B-cell survival. We found that the JAK2/FLT3 inhibitor pacritinib suppresses CSF-1R signaling, thereby preventing the generation and survival of NLC and impairs CLL B-cell viability.

CONCLUSIONS

CSF-1R is a novel therapeutic target that may be exploited in lymphoproliferative disorders, like CLL, that are dependent upon lymphoma-associated macrophages. Clin Cancer Res; 22(24); 6118-28. ©2016 AACR.

Comparative study of joint analysis of microarray gene expression data in survival prediction and risk assessment of breast cancer patients

Microarray gene expression data sets are jointly analyzed to increase statistical power. They could either be merged together or analyzed by meta-analysis. For a given ensemble of data sets, it cannot be foreseen which of these paradigms, merging or meta-analysis, works better. In this article, three joint analysis methods, Z -score normalization, ComBat and the inverse normal method (meta-analysis) were selected for survival prognosis and risk assessment of breast cancer patients. The methods were applied to eight microarray gene expression data sets, totaling 1324 patients with two clinical endpoints, overall survival and relapse-free survival. The performance derived from the joint analysis methods was evaluated using Cox regression for survival analysis and independent validation used as bias estimation. Overall, Z -score normalization had a better performance than ComBat and meta-analysis. Higher Area Under the Receiver Operating Characteristic curve and hazard ratio were also obtained when independent validation was used as bias estimation. With a lower time and memory complexity, Z -score normalization is a simple method for joint analysis of microarray gene expression data sets. The derived findings suggest further assessment of this method in future survival prediction and cancer classification applications.

Methods: for studying pharmacogenetic profiles of combination chemotherapeutic drugs

INTRODUCTION

Molecular and genetic analysis of tumors and individuals has led to patient-centered therapies, through the discovery and identification of genetic markers predictive of drug efficacy and tolerability. Present therapies often include a combination of synergic drugs, each of them directed against different targets. Therefore, the pharmacogenetic profiling of tumor masses and patients is becoming a challenge, and several questions may arise when planning a translational study.

AREAS COVERED

The review presents the different techniques used to stratify oncology patients and to tailor antineoplastic treatments according to individual pharmacogenetic profiling. The advantages of these methodologies are discussed as well as current limits.

EXPERT OPINION

Facing the rapid technological evolution for genetic analyses, the most pressing issues are the choice of appropriate strategies (i.e., from gene candidate up to next-generation sequencing) and the possibility to replicate study results for their final validation. It is likely that the latter will be the major obstacle in the future. However, the present landscape is opening up new possibilities, overcoming those hurdles that have limited result translation into clinical settings for years.

The expression of RUNDC3B is associated with promoter methylation in lymphoid malignancies

DNA methylation is an epigenetic modification that plays an important role in the regulation of gene expression. The function of RUNDC3B has yet to be determined, although its dysregulated expression has been associated with malignant potential of both breast and lung carcinoma. To elucidate the potential of using DNA methylation in RUNDC3B as a biomarker in lymphoid malignancies, the methylation status of six regions spanning the CpG island in the promoter region of RUNDC3B was determined in cancer cell lines. Lymphoid malignancies were found to have more prominent methylation and did not express RUNDC3B compared with myeloid malignancies and solid tumours, supporting the potential use of DNA methylation in this region as a biomarker for lymphoid malignancies. RUNDC3B contains a RUN domain in its N‐terminal region that mediates interaction with Rap2, an important component of the mitogen‐activated protein kinase (MAPK) cascade, which regulates cellular proliferation and differentiation. The protein sequence of RUNDC3B also contains characteristic binding sites for MAPK intermediates. Therefore, it is possible that RUNDC3B serves as a mediator between Rap2 and the MAPK signalling cascade. Three genes with MAPK‐inducible expression were downregulated in a methylated leukaemia cell line (HSPA5, Jun and Fos). Jun and Fos combine to form the activating protein 1 transcription factor, and loss of this factor is associated with the dysregulation of genes involved in differentiation and proliferation. We hypothesize that the loss of RUNDC3B secondary to aberrant hypermethylation of the early growth response 3 transcription factor binding site results in dysregulated MAPK signalling and carcinogenesis in lymphoid malignancies. © 2015 The Authors. Hematological Oncology published by John Wiley & Sons Ltd

Pathway activity inference for multiclass disease classification through a mathematical programming optimisation framework

BACKGROUND

Applying machine learning methods on microarray gene expression profiles for disease classification problems is a popular method to derive biomarkers, i.e. sets of genes that can predict disease state or outcome. Traditional approaches where expression of genes were treated independently suffer from low prediction accuracy and difficulty of biological interpretation. Current research efforts focus on integrating information on protein interactions through biochemical pathway datasets with expression profiles to propose pathway-based classifiers that can enhance disease diagnosis and prognosis. As most of the pathway activity inference methods in literature are either unsupervised or applied on two-class datasets, there is good scope to address such limitations by proposing novel methodologies.

RESULTS

A supervised multiclass pathway activity inference method using optimisation techniques is reported. For each pathway expression dataset, patterns of its constituent genes are summarised into one composite feature, termed pathway activity, and a novel mathematical programming model is proposed to infer this feature as a weighted linear summation of expression of its constituent genes. Gene weights are determined by the optimisation model, in a way that the resulting pathway activity has the optimal discriminative power with regards to disease phenotypes. Classification is then performed on the resulting low-dimensional pathway activity profile.

CONCLUSIONS

The model was evaluated through a variety of published gene expression profiles that cover different types of disease. We show that not only does it improve classification accuracy, but it can also perform well in multiclass disease datasets, a limitation of other approaches from the literature. Desirable features of the model include the ability to control the maximum number of genes that may participate in determining pathway activity, which may be pre-specified by the user. Overall, this work highlights the potential of building pathway-based multi-phenotype classifiers for accurate disease diagnosis and prognosis problems.

A three step network based approach (TSNBA) to finding disease molecular signature and key regulators: a case study of IL-1 and TNF-alpha stimulated inflammation

A disease molecular signature is a set of biomolecular features that are prognostic of clinical phenotypes and indicative of underlying pathology. It is of great importance to develop computational approaches for finding more relevant molecular signatures. Based upon the hypothesis that various components in a molecular signature are more likely to share similar patterns, we introduced a novel three step network based approach (TSNBA) to identify the molecular signature and key pathological regulators. Protein-protein interaction (PPI) network and ranking algorithm were integrated in the first step to find pathology related proteins with high accuracy. It was followed by the second step to further screen with co-expression patterns for better pathology enrichment. Context likelihood of relatedness (CLR) algorithm was used in the third step to infer gene regulatory networks and identify key transcription regulators. We applied this approach to study IL-1 (interleukin-1) and TNF-alpha (tumor necrosis factor-alpha) stimulated inflammation. TSNBA identified inflammatory signature with high accuracy and outperformed 5 competing methods namely fold change, degree, interconnectivity, neighborhood score and network propagation based approaches. The best molecular signature, with 80% (40/50) confirmed inflammatory genes, was used to predict inflammation related genes. As a result, 8 out of 10 predicted inflammation genes that were not included in the benchmark Entrez Gene database were validated by literature evidence. Furthermore, 23 of the 32 predicted inflammation regulators were validated by literature evidence. The rest 9 were also validated with TF (transcription factor) binding site analysis. In conclusion, we developed an efficient strategy for disease molecular signature finding and key pathological regulator identification.

Transcriptome characterization by RNA sequencing identifies a major molecular and clinical subdivision in chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) has heterogeneous clinical and biological behavior. Whole-genome and -exome sequencing has contributed to the characterization of the mutational spectrum of the disease, but the underlying transcriptional profile is still poorly understood. We have performed deep RNA sequencing in different subpopulations of normal B-lymphocytes and CLL cells from a cohort of 98 patients, and characterized the CLL transcriptional landscape with unprecedented resolution. We detected thousands of transcriptional elements differentially expressed between the CLL and normal B cells, including protein-coding genes, noncoding RNAs, and pseudogenes. Transposable elements are globally derepressed in CLL cells. In addition, two thousand genes—most of which are not differentially expressed—exhibit CLL-specific splicing patterns. Genes involved in metabolic pathways showed higher expression in CLL, while genes related to spliceosome, proteasome, and ribosome were among the most down-regulated in CLL. Clustering of the CLL samples according to RNA-seq derived gene expression levels unveiled two robust molecular subgroups, C1 and C2. C1/C2 subgroups and the mutational status of the immunoglobulin heavy variable (IGHV) region were the only independent variables in predicting time to treatment in a multivariate analysis with main clinico-biological features. This subdivision was validated in an independent cohort of patients monitored through DNA microarrays. Further analysis shows that B-cell receptor (BCR) activation in the microenvironment of the lymph node may be at the origin of the C1/C2 differences.

Multi-study integration of brain cancer transcriptomes reveals organ-level molecular signatures

We utilized abundant transcriptomic data for the primary classes of brain cancers to study the feasibility of separating all of these diseases simultaneously based on molecular data alone. These signatures were based on a new method reported herein – Identification of Structured Signatures and Classifiers (ISSAC) – that resulted in a brain cancer marker panel of 44 unique genes. Many of these genes have established relevance to the brain cancers examined herein, with others having known roles in cancer biology. Analyses on large-scale data from multiple sources must deal with significant challenges associated with heterogeneity between different published studies, for it was observed that the variation among individual studies often had a larger effect on the transcriptome than did phenotype differences, as is typical. For this reason, we restricted ourselves to studying only cases where we had at least two independent studies performed for each phenotype, and also reprocessed all the raw data from the studies using a unified pre-processing pipeline. We found that learning signatures across multiple datasets greatly enhanced reproducibility and accuracy in predictive performance on truly independent validation sets, even when keeping the size of the training set the same. This was most likely due to the meta-signature encompassing more of the heterogeneity across different sources and conditions, while amplifying signal from the repeated global characteristics of the phenotype. When molecular signatures of brain cancers were constructed from all currently available microarray data, 90% phenotype prediction accuracy, or the accuracy of identifying a particular brain cancer from the background of all phenotypes, was found. Looking forward, we discuss our approach in the context of the eventual development of organ-specific molecular signatures from peripheral fluids such as the blood.

From a multi-study, integrated transcriptomic dataset, we identified a marker panel for differentiating major human brain cancers at the gene-expression level. The ISSAC molecular signatures for brain cancers, composed of 44 unique genes, are based on comparing expression levels of pairs of genes, and phenotype prediction follows a diagnostic hierarchy. We found that sufficient dataset integration across multiple studies greatly enhanced diagnostic performance on truly independent validation sets, whereas signatures learned from only one dataset typically led to high error rate. Molecular signatures of brain cancers, when obtained using all currently available gene-expression data, achieved 90% phenotype prediction accuracy. Thus, our integrative approach holds significant promise for developing organ-level, comprehensive, molecular signatures of disease.

Clinical and biological relevance of genomic heterogeneity in chronic lymphocytic leukemia

BACKGROUND

Chronic lymphocytic leukemia (CLL) is typically regarded as an indolent B-cell malignancy. However, there is wide variability with regards to need for therapy, time to progressive disease, and treatment response. This clinical variability is due, in part, to biological heterogeneity between individual patients' leukemias. While much has been learned about this biological variation using genomic approaches, it is unclear whether such efforts have sufficiently evaluated biological and clinical heterogeneity in CLL.

METHODS

To study the extent of genomic variability in CLL and the biological and clinical attributes of genomic classification in CLL, we evaluated 893 unique CLL samples from fifteen publicly available gene expression profiling datasets. We used unsupervised approaches to divide the data into subgroups, evaluated the biological pathways and genetic aberrations that were associated with the subgroups, and compared prognostic and clinical outcome data between the subgroups.

RESULTS

Using an unsupervised approach, we determined that approximately 600 CLL samples are needed to define the spectrum of diversity in CLL genomic expression. We identified seven genomically-defined CLL subgroups that have distinct biological properties, are associated with specific chromosomal deletions and amplifications, and have marked differences in molecular prognostic markers and clinical outcomes.

CONCLUSIONS

Our results indicate that investigations focusing on small numbers of patient samples likely provide a biased outlook on CLL biology. These findings may have important implications in identifying patients who should be treated with specific targeted therapies, which could have efficacy against CLL cells that rely on specific biological pathways.

Network information improves cancer outcome prediction

Disease progression in cancer can vary substantially between patients. Yet, patients often receive the same treatment. Recently, there has been much work on predicting disease progression and patient outcome variables from gene expression in order to personalize treatment options. Despite first diagnostic kits in the market, there are open problems such as the choice of random gene signatures or noisy expression data. One approach to deal with these two problems employs protein-protein interaction networks and ranks genes using the random surfer model of Google's PageRank algorithm. In this work, we created a benchmark dataset collection comprising 25 cancer outcome prediction datasets from literature and systematically evaluated the use of networks and a PageRank derivative, NetRank, for signature identification. We show that the NetRank performs significantly better than classical methods such as fold change or t-test. Despite an order of magnitude difference in network size, a regulatory and protein-protein interaction network perform equally well. Experimental evaluation on cancer outcome prediction in all of the 25 underlying datasets suggests that the network-based methodology identifies highly overlapping signatures over all cancer types, in contrast to classical methods that fail to identify highly common gene sets across the same cancer types. Integration of network information into gene expression analysis allows the identification of more reliable and accurate biomarkers and provides a deeper understanding of processes occurring in cancer development and progression.

Gene signature critical to cancer phenotype as a paradigm for anticancer drug discovery

Malignant cell transformation commonly results in the deregulation of thousands of cellular genes, an observation that suggests a complex biological process and an inherently challenging scenario for the development of effective cancer interventions. To better define the genes/pathways essential to regulating the malignant phenotype, we recently described a novel strategy based on the cooperative nature of carcinogenesis that focuses on genes synergistically deregulated in response to cooperating oncogenic mutations. These so-called 'cooperation response genes' (CRGs) are highly enriched for genes critical for the cancer phenotype, thereby suggesting their causal role in the malignant state. Here, we show that CRGs have an essential role in drug-mediated anticancer activity and that anticancer agents can be identified through their ability to antagonize the CRG expression profile. These findings provide proof-of-concept for the use of the CRG signature as a novel means of drug discovery with relevance to underlying anticancer drug mechanisms.

Molecular signatures from omics data: from chaos to consensus

In the past 15 years, new "omics" technologies have made it possible to obtain high-resolution molecular snapshots of organisms, tissues, and even individual cells at various disease states and experimental conditions. It is hoped that these developments will usher in a new era of personalized medicine in which an individual's molecular measurements are used to diagnose disease, guide therapy, and perform other tasks more accurately and effectively than is possible using standard approaches. There now exists a vast literature of reported "molecular signatures". However, despite some notable exceptions, many of these signatures have suffered from limited reproducibility in independent datasets, insufficient sensitivity or specificity to meet clinical needs, or other challenges. In this paper, we discuss the process of molecular signature discovery on the basis of omics data. In particular, we highlight potential pitfalls in the discovery process, as well as strategies that can be used to increase the odds of successful discovery. Despite the difficulties that have plagued the field of molecular signature discovery, we remain optimistic about the potential to harness the vast amounts of available omics data in order to substantially impact clinical practice.

Subnetwork-based analysis of chronic lymphocytic leukemia identifies pathways that associate with disease progression

The clinical course of patients with chronic lymphocytic leukemia (CLL) is heterogeneous. Several prognostic factors have been identified that can stratify patients into groups that differ in their relative tendency for disease progression and/or survival. Here, we pursued a subnetwork-based analysis of gene expression profiles to discriminate between groups of patients with disparate risks for CLL progression. From an initial cohort of 130 patients, we identified 38 prognostic subnetworks that could predict the relative risk for disease progression requiring therapy from the time of sample collection, more accurately than established markers. The prognostic power of these subnetworks then was validated on 2 other cohorts of patients. We noted reduced divergence in gene expression between leukemia cells of CLL patients classified at diagnosis with aggressive versus indolent disease over time. The predictive subnetworks vary in levels of expression over time but exhibit increased similarity at later time points before therapy, suggesting that degenerate pathways apparently converge into common pathways that are associated with disease progression. As such, these results have implications for understanding cancer evolution and for the development of novel treatment strategies for patients with CLL.

Single sample expression-anchored mechanisms predict survival in head and neck cancer

Gene expression signatures that are predictive of therapeutic response or prognosis are increasingly useful in clinical care; however, mechanistic (and intuitive) interpretation of expression arrays remains an unmet challenge. Additionally, there is surprisingly little gene overlap among distinct clinically validated expression signatures. These “causality challenges” hinder the adoption of signatures as compared to functionally well-characterized single gene biomarkers. To increase the utility of multi-gene signatures in survival studies, we developed a novel approach to generate “personal mechanism signatures” of molecular pathways and functions from gene expression arrays. FAIME, the Functional Analysis of Individual Microarray Expression, computes mechanism scores using rank-weighted gene expression of an individual sample. By comparing head and neck squamous cell carcinoma (HNSCC) samples with non-tumor control tissues, the precision and recall of deregulated FAIME-derived mechanisms of pathways and molecular functions are comparable to those produced by conventional cohort-wide methods (e.g. GSEA). The overlap of “Oncogenic FAIME Features of HNSCC” (statistically significant and differentially regulated FAIME-derived genesets representing GO functions or KEGG pathways derived from HNSCC tissue) among three distinct HNSCC datasets (pathways:46%, p<0.001) is more significant than the gene overlap (genes:4%). These Oncogenic FAIME Features of HNSCC can accurately discriminate tumors from control tissues in two additional HNSCC datasets (n = 35 and 91, F-accuracy = 100% and 97%, empirical p<0.001, area under the receiver operating characteristic curves = 99% and 92%), and stratify recurrence-free survival in patients from two independent studies (p = 0.0018 and p = 0.032, log-rank). Previous approaches depending on group assignment of individual samples before selecting features or learning a classifier are limited by design to discrete-class prediction. In contrast, FAIME calculates mechanism profiles for individual patients without requiring group assignment in validation sets. FAIME is more amenable for clinical deployment since it translates the gene-level measurements of each given sample into pathways and molecular function profiles that can be applied to analyze continuous phenotypes in clinical outcome studies (e.g. survival time, tumor volume).

Clinical utilization of multi-gene expression signatures that are predictive of therapeutic response has been steadily increasing, however, interpretation of such results remains challenging because multi-gene signatures, generated from analyzing different patient cohorts, tend to be equally predictive but contain minimal overlap. Whereas pathway-level analyses of expression arrays show promise for generating clinically meaningful mechanistic signatures, current approaches do not permit single-patient based analyses that are independent of cross-group calculations. To bridge the gap between deterministic biological mechanisms of single-gene biomarkers and the statistical predictive power of multi-gene signatures that are disconnected from mechanisms, we developed FAIME, a novel method that transforms microarray gene expression data into individualized patient profiles of molecular mechanisms. We have validated its capability for predicting clinical outcomes, including cancer patient samples derived from six different clinical trial cohorts of head and neck cancers. This method provides opportunities to harness an untapped resource for personal genomics: clinical evaluation and testing of individually interpretable mechanistic profiles derived from gene expression arrays.

Identification of predictive markers of response to the MEK1/2 inhibitor selumetinib (AZD6244) in K-ras-mutated colorectal cancer

Mutant K-ras activity leads to the activation of the RAS/RAF/MEK/ERK pathway in approximately 44% of colorectal cancer (CRC) tumors. Accordingly, several inhibitors of the MEK pathway are under clinical evaluation in several malignancies including CRC. The aim of this study was to develop and characterize predictive biomarkers of response to the MEK1/2 inhibitor AZD6244 in CRC in order to maximize the clinical utility of this agent. Twenty-seven human CRC cell lines were exposed to AZD6244 and classified according to the IC(50) value as sensitive (≤ 0.1 μmol/L) or resistant (>1 μmol/L). All cell lines were subjected to immunoblotting for effector proteins, K-ras/BRAF mutation status, and baseline gene array analysis. Further testing was done in cell line xenografts and K-ras mutant CRC human explants models to develop a predictive genomic classifier for AZD6244. The most sensitive and resistant cell lines were subjected to differential gene array and pathway analyses. Members of the Wnt signaling pathway were highly overexpressed in cell lines resistant to AZD6244 and seem to be functionally involved in mediating resistance by shRNA knockdown studies. Baseline gene array data from CRC cell lines and xenografts were used to develop a k-top scoring pair (k-TSP) classifier, which predicted with 71% accuracy which of a test set of patient-derived K-ras mutant CRC explants would respond to AZD6244, providing the basis for a patient-selective clinical trial. These results also indicate that resistance to AZD6244 may be mediated, in part, by the upregulation of the Wnt pathway, suggesting potential rational combination partners for AZD6244 in CRC.

The basis and rational use of molecular genetic testing in mature B-cell lymphomas

An increasing number of neoplasms are associated with variably specific genetic abnormalities. This is best exemplified by hematological malignancies, in which there is a growing list of entities that are defined by their genetic lesion(s); this is not (yet) the case in mature B-cell lymphomas. However, enhanced insights into the pathogenesis of this large and diverse group of lymphomas have emerged with the ongoing unraveling of a plethora of fascinating genetic abnormalities. The purpose of this review is to synthesize well-recognized data and nascent discoveries in our understanding of the genetic basis of a spectrum of mature B-cell lymphomas, and how this may be applied to contemporary clinical practice. Despite the explosion of new and exciting knowledge in this arena, with the potential for enhanced diagnostic and prognostic strategies, it is essential to remain cognizant of the limitations (and complexity) of genetic investigations, so that assays can be developed and used both judiciously and rationally.