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.

Epidermal homeostasis is tightly controlled by a balancing act of self-renewal or terminal differentiation of proliferating basal keratinocytes. An increase in DNA content as a consequence of a mitotic block is a recognized mechanism underlying keratinocyte differentiation, but the molecular mechanisms involved in this process are not yet fully understood. Using cultured primary keratinocytes, here we report that the expression of the mammalian microtubule and kinetochore-associated protein Clasp2 is intimately associated with the basal proliferative makeup of keratinocytes, and its deficiency leads to premature differentiation. Clasp2-deficient keratinocytes exhibit increased centrosomal numbers and numerous mitotic alterations, including multipolar spindles and chromosomal misalignments that overall result in mitotic stress and a high DNA content. Such mitotic block prompts premature keratinocyte differentiation in a p53-dependent manner in the absence of cell death. Our findings reveal a new role for Clasp2 in governing keratinocyte undifferentiated features and highlight the presence of surveillance mechanisms that prevent cell cycle entry in cells that have alterations in the DNA content.

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.

4-Hydroxy-tamoxifen (4OHT) triggers Cre-mediated K-Ras removal in [H-Ras-/-; N-Ras-/-; K-Ras lox/lox; RERT ert/ert] fibroblasts, generating growth-arrested "Rasless" MEFs which are able to recover their proliferative ability after ectopic expression of Ras oncoproteins or constitutively active BRAF or MEK1.

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.

HMG-box proteins, including Sox/SRY (Sox) and TCF/LEF1 (TCF) family members, bind DNA via their HMG-box. This binding, however, is relatively weak and both Sox and TCF factors employ distinct mechanisms for enhancing their affinity and specificity for DNA. Here we report that Capicua (CIC), an HMG-box transcriptional repressor involved in Ras/MAPK signaling and cancer progression, employs an additional distinct mode of DNA binding that enables selective recognition of its targets. We find that, contrary to previous assumptions, the HMG-box of CIC does not bind DNA alone but instead requires a distant motif (referred to as C1) present at the C-terminus of all CIC proteins. The HMG-box and C1 domains are both necessary for binding specific TGAATGAA-like sites, do not function via dimerization, and are active in the absence of cofactors, suggesting that they form a bipartite structure for sequence-specific binding to DNA. We demonstrate that this binding mechanism operates throughout Drosophila development and in human cells, ensuring specific regulation of multiple CIC targets. It thus appears that HMG-box proteins generally depend on auxiliary DNA binding mechanisms for regulating their appropriate genomic targets, but that each sub-family has evolved unique strategies for this purpose. Finally, the key role of C1 in DNA binding also explains the fact that this domain is a hotspot for inactivating mutations in oligodendroglioma and other tumors, while being preserved in oncogenic CIC-DUX4 fusion chimeras associated to Ewing-like sarcomas.

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.

Drug combinations are key to circumvent resistance mechanisms compromising response to single anti-cancer targeted therapies. The implementation of combinatorial approaches involving MEK1/2 or KRASG12C inhibitors in the context of KRAS-mutated lung cancers focuses fundamentally on targeting KRAS proximal activators or effectors. However, the antitumor effect is highly determined by compensatory mechanisms arising in defined cell types or tumor subgroups. A potential strategy to find drug combinations targeting a larger fraction of KRAS-mutated lung cancers may capitalize on the common, distal gene expression output elicited by oncogenic KRAS. By integrating a signature-driven drug repurposing approach with a pairwise pharmacological screen, here we show synergistic drug combinations consisting of multi-tyrosine kinase PKC inhibitors together with MEK1/2 or KRASG12C inhibitors. Such combinations elicit a cytotoxic response in both in vitro and in vivo models, which in part involves inhibition of the PKC inhibitor target AURKB. Proteome profiling links dysregulation of MYC expression to the effect of both PKC inhibitor-based drug combinations. Furthermore, MYC overexpression appears as a resistance mechanism to MEK1/2 and KRASG12C inhibitors. Our study provides a rational framework for selecting drugs entering combinatorial strategies and unveils MEK1/2- and KRASG12C-based therapies for lung cancer.

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.

Infections are a significant cause of morbidity and mortality in patients with chronic lymphocytic leukemia (CLL). The pathogenesis of infections is multifactorial and includes hypogammaglobulinemia, conventional therapy with alkylating drugs, and recently, purine analogs and mAb-associated T cells. Patients without these risk factors also suffer from infections, although the mechanism remains unknown. In a cohort of 70 patients with CLL, we demonstrated that their monocytes were locked into a refractory state and were unable to mount a classic inflammatory response to pathogens. In addition, they exhibited the primary features of endotoxin tolerance, including low cytokine production, high phagocytic activity, and impaired Ag presentation. The involvement of miR-146a in this phenomenon was suspected. We found miR-146a target genes, such as IRAK1 and TRAF6, were manifestly downregulated. Our study provides a new explanation for infections in patients with CLL and describes a cross-tolerance between endotoxins and tumors.

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 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.

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 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.

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.

The ongoing COVID-19 pandemic is one of the biggest health challenges of recent decades. Among the causes of mortality triggered by SARS-CoV-2 infection, the development of an inflammatory "cytokine storm" (CS) plays a determinant role. Here, we used transcriptomic data from the bronchoalveolar lavage fluid (BALF) of COVID-19 patients undergoing a CS to obtain gene-signatures associated to this pathology. Using these signatures, we interrogated the Connectivity Map (CMap) dataset that contains the effects of over 5000 small molecules on the transcriptome of human cell lines, and looked for molecules which effects on transcription mimic or oppose those of the CS. As expected, molecules that potentiate immune responses such as PKC activators are predicted to worsen the CS. In addition, we identified the negative regulation of female hormones among pathways potentially aggravating the CS, which helps to understand the gender-related differences in COVID-19 mortality. Regarding drugs potentially counteracting the CS, we identified glucocorticoids as a top hit, which validates our approach as this is the primary treatment for this pathology. Interestingly, our analysis also reveals a potential effect of MEK inhibitors in reverting the COVID-19 CS, which is supported by in vitro data that confirms the anti-inflammatory properties of these compounds.