Small but mighty: microexons in glucose homeostasis.

Many molecular mechanisms underlying blood glucose homeostasis remain elusive. Juan-Mateu et al. find that pancreatic islet cells utilize a regulatory program, originally identified in neurons, that involves alternative splicing of microexons in genes important for insulin secretion or diabetes risk.

Learning chemical sensitivity reveals mechanisms of cellular response.

Chemical probes interrogate disease mechanisms at the molecular level by linking genetic changes to observable traits. However, comprehensive chemical screens in diverse biological models are impractical. To address this challenge, we develop ChemProbe, a model that predicts cellular sensitivity to hundreds of molecular probes and drugs by learning to combine transcriptomes and chemical structures. Using ChemProbe, we infer the chemical sensitivity of cancer cell lines and tumor samples and analyze how the model makes predictions. We retrospectively evaluate drug response predictions for precision breast cancer treatment and prospectively validate chemical sensitivity predictions in new cellular models, including a genetically modified cell line. Our model interpretation analysis identifies transcriptome features reflecting compound targets and protein network modules, identifying genes that drive ferroptosis. ChemProbe is an interpretable in silico screening tool that allows researchers to measure cellular response to diverse compounds, facilitating research into molecular mechanisms of chemical sensitivity.

Functional microRNA-targeting drug discovery by graph-based deep learning.

MicroRNAs are recognized as key drivers in many cancers but targeting them with small molecules remains a challenge. We present RiboStrike, a deep-learning framework that identifies small molecules against specific microRNAs. To demonstrate its capabilities, we applied it to microRNA-21 (miR-21), a known driver of breast cancer. To ensure selectivity toward miR-21, we performed counter-screens against miR-122 and DICER. Auxiliary models were used to evaluate toxicity and rank the candidates. Learning from various datasets, we screened a pool of nine million molecules and identified eight, three of which showed anti-miR-21 activity in both reporter assays and RNA sequencing experiments. Target selectivity of these compounds was assessed using microRNA profiling and RNA sequencing analysis. The top candidate was tested in a xenograft mouse model of breast cancer metastasis, demonstrating a significant reduction in lung metastases. These results demonstrate RiboStrike's ability to nominate compounds that target the activity of miRNAs in cancer.

Integrative identification of non-coding regulatory regions driving metastatic prostate cancer.

Large-scale sequencing efforts have been undertaken to understand the mutational landscape of the coding genome. However, the vast majority of variants occur within non-coding genomic regions. We designed an integrative computational and experimental framework to identify recurrently mutated non-coding regulatory regions that drive tumor progression. Applying this framework to sequencing data from a large prostate cancer patient cohort revealed a large set of candidate drivers. We used (1) in silico analyses, (2) massively parallel reporter assays, and (3) in vivo CRISPR interference screens to systematically validate metastatic castration-resistant prostate cancer (mCRPC) drivers. One identified enhancer region, GH22I030351, acts on a bidirectional promoter to simultaneously modulate expression of the U2-associated splicing factor SF3A1 and chromosomal protein CCDC157. SF3A1 and CCDC157 promote tumor growth in vivo. We nominated a number of transcription factors, notably SOX6, to regulate expression of SF3A1 and CCDC157. Our integrative approach enables the systematic detection of non-coding regulatory regions that drive human cancers.

Systematic annotation of orphan RNAs reveals blood-accessible molecular barcodes of cancer identity and cancer-emergent oncogenic drivers.

From extrachromosomal DNA to neo-peptides, the broad reprogramming of the cancer genome leads to the emergence of molecules that are specific to the cancer state. We recently described orphan non-coding RNAs (oncRNAs) as a class of cancer-specific small RNAs with the potential to play functional roles in breast cancer progression1. Here, we report a systematic and comprehensive search to identify, annotate, and characterize cancer-emergent oncRNAs across 32 tumor types. We also leverage large-scale in vivo genetic screens in xenografted mice to functionally identify driver oncRNAs in multiple tumor types. We have not only discovered a large repertoire of oncRNAs, but also found that their presence and absence represent a digital molecular barcode that faithfully captures the types and subtypes of cancer. Importantly, we discovered that this molecular barcode is partially accessible from the cell-free space as some oncRNAs are secreted by cancer cells. In a large retrospective study across 192 breast cancer patients, we showed that oncRNAs can be reliably detected in the blood and that changes in the cell-free oncRNA burden captures both short-term and long-term clinical outcomes upon completion of a neoadjuvant chemotherapy regimen. Together, our findings establish oncRNAs as an emergent class of cancer-specific non-coding RNAs with potential roles in tumor progression and clinical utility in liquid biopsies and disease monitoring.

Systematic identification of post-transcriptional regulatory modules.

In our cells, a limited number of RNA binding proteins (RBPs) are responsible for all aspects of RNA metabolism across the entire transcriptome. To accomplish this, RBPs form regulatory units that act on specific target regulons. However, the landscape of RBP combinatorial interactions remains poorly explored. Here, we perform a systematic annotation of RBP combinatorial interactions via multimodal data integration. We build a large-scale map of RBP protein neighborhoods by generating in vivo proximity-dependent biotinylation datasets of 50 human RBPs. In parallel, we use CRISPR interference with single-cell readout to capture transcriptomic changes upon RBP knockdowns. By combining these physical and functional interaction readouts, along with the atlas of RBP mRNA targets from eCLIP assays, we generate an integrated map of functional RBP interactions. We then use this map to match RBPs to their context-specific functions and validate the predicted functions biochemically for four RBPs. This study provides a detailed map of RBP interactions and deconvolves them into distinct regulatory modules with annotated functions and target regulons. This multimodal and integrative framework provides a principled approach for studying post-transcriptional regulatory processes and enriches our understanding of their underlying mechanisms.

Learning chemical sensitivity reveals mechanisms of cellular response.

Chemical probes interrogate disease mechanisms at the molecular level by linking genetic changes to observable traits. However, comprehensive chemical screens in diverse biological models are impractical. To address this challenge, we developed ChemProbe, a model that predicts cellular sensitivity to hundreds of molecular probes and drugs by learning to combine transcriptomes and chemical structures. Using ChemProbe, we inferred the chemical sensitivity of cancer cell lines and tumor samples and analyzed how the model makes predictions. We retrospectively evaluated drug response predictions for precision breast cancer treatment and prospectively validated chemical sensitivity predictions in new cellular models, including a genetically modified cell line. Our model interpretation analysis identified transcriptome features reflecting compound targets and protein network modules, identifying genes that drive ferroptosis. ChemProbe is an interpretable in silico screening tool that allows researchers to measure cellular response to diverse compounds, facilitating research into molecular mechanisms of chemical sensitivity.

Functional microRNA-Targeting Drug Discovery by Graph-Based Deep Learning.

MicroRNAs are recognized as key drivers in many cancers, but targeting them with small molecules remains a challenge. We present RiboStrike, a deep learning framework that identifies small molecules against specific microRNAs. To demonstrate its capabilities, we applied it to microRNA-21 (miR-21), a known driver of breast cancer. To ensure the selected molecules only targeted miR-21 and not other microRNAs, we also performed a counter-screen against DICER, an enzyme involved in microRNA biogenesis. Additionally, we used auxiliary models to evaluate toxicity and select the best candidates. Using datasets from various sources, we screened a pool of nine million molecules and identified eight, three of which showed anti-miR-21 activity in both reporter assays and RNA sequencing experiments. One of these was also tested in mouse models of breast cancer, resulting in a significant reduction of lung metastases. These results demonstrate RiboStrike’s ability to effectively screen for microRNA-targeting compounds in cancer.

Integrative identification of non-coding regulatory regions driving metastatic prostate cancer.

Large-scale sequencing efforts of thousands of tumor samples have been undertaken to understand the mutational landscape of the coding genome. However, the vast majority of germline and somatic variants occur within non-coding portions of the genome. These genomic regions do not directly encode for specific proteins, but can play key roles in cancer progression, for example by driving aberrant gene expression control. Here, we designed an integrative computational and experimental framework to identify recurrently mutated non-coding regulatory regions that drive tumor progression. Application of this approach to whole-genome sequencing (WGS) data from a large cohort of metastatic castration-resistant prostate cancer (mCRPC) revealed a large set of recurrently mutated regions. We used (i) in silico prioritization of functional non-coding mutations, (ii) massively parallel reporter assays, and (iii) in vivo CRISPR-interference (CRISPRi) screens in xenografted mice to systematically identify and validate driver regulatory regions that drive mCRPC. We discovered that one of these enhancer regions, GH22I030351, acts on a bidirectional promoter to simultaneously modulate expression of U2-associated splicing factor SF3A1 and chromosomal protein CCDC157. We found that both SF3A1 and CCDC157 are promoters of tumor growth in xenograft models of prostate cancer. We nominated a number of transcription factors, including SOX6, to be responsible for higher expression of SF3A1 and CCDC157. Collectively, we have established and confirmed an integrative computational and experimental approach that enables the systematic detection of non-coding regulatory regions that drive the progression of human cancers.

An mRNA processing pathway suppresses metastasis by governing translational control from the nucleus.

Cancer cells often co-opt post-transcriptional regulatory mechanisms to achieve pathologic expression of gene networks that drive metastasis. Translational control is a major regulatory hub in oncogenesis; however, its effects on cancer progression remain poorly understood. Here, to address this, we used ribosome profiling to compare genome-wide translation efficiencies of poorly and highly metastatic breast cancer cells and patient-derived xenografts. We developed dedicated regression-based methods to analyse ribosome profiling and alternative polyadenylation data, and identified heterogeneous nuclear ribonucleoprotein C (HNRNPC) as a translational controller of a specific mRNA regulon. We found that HNRNPC is downregulated in highly metastatic cells, which causes HNRNPC-bound mRNAs to undergo 3' untranslated region lengthening and, subsequently, translational repression. We showed that modulating HNRNPC expression impacts the metastatic capacity of breast cancer cells in xenograft mouse models. In addition, the reduced expression of HNRNPC and its regulon is associated with the worse prognosis in breast cancer patient cohorts.

A sense-antisense RNA interaction promotes breast cancer metastasis via regulation of NQO1 expression.

Antisense RNAs are ubiquitous in human cells, yet their role is largely unexplored. Here we profiled antisense RNAs in the MDA-MB-231 breast cancer cell line and its highly lung metastatic derivative. We identified one antisense RNA that drives cancer progression by upregulating the redox enzyme NADPH quinone dehydrogenase 1 (NQO1), and named it NQO1-AS. Knockdown of either NQO1 or NQO1-AS reduced lung colonization in a mouse model, and investigation into the role of NQO1 indicated that it is broadly protective against oxidative damage and ferroptosis. Breast cancer cells in the lung are dependent on this pathway, and this dependence can be exploited therapeutically by inducing ferroptosis while inhibiting NQO1. Together, our findings establish a role for NQO1-AS in the progression of breast cancer by regulating its sense mRNA post-transcriptionally. Because breast cancer predominantly affects females, the disease models used in this study are of female origin and the results are primarily applicable to females.

Generation of Schwann cell derived melanocytes from hPSCs identifies pro-metastatic factors in melanoma.

The neural crest (NC) is highly multipotent and generates diverse lineages in the developing embryo. However, spatiotemporally distinct NC populations display differences in fate potential, such as increased gliogenic and parasympathetic potential from later migrating, nerve-associated Schwann cell precursors (SCPs). Interestingly, while melanogenic potential is shared by both early migrating NC and SCPs, differences in melanocyte identity resulting from differentiation through these temporally distinct progenitors have not been determined. Here, we leverage a human pluripotent stem cell (hPSC) model of NC temporal patterning to comprehensively characterize human NC heterogeneity, fate bias, and lineage development. We captured the transition of NC differentiation between temporally and transcriptionally distinct melanogenic progenitors and identified modules of candidate transcription factor and signaling activity associated with this transition. For the first time, we established a protocol for the directed differentiation of melanocytes from hPSCs through a SCP intermediate, termed trajectory 2 (T2) melanocytes. Leveraging an existing protocol for differentiating early NC-derived melanocytes, termed trajectory 1 (T1), we performed the first comprehensive comparison of transcriptional and functional differences between these distinct melanocyte populations, revealing differences in pigmentation and unique expression of transcription factors, ligands, receptors and surface markers. We found a significant link between the T2 melanocyte transcriptional signature and decreased survival in melanoma patients in the cancer genome atlas (TCGA). We performed an in vivo CRISPRi screen of T1 and T2 melanocyte signature genes in a human melanoma cell line and discovered several T2-specific markers that promote lung metastasis in mice. We further demonstrated that one of these factors, SNRPB, regulates the splicing of transcripts involved in metastasis relevant functions such as migration, cell adhesion and proliferation. Overall, this study identifies distinct developmental trajectories as a source of diversity in melanocytes and implicates the unique molecular signature of SCP-derived melanocytes in metastatic melanoma.

Postmortem Human Dura Mater Cells Exhibit Phenotypic, Transcriptomic and Genetic Abnormalities that Impact their Use for Disease Modeling.

Patient-derived cells hold great promise for precision medicine approaches in human health. Human dermal fibroblasts have been a major source of cells for reprogramming and differentiating into specific cell types for disease modeling. Postmortem human dura mater has been suggested as a primary source of fibroblasts for in vitro modeling of neurodegenerative diseases. Although fibroblast-like cells from human and mouse dura mater have been previously described, their utility for reprogramming and direct differentiation protocols has not been fully established. In this study, cells derived from postmortem dura mater are directly compared to those from dermal biopsies of living subjects. In two instances, we have isolated and compared dermal and dural cell lines from the same subject. Notably, striking differences were observed between cells of dermal and dural origin. Compared to dermal fibroblasts, postmortem dura mater-derived cells demonstrated different morphology, slower growth rates, and a higher rate of karyotype abnormality. Dura mater-derived cells also failed to express fibroblast protein markers. When dermal fibroblasts and dura mater-derived cells from the same subject were compared, they exhibited highly divergent gene expression profiles that suggest dura mater cells originated from a mixed mural lineage. Given their postmortem origin, somatic mutation signatures of dura mater-derived cells were assessed and suggest defective DNA damage repair. This study argues for rigorous karyotyping of postmortem derived cell lines and highlights limitations of postmortem human dura mater-derived cells for modeling normal biology or disease-associated pathobiology.

A prometastatic splicing program regulated by SNRPA1 interactions with structured RNA elements.

Aberrant alternative splicing is a hallmark of cancer, yet the underlying regulatory programs that control this process remain largely unknown. Here, we report a systematic effort to decipher the RNA structural code that shapes pathological splicing during breast cancer metastasis. We discovered a previously unknown structural splicing enhancer that is enriched near cassette exons with increased inclusion in highly metastatic cells. We show that the spliceosomal protein small nuclear ribonucleoprotein polypeptide A' (SNRPA1) interacts with these enhancers to promote cassette exon inclusion. This interaction enhances metastatic lung colonization and cancer cell invasion, in part through SNRPA1-mediated regulation of PLEC alternative splicing, which can be counteracted by splicing modulating morpholinos. Our findings establish a noncanonical regulatory role for SNRPA1 as a prometastatic splicing enhancer in breast cancer.

RBMS1 Suppresses Colon Cancer Metastasis through Targeted Stabilization of Its mRNA Regulon.

Identifying master regulators that drive pathologic gene expression is a key challenge in precision oncology. Here, we have developed an analytic framework, named PRADA, that identifies oncogenic RNA-binding proteins through the systematic detection of coordinated changes in their target regulons. Application of this approach to data collected from clinical samples, patient-derived xenografts, and cell line models of colon cancer metastasis revealed the RNA-binding protein RBMS1 as a suppressor of colon cancer progression. We observed that silencing RBMS1 results in increased metastatic capacity in xenograft mouse models, and that restoring its expression blunts metastatic liver colonization. We have found that RBMS1 functions as a posttranscriptional regulator of RNA stability by directly binding its target mRNAs. Together, our findings establish a role for RBMS1 as a previously unknown regulator of RNA stability and as a suppressor of colon cancer metastasis with clinical utility for risk stratification of patients. SIGNIFICANCE: By applying a new analytic approach to transcriptomic data from clinical samples and models of colon cancer progression, we have identified RBMS1 as a suppressor of metastasis and as a post-transcriptional regulator of RNA stability. Notably, RBMS1 silencing and downregulation of its targets are negatively associated with patient survival.See related commentary by Carter, p. 1261.This article is highlighted in the In This Issue feature, p. 1241.