Using multi-scale genomics to associate poorly annotated genes with rare diseases.

BACKGROUND: Next-generation sequencing (NGS) has significantly transformed the landscape of identifying disease-causing genes associated with genetic disorders. However, a substantial portion of sequenced patients remains undiagnosed. This may be attributed not only to the challenges posed by harder-to-detect variants, such as non-coding and structural variations but also to the existence of variants in genes not previously associated with the patient's clinical phenotype. This study introduces EvORanker, an algorithm that integrates unbiased data from 1,028 eukaryotic genomes to link mutated genes to clinical phenotypes. METHODS: EvORanker utilizes clinical data, multi-scale phylogenetic profiling, and other omics data to prioritize disease-associated genes. It was evaluated on solved exomes and simulated genomes, compared with existing methods, and applied to 6260 knockout genes with mouse phenotypes lacking human associations. Additionally, EvORanker was made accessible as a user-friendly web tool. RESULTS: In the analyzed exomic cohort, EvORanker accurately identified the "true" disease gene as the top candidate in 69% of cases and within the top 5 candidates in 95% of cases, consistent with results from the simulated dataset. Notably, EvORanker outperformed existing methods, particularly for poorly annotated genes. In the case of the 6260 knockout genes with mouse phenotypes, EvORanker linked 41% of these genes to observed human disease phenotypes. Furthermore, in two unsolved cases, EvORanker successfully identified DLGAP2 and LPCAT3 as disease candidates for previously uncharacterized genetic syndromes. CONCLUSIONS: We highlight clade-based phylogenetic profiling as a powerful systematic approach for prioritizing potential disease genes. Our study showcases the efficacy of EvORanker in associating poorly annotated genes to disease phenotypes observed in patients. The EvORanker server is freely available at https://ccanavati.shinyapps.io/EvORanker/ .

Systematic proximal mapping of the classical RAD51 paralogs unravel functionally and clinically relevant interactors for genome stability.

Homologous recombination (HR) plays an essential role in the maintenance of genome stability by promoting the repair of cytotoxic DNA double strand breaks (DSBs). More recently, the HR pathway has emerged as a core component of the response to replication stress, in part by protecting stalled replication forks from nucleolytic degradation. In that regard, the mammalian RAD51 paralogs (RAD51B, RAD51C, RAD51D, XRCC2, and XRCC3) have been involved in both HR-mediated DNA repair and collapsed replication fork resolution. Still, it remains largely obscure how they participate in both processes, thereby maintaining genome stability and preventing cancer development. To gain better insight into their contribution in cellulo, we mapped the proximal interactome of the classical RAD51 paralogs using the BioID approach. Aside from identifying the well-established BCDX2 and CX3 sub-complexes, the spliceosome machinery emerged as an integral component of our proximal mapping, suggesting a crosstalk between this pathway and the RAD51 paralogs. Furthermore, we noticed that factors involved RNA metabolic pathways are significantly modulated within the BioID of the classical RAD51 paralogs upon exposure to hydroxyurea (HU), pointing towards a direct contribution of RNA processing during replication stress. Importantly, several members of these pathways have prognostic potential in breast cancer (BC), where their RNA expression correlates with poorer patient outcome. Collectively, this study uncovers novel functionally relevant partners of the different RAD51 paralogs in the maintenance of genome stability that could be used as biomarkers for the prognosis of BC.

Multi-omics data integration analysis identifies the spliceosome as a key regulator of DNA double-strand break repair.

DNA repair by homologous recombination (HR) is critical for the maintenance of genome stability. Germline and somatic mutations in HR genes have been associated with an increased risk of developing breast (BC) and ovarian cancers (OvC). However, the extent of factors and pathways that are functionally linked to HR with clinical relevance for BC and OvC remains unclear. To gain a broader understanding of this pathway, we used multi-omics datasets coupled with machine learning to identify genes that are associated with HR and to predict their sub-function. Specifically, we integrated our phylogenetic-based co-evolution approach (CladePP) with 23 distinct genetic and proteomic screens that monitored, directly or indirectly, DNA repair by HR. This omics data integration analysis yielded a new database (HRbase) that contains a list of 464 predictions, including 76 gold standard HR genes. Interestingly, the spliceosome machinery emerged as one major pathway with significant cross-platform interactions with the HR pathway. We functionally validated 6 spliceosome factors, including the RNA helicase SNRNP200 and its co-factor SNW1. Importantly, their RNA expression correlated with BC/OvC patient outcome. Altogether, we identified novel clinically relevant DNA repair factors and delineated their specific sub-function by machine learning. Our results, supported by evolutionary and multi-omics analyses, suggest that the spliceosome machinery plays an important role during the repair of DNA double-strand breaks (DSBs).

CladeOScope: functional interactions through the prism of clade-wise co-evolution.

Mapping co-evolved genes via phylogenetic profiling (PP) is a powerful approach to uncover functional interactions between genes and to associate them with pathways. Despite many successful endeavors, the understanding of co-evolutionary signals in eukaryotes remains partial. Our hypothesis is that 'Clades', branches of the tree of life (e.g. primates and mammals), encompass signals that cannot be detected by PP using all eukaryotes. As such, integrating information from different clades should reveal local co-evolution signals and improve function prediction. Accordingly, we analyzed 1028 genomes in 66 clades and demonstrated that the co-evolutionary signal was scattered across clades. We showed that functionally related genes are frequently co-evolved in only parts of the eukaryotic tree and that clades are complementary in detecting functional interactions within pathways. We examined the non-homologous end joining pathway and the UFM1 ubiquitin-like protein pathway and showed that both demonstrated distinguished co-evolution patterns in specific clades. Our research offers a different way to look at co-evolution across eukaryotes and points to the importance of modular co-evolution analysis. We developed the 'CladeOScope' PP method to integrate information from 16 clades across over 1000 eukaryotic genomes and is accessible via an easy to use web server at http://cladeoscope.cs.huji.ac.il.

Reversal of diet-induced hepatic steatosis by peripheral CB1 receptor blockade in mice is p53/miRNA-22/SIRT1/PPARα dependent.

OBJECTIVE: The endocannabinoid (eCB) system is increasingly recognized as being crucially important in obesity-related hepatic steatosis. By activating the hepatic cannabinoid-1 receptor (CB1R), eCBs modulate lipogenesis and fatty acid oxidation. However, the underlying molecular mechanisms are largely unknown. METHODS: We combined unbiased bioinformatics techniques, mouse genetic manipulations, multiple pharmacological, molecular, and cellular biology approaches, and genomic sequencing to systematically decipher the role of the hepatic CB1R in modulating fat utilization in the liver and explored the downstream molecular mechanisms. RESULTS: Using an unbiased normalized phylogenetic profiling analysis, we found that the CB1R evolutionarily coevolves with peroxisome proliferator-activated receptor-alpha (PPARα), a key regulator of hepatic lipid metabolism. In diet-induced obese (DIO) mice, peripheral CB1R blockade (using AM6545) induced the reversal of hepatic steatosis and improved liver injury in WT, but not in PPARα-/- mice. The antisteatotic effect mediated by AM6545 in WT DIO mice was accompanied by increased hepatic expression and activity of PPARα as well as elevated hepatic levels of the PPARα-activating eCB-like molecules oleoylethanolamide and palmitoylethanolamide. Moreover, AM6545 was unable to rescue hepatic steatosis in DIO mice lacking liver sirtuin 1 (SIRT1), an upstream regulator of PPARα. Both of these signaling molecules were modulated by the CB1R as measured in hepatocytes exposed to lipotoxic conditions or treated with CB1R agonists in the absence/presence of AM6545. Furthermore, using microRNA transcriptomic profiling, we found that the CB1R regulated the hepatic expression, acetylation, and transcriptional activity of p53, resulting in the enhanced expression of miR-22, which was found to specifically target SIRT1 and PPARα. CONCLUSIONS: We provide strong evidence for a functional role of the p53/miR-22/SIRT1/PPARα signaling pathway in potentially mediating the antisteatotic effect of peripherally restricted CB1R blockade.

Titanium Tackles the Endoplasmic Reticulum: A First Genomic Study on a Titanium Anticancer Metallodrug.

PhenolaTi is an advanced non-toxic anticancer chemotherapy; this inert bis(phenolato)bis(alkoxo) Ti(IV) complex demonstrates the intriguing combination of high and wide efficacy with no detected toxicity in animals. Here we unravel the cellular pathways involved in its mechanism of action by a first genome study on Ti(IV)-treated cells, using an attuned RNA sequencing-based available technology. First, phenolaTi induced apoptosis and cell-cycle arrest at the G2/M phase in MCF7 cells. Second, the transcriptome of the treated cells was analyzed, identifying alterations in pathways relating to protein translation, DNA damage, and mitochondrial eruption. Unlike for common metallodrugs, electrophoresis assay showed no inhibition of DNA polymerase activity. Reduced in vitro cytotoxicity with added endoplasmic reticulum (ER) stress inhibitor supported the ER as a putative cellular target. Altogether, this paper reveals a distinct ER-related mechanism by the Ti(IV) anticancer coordination complex, paving the way for wider applicability of related techniques in mechanistic analyses of metallodrugs.

Optimization of co-evolution analysis through phylogenetic profiling reveals pathway-specific signals.

SUMMARY: The exponential growth in available genomic data is expected to reach full sequencing of a million genomes in the coming decade. Improving and developing methods to analyze these genomes and to reveal their utility is of major interest in a wide variety of fields, such as comparative and functional genomics, evolution and bioinformatics. Phylogenetic profiling is an established method for predicting functional interactions between proteins based on similarities in their evolutionary patterns across species. Proteins that function together (i.e. generate complexes, interact in the same pathways or improve adaptation to environmental niches) tend to show coordinated evolution across the tree of life. The normalized phylogenetic profiling (NPP) method takes into account minute changes in proteins across species to identify protein co-evolution. Despite the success of this method, it is still not clear what set of parameters is required for optimal use of co-evolution in predicting functional interactions. Moreover, it is not clear if pathway evolution or function should direct parameter choice. Here, we create a reliable and usable NPP construction pipeline. We explore the effect of parameter selection on functional interaction prediction using NPP from 1028 genomes, both separately and in various value combinations. We identify several parameter sets that optimize performance for pathways with certain biological annotation. This work reveals the importance of choosing the right parameters for optimized function prediction based on a biological context. AVAILABILITY AND IMPLEMENTATION: Source code and documentation are available on GitHub: https://github.com/iditam/CompareNPPs. CONTACT: yuvaltab@ekmd.huji.ac.il. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Author Correction: Clinical Implications of Sub-grouping HER2 Positive Tumors by Amplicon Structure and Co-amplified Genes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Clinical Implications of Sub-grouping HER2 Positive Tumors by Amplicon Structure and Co-amplified Genes.

ERBB2 amplification is a prognostic marker for aggressive tumors and a predictive marker for prolonged survival following treatment with HER2 inhibitors. We attempt to sub-group HER2+ tumors based on amplicon structures and co-amplified genes. We examined five HER2+ cell lines, three HER2+ xenographs and 57 HER2+ tumor tissues. ERBB2 amplification was analyzed using digital droplet PCR and low coverage whole genome sequencing. In some HER2+ tumors PPM1D, that encodes WIP1, is co-amplified. Cell lines were treated with HER2 and WIP1 inhibitors. We find that inverted duplication is the amplicon structure in the majority of HER2+ tumors. In patients suffering from an early stage disease the ERBB2 amplicon is composed of a single segment while in patients suffering from advanced cancer the amplicon is composed of several different segments. We find robust WIP1 inhibition in some HER2+ PPM1D amplified cell lines. Sub-grouping HER2+ tumors using low coverage whole genome sequencing identifies inverted duplications as the main amplicon structure and based on the number of segments, differentiates between local and advanced tumors. In addition, we found that we could determine if a tumor is a recurrent tumor or second primary tumor and identify co-amplified oncogenes that may serve as targets for therapy.

Mapping global and local coevolution across 600 species to identify novel homologous recombination repair genes.

The homologous recombination repair (HRR) pathway repairs DNA double-strand breaks in an error-free manner. Mutations in HRR genes can result in increased mutation rate and genomic rearrangements, and are associated with numerous genetic disorders and cancer. Despite intensive research, the HRR pathway is not yet fully mapped. Phylogenetic profiling analysis, which detects functional linkage between genes using coevolution, is a powerful approach to identify factors in many pathways. Nevertheless, phylogenetic profiling has limited predictive power when analyzing pathways with complex evolutionary dynamics such as the HRR. To map novel HRR genes systematically, we developed clade phylogenetic profiling (CladePP). CladePP detects local coevolution across hundreds of genomes and points to the evolutionary scale (e.g., mammals, vertebrates, animals, plants) at which coevolution occurred. We found that multiscale coevolution analysis is significantly more biologically relevant and sensitive to detect gene function. By using CladePP, we identified dozens of unrecognized genes that coevolved with the HRR pathway, either globally across all eukaryotes or locally in different clades. We validated eight genes in functional biological assays to have a role in DNA repair at both the cellular and organismal levels. These genes are expected to play a role in the HRR pathway and might lead to a better understanding of missing heredity in HRR-associated cancers (e.g., heredity breast and ovarian cancer). Our platform presents an innovative approach to predict gene function, identify novel factors related to different diseases and pathways, and characterize gene evolution.

Sex with no regrets: How sexual reproduction uses a no regret learning algorithm for evolutionary advantage.

The question of 'why sex' has long been a puzzle. The randomness of recombination, which potentially produces low fitness progeny, contradicts notions of fitness landscape hill climbing. We use the concept of evolution as an algorithm for learning unpredictable environments to provide a possible answer. While sex and asex both implement similar machine learning no-regret algorithms in the context of random samples that are small relative to a vast genotype space, the algorithm of sex constitutes a more efficient goal-directed walk through this space. Simulations indicate this gives sex an evolutionary advantage, even in stable, unchanging environments. Asexual populations rapidly reach a fitness plateau, but the learning aspect of the no-regret algorithm most often eventually boosts the fitness of sexual populations past the maximal viability of corresponding asexual populations. In this light, the randomness of sexual recombination is not a hindrance but a crucial component of the 'sampling for learning' algorithm of sexual reproduction.