Prediction of Metabolic Profiles from Transcriptomics Data in Human Cancer Cell Lines

The Metabolome and Transcriptome are mutually communicating within cancer cells, and this interplay is translated into the existence of quantifiable correlation structures between gene expression and metabolite abundance levels. Studying these correlations could provide a novel venue of understanding cancer and the discovery of novel biomarkers and pharmacological strategies, as well as laying the foundation for the prediction of metabolite quantities by leveraging information from the more widespread transcriptomics data. In the current paper, we investigate the correlation between gene expression and metabolite levels in the Cancer Cell Line Encyclopedia dataset, building a direct correlation network between the two molecular ensembles. We show that a metabolite/transcript correlation network can be used to predict metabolite levels in different samples and datasets, such as the NCI-60 cancer cell line dataset, both on a sample-by-sample basis and in differential contrasts. We also show that metabolite levels can be predicted in principle on any sample and dataset for which transcriptomics data are available, such as the Cancer Genome Atlas (TCGA)

The Transcriptome of SH-SY5Y at Single-Cell Resolution: A CITE-Seq Data Analysis Workflow

Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE-seq) is a recently established multimodal single cell analysis technique combining the immunophenotyping capabilities of antibody labeling and cell sorting with the resolution of single-cell RNA sequencing (scRNA-seq). By simply adding a 12-bp nucleotide barcode to antibodies (cell hashing), CITE-seq can be used to sequence antibody-bound tags alongside the cellular mRNA, thus reducing costs of scRNA-seq by performing it at the same time on multiple barcoded samples in a single run. Here, we illustrate an ideal CITE-seq data analysis workflow by characterizing the transcriptome of SH-SY5Y neuroblastoma cell line, a widely used model to study neuronal function and differentiation. We obtained transcriptomes from a total of 2879 single cells, measuring an average of 1600 genes/cell. Along with standard scRNA-seq data handling procedures, such as quality checks and cell filtering procedures, we performed exploratory analyses to identify most stable genes to be possibly used as reference housekeeping genes in qPCR experiments. We also illustrate how to use some popular R packages to investigate cell heterogeneity in scRNA-seq data, namely Seurat, Monocle, and slalom. Both the CITE-seq dataset and the code used to analyze it are freely shared and fully reusable for future research

Single-cell gene network analysis and transcriptional landscape of MYCN-amplified neuroblastoma cell lines

Neuroblastoma (NBL) is a pediatric cancer responsible for more than 15% of cancer deaths in children, with 800 new cases each year in the United States alone. Genomic amplification of the MYC oncogene family member MYCN characterizes a subset of high-risk pediatric neuroblastomas. Several cellular models have been implemented to study this disease over the years. Two of these, SK-N-BE-2-C (BE2C) and Kelly, are amongst the most used worldwide as models of MYCN-Amplified human NBL. Here, we provide a transcriptome-wide quantitative measurement of gene expression and transcriptional network activity in BE2C and Kelly cell lines at an unprecedented single-cell resolution. We obtained 1105 Kelly and 962 BE2C unsynchronized cells, with an average number of mapped reads/cell of roughly 38,000. The single-cell data recapitulate gene expression signatures previously generated from bulk RNA-Seq. We highlight low variance for commonly used housekeeping genes between different cells (ACTB, B2M and GAPDH), while showing higher than expected variance for metallothionein transcripts in Kelly cells. The high number of samples, despite the relatively low read coverage of single cells, allowed for robust pathway enrichment analysis and master regulator analysis (MRA), both of which highlight the more mesenchymal nature of BE2C cells as compared to Kelly cells, and the upregulation of TWIST1 and DNAJC1 transcriptional networks. We further defined master regulators at the single cell level and showed that MYCN is not constantly active or expressed within Kelly and BE2C cells, independently of cell cycle phase. The dataset, alongside a detailed and commented programming protocol to analyze it, is fully shared and reusable

A History of Repeated Alcohol Intoxication Promotes Cognitive Impairment and Gene Expression Signatures of Disease Progression in the 3xTg Mouse Model of Alzheimer’s Disease

The impact of alcohol abuse on Alzheimer’s disease (AD) is poorly understood. Here, we show that the onset of neurocognitive impairment in a mouse model of AD is hastened by repeated alcohol intoxication through exposure to alcohol vapor, and we provide a comprehensive gene expression dataset of the prefrontal cortex by the single-nucleus RNA sequencing of 113,242 cells. We observed a broad dysregulation of gene expression that involves neuronal excitability, neurodegeneration, and inflammation, including interferon genes. Several genes previously associated with AD in humans by genome-wide association studies were differentially regulated in specific neuronal populations. The gene expression signatures of AD mice with a history of alcohol intoxication showed greater similarity to the signatures of older AD mice with advanced disease and cognitive impairment than did the gene expression signatures of AD mice not exposed to alcohol, suggesting that alcohol promotes transcriptional changes consistent with AD progression. Our gene expression dataset at the single-cell level provides a unique re-source for investigations of the molecular bases of the detrimental role of excessive alcohol intake in AD

Gene Co-expression Analysis Identifies Histone Deacetylase 5 and 9 Expression in Midbrain Dopamine Neurons and as Regulators of Neurite Growth via Bone Morphogenetic Protein Signaling

Parkinson’s disease is characterized by the intracellular accumulation of α-synuclein which has been linked to early dopaminergic axonal degeneration. Identifying druggable targets that can promote axonal growth in cells overexpressing α-synuclein is important in order to develop strategies for early intervention. Class-IIa histone deacetylases (HDACs) have previously emerged as druggable targets, however, it is not known which specific class-IIa HDACs should be targeted to promote neurite growth in dopaminergic neurons. To provide insight into this, we used gene co-expression analysis to identify which, if any, of the class-IIa HDACs had a positive correlation with markers of dopaminergic neurons in the human substantia nigra. This revealed that two histone deacetylases, HDAC5 and HDAC9, are co-expressed with TH, GIRK2 and ALDH1A1 in the human SN. We further found that HDAC5 and HDAC9 are expressed in dopaminergic neurons in the adult mouse substantia nigra. We show that siRNAs targeting HDAC5 or HDAC9 can promote neurite growth in SH-SY5Y cells, and that their pharmacological inhibition, using the drug MC1568, promoted neurite growth in cultured rat dopaminergic neurons. Moreover, MC1568 treatment upregulated the expression of the neurotrophic factor, BMP2, and its downstream transcription factor, SMAD1. In addition, MC1568 or siRNAs targeting HDAC5 or HDAC9 led to an increase in Smad-dependent GFP expression in a reporter assay. Furthermore, MC1568 treatment of cultured rat dopaminergic neurons increased cellular levels of phosphorylated Smad1, which was prevented by the BMP receptor inhibitor, dorsomorphin. Dorsomorphin treatment prevented the neurite growth-promoting effects of siRNAs targeting HDAC5, as did overexpression of dominant-negative Smad4 or of the inhibitory Smad7, demonstrating a functional link to BMP signaling. Supplementation with BMP2 prevented the neurite growth-inhibitory effects of nuclear-restricted HDAC5. Finally, we report that siRNAs targeting HDAC5 or HDAC9 promoted neurite growth in cells overexpressing wild-type or A53T-α-synuclein and that MC1568 protected cultured rat dopaminergic neurons against the neurotoxin, MPP+. These findings establish HDAC5 and HDAC9 as novel regulators of BMP-Smad signaling, that additionally may be therapeutic targets worthy of further exploration in iPSC-derived human DA neurons and in vivo models of Parkinson’s disease

The inhibition of the highly expressed miR-221 and miR-222 impairs the growth of prostate carcinoma xenografts in mice

MiR-221 and miR-222 are two highly homologous microRNAs whose upregulation has been recently described in several types of human tumors, for some of which their oncogenic role was explained by the discovery of their target p27, a key cell cycle regulator. We previously showed this regulatory relationship in prostate carcinoma cell lines in vitro, underlying the role of miR-221/222 as inducers of proliferation and tumorigenicity

NF-kB and c-Jun induce the expression of the oncogenic miR-221 and miR-222 in prostate carcinoma and glioblastoma cells

MicroRNAs (miRNAs) are potent negative regulators of gene expression involved in all aspects of cell biology. They finely modulate virtually all physiological pathways in metazoans, and are deeply implicated in all main pathologies, among which cancer. Mir-221 and miR-222, two closely related miRNAs encoded in cluster from a genomic region on chromosome X, are strongly upregulated in several forms of human tumours. In this work, we report that the ectopic modulation of NF-kB modifies miR-221/222 expression in prostate carcinoma and glioblastoma cell lines, where we had previously shown their oncogenic activity. We identify two separate distal regions upstream of miR-221/222 promoter which are bound by the NF-kB subunit p65 and drive efficient transcription in luciferase reporter assays; consistently, the site-directed mutagenesis disrupting p65 binding sites or the ectopical inhibition of NF-kB activity significantly reduce luciferase activity. In the most distal enhancer region, we also define a binding site for c-Jun, and we show that the binding of this factor cooperates with that of p65, fully accounting for the observed upregulation of miR-221/222. Thus our work uncovers an additional mechanism through which NF-kB and c-Jun, two transcription factors deeply involved in cancer onset and progression, contribute to oncogenesis, by inducing miR-221/222 transcription

Sam68 splicing regulation contributes to motor unit establishment in the postnatal skeletal muscle

RNA-binding proteins orchestrate the composite life of RNA molecules and impact most physiological processes, thus underlying complex phenotypes. The RNA-binding protein Sam68 regulates differentiation processes by modulating splicing, polyadenylation, and stability of select transcripts. Herein, we found that Sam68-/- mice display altered regulation of alternative splicing in the spinal cord of key target genes involved in synaptic functions. Analysis of the motor units revealed that Sam68 ablation impairs the establishment of neuromuscular junctions and causes progressive loss of motor neurons in the spinal cord. Importantly, alterations of neuromuscular junction morphology and properties in Sam68-/- mice correlate with defects in muscle and motor unit integrity. Sam68-/- muscles display defects in postnatal development, with manifest signs of atrophy. Furthermore, fast-twitch muscles in Sam68-/- mice show structural features typical of slow-twitch muscles, suggesting alterations in the metabolic and functional properties of myofibers. Collectively, our data identify a key role for Sam68 in muscle development and suggest that proper establishment of motor units requires timely expression of synaptic splice variants

RGS4 negatively modulates Nociceptin/Orphanin FQ opioid receptor signaling: implication for L-Dopa induced dyskinesia

Background and purpose Regulator of G-protein signal 4 (RGS4) is a signal transduction protein that accelerates intrinsic GTPase activity of Gαi/o and Gαq subunits, suppressing GPCR signaling. Here we investigate whether RGS4 modulates nociceptin/orphanin FQ (N/OFQ) opioid (NOP) receptor signaling and this modulation has relevance for L-Dopa-induced dyskinesia. Experimental approach HEK293T cells transfected with NOP, NOP/RGS4 or NOP/RGS19 were challenged with N/OFQ and the small molecule NOP agonist AT-403, using D1-stimulated cAMP levels as a readout. Primary rat striatal neurons and adult mouse striatal slices were challenged with N/OFQ or AT-403 in the presence of the experimental RGS4 chemical probe, CCG-203920, and D1-stimulated cAMP or phosphorylated extracellular signal regulated kinase 1/2 (pERK) responses were monitored. In vivo, CCG-203920 was co-administered with AT-403 and L-Dopa to 6-hydroxydopamine hemilesioned rats, and dyskinetic movements, striatal biochemical correlates of dyskinesia (pERK and pGluR1 levels) and striatal RGS4 levels were measured. Key results RGS4 expression reduced NOFQ and AT-403 potency and efficacy in HEK293T cells. CCG-203920 increased N/OFQ potency in primary rat striatal neurons, and potentiated AT-403 response in mouse striatal slices. CCG-203920 enhanced AT-403 mediated inhibition of dyskinesia and its biochemical correlates, without compromising its motor-improving effects. Unilateral dopamine depletion caused bilateral reduction of RGS4 levels, which was reversed by L-Dopa. L-Dopa acutely upregulated RGS4 in the lesioned striatum. Conclusions and Implications RGS4 physiologically inhibits NOP receptor signaling. CCG-203920 enhanced NOP responses and improved the antidyskinetic potential of NOP receptor agonists, mitigating the effects of striatal RGS4 upregulation occurring during dyskinesia expression

Profiling cytotoxic microRNAs in pediatric and adult glioblastoma cells by high-content screening, identification, and validation of miR-1300

MicroRNAs play an important role in the regulation of mRNA translation and have therapeutic potential in cancer and other diseases. To profile the landscape of microRNAs with significant cytotoxicity in the context of glioblastoma (GBM), we performed a high-throughput screen in adult and pediatric GBM cells using a synthetic oligonucleotide library representing all known human microRNAs. Bioinformatics analysis was used to refine this list and the top seven microRNAs were validated in a larger panel of GBM cells using state-of-the-art in vitro assays. The cytotoxic effect of our most relevant candidate was assessed in a preclinical model. Our screen identified ~100 significantly cytotoxic microRNAs with 70% concordance between cell lines. MicroRNA-1300 (miR-1300) was the most potent and robust candidate. We observed a striking binucleated phenotype in miR-1300 transfected cells due to cytokinesis failure followed by apoptosis. This was also observed in two stem-like patient-derived cultures. We identified the physiological role of miR-1300 as a regulator of endomitosis in megakaryocyte differentiation where blockade of cytokinesis is an essential step. In GBM cells, where miR-1300 is normally not expressed, the oncogene Epithelial Cell Transforming 2 (ECT2) was validated as a direct key target. ECT2 siRNA phenocopied the effects of miR-1300, and ECT2 overexpression led to rescue of miR-1300 induced binucleation. We showed that ectopic expression of miR-1300 led to decreased tumor growth in an orthotopic GBM model. Our screen provides a resource for the neuro-oncology community and identified miR-1300 as a novel regulator of endomitosis with translatable potential for therapeutic application