Lessons from a decade of integrating cancer copy number alterations with gene expression profiles

Fundamental immune-oncogenicity trade-offs define driver mutation fitness

Missense driver mutations in cancer are concentrated in a few hotspots1. Various mechanisms have been proposed to explain this skew, including biased mutational processes2, phenotypic differences3-6 and immunoediting of neoantigens7,8; however, to our knowledge, no existing model weighs the relative contribution of these features to tumour evolution. We propose a unified theoretical 'free fitness' framework that parsimoniously integrates multimodal genomic, epigenetic, transcriptomic and proteomic data into a biophysical model of the rate-limiting processes underlying the fitness advantage conferred on cancer cells by driver gene mutations. Focusing on TP53, the most mutated gene in cancer1, we present an inference of mutant p53 concentration and demonstrate that TP53 hotspot mutations optimally solve an evolutionary trade-off between oncogenic potential and neoantigen immunogenicity. Our model anticipates patient survival in The Cancer Genome Atlas and patients with lung cancer treated with immunotherapy as well as the age of tumour onset in germline carriers of TP53 variants. The predicted differential immunogenicity between hotspot mutations was validated experimentally in patients with cancer and in a unique large dataset of healthy individuals. Our data indicate that immune selective pressure on TP53 mutations has a smaller role in non-cancerous lesions than in tumours, suggesting that targeted immunotherapy may offer an early prophylactic opportunity for the former. Determining the relative contribution of immunogenicity and oncogenic function to the selective advantage of hotspot mutations thus has important implications for both precision immunotherapies and our understanding of tumour evolution.

Obscurin: A multitasking giant in the fight against cancer

Giant obscurins (720-870 kDa), encoded by OBSCN, were originally discovered in striated muscles as cytoskeletal proteins with scaffolding and regulatory roles. Recently though, they have risen to the spotlight as key players in cancer development and progression. Herein, we provide a timely prudent synopsis of the expanse of OBSCN mutations across 16 cancer types. Given the extensive work on OBSCN's role in breast epithelium, we summarize functional studies implicating obscurins as potent tumor suppressors in breast cancer and delve into an in silico analysis of its mutational profile and epigenetic (de)regulation using different dataset platforms and sophisticated computational tools. Lastly, we formally describe the OBSCN-Antisense-RNA-1 gene, which belongs to the long non-coding RNA family and discuss its potential role in modulating OBSCN expression in breast cancer. Collectively, we highlight the escalating involvement of obscurins in cancer biology and outline novel potential mechanisms of OBSCN (de)regulation that warrant further investigation.

A Bayesian two-way latent structure model for genomic data integration reveals few pan-genomic cluster subtypes in a breast cancer cohort

MOTIVATION

Unsupervised clustering is important in disease subtyping, among having other genomic applications. As genomic data has become more multifaceted, how to cluster across data sources for more precise subtyping is an ever more important area of research. Many of the methods proposed so far, including iCluster and Cluster of Cluster Assignments (COCAs), make an unreasonable assumption of a common clustering across all data sources, and those that do not are fewer and tend to be computationally intensive.

RESULTS

We propose a Bayesian parametric model for integrative, unsupervised clustering across data sources. In our two-way latent structure model, samples are clustered in relation to each specific data source, distinguishing it from methods like COCAs and iCluster, but cluster labels have across-dataset meaning, allowing cluster information to be shared between data sources. A common scaling across data sources is not required, and inference is obtained by a Gibbs Sampler, which we improve with a warm start strategy and modified density functions to robustify and speed convergence. Posterior interpretation allows for inference on common clusterings occurring among subsets of data sources. An interesting statistical formulation of the model results in sampling from closed-form posteriors despite incorporation of a complex latent structure. We fit the model with Gaussian and more general densities, which influences the degree of across-dataset cluster label sharing. Uniquely among integrative clustering models, our formulation makes no nestedness assumptions of samples across data sources so that a sample missing data from one genomic source can be clustered according to its existing data sources. We apply our model to a Norwegian breast cancer cohort of ductal carcinoma in situ and invasive tumors, comprised of somatic copy-number alteration, methylation and expression datasets. We find enrichment in the Her2 subtype and ductal carcinoma among those observations exhibiting greater cluster correspondence across expression and CNA data. In general, there are few pan-genomic clusterings, suggesting that models assuming a common clustering across genomic data sources might yield misleading results.

AVAILABILITY AND IMPLEMENTATION

The model is implemented in an R package called twl ('two-way latent'), available on CRAN. Data for analysis are available within the R package.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A meta-analysis of multiple matched aCGH/expression cancer datasets reveals regulatory relationships and pathway enrichment of potential oncogenes

The copy numbers of genes in cancer samples are often highly disrupted and form a natural amplification/deletion experiment encompassing multiple genes. Matched array comparative genomics and transcriptomics datasets from such samples can be used to predict inter-chromosomal gene regulatory relationships. Previously we published the database METAMATCHED, comprising the results from such an analysis of a large number of publically available cancer datasets. Here we investigate genes in the database which are unusual in that their copy number exhibits consistent heterogeneous disruption in a high proportion of the cancer datasets. We assess the potential relevance of these genes to the pathology of the cancer samples, in light of their predicted regulatory relationships and enriched biological pathways. A network-based method was used to identify enriched pathways from the genes’ inferred targets. The analysis predicts both known and new regulator-target interactions and pathway memberships. We examine examples in detail, in particular the gene POGZ, which is disrupted in many of the cancer datasets and has an unusually large number of predicted targets, from which the network analysis predicts membership of cancer related pathways. The results suggest close involvement in known cancer pathways of genes exhibiting consistent heterogeneous copy number disruption. Further experimental work would clarify their relevance to tumor biology. The results of the analysis presented in the database METAMATCHED, and included here as an R archive file, constitute a large number of predicted regulatory relationships and pathway memberships which we anticipate will be useful in informing such experiments.

MVisAGe Identifies Concordant and Discordant Genomic Alterations of Driver Genes in Squamous Tumors

Integrated analyses of multiple genomic datatypes are now common in cancer profiling studies. Such data present opportunities for numerous computational experiments, yet analytic pipelines are limited. Tools such as the cBioPortal and Regulome Explorer, although useful, are not easy to access programmatically or to implement locally. Here, we introduce the MVisAGe R package, which allows users to quantify gene-level associations between two genomic datatypes to investigate the effect of genomic alterations (e.g., DNA copy number changes on gene expression). Visualizing Pearson/Spearman correlation coefficients according to the genomic positions of the underlying genes provides a powerful yet novel tool for conducting exploratory analyses. We demonstrate its utility by analyzing three publicly available cancer datasets. Our approach highlights canonical oncogenes in chr11q13 that displayed the strongest associations between expression and copy number, including CCND1 and CTTN, genes not identified by copy number analysis in the primary reports. We demonstrate highly concordant usage of shared oncogenes on chr3q, yet strikingly diverse oncogene usage on chr11q as a function of HPV infection status. Regions of chr19 that display remarkable associations between methylation and gene expression were identified, as were previously unreported miRNA-gene expression associations that may contribute to the epithelial-to-mesenchymal transition.Significance: This study presents an important bioinformatics tool that will enable integrated analyses of multiple genomic datatypes. Cancer Res; 78(12); 3375-85. ©2018 AACR.

Differentially expressed microRNAs in lung adenocarcinoma invert effects of copy number aberrations of prognostic genes

In many cancers, significantly down- or upregulated genes are found within chromosomal regions with DNA copy number alteration opposite to the expression changes. Generally, this paradox has been overlooked as noise, but can potentially be a consequence of interference of epigenetic regulatory mechanisms, including microRNA-mediated control of mRNA levels.

To explore potential associations between microRNAs and paradoxes in non-small-cell lung cancer (NSCLC) we curated and analyzed lung adenocarcinoma (LUAD) data, comprising gene expressions, copy number aberrations (CNAs) and microRNA expressions. We integrated data from 1,062 tumor samples and 241 normal lung samples, including newly-generated array comparative genomic hybridization (aCGH) data from 63 LUAD samples.

We identified 85 “paradoxical” genes whose differential expression consistently contrasted with aberrations of their copy numbers. Paradoxical status of 70 out of 85 genes was validated on sample-wise basis using The Cancer Genome Atlas (TCGA) LUAD data. Of these, 41 genes are prognostic and form a clinically relevant signature, which we validated on three independent datasets. By meta-analysis of results from 9 LUAD microRNA expression studies we identified 24 consistently-deregulated microRNAs. Using TCGA-LUAD data we showed that deregulation of 19 of these microRNAs explains differential expression of the paradoxical genes.

Our results show that deregulation of paradoxical genes is crucial in LUAD and their expression pattern is maintained epigenetically, defying gene copy number status.

MicroRNA Networks in Breast Cancer Cells

A variety of molecular techniques can be used in order to unravel the molecular composition of cells. In particular, the microarray technology has been used to identify novel biomarkers that may be useful in the diagnosis, prognosis, or treatment of cancer. The microarray technology is ideal for biomarker discovery as it allows for the screening of a large number of molecules at once. In this review, we focus on microRNAs (miRNAs) which are key molecules in cells and regulate gene expression post-transcriptionally. miRNAs are small, single-stranded RNA molecules that bind to complementary mRNAs. Binding of miRNAs to mRNAs leads either to degradation, or translational inhibition of the target mRNA. Roughly one third of all the mRNAs are postulated to be regulated by miRNAs. miRNAs are known to be deregulated in different types of cancer, including breast cancer, and it has been demonstrated that deregulation of several miRNAs can be used as biological markers in cancer. miRNA expression can for example discriminate between normal, benign and malignant breast tissue, and between different breast cancer subtypes.In the post-genomic era, an important task of molecular biology is to understand gene regulation in the context of biological networks. Because miRNAs have such a pronounced role in cells, it is pivotal to understand the mechanisms that underlie their control, and to identify how miRNAs influence cancer development and progression.

Driver pattern identification over the gene co-expression of drug response in ovarian cancer by integrating high throughput genomics data

Multiple types of high throughput genomics data create a potential opportunity to identify driver patterns in ovarian cancer, which will acquire some novel and clinical biomarkers for appropriate diagnosis and treatment to cancer patients. To identify candidate driver genes and the corresponding driving patterns for resistant and sensitive tumors from the heterogeneous data, we combined gene co-expression modules with mutation modulators and proposed the method to identify driver patterns. Firstly, co-expression network analysis is applied to explore gene modules for gene expression profiles through weighted correlation network analysis (WGCNA). Secondly, mutation matrix is generated by integrating the CNV data and somatic mutation data, and a mutation network is constructed from the mutation matrix. Thirdly, candidate modulators are selected from significant genes by clustering vertexs of the mutation network. Finally, a regression tree model is utilized for module network learning, in which the obtained gene modules and candidate modulators are trained for the driving pattern identification and modulators regulatory exploration. Many identified candidate modulators are known to be involved in biological meaningful processes associated with ovarian cancer, such as CCL11, CCL16, CCL18, CCL23, CCL8, CCL5, APOB, BRCA1, SLC18A1, FGF22, GADD45B, GNA15, GNA11, and so on.

Driver pattern identification over the gene co-expression of drug response in ovarian cancer by integrating high throughput genomics data

Multiple types of high throughput genomics data create a potential opportunity to identify driver patterns in ovarian cancer, which will acquire some novel and clinical biomarkers for appropriate diagnosis and treatment to cancer patients. To identify candidate driver genes and the corresponding driving patterns for resistant and sensitive tumors from the heterogeneous data, we combined gene co-expression modules with mutation modulators and proposed the method to identify driver patterns. Firstly, co-expression network analysis is applied to explore gene modules for gene expression profiles through weighted correlation network analysis (WGCNA). Secondly, mutation matrix is generated by integrating the CNV data and somatic mutation data, and a mutation network is constructed from the mutation matrix. Thirdly, candidate modulators are selected from significant genes by clustering vertexs of the mutation network. Finally, a regression tree model is utilized for module network learning, in which the obtained gene modules and candidate modulators are trained for the driving pattern identification and modulators regulatory exploration. Many identified candidate modulators are known to be involved in biological meaningful processes associated with ovarian cancer, such as CCL11, CCL16, CCL18, CCL23, CCL8, CCL5, APOB, BRCA1, SLC18A1, FGF22, GADD45B, GNA15, GNA11, and so on.

Driver pattern identification over the gene co-expression of drug response in ovarian cancer by integrating high throughput genomics data

Multiple types of high throughput genomics data create a potential opportunity to identify driver patterns in ovarian cancer, which will acquire some novel and clinical biomarkers for appropriate diagnosis and treatment to cancer patients. To identify candidate driver genes and the corresponding driving patterns for resistant and sensitive tumors from the heterogeneous data, we combined gene co-expression modules with mutation modulators and proposed the method to identify driver patterns. Firstly, co-expression network analysis is applied to explore gene modules for gene expression profiles through weighted correlation network analysis (WGCNA). Secondly, mutation matrix is generated by integrating the CNV data and somatic mutation data, and a mutation network is constructed from the mutation matrix. Thirdly, candidate modulators are selected from significant genes by clustering vertexs of the mutation network. Finally, a regression tree model is utilized for module network learning, in which the obtained gene modules and candidate modulators are trained for the driving pattern identification and modulators regulatory exploration. Many identified candidate modulators are known to be involved in biological meaningful processes associated with ovarian cancer, such as CCL11, CCL16, CCL18, CCL23, CCL8, CCL5, APOB, BRCA1, SLC18A1, FGF22, GADD45B, GNA15, GNA11, and so on.

The functional consequences and prognostic value of dosage sensitivity in ovarian cancer

Copy number alteration (CNA) represents an important class of genetic variations that may contribute to tumorigenesis, tumor growth and metastatic spread. CNA can directly affect the expression of genes within the CNA regions; however, genes within the CNA regions exhibit heterogeneity in gene dosage sensitivity. In this study, a computational framework was built to identify 1170 dosage-sensitive genes (DSGs) and 1215 dosage-resistant genes (DRGs) that were related to ovarian serous cystadenocarcinoma (OV) through the association between CNA and gene expression. To analyze the different functions of the genes within the two groups, the functional annotation results indicated that DRGs were involved in cancer-related processes like immune response, cell death and apoptosis, while DSGs were enriched in essential processes like the cell cycle and the DNA metabolic process. Meanwhile, two three-dimensional regulatory networks for differentially expressed miRNAs, differentially expressed transcription factors (TFs) and DSGs or DRGs were constructed based on feed-forward loops. We identified key regulators (such as miR-16-5p, miR-98-5p, MYB and HOXA5) and cancer prognosis-related network motifs (such as miR-98-5p-HOXA5-TP53 and miR-16-5p-MYB-IGF1R) after the analysis of network topological features. Our results lead us to speculate that these genes and associated regulators may be potential mechanistic biomarkers for tumorigenesis and progression of cancer. Research on the network characteristics and the role of feed-forward loops in OV tumorigenesis and development could lead to feasible suggestions for the prevention and early diagnosis of OV, which will shed light on understanding the functional mechanism of CNA in cancer.

iGC-an integrated analysis package of gene expression and copy number alteration

BACKGROUND

With the advancement in high-throughput technologies, researchers can simultaneously investigate gene expression and copy number alteration (CNA) data from individual patients at a lower cost. Traditional analysis methods analyze each type of data individually and integrate their results using Venn diagrams. Challenges arise, however, when the results are irreproducible and inconsistent across multiple platforms. To address these issues, one possible approach is to concurrently analyze both gene expression profiling and CNAs in the same individual.

RESULTS

We have developed an open-source R/Bioconductor package (iGC). Multiple input formats are supported and users can define their own criteria for identifying differentially expressed genes driven by CNAs. The analysis of two real microarray datasets demonstrated that the CNA-driven genes identified by the iGC package showed significantly higher Pearson correlation coefficients with their gene expression levels and copy numbers than those genes located in a genomic region with CNA. Compared with the Venn diagram approach, the iGC package showed better performance.

CONCLUSION

The iGC package is effective and useful for identifying CNA-driven genes. By simultaneously considering both comparative genomic and transcriptomic data, it can provide better understanding of biological and medical questions. The iGC package's source code and manual are freely available at https://www.bioconductor.org/packages/release/bioc/html/iGC.html .

Deletions of multidrug resistance gene loci in breast cancer leads to the down-regulation of its expression and predict tumor response to neoadjuvant chemotherapy

Neoadjuvant chemotherapy (NAC) is intensively used for the treatment of primary breast cancer. In our previous studies, we reported that clinical tumor response to NAC is associated with the change of multidrug resistance (MDR) gene expression in tumors after chemotherapy. In this study we performed a combined analysis of MDR gene locus deletions in tumor DNA, MDR gene expression and clinical response to NAC in 73 BC patients. Copy number variations (CNVs) in biopsy specimens were tested using high-density microarray platform CytoScanTM HD Array (Affymetrix, USA). 75%–100% persons having deletions of MDR gene loci demonstrated the down-regulation of MDR gene expression. Expression of MDR genes was 2–8 times lower in patients with deletion than in patients having no deletion only in post-NAC tumors samples but not in tumor tissue before chemotherapy. All patients with deletions of ABCB1 ABCB 3 ABCC5 gene loci – 7q21.1, 6p21.32, 3q27 correspondingly, and most patients having deletions in ABCC1 (16p13.1), ABCC2 (10q24), ABCG1 (21q22.3), ABCG2 (4q22.1), responded favorably to NAC. The analysis of all CNVs, including both amplification and deletion showed that the frequency of 13q14.2 deletion was 85% among patients bearing tumor with the deletion at least in one MDR gene locus versus 9% in patients with no deletions. Differences in the frequency of 13q14.2 deletions between the two groups were statistically significant (p = 2.03 ×10⁻¹¹, Fisher test, Bonferroni-adjusted p = 1.73 × 10⁻⁸). In conclusion, our study for the first time demonstrates that deletion MDR gene loci can be used as predictive marker for tumor response to NAC.

FocalScan: Scanning for altered genes in cancer based on coordinated DNA and RNA change

Somatic genomic copy-number alterations can lead to transcriptional activation or inactivation of tumor driver or suppressor genes, contributing to the malignant properties of cancer cells. Selection for such events may manifest as recurrent amplifications or deletions of size-limited (focal) regions. While methods have been developed to identify such focal regions, finding the exact targeted genes remains a challenge. Algorithms are also available that integrate copy number and RNA expression data, to aid in identifying individual targeted genes, but specificity is lacking. Here, we describe FocalScan, a tool designed to simultaneously uncover patterns of focal copy number alteration and coordinated expression change, thus combining both principles. The method outputs a ranking of tentative cancer drivers or suppressors. FocalScan works with RNA-seq data, and unlike other tools it can scan the genome unaided by a gene annotation, enabling identification of novel putatively functional elements including lncRNAs. Application on a breast cancer data set suggests considerably better performance than other DNA/RNA integration tools.

A comparative analysis reveals the dosage sensitivity and regulatory patterns of lncRNA in prostate cancer

Although the key roles of long non-coding RNAs (lncRNAs) in multiple diseases are well documented, the relationship between the lncRNA copy number and expression is unknown. Here, we present a comprehensive study that demonstrates the impact of miRNA-TF co-regulatory motifs on the dosage effects of protein-coding genes (PCGs) and lncRNAs in prostate cancer. By integrating copy number profiles, expression profiles and regulatory relationships with miRNAs and transcription factors, we reveal that lncRNAs and PCGs correlate with distinct dosage sensitivity and regulatory pattern characteristics. We also show that lncRNAs from different genomic regions have different features. Using a custom-built framework, we identified 24 candidate prostate cancer-related lncRNAs based on the known properties of established prostate-related lncRNAs. Our method will enable the identification of cancer-related lncRNAs, which will provide new insights into cancer lncRNA biology.

Identifying candidate drivers of drug response in heterogeneous cancer by mining high throughput genomics data

BACKGROUND

With advances in technologies, huge amounts of multiple types of high-throughput genomics data are available. These data have tremendous potential to identify new and clinically valuable biomarkers to guide the diagnosis, assessment of prognosis, and treatment of complex diseases, such as cancer. Integrating, analyzing, and interpreting big and noisy genomics data to obtain biologically meaningful results, however, remains highly challenging. Mining genomics datasets by utilizing advanced computational methods can help to address these issues.

RESULTS

To facilitate the identification of a short list of biologically meaningful genes as candidate drivers of anti-cancer drug resistance from an enormous amount of heterogeneous data, we employed statistical machine-learning techniques and integrated genomics datasets. We developed a computational method that integrates gene expression, somatic mutation, and copy number aberration data of sensitive and resistant tumors. In this method, an integrative method based on module network analysis is applied to identify potential driver genes. This is followed by cross-validation and a comparison of the results of sensitive and resistance groups to obtain the final list of candidate biomarkers. We applied this method to the ovarian cancer data from the cancer genome atlas. The final result contains biologically relevant genes, such as COL11A1, which has been reported as a cis-platinum resistant biomarker for epithelial ovarian carcinoma in several recent studies.

CONCLUSIONS

The described method yields a short list of aberrant genes that also control the expression of their co-regulated genes. The results suggest that the unbiased data driven computational method can identify biologically relevant candidate biomarkers. It can be utilized in a wide range of applications that compare two conditions with highly heterogeneous datasets.

Transcriptome Analysis in Domesticated Species: Challenges and Strategies

Domesticated species occupy a special place in the human world due to their economic and cultural value. In the era of genomic research, domesticated species provide unique advantages for investigation of diseases and complex phenotypes. RNA sequencing, or RNA-seq, has recently emerged as a new approach for studying transcriptional activity of the whole genome, changing the focus from individual genes to gene networks. RNA-seq analysis in domesticated species may complement genome-wide association studies of complex traits with economic importance or direct relevance to biomedical research. However, RNA-seq studies are more challenging in domesticated species than in model organisms. These challenges are at least in part associated with the lack of quality genome assemblies for some domesticated species and the absence of genome assemblies for others. In this review, we discuss strategies for analyzing RNA-seq data, focusing particularly on questions and examples relevant to domesticated species.

Analysing multiple types of molecular profiles simultaneously: connecting the needles in the haystack

Background

It has been shown that a random-effects framework can be used to test the association between a gene’s expression level and the number of DNA copies of a set of genes. This gene-set modelling framework was later applied to find associations between mRNA expression and microRNA expression, by defining the gene sets using target prediction information.

Methods and results

Here, we extend the model introduced by Menezes et al. 2009 to consider the effect of not just copy number, but also of other molecular profiles such as methylation changes and loss-of-heterozigosity (LOH), on gene expression levels. We will consider again sets of measurements, to improve robustness of results and increase the power to find associations. Our approach can be used genome-wide to find associations and yields a test to help separate true associations from noise.

We apply our method to colon and to breast cancer samples, for which genome-wide copy number, methylation and gene expression profiles are available. Our findings include interesting gene expression-regulating mechanisms, which may involve only one of copy number or methylation, or both for the same samples. We even are able to find effects due to different molecular mechanisms in different samples.

Conclusions

Our method can equally well be applied to cases where other types of molecular (high-dimensional) data are collected, such as LOH, SNP genotype and microRNA expression data. Computationally efficient, it represents a flexible and powerful tool to study associations between high-dimensional datasets. The method is freely available via the SIM BioConductor package.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0926-8) contains supplementary material, which is available to authorized users.

Understanding the functional impact of copy number alterations in breast cancer using a network modeling approach

We apply a network approach to identify genes associated incisor intranswith copy-number alterations in breast cancer pathogenesis.

Downregulation of transient receptor potential cation channel, subfamily C, member 1 contributes to drug resistance and high histological grade in ovarian cancer

Transient receptor potential cation channel, subfamily C, member 1 (TRPC1) participates in many physiological functions but has also been implicated in cancer development. However, little is known about the role of TRPC1 in ovarian cancer (OC), including the drug resistance of these tumors. In the present study, a significant and consistent downregulation of TRPC1 in drug-resistant OC tissues/cells was determined using real-time quantitative polymerase chain reaction assays and the microarrays deposited in Oncomine and Gene Expression Omnibus (GEO) profiles. Protein/gene-protein/gene and protein-chemical interactions indicated that TRPC1 interacts with 14 proteins/genes and 6 chemicals, all of which are involved in the regulation of drug resistance in OC. Biological process annotation of TRPC1, OC, and drug resistance indicated a role for TRPC1 in drug-resistance-related functions in OC, mainly via the cell cycle, gene expression and cell growth and cell death. Analysis of mRNA-microRNA interactions showed that 8 out of 11 major pathways enriched from 38 predominant microRNAs targeting TRPC1 were involved in the regulation of drug resistance in OC, and 8 out of these top 10 microRNAs were implicated in the drug resistance in ovarian and other cancers. In a clinical analysis using data obtained from The Cancer Genome Atlas project (TCGA) cohort on 341 OC patients, TRPC1 expression was found to differ significantly between grade 2 and grade 3 tumors, with low-level expression correlating with higher tumor grade. This is the first report to show a potential association between the downregulation of TRPC1 and both drug resistance and high histological tumor grade in OC. Our results provide the basis for further investigations of the drug-resistance-related functions of TRPC1 in OC and other forms of cancer.

Investigating inter-chromosomal regulatory relationships through a comprehensive meta-analysis of matched copy number and transcriptomics data sets

Background

Gene regulatory relationships can be inferred using matched array comparative genomics and transcriptomics data sets from cancer samples. The way in which copy numbers of genes in cancer samples are often greatly disrupted works like a natural gene amplification/deletion experiment. There are now a large number of such data sets publicly available making a meta-analysis of the data possible.

Results

We infer inter-chromosomal acting gene regulatory relationships from a meta-analysis of 31 publicly available matched array comparative genomics and transcriptomics data sets in humans. We obtained statistically significant predictions of target genes for 1430 potential regulatory genes. The regulatory relationships being inferred are either direct relationships, of a transcription factor on its target, or indirect ones, through pathways containing intermediate steps. We analyse the predictions in terms of cocitations, both publications which cite a regulator with any of its inferred targets and cocitations of any genes in a target list.

Conclusions

The most striking observation from the results is the greater number of inter-chromosomal regulatory relationships involving repression compared to those involving activation. The complete results of the meta-analysis are presented in the database METAMATCHED. We anticipate that the predictions contained in the database will be useful in informing experiments and in helping to construct networks of regulatory relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2100-5) contains supplementary material, which is available to authorized users.

Systematic analysis of somatic mutations impacting gene expression in 12 tumour types

We present a novel hierarchical Bayes statistical model, xseq, to systematically quantify the impact of somatic mutations on expression profiles. We establish the theoretical framework and robust inference characteristics of the method using computational benchmarking. We then use xseq to analyse thousands of tumour data sets available through The Cancer Genome Atlas, to systematically quantify somatic mutations impacting expression profiles. We identify 30 novel cis-effect tumour suppressor gene candidates, enriched in loss-of-function mutations and biallelic inactivation. Analysis of trans-effects of mutations and copy number alterations with xseq identifies mutations in 150 genes impacting expression networks, with 89 novel predictions. We reveal two important novel characteristics of mutation impact on expression: (1) patients harbouring known driver mutations exhibit different downstream gene expression consequences; (2) expression patterns for some mutations are stable across tumour types. These results have critical implications for identification and interpretation of mutations with consequent impact on transcription in cancer.

Assessing functional impact of mutations in cancer on gene expression can improve our understanding of cancer biology and may identify potential therapeutic targets. Here, Ding et al. describe a novel statistical model named xseq for a systematic survey of how mutations impact transcriptome landscapes across 12 different tumour types.

Integrating genetics and epigenetics in breast cancer: biological insights, experimental, computational methods and therapeutic potential

BACKGROUND

Development of human cancer can proceed through the accumulation of different genetic changes affecting the structure and function of the genome. Combined analyses of molecular data at multiple levels, such as DNA copy-number alteration, mRNA and miRNA expression, can clarify biological functions and pathways deregulated in cancer. The integrative methods that are used to investigate these data involve different fields, including biology, bioinformatics, and statistics.

RESULTS

These methodologies are presented in this review, and their implementation in breast cancer is discussed with a focus on integration strategies. We report current applications, recent studies and interesting results leading to the identification of candidate biomarkers for diagnosis, prognosis, and therapy in breast cancer by using both individual and combined analyses.

CONCLUSION

This review presents a state of art of the role of different technologies in breast cancer based on the integration of genetics and epigenetics, and shares some issues related to the new opportunities and challenges offered by the application of such integrative approaches.

An integrated approach to reveal miRNAs' impacts on the functional consequence of copy number alterations in cancer

Copy number alteration (CNA) is known to induce gene expression changes mainly through dosage effect, and therefore affect the initiation and progression of tumor. However, tumor samples exhibit heterogeneity in gene dosage sensitivity due to the complicated mechanisms of transcriptional regulation. Currently, no high-throughput method has been available for identifying the regulatory factors affecting the functional consequences of CNA, and determining their effects on cancer. In view of the important regulatory role of miRNA, we investigated the influence of miRNAs on the dosage sensitivities of genes within the CNA regions. By integrating copy number, mRNA expression, miRNA expression profiles of three kinds of cancer, we observed a tendency for high dosage-sensitivity genes to be more targeted by miRNAs in cancer, and identified the miRNAs regulating the dosage sensitivity of amplified/deleted target genes. The results show that miRNAs can modulate oncogenic biological functions by regulating the genes within the CNA regions, and thus play a role as a trigger or balancer in cancer, affecting cancer processes, even survival. This work provided a framework for analyzing the regulation of dosage effect, which will shed a light on understanding the oncogenic and tumor suppressive mechanisms of CNA. Besides, new cancer-related miRNAs were identified.

Discovering key regulatory mechanisms from single-factor and multi-factor regulations in glioblastoma utilizing multi-dimensional data

Glioblastoma (GBM) is the most common malignant brain cancer in adults. Investigating the regulatory mechanisms underlying GBM is effective for the in-depth study of GBM. The Cancer Genome Atlas (TCGA) project is producing large-scale data and makes the comprehensive study of the diverse regulatory mechanisms underlying GBM possible. Although there have been research studies on GBM with large-scale data, distinguishing different regulatory mechanisms and identifying the key regulation types remain challenging. In this study, we integrated multi-dimensional data of differentially expressed genes in GBM: copy number variation (CNV), gene expression, miRNA expression and methylation, by performing partial correlation analysis with the Lasso technique. Our results showed that there were single-factor and multi-factor regulatory mechanisms in GBM. In further analysis of the regulation subtypes, we discovered that single-factor and multi-factor regulations are potentially distinct in functionality. Moreover, multi-factor regulations especially the key regulation subtypes may be more relevant to GBM and affect many GBM-related genes such as ERBB2 and MAPK1. This study not only verifies the utility of multi-dimensional data integration into GBM research but also distinguishes the key multi-factor regulatory subtypes that may drive pathogenesis of GBM from various regulatory mechanisms.

ICan: an integrated co-alteration network to identify ovarian cancer-related genes

BACKGROUND

Over the last decade, an increasing number of integrative studies on cancer-related genes have been published. Integrative analyses aim to overcome the limitation of a single data type, and provide a more complete view of carcinogenesis. The vast majority of these studies used sample-matched data of gene expression and copy number to investigate the impact of copy number alteration on gene expression, and to predict and prioritize candidate oncogenes and tumor suppressor genes. However, correlations between genes were neglected in these studies. Our work aimed to evaluate the co-alteration of copy number, methylation and expression, allowing us to identify cancer-related genes and essential functional modules in cancer.

RESULTS

We built the Integrated Co-alteration network (ICan) based on multi-omics data, and analyzed the network to uncover cancer-related genes. After comparison with random networks, we identified 155 ovarian cancer-related genes, including well-known (TP53, BRCA1, RB1 and PTEN) and also novel cancer-related genes, such as PDPN and EphA2. We compared the results with a conventional method: CNAmet, and obtained a significantly better area under the curve value (ICan: 0.8179, CNAmet: 0.5183).

CONCLUSION

In this paper, we describe a framework to find cancer-related genes based on an Integrated Co-alteration network. Our results proved that ICan could precisely identify candidate cancer genes and provide increased mechanistic understanding of carcinogenesis. This work suggested a new research direction for biological network analyses involving multi-omics data.

Evaluation of copy number variation and gene expression in neurofibromatosis type-1-associated malignant peripheral nerve sheath tumours

BACKGROUND

Neurofibromatosis type-1 (NF1) is a complex neurogenetic disorder characterised by the development of benign and malignant tumours of the peripheral nerve sheath (MPNSTs). Whilst biallelic NF1 gene inactivation contributes to benign tumour formation, additional cellular changes in gene structure and/or expression are required to induce malignant transformation. Although few molecular profiling studies have been performed on the process of progression of pre-existing plexiform neurofibromas to MPNSTs, the integrated analysis of copy number alterations (CNAs) and gene expression is likely to be key to understanding the molecular mechanisms underlying NF1-MPNST tumorigenesis. In a pilot study, we employed this approach to identify genes differentially expressed between benign and malignant NF1 tumours.

RESULTS

SPP1 (osteopontin) was the most differentially expressed gene (85-fold increase in expression), compared to benign plexiform neurofibromas. Short hairpin RNA (shRNA) knockdown of SPP1 in NF1-MPNST cells reduced tumour spheroid size, wound healing and invasion in four different MPNST cell lines. Seventy-six genes were found to exhibit concordance between CNA and gene expression level.

CONCLUSIONS

Pathway analysis of these genes suggested that glutathione metabolism and Wnt signalling may be specifically involved in NF1-MPNST development. SPP1 is associated with malignant transformation in NF1-associated MPNSTs and could prove to be an important target for therapeutic intervention.

A meta-analysis of multiple matched copy number and transcriptomics data sets for inferring gene regulatory relationships

Inferring gene regulatory relationships from observational data is challenging. Manipulation and intervention is often required to unravel causal relationships unambiguously. However, gene copy number changes, as they frequently occur in cancer cells, might be considered natural manipulation experiments on gene expression. An increasing number of data sets on matched array comparative genomic hybridisation and transcriptomics experiments from a variety of cancer pathologies are becoming publicly available. Here we explore the potential of a meta-analysis of thirty such data sets. The aim of our analysis was to assess the potential of in silico inference of trans-acting gene regulatory relationships from this type of data. We found sufficient correlation signal in the data to infer gene regulatory relationships, with interesting similarities between data sets. A number of genes had highly correlated copy number and expression changes in many of the data sets and we present predicted potential trans-acted regulatory relationships for each of these genes. The study also investigates to what extent heterogeneity between cell types and between pathologies determines the number of statistically significant predictions available from a meta-analysis of experiments.

Identification of metastasis-associated genes in colorectal cancer through an integrated genomic and transcriptomic analysis

OBJECTIVE

Identification of colorectal cancer (CRC) metastasis genes is one of the most important issues in CRC research. For the purpose of mining CRC metastasis-associated genes, an integrated analysis of microarray data was presented, by combined with evidence acquired from comparative genomic hybridization (CGH) data.

METHODS

Gene expression profile data of CRC samples were obtained at Gene Expression Omnibus (GEO) website. The 15 important chromosomal aberration sites detected by using CGH technology were used for integrated genomic and transcriptomic analysis. Significant Analysis of Microarray (SAM) was used to detect significantly differentially expressed genes across the whole genome. The overlapping genes were selected in their corresponding chromosomal aberration regions, and analyzed by using the Database for Annotation, Visualization and Integrated Discovery (DAVID). Finally, SVM-T-RFE gene selection algorithm was applied to identify metastasis-associated genes in CRC.

RESULTS

A minimum gene set was obtained with the minimum number [14] of genes, and the highest classification accuracy (100%) in both PRI and META datasets. A fraction of selected genes are associated with CRC or its metastasis.

CONCLUSIONS

Our results demonstrated that integration analysis is an effective strategy for mining cancer-associated genes.

The shaping and functional consequences of the dosage effect landscape in multiple myeloma

Background

Multiple myeloma (MM) is a malignant proliferation of plasma B cells. Based on recurrent aneuploidy such as copy number alterations (CNAs), myeloma is divided into two subtypes with different CNA patterns and patient survival outcomes. How aneuploidy events arise, and whether they contribute to cancer cell evolution are actively studied. The large amount of transcriptomic changes resultant of CNAs (dosage effect) pose big challenges for identifying functional consequences of CNAs in myeloma in terms of specific driver genes and pathways. In this study, we hypothesize that gene-wise dosage effect varies as a result from complex regulatory networks that translate the impact of CNAs to gene expression, and studying this variation can provide insights into functional effects of CNAs.

Results

We propose gene-wise dosage effect score and genome-wide karyotype plot as tools to measure and visualize concordant copy number and expression changes across cancer samples. We find that dosage effect in myeloma is widespread yet variable, and it is correlated with gene expression level and CNA frequencies in different chromosomes. Our analysis suggests that despite the enrichment of differentially expressed genes between hyperdiploid MM and non-hyperdiploid MM in the trisomy chromosomes, the chromosomal proportion of dosage sensitive genes is higher in the non-trisomy chromosomes. Dosage-sensitive genes are enriched by genes with protein translation and localization functions, and dosage resistant genes are enriched by apoptosis genes. These results point to future studies on differential dosage sensitivity and resistance of pro- and anti-proliferation pathways and their variation across patients as therapeutic targets and prognosis markers.

Conclusions

Our findings support the hypothesis that recurrent CNAs in myeloma are selected by their functional consequences. The novel dosage effect score defined in this work will facilitate integration of copy number and expression data for identifying driver genes in cancer genomics studies. The accompanying R code is available at http://www.canevolve.org/dosageEffect/.

canEvolve: a web portal for integrative oncogenomics

Background & Objective

Genome-wide profiles of tumors obtained using functional genomics platforms are being deposited to the public repositories at an astronomical scale, as a result of focused efforts by individual laboratories and large projects such as the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium. Consequently, there is an urgent need for reliable tools that integrate and interpret these data in light of current knowledge and disseminate results to biomedical researchers in a user-friendly manner. We have built the canEvolve web portal to meet this need.

Results

canEvolve query functionalities are designed to fulfill most frequent analysis needs of cancer researchers with a view to generate novel hypotheses. canEvolve stores gene, microRNA (miRNA) and protein expression profiles, copy number alterations for multiple cancer types, and protein-protein interaction information. canEvolve allows querying of results of primary analysis, integrative analysis and network analysis of oncogenomics data. The querying for primary analysis includes differential gene and miRNA expression as well as changes in gene copy number measured with SNP microarrays. canEvolve provides results of integrative analysis of gene expression profiles with copy number alterations and with miRNA profiles as well as generalized integrative analysis using gene set enrichment analysis. The network analysis capability includes storage and visualization of gene co-expression, inferred gene regulatory networks and protein-protein interaction information. Finally, canEvolve provides correlations between gene expression and clinical outcomes in terms of univariate survival analysis.

Conclusion

At present canEvolve provides different types of information extracted from 90 cancer genomics studies comprising of more than 10,000 patients. The presence of multiple data types, novel integrative analysis for identifying regulators of oncogenesis, network analysis and ability to query gene lists/pathways are distinctive features of canEvolve. canEvolve will facilitate integrative and meta-analysis of oncogenomics datasets.

Availability

The canEvolve web portal is available at http://www.canevolve.org/.

Genome Fusion Detection: a novel method to detect fusion genes from SNP-array data

MOTIVATION

Fusion genes result from genomic rearrangements, such as deletions, amplifications and translocations. Such rearrangements can also frequently be observed in cancer and have been postulated as driving event in cancer development. to detect them, one needs to analyze the transition region of two segments with different copy number, the location where fusions are known to occur. Finding fusion genes is essential to understanding cancer development and may lead to new therapeutic approaches.

RESULTS

Here we present a novel method, the Genomic Fusion Detection algorithm, to predict fusion genes on a genomic level based on SNP-array data. This algorithm detects genes at the transition region of segments with copy number variation. With the application of defined constraints, certain properties of the detected genes are evaluated to predict whether they may be fused. We evaluated our prediction by calculating the observed frequency of known fusions in both primary cancers and cell lines. We tested a set of cell lines positive for the BCR-ABL1 fusion and prostate cancers positive for the TMPRSS2-ERG fusion. We could detect the fusions in all positive cell lines, but not in the negative controls.

integIRTy: a method to identify genes altered in cancer by accounting for multiple mechanisms of regulation using item response theory

MOTIVATION

Identifying genes altered in cancer plays a crucial role in both understanding the mechanism of carcinogenesis and developing novel therapeutics. It is known that there are various mechanisms of regulation that can lead to gene dysfunction, including copy number change, methylation, abnormal expression, mutation and so on. Nowadays, all these types of alterations can be simultaneously interrogated by different types of assays. Although many methods have been proposed to identify altered genes from a single assay, there is no method that can deal with multiple assays accounting for different alteration types systematically.

RESULTS

In this article, we propose a novel method, integration using item response theory (integIRTy), to identify altered genes by using item response theory that allows integrated analysis of multiple high-throughput assays. When applied to a single assay, the proposed method is more robust and reliable than conventional methods such as Student's t-test or the Wilcoxon rank-sum test. When used to integrate multiple assays, integIRTy can identify novel-altered genes that cannot be found by looking at individual assay separately. We applied integIRTy to three public cancer datasets (ovarian carcinoma, breast cancer, glioblastoma) for cross-assay type integration which all show encouraging results.

AVAILABILITY AND IMPLEMENTATION

The R package integIRTy is available at the web site http://bioinformatics.mdanderson.org/main/OOMPA:Overview.

CONTACT

kcoombes@mdanderson.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recurrent transcriptional clusters in the genome of mouse pluripotent stem cells

A number of studies have shown that transcriptome analysis in terms of chromosomal location can reveal regions of non-random transcriptional activity within the genome. Genomic clusters of differentially expressed genes can identify genomic patterns of structural organization, underlying copy number variations or long-range epigenetic regulation such as X-chromosome inactivation. Here we apply an integrative bioinformatics analysis to a collection of 315 freely available mouse pluripotent stem cell samples to discover transcriptional clusters in the genome. We show that over half of the analysed samples (56.83%) carry whole or partial-chromosome spanning clusters which recur in genomic regions previously implicated in chromosomal imbalances. Strikingly, we found that the presence of such large-clusters is linked to the differential expression of a limited number of genes, common to all samples carrying clusters irrespectively of the chromosome where the cluster is found. We have used these genes to train and test classification models that can predict samples that carry large-scale clusters on any chromosome with over 90% accuracy. Our findings suggest that there is a common downstream activation in these cells that affects a limited number of nodes. We propose that this effect is linked to selective advantage and identify potential driver genes.

Cancer gene prioritization by integrative analysis of mRNA expression and DNA copy number data: a comparative review

A variety of genome-wide profiling techniques are available to investigate complementary aspects of genome structure and function. Integrative analysis of heterogeneous data sources can reveal higher level interactions that cannot be detected based on individual observations. A standard integration task in cancer studies is to identify altered genomic regions that induce changes in the expression of the associated genes based on joint analysis of genome-wide gene expression and copy number profiling measurements. In this review, we highlight common approaches to genomic data integration and provide a transparent benchmarking procedure to quantitatively compare method performances in cancer gene prioritization. Algorithms, data sets and benchmarking results are available at http://intcomp.r-forge.r-project.org.