Detecting cancer gene networks characterized by recurrent genomic alterations in a population

High resolution, system-wide characterizations have demonstrated the capacity to identify genomic regions that undergo genomic aberrations. Such research efforts often aim at associating these regions with disease etiology and outcome. Identifying the corresponding biologic processes that are responsible for disease and its outcome remains challenging. Using novel analytic methods that utilize the structure of biologic networks, we are able to identify the specific networks that are highly significantly, nonrandomly altered by regions of copy number amplification observed in a systems-wide analysis. We demonstrate this method in breast cancer, where the state of a subset of the pathways identified through these regions is shown to be highly associated with disease survival and recurrence.     

Elucidating cis- and trans-regulatory variation using genetical genomics

Genetic variation, through its effects on gene regulatory networks, has an important role in determining biologic phenotypes and disease susceptibility. Recent studies using segregating populations have characterized variation in gene expression resulting from genetic differences and have used well-established assays to identify many cis-acting polymorphisms. Although trans effects have been shown to be globally important, their individual effects are small, such that their genome-wide identification remains a challenge. Two recent studies using approaches capitalizing on intercrossed mouse strains and novel computational methods have yielded examples of trans-regulatory networks.     

Principles and strategies for developing network models in cancer

The flood of genome-wide data generated by high-throughput technologies currently provides biologists with an unprecedented opportunity: to manipulate, query, and reconstruct functional molecular networks of cells. Here, we outline three underlying principles and six strategies to infer network models from genomic data. Then, using cancer as an example, we describe experimental and computational approaches to infer "differential" networks that can identify genes and processes driving disease phenotypes. In conclusion, we discuss how a network-level understanding of cancer can be used to predict drug response and guide therapeutics.     

Design and analysis issues in genome-wide somatic mutation studies of cancer

The availability of the human genome sequence and progress in sequencing and bioinformatic technologies have enabled genome-wide investigation of somatic mutations in human cancers. This article briefly reviews challenges arising in the statistical analysis of mutational data of this kind. A first challenge is that of designing studies that efficiently allocate sequencing resources. We show that this can be addressed by two-stage designs and demonstrate via simulations that even relatively small studies can produce lists of candidate cancer genes that are highly informative for future research efforts. A second challenge is to distinguish mutated genes that are selected for by cancer (drivers) from mutated genes that have no role in the development of cancer and simply happened to mutate (passengers). We suggest that this question is best approached as a classification problem and discuss some of the difficulties of more traditional testing-based approaches. A third challenge is to identify biologic processes affected by the driver genes. This can be pursued by gene set analyses. These can reliably identify functional groups and pathways that are enriched for mutated genes even when the individual genes involved in those pathways or sets are not mutated at sufficient frequencies to provide conclusive evidence as drivers.     

Pathway-based genetic analysis of preterm birth

Preterm birth in the United States is now 12%. Multiple genes, gene networks, and variants have been associated with this disease. Using a custom database for preterm birth (dbPTB) with a refined set of genes extensively curated from literature and biological databases, we analyzed GWAS of preterm birth for complete genotype data on nearly 2000 preterm and term mothers. We used both the curated genes and a genome-wide approach to carry out a pathway-based analysis. There were 19 significant pathways, which withstood FDR correction for multiple testing that were identified using both the curated genes and the genome-wide approach. The analysis based on the curated genes was more significant than genome-wide in 15 out of 19 pathways. This approach demonstrates the use of a validated set of genes, in the analysis of otherwise unsuccessful GWAS data, to identify gene–gene interactions in a way that enhances statistical power and discovery.     

Genotype-phenotype mapping in a post-GWAS world

Understanding how metabolic reactions, cell signaling, and developmental pathways translate the genome of an organism into its phenotype is a grand challenge in biology. Genome-wide association studies (GWAS) statistically connect genotypes to phenotypes, without any recourse to known molecular interactions, whereas a molecular biology approach directly ties gene function to phenotype through gene regulatory networks (GRNs). Using natural variation in allele-specific expression, GWAS and GRN approaches can be merged into a single framework via structural equation modeling (SEM). This approach leverages the myriad of polymorphisms in natural populations to elucidate and quantitate the molecular pathways that underlie phenotypic variation. The SEM framework can be used to quantitate a GRN, evaluate its consistency across environments or sexes, identify the differences in GRNs between species, and annotate GRNs de novo in non-model organisms.     

Circulating tumor cells: highlight on practical implications

Circulating tumor cells (CTCs) are cells of presumed epithelial origin, whose prognostic and predictive value in metastatic cancer patients has recently been demonstrated. To date, the count of CTCs through the CellSearch? system represents a valid approach for monitoring disease status in patients with metastatic colorectal, breast, and prostate cancer; in these cancer types, a rise in the CTC count at any time during treatment predicts a poor outcome. Nevertheless, the clinical utility of monitoring CTC counts remains controversial, and what to do when CTC counts rise during therapy still remains an unanswered question. In this report, we suggest how to integrate CTC counts with their molecular characterization to better translate biologic information obtained on CTCs into daily clinical practice.     

Broad network-based predictability of <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>gene loss-of-function phenotypes

We demonstrate that loss-of-function yeast phenotypes are predictable by guilt-by-association in functional gene networks. Testing 1,102 loss-of-function phenotypes from genome-wide assays of yeast reveals predictability of diverse phenotypes, spanning cellular morphology, growth, metabolism, and quantitative cell shape features. We apply the method to extend a genome-wide screen by predicting, then verifying, genes whose disruption elongates yeast cells, and to predict human disease genes. To facilitate network-guided screens, a web server is available .     

Comparative proteomic analysis of cysteine oxidation in colorectal cancer patients

Oxidative stress promotes damage to cellular proteins, lipids, membranes and DNA, and plays a key role in the development of cancer. Reactive oxygen species disrupt redox homeostasis and promote tumor formation by initiating aberrant activation of signaling pathways that lead to tumorigenesis. We used shotgun proteomics to identify proteins containing oxidation-sensitive cysteines in tissue specimens from colorectal cancer patients. We then compared the patterns of cysteine oxidation in the membrane fractions between the tumor and non-tumor tissues. Using nano-UPLC-MS^E proteomics, we identified 31 proteins containing 37 oxidation-sensitive cysteines. These proteins were observed with IAM-binding cysteines in non-tumoral region more than tumoral region of CRC patients. Then using the Ingenuity pathway program, we evaluated the cellular canonical networks connecting those proteins. Within the networks, proteins with multiple connections were related with organ morphology, cellular metabolism, and various disorders. We have thus identified networks of proteins whose redox status is altered by oxidative stress, perhaps leading to changes in cellular functionality that promotes tumorigenesis.     

Differences in the Tumor Microenvironment between African-American and European-American Breast Cancer Patients

Background

African-American breast cancer patients experience higher mortality rates than European-American patients despite having a lower incidence of the disease. We tested the hypothesis that intrinsic differences in the tumor biology may contribute to this cancer health disparity.

Methods and Results

Using laser capture microdissection, we examined genome-wide mRNA expression specific to tumor epithelium and tumor stroma in 18 African-American and 17 European-American patients. Numerous genes were differentially expressed between these two patient groups and a two-gene signature in the tumor epithelium distinguished between them. To identify the biological processes in tumors that are different by race/ethnicity, Gene Ontology and disease association analyses were performed. Several biological processes were identified which may contribute to enhanced disease aggressiveness in African-American patients, including angiogenesis and chemotaxis. African-American tumors also contained a prominent interferon signature. The role of angiogenesis in the tumor biology of African-Americans was further investigated by examining the extent of vascularization and macrophage infiltration in an expanded set of 248 breast tumors. Immunohistochemistry revealed that microvessel density and macrophage infiltration is higher in tumors of African-Americans than in tumors of European-Americans. Lastly, using an in silico approach, we explored the potential of tailored treatment options for African-American patients based on their gene expression profile. This exploratory approach generated lists of therapeutics that may have specific antagonistic activity against tumors of African-American patients, e.g., sirolimus, resveratrol, and chlorpromazine in estrogen receptor-negative tumors.

Conclusions

The gene expression profiles of breast tumors indicate that differences in tumor biology may exist between African-American and European-American patients beyond the knowledge of current markers. Notably, pathways related to tumor angiogenesis and chemotaxis could be functionally different in these two patient groups.     

Exploiting Protein-Protein Interaction Networks for Genome-Wide Disease-Gene Prioritization

Complex genetic disorders often involve products of multiple genes acting cooperatively. Hence, the pathophenotype is the outcome of the perturbations in the underlying pathways, where gene products cooperate through various mechanisms such as protein-protein interactions. Pinpointing the decisive elements of such disease pathways is still challenging. Over the last years, computational approaches exploiting interaction network topology have been successfully applied to prioritize individual genes involved in diseases. Although linkage intervals provide a list of disease-gene candidates, recent genome-wide studies demonstrate that genes not associated with any known linkage interval may also contribute to the disease phenotype. Network based prioritization methods help highlighting such associations. Still, there is a need for robust methods that capture the interplay among disease-associated genes mediated by the topology of the network. Here, we propose a genome-wide network-based prioritization framework named GUILD. This framework implements four network-based disease-gene prioritization algorithms. We analyze the performance of these algorithms in dozens of disease phenotypes. The algorithms in GUILD are compared to state-of-the-art network topology based algorithms for prioritization of genes. As a proof of principle, we investigate top-ranking genes in Alzheimer''s disease (AD), diabetes and AIDS using disease-gene associations from various sources. We show that GUILD is able to significantly highlight disease-gene associations that are not used a priori. Our findings suggest that GUILD helps to identify genes implicated in the pathology of human disorders independent of the loci associated with the disorders.     

A Multi-Method Approach for Proteomic Network Inference in 11 Human Cancers

Protein expression and post-translational modification levels are tightly regulated in neoplastic cells to maintain cellular processes known as ‘cancer hallmarks’. The first Pan-Cancer initiative of The Cancer Genome Atlas (TCGA) Research Network has aggregated protein expression profiles for 3,467 patient samples from 11 tumor types using the antibody based reverse phase protein array (RPPA) technology. The resultant proteomic data can be utilized to computationally infer protein-protein interaction (PPI) networks and to study the commonalities and differences across tumor types. In this study, we compare the performance of 13 established network inference methods in their capacity to retrieve the curated Pathway Commons interactions from RPPA data. We observe that no single method has the best performance in all tumor types, but a group of six methods, including diverse techniques such as correlation, mutual information, and regression, consistently rank highly among the tested methods. We utilize the high performing methods to obtain a consensus network; and identify four robust and densely connected modules that reveal biological processes as well as suggest antibody–related technical biases. Mapping the consensus network interactions to Reactome gene lists confirms the pan-cancer importance of signal transduction pathways, innate and adaptive immune signaling, cell cycle, metabolism, and DNA repair; and also suggests several biological processes that may be specific to a subset of tumor types. Our results illustrate the utility of the RPPA platform as a tool to study proteomic networks in cancer.     

MiRNA Profiles in Lymphoblastoid Cell Lines of Finnish Prostate Cancer Families

BackgroundHeritable factors are evidently involved in prostate cancer (PrCa) carcinogenesis, but currently, genetic markers are not routinely used in screening or diagnostics of the disease. More precise information is needed for making treatment decisions to distinguish aggressive cases from indolent disease, for which heritable factors could be a useful tool. The genetic makeup of PrCa has only recently begun to be unravelled through large-scale genome-wide association studies (GWAS). The thus far identified Single Nucleotide Polymorphisms (SNPs) explain, however, only a fraction of familial clustering. Moreover, the known risk SNPs are not associated with the clinical outcome of the disease, such as aggressive or metastasised disease, and therefore cannot be used to predict the prognosis. Annotating the SNPs with deep clinical data together with miRNA expression profiles can improve the understanding of the underlying mechanisms of different phenotypes of prostate cancer.ResultsIn this study microRNA (miRNA) profiles were studied as potential biomarkers to predict the disease outcome. The study subjects were from Finnish high risk prostate cancer families. To identify potential biomarkers we combined a novel non-parametrical test with an importance measure provided from a Random Forest classifier. This combination delivered a set of nine miRNAs that was able to separate cases from controls. The detected miRNA expression profiles could predict the development of the disease years before the actual PrCa diagnosis or detect the existence of other cancers in the studied individuals. Furthermore, using an expression Quantitative Trait Loci (eQTL) analysis, regulatory SNPs for miRNA miR-483-3p that were also directly associated with PrCa were found.ConclusionBased on our findings, we suggest that blood-based miRNA expression profiling can be used in the diagnosis and maybe even prognosis of the disease. In the future, miRNA profiling could possibly be used in targeted screening, together with Prostate Specific Antigene (PSA) testing, to identify men with an elevated PrCa risk.     

Incorporating regulatory interactions into gene-set analyses for GWAS data: A controlled analysis with the MAGMA tool

To date, genome-wide association studies have identified thousands of statistically-significant associations between genetic variants, and phenotypes related to a myriad of traits and diseases. A key goal for human-genetics research is to translate these associations into functional mechanisms. Popular gene-set analysis tools, like MAGMA, map variants to genes they might affect, and then integrate genome-wide association study data (that is, variant-level associations for a phenotype) to score genes for association with a phenotype. Gene scores are subsequently used in competitive gene-set analyses to identify biological processes that are enriched for phenotype association. By default, variants are mapped to genes in their proximity. However, many variants that affect phenotypes are thought to act at regulatory elements, which can be hundreds of kilobases away from their target genes. Thus, we explored the idea of augmenting a proximity-based mapping scheme with publicly-available datasets of regulatory interactions. We used MAGMA to analyze genome-wide association study data for ten different phenotypes, and evaluated the effects of augmentation by comparing numbers, and identities, of genes and gene sets detected as statistically significant between mappings. We detected several pitfalls and confounders of such “augmented analyses”, and introduced ways to control for them. Using these controls, we demonstrated that augmentation with datasets of regulatory interactions only occasionally strengthened the enrichment for phenotype association amongst (biologically-relevant) gene sets for different phenotypes. Still, in such cases, genes and regulatory elements responsible for the improvement could be pinpointed. For instance, using brain regulatory-interactions for augmentation, we were able to implicate two acetylcholine receptor subunits involved in post-synaptic chemical transmission, namely CHRNB2 and CHRNE, in schizophrenia. Collectively, our study presents a critical approach for integrating regulatory interactions into gene-set analyses for genome-wide association study data, by introducing various controls to distinguish genuine results from spurious discoveries.     

Mapping of signaling networks through synthetic genetic interaction analysis by RNAi

The analysis of synthetic genetic interaction networks can reveal how biological systems achieve a high level of complexity with a limited repertoire of components. Studies in yeast and bacteria have taken advantage of collections of deletion strains to construct matrices of quantitative interaction profiles and infer gene function. Yet comparable approaches in higher organisms have been difficult to implement in a robust manner. Here we report a method to identify genetic interactions in tissue culture cells through RNAi. By performing more than 70,000 pairwise perturbations of signaling factors, we identified >600 interactions affecting different quantitative phenotypes of Drosophila melanogaster cells. Computational analysis of this interaction matrix allowed us to reconstruct signaling pathways and identify a conserved regulator of Ras-MAPK signaling. Large-scale genetic interaction mapping by RNAi is a versatile, scalable approach for revealing gene function and the connectivity of cellular networks.     

Target Inhibition Networks: Predicting Selective Combinations of Druggable Targets to Block Cancer Survival Pathways

A recent trend in drug development is to identify drug combinations or multi-target agents that effectively modify multiple nodes of disease-associated networks. Such polypharmacological effects may reduce the risk of emerging drug resistance by means of attacking the disease networks through synergistic and synthetic lethal interactions. However, due to the exponentially increasing number of potential drug and target combinations, systematic approaches are needed for prioritizing the most potent multi-target alternatives on a global network level. We took a functional systems pharmacology approach toward the identification of selective target combinations for specific cancer cells by combining large-scale screening data on drug treatment efficacies and drug-target binding affinities. Our model-based prediction approach, named TIMMA, takes advantage of the polypharmacological effects of drugs and infers combinatorial drug efficacies through system-level target inhibition networks. Case studies in MCF-7 and MDA-MB-231 breast cancer and BxPC-3 pancreatic cancer cells demonstrated how the target inhibition modeling allows systematic exploration of functional interactions between drugs and their targets to maximally inhibit multiple survival pathways in a given cancer type. The TIMMA prediction results were experimentally validated by means of systematic siRNA-mediated silencing of the selected targets and their pairwise combinations, showing increased ability to identify not only such druggable kinase targets that are essential for cancer survival either individually or in combination, but also synergistic interactions indicative of non-additive drug efficacies. These system-level analyses were enabled by a novel model construction method utilizing maximization and minimization rules, as well as a model selection algorithm based on sequential forward floating search. Compared with an existing computational solution, TIMMA showed both enhanced prediction accuracies in cross validation as well as significant reduction in computation times. Such cost-effective computational-experimental design strategies have the potential to greatly speed-up the drug testing efforts by prioritizing those interventions and interactions warranting further study in individual cancer cases.     

A Predictive Phosphorylation Signature of Lung Cancer

Background

Aberrant activation of signaling pathways drives many of the fundamental biological processes that accompany tumor initiation and progression. Inappropriate phosphorylation of intermediates in these signaling pathways are a frequently observed molecular lesion that accompanies the undesirable activation or repression of pro- and anti-oncogenic pathways. Therefore, methods which directly query signaling pathway activation via phosphorylation assays in individual cancer biopsies are expected to provide important insights into the molecular “logic” that distinguishes cancer and normal tissue on one hand, and enables personalized intervention strategies on the other.

Results

We first document the largest available set of tyrosine phosphorylation sites that are, individually, differentially phosphorylated in lung cancer, thus providing an immediate set of drug targets. Next, we develop a novel computational methodology to identify pathways whose phosphorylation activity is strongly correlated with the lung cancer phenotype. Finally, we demonstrate the feasibility of classifying lung cancers based on multi-variate phosphorylation signatures.

Conclusions

Highly predictive and biologically transparent phosphorylation signatures of lung cancer provide evidence for the existence of a robust set of phosphorylation mechanisms (captured by the signatures) present in the majority of lung cancers, and that reliably distinguish each lung cancer from normal. This approach should improve our understanding of cancer and help guide its treatment, since the phosphorylation signatures highlight proteins and pathways whose phosphorylation should be inhibited in order to prevent unregulated proliferation.     

Extending pathways and processes using molecular interaction networks to analyse cancer genome data

Background  

Cellular processes and pathways, whose deregulation may contribute to the development of cancers, are often represented as cascades of proteins transmitting a signal from the cell surface to the nucleus. However, recent functional genomic experiments have identified thousands of interactions for the signalling canonical proteins, challenging the traditional view of pathways as independent functional entities. Combining information from pathway databases and interaction networks obtained from functional genomic experiments is therefore a promising strategy to obtain more robust pathway and process representations, facilitating the study of cancer-related pathways.     

Simultaneous analysis of all SNPs in genome-wide and re-sequencing association studies

Testing one SNP at a time does not fully realise the potential of genome-wide association studies to identify multiple causal variants, which is a plausible scenario for many complex diseases. We show that simultaneous analysis of the entire set of SNPs from a genome-wide study to identify the subset that best predicts disease outcome is now feasible, thanks to developments in stochastic search methods. We used a Bayesian-inspired penalised maximum likelihood approach in which every SNP can be considered for additive, dominant, and recessive contributions to disease risk. Posterior mode estimates were obtained for regression coefficients that were each assigned a prior with a sharp mode at zero. A non-zero coefficient estimate was interpreted as corresponding to a significant SNP. We investigated two prior distributions and show that the normal-exponential-gamma prior leads to improved SNP selection in comparison with single-SNP tests. We also derived an explicit approximation for type-I error that avoids the need to use permutation procedures. As well as genome-wide analyses, our method is well-suited to fine mapping with very dense SNP sets obtained from re-sequencing and/or imputation. It can accommodate quantitative as well as case-control phenotypes, covariate adjustment, and can be extended to search for interactions. Here, we demonstrate the power and empirical type-I error of our approach using simulated case-control data sets of up to 500 K SNPs, a real genome-wide data set of 300 K SNPs, and a sequence-based dataset, each of which can be analysed in a few hours on a desktop workstation.