Dissecting celastrol with machine learning to unveil dark pharmacology

Coallescing bespoke machine learning and bioinformatics analyses with cell-based assays we unveil pharmacology of celastrol. Celastrol is a direct modulator of the progesterone and cannabinoid receptors, and its effects correlate with the antiproliferative activity. We demonstrate how in silico methods may drive systems biology studies for natural products.We thank the Royal Society (URF\R\180019), H2020 (ERC StG, GA no. 676832, 743640 ad 807281) and FCT Portugal (IF/00624/2015 and 02/SAICT/2017, Grant 28333) for funding

Androgen-regulated transcription of ESRP2 drives alternative splicing patterns in prostate cancer

Prostate is the most frequent cancer in men. Prostate cancer progression is driven by androgen steroid hormones, and delayed by androgen deprivation therapy (ADT). Androgens control transcription by stimulating androgen receptor (AR) activity, yet also control pre-mRNA splicing through less clear mechanisms. Here we find androgens regulate splicing through AR-mediated transcriptional control of the epithelial-specific splicing regulator ESRP2. Both ESRP2 and its close paralog ESRP1 are highly expressed in primary prostate cancer. Androgen stimulation induces splicing switches in many endogenous ESRP2-controlled mRNA isoforms, including splicing switches correlating with disease progression. ESRP2 expression in clinical prostate cancer is repressed by ADT, which may thus inadvertently dampen epithelial splice programmes. Supporting this, treatment with the AR antagonist bicalutamide (Casodex) induced mesenchymal splicing patterns of genes including FLNB and CTNND1. Our data reveals a new mechanism of splicing control in prostate cancer with important implications for disease progression.This article is freely available via Open Access. Click on the Publisher URL to access the full-text via the publisher's site

Widespread sex differences in gene expression and splicing in the adult human brain

There is strong evidence to show that men and women differ in terms of neurodevelopment, neurochemistry and susceptibility to neurodegenerative and neuropsychiatric disease. The molecular basis of these differences remains unclear. Progress in this field has been hampered by the lack of genome-wide information on sex differences in gene expression and in particular splicing in the human brain. Here we address this issue by using post-mortem adult human brain and spinal cord samples originating from 137 neuropathologically confirmed control individuals to study whole-genome gene expression and splicing in 12 CNS regions. We show that sex differences in gene expression and splicing are widespread in adult human brain, being detectable in all major brain regions and involving 2.5% of all expressed genes. We give examples of genes where sex-biased expression is both disease-relevant and likely to have functional consequences, and provide evidence suggesting that sex biases in expression may reflect sex-biased gene regulatory structures

Cis and Trans Effects of Human Genomic Variants on Gene Expression

This work was funded by the Louis-Jeantet Foundation (http://www.jeantet.ch/), the European Research Council (Grant ID: 260927 http://erc.europa.eu/), the Swiss National Foundation (Grant ID: 130342 http://www.snf.ch), NCCR Frontiers In Genetics (http://www.frontiers-in-genetics.org), the UK Medical Research Council (http://www.mrc.ac.uk) and the Wellcome Trust (Grant ID: 092731).

The Choice of the Filtering Method in Microarrays Affects the Inference Regarding Dosage Compensation of the Active X-Chromosome

The hypothesis of dosage compensation of genes of the X chromosome, supported by previous microarray studies, was recently challenged by RNA-sequencing data. It was suggested that microarray studies were biased toward an over-estimation of X-linked expression levels as a consequence of the filtering of genes below the detection threshold of microarrays.To investigate this hypothesis, we used microarray expression data from circulating monocytes in 1,467 individuals. In total, 25,349 and 1,156 probes were unambiguously assigned to autosomes and the X chromosome, respectively. Globally, there was a clear shift of X-linked expressions toward lower levels than autosomes. We compared the ratio of expression levels of X-linked to autosomal transcripts (X∶AA) using two different filtering methods: 1. gene expressions were filtered out using a detection threshold irrespective of gene chromosomal location (the standard method in microarrays); 2. equal proportions of genes were filtered out separately on the X and on autosomes. For a wide range of filtering proportions, the X∶AA ratio estimated with the first method was not significantly different from 1, the value expected if dosage compensation was achieved, whereas it was significantly lower than 1 with the second method, leading to the rejection of the hypothesis of dosage compensation. We further showed in simulated data that the choice of the most appropriate method was dependent on biological assumptions regarding the proportion of actively expressed genes on the X chromosome comparative to the autosomes and the extent of dosage compensation.This study shows that the method used for filtering out lowly expressed genes in microarrays may have a major impact according to the hypothesis investigated. The hypothesis of dosage compensation of X-linked genes cannot be firmly accepted or rejected using microarray-based data

BeadArray Expression Analysis Using Bioconductor

Illumina whole-genome expression BeadArrays are a popular choice in gene profiling studies. Aside from the vendor-provided software tools for analyzing BeadArray expression data (GenomeStudio/BeadStudio), there exists a comprehensive set of open-source analysis tools in the Bioconductor project, many of which have been tailored to exploit the unique properties of this platform. In this article, we explore a number of these software packages and demonstrate how to perform a complete analysis of BeadArray data in various formats. The key steps of importing data, performing quality assessments, preprocessing, and annotation in the common setting of assessing differential expression in designed experiments will be covered

The fitness cost of mis-splicing is the main determinant of alternative splicing patterns

Background Most eukaryotic genes are subject to alternative splicing (AS), which may contribute to the production of protein variants or to the regulation of gene expression via nonsense-mediated messenger RNA (mRNA) decay (NMD). However, a fraction of splice variants might correspond to spurious transcripts and the question of the relative proportion of splicing errors to functional splice variants remains highly debated. Results We propose a test to quantify the fraction of AS events corresponding to errors. This test is based on the fact that the fitness cost of splicing errors increases with the number of introns in a gene and with expression level. We analyzed the transcriptome of the intron-rich eukaryote Paramecium tetraurelia. We show that in both normal and in NMD-deficient cells, AS rates strongly decrease with increasing expression level and with increasing number of introns. This relationship is observed for AS events that are detectable by NMD as well as for those that are not, which invalidates the hypothesis of a link with the regulation of gene expression. Our results show that in genes with a median expression level, 92–98% of observed splice variants correspond to errors. We observed the same patterns in human transcriptomes and we further show that AS rates correlate with the fitness cost of splicing errors. Conclusions These observations indicate that genes under weaker selective pressure accumulate more maladaptive substitutions and are more prone to splicing errors. Thus, to a large extent, patterns of gene expression variants simply reflect the balance between selection, mutation, and drift

Direct integration of intensity-level data from Affymetrix and Illumina microarrays improves statistical power for robust reanalysis

<p>Abstract</p> <p>Background</p> <p>Affymetrix GeneChips and Illumina BeadArrays are the most widely used commercial single channel gene expression microarrays. Public data repositories are an extremely valuable resource, providing array-derived gene expression measurements from many thousands of experiments. Unfortunately many of these studies are underpowered and it is desirable to improve power by combining data from more than one study; we sought to determine whether platform-specific bias precludes direct integration of probe intensity signals for combined reanalysis.</p> <p>Results</p> <p>Using Affymetrix and Illumina data from the microarray quality control project, from our own clinical samples, and from additional publicly available datasets we evaluated several approaches to directly integrate intensity level expression data from the two platforms. After mapping probe sequences to Ensembl genes we demonstrate that, ComBat and cross platform normalisation (XPN), significantly outperform mean-centering and distance-weighted discrimination (DWD) in terms of minimising inter-platform variance. In particular we observed that DWD, a popular method used in a number of previous studies, removed systematic bias at the expense of genuine biological variability, potentially reducing legitimate biological differences from integrated datasets.</p> <p>Conclusion</p> <p>Normalised and batch-corrected intensity-level data from Affymetrix and Illumina microarrays can be directly combined to generate biologically meaningful results with improved statistical power for robust, integrated reanalysis.</p

An shRNA-Based Screen of Splicing Regulators Identifies SFRS3 as a Negative Regulator of IL-1β Secretion

The generation of diversity and plasticity of transcriptional programs are key components of effective vertebrate immune responses. The role of Alternative Splicing has been recognized, but it is underappreciated and poorly understood as a critical mechanism for the regulation and fine-tuning of physiological immune responses. Here we report the generation of loss-of-function phenotypes for a large collection of genes known or predicted to be involved in the splicing reaction and the identification of 19 novel regulators of IL-1β secretion in response to E. coli challenge of THP-1 cells. Twelve of these genes are required for IL-1β secretion, while seven are negative regulators of this process. Silencing of SFRS3 increased IL-1β secretion due to elevation of IL-1β and caspase-1 mRNA in addition to active caspase-1 levels. This study points to the relevance of splicing in the regulation of auto-inflammatory diseases