Lukas M. Simon

The spliceosome is a therapeutic vulnerability in MYC-driven cancer

MYC (also known as c-MYC) overexpression or hyperactivation is one of the most common drivers of human cancer. Despite intensive study, the MYC oncogene remains recalcitrant to therapeutic inhibition. MYC is a transcription factor, and many of its pro-tumorigenic functions have been attributed to its ability to regulate gene expression programs. Notably, oncogenic MYC activation has also been shown to increase total RNA and protein production in many tissue and disease contexts. While such increases in RNA and protein production may endow cancer cells with pro-tumour hallmarks, this increase in synthesis may also generate new or heightened burden on MYC-driven cancer cells to process these macromolecules properly. Here we discover that the spliceosome is a new target of oncogenic stress in MYC-driven cancers. We identify BUD31 as a MYC-synthetic lethal gene in human mammary epithelial cells, and demonstrate that BUD31 is a component of the core spliceosome required for its assembly and catalytic activity. Core spliceosomal factors (such as SF3B1 and U2AF1) associated with BUD31 are also required to tolerate oncogenic MYC. Notably, MYC hyperactivation induces an increase in total precursor messenger RNA synthesis, suggesting an increased burden on the core spliceosome to process pre-mRNA. In contrast to normal cells, partial inhibition of the spliceosome in MYC-hyperactivated cells leads to global intron retention, widespread defects in pre-mRNA maturation, and deregulation of many essential cell processes. Notably, genetic or pharmacological inhibition of the spliceosome in vivo impairs survival, tumorigenicity and metastatic proclivity of MYC-dependent breast cancers. Collectively, these data suggest that oncogenic MYC confers a collateral stress on splicing, and that components of the spliceosome may be therapeutic entry points for aggressive MYC-driven cancers.

Racial differences in human platelet PAR4 reactivity reflect expression of PCTP and miR-376c

Racial differences in the pathophysiology of atherothrombosis are poorly understood. We explored the function and transcriptome of platelets in healthy black (n = 70) and white (n = 84) subjects. Platelet aggregation and calcium mobilization induced by the PAR4 thrombin receptor were significantly greater in black subjects. Numerous differentially expressed RNAs were associated with both race and PAR4 reactivity, including PCTP (encoding phosphatidylcholine transfer protein), and platelets from black subjects expressed higher levels of PC-TP protein. PC-TP inhibition or depletion blocked PAR4- but not PAR1-mediated activation of platelets and megakaryocytic cell lines. miR-376c levels were differentially expressed by race and PAR4 reactivity and were inversely correlated with PCTP mRNA levels, PC-TP protein levels and PAR4 reactivity. miR-376c regulated the expression of PC-TP in human megakaryocytes. A disproportionately high number of microRNAs that were differentially expressed by race and PAR4 reactivity, including miR-376c, are encoded in the DLK1-DIO3 locus and were expressed at lower levels in platelets from black subjects. These results suggest that PC-TP contributes to the racial difference in PAR4-mediated platelet activation, indicate a genomic contribution to platelet function that differs by race and emphasize a need to consider the effects of race when developing anti-thrombotic drugs.

Human platelet microRNA-mRNA networks associated with age and gender revealed by integrated plateletomics.

There is little data considering relationships among human RNA, demographic variables, and primary human cell physiology. The platelet RNA and expression-1 study measured platelet aggregation to arachidonic acid, ADP, protease-activated receptor (PAR) 1 activation peptide (PAR1-AP), and PAR4-AP, as well as mRNA and microRNA (miRNA) levels in platelets from 84 white and 70 black healthy subjects. A total of 5911 uniquely mapped mRNAs and 181 miRNAs were commonly expressed and validated in a separate cohort. One hundred twenty-nine mRNAs and 15 miRNAs were differentially expressed (DE) by age, and targets of these miRNAs were over-represented among these mRNAs. Fifty-four mRNAs and 9 miRNAs were DE by gender. Networks of miRNAs targeting mRNAs, both DE by age and gender, were identified. The inverse relationship in these RNA pairs suggests miRNAs regulate mRNA levels on aging and between genders. A simple, interactive public web tool (www.plateletomics.com) was developed that permits queries of RNA levels and associations among RNA, platelet aggregation and demographic variables. Access to these data will facilitate discovery of mechanisms of miRNA regulation of gene expression. These results provide new insights into aging and gender, and future platelet RNA association studies must account for age and gender.

Common variants in the human platelet PAR4 thrombin receptor alter platelet function and differ by race

Human platelets express 2 thrombin receptors: protease-activated receptor (PAR)-1 and PAR4. Recently, we reported 3.7-fold increased PAR4-mediated aggregation kinetics in platelets from black subjects compared with white subjects. We now show that platelets from blacks (n = 70) express 14% more PAR4 protein than those from whites (n = 84), but this difference is not associated with platelet PAR4 function. Quantitative trait locus analysis identified 3 common single nucleotide polymorphisms in the PAR4 gene (F2RL3) associated with PAR4-induced platelet aggregation. Among these single nucleotide polymorphisms, rs773902 determines whether residue 120 in transmembrane domain 2 is an alanine (Ala) or threonine (Thr). Compared with the Ala120 variant, Thr120 was more common in black subjects than in white subjects (63% vs 19%), was associated with higher PAR4-induced human platelet aggregation and Ca2+ flux, and generated greater inositol 1,4,5-triphosphate in transfected cells. A second, less frequent F2RL3 variant, Phe296Val, was only observed in blacks and abolished the enhanced PAR4-induced platelet aggregation and 1,4,5-triphosphate generation associated with PAR4-Thr120. PAR4 genotype did not affect vorapaxar inhibition of platelet PAR1 function, but a strong pharmacogenetic effect was observed with the PAR4-specific antagonist YD-3 [1-benzyl-3(ethoxycarbonylphenyl)-indazole]. These findings may have an important pharmacogenetic effect on the development of new PAR antagonists.

A paradox of transcriptional and functional innate interferon responses of human intestinal enteroids to enteric virus infection

Significance Understanding host–enteric virus interactions has been limited by the inability to culture nontransformed small intestinal epithelial cells and to infect animal models with human viruses. We report epithelial responses in human small intestinal enteroid cultures from different individuals following infection with human rotavirus (HRV), a model enteric pathogen. RNA-sequencing and functional assays revealed type III IFN as the dominant transcriptional response that activates interferon-stimulated genes, but antagonism of the IFN response negates restriction of HRV replication. Exogenously added IFNs reduce HRV replication, with type I IFN being most effective. This highlights a paradox between the strong type III transcriptional response and the weaker functional role of type III IFN in human enteric viral restriction in human small intestinal cultures. The intestinal epithelium can limit enteric pathogens by producing antiviral cytokines, such as IFNs. Type I IFN (IFN-α/β) and type III IFN (IFN-λ) function at the epithelial level, and their respective efficacies depend on the specific pathogen and site of infection. However, the roles of type I and type III IFN in restricting human enteric viruses are poorly characterized as a result of the difficulties in cultivating these viruses in vitro and directly obtaining control and infected small intestinal human tissue. We infected nontransformed human intestinal enteroid cultures from multiple individuals with human rotavirus (HRV) and assessed the host epithelial response by using RNA-sequencing and functional assays. The dominant transcriptional pathway induced by HRV infection is a type III IFN-regulated response. Early after HRV infection, low levels of type III IFN protein activate IFN-stimulated genes. However, this endogenous response does not restrict HRV replication because replication-competent HRV antagonizes the type III IFN response at pre- and posttranscriptional levels. In contrast, exogenous IFN treatment restricts HRV replication, with type I IFN being more potent than type III IFN, suggesting that extraepithelial sources of type I IFN may be the critical IFN for limiting enteric virus replication in the human intestine.

CELF1 is a central node in post-transcriptional regulatory programmes underlying EMT

The importance of translational regulation in tumour biology is increasingly appreciated. Here, we leverage polyribosomal profiling to prospectively define translational regulatory programs underlying epithelial-to-mesenchymal transition (EMT) in breast epithelial cells. We identify a group of ten translationally regulated drivers of EMT sharing a common GU-rich cis-element within the 3′-untranslated region (3′-UTR) of their mRNA. These cis-elements, necessary for the regulatory activity imparted by these 3′-UTRs, are directly bound by the CELF1 protein, which itself is regulated post-translationally during the EMT program. CELF1 is necessary and sufficient for both mesenchymal transition and metastatic colonization, and CELF1 protein, but not mRNA, is significantly overexpressed in human breast cancer tissues. Our data present an 11-component genetic pathway, invisible to transcriptional profiling approaches, in which the CELF1 protein functions as a central node controlling translational activation of genes driving EMT and ultimately tumour progression.

EMT cells increase breast cancer metastasis via paracrine GLI activation in neighbouring tumour cells

Recent fate-mapping studies concluded that EMT is not required for metastasis of carcinomas. Here we challenge this conclusion by showing that these studies failed to account for possible crosstalk between EMT and non-EMT cells that promotes dissemination of non-EMT cells. In breast cancer models, EMT cells induce increased metastasis of weakly metastatic, non-EMT tumour cells in a paracrine manner, in part by non-cell autonomous activation of the GLI transcription factor. Treatment with GANT61, a GLI1/2 inhibitor, but not with IPI 926, a Smoothened inhibitor, blocks this effect and inhibits growth in PDX models. In human breast tumours, the EMT-transcription factors strongly correlate with activated Hedgehog/GLI signalling but not with the Hh ligands. Our findings indicate that EMT contributes to metastasis via non-cell autonomous effects that activate the Hh pathway. Although all Hh inhibitors may act against tumours with canonical Hh/GLI signalling, only GLI inhibitors would act against non-canonical EMT-induced GLI activation.

GATA2/3-TFAP2A/C transcription factor network couples human pluripotent stem cell differentiation to trophectoderm with repression of pluripotency

Significance This study provides a mechanistic explanation for the differentiation of trophoblasts from human pluripotent stem cells, a process relying on BMP morphogens. We found that a network of the transcription factors GATA2, GATA3, TFAP2A, and TFAP2C regulates early trophoblast progenitor specification by activating placental genes and inhibiting the pluripotency gene OCT4, thus acting to couple trophoblast specification with exit from pluripotency. To demonstrate the relevance of our findings in vivo, we show that down-regulating GATA3 in primate embryos prevents trophectoderm specification. In addition, we present a genome-wide analysis of active and inactive chromatin during trophoblast progenitor specification. These results provide a basis to guide investigations of human trophectoderm development. To elucidate the molecular basis of BMP4-induced differentiation of human pluripotent stem cells (PSCs) toward progeny with trophectoderm characteristics, we produced transcriptome, epigenome H3K4me3, H3K27me3, and CpG methylation maps of trophoblast progenitors, purified using the surface marker APA. We combined them with the temporally resolved transcriptome of the preprogenitor phase and of single APA+ cells. This revealed a circuit of bivalent TFAP2A, TFAP2C, GATA2, and GATA3 transcription factors, coined collectively the “trophectoderm four” (TEtra), which are also present in human trophectoderm in vivo. At the onset of differentiation, the TEtra factors occupy multiple sites in epigenetically inactive placental genes and in OCT4. Functional manipulation of GATA3 and TFAP2A indicated that they directly couple trophoblast-specific gene induction with suppression of pluripotency. In accordance, knocking down GATA3 in primate embryos resulted in a failure to form trophectoderm. The discovery of the TEtra circuit indicates how trophectoderm commitment is regulated in human embryogenesis.

Single cell RNA-seq denoising using a deep count autoencoder

Single-cell RNA sequencing (scRNA-seq) has enabled researchers to study gene expression at a cellular resolution. However, noise due to amplification and dropout may obstruct analyses, so scalable denoising methods for increasingly large but sparse scRNAseq data are needed. We propose a deep count autoencoder network (DCA) to denoise scRNA-seq datasets. DCA takes the count distribution, overdispersion and sparsity of the data into account using a zero-inflated negative binomial noise model, and nonlinear gene-gene or gene-dispersion interactions are captured. Our method scales linearly with the number of cells and can therefore be applied to datasets of millions of cells. We demonstrate that DCA denoising improves a diverse set of typical scRNA-seq data analyses using simulated and real datasets. DCA outperforms existing methods for data imputation in quality and speed, enhancing biological discovery.

MetaMap, an interactive webtool for the exploration of metatranscriptomic reads in human disease-related RNA-seq data

Motivation The MetaMap resource contains metatranscriptomic expression data from screening >17,000 RNA-seq samples from >400 archived human disease-related studies for viral and microbial reads, so-called “metafeatures”. However, navigating this set of large and heterogeneous data is challenging, especially for researchers without bioinformatic expertise. Therefore, a user-friendly interface is needed that allows users to visualize and statistically analyse the data. Results We developed an interactive frontend to facilitate the exploration of the MetaMap resource. The webtool allows users to query the resource by searching study abstracts for keywords or browsing expression patterns for specific metafeatures. Moreover, users can manually define sample groupings or use the existing annotation for downstream analysis. The web tool provides a large variety of analyses and visualizations including dimension reduction, differential abundance analysis and Krona visualizations. The MetaMap webtool represents a valuable resource for hypothesis generation regarding the impact of the microbiome in human disease. Availability The presented web tool can be accessed at https://github.com/theislab/MetaMap