Revealing disease-associated pathways by network integration of untargeted metabolomics.

Uncovering the molecular context of dysregulated metabolites is crucial to understand pathogenic pathways. However, their system-level analysis has been limited owing to challenges in global metabolite identification. Most metabolite features detected by untargeted metabolomics carried out by liquid-chromatography-mass spectrometry cannot be uniquely identified without additional, time-consuming experiments. We report a network-based approach, prize-collecting Steiner forest algorithm for integrative analysis of untargeted metabolomics (PIUMet), that infers molecular pathways and components via integrative analysis of metabolite features, without requiring their identification. We demonstrated PIUMet by analyzing changes in metabolism of sphingolipids, fatty acids and steroids in a Huntington's disease model. Additionally, PIUMet enabled us to elucidate putative identities of altered metabolite features in diseased cells, and infer experimentally undetected, disease-associated metabolites and dysregulated proteins. Finally, we established PIUMet's ability for integrative analysis of untargeted metabolomics data with proteomics data, demonstrating that this approach elicits disease-associated metabolites and proteins that cannot be inferred by individual analysis of these data.

Multiplexed barcoded CRISPR-Cas9 screening enabled by CombiGEM.

The orchestrated action of genes controls complex biological phenotypes, yet the systematic discovery of gene and drug combinations that modulate these phenotypes in human cells is labor intensive and challenging to scale. Here, we created a platform for the massively parallel screening of barcoded combinatorial gene perturbations in human cells and translated these hits into effective drug combinations. This technology leverages the simplicity of the CRISPR-Cas9 system for multiplexed targeting of specific genomic loci and the versatility of combinatorial genetics en masse (CombiGEM) to rapidly assemble barcoded combinatorial genetic libraries that can be tracked with high-throughput sequencing. We applied CombiGEM-CRISPR to create a library of 23,409 barcoded dual guide-RNA (gRNA) combinations and then perform a high-throughput pooled screen to identify gene pairs that inhibited ovarian cancer cell growth when they were targeted. We validated the growth-inhibiting effects of specific gene sets, including epigenetic regulators KDM4C/BRD4 and KDM6B/BRD4, via individual assays with CRISPR-Cas-based knockouts and RNA-interference-based knockdowns. We also tested small-molecule drug pairs directed against our pairwise hits and showed that they exerted synergistic antiproliferative effects against ovarian cancer cells. We envision that the CombiGEM-CRISPR platform will be applicable to a broad range of biological settings and will accelerate the systematic identification of genetic combinations and their translation into novel drug combinations that modulate complex human disease phenotypes.

Hydrogen peroxide-mediated induction of SOD1 gene transcription is independent from Nrf2 in a cellular model of neurodegeneration.

BACKGROUND: It is still unclear whether oxidative stress (OS) is a disease consequence or is directly involved in the etiology of neurodegenerative disorders (NDs) onset and/or progression; however, many of these conditions are associated with increased levels of oxidation markers and damaged cell components. Previously we demonstrated the accumulation of reactive oxygen species (ROS) and increased SOD1 gene expression in H2O2-treated SH-SY5Y cells, recapitulating pathological features of Amyotrophic Lateral Sclerosis (ALS). Since we observed a post-transcriptional regulation of SOD1 gene in this cellular model, we investigated the transcriptional regulation of SOD1 mRNA under oxidative stress (OS). RESULTS: In response to H2O2 treatment, PolII increased its association to SOD1 promoter. Electrophoretic mobility shift assays and mass spectrometry analyses on SOD1 promoter highlighted the formation of a transcriptional complex bound to the ARE sequences. Western Blotting experiments showed that in our in vitro model, H2O2 exposure increases Nrf2 expression in the nuclear fraction while immunoprecipitation confirmed its phosphorylation and release from Keap1 inhibition. However, H2O2 treatment did not modify Nrf2 binding on SOD1 promoter, which seems to be regulated by different transcription factors (TFs). CONCLUSIONS: Although our data suggest that SOD1 is transcriptionally regulated in response to OS, Nrf2 does not appear to associate with SOD1 promoter in this cellular model of neurodegeneration. Our results open new perspectives in the comprehension of two key antioxidant pathways involved in neurodegenerative disorders.

Posttranscriptional regulation of SOD1 gene expression under oxidative stress: Potential role of ELAV proteins in sporadic ALS.

Increased levels of SOD1 mRNA have been observed in sporadic ALS patients (SALS) compared to controls. Hence, the understanding of the mechanisms by which SOD1 gene expression is modulated may shed new light on SOD1 involvement in ALS. Of interest, some adenine/uracil-rich elements (AREs) in SOD1 3'-untranslated region have been identified. These sequences represent the docking sites for several RNA-binding proteins such as ELAV proteins (ELAVs), positive regulators of gene expression. We first investigated in SH-SY5Y cells whether SOD1 mRNA represents a target of ELAVs. Results from RNA Electrophoretic Mobility Shift and RNA-immunoprecipitation assays showed a molecular interaction between ELAVs and SOD1 mRNA. We also observed that the treatment with H2O2 induced a significant increase of the amount of SOD1 mRNA bound by ELAVs and an up-regulation of SOD1 protein levels. We found a specific increase in ELAV/HuR phosphorylation, suggesting activation of this protein, in peripheral blood mononuclear cells from SALS patients compared to controls. Finally, we found increased levels of ELAV proteins in the motor cortex and spinal cord from SALS patients compared to controls, in parallel with SOD1 up-regulation in the same areas. This study suggests, for the first time, that ELAVs are involved in the regulation of SOD1 gene expression at post-transcriptional level and that these proteins are more activated in ALS pathology. The link between ELAVs and SOD1 may open novel perspectives for ALS research, paving the way for new therapeutic options.

Ubiquilin 2 mutations in Italian patients with amyotrophic lateral sclerosis and frontotemporal dementia.

OBJECTIVES: Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease mainly involving cortical and spinal motor neurones. Molecular studies have recently identified different mutations in the  ubiquilin-2 (UBQLN2) gene as causative of a familial form of X-linked ALS, 90% penetrant in women. The aim of our study was to analyse the UBQLN2 gene in a large cohort of patients with familial (FALS) and sporadic (SALS) amyotrophic lateral sclerosis, with or without frontotemporal dementia (FTD), and in patients with FTD. METHODS: We analysed the UBQLN2 gene in 819 SALS cases, 226 FALS cases, 53 ALS-FTD patients, and 63 patients with a clinical record of FTD. Molecular analysis of the entire coding sequence was carried out in all FALS and ALS-FTD patients, while SALS and FTD patients were analysed specifically for the genomic region coding for the PXX repeat tract. Healthy controls were 845 anonymous blood donors and were screened for the PXX repeat region only. RESULTS: We found five different variants in the UBQLN2 gene in five unrelated ALS patients. Three variants, including two novel ones, involved a proline residue in the PXX repeat region and were found in three FALS cases. The other two were novel variants, identified in one FALS and one SALS patient. None of these variants was present in controls, while one control carried a new heterozygous variant. CONCLUSIONS: Our data support the role of the UBQLN2 gene in the pathogenesis of FALS, being conversely a rare genetic cause in SALS even when complicated by FTD.

Prediction of Immunotherapy Response in Melanoma through Combined Modeling of Neoantigen Burden and Immune-Related Resistance Mechanisms.

PURPOSE: While immune checkpoint blockade (ICB) has become a pillar of cancer treatment, biomarkers that consistently predict patient response remain elusive due to the complex mechanisms driving immune response to tumors. We hypothesized that a multi-dimensional approach modeling both tumor and immune-related molecular mechanisms would better predict ICB response than simpler mutation-focused biomarkers, such as tumor mutational burden (TMB). EXPERIMENTAL DESIGN: Tumors from a cohort of patients with late-stage melanoma (n = 51) were profiled using an immune-enhanced exome and transcriptome platform. We demonstrate increasing predictive power with deeper modeling of neoantigens and immune-related resistance mechanisms to ICB. RESULTS: Our neoantigen burden score, which integrates both exome and transcriptome features, more significantly stratified responders and nonresponders (P = 0.016) than TMB alone (P = 0.049). Extension of this model to include immune-related resistance mechanisms affecting the antigen presentation machinery, such as HLA allele-specific LOH, resulted in a composite neoantigen presentation score (NEOPS) that demonstrated further increased association with therapy response (P = 0.002). CONCLUSIONS: NEOPS proved the statistically strongest biomarker compared with all single-gene biomarkers, expression signatures, and TMB biomarkers evaluated in this cohort. Subsequent confirmation of these findings in an independent cohort of patients (n = 110) suggests that NEOPS is a robust, novel biomarker of ICB response in melanoma.

An integrated multi-omic analysis of iPSC-derived motor neurons from C9ORF72 ALS patients.

Neurodegenerative diseases are challenging for systems biology because of the lack of reliable animal models or patient samples at early disease stages. Induced pluripotent stem cells (iPSCs) could address these challenges. We investigated DNA, RNA, epigenetics, and proteins in iPSC-derived motor neurons from patients with ALS carrying hexanucleotide expansions in C9ORF72. Using integrative computational methods combining all omics datasets, we identified novel and known dysregulated pathways. We used a C9ORF72 Drosophila model to distinguish pathways contributing to disease phenotypes from compensatory ones and confirmed alterations in some pathways in postmortem spinal cord tissue of patients with ALS. A different differentiation protocol was used to derive a separate set of C9ORF72 and control motor neurons. Many individual -omics differed by protocol, but some core dysregulated pathways were consistent. This strategy of analyzing patient-specific neurons provides disease-related outcomes with small numbers of heterogeneous lines and reduces variation from single-omics to elucidate network-based signatures.

A Multi-Omics Interpretable Machine Learning Model Reveals Modes of Action of Small Molecules.

High-throughput screening and gene signature analyses frequently identify lead therapeutic compounds with unknown modes of action (MoAs), and the resulting uncertainties can lead to the failure of clinical trials. We developed an approach for uncovering MoAs through an interpretable machine learning model of transcriptomics, epigenomics, metabolomics, and proteomics. Examining compounds with beneficial effects in models of Huntington's Disease, we found common MoAs for compounds with unrelated structures, connectivity scores, and binding targets. The approach also predicted highly divergent MoAs for two FDA-approved antihistamines. We experimentally validated these effects, demonstrating that one antihistamine activates autophagy, while the other targets bioenergetics. The use of multiple omics was essential, as some MoAs were virtually undetectable in specific assays. Our approach does not require reference compounds or large databases of experimental data in related systems and thus can be applied to the study of agents with uncharacterized MoAs and to rare or understudied diseases.

The Library of Integrated Network-Based Cellular Signatures NIH Program: System-Level Cataloging of Human Cells Response to Perturbations.

The Library of Integrated Network-Based Cellular Signatures (LINCS) is an NIH Common Fund program that catalogs how human cells globally respond to chemical, genetic, and disease perturbations. Resources generated by LINCS include experimental and computational methods, visualization tools, molecular and imaging data, and signatures. By assembling an integrated picture of the range of responses of human cells exposed to many perturbations, the LINCS program aims to better understand human disease and to advance the development of new therapies. Perturbations under study include drugs, genetic perturbations, tissue micro-environments, antibodies, and disease-causing mutations. Responses to perturbations are measured by transcript profiling, mass spectrometry, cell imaging, and biochemical methods, among other assays. The LINCS program focuses on cellular physiology shared among tissues and cell types relevant to an array of diseases, including cancer, heart disease, and neurodegenerative disorders. This Perspective describes LINCS technologies, datasets, tools, and approaches to data accessibility and reusability.

Identifying therapeutic targets by combining transcriptional data with ordinal clinical measurements.

The immense and growing repositories of transcriptional data may contain critical insights for developing new therapies. Current approaches to mining these data largely rely on binary classifications of disease vs. control, and are not able to incorporate measures of disease severity. We report an analytical approach to integrate ordinal clinical information with transcriptomics. We apply this method to public data for a large cohort of Huntington's disease patients and controls, identifying and prioritizing phenotype-associated genes. We verify the role of a high-ranked gene in dysregulation of sphingolipid metabolism in the disease and demonstrate that inhibiting the enzyme, sphingosine-1-phosphate lyase 1 (SPL), has neuroprotective effects in Huntington's disease models. Finally, we show that one consequence of inhibiting SPL is intracellular inhibition of histone deacetylases, thus linking our observations in sphingolipid metabolism to a well-characterized Huntington's disease pathway. Our approach is easily applied to any data with ordinal clinical measurements, and may deepen our understanding of disease processes.Identifying gene subsets affecting disease phenotypes from transcriptome data is challenge. Here, the authors develop a method that combines transcriptional data with disease ordinal clinical measurements to discover a sphingolipid metabolism regulator involving in Huntington's disease progression.

Cell freezing protocol suitable for ATAC-Seq on motor neurons derived from human induced pluripotent stem cells.

In recent years, the assay for transposase-accessible chromatin using sequencing (ATAC-Seq) has become a fundamental tool of epigenomic research. However, it is difficult to perform this technique on frozen samples because freezing cells before extracting nuclei can impair nuclear integrity and alter chromatin structure, especially in fragile cells such as neurons. Our aim was to develop a protocol for freezing neuronal cells that is compatible with ATAC-Seq; we focused on a disease-relevant cell type, namely motor neurons differentiated from induced pluripotent stem cells (iMNs) from a patient affected by spinal muscular atrophy. We found that while flash-frozen iMNs are not suitable for ATAC-Seq, the assay is successful with slow-cooled cryopreserved cells. Using this method, we were able to isolate high quality, intact nuclei, and we verified that epigenetic results from fresh and cryopreserved iMNs quantitatively agree.

Elucidating MicroRNA Regulatory Networks Using Transcriptional, Post-transcriptional, and Histone Modification Measurements.

MicroRNAs (miRNAs) regulate diverse biological processes by repressing mRNAs, but their modest effects on direct targets, together with their participation in larger regulatory networks, make it challenging to delineate miRNA-mediated effects. Here, we describe an approach to characterizing miRNA-regulatory networks by systematically profiling transcriptional, post-transcriptional and epigenetic activity in a pair of isogenic murine fibroblast cell lines with and without Dicer expression. By RNA sequencing (RNA-seq) and CLIP (crosslinking followed by immunoprecipitation) sequencing (CLIP-seq), we found that most of the changes induced by global miRNA loss occur at the level of transcription. We then introduced a network modeling approach that integrated these data with epigenetic data to identify specific miRNA-regulated transcription factors that explain the impact of miRNA perturbation on gene expression. In total, we demonstrate that combining multiple genome-wide datasets spanning diverse regulatory modes enables accurate delineation of the downstream miRNA-regulated transcriptional network and establishes a model for studying similar networks in other systems.

Network Modeling Identifies Patient-specific Pathways in Glioblastoma.

Glioblastoma is the most aggressive type of malignant human brain tumor. Molecular profiling experiments have revealed that these tumors are extremely heterogeneous. This heterogeneity is one of the principal challenges for developing targeted therapies. We hypothesize that despite the diverse molecular profiles, it might still be possible to identify common signaling changes that could be targeted in some or all tumors. Using a network modeling approach, we reconstruct the altered signaling pathways from tumor-specific phosphoproteomic data and known protein-protein interactions. We then develop a network-based strategy for identifying tumor specific proteins and pathways that were predicted by the models but not directly observed in the experiments. Among these hidden targets, we show that the ERK activator kinase1 (MEK1) displays increased phosphorylation in all tumors. By contrast, protein numb homolog (NUMB) is present only in the subset of the tumors that are the most invasive. Additionally, increased S100A4 is associated with only one of the tumors. Overall, our results demonstrate that despite the heterogeneity of the proteomic data, network models can identify common or tumor specific pathway-level changes. These results represent an important proof of principle that can improve the target selection process for tumor specific treatments.

SOD1 and DJ-1 converge at Nrf2 pathway: a clue for antioxidant therapeutic potential in neurodegeneration.

Neurodegenerative diseases share diverse pathological features and among these oxidative stress (OS) plays a leading role. Impaired activity and reduced expression of antioxidant proteins have been reported as common events in several aging-associated disorders. In this review paper, we first provide an overview of the involvement of reactive oxygen species- (ROS-) induced oxidative damage in Parkinson's disease (PD) and amyotrophic lateral sclerosis (ALS). Subsequently, we focus on DJ-1 and SOD1 proteins, which are involved in PD and ALS and also exert a prominent role in the interaction between redox homeostasis and neurodegeneration. Interestingly, recent studies demonstrated that DJ-1 and SOD1 are both tightly connected with Nrf2 protein, a transcriptional factor and master regulator of the expression of many antioxidant/detoxification genes. Nrf2 is emerging as a key neuroprotective protein in neurodegenerative diseases, since it helps neuronal cells to cope with toxic insults and OS. We herein summarize the recent literature providing a detailed picture of the promising therapeutic efficacy of Nrf2 natural and synthetic inducers as disease-modifying molecules for the treatment of neurodegenerative diseases.

Regulation of FMO and PON detoxication systems in ALS human tissues.

Amyotrophic lateral sclerosis (ALS) is an adult-onset, progressive, and fatal neurodegenerative disease with unknown etiology. Recent evidence suggests an association between the exposure to toxic environmental factors and sporadic ALS. The flavin-containing monooxygenases (FMOs) and paraoxonase (PONs) genes encode enzymes involved in xenobiotic detoxication and are associated with ALS. FMO and PON gene expression has been examined in the human central nervous system including human brain subregions defined as the spinal cord, medulla, and cerebral cortex and in the peripheral tissues (lymphocytes, fibroblasts) in ALS patients and normal control subjects. FMO expression was generally higher in tissues from ALS subjects than in control tissues, with the largest increases in FMO expression detected in the spinal cord. In peripheral tissues, the FMO mRNA level was found to be lower compared with FMO expression in brain tissue, and no differences were detected between ALS patients and the control tissue. FMO and PON gene expression was low in peripheral tissues. In contrast to FMO5 expression, the PON2 gene was down-regulated in ALS patients compared to the controls. Because FMO and PON are involved in the detoxication processes and their functional activity to bioactivate chemicals to toxins has been documented, the data herein suggest that environmental toxin exposure may play a role in a subset of individuals who contract ALS by altering FMO and PON gene expression. Although the precise pathogenic link is presently unknown, these findings suggest a role at FMO and PON genes in the development of ALS.

Altered intracellular localization of SOD1 in leukocytes from patients with sporadic amyotrophic lateral sclerosis.

Several lines of evidence support the hypothesis of a toxic role played by wild type SOD1 (WT-SOD1) in the pathogenesis of sporadic amyotrophic lateral sclerosis (SALS). In this study we investigated both distribution and expression profile of WT-SOD1 in leukocytes from 19 SALS patients and 17 healthy individuals. Immunofluorescence experiments by confocal microscopy showed that SOD1 accumulates in the nuclear compartment in a group of SALS subjects. These results were also confirmed by western blot carried out on soluble nuclear and cytoplasmic fractions, with increased nuclear SOD1 level (p<0.05). In addition, we observed the presence of cytoplasmic SOD1 aggregates in agreement with an increased amount of the protein recovered by the insoluble fraction. A further confirmation of the overall increased level of SOD1 has been obtained from single cells analysis using flow cytometry as cells from SALS patients showed an higher SOD1 protein content (p<0.05). These findings add further evidence to the hypothesis of an altered WT-SOD1 expression profile in peripheral blood mononuclear cells (PBMCs) from patients with ALS suggesting that WT-SOD1 species with different degrees of solubility could be involved in the pathogenesis of the disease.

Analysis of hnRNPA1, A2/B1, and A3 genes in patients with amyotrophic lateral sclerosis.

Mutations in the prion-like domain (PrLD) of hnRNPA1 and A2/B1 genes were recently identified in 2 families with inclusion body myopathy associated with Paget disease of bone, frontotemporal dementia (FTD), and amyotrophic lateral sclerosis, and in ALS patients. These mutations were shown to increase the propensity of hnRNPA1 and A2/B1 proteins, which are TDP-43-binding partners, to self-aggregate. hnRNPA3 protein contains a similar PrLD and was recently described in the p62-positive/TDP-43-negative inclusions in affected tissues of C9orf72-mutated ALS/FTD patients. We screened hnRNPA1, A2/B1, and A3 genes in a cohort of 113 familial ALS (FALS) individuals without mutations in other known ALS-causative genes. We extended our analysis to 108 FALS with mutations in other ALS-associated genes and to 622 sporadic cases by screening specifically the PrLDs of hnRNPA1, A2/B1, and A3. We failed to find variants in each cohort. Our results suggest that mutations in hnRNPA1, A2/B1, and A3 genes are a rare finding in ALS.

From Transcriptome to Noncoding RNAs: Implications in ALS Mechanism.

In the last years, numerous studies have focused on understanding the metabolism of RNA and its implication in disease processes but abnormal RNA metabolism is still unknown. RNA plays a central role in translating genetic information into proteins and in many other catalytic and regulatory tasks. Recent advances in the study of RNA metabolism revealed complex pathways for the generation and maintenance of functional RNA in amyotrophic lateral sclerosis (ALS). Interestingly, perturbations in RNA processing have been described in ALS at various levels such as gene transcription, mRNA stabilization, transport, and translational regulations. In this paper, we will discuss the alteration of RNA profile in ALS disease, starting from transcription, the first step leading to gene expression, through the posttranscriptional regulation, including RNA/DNA binding proteins and aberrant exon splicing to protein noncoding RNAs, as lncRNA and microRNA.

SOD1 Transcriptional and Posttranscriptional Regulation and Its Potential Implications in ALS.

Copper-zinc superoxide dismutase (SOD1) is a detoxifying enzyme localized in the cytosol, nucleus, peroxisomes, and mitochondria. The discovery that mutations in SOD1 gene cause a subset of familial amyotrophic lateral sclerosis (FALS) has attracted great attention, and studies to date have been mainly focused on discovering mutations in the coding region and investigation at protein level. Considering that changes in SOD1 mRNA levels have been associated with sporadic ALS (SALS), a molecular understanding of the processes involved in the regulation of SOD1 gene expression could not only unravel novel regulatory pathways that may govern cellular phenotypes and changes in diseases but also might reveal therapeutic targets and treatments. This review seeks to provide an overview of SOD1 gene structure and of the processes through which SOD1 transcription is controlled. Furthermore, we emphasize the importance to focus future researches on investigating posttranscriptional mechanisms and their relevance to ALS.

Lack of association of PON polymorphisms with sporadic ALS in an Italian population.

Paraoxonase (PON) gene polymorphisms have been associated with susceptibility to sporadic amyotrophic lateral sclerosis (ALS). We have investigated the role of the previously associated single nucleotide polymorphisms rs854560, rs662, and rs6954345 in 350 ALS patients and 376 matched controls from Italy. No significant association was observed at genotype and haplotype level. Our data suggest that PON polymorphisms are not involved in ALS pathogenesis in an Italian population.