Individualized identification of disease-associated pathways with disrupted coordination of gene expression

Metabolism-Related Gene Pairs to Predict the Clinical Outcome and Molecular Characteristics of Early Hepatocellular Carcinoma

Simple Summary

After surgery, about 60–70% of early hepatocellular carcinoma patients suffer from relapse within 5 years, hindering long-term survival. Clinical and pathologic features cannot provide an accurate evaluation. We aimed to construct a stratification model from the metabolic aspect to predict the clinical outcome and reveal the molecular characteristics of different prognostic subgroups. An individualized metabolic signature of 10 gene pairs was developed from 250 early HCCs and validated in 311 samples from different datasets. The signature stratified early HCC cases one-by-one into two risk groups with different survival rates. The molecular characteristics of the two risk groups were analyzed by multi-omics data. The relationships with proliferation, immunity, and drug benefits were summarized. The signature was further validated in 47 institutional transcriptomic HCC samples and 101 public proteomic samples.

Abstract

Recurrence is the main factor affecting the prognosis of early hepatocellular carcinoma (HCC), which is not accurately evaluated by clinical indicators. The metabolic heterogeneity of HCC hints at the possibility of constructing a stratification model to predict the clinical outcome. On the basis of the relative expression orderings of 2939 metabolism-related genes, an individualized signature with 10 metabolism-related gene pairs (10-GPS) was developed from 250 early HCC samples in the discovery datasets, which stratified HCC patients into the high- and low-risk subgroups with significantly different survival rates. The 10-GPS was validated in 311 public transcriptomic samples from two independent validation datasets. A nomogram that included the 10-GPS, age, gender, and stage was constructed for eventual clinical evaluation. The low-risk group was characterized by lower proliferation, higher metabolism, increased activated immune microenvironment, and lower TIDE scores, suggesting a better response to immunotherapy. The high-risk group displayed hypomethylation, higher copy number alterations, mutations, and more overexpression of immune-checkpoint genes, which might jointly lead to poor outcomes. The prognostic accuracy of the 10-GPS was further validated in 47 institutional transcriptomic samples and 101 public proteomic samples. In conclusion, the 10-GPS is a robust predictor of the clinical outcome for early HCC patients and could help evaluate prognosis and characterize molecular heterogeneity.

Individualized pathway activity algorithm identifies oncogenic pathways in pan-cancer analysis

Background

Accumulative evidences have shown that dysregulation of biological pathways contributed to the initiation and progression of malignant tumours. Several methods for pathway activity measurement have been proposed, but they are restricted to making comparisons between groups or sensitive to experimental batch effects.

Methods

We introduced a novel method for individualized pathway activity measurement (IPAM) that is based on the ranking of gene expression levels in individual sample. Taking advantage of IPAM, we calculated the pathway activity of 318 pathways from KEGG database in the 10528 tumour/normal samples of 33 cancer types from TCGA to identify characteristic dysregulated pathways among different cancer types.

Findings

IPAM precisely quantified the level of activity of each pathway in pan-cancer analysis and exhibited better performance in cancer classification and prognosis prediction over five widely used tools. The average ROC-AUC of cancer diagnostic model using tumour-educated platelets (TEPs) reached 92.84%, suggesting the potential of our algorithm in early diagnosis of cancer. We identified several pathways significantly deregulated and associated with patient survival in a large fraction of cancer types, such as tyrosine metabolism, fatty acid degradation, cell cycle, p53 signalling pathway and DNA replication. We also confirmed the dominant role of metabolic pathways in cancer pathway dysregulation and identified the driving factors of specific pathway dysregulation, such as PPARA for branched-chain amino acid metabolism and NR1I2, NR1I3 for fatty acid metabolism.

Interpretation

Our study will provide novel clues for understanding the pathological mechanisms of cancer, ultimately paving the way for personalized medicine of cancer.

Funding

A full list of funding can be found in the Acknowledgements section.

SEPA: Signalling entropy-based algorithm to evaluate personalized pathway activation for survival analysis on pan-cancer data

MOTIVATION

Biomarkers with prognostic ability and biological interpretability can be used to support decision-making in the survival analysis. Genes usually form functional modules to play synergistic roles, such as pathways. Predicting significant features from the functional level can effectively reduce the adverse effects of heterogeneity and obtain more reproducible and interpretable biomarkers. Personalized pathway activation inference can quantify the dysregulation of essential pathways involved in the initiation and progression of cancers, and can contribute to the development of personalized medical treatments.

RESULTS

In this study, we propose a novel method to evaluate personalized pathway activation based on signalling entropy for survival analysis (SEPA), which is a new attempt to introduce the information-theoretic entropy in generating pathway representation for each patient. SEPA effectively integrates pathway-level information into gene expression data, converting the high-dimensional gene expression data into the low-dimensional biological pathway activation scores. SEPA shows its classification power on the prognostic pan-cancer genomic data, and the potential pathway markers identified based on SEPA have statistical significance in the discrimination of high-risk and low-risk cohorts and are likely to be associated with the initiation and progress of cancers. The results show that SEPA scores can be used as an indicator to precisely distinguish cancer patients with different clinical outcomes, and identify important pathway features with strong discriminative power and biological interpretability.

AVAILABILITY

The MATLAB-package for SEPA is freely available from https://github.com/xingyili/SEPA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

IndGOterm: a qualitative method for the identification of individually dysregulated GO terms in cancer

Individual pathway analysis can dissect heterogeneities among different cancer patients and provide efficient guidelines for individualized therapy. However, the existence of the batch effect brings extensive limitations for the application of many individual methods for pathway analysis. Previously, researchers proposed that methods based on within-sample relative expression ordering (REO) of the genes are notably insensitive to ‘batch effects’. In this article, we focus on the Gene Ontology (GO) database and propose an individual qualitative GO term analysis method (IndGOterm) based on the REO of genes. Compared with some current widely used single-sample enrichment analysis methods, such as ssGSEA and GSVA, IndGOterm has a predominance of ignoring the batch effects caused by diverse technologies. Through the survival and drug responses analysis, we found IndGOterm could capture more terms connected to cancer than other single-sample enrichment analysis methods. Furthermore, through the application of IndGOterm, we found some terms that present different dysregulation models that manifest heterogenetic in homologous patients. Collectively, these results attested that IndGOterm could capture useful information from patients and be a useful tool to reveal the intrinsic characteristic of cancer. An open-source R statistical analysis package ‘IndGOterm’ is available at https://github.com/robert19960424/IndGOterm.

Identifying individualized risk subpathways reveals pan-cancer molecular classification based on multi-omics data

Cancer is a highly heterogeneous disease with different functional disorders among individuals. The initiation and progression of cancer is usually related to dysregulation of local regions within pathways. Identification of individualized risk pathways is crucial for revealing the mechanisms of tumorigenesis and heterogeneity. However, approach that focused on mining patient-specific risk subpathway regions is still lacking. Here, we developed an individualized cancer risk subpathway identification method that was referred as InCRiS by integrating multi-omics data. Then, the method was applied to nearly 3000 samples across 9 TCGA cancer types and its robustness and reliability were evaluated. Dissecting dysregulated subpathways in these tumor samples revealed several key regions which may play oncogenic roles in cancer. The construction of risk subpathway dysregulation profile of pan-cancers revealed 11 pan-cancer molecular classification (InCRiS subtypes) with significantly different clinical outcomes. Moreover, subpathway regions that tend to be disordered in individuals of specific subtypes were examined for understanding the pathogenesis of tumor and some key genes such as CTNNB1, EP300 and PRKDC were nominated in different subtypes. In summary, the proposed method and resulting data presented useful resources to study the mechanism of tumor heterogeneity and to discovery novel therapeutic targets for precise treatment of cancer.

Highly accurate diagnosis of papillary thyroid carcinomas based on personalized pathways coupled with machine learning

Thyroid nodules are neoplasms commonly found among adults, with papillary thyroid carcinoma (PTC) being the most prevalent malignancy. However, current diagnostic methods often subject patients to unnecessary surgical burden. In this study, we developed and validated an automated, highly accurate multi-study-derived diagnostic model for PTCs using personalized biological pathways coupled with a sophisticated machine learning algorithm. Surprisingly, the algorithm achieved near-perfect performance in discriminating PTCs from non-tumoral thyroid samples with an overall cross-study-validated area under the receiver operating characteristic curve (AUROC) of 0.999 (95% confidence interval [CI]: 0.995-1) and a Brier score of 0.013 on three independent development cohorts. In addition, the algorithm showed excellent generalizability and transferability on two large-scale external blind PTC cohorts consisting of The Cancer Genome Atlas (TCGA), which is the largest genomic PTC cohort studied to date, and the post-Chernobyl cohort, which includes PTCs reported after exposure to radiation from the Chernobyl accident. When applied to the TCGA cohort, the model yielded an AUROC of 0.969 (95% CI: 0.950-0.987) and a Brier score of 0.109. On the post-Chernobyl cohort, it yielded an AUROC of 0.962 (95% CI: 0.918-1) and a Brier score of 0.073. This algorithm also is robust against other various types of clinical scenarios, discriminating malignant from benign lesions as well as clinically aggressive thyroid cancer with poor prognosis from indolent ones. Furthermore, we discovered novel pathway alterations and prognostic signatures for PTC, which can provide directions for follow-up studies.

GMIEC: a shiny application for the identification of gene-targeted drugs for precision medicine

Abstract

Background

Precision medicine is a medical approach that takes into account individual genetic variability and often requires Next Generation Sequencing data in order to predict new treatments. Here we present GMIEC, Genomic Modules Identification et Characterization for genomics medicine, an application that is able to identify specific drugs at the level of single patient integrating multi-omics data such as RNA-sequencing, copy-number variation, methylation, Chromatin Immuno-Precipitation and Exome/Whole Genome sequencing. It is also possible to include clinical data related to each patient. GMIEC has been developed as a web-based R-Shiny platform and gives as output a table easy to use and explore.

Results

We present GMIEC, a Shiny application for genomics medicine. The tool allows the users the integration of two or more multiple omics datasets (e.g. gene-expression, copy-number), at sample level, to identify groups of genes that share common genomic and corresponding drugs. We demonstrate the characteristics of our application by using it to analyze a prostate cancer data set.

Conclusions

GMIEC provides a simple interface for genomics medicine. GMIEC was develop with Shiny to provide an application that does not require advanced programming skills. GMIEC consists of three sub-application for the analysis (GMIEC-AN), the visualization (GMIEC-VIS) and the exploration of results (GMIEC-RES). GMIEC is an open source software and is available at https://github.com/guidmt/GMIEC-shiny

Developing a 'personalome' for precision medicine: emerging methods that compute interpretable effect sizes from single-subject transcriptomes

The development of computational methods capable of analyzing -omics data at the individual level is critical for the success of precision medicine. Although unprecedented opportunities now exist to gather data on an individual's -omics profile ('personalome'), interpreting and extracting meaningful information from single-subject -omics remain underdeveloped, particularly for quantitative non-sequence measurements, including complete transcriptome or proteome expression and metabolite abundance. Conventional bioinformatics approaches have largely been designed for making population-level inferences about 'average' disease processes; thus, they may not adequately capture and describe individual variability. Novel approaches intended to exploit a variety of -omics data are required for identifying individualized signals for meaningful interpretation. In this review-intended for biomedical researchers, computational biologists and bioinformaticians-we survey emerging computational and translational informatics methods capable of constructing a single subject's 'personalome' for predicting clinical outcomes or therapeutic responses, with an emphasis on methods that provide interpretable readouts. Key points: (i) the single-subject analytics of the transcriptome shows the greatest development to date and, (ii) the methods were all validated in simulations, cross-validations or independent retrospective data sets. This survey uncovers a growing field that offers numerous opportunities for the development of novel validation methods and opens the door for future studies focusing on the interpretation of comprehensive 'personalomes' through the integration of multiple -omics, providing valuable insights into individual patient outcomes and treatments.

Impact of Colon-Specific DNA Methylation-Regulated Gene Modules on Colorectal Cancer Patient Survival

Background

Colorectal cancer (CRC) is the second most commonly diagnosed cancer in females and the third in males worldwide. Although existing evidence explained some critical functions of the single genetic abnormality in the pathogenesis of CRC, the function of interactors involved in the colon-specific regulatory network, especially DNA methylation regulated network is still poorly understood.

Material/Methods

In this work, matched gene expression and DNA methylation samples of CRC patients were retrieved. Differential gene expression and methylation analyses were performed. In addition, gene expression and DNA methylation were integrated into a colon-specific regulatory gene network, detecting the epigenetically regulated gene modules which drive CRC through an underlying epigenetic mechanism. Finally, the colon-specific DNA methylation-regulated gene modules were validated using an independent set of CRC patients.

Results

Differential gene expression analysis demonstrated the upregulation of the cell cycle and DNA replication and downregulation of cGMP-PKG signaling pathway and calcium signaling pathway in CRC. Differentially methylated regions (DMRs) showed the different levels of methylation in promoters, CpG islands, and genes in CRC. In addition, gene expression and DNA methylation were integrated into a colon-specific regulatory gene network, detecting 8 epigenetically regulated gene modules which drive CRC through an underlying epigenetic mechanism. Interestingly, 2 of the colon-specific DNA methylation-regulated gene modules showed a significant predictive ability for the survival of an independent set of CRC patients.

Conclusions

The results of this study could open a new era and aid the development of novel therapeutic targets for the treatment of CRC patient.

Comprehensive and critical evaluation of individualized pathway activity measurement tools on pan-cancer data

Motivation

: Biological pathways are extensively used for the analysis of transcriptome data to characterize biological mechanisms underlying various phenotypes. There are a number of computational tools that summarize transcriptome data at the pathway level. However, there is no comparative study on how well these tools produce useful information at the cohort level, enabling comparison of many samples or patients.

Results

: In this study, we systematically compared and evaluated 13 different pathway activity inference tools based on 5 comparison criteria using pan-cancer data set. This study has two major contributions. First, our study provides a comprehensive survey on computational techniques used by existing pathway activity inference tools. The tools use different strategies and assume different requirements on data: input transformation, use of labels, necessity of cohort-level input data, use of gene relations and scoring metric. Second, we performed extensive evaluations on the performance of these tools. Because different tools use different methods to map samples to the pathway dimension, the tools are evaluated at the pathway level using five comparison criteria. Starting from measuring how well a tool maintains the characteristics of original gene expression values, robustness was also investigated by adding noise into gene expression data. Classification tasks on three clinical variables (tumor versus normal, survival and cancer subtypes) were performed to evaluate the utility of tools for their clinical applications. In addition, the inferred activity values were compared between the tools to see how similar they are along with the scoring schemes they use.

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.

Evaluating hepatocellular carcinoma cell lines for tumour samples using within-sample relative expression orderings of genes

BACKGROUND & AIMS

Concerns are raised about the representativeness of cell lines for tumours due to the culture environment and misidentification. Liver is a major metastatic destination of many cancers, which might further confuse the origin of hepatocellular carcinoma cell lines. Therefore, it is of crucial importance to understand how well they can represent hepatocellular carcinoma.

METHODS

The HCC-specific gene pairs with highly stable relative expression orderings in more than 99% of hepatocellular carcinoma but with reversed relative expression orderings in at least 99% of one of the six types of cancer, colorectal carcinoma, breast carcinoma, non-small-cell lung cancer, gastric carcinoma, pancreatic carcinoma and ovarian carcinoma, were identified.

RESULTS

With the simple majority rule, the HCC-specific relative expression orderings from comparisons with colorectal carcinoma and breast carcinoma could exactly discriminate primary hepatocellular carcinoma samples from both primary colorectal carcinoma and breast carcinoma samples. Especially, they correctly classified more than 90% of liver metastatic samples from colorectal carcinoma and breast carcinoma to their original tumours. Finally, using these HCC-specific relative expression orderings from comparisons with six cancer types, we identified eight of 24 hepatocellular carcinoma cell lines in the Cancer Cell Line Encyclopedia (Huh-7, Huh-1, HepG2, Hep3B, JHH-5, JHH-7, C3A and Alexander cells) that are highly representative of hepatocellular carcinoma. Evaluated with a REOs-based prognostic signature for hepatocellular carcinoma, all these eight cell lines showed the same metastatic properties of the high-risk metastatic hepatocellular carcinoma tissues.

CONCLUSIONS

Caution should be taken for using hepatocellular carcinoma cell lines. Our results should be helpful to select proper hepatocellular carcinoma cell lines for biological experiments.

Differential expression analysis at the individual level reveals a lncRNA prognostic signature for lung adenocarcinoma

BACKGROUND

Deregulations of long non-coding RNAs (lncRNAs) have been implicated in cancer initiation and progression. Current methods can only capture differential expression of lncRNAs at the population level and ignore the heterogeneous expression of lncRNAs in individual patients.

METHODS

We propose a method (LncRIndiv) to identify differentially expressed (DE) lncRNAs in individual cancer patients by exploiting the disrupted ordering of expression levels of lncRNAs in each disease sample in comparison with stable normal ordering. LncRIndiv was applied to lncRNA expression profiles of lung adenocarcinoma (LUAD). Based on the expression profile of LUAD individual-level DE lncRNAs, we used a forward selection procedure to identify prognostic signature for stage I-II LUAD patients without adjuvant therapy.

RESULTS

In both simulated data and real pair-wise cancer and normal sample data, LncRIndiv method showed good performance. Based on the individual-level DE lncRNAs, we developed a robust prognostic signature consisting of two lncRNA (C1orf132 and TMPO-AS1) for stage I-II LUAD patients without adjuvant therapy (P = 3.06 × 10^-6, log-rank test), which was confirmed in two independent datasets of GSE50081 (P = 1.82 × 10^-2, log-rank test) and GSE31210 (P = 7.43 × 10^-4, log-rank test) after adjusting other clinical factors such as smoking status and stages. Pathway analysis showed that TMPO-AS1 and C1orf132 could affect the prognosis of LUAD patients through regulating cell cycle and cell adhesion.

CONCLUSIONS

LncRIndiv can successfully detect DE lncRNAs in individuals and be applied to identify prognostic signature for LUAD patients.

Testing for differentially expressed genetic pathways with single-subject N-of-1 data in the presence of inter-gene correlation

Modern precision medicine increasingly relies on molecular data analytics, wherein development of interpretable single-subject ("N-of-1") signals is a challenging goal. A previously developed global framework, N-of-1- pathways, employs single-subject gene expression data to identify differentially expressed gene set pathways in an individual patient. Unfortunately, the limited amount of data within the single-subject, N-of-1 setting makes construction of suitable statistical inferences for identifying differentially expressed gene set pathways difficult, especially when non-trivial inter-gene correlation is present. We propose a method that exploits external information on gene expression correlations to cluster positively co-expressed genes within pathways, then assesses differential expression across the clusters within a pathway. A simulation study illustrates that the cluster-based approach exhibits satisfactory false-positive error control and reasonable power to detect differentially expressed gene set pathways. An example with a single N-of-1 patient's triple negative breast cancer data illustrates use of the methodology.

Identifying disease-associated pathways in one-phenotype data based on reversal gene expression orderings

Due to the invasiveness nature of tissue biopsy, it is common that investigators cannot collect sufficient normal controls for comparison with diseased samples. We developed a pathway enrichment tool, DRFunc, to detect significantly disease-disrupted pathways by incorporating normal controls from other experiments. The method was validated using both microarray and RNA-seq expression data for different cancers. The high concordant differentially ranked (DR) gene pairs were identified between cases and controls from different independent datasets. The DR gene pairs were used in the DRFunc algorithm to detect significantly disrupted pathways in one-phenotype expression data by combing controls from other studies. The DRFunc algorithm was exemplified by the detection of significant pathways in glioblastoma samples. The algorithm can also be used to detect altered pathways in the datasets with weak expression signals, as shown by the analysis on the expression data of chemotherapy-treated breast cancer samples.

Individualized analysis reveals CpG sites with methylation aberrations in almost all lung adenocarcinoma tissues

BACKGROUND

Due to the heterogeneity of cancer, identifying differentially methylated (DM) CpG sites between a set of cancer samples and a set of normal samples cannot tell us which patients have methylation aberrations in a particular DM CpG site.

METHODS

We firstly showed that the relative methylation-level orderings (RMOs) of CpG sites within individual normal lung tissues are highly stable but widely disrupted in lung adenocarcinoma tissues. This finding provides the basis of using the RankComp algorithm, previously developed for differential gene expression analysis at the individual level, to identify DM CpG sites in each cancer tissue compared with its own normal state. Briefly, through comparing with the highly stable normal RMOs predetermined in a large collection of samples for normal lung tissues, the algorithm finds those CpG sites whose hyper- or hypo-methylations may lead to the disrupted RMOs of CpG site pairs within a disease sample based on Fisher's exact test.

RESULTS

Evaluated in 59 lung adenocarcinoma tissues with paired adjacent normal tissues, RankComp reached an average precision of 94.26% for individual-level DM CpG sites. Then, after identifying DM CpG sites in each of the 539 lung adenocarcinoma samples from TCGA, we found five and 44 CpG sites hypermethylated and hypomethylated in above 90% of the disease samples, respectively. These findings were validated in 140 publicly available and eight additionally measured paired cancer-normal samples. Gene expression analysis revealed that four of the five genes, HOXA9, TAL1, ATP8A2, ENG and SPARCL1, each harboring one of the five frequently hypermethylated CpG sites within its promoters, were also frequently down-regulated in lung adenocarcinoma.

CONCLUSIONS

The common DNA methylation aberrations in lung adenocarcinoma tissues may be important for lung adenocarcinoma diagnosis and therapy.

Differential expression analysis for individual cancer samples based on robust within-sample relative gene expression orderings across multiple profiling platforms

The highly stable within-sample relative expression orderings (REOs) of gene pairs in a particular type of human normal tissue are widely reversed in the cancer condition. Based on this finding, we have recently proposed an algorithm named RankComp to detect differentially expressed genes (DEGs) for individual disease samples measured by a particular platform. In this paper, with 461 normal lung tissue samples separately measured by four commonly used platforms, we demonstrated that tens of millions of gene pairs with significantly stable REOs in normal lung tissue can be consistently detected in samples measured by different platforms. However, about 20% of stable REOs commonly detected by two different platforms (e.g., Affymetrix and Illumina platforms) showed inconsistent REO patterns due to the differences in probe design principles. Based on the significantly stable REOs (FDR<0.01) for normal lung tissue consistently detected by the four platforms, which tended to have large rank differences, RankComp detected averagely 1184, 1335 and 1116 DEGs per sample with averagely 96.51%, 95.95% and 94.78% precisions in three evaluation datasets with 25, 57 and 58 paired lung cancer and normal samples, respectively. Individualized pathway analysis revealed some common and subtype-specific functional mechanisms of lung cancer. Similar results were observed for colorectal cancer. In conclusion, based on the cross-platform significantly stable REOs for a particular normal tissue, differentially expressed genes and pathways in any disease sample measured by any of the platforms can be readily and accurately detected, which could be further exploited for dissecting the heterogeneity of cancer.

Identification of differentially expressed miRNAs in individual breast cancer patient and application in personalized medicine

MicroRNAs (miRNAs) have key roles in breast cancer progression, and their expression levels are heterogeneous across individual breast cancer patients. Traditional methods aim to identify differentially expressed miRNAs in populations rather than in individuals and are affected by the expression intensities of miRNAs in different experimental batches or platforms. Thus it is urgent to conduct miRNA differential expression analysis at an individual level for further personalized medicine research. We proposed a straightforward method to determine the differential expression of each miRNA in an individual patient by utilizing the reversal expression order of miRNA pairs between two conditions (cancer and normal tissue). We applied our method to breast cancer miRNA expression profiles from The Cancer Genome Atlas and two other independent data sets. In total, 292 miRNAs were differentially expressed in individual breast cancer patients. Using the differential expression profile of miRNAs in individual patients, we found that the deregulations of miRNA tend to occur in specific breast cancer subtypes. We investigated the coordination effect between the miRNA and its target, based on the hypothesis that one gene function can be changed by copy number alterations of the corresponding gene or deregulation of the miRNA. We revealed that patients exhibiting an upregulation of hsa-miR-92b and patients with deletions of PTEN did not tend to overlap, and hsa-miR-92b and PTEN coordinately regulated the pathway of 'cell cycle' and so on. Moreover, we discovered a new prognostic signature, hsa-miR-29c, whose downregulation was associated with poor survival of breast cancer patients.

Personalised pathway analysis reveals association between DNA repair pathway dysregulation and chromosomal instability in sporadic breast cancer

The Homologous Recombination (HR) pathway is crucial for the repair of DNA double-strand breaks (DSBs) generated during DNA replication. Defects in HR repair have been linked to the initiation and development of a wide variety of human malignancies, and exploited in chemical, radiological and targeted therapies. In this study, we performed a personalised pathway analysis independently for four large sporadic breast cancer cohorts to investigate the status of HR pathway dysregulation in individual sporadic breast tumours, its association with HR repair deficiency and its impact on tumour characteristics. Specifically, we first manually curated a list of HR genes according to our recent review on this pathway (Liu et al., 2014), and then applied a personalised pathway analysis method named Pathifier (Drier et al., 2013) on the expression levels of the curated genes to obtain an HR score quantifying HR pathway dysregulation in individual tumours. Based on the score, we observed a great diversity in HR dysregulation between and within gene expression-based breast cancer subtypes, and by using two published HR-defect signatures, we found HR pathway dysregulation reflects HR repair deficiency. Furthermore, we identified a novel association between HR pathway dysregulation and chromosomal instability (CIN) in sporadic breast cancer. Although CIN has long been considered as a hallmark of most solid tumours, with recent extensive studies highlighting its importance in tumour evolution and drug resistance, the molecular basis of CIN in sporadic cancers remains poorly understood. Our results imply that HR pathway dysregulation might contribute to CIN in sporadic breast cancer.