Seeking unique and common biological themes in multiple gene lists or datasets: pathway pattern extraction pipeline for pathway-level comparative analysis

The Metabolic Landscape of RAS-Driven Cancers from biology to therapy

Our understanding of how the RAS protein family, and in particular mutant KRAS promote metabolic dysregulation in cancer cells has advanced significantly over the last decade. In this Review, we discuss the metabolic reprogramming mediated by oncogenic RAS in cancer, and elucidating the underlying mechanisms could translate to novel therapeutic opportunities to target metabolic vulnerabilities in RAS-driven cancers.

ssGSEA score-based Ras dependency indexes derived from gene expression data reveal potential Ras addiction mechanisms with possible clinical implications

For nearly a decade, the difficulties associated with both the determination and reproducibility of Ras-dependency indexes (RDIs) have limited their application and further delineation of the biology underlying Ras dependency. In this report, we describe the application of a computational single sample gene set enrichment analysis (ssGSEA) method to derive RDIs with gene expression data. The computationally derived RDIs across the Cancer Cell Line Encyclopedia (CCLE) cell lines show excellent agreement with the experimentally derived values and high correlation with a previous in-house siRNA effector node (siREN) study and external studies. Using EMT signature-derived RDIs and data from cell lines representing the extremes in RAS dependency, we identified enriched pathways distinguishing these classes, including the Fas signaling pathway and a putative Ras-independent pathway first identified in NK cells. Importantly, extension of the method to patient samples from The Cancer Genome Atlas (TCGA) showed the same consensus differential expression patterns for these two pathways across multiple tissue types. Last, the computational RDIs displayed a significant association with TCGA cancer patients' survival outcomes. Together, these lines of evidence confirm that our computationally derived RDIs faithfully represent a measure of Ras dependency in both cancer cell lines and patient samples. The application of such computational RDIs can provide insights into Ras biology and potential clinical applications.

MicroRNA-10b inhibition reduces E2F1-mediated transcription and miR-15/16 activity in glioblastoma

MicroRNA-10b (miR-10b) is commonly elevated in glioblastoma (GBM), while not expressed in normal brain tissues. Targeted inhibition of miR-10b has pleiotropic effects on GBM derived cell lines, it reduces GBM growth in animal models, but does not affect normal neurons and astrocytes. This data raises the possibility of developing miR-10b-targeting GBM therapy. However, the mechanisms contributing to miR-10b-mediated glioma cell survival and proliferation are unexplored. We found that inhibition of miR-10b has distinct effects on specific glioma cell lines. In cells expressing high levels of tumor suppressor p21WAF1/Cip1, it represses E2F1-mediated transcription, leading to down-regulation of multiple E2F1 target genes encoding for S-phase specific proteins, epigenetic modulators, and miRNAs (e.g. miR-15/16), and thereby stalling progression through the S-phase of cell cycle. Subsequently, miR-15/16 activities are reduced and many of their direct targets are de-repressed, including ubiquitin ligase FBXW7 that destabilizes Cyclin E. Conversely, GBM cells expressing low p21 level, or after p21 knock-down, exhibit weaker or no E2F1 response to miR-10b inhibition. Comparative analysis of The Cancer Genome Atlas revealed a strong correlation between miR-10b and multiple E2F target genes in GBM and low-grade glioma. Taken together, these findings indicate that miR-10b regulates E2F1-mediated transcription in GBM, in a p21-dependent fashion.

The Cancer Genome Atlas Analysis Predicts MicroRNA for Targeting Cancer Growth and Vascularization in Glioblastoma

Using in silico analysis of The Cancer Genome Atlas (TCGA), we identified microRNAs associated with glioblastoma (GBM) survival, and predicted their functions in glioma growth and progression. Inhibition of two "risky" miRNAs, miR-148a and miR-31, in orthotopic xenograft GBM mouse models suppressed tumor growth and thereby prolonged animal survival. Intracranial tumors treated with uncomplexed miR-148a and miR-31 antagomirs exhibited reduced proliferation, stem cell depletion, and normalized tumor vasculature. Growth-promoting functions of these two miRNAs were, in part, mediated by the common target, the factor inhibiting hypoxia-inducible factor 1 (FIH1), and the downstream pathways involving hypoxia-inducible factor HIF1α and Notch signaling. Therefore, miR-31 and miR-148a regulate glioma growth by maintaining tumor stem cells and their niche, and providing the tumor a way to activate angiogenesis even in a normoxic environment. This is the first study that demonstrates intratumoral uptake and growth-inhibiting effects of uncomplexed antagomirs in orthotopic glioma.

ChiNet uncovers rewired transcription subnetworks in tolerant yeast for advanced biofuels conversion

Analysis of rewired upstream subnetworks impacting downstream differential gene expression aids the delineation of evolving molecular mechanisms. Cumulative statistics based on conventional differential correlation are limited for subnetwork rewiring analysis since rewiring is not necessarily equivalent to change in correlation coefficients. Here we present a computational method ChiNet to quantify subnetwork rewiring by statistical heterogeneity that enables detection of potential genotype changes causing altered transcription regulation in evolving organisms. Given a differentially expressed downstream gene set, ChiNet backtracks a rewired upstream subnetwork from a super-network including gene interactions known to occur under various molecular contexts. We benchmarked ChiNet for its high accuracy in distinguishing rewired artificial subnetworks, in silico yeast transcription-metabolic subnetworks, and rewired transcription subnetworks for Candida albicans versus Saccharomyces cerevisiae, against two differential-correlation based subnetwork rewiring approaches. Then, using transcriptome data from tolerant S. cerevisiae strain NRRL Y-50049 and a wild-type intolerant strain, ChiNet identified 44 metabolic pathways affected by rewired transcription subnetworks anchored to major adaptively activated transcription factor genes YAP1, RPN4, SFP1 and ROX1, in response to toxic chemical challenges involved in lignocellulose-to-biofuels conversion. These findings support the use of ChiNet in rewiring analysis of subnetworks where differential interaction patterns resulting from divergent nonlinear dynamics abound.

miR-126 contributes to Parkinson's disease by dysregulating the insulin-like growth factor/phosphoinositide 3-kinase signaling

Dopamine (DA) neurons in sporadic Parkinson's disease (PD) display dysregulated gene expression networks and signaling pathways that are implicated in PD pathogenesis. Micro (mi)RNAs are regulators of gene expression, which could be involved in neurodegenerative diseases. We determined the miRNA profiles in laser microdissected DA neurons from postmortem sporadic PD patients' brains and age-matched controls. DA neurons had a distinctive miRNA signature and a set of miRNAs was dysregulated in PD. Bioinformatics analysis provided evidence for correlations of miRNAs with signaling pathways relevant to PD, including an association of miR-126 with insulin/IGF-1/PI3K signaling. In DA neuronal cell systems, enhanced expression of miR-126 impaired IGF-1 signaling and increased vulnerability to the neurotoxin 6-OHDA by downregulating factors in IGF-1/PI3K signaling, including its targets p85β, IRS-1, and SPRED1. Blocking of miR-126 function increased IGF-1 trophism and neuroprotection to 6-OHDA. Our data imply that elevated levels of miR-126 may play a functional role in DA neurons and in PD pathogenesis by downregulating IGF-1/PI3K/AKT signaling and that its inhibition could be a mechanism of neuroprotection.

An integrated understanding of the physiological response to elevated extracellular phosphate

Recent studies have suggested that changes in serum phosphate levels influence pathological states associated with aging such as cancer, bone metabolism, and cardiovascular function, even in individuals with normal renal function. The causes are only beginning to be elucidated but are likely a combination of endocrine, paracrine, autocrine, and cell autonomous effects. We have used an integrated quantitative biology approach, combining transcriptomics and proteomics to define a multi-phase, extracellular phosphate-induced, signaling network in pre-osteoblasts as well as primary human and mouse mesenchymal stromal cells. We identified a rapid mitogenic response stimulated by elevated phosphate that results in the induction of immediate early genes including c-fos. The mechanism of activation requires FGF receptor signaling followed by stimulation of N-Ras and activation of AP-1 and serum response elements. A distinct long-term response also requires FGF receptor signaling and results in N-Ras activation and expression of genes and secretion of proteins involved in matrix regulation, calcification, and angiogenesis. The late response is synergistically enhanced by addition of FGF23 peptide. The intermediate phase results in increased oxidative phosphorylation and ATP production and is necessary for the late response providing a functional link between the phases. Collectively, the results define elevated phosphate, as a mitogen and define specific mechanisms by which phosphate stimulates proliferation and matrix regulation. Our approach provides a comprehensive understanding of the cellular response to elevated extracellular phosphate, functionally connecting temporally coordinated signaling, transcriptional, and metabolic events with changes in long-term cell behavior.

Guanine holes are prominent targets for mutation in cancer and inherited disease

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G • C bp in the context of all 64 5'-NGNN-3' motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

IL-1R-MyD88 signaling in keratinocyte transformation and carcinogenesis

Constitutively active RAS plays a central role in the development of human cancer and is sufficient to induce tumors in two-stage skin carcinogenesis. RAS-mediated tumor formation is commonly associated with up-regulation of cytokines and chemokines that mediate an inflammatory response considered relevant to oncogenesis. In this study, we report that mice lacking IL-1R or MyD88 are less sensitive to topical skin carcinogenesis than their respective wild-type (WT) controls. MyD88(-/-) or IL-1R(-/-) keratinocytes expressing oncogenic RAS are hyperproliferative and fail to up-regulate proinflammatory genes or down-regulate differentiation markers characteristic of RAS-expressing WT keratinocytes. Although RAS-expressing MyD88(-/-) keratinocytes form only a few small tumors in orthotopic grafts, IL-1R-deficient RAS-expressing keratinocytes retain the ability to form tumors in orthotopic grafts. Using both genetic and pharmacological approaches, we find that the differentiation and proinflammatory effects of oncogenic RAS in keratinocytes require the establishment of an autocrine loop through IL-1α, IL-1R, and MyD88 leading to phosphorylation of IκBα and NF-κB activation. Blocking IL-1α-mediated NF-κB activation in RAS-expressing WT keratinocytes reverses the differentiation defect and inhibits proinflammatory gene expression. Collectively, these results demonstrate that MyD88 exerts a cell-intrinsic function in RAS-mediated transformation of keratinocytes.

MicroRNA-1 is a candidate tumor suppressor and prognostic marker in human prostate cancer

We previously reported that miR-1 is among the most consistently down-regulated miRs in primary human prostate tumors. In this follow-up study, we further corroborated this finding in an independent data set and made the novel observation that miR-1 expression is further reduced in distant metastasis and is a candidate predictor of disease recurrence. Moreover, we performed in vitro experiments to explore the tumor suppressor function of miR-1. Cell-based assays showed that miR-1 is epigenetically silenced in human prostate cancer. Overexpression of miR-1 in these cells led to growth inhibition and down-regulation of genes in pathways regulating cell cycle progression, mitosis, DNA replication/repair and actin dynamics. This observation was further corroborated with protein expression analysis and 3'-UTR-based reporter assays, indicating that genes in these pathways are either direct or indirect targets of miR-1. A gene set enrichment analysis revealed that the miR-1-mediated tumor suppressor effects are globally similar to those of histone deacetylase inhibitors. Lastly, we obtained preliminary evidence that miR-1 alters the cellular organization of F-actin and inhibits tumor cell invasion and filipodia formation. In conclusion, our findings indicate that miR-1 acts as a tumor suppressor in prostate cancer by influencing multiple cancer-related processes and by inhibiting cell proliferation and motility.

Evidence for gender-specific transcriptional profiles of nigral dopamine neurons in Parkinson disease

Background

Epidemiological data suggest that the male gender is one of the risks factors for the development of Parkinson Disease (PD). Also, differences in the clinical manifestation and the course of PD have been observed between males and females. However, little is known about the molecular aspects underlying gender-specificity in PD. To address this issue, we determined the gene expression profiles of male and female dopamine (DA) neurons in sporadic PD.

Methodology/Principal Findings

We analyzed Affymetrix-based microarrays on laser microdissected DA neurons from postmortem brains of sporadic PD patients and age-matched controls across genders. Pathway enrichment demonstrated that major cellular pathways involved in PD pathogenesis showed different patterns of deregulation between males and females with more prominent downregulation of genes related to oxidative phosporylation, apoptosis, synaptic transmission and transmission of nerve impulse in the male population. In addition, we found upregulation of gene products for metabolic processes and mitochondrial energy consumption in the age-matched male control neurons. On the single cell level, selected data validation using quantitative Real-Time (qRT)-PCR was consistent with microarray raw data and supported some of the observations from data analysis.

Conclusions/Significance

On the molecular level, our results provide evidence that the expression profiles of aged normal and PD midbrain DA neurons are gender-specific. The observed differences in the expression profiles suggest a disease bias of the male gender, which could be in concordance with clinical observations that the male gender represents a risk factor for sporadic PD. Validation of gene expression by qRT-PCR supported the microarray results, but also pointed to several caveats involved in data interpretation.

STARNET 2: a web-based tool for accelerating discovery of gene regulatory networks using microarray co-expression data

Background

Although expression microarrays have become a standard tool used by biologists, analysis of data produced by microarray experiments may still present challenges. Comparison of data from different platforms, organisms, and labs may involve complicated data processing, and inferring relationships between genes remains difficult.

Results

STARNET 2 is a new web-based tool that allows post hoc visual analysis of correlations that are derived from expression microarray data. STARNET 2 facilitates user discovery of putative gene regulatory networks in a variety of species (human, rat, mouse, chicken, zebrafish, Drosophila, C. elegans, S. cerevisiae, Arabidopsis and rice) by graphing networks of genes that are closely co-expressed across a large heterogeneous set of preselected microarray experiments. For each of the represented organisms, raw microarray data were retrieved from NCBI's Gene Expression Omnibus for a selected Affymetrix platform. All pairwise Pearson correlation coefficients were computed for expression profiles measured on each platform, respectively. These precompiled results were stored in a MySQL database, and supplemented by additional data retrieved from NCBI. A web-based tool allows user-specified queries of the database, centered at a gene of interest. The result of a query includes graphs of correlation networks, graphs of known interactions involving genes and gene products that are present in the correlation networks, and initial statistical analyses. Two analyses may be performed in parallel to compare networks, which is facilitated by the new HEATSEEKER module.

Conclusion

STARNET 2 is a useful tool for developing new hypotheses about regulatory relationships between genes and gene products, and has coverage for 10 species. Interpretation of the correlation networks is supported with a database of previously documented interactions, a test for enrichment of Gene Ontology terms, and heat maps of correlation distances that may be used to compare two networks. The list of genes in a STARNET network may be useful in developing a list of candidate genes to use for the inference of causal networks. The tool is freely available at , and does not require user registration.