Recent observations point towards the existence of a large number of neighborhoods composed of functionally-related gene modules that lie together in the genome. This local component in the distribution of the functionality across chromosomes is probably affecting the own chromosomal architecture by limiting the possibilities in which genes can be arranged and distributed across the genome. As a direct consequence of this fact it is therefore presumable that diseases such as cancer, harboring DNA copy number alterations (CNAs), will have a symptomatology strongly dependent on modules of functionally-related genes rather than on a unique "important" gene.

A quarter of all solid tumors harbor KRAS oncogenes. Yet, no selective drugs have been approved to treat these malignancies. Genetic interrogation of the MAPK pathway revealed that systemic ablation of MEK or ERK kinases in adult mice prevent tumor development but are unacceptably toxic. Here, we demonstrate that ablation of c-RAF expression in advanced tumors driven by KrasG12V/Trp53 mutations leads to significant tumor regression with no detectable appearance of resistance mechanisms. Tumor regression results from massive apoptosis. Importantly, systemic abrogation of c-RAF expression does not inhibit canonical MAPK signaling, hence, resulting in limited toxicities. These results are of significant relevance for the design of therapeutic strategies to treat K-RAS mutant cancers.

In silico drug prescription tools for precision cancer medicine can match molecular alterations with tailored candidate treatments. These methodologies require large and well-annotated datasets to systematically evaluate their performance, but this is currently constrained by the lack of complete patient clinicopathological data. Moreover, in silico drug prescription performance could be improved by integrating additional tumour information layers like intra-tumour heterogeneity (ITH) which has been related to drug response and tumour progression. PanDrugs is an in silico drug prescription method which prioritizes anticancer drugs combining both biological and clinical evidence. We have systematically evaluated PanDrugs in the Genomic Data Commons repository (GDC). Our results showed that PanDrugs is able to establish an a priori stratification of cancer patients treated with Epidermal Growth Factor Receptor (EGFR) inhibitors. Patients labelled as responders according to PanDrugs predictions showed a significantly increased overall survival (OS) compared to non-responders. PanDrugs was also able to suggest alternative tailored treatments for non-responder patients. Additionally, PanDrugs usefulness was assessed considering spatial and temporal ITH in cancer patients and showed that ITH can be approached therapeutically proposing drugs or combinations potentially capable of targeting the clonal diversity. In summary, this study is a proof of concept where PanDrugs predictions have been correlated to OS and can be useful to manage ITH in patients while increasing therapeutic options and demonstrating its clinical utility.

FBXW7 is one of the most frequently mutated tumor suppressors, deficiency of which has been associated with resistance to some anticancer therapies. Through bioinformatics and genome-wide CRISPR screens, we here reveal that FBXW7 deficiency leads to multidrug resistance (MDR). Proteomic analyses found an upregulation of mitochondrial factors as a hallmark of FBXW7 deficiency, which has been previously linked to chemotherapy resistance. Despite this increased expression of mitochondrial factors, functional analyses revealed that mitochondria are under stress, and genetic or chemical targeting of mitochondria is preferentially toxic for FBXW7-deficient cells. Mechanistically, the toxicity of therapies targeting mitochondrial translation such as the antibiotic tigecycline relates to the activation of the integrated stress response (ISR) in a GCN2 kinase-dependent manner. Furthermore, the discovery of additional drugs that are toxic for FBXW7-deficient cells showed that all of them unexpectedly activate a GCN2-dependent ISR regardless of their accepted mechanism of action. Our study reveals that while one of the most frequent mutations in cancer reduces the sensitivity to the vast majority of available therapies, it renders cells vulnerable to ISR-activating drugs.

The role of heterochromatin in cell fate specification during development is unclear. We demonstrate that loss of the lysine 9 of histone H3 (H3K9) methyltransferase G9a in the mammary epithelium results in de novo chromatin opening, aberrant formation of the mammary ductal tree, impaired stem cell potential, disrupted intraductal polarity, and loss of tissue function. G9a loss derepresses long terminal repeat (LTR) retroviral sequences (predominantly the ERVK family). Transcriptionally activated endogenous retroviruses generate double-stranded DNA (dsDNA) that triggers an antiviral innate immune response, and knockdown of the cytosolic dsDNA sensor Aim2 in G9a knockout (G9acKO) mammary epithelium rescues mammary ductal invasion. Mammary stem cell transplantation into immunocompromised or G9acKO-conditioned hosts shows partial dependence of the G9acKO mammary morphological defects on the inflammatory milieu of the host mammary fat pad. Thus, altering the chromatin accessibility of retroviral elements disrupts mammary gland development and stem cell activity through both cell-autonomous and non-autonomous mechanisms.

Gene Expression Profile Analysis Suite (GEPAS) is one of the most complete and extensively used web-based packages for microarray data analysis. During its more than 5 years of activity it has continuously been updated to keep pace with the state-of-the-art in the changing microarray data analysis arena. GEPAS offers diverse analysis options that include well established as well as novel algorithms for normalization, gene selection, class prediction, clustering and functional profiling of the experiment. New options for time-course (or dose-response) experiments, microarray-based class prediction, new clustering methods and new tests for differential expression have been included. The new pipeliner module allows automating the execution of sequential analysis steps by means of a simple but powerful graphic interface. An extensive re-engineering of GEPAS has been carried out which includes the use of web services and Web 2.0 technology features, a new user interface with persistent sessions and a new extended database of gene identifiers. GEPAS is nowadays the most quoted web tool in its field and it is extensively used by researchers of many countries and its records indicate an average usage rate of 500 experiments per day. GEPAS, is available at http://www.gepas.org.

The Gene Expression Profile Analysis Suite (GEPAS) has been running for more than four years. During this time it has evolved to keep pace with the new interests and trends in the still changing world of microarray data analysis. GEPAS has been designed to provide an intuitive although powerful web-based interface that offers diverse analysis options from the early step of preprocessing (normalization of Affymetrix and two-colour microarray experiments and other preprocessing options), to the final step of the functional annotation of the experiment (using Gene Ontology, pathways, PubMed abstracts etc.), and include different possibilities for clustering, gene selection, class prediction and array-comparative genomic hybridization management. GEPAS is extensively used by researchers of many countries and its records indicate an average usage rate of 400 experiments per day. The web-based pipeline for microarray gene expression data, GEPAS, is available at http://www.gepas.org.

Epidemiological studies indicate that patients suffering from Alzheimer's disease have a lower risk of developing lung cancer, and suggest a higher risk of developing glioblastoma. Here we explore the molecular scenarios that might underlie direct and inverse co-morbidities between these diseases. Transcriptomic meta-analyses reveal significant numbers of genes with inverse patterns of expression in Alzheimer's disease and lung cancer, and with similar patterns of expression in Alzheimer's disease and glioblastoma. These observations support the existence of molecular substrates that could at least partially account for these direct and inverse co-morbidity relationships. A functional analysis of the sets of deregulated genes points to the immune system, up-regulated in both Alzheimer's disease and glioblastoma, as a potential link between these two diseases. Mitochondrial metabolism is regulated oppositely in Alzheimer's disease and lung cancer, indicating that it may be involved in the inverse co-morbidity between these diseases. Finally, oxidative phosphorylation is a good candidate to play a dual role by decreasing or increasing the risk of lung cancer and glioblastoma in Alzheimer's disease.

A high percentage of patients with brain metastases frequently develop neurocognitive symptoms; however, understanding how brain metastasis co-opts the function of neuronal circuits beyond a tumor mass effect remains unknown. We report a comprehensive multidimensional modeling of brain functional analyses in the context of brain metastasis. By testing different preclinical models of brain metastasis from various primary sources and oncogenic profiles, we dissociated the heterogeneous impact on local field potential oscillatory activity from cortical and hippocampal areas that we detected from the homogeneous inter-model tumor size or glial response. In contrast, we report a potential underlying molecular program responsible for impairing neuronal crosstalk by scoring the transcriptomic and mutational profiles in a model-specific manner. Additionally, measurement of various brain activity readouts matched with machine learning strategies confirmed model-specific alterations that could help predict the presence and subtype of metastasis.

Chronic lymphocytic leukaemia is the most prevalent leukaemia in Western countries. It is an incurable disease characterized by a highly variable clinical course. Chronic lymphocytic leukaemia is an ideal model for studying clonal heterogeneity and dynamics during cancer progression, response to therapy and/or relapse because the disease usually develops over several years. Here we report an analysis by deep sequencing of sequential samples taken at different times from the affected organs of two patients with 12- and 7-year disease courses, respectively. One of the patients followed a linear pattern of clonal evolution, acquiring and selecting new mutations in response to salvage therapy and/or allogeneic transplantation, while the other suffered loss of cellular tumoral clones during progression and histological transformation.

The Gene Expression Profile Analysis Suite, GEPAS, has been running for more than three years. With >76,000 experiments analysed during the last year and a daily average of almost 300 analyses, GEPAS can be considered a well-established and widely used platform for gene expression microarray data analysis. GEPAS is oriented to the analysis of whole series of experiments. Its design and development have been driven by the demands of the biomedical community, probably the most active collective in the field of microarray users. Although clustering methods have obviously been implemented in GEPAS, our interest has focused more on methods for finding genes differentially expressed among distinct classes of experiments or correlated to diverse clinical outcomes, as well as on building predictors. There is also a great interest in CGH-arrays which fostered the development of the corresponding tool in GEPAS: InSilicoCGH. Much effort has been invested in GEPAS for developing and implementing efficient methods for functional annotation of experiments in the proper statistical framework. Thus, the popular FatiGO has expanded to a suite of programs for functional annotation of experiments, including information on transcription factor binding sites, chromosomal location and tissues. The web-based pipeline for microarray gene expression data, GEPAS, is available at http://www.gepas.org.

We present a new version of Babelomics, a complete suite of web tools for functional analysis of genome-scale experiments, with new and improved tools. New functionally relevant terms have been included such as CisRed motifs or bioentities obtained by text-mining procedures. An improved indexing has considerably speeded up several of the modules. An improved version of the FatiScan method for studying the coordinate behaviour of groups of functionally related genes is presented, along with a similar tool, the Gene Set Enrichment Analysis. Babelomics is now more oriented to test systems biology inspired hypotheses. Babelomics can be found at http://www.babelomics.org.

Large-sequencing cancer genome projects have shown that tumors have thousands of molecular alterations and their frequency is highly heterogeneous. In such scenarios, physicians and oncologists routinely face lists of cancer genomic alterations where only a minority of them are relevant biomarkers to drive clinical decision-making. For this reason, the medical community agrees on the urgent need of methodologies to establish the relevance of tumor alterations, assisting in genomic profile interpretation, and, more importantly, to prioritize those that could be clinically actionable for cancer therapy.

Functional profiling methods have been extensively used in the context of high-throughput experiments and, in particular, in microarray data analysis. Such methods use available biological information to define different types of functional gene modules (e.g. gene ontology -GO-, KEGG pathways, etc.) whose representation in a pre-defined list of genes is further studied. In the most popular type of microarray experimental designs (e.g. up- or down-regulated genes, clusters of co-expressing genes, etc.) or in other genomic experiments (e.g. Chip-on-chip, epigenomics, etc.) these lists are composed by genes with a high degree of co-expression. Therefore, an implicit assumption in the application of functional profiling methods within this context is that the genes corresponding to the modules tested are effectively defining sets of co-expressing genes. Nevertheless not all the functional modules are biologically coherent entities in terms of co-expression, which will eventually hinder its detection with conventional methods of functional enrichment.

Contrarily to the traditional view in which only one or a few key genes were supposed to be the causative factors of diseases, we discuss the importance of considering groups of functionally related genes in the study of pathologies characterised by chromosomal copy number alterations. Recent observations have reported the existence of regions in higher eukaryotic chromosomes (including humans) containing genes of related function that show a high degree of coregulation. Copy number alterations will consequently affect to clusters of functionally related genes, which will be the final causative agents of the diseased phenotype, in many cases. Therefore, we propose that the functional profiling of the regions affected by copy number alterations must be an important aspect to take into account in the understanding of this type of pathologies. To illustrate this, we present an integrated study of DNA copy number variations, gene expression along with the functional profiling of chromosomal regions in a case of multiple myeloma.

We present the ISACGH, a web-based system that allows for the combination of genomic data with gene expression values and provides different options for functional profiling of the regions found. Several visualization options offer a convenient representation of the results. Different efficient methods for accurate estimation of genomic copy number from array-CGH hybridization data have been included in the program. Moreover, the connection to the gene expression analysis package GEPAS allows the use of different facilities for data pre-processing and analysis. A DAS server allows exporting the results to the Ensembl viewer where contextual genomic information can be obtained. The program is freely available at: http://isacgh.bioinfo.cipf.es or within http://www.gepas.org.

We present a new version of Babelomics, a complete suite of web tools for the functional profiling of genome scale experiments, with new and improved methods as well as more types of functional definitions. Babelomics includes different flavours of conventional functional enrichment methods as well as more advanced gene set analysis methods that makes it a unique tool among the similar resources available. In addition to the well-known functional definitions (GO, KEGG), Babelomics includes new ones such as Biocarta pathways or text mining-derived functional terms. Regulatory modules implemented include transcriptional control (Transfac, CisRed) and other levels of regulation such as miRNA-mediated interference. Moreover, Babelomics allows for sub-selection of terms in order to test more focused hypothesis. Also gene annotation correspondence tables can be imported, which allows testing with user-defined functional modules. Finally, a tool for the 'de novo' functional annotation of sequences has been included in the system. This allows using yet unannotated organisms in the program. Babelomics has been extensively re-engineered and now it includes the use of web services and Web 2.0 technology features, a new user interface with persistent sessions and a new extended database of gene identifiers. Babelomics is available at http://www.babelomics.org.

Bone marrow (BM) stromal cells provide the regulatory framework for hematopoiesis and contribute to developmental stage-specific niches, such as those preserving hematopoietic stem cells. Despite advances in our understanding of stromal function, little is known about the transcriptional changes that this compartment undergoes throughout lifespan and during adaptation to stress. Using RNA sequencing, we perform transcriptional analyses of four principal stromal subsets, namely CXCL12-abundant reticular, platelet-derived growth factor receptor (PDGFR)-α+Sca1+, sinusoidal, and arterial endothelial cells, from early postnatal, adult, and aged mice. Our data reveal (1) molecular fingerprints defining cell-specific anatomical and functional features, (2) a radical reprogramming of pro-hematopoietic, immune, and matrisomic transcriptional programs during the transition from juvenile stages to adulthood, and (3) the aging-driven progressive upregulation of pro-inflammatory gene expression in stroma. We further demonstrate that transcriptomic pathways elicited in vivo by prototypic microbial molecules are largely recapitulated during aging, thereby supporting the inflammatory basis of age-related adaptations of BM hematopoietic function.

Acute T-cell lymphoblastic leukaemia (T-ALL) is an aggressive disorder derived from immature thymocytes. The variability observed in clinical responses on this type of tumours to treatments, the high toxicity of current protocols and the poor prognosis of patients with relapse or refractory make it urgent to find less toxic and more effective therapies in the context of a personalized medicine of precision.

Five-year survival for pancreatic ductal adenocarcinoma (PDAC) patients remains below 7% due to the lack of effective treatments. Here, we report that combined ablation of EGFR and c-RAF expression results in complete regression of a significant percentage of PDAC tumors driven by Kras/Trp53 mutations in genetically engineered mice. Moreover, systemic elimination of these targets induces toxicities that are well tolerated. Response to this targeted therapy correlates with transcriptional profiles that resemble those observed in human PDACs. Finally, inhibition of EGFR and c-RAF expression effectively blocked tumor progression in nine independent patient-derived xenografts carrying KRAS and TP53 mutations. These results open the door to the development of targeted therapies for PDAC patients.