Spliceosome-targeted therapies trigger an antiviral immune response in triple-negative breast cancer.

Many oncogenic insults deregulate RNA splicing, often leading to hypersensitivity of tumors to spliceosome-targeted therapies (STTs). However, the mechanisms by which STTs selectively kill cancers remain largely unknown. Herein, we discover that mis-spliced RNA itself is a molecular trigger for tumor killing through viral mimicry. In MYC-driven triple-negative breast cancer, STTs cause widespread cytoplasmic accumulation of mis-spliced mRNAs, many of which form double-stranded structures. Double-stranded RNA (dsRNA)-binding proteins recognize these endogenous dsRNAs, triggering antiviral signaling and extrinsic apoptosis. In immune-competent models of breast cancer, STTs cause tumor cell-intrinsic antiviral signaling, downstream adaptive immune signaling, and tumor cell death. Furthermore, RNA mis-splicing in human breast cancers correlates with innate and adaptive immune signatures, especially in MYC-amplified tumors that are typically immune cold. These findings indicate that dsRNA-sensing pathways respond to global aberrations of RNA splicing in cancer and provoke the hypothesis that STTs may provide unexplored strategies to activate anti-tumor immune pathways.

Gene expression imputation and cell-type deconvolution in human brain with spatiotemporal precision and its implications for brain-related disorders.

As the most complex organ of the human body, the brain is composed of diverse regions, each consisting of distinct cell types and their respective cellular interactions. Human brain development involves a finely tuned cascade of interactive events. These include spatiotemporal gene expression changes and dynamic alterations in cell-type composition. However, our understanding of this process is still largely incomplete owing to the difficulty of brain spatiotemporal transcriptome collection. In this study, we developed a tensor-based approach to impute gene expression on a transcriptome-wide level. After rigorous computational benchmarking, we applied our approach to infer missing data points in the widely used BrainSpan resource and completed the entire grid of spatiotemporal transcriptomics. Next, we conducted deconvolutional analyses to comprehensively characterize major cell-type dynamics across the entire BrainSpan resource to estimate the cellular temporal changes and distinct neocortical areas across development. Moreover, integration of these results with GWAS summary statistics for 13 brain-associated traits revealed multiple novel trait-cell-type associations and trait-spatiotemporal relationships. In summary, our imputed BrainSpan transcriptomic data provide a valuable resource for the research community and our findings help further studies of the transcriptional and cellular dynamics of the human brain and related diseases.

Single-Cell RNA Sequencing Uncovers Paracrine Functions of the Epicardial-Derived Cells in Arrhythmogenic Cardiomyopathy.

BACKGROUND: Arrhythmogenic cardiomyopathy (ACM) manifests with sudden death, arrhythmias, heart failure, apoptosis, and myocardial fibro-adipogenesis. The phenotype typically starts at the epicardium and advances transmurally. Mutations in genes encoding desmosome proteins, including DSP (desmoplakin), are major causes of ACM. METHODS: To delineate contributions of the epicardium to the pathogenesis of ACM, the Dsp allele was conditionally deleted in the epicardial cells in mice upon expression of tamoxifen-inducible Cre from the Wt1 locus. Wild type (WT) and Wt1-CreERT2:DspW/F were crossed to Rosa26mT/mG (R26mT/mG) dual reporter mice to tag the epicardial-derived cells with the EGFP (enhanced green fluorescent protein) reporter protein. Tagged epicardial-derived cells from adult Wt1-CreERT2:R26mT/mG and Wt1-CreERT2: R26mT/mG:DspW/F mouse hearts were isolated by fluorescence-activated cell staining and sequenced by single-cell RNA sequencing. RESULTS: WT1 (Wilms tumor 1) expression was progressively restricted postnatally and was exclusive to the epicardium by postnatal day 21. Expression of Dsp was reduced in the epicardial cells but not in cardiac myocytes in the Wt1-CreERT2:DspW/F mice. The Wt1-CreERT2:DspW/F mice exhibited premature death, cardiac dysfunction, arrhythmias, myocardial fibro-adipogenesis, and apoptosis. Single-cell RNA sequencing of ≈18 000 EGFP-tagged epicardial-derived cells identified genotype-independent clusters of endothelial cells, fibroblasts, epithelial cells, and a very small cluster of cardiac myocytes, which were confirmed on coimmunofluorescence staining of the myocardial sections. Differentially expressed genes between the paired clusters in the 2 genotypes predicted activation of the inflammatory and mitotic pathways-including the TGFß1 (transforming growth factor ß1) and fibroblast growth factors-in the epicardial-derived fibroblast and epithelial clusters, but predicted their suppression in the endothelial cell cluster. The findings were corroborated by analysis of gene expression in the pooled RNA-sequencing data, which identified predominant dysregulation of genes involved in epithelial-mesenchymal transition, and dysregulation of 146 genes encoding the secreted proteins (secretome), including genes in the TGFß1 pathway. Activation of the TGFß1 and its colocalization with fibrosis in the Wt1-CreERT2:R26mT/mG:DspW/F mouse heart was validated by complementary methods. CONCLUSIONS: Epicardial-derived cardiac fibroblasts and epithelial cells express paracrine factors, including TGFß1 and fibroblast growth factors, which mediate epithelial-mesenchymal transition, and contribute to the pathogenesis of myocardial fibrosis, apoptosis, arrhythmias, and cardiac dysfunction in a mouse model of ACM. The findings uncover contributions of the epicardial-derived cells to the pathogenesis of ACM.

Longitudinal RNA-Seq Analysis of the Repeatability of Gene Expression and Splicing in Human Platelets Identifies a Platelet SELP Splice QTL.

RATIONALE: Longitudinal studies are required to distinguish within versus between-individual variation and repeatability of gene expression. They are uniquely positioned to decipher genetic signal from environmental noise, with potential application to gene variant and expression studies. However, longitudinal analyses of gene expression in healthy individuals-especially with regards to alternative splicing-are lacking for most primary cell types, including platelets. OBJECTIVE: To assess repeatability of gene expression and splicing in platelets and use repeatability to identify novel platelet expression quantitative trait loci (QTLs) and splice QTLs. METHODS AND RESULTS: We sequenced the transcriptome of platelets isolated repeatedly up to 4 years from healthy individuals. We examined within and between individual variation and repeatability of platelet RNA expression and exon skipping, a readily measured alternative splicing event. We find that platelet gene expression is generally stable between and within-individuals over time-with the exception of a subset of genes enriched for the inflammation gene ontology. We show an enrichment among repeatable genes for associations with heritable traits, including known and novel platelet expression QTLs. Several exon skipping events were also highly repeatable, suggesting heritable patterns of splicing in platelets. One of the most repeatable was exon 14 skipping of SELP. Accordingly, we identify rs6128 as a platelet splice QTL and define an rs6128-dependent association between SELP exon 14 skipping and race. In vitro experiments demonstrate that this single nucleotide variant directly affects exon 14 skipping and changes the ratio of transmembrane versus soluble P-selectin protein production. CONCLUSIONS: We conclude that the platelet transcriptome is generally stable over 4 years. We demonstrate the use of repeatability of gene expression and splicing to identify novel platelet expression QTLs and splice QTLs. rs6128 is a platelet splice QTL that alters SELP exon 14 skipping and soluble versus transmembrane P-selectin protein production.

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.

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

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.

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

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.

The Human Lung Cell Atlas: A High-Resolution Reference Map of the Human Lung in Health and Disease.

Lung disease accounts for every sixth death globally. Profiling the molecular state of all lung cell types in health and disease is currently revolutionizing the identification of disease mechanisms and will aid the design of novel diagnostic and personalized therapeutic regimens. Recent progress in high-throughput techniques for single-cell genomic and transcriptomic analyses has opened up new possibilities to study individual cells within a tissue, classify these into cell types, and characterize variations in their molecular profiles as a function of genetics, environment, cell-cell interactions, developmental processes, aging, or disease. Integration of these cell state definitions with spatial information allows the in-depth molecular description of cellular neighborhoods and tissue microenvironments, including the tissue resident structural and immune cells, the tissue matrix, and the microbiome. The Human Cell Atlas consortium aims to characterize all cells in the healthy human body and has prioritized lung tissue as one of the flagship projects. Here, we present the rationale, the approach, and the expected impact of a Human Lung Cell Atlas.

Integrative Multi-omic Analysis of Human Platelet eQTLs Reveals Alternative Start Site in Mitofusin 2.

Platelets play a central role in ischemic cardiovascular events. Cardiovascular disease (CVD) is a major cause of death worldwide. Numerous genome-wide association studies (GWASs) have identified loci associated with CVD risk. However, our understanding of how these variants contribute to disease is limited. Using data from the platelet RNA and expression 1 (PRAX1) study, we analyzed cis expression quantitative trait loci (eQTLs) in platelets from 154 normal human subjects. We confirmed these results in silico by performing allele-specific expression (ASE) analysis, which demonstrated that the allelic directionality of eQTLs and ASE patterns correlate significantly. Comparison of platelet eQTLs with data from the Genotype-Tissue Expression (GTEx) project revealed that a number of platelet eQTLs are platelet specific and that platelet eQTL peaks localize to the gene body at a higher rate than eQTLs from other tissues. Upon integration with data from previously published GWASs, we found that the trait-associated variant rs1474868 coincides with the eQTL peak for mitofusin 2 (MFN2). Additional experimental and computational analyses revealed that this eQTL is linked to an unannotated alternate MFN2 start site preferentially expressed in platelets. Integration of phenotype data from the PRAX1 study showed that MFN2 expression levels were significantly associated with platelet count. This study links the variant rs1474868 to a platelet-specific regulatory role for MFN2 and demonstrates the utility of integrating multi-omic data with eQTL analysis in disease-relevant tissues for interpreting GWAS results.

Effects of genetic variation in protease activated receptor 4 after an acute coronary syndrome: Analysis from the TRACER trial.

Variation in platelet response to thrombin may affect the safety and efficacy of PAR antagonism. The Thr120 variant of the common single nucleotide polymorphism (SNP) rs773902 in the protease-activated receptor (PAR) 4 gene is associated with higher platelet aggregation compared to the Ala120 variant. We investigated the relationship between the rs773902 SNP with major bleeding and ischemic events, safety, and efficacy of PAR1 inhibition in 6177 NSTE ACS patients in the TRACER trial. There was a lower rate of GUSTO moderate/severe bleeding in patients with the Thr120 variant. The difference was driven by a lower rate in the smaller homozygous group (recessive model, HR 0.13 [0.02-0.92] P = 0.042). No significant differences were observed in the ischemic outcomes. The excess in bleeding observed with PAR1 inhibition was attenuated in patients with the Thr120 variant, but the interactions were not statistically significant. In summary, lower major bleeding rates were observed in the overall TRACER cohort with the hyperreactive PAR4 Thr120 variant. The increase in bleeding with vorapaxar was attenuated with the Thr120 variant, but we could not demonstrate an interaction with PAR1 inhibition. These findings warrant further exploration, including those of African ancestry where the A allele (Thr120) frequency is ~65%.