MicroRNAs (miRNAs) are small regulatory RNAs that act by blocking the translation and increasing the degradation of target transcripts. MiRNAs play a critical role in many biological processes including development and differentiation and many studies have shown that major changes in miRNA levels occur in cancer. Since miRNAs degrade target messages, we used this property to develop a novel computational method aimed at determining the actual biological activity of miRNAs using variations in gene expression. Using the method described here, we quantified miRNA activity in papillary thyroid carcinoma and breast cancer, and found a strong and distinctive signal of increased global miRNA activity, embedded in the pertaining gene expression measurements. Interestingly, we found that in these two cancers, miRNA activity is globally increased, and is associated with a global downregulation of miRNA target genes. This downregulation of miRNA regulated genes is particularly noticeable for genes carrying multiple target sites for miRNAs. Among the miRNA-repressed genes, we found a significant enrichment of known tumor suppressors, thereby suggesting that the increased miRNA activity was indeed tumorigenic.

It is increasingly clear that transcription factors play versatile roles in turning genes "on" or "off" depending on cellular context via the various transcription complexes they form. This poses a major challenge in unraveling combinatorial transcription complex codes. Here we use the powerful genetics of Drosophila combined with microarray and bioinformatics analyses to tackle this challenge. The nuclear adaptor CHIP/LDB is a major developmental regulator capable of forming tissue-specific transcription complexes with various types of transcription factors and cofactors, making it a valuable model to study the intricacies of gene regulation. To date only few CHIP/LDB complexes target genes have been identified, and possible tissue-dependent crosstalk between these complexes has not been rigorously explored. SSDP proteins protect CHIP/LDB complexes from proteasome dependent degradation and are rate-limiting cofactors for these complexes. By using mutations in SSDP, we identified 189 down-stream targets of CHIP/LDB and show that these genes are enriched for the binding sites of APTEROUS (AP) and PANNIER (PNR), two well studied transcription factors associated with CHIP/LDB complexes. We performed extensive genetic screens and identified target genes that genetically interact with components of CHIP/LDB complexes in directing the development of the wings (28 genes) and thoracic bristles (23 genes). Moreover, by in vivo RNAi silencing we uncovered novel roles for two of the target genes, xbp1 and Gs-alpha, in early development of these structures. Taken together, our results suggest that loss of SSDP disrupts the normal balance between the CHIP-AP and the CHIP-PNR transcription complexes, resulting in down-regulation of CHIP-AP target genes and the concomitant up-regulation of CHIP-PNR target genes. Understanding the combinatorial nature of transcription complexes as presented here is crucial to the study of transcription regulation of gene batteries required for development.

The putative link between gene expression of brain regions and their neural connectivity patterns is a fundamental question in neuroscience. Here this question is addressed in the first large scale study of a prototypical mammalian rodent brain, using a combination of rat brain regional connectivity data with gene expression of the mouse brain. Remarkably, even though this study uses data from two different rodent species (due to the data limitations), we still find that the connectivity of the majority of brain regions is highly predictable from their gene expression levels-the outgoing (incoming) connectivity is successfully predicted for 73% (56%) of brain regions, with an overall fairly marked accuracy level of 0.79 (0.83). Many genes are found to play a part in predicting both the incoming and outgoing connectivity (241 out of the 500 top selected genes, p-value<1e-5). Reassuringly, the genes previously known from the literature to be involved in axon guidance do carry significant information about regional brain connectivity. Surveying the genes known to be associated with the pathogenesis of several brain disorders, we find that those associated with schizophrenia, autism and attention deficit disorder are the most highly enriched in the connectivity-related genes identified here. Finally, we find that the profile of functional annotation groups that are associated with regional connectivity in the rodent is significantly correlated with the annotation profile of genes previously found to determine neural connectivity in C. elegans (Pearson correlation of 0.24, p<1e-6 for the outgoing connections and 0.27, p<1e-5 for the incoming). Overall, the association between connectivity and gene expression in a specific extant rodent species' brain is likely to be even stronger than found here, given the limitations of current data.

A fundamental challenge in human health is the identification of disease-causing genes. Recently, several studies have tackled this challenge via a network-based approach, motivated by the observation that genes causing the same or similar diseases tend to lie close to one another in a network of protein-protein or functional interactions. However, most of these approaches use only local network information in the inference process and are restricted to inferring single gene associations. Here, we provide a global, network-based method for prioritizing disease genes and inferring protein complex associations, which we call PRINCE. The method is based on formulating constraints on the prioritization function that relate to its smoothness over the network and usage of prior information. We exploit this function to predict not only genes but also protein complex associations with a disease of interest. We test our method on gene-disease association data, evaluating both the prioritization achieved and the protein complexes inferred. We show that our method outperforms extant approaches in both tasks. Using data on 1,369 diseases from the OMIM knowledgebase, our method is able (in a cross validation setting) to rank the true causal gene first for 34% of the diseases, and infer 139 disease-related complexes that are highly coherent in terms of the function, expression and conservation of their member proteins. Importantly, we apply our method to study three multi-factorial diseases for which some causal genes have been found already: prostate cancer, alzheimer and type 2 diabetes mellitus. PRINCE's predictions for these diseases highly match the known literature, suggesting several novel causal genes and protein complexes for further investigation.

Synonymous mutations are considered to be "silent" as they do not affect protein sequence. However, different silent codons have different translation efficiency (TE), which raises the question to what extent such mutations are really neutral. We perform the first genome-wide study of natural selection operating on TE in recent human evolution, surveying 13,798 synonymous single nucleotide polymorphisms (SNPs) in 1,198 unrelated individuals from 11 populations. We find evidence for both negative and positive selection on TE, as measured based on differentiation in allele frequencies between populations. Notably, the likelihood of an SNP to be targeted by positive or negative selection is correlated with the magnitude of its effect on the TE of the corresponding protein. Furthermore, negative selection acting against changes in TE is more marked in highly expressed genes, highly interacting proteins, complex members, and regulatory genes. It is also more common in functional regions and in the initial segments of highly expressed genes. Positive selection targeting sites with a large effect on TE is stronger in lowly interacting proteins and in regulatory genes. Similarly, essential genes are enriched for negative TE selection while underrepresented for positive TE selection. Taken together, these results point to the significant role of TE as a selective force operating in humans and hence underscore the importance of considering silent SNPs in interpreting associations with complex human diseases. Testifying to this potential, we describe two synonymous SNPs that may have clinical implications in phenylketonuria and in Best's macular dystrophy due to TE differences between alleles.

There is a growing demand for genome-scale metabolic reconstructions for plants, fueled by the need to understand the metabolic basis of crop yield and by progress in genome and transcriptome sequencing. Methods are also required to enable the interpretation of plant transcriptome data to study how cellular metabolic activity varies under different growth conditions or even within different organs, tissues, and developmental stages. Such methods depend extensively on the accuracy with which genes have been mapped to the biochemical reactions in the plant metabolic pathways. Errors in these mappings lead to metabolic reconstructions with an inflated number of reactions and possible generation of unreliable metabolic phenotype predictions. Here we introduce a new evidence-based genome-scale metabolic reconstruction of maize, with significant improvements in the quality of the gene-reaction associations included within our model. We also present a new approach for applying our model to predict active metabolic genes based on transcriptome data. This method includes a minimal set of reactions associated with low expression genes to enable activity of a maximum number of reactions associated with high expression genes. We apply this method to construct an organ-specific model for the maize leaf, and tissue specific models for maize embryo and endosperm cells. We validate our models using fluxomics data for the endosperm and embryo, demonstrating an improved capacity of our models to fit the available fluxomics data. All models are publicly available via the DOE Systems Biology Knowledgebase and PlantSEED, and our new method is generally applicable for analysis transcript profiles from any plant, paving the way for further in silico studies with a wide variety of plant genomes.

Growth rate has long been considered one of the most valuable phenotypes that can be measured in cells. Aside from being highly accessible and informative in laboratory cultures, maximal growth rate is often a prime determinant of cellular fitness, and predicting phenotypes that underlie fitness is key to both understanding and manipulating life. Despite this, current methods for predicting microbial fitness typically focus on yields [e.g., predictions of biomass yield using GEnome-scale metabolic Models (GEMs)] or notably require many empirical kinetic constants or substrate uptake rates, which render these methods ineffective in cases where fitness derives most directly from growth rate. Here we present a new method for predicting cellular growth rate, termed SUMEX, which does not require any empirical variables apart from a metabolic network (i.e., a GEM) and the growth medium. SUMEX is calculated by maximizing the SUM of molar EXchange fluxes (hence SUMEX) in a genome-scale metabolic model. SUMEX successfully predicts relative microbial growth rates across species, environments, and genetic conditions, outperforming traditional cellular objectives (most notably, the convention assuming biomass maximization). The success of SUMEX suggests that the ability of a cell to catabolize substrates and produce a strong proton gradient enables fast cell growth. Easily applicable heuristics for predicting growth rate, such as what we demonstrate with SUMEX, may contribute to numerous medical and biotechnological goals, ranging from the engineering of faster-growing industrial strains, modeling of mixed ecological communities, and the inhibition of cancer growth.

Alterations in gene expression resulting from Alzheimer's disease have received considerable attention in recent years. Although expression has been investigated separately in whole brain tissue, in astrocytes and in neurons, a rigorous comparative study quantifying the relative utility of these sources in predicting the progression of Alzheimer's disease has been lacking. Here we analyze gene expression from neurons, astrocytes and whole tissues across different brain regions, and compare their ability to predict Alzheimer's disease progression by building pertaining classification models based on gene expression sets annotated to different biological processes. Remarkably, we find that predictions based on neuronal gene expression are significantly more accurate than those based on astrocyte or whole tissue expression. The findings explicate the central role of neurons, particularly as compared to glial cells, in the pathogenesis of Alzheimer's disease, and emphasize the importance of measuring gene expression in the most relevant (pathogenically 'proximal') single cell types.

Recent technological advances have enabled high-throughput measurements of protein-protein interactions in the cell, producing large protein interaction networks for various species at an ever-growing pace. However, common technologies like yeast two-hybrid may experience high rates of false positive detection. To combat false positive discoveries, a number of different methods have been recently developed that associate confidence scores with protein interactions. Here, we perform a rigorous comparative analysis and performance assessment among these different methods.

Metabolic alterations play an important role in cancer and yet, few metabolic cancer driver genes are known. Here we perform a combined genomic and metabolic modeling analysis searching for metabolic drivers of colorectal cancer. Our analysis predicts FUT9, which catalyzes the biosynthesis of Ley glycolipids, as a driver of advanced-stage colon cancer. Experimental testing reveals FUT9's complex dual role; while its knockdown enhances proliferation and migration in monolayers, it suppresses colon cancer cells expansion in tumorspheres and inhibits tumor development in a mouse xenograft models. These results suggest that FUT9's inhibition may attenuate tumor-initiating cells (TICs) that are known to dominate tumorspheres and early tumor growth, but promote bulk tumor cells. In agreement, we find that FUT9 silencing decreases the expression of the colorectal cancer TIC marker CD44 and the level of the OCT4 transcription factor, which is known to support cancer stemness. Beyond its current application, this work presents a novel genomic and metabolic modeling computational approach that can facilitate the systematic discovery of metabolic driver genes in other types of cancer.

Disease epidemiology during ageing shows a transition from cancer to degenerative chronic disorders as dominant contributors to mortality in the old. Nevertheless, it has remained unclear to what extent molecular signatures of ageing reflect this phenomenon. Here we report on the identification of a conserved transcriptomic signature of ageing based on gene expression data from four vertebrate species across four tissues. We find that ageing-associated transcriptomic changes follow trajectories similar to the transcriptional alterations observed in degenerative ageing diseases but are in opposite direction to the transcriptomic alterations observed in cancer. We confirm the existence of a similar antagonism on the genomic level, where a majority of shared risk alleles which increase the risk of cancer decrease the risk of chronic degenerative disorders and vice versa. These results reveal a fundamental trade-off between cancer and degenerative ageing diseases that sheds light on the pronounced shift in their epidemiology during ageing.

Identification of targeted therapies for TNBC is an urgent medical need. Using a drug combination screen reliant on synthetic lethal interactions, we identified clinically relevant combination therapies for different TNBC subtypes. Two drug combinations targeting the BET family were further explored. The first, targeting BET and CXCR2, is specific for mesenchymal TNBC and induces apoptosis, whereas the second, targeting BET and the proteasome, is effective for major TNBC subtypes and triggers ferroptosis. Ferroptosis was induced at low drug doses and was associated with increased cellular iron and decreased glutathione levels, concomitant with reduced levels of GPX4 and key glutathione biosynthesis genes. Further functional studies, analysis of clinical datasets and breast cancer specimens revealed a unique vulnerability of TNBC to ferroptosis inducers, enrichment of ferroptosis gene signature, and differential expression of key proteins that increase labile iron and decrease glutathione levels. This study identified potent combination therapies for TNBC and unveiled ferroptosis as a promising therapeutic strategy.

The pathological hallmarks of Alzheimer's disease (AD) are pathogenic oligomers and fibrils of misfolded amyloidogenic proteins (e.g., β-amyloid and hyper-phosphorylated tau in AD), which cause progressive loss of neurons in the brain and nervous system. Although deviations from normal protein glycosylation have been documented in AD, their role in disease pathology has been barely explored. Here our analysis of available expression data sets indicates that many glycosylation-related genes are differentially expressed in brains of AD patients compared with healthy controls. The robust differences found enabled us to predict the occurrence of AD with remarkable accuracy in a test cohort and identify a set of key genes whose expression determines this classification. We then studied in vivo the effect of reducing expression of homologs of 6 of these genes in transgenic Drosophila overexpressing human tau, a well-established invertebrate AD model. These experiments have led to the identification of glycosylation genes that may augment or ameliorate tauopathy phenotypes. Our results indicate that OstDelta, l(2)not and beta4GalT7 are tauopathy suppressors, whereas pgnat5 and CG33303 are enhancers, of tauopathy. These results suggest that specific alterations in protein glycosylation may play a causal role in AD etiology.

Modifier genes are believed to account for the clinical variability observed in many Mendelian disorders, but their identification remains challenging due to the limited availability of genomics data from large patient cohorts. Here, we present GENDULF (GENetic moDULators identiFication), one of the first methods to facilitate prediction of disease modifiers using healthy and diseased tissue gene expression data. GENDULF is designed for monogenic diseases in which the mechanism is loss of function leading to reduced expression of the mutated gene. When applied to cystic fibrosis, GENDULF successfully identifies multiple, previously established disease modifiers, including EHF, SLC6A14, and CLCA1. It is then utilized in spinal muscular atrophy (SMA) and predicts U2AF1 as a modifier whose low expression correlates with higher SMN2 pre-mRNA exon 7 retention. Indeed, knockdown of U2AF1 in SMA patient-derived cells leads to increased full-length SMN2 transcript and SMN protein expression. Taking advantage of the increasing availability of transcriptomic data, GENDULF is a novel addition to existing strategies for prediction of genetic disease modifiers, providing insights into disease pathogenesis and uncovering novel therapeutic targets.

With most cancer-related deaths resulting from metastasis, the development of new therapeutic approaches against metastatic colorectal cancer (mCRC) is essential to increasing patient survival. The metabolic adaptations that support mCRC remain undefined and their elucidation is crucial to identify potential therapeutic targets. Here, we employed a strategy for the rational identification of targetable metabolic vulnerabilities. This strategy involved first a thorough metabolic characterisation of same-patient-derived cell lines from primary colon adenocarcinoma (SW480), its lymph node metastasis (SW620) and a liver metastatic derivative (SW620-LiM2), and second, using a novel multi-omics integration workflow, identification of metabolic vulnerabilities specific to the metastatic cell lines. We discovered that the metastatic cell lines are selectively vulnerable to the inhibition of cystine import and folate metabolism, two key pathways in redox homeostasis. Specifically, we identified the system xCT and MTHFD1 genes as potential therapeutic targets, both individually and combined, for combating mCRC.

Most patients with advanced cancer eventually acquire resistance to targeted therapies, spurring extensive efforts to identify molecular events mediating therapy resistance. Many of these events involve synthetic rescue (SR) interactions, where the reduction in cancer cell viability caused by targeted gene inactivation is rescued by an adaptive alteration of another gene (the rescuer). Here, we perform a genome-wide in silico prediction of SR rescuer genes by analyzing tumor transcriptomics and survival data of 10,000 TCGA cancer patients. Predicted SR interactions are validated in new experimental screens. We show that SR interactions can successfully predict cancer patients' response and emerging resistance. Inhibiting predicted rescuer genes sensitizes resistant cancer cells to therapies synergistically, providing initial leads for developing combinatorial approaches to overcome resistance proactively. Finally, we show that the SR analysis of melanoma patients successfully identifies known mediators of resistance to immunotherapy and predicts novel rescuers.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 2019 has triggered an ongoing global pandemic of the severe pneumonia-like disease coronavirus disease 2019 (COVID-19)1. The development of a vaccine is likely to take at least 12-18 months, and the typical timeline for approval of a new antiviral therapeutic agent can exceed 10 years. Thus, repurposing of known drugs could substantially accelerate the deployment of new therapies for COVID-19. Here we profiled a library of drugs encompassing approximately 12,000 clinical-stage or Food and Drug Administration (FDA)-approved small molecules to identify candidate therapeutic drugs for COVID-19. We report the identification of 100 molecules that inhibit viral replication of SARS-CoV-2, including 21 drugs that exhibit dose-response relationships. Of these, thirteen were found to harbour effective concentrations commensurate with probable achievable therapeutic doses in patients, including the PIKfyve kinase inhibitor apilimod2-4 and the cysteine protease inhibitors MDL-28170, Z LVG CHN2, VBY-825 and ONO 5334. Notably, MDL-28170, ONO 5334 and apilimod were found to antagonize viral replication in human pneumocyte-like cells derived from induced pluripotent stem cells, and apilimod also demonstrated antiviral efficacy in a primary human lung explant model. Since most of the molecules identified in this study have already advanced into the clinic, their known pharmacological and human safety profiles will enable accelerated preclinical and clinical evaluation of these drugs for the treatment of COVID-19.

Precision oncology is gradually advancing into mainstream clinical practice, demonstrating significant survival benefits. However, eligibility and response rates remain limited in many cases, calling for better predictive biomarkers.

Cytotoxic T lymphocytes (CTL) play a crucial role in anti-cancer immunity. Progression of CTL to exhausted T lymphocytes (ETL) that overexpress inhibitory receptors can substantially decrease effector cytokines production and diminish cytolytic activity in tumor microenvironment (TME). However, while the activity levels of CTL and ETL are considered important determinants of immune checkpoint inhibitors (ICIs) response, it has been repeatedly observed that their predictive power of the latter is quite limited. Studying this conundrum on a large scale across the TCGA cohort, we find that ETL and CTL activity (estimated based on conventional gene signatures in the bulk tumor expression) is strongly positively correlated in most cancer types. We hypothesized that the limited predictive power of CTL activity might result from the high concordance of CTL and ETL activities, which mutually cancels out their individual antagonistic effects on ICI response.

Malignant mesothelioma is an aggressive cancer with limited treatment options and poor prognosis. A better understanding of mesothelioma genomics and transcriptomics could advance therapies. Here, we present a mesothelioma cohort of 122 patients along with their germline and tumor whole-exome and tumor RNA sequencing data as well as phenotypic and drug response information. We identify a 48-gene prognostic signature that is highly predictive of mesothelioma patient survival, including CCNB1, the expression of which is highly predictive of patient survival on its own. In addition, we analyze the transcriptomics data to study the tumor immune microenvironment and identify synthetic-lethality-based signatures predictive of response to therapy. This germline and somatic whole-exome sequencing as well as transcriptomics data from the same patient are a valuable resource to address important biological questions, including prognostic biomarkers and determinants of treatment response in mesothelioma.