MYC Induces Oncogenic Stress through RNA Decay and Ribonucleotide Catabolism in Breast Cancer.

Upregulation of MYC is a hallmark of cancer, wherein MYC drives oncogenic gene expression and elevates total RNA synthesis across cancer cell transcriptomes. Although this transcriptional anabolism fuels cancer growth and survival, the consequences and metabolic stresses induced by excess cellular RNA are poorly understood. Herein, we discover that RNA degradation and downstream ribonucleotide catabolism is a novel mechanism of MYC-induced cancer cell death. Combining genetics and metabolomics, we find that MYC increases RNA decay through the cytoplasmic exosome, resulting in the accumulation of cytotoxic RNA catabolites and reactive oxygen species. Notably, tumor-derived exosome mutations abrogate MYC-induced cell death, suggesting excess RNA decay may be toxic to human cancers. In agreement, purine salvage acts as a compensatory pathway that mitigates MYC-induced ribonucleotide catabolism, and inhibitors of purine salvage impair MYC+ tumor progression. Together, these data suggest that MYC-induced RNA decay is an oncogenic stress that can be exploited therapeutically. Significance: MYC is the most common oncogenic driver of poor-prognosis cancers but has been recalcitrant to therapeutic inhibition. We discovered a new vulnerability in MYC+ cancer where MYC induces cell death through excess RNA decay. Therapeutics that exacerbate downstream ribonucleotide catabolism provide a therapeutically tractable approach to TNBC (Triple-negative Breast Cancer) and other MYC-driven cancers.

Leveraging Integrated RNA Sequencing to Decipher Adrenomedullin's Protective Mechanisms in Experimental Bronchopulmonary Dysplasia.

Bronchopulmonary dysplasia (BPD) is a chronic lung disease commonly affecting premature infants, with limited therapeutic options and increased long-term consequences. Adrenomedullin (Adm), a proangiogenic peptide hormone, has been found to protect rodents against experimental BPD. This study aims to elucidate the molecular and cellular mechanisms through which Adm influences BPD pathogenesis using a lipopolysaccharide (LPS)-induced model of experimental BPD in mice. Bulk RNA sequencing of Adm-sufficient (wild-type or Adm+/+) and Adm-haplodeficient (Adm+/-) mice lungs, integrated with single-cell RNA sequencing data, revealed distinct gene expression patterns and cell type alterations associated with Adm deficiency and LPS exposure. Notably, computational integration with cell atlas data revealed that Adm-haplodeficient mouse lungs exhibited gene expression signatures characteristic of increased inflammation, natural killer (NK) cell frequency, and decreased endothelial cell and type II pneumocyte frequency. Furthermore, in silico human BPD patient data analysis supported our cell type frequency finding, highlighting elevated NK cells in BPD infants. These results underscore the protective role of Adm in experimental BPD and emphasize that it is a potential therapeutic target for BPD infants with an inflammatory phenotype.

Let-7 restrains an oncogenic epigenetic circuit in AT2 cells to prevent ectopic formation of fibrogenic transitional cell intermediates and pulmonary fibrosis.

Analysis of lung alveolar type 2 (AT2) progenitor stem cells has highlighted fundamental mechanisms that direct their differentiation into alveolar type 1 cells (AT1s) in lung repair and disease. However, microRNA (miRNA) mediated post-transcriptional mechanisms which govern this nexus remain understudied. We show here that the let-7 miRNA family serves a homeostatic role in governance of AT2 quiescence, specifically by preventing the uncontrolled accumulation of AT2 transitional cells and by promoting AT1 differentiation to safeguard the lung from spontaneous alveolar destruction and fibrosis. Using mice and organoid models with genetic ablation of let-7a1/let-7f1/let-7d cluster (let-7afd) in AT2 cells, we demonstrate prevents AT1 differentiation and results in aberrant accumulation of AT2 transitional cells in progressive pulmonary fibrosis. Integration of enhanced AGO2 UV-crosslinking and immunoprecipitation sequencing (AGO2-eCLIP) with RNA-sequencing from AT2 cells uncovered the induction of direct targets of let-7 in an oncogene feed-forward regulatory network including BACH1/EZH2 which drives an aberrant fibrotic cascade. Additional analyses by CUT&RUN-sequencing revealed loss of let-7afd hampers AT1 differentiation by eliciting aberrant histone EZH2 methylation which prevents the exit of AT2 transitional cells into terminal AT1s. This study identifies let-7 as a key gatekeeper of post-transcriptional and epigenetic chromatin signals to prevent AT2-driven pulmonary fibrosis.

The HSP90-MYC-CDK9 network drives therapeutic resistance in mantle cell lymphoma.

Brexucabtagene autoleucel CAR-T therapy is highly efficacious in overcoming resistance to Bruton's tyrosine kinase inhibitors (BTKi) in mantle cell lymphoma. However, many patients relapse post CAR-T therapy with dismal outcomes. To dissect the underlying mechanisms of sequential resistance to BTKi and CAR-T therapy, we performed single-cell RNA sequencing analysis for 66 samples from 25 patients treated with BTKi and/or CAR-T therapy and conducted in-depth bioinformatics™ analysis. Our analysis revealed that MYC activity progressively increased with sequential resistance. HSP90AB1 (Heat shock protein 90 alpha family class B member 1), a MYC target, was identified as early driver of CAR-T resistance. CDK9 (Cyclin-dependent kinase 9), another MYC target, was significantly upregulated in Dual-R samples. Both HSP90AB1 and CDK9 expression were correlated with MYC activity levels. Pharmaceutical co-targeting of HSP90 and CDK9 synergistically diminished MYC activity, leading to potent anti-MCL activity. Collectively, our study revealed that HSP90-MYC-CDK9 network is the primary driving force of therapeutic resistance.

Investigating Cellular Trajectories in the Severity of COVID-19 and Their Transcriptional Programs Using Machine Learning Approaches.

Single-cell RNA sequencing of the bronchoalveolar lavage fluid (BALF) samples from COVID-19 patients has enabled us to examine gene expression changes of human tissue in response to the SARS-CoV-2 virus infection. However, the underlying mechanisms of COVID-19 pathogenesis at single-cell resolution, its transcriptional drivers, and dynamics require further investigation. In this study, we applied machine learning algorithms to infer the trajectories of cellular changes and identify their transcriptional programs. Our study generated cellular trajectories that show the COVID-19 pathogenesis of healthy-to-moderate and healthy-to-severe on macrophages and T cells, and we observed more diverse trajectories in macrophages compared to T cells. Furthermore, our deep-learning algorithm DrivAER identified several pathways (e.g., xenobiotic pathway and complement pathway) and transcription factors (e.g., MITF and GATA3) that could be potential drivers of the transcriptomic changes for COVID-19 pathogenesis and the markers of the COVID-19 severity. Moreover, macrophages-related functions corresponded more to the disease severity compared to T cells-related functions. Our findings more proficiently dissected the transcriptomic changes leading to the severity of a COVID-19 infection.

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.

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.

Charting Extracellular Transcriptomes in The Human Biofluid RNA Atlas.

Extracellular RNAs present in biofluids have emerged as potential biomarkers for disease. Where most studies focus on blood-derived fluids, other biofluids may be more informative. We present an atlas of messenger, circular, and small RNA transcriptomes of a comprehensive collection of 20 human biofluids. By means of synthetic spike-in controls, we compare RNA content across biofluids, revealing a 10,000-fold difference in concentration. The circular RNA fraction is increased in most biofluids compared to tissues. Each biofluid transcriptome is enriched for RNA molecules derived from specific tissues and cell types. Our atlas enables an informed selection of the most relevant biofluid to monitor particular diseases. To verify the biomarker potential in these biofluids, four validation cohorts representing a broad spectrum of diseases were profiled, revealing numerous differential RNAs between case and control subjects. Spike-normalized data are publicly available in the R2 web portal for further exploration.

Alveolar regeneration through a Krt8+ transitional stem cell state that persists in human lung fibrosis.

The cell type specific sequences of transcriptional programs during lung regeneration have remained elusive. Using time-series single cell RNA-seq of the bleomycin lung injury model, we resolved transcriptional dynamics for 28 cell types. Trajectory modeling together with lineage tracing revealed that airway and alveolar stem cells converge on a unique Krt8 + transitional stem cell state during alveolar regeneration. These cells have squamous morphology, feature p53 and NFkB activation and display transcriptional features of cellular senescence. The Krt8+ state appears in several independent models of lung injury and persists in human lung fibrosis, creating a distinct cell-cell communication network with mesenchyme and macrophages during repair. We generated a model of gene regulatory programs leading to Krt8+ transitional cells and their terminal differentiation to alveolar type-1 cells. We propose that in lung fibrosis, perturbed molecular checkpoints on the way to terminal differentiation can cause aberrant persistence of regenerative intermediate stem cell states.

Targeting the Mevalonate Pathway to Overcome Acquired Anti-HER2 Treatment Resistance in Breast Cancer.

Despite effective strategies, resistance in HER2+ breast cancer remains a challenge. While the mevalonate pathway (MVA) is suggested to promote cell growth and survival, including in HER2+ models, its potential role in resistance to HER2-targeted therapy is unknown. Parental HER2+ breast cancer cells and their lapatinib-resistant and lapatinib + trastuzumab-resistant derivatives were used for this study. MVA activity was found to be increased in lapatinib-resistant and lapatinib + trastuzumab-resistant cells. Specific blockade of this pathway with lipophilic but not hydrophilic statins and with the N-bisphosphonate zoledronic acid led to apoptosis and substantial growth inhibition of R cells. Inhibition was rescued by mevalonate or the intermediate metabolites farnesyl pyrophosphate or geranylgeranyl pyrophosphate, but not cholesterol. Activated Yes-associated protein (YAP)/transcriptional coactivator with PDZ-binding motif (TAZ) and mTORC1 signaling, and their downstream target gene product Survivin, were inhibited by MVA blockade, especially in the lapatinib-resistant/lapatinib + trastuzumab-resistant models. Overexpression of constitutively active YAP rescued Survivin and phosphorylated-S6 levels, despite blockade of the MVA. These results suggest that the MVA provides alternative signaling leading to cell survival and resistance by activating YAP/TAZ-mTORC1-Survivin signaling when HER2 is blocked, suggesting novel therapeutic targets. MVA inhibitors including lipophilic statins and N-bisphosphonates may circumvent resistance to anti-HER2 therapy warranting further clinical investigation. IMPLICATIONS: The MVA was found to constitute an escape mechanism of survival and growth in HER2+ breast cancer models resistant to anti-HER2 therapies. MVA inhibitors such as simvastatin and zoledronic acid are potential therapeutic agents to resensitize the tumors that depend on the MVA to progress on anti-HER2 therapies.

A cellular census of human lungs identifies novel cell states in health and in asthma.

Human lungs enable efficient gas exchange and form an interface with the environment, which depends on mucosal immunity for protection against infectious agents. Tightly controlled interactions between structural and immune cells are required to maintain lung homeostasis. Here, we use single-cell transcriptomics to chart the cellular landscape of upper and lower airways and lung parenchyma in healthy lungs, and lower airways in asthmatic lungs. We report location-dependent airway epithelial cell states and a novel subset of tissue-resident memory T cells. In the lower airways of patients with asthma, mucous cell hyperplasia is shown to stem from a novel mucous ciliated cell state, as well as goblet cell hyperplasia. We report the presence of pathogenic effector type 2 helper T cells (TH2) in asthmatic lungs and find evidence for type 2 cytokines in maintaining the altered epithelial cell states. Unbiased analysis of cell-cell interactions identifies a shift from airway structural cell communication in healthy lungs to a TH2-dominated interactome in asthmatic lungs.

PAGA: graph abstraction reconciles clustering with trajectory inference through a topology preserving map of single cells.

Single-cell RNA-seq quantifies biological heterogeneity across both discrete cell types and continuous cell transitions. Partition-based graph abstraction (PAGA) provides an interpretable graph-like map of the arising data manifold, based on estimating connectivity of manifold partitions ( https://github.com/theislab/paga ). PAGA maps preserve the global topology of data, allow analyzing data at different resolutions, and result in much higher computational efficiency of the typical exploratory data analysis workflow. We demonstrate the method by inferring structure-rich cell maps with consistent topology across four hematopoietic datasets, adult planaria and the zebrafish embryo and benchmark computational performance on one million neurons.

An atlas of the aging lung mapped by single cell transcriptomics and deep tissue proteomics.

Aging promotes lung function decline and susceptibility to chronic lung diseases, which are the third leading cause of death worldwide. Here, we use single cell transcriptomics and mass spectrometry-based proteomics to quantify changes in cellular activity states across 30 cell types and chart the lung proteome of young and old mice. We show that aging leads to increased transcriptional noise, indicating deregulated epigenetic control. We observe cell type-specific effects of aging, uncovering increased cholesterol biosynthesis in type-2 pneumocytes and lipofibroblasts and altered relative frequency of airway epithelial cells as hallmarks of lung aging. Proteomic profiling reveals extracellular matrix remodeling in old mice, including increased collagen IV and XVI and decreased Fraser syndrome complex proteins and collagen XIV. Computational integration of the aging proteome with the single cell transcriptomes predicts the cellular source of regulated proteins and creates an unbiased reference map of the aging lung.

IRE1α-XBP1s pathway promotes prostate cancer by activating c-MYC signaling.

Activation of endoplasmic reticulum (ER) stress/the unfolded protein response (UPR) has been linked to cancer, but the molecular mechanisms are poorly understood and there is a paucity of reagents to translate this for cancer therapy. Here, we report that an IRE1α RNase-specific inhibitor, MKC8866, strongly inhibits prostate cancer (PCa) tumor growth as monotherapy in multiple preclinical models in mice and shows synergistic antitumor effects with current PCa drugs. Interestingly, global transcriptomic analysis reveal that IRE1α-XBP1s pathway activity is required for c-MYC signaling, one of the most highly activated oncogenic pathways in PCa. XBP1s is necessary for optimal c-MYC mRNA and protein expression, establishing, for the first time, a direct link between UPR and oncogene activation. In addition, an XBP1-specific gene expression signature is strongly associated with PCa prognosis. Our data establish IRE1α-XBP1s signaling as a central pathway in PCa and indicate that its targeting may offer novel treatment strategies.

Publisher Correction: EMT cells increase breast cancer metastasis via paracrine GLI activation in neighbouring tumour cells.

This Article contains an error in Figure 2. In panel a, the second lane of the western blot should have been labelled 'siNT'. A correct version of Figure 2a appears in the Author Correction associated with this Article; the error has not been fixed in the original Article.

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.

SNP in human ARHGEF3 promoter is associated with DNase hypersensitivity, transcript level and platelet function, and Arhgef3 KO mice have increased mean platelet volume.

Genome-wide association studies have identified a genetic variant at 3p14.3 (SNP rs1354034) that strongly associates with platelet number and mean platelet volume in humans. While originally proposed to be intronic, analysis of mRNA expression in primary human hematopoietic subpopulations reveals that this SNP is located directly upstream of the predominantly expressed ARHGEF3 isoform in megakaryocytes (MK). We found that ARHGEF3, which encodes a Rho guanine exchange factor, is dramatically upregulated during both human and murine MK maturation. We show that the SNP (rs1354034) is located in a DNase I hypersensitive region in human MKs and is an expression quantitative locus (eQTL) associated with ARHGEF3 expression level in human platelets, suggesting that it may be the causal SNP that accounts for the variations observed in human platelet traits and ARHGEF3 expression. In vitro human platelet activation assays revealed that rs1354034 is highly correlated with human platelet activation by ADP. In order to test whether ARHGEF3 plays a role in MK development and/or platelet function, we developed an Arhgef3 KO/LacZ reporter mouse model. Reflecting changes in gene expression, LacZ expression increases during MK maturation in these mice. Although Arhgef3 KO mice have significantly larger platelets, loss of Arhgef3 does not affect baseline MK or platelets nor does it affect platelet function or platelet recovery in response to antibody-mediated platelet depletion compared to littermate controls. In summary, our data suggest that modulation of ARHGEF3 gene expression in humans with a promoter-localized SNP plays a role in human MKs and human platelet function-a finding resulting from the biological follow-up of human genetic studies. Arhgef3 KO mice partially recapitulate the human phenotype.

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.

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.