Mucosal Genes Encoding Clock, Inflammation and Their Mutual Regulators Are Disrupted in Pediatric Patients with Active Ulcerative Colitis.

Patients with active ulcerative colitis (UC) display a misalignment of the circadian clock, which plays a vital role in various immune functions. Our aim was to characterize the expression of clock and inflammation genes, and their mutual regulatory genes in treatment-naïve pediatric patients with UC. Using the Inflammatory Bowel Disease Transcriptome and Metatranscriptome Meta-Analysis (IBD TaMMA) platform and R algorithms, we analyzed rectal biopsy transcriptomic data from two cohorts (206 patients with UC vs. 20 healthy controls from the GSE-109142 study, and 43 patients with UC vs. 55 healthy controls from the GSE-117993 study). We compared gene expression levels and correlation of clock genes (BMAL1, CLOCK, PER1, PER2, CRY1, CRY2), inflammatory genes (IκB, IL10, NFκB1, NFκB2, IL6, TNFα) and their mutual regulatory genes (RORα, RORγ, REV-ERBα, PGC1α, PPARα, PPARγ, AMPK, SIRT1) in patients with active UC and healthy controls. The clock genes BMAL1, CLOCK, PER1 and CRY1 and the inflammatory genes IκB, IL10, NFκB1, NFκB2, IL6 and TNFα were significantly upregulated in patients with active UC. The genes encoding the mutual regulators RORα, RORγ, PGC1α, PPARα and PPARγ were significantly downregulated in patients with UC. A uniform pattern of gene expression was found in healthy controls compared to the highly variable expression pattern in patients with UC. Among the healthy controls, inflammatory genes were positively correlated with clock genes and they all showed reduced expression. The difference in gene expression levels was associated with disease severity and endoscopic score but not with histological score. In patients with active UC, clock gene disruption is associated with abnormal mucosal immune response. Disrupted expression of genes encoding clock, inflammation and their mutual regulators together may play a role in active UC.

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

Developmental and temporal changes in petunia petal transcriptome reveal scent-repressing plant-specific RING-kinase-WD40 protein.

In moth-pollinated petunias, production of floral volatiles initiates when the flower opens and occurs rhythmically during the day, for optimal flower-pollinator interaction. To characterize the developmental transcriptomic response to time of day, we generated RNA-Seq databases for corollas of floral buds and mature flowers in the morning and in the evening. Around 70% of transcripts accumulating in petals demonstrated significant changes in expression levels in response to the flowers' transition from a 4.5-cm bud to a flower 1 day postanthesis (1DPA). Overall, 44% of the petal transcripts were differentially expressed in the morning vs. evening. Morning/evening changes were affected by flower developmental stage, with a 2.5-fold larger transcriptomic response to daytime in 1DPA flowers compared to buds. Analyzed genes known to encode enzymes in volatile organic compound biosynthesis were upregulated in 1DPA flowers vs. buds-in parallel with the activation of scent production. Based on analysis of global changes in the petal transcriptome, PhWD2 was identified as a putative scent-related factor. PhWD2 is a protein that is uniquely present in plants and has a three-domain structure: RING-kinase-WD40. Suppression of PhWD2 (termed UPPER - Unique Plant PhEnylpropanoid Regulator) resulted in a significant increase in the levels of volatiles emitted from and accumulated in internal pools, suggesting that it is a negative regulator of petunia floral scent production.

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.

Transcriptome changes in DM1 patients' tissues are governed by the RNA interference pathway.

Myotonic dystrophy type 1 (DM1) is a multisystemic disease caused by pathogenic expansions of CTG repeats. The expanded repeats are transcribed to long RNA and induce cellular toxicity. Recent studies suggest that the CUG repeats are processed by the RNA interference (RNAi) pathway to generate small interfering repeated RNA (siRNA). However, the effects of the CTG repeat-derived siRNAs remain unclear. We hypothesize that the RNAi machinery in DM1 patients generates distinct gene expression patterns that determine the disease phenotype in the individual patient. The abundance of genes with complementary repeats that are targeted by siRNAs in each tissue determines the way that the tissue is affected in DM1. We integrated and analyzed published transcriptome data from muscle, heart, and brain biopsies of DM1 patients, and revealed shared, characteristic changes that correlated with disease phenotype. These signatures are overrepresented by genes and transcription factors bearing endogenous CTG/CAG repeats and are governed by aberrant activity of the RNAi machinery, miRNAs, and a specific gain-of-function of the CTG repeats. Computational analysis of the DM1 transcriptome enhances our understanding of the complex pathophysiology of the disease and may reveal a path for cure.

An Exercise-Induced Metabolic Shield in Distant Organs Blocks Cancer Progression and Metastatic Dissemination.

Exercise prevents cancer incidence and recurrence, yet the underlying mechanism behind this relationship remains mostly unknown. Here we report that exercise induces the metabolic reprogramming of internal organs that increases nutrient demand and protects against metastatic colonization by limiting nutrient availability to the tumor, generating an exercise-induced metabolic shield. Proteomic and ex vivo metabolic capacity analyses of murine internal organs revealed that exercise induces catabolic processes, glucose uptake, mitochondrial activity, and GLUT expression. Proteomic analysis of routinely active human subject plasma demonstrated increased carbohydrate utilization following exercise. Epidemiologic data from a 20-year prospective study of a large human cohort of initially cancer-free participants revealed that exercise prior to cancer initiation had a modest impact on cancer incidence in low metastatic stages but significantly reduced the likelihood of highly metastatic cancer. In three models of melanoma in mice, exercise prior to cancer injection significantly protected against metastases in distant organs. The protective effects of exercise were dependent on mTOR activity, and inhibition of the mTOR pathway with rapamycin treatment ex vivo reversed the exercise-induced metabolic shield. Under limited glucose conditions, active stroma consumed significantly more glucose at the expense of the tumor. Collectively, these data suggest a clash between the metabolic plasticity of cancer and exercise-induced metabolic reprogramming of the stroma, raising an opportunity to block metastasis by challenging the metabolic needs of the tumor. SIGNIFICANCE: Exercise protects against cancer progression and metastasis by inducing a high nutrient demand in internal organs, indicating that reducing nutrient availability to tumor cells represents a potential strategy to prevent metastasis. See related commentary by Zerhouni and Piskounova, p. 4124.

Cross-species identification of cancer resistance-associated genes that may mediate human cancer risk.

Cancer is a predominant disease across animals. We applied a comparative genomics approach to systematically characterize genes whose conservation levels correlate positively (PC) or negatively (NC) with cancer resistance estimates across 193 vertebrates. Pathway analysis reveals that NC genes are enriched for metabolic functions and PC genes in cell cycle regulation, DNA repair, and immune response, pointing to their corresponding roles in mediating cancer risk. We find that PC genes are less tolerant to loss-of-function (LoF) mutations, are enriched in cancer driver genes, and are associated with germline mutations that increase human cancer risk. Their relevance to cancer risk is further supported via the analysis of mouse functional genomics and cancer mortality of zoo mammals' data. In sum, our study describes a cross-species genomic analysis pointing to candidate genes that may mediate human cancer risk.

Asymmetric inheritance of RNA toxicity in C. elegans expressing CTG repeats.

Nucleotide repeat expansions are a hallmark of over 40 neurodegenerative diseases and cause RNA toxicity and multisystemic symptoms that worsen with age. Through an unclear mechanism, RNA toxicity can trigger severe disease manifestation in infants if the repeats are inherited from their mother. Here we use Caenorhabditis elegans bearing expanded CUG repeats to show that this asymmetric intergenerational inheritance of toxicity contributes to disease pathogenesis. In addition, we show that this mechanism is dependent on small RNA pathways with maternal repeat-derived small RNAs causing transcriptomic changes in the offspring, reduced motility, and shortened lifespan. We rescued the toxicity phenotypes in the offspring by perturbing the RNAi machinery in the affected hermaphrodites. This points to a novel mechanism linking maternal bias and the RNAi machinery and suggests that toxic RNA is transmitted to offspring, causing disease phenotypes through intergenerational epigenetic inheritance.

DEPCOD: a tool to detect and visualize co-evolution of protein domains.

Proteins with similar phylogenetic patterns of conservation or loss across evolutionary taxa are strong candidates to work in the same cellular pathways or engage in physical or functional interactions. Our previously published tools implemented our method of normalized phylogenetic sequence profiling to detect functional associations between non-homologous proteins. However, many proteins consist of multiple protein domains subjected to different selective pressures, so using protein domain as the unit of analysis improves the detection of similar phylogenetic patterns. Here we analyze sequence conservation patterns across the whole tree of life for every protein domain from a set of widely studied organisms. The resulting new interactive webserver, DEPCOD (DEtection of Phylogenetically COrrelated Domains), performs searches with either a selected pre-defined protein domain or a user-supplied sequence as a query to detect other domains from the same organism that have similar conservation patterns. Top similarities on two evolutionary scales (the whole tree of life or eukaryotic genomes) are displayed along with known protein interactions and shared complexes, pathway enrichment among the hits, and detailed visualization of sources of detected similarities. DEPCOD reveals functional relationships between often non-homologous domains that could not be detected using whole-protein sequences. The web server is accessible at http://genetics.mgh.harvard.edu/DEPCOD.

Machine-learning of complex evolutionary signals improves classification of SNVs.

Conservation is a strong predictor for the pathogenicity of single-nucleotide variants (SNVs). However, some positions that present complex conservation patterns across vertebrates stray from this paradigm. Here, we analyzed the association between complex conservation patterns and the pathogenicity of SNVs in the 115 disease-genes that had sufficient variant data. We show that conservation is not a one-rule-fits-all solution since its accuracy highly depends on the analyzed set of species and genes. For example, pairwise comparisons between the human and 99 vertebrate species showed that species differ in their ability to predict the clinical outcomes of variants among different genes using conservation. Furthermore, certain genes were less amenable for conservation-based variant prediction, while others demonstrated species that optimize prediction. These insights led to developing EvoDiagnostics, which uses the conservation against each species as a feature within a random-forest machine-learning classification algorithm. EvoDiagnostics outperformed traditional conservation algorithms, deep-learning based methods and most ensemble tools in every prediction-task, highlighting the strength of optimizing conservation analysis per-species and per-gene. Overall, we suggest a new and a more biologically relevant approach for analyzing conservation, which improves prediction of variant pathogenicity.

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

Aurintricarboxylic Acid Decreases RNA Toxicity in a C. elegans Model of Repeat Expansions.

Pathologic expansions of DNA nucleotide tandem repeats may generate toxic RNA that triggers disease phenotypes. RNA toxicity is the hallmark of multiple expansion repeat disorders, including myotonic dystrophy type 1 (DM1). To date, there are no available disease-modifying therapies for DM1. Our aim was to use drug repositioning to ameliorate the phenotype of affected individuals in a nematode model of DM1. As the RNA interference pathway plays a key role in mediating RNA toxicity, we investigated the effect of aurintricarboxylic acid. We demonstrated that by perturbing the RNA interference machinery using aurintricarboxylic acid, we could annihilate the RNA toxicity and ameliorate the phenotype. As our approach targets a universal disease mechanism, it is potentially relevant for more expansion repeat disorders.

Co-evolution based machine-learning for predicting functional interactions between human genes.

Over the next decade, more than a million eukaryotic species are expected to be fully sequenced. This has the potential to improve our understanding of genotype and phenotype crosstalk, gene function and interactions, and answer evolutionary questions. Here, we develop a machine-learning approach for utilizing phylogenetic profiles across 1154 eukaryotic species. This method integrates co-evolution across eukaryotic clades to predict functional interactions between human genes and the context for these interactions. We benchmark our approach showing a 14% performance increase (auROC) compared to previous methods. Using this approach, we predict functional annotations for less studied genes. We focus on DNA repair and verify that 9 of the top 50 predicted genes have been identified elsewhere, with others previously prioritized by high-throughput screens. Overall, our approach enables better annotation of function and functional interactions and facilitates the understanding of evolutionary processes underlying co-evolution. The manuscript is accompanied by a webserver available at: https://mlpp.cs.huji.ac.il .

Expanding the MECP2 network using comparative genomics reveals potential therapeutic targets for Rett syndrome.

Inactivating mutations in the Methyl-CpG Binding Protein 2 (MECP2) gene are the main cause of Rett syndrome (RTT). Despite extensive research into MECP2 function, no treatments for RTT are currently available. Here, we used an evolutionary genomics approach to construct an unbiased MECP2 gene network, using 1028 eukaryotic genomes to prioritize proteins with strong co-evolutionary signatures with MECP2. Focusing on proteins targeted by FDA-approved drugs led to three promising targets, two of which were previously linked to MECP2 function (IRAK, KEAP1) and one that was not (EPOR). The drugs targeting these three proteins (Pacritinib, DMF, and EPO) were able to rescue different phenotypes of MECP2 inactivation in cultured human neural cell types, and appeared to converge on Nuclear Factor Kappa B (NF-κB) signaling in inflammation. This study highlights the potential of comparative genomics to accelerate drug discovery, and yields potential new avenues for the treatment of RTT.

Current gene panels account for nearly all homologous recombination repair-associated multiple-case breast cancer families.

It was hypothesized that variants in underexplored homologous recombination repair (HR) genes could explain unsolved multiple-case breast cancer (BC) families. We investigated HR deficiency (HRD)-associated mutational signatures and second hits in tumor DNA from familial BC cases. No candidates genes were associated with HRD in 38 probands previously tested negative with gene panels. We conclude it is unlikely that unknown HRD-associated genes explain a large fraction of unsolved familial BC.

PAX8 plays an essential antiapoptotic role in uterine serous papillary cancer.

Endometrial carcinoma (EC) is the fourth-most common cancer in women in the United States, and generally carries a favorable prognosis. However, about 10% of EC patients have a rare and aggressive form, uterine serous papillary carcinoma (USPC), which carries a much higher mortality rate. The developmental transcription factor PAX8 is expressed in nearly 100% of USPCs. We show that PAX8 plays a critical antiapoptotic role in USPC and this role is established via transcriptional activation of two aberrant signaling pathways. First, PAX8 positively regulates mutated p53, and missense p53 mutations have an oncogenic gain of function effect. Second, PAX8 directly transcriptionally regulates p21, in a p53-independent manner, and p21 acquires a growth promoting role that is mediated via cytoplasmic localization of the protein. We propose that mutated p53 and cytoplasmic p21 can independently mediate the pro-proliferative role of PAX8 in USPC. In addition, we performed a genome-wide transcriptome analysis to detect pathways that are regulated by PAX8, and propose that metabolism and HIF-1alpha -related pathways are potential candidates for mediating the role of PAX8 in USPC. Taken together our findings demonstrate for the first time that PAX8 is an essential lineage marker in USPC, and suggest its mechanism of action.

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.

Spatial patterning of scent in petunia corolla is discriminated by bees and involves the ABCG1 transporter.

Floral guides are patterned cues that direct the pollinator to the plant reproductive organs. The spatial distribution of showy visual and olfactory traits allows efficient plant-pollinator interactions. Data on the mechanisms underlying floral volatile patterns or their interactions with pollinators are lacking. Here we characterize the spatial emission patterns of volatiles from the corolla of the model plant Petunia × hybrida and reveal the ability of honeybees to distinguish these patterns. Along the adaxial epidermis, in correlation with cell density, the petal base adjacent to reproductive organs emitted significantly higher levels of volatiles than the distal petal rim. Volatile emission could also be differentiated between the two epidermal surfaces: emission from the adaxial side was significantly higher than that from the abaxial side. Similar emission patterns were also observed in other petunias, Dianthus caryophyllus (carnation) and Argyranthemum frutescens (Marguerite daisy). Analyses of transcripts involved in volatile production/emission revealed lower levels of the plasma-membrane transporter ABCG1 in the abaxial versus adaxial epidermis. Transient overexpression of ABCG1 enhanced emission from the abaxial epidermis to the level of the adaxial epidermis, suggesting its involvement in spatial emission patterns in the epidermal layers. Proboscis extension response experiments showed that differences in emission levels along the adaxial epidermis, that is, petal base versus rim, detected by GC-MS are also discernible by honeybees.

Alternative Splicing of the Inhibitory Immune Checkpoint Receptor SLAMF6 Generates a Dominant Positive Form, Boosting T-cell Effector Functions.

SLAMF6 is a homotypic receptor of the Ig-superfamily associated with progenitor-exhausted T cells. Here we show that in humans, SLAMF6 has three splice isoforms involving its V-domain. Although the canonical receptor inhibited T-cell activation through SAP recruitment, the short isoform SLAMF6Δ17-65 had a strong agonistic effect. The costimulatory action depended on protein phosphatase SHP1 and led to a cytotoxic molecular profile mediated by the expression of TBX21 and RUNX3. Patients treated with immune checkpoint blockade showed a shift toward SLAMF6Δ17-65 in peripheral blood T cells. We developed splice-switching antisense oligonucleotides (ASO) designed to target the relevant SLAMF6 splice junction. Our ASOs enhanced SLAMF6Δ17-65 expression in human tumor-infiltrating lymphocytes and improved their capacity to inhibit human melanoma in mice. The yin-yang relationship of SLAMF6 splice isoforms may represent a balancing mechanism that could be exploited to improve cancer immunotherapy.