A systematic search for RNA structural switches across the human transcriptome.

RNA structural switches are key regulators of gene expression in bacteria, but their characterization in Metazoa remains limited. Here, we present SwitchSeeker, a comprehensive computational and experimental approach for systematic identification of functional RNA structural switches. We applied SwitchSeeker to the human transcriptome and identified 245 putative RNA switches. To validate our approach, we characterized a previously unknown RNA switch in the 3' untranslated region of the RORC (RAR-related orphan receptor C) transcript. In vivo dimethyl sulfate (DMS) mutational profiling with sequencing (DMS-MaPseq), coupled with cryogenic electron microscopy, confirmed its existence as two alternative structural conformations. Furthermore, we used genome-scale CRISPR screens to identify trans factors that regulate gene expression through this RNA structural switch. We found that nonsense-mediated messenger RNA decay acts on this element in a conformation-specific manner. SwitchSeeker provides an unbiased, experimentally driven method for discovering RNA structural switches that shape the eukaryotic gene expression landscape.

High Throughput Variant Libraries and Machine Learning Yield Design Rules for Retron Gene Editors.

The bacterial retron reverse transcriptase system has served as an intracellular factory for single-stranded DNA in many biotechnological applications. In these technologies, a natural retron non-coding RNA (ncRNA) is modified to encode a template for the production of custom DNA sequences by reverse transcription. The efficiency of reverse transcription is a major limiting step for retron technologies, but we lack systematic knowledge of how to improve or maintain reverse transcription efficiency while changing the retron sequence for custom DNA production. Here, we test thousands of different modifications to the retron-Eco1 ncRNA and measure DNA production in pooled variant library experiments, identifying regions of the ncRNA that are tolerant and intolerant to modification. We apply this new information to a specific application: the use of the retron to produce a precise genome editing donor in combination with a CRISPR-Cas9 RNA-guided nuclease (an editron). We use high-throughput libraries in S. cerevisiae to additionally define design rules for editrons. We extend our new knowledge of retron DNA production and editron design rules to human genome editing to achieve the highest efficiency retron-Eco1 editrons to date.

Accurate isoform quantification by joint short- and long-read RNA-sequencing.

Accurate quantification of transcript isoforms is crucial for understanding gene regulation, functional diversity, and cellular behavior. Existing RNA sequencing methods have significant limitations: short-read (SR) sequencing provides high depth but struggles with isoform deconvolution, whereas long-read (LR) sequencing offers isoform resolution at the cost of lower depth, higher noise, and technical biases. Addressing this gap, we introduce Multi-Platform Aggregation and Quantification of Transcripts (MPAQT), a generative model that combines the complementary strengths of different sequencing platforms to achieve state-of-the-art isoform-resolved transcript quantification, as demonstrated by extensive simulations and experimental benchmarks. By applying MPAQT to an in vitro model of human embryonic stem cell differentiation into cortical neurons, followed by machine learning-based modeling of transcript abundances, we show that untranslated regions (UTRs) are major determinants of isoform proportion and exon usage; this effect is mediated through isoform-specific sequence features embedded in UTRs, which likely interact with RNA-binding proteins that modulate mRNA stability. These findings highlight MPAQT's potential to enhance our understanding of transcriptomic complexity and underline the role of splicing-independent post-transcriptional mechanisms in shaping the isoform and exon usage landscape of the cell.

Stem-cell states converge in multistage cutaneous squamous cell carcinoma development.

Stem cells play a critical role in cancer development by contributing to cell heterogeneity, lineage plasticity, and drug resistance. We created gene expression networks from hundreds of mouse tissue samples (both normal and tumor) and integrated these with lineage tracing and single-cell RNA-seq, to identify convergence of cell states in premalignant tumor cells expressing markers of lineage plasticity and drug resistance. Two of these cell states representing multilineage plasticity or proliferation were inversely correlated, suggesting a mutually exclusive relationship. Treatment of carcinomas in vivo with chemotherapy repressed the proliferative state and activated multilineage plasticity whereas inhibition of differentiation repressed plasticity and potentiated responses to cell cycle inhibitors. Manipulation of this cell state transition point may provide a source of potential combinatorial targets for cancer therapy.

Direct transposition of native DNA for sensitive multimodal single-molecule sequencing.

Concurrent readout of sequence and base modifications from long unamplified DNA templates by Pacific Biosciences of California (PacBio) single-molecule sequencing requires large amounts of input material. Here we adapt Tn5 transposition to introduce hairpin oligonucleotides and fragment (tagment) limiting quantities of DNA for generating PacBio-compatible circular molecules. We developed two methods that implement tagmentation and use 90-99% less input than current protocols: (1) single-molecule real-time sequencing by tagmentation (SMRT-Tag), which allows detection of genetic variation and CpG methylation; and (2) single-molecule adenine-methylated oligonucleosome sequencing assay by tagmentation (SAMOSA-Tag), which uses exogenous adenine methylation to add a third channel for probing chromatin accessibility. SMRT-Tag of 40 ng or more human DNA (approximately 7,000 cell equivalents) yielded data comparable to gold standard whole-genome and bisulfite sequencing. SAMOSA-Tag of 30,000-50,000 nuclei resolved single-fiber chromatin structure, CTCF binding and DNA methylation in patient-derived prostate cancer xenografts and uncovered metastasis-associated global epigenome disorganization. Tagmentation thus promises to enable sensitive, scalable and multimodal single-molecule genomics for diverse basic and clinical applications.

Identification of the extracellular membrane protein ENPP3 as a major cGAMP hydrolase and innate immune checkpoint.

2'3'-Cyclic guanosine monophosphate (GMP)-AMP (cGAMP) is a second messenger synthesized upon detection of cytosolic double-stranded DNA (dsDNA) and passed between cells to facilitate downstream immune signaling. Ectonucleotide pyrophosphatase phosphodiesterase I (ENPP1), an extracellular enzyme, was the only metazoan hydrolase known to regulate cGAMP levels to dampen anti-cancer immunity. Here, we uncover ENPP3 as the second and likely the only other metazoan cGAMP hydrolase under homeostatic conditions. ENPP3 has a tissue expression pattern distinct from ENPP1's and accounts for all cGAMP hydrolysis activity in ENPP1-deficient mice. Importantly, we also show that, as with ENPP1, selectively abolishing ENPP3's cGAMP hydrolysis activity results in diminished cancer growth and metastasis of certain tumor types in a stimulator of interferon genes (STING)-dependent manner. Both ENPP1 and ENPP3 are extracellular enzymes, suggesting the dominant role that extracellular cGAMP must play as a mediator of cell-cell innate immune communication. Our work demonstrates that ENPP1 and ENPP3 non-redundantly dampen extracellular cGAMP-STING signaling, pointing to ENPP3 as a target for cancer immunotherapy.

Focal adhesion kinase-YAP signaling axis drives drug-tolerant persister cells and residual disease in lung cancer.

Targeted therapy is effective in many tumor types including lung cancer, the leading cause of cancer mortality. Paradigm defining examples are targeted therapies directed against non-small cell lung cancer (NSCLC) subtypes with oncogenic alterations in EGFR, ALK and KRAS. The success of targeted therapy is limited by drug-tolerant persister cells (DTPs) which withstand and adapt to treatment and comprise the residual disease state that is typical during treatment with clinical targeted therapies. Here, we integrate studies in patient-derived and immunocompetent lung cancer models and clinical specimens obtained from patients on targeted therapy to uncover a focal adhesion kinase (FAK)-YAP signaling axis that promotes residual disease during oncogenic EGFR-, ALK-, and KRAS-targeted therapies. FAK-YAP signaling inhibition combined with the primary targeted therapy suppressed residual drug-tolerant cells and enhanced tumor responses. This study unveils a FAK-YAP signaling module that promotes residual disease in lung cancer and mechanism-based therapeutic strategies to improve tumor response.

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.

Revealing the grammar of small RNA secretion using interpretable machine learning.

Small non-coding RNAs can be secreted through a variety of mechanisms, including exosomal sorting, in small extracellular vesicles, and within lipoprotein complexes. However, the mechanisms that govern their sorting and secretion are not well understood. Here, we present ExoGRU, a machine learning model that predicts small RNA secretion probabilities from primary RNA sequences. We experimentally validated the performance of this model through ExoGRU-guided mutagenesis and synthetic RNA sequence analysis. Additionally, we used ExoGRU to reveal cis and trans factors that underlie small RNA secretion, including known and novel RNA-binding proteins (RBPs), e.g., YBX1, HNRNPA2B1, and RBM24. We also developed a novel technique called exoCLIP, which reveals the RNA interactome of RBPs within the cell-free space. Together, our results demonstrate the power of machine learning in revealing novel biological mechanisms. In addition to providing deeper insight into small RNA secretion, this knowledge can be leveraged in therapeutic and synthetic biology applications.

Translational adaptation in breast cancer metastasis and emerging therapeutic opportunities.

Breast cancer's tendency to metastasize poses a critical barrier to effective treatment, making it a leading cause of mortality among women worldwide. A growing body of evidence is showing that translational adaptation is emerging as a key mechanism enabling cancer cells to thrive in the dynamic tumor microenvironment (TME). Here, we systematically summarize how breast cancer cells utilize translational adaptation to drive metastasis, highlighting the intricate regulation by specific translation machinery and mRNA attributes such as sequences and structures, along with the involvement of tRNAs and other trans-acting RNAs. We provide an overview of the latest findings and emerging concepts in this area, discussing their potential implications for therapeutic strategies in breast cancer.

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.

Autism genes converge on microtubule biology and RNA-binding proteins during excitatory neurogenesis.

Recent studies have identified over one hundred high-confidence (hc) autism spectrum disorder (ASD) genes. Systems biological and functional analyses on smaller subsets of these genes have consistently implicated excitatory neurogenesis. However, the extent to which the broader set of hcASD genes are involved in this process has not been explored systematically nor have the biological pathways underlying this convergence been identified. Here, we leveraged CROP-Seq to repress 87 hcASD genes in a human in vitro model of cortical neurogenesis. We identified 17 hcASD genes whose repression significantly alters developmental trajectory and results in a common cellular state characterized by disruptions in proliferation, differentiation, cell cycle, microtubule biology, and RNA-binding proteins (RBPs). We also characterized over 3,000 differentially expressed genes, 286 of which had expression profiles correlated with changes in developmental trajectory. Overall, we uncovered transcriptional disruptions downstream of hcASD gene perturbations, correlated these disruptions with distinct differentiation phenotypes, and reinforced neurogenesis, microtubule biology, and RBPs as convergent points of disruption in ASD.

ENPP1 is an innate immune checkpoint of the anticancer cGAMP-STING pathway in breast cancer.

Ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1) expression correlates with poor prognosis in many cancers, and we previously discovered that ENPP1 is the dominant hydrolase of extracellular cGAMP: a cancer-cell-produced immunotransmitter that activates the anticancer stimulator of interferon genes (STING) pathway. However, ENPP1 has other catalytic activities and the molecular and cellular mechanisms contributing to its tumorigenic effects remain unclear. Here, using single-cell RNA-seq, we show that ENPP1 in both cancer and normal tissues drives primary breast tumor growth and metastasis by dampening extracellular 2'3'-cyclic-GMP-AMP (cGAMP)-STING-mediated antitumoral immunity. ENPP1 loss-of-function in both cancer cells and normal tissues slowed primary tumor growth and abolished metastasis. Selectively abolishing the cGAMP hydrolysis activity of ENPP1 phenocopied ENPP1 knockout in a STING-dependent manner, demonstrating that restoration of paracrine cGAMP-STING signaling is the dominant anti-cancer mechanism of ENPP1 inhibition. Finally, ENPP1 expression in breast tumors deterministically predicated whether patients would remain free of distant metastasis after pembrolizumab (anti-PD-1) treatment followed by surgery. Altogether, ENPP1 blockade represents a strategy to exploit cancer-produced extracellular cGAMP for controlled local activation of STING and is therefore a promising therapeutic approach against breast cancer.

RPL22 is a tumor suppressor in MSI-high cancers and a key splicing regulator of MDM4.

Microsatellite instability high (MSI-H) tumors are malignant tumors that, despite harboring a high mutational burden, often have intact TP53. One of the most frequent mutations in MSI-H tumors is a frameshift mutation in RPL22, a ribosomal protein. Here, we identified RPL22 as a modulator of MDM4 splicing through an alternative splicing switch in exon 6. RPL22 loss increases MDM4 exon 6 inclusion, cell proliferation, and augments resistance to the MDM inhibitor Nutlin-3a. RPL22 represses expression of its paralog, RPL22L1, by mediating the splicing of a cryptic exon corresponding to a truncated transcript. Therefore, damaging mutations in RPL22 drive oncogenic MDM4 induction and reveal a common splicing circuit in MSI-H tumors that may inform therapeutic targeting of the MDM4-p53 axis and oncogenic RPL22L1 induction.

Nucleosome density shapes kilobase-scale regulation by a mammalian chromatin remodeler.

Nearly all essential nuclear processes act on DNA packaged into arrays of nucleosomes. However, our understanding of how these processes (for example, DNA replication, RNA transcription, chromatin extrusion and nucleosome remodeling) occur on individual chromatin arrays remains unresolved. Here, to address this deficit, we present SAMOSA-ChAAT: a massively multiplex single-molecule footprinting approach to map the primary structure of individual, reconstituted chromatin templates subject to virtually any chromatin-associated reaction. We apply this method to distinguish between competing models for chromatin remodeling by the essential imitation switch (ISWI) ATPase SNF2h: nucleosome-density-dependent spacing versus fixed-linker-length nucleosome clamping. First, we perform in vivo single-molecule nucleosome footprinting in murine embryonic stem cells, to discover that ISWI-catalyzed nucleosome spacing correlates with the underlying nucleosome density of specific epigenomic domains. To establish causality, we apply SAMOSA-ChAAT to quantify the activities of ISWI ATPase SNF2h and its parent complex ACF on reconstituted nucleosomal arrays of varying nucleosome density, at single-molecule resolution. We demonstrate that ISWI remodelers operate as density-dependent, length-sensing nucleosome sliders, whose ability to program DNA accessibility is dictated by single-molecule nucleosome density. We propose that the long-observed, context-specific regulatory effects of ISWI complexes can be explained in part by the sensing of nucleosome density within epigenomic domains. More generally, our approach promises molecule-precise views of the essential processes that shape nuclear physiology.

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.

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.

Computational drug repositioning for the identification of new agents to sensitize drug-resistant breast tumors across treatments and receptor subtypes.

Introduction: Drug resistance is a major obstacle in cancer treatment and can involve a variety of different factors. Identifying effective therapies for drug resistant tumors is integral for improving patient outcomes. Methods: In this study, we applied a computational drug repositioning approach to identify potential agents to sensitize primary drug resistant breast cancers. We extracted drug resistance profiles from the I-SPY 2 TRIAL, a neoadjuvant trial for early stage breast cancer, by comparing gene expression profiles of responder and non-responder patients stratified into treatments within HR/HER2 receptor subtypes, yielding 17 treatment-subtype pairs. We then used a rank-based pattern-matching strategy to identify compounds in the Connectivity Map, a database of cell line derived drug perturbation profiles, that can reverse these signatures in a breast cancer cell line. We hypothesize that reversing these drug resistance signatures will sensitize tumors to treatment and prolong survival. Results: We found that few individual genes are shared among the drug resistance profiles of different agents. At the pathway level, however, we found enrichment of immune pathways in the responders in 8 treatments within the HR+HER2+, HR+HER2-, and HR-HER2- receptor subtypes. We also found enrichment of estrogen response pathways in the non-responders in 10 treatments primarily within the hormone receptor positive subtypes. Although most of our drug predictions are unique to treatment arms and receptor subtypes, our drug repositioning pipeline identified the estrogen receptor antagonist fulvestrant as a compound that can potentially reverse resistance across 13/17 of the treatments and receptor subtypes including HR+ and triple negative. While fulvestrant showed limited efficacy when tested in a panel of 5 paclitaxel resistant breast cancer cell lines, it did increase drug response in combination with paclitaxel in HCC-1937, a triple negative breast cancer cell line. Conclusion: We applied a computational drug repurposing approach to identify potential agents to sensitize drug resistant breast cancers in the I-SPY 2 TRIAL. We identified fulvestrant as a potential drug hit and showed that it increased response in a paclitaxel-resistant triple negative breast cancer cell line, HCC-1937, when treated in combination with paclitaxel.

ENPP1 is an innate immune checkpoint of the anticancer cGAMP-STING pathway.

ENPP1 expression correlates with poor prognosis in many cancers, and we previously discovered that ENPP1 is the dominant hydrolase of extracellular cGAMP: a cancer-cell-produced immunotransmitter that activates the anticancer STING pathway. However, ENPP1 has other catalytic activities and the molecular and cellular mechanisms contributing to its tumorigenic effects remain unclear. Here, using single cell RNA-seq (scRNA-seq), we show that ENPP1 overexpression drives primary breast tumor growth and metastasis by synergistically dampening extracellular cGAMP-STING mediated antitumoral immunity and activating immunosuppressive extracellular adenosine (eADO) signaling. In addition to cancer cells, stromal and immune cells in the tumor microenvironment (TME) also express ENPP1 that restrains their response to tumor-derived cGAMP. Enpp1 loss-of-function in both cancer cells and normal tissues slowed primary tumor initiation and growth and prevented metastasis in an extracellular cGAMP- and STING-dependent manner. Selectively abolishing the cGAMP hydrolysis activity of ENPP1 phenocopied total ENPP1 knockout, demonstrating that restoration of paracrine cGAMP-STING signaling is the dominant anti-cancer mechanism of ENPP1 inhibition. Strikingly, we find that breast cancer patients with low ENPP1 expression have significantly higher immune infiltration and improved response to therapeutics impacting cancer immunity upstream or downstream of the cGAMP-STING pathway, like PARP inhibitors and anti-PD1. Altogether, selective inhibition of ENPP1's cGAMP hydrolase activity alleviates an innate immune checkpoint to boost cancer immunity and is therefore a promising therapeutic approach against breast cancer that may synergize with other cancer immunotherapies.

A bi-steric mTORC1-selective inhibitor overcomes drug resistance in breast cancer.

Activation of the PI3K-mTOR pathway is central to breast cancer pathogenesis including resistance to many targeted therapies. The mTOR kinase forms two distinct complexes, mTORC1 and mTORC2, and understanding which is required for the survival of malignant cells has been limited by tools to selectively and completely impair either subcomplex. To address this, we used RMC-6272, a bi-steric molecule with a rapamycin-like moiety linked to an mTOR active-site inhibitor that displays >25-fold selectivity for mTORC1 over mTORC2 substrates. Complete suppression of mTORC1 by RMC-6272 causes apoptosis in ER+/HER2- breast cancer cell lines, particularly in those that harbor mutations in PIK3CA or PTEN, due to inhibition of the rapamycin resistant, mTORC1 substrate 4EBP1 and reduction of the pro-survival protein MCL1. RMC-6272 reduced translation of ribosomal mRNAs, MYC target genes, and components of the CDK4/6 pathway, suggesting enhanced impairment of oncogenic pathways compared to the partial mTORC1 inhibitor everolimus. RMC-6272 maintained efficacy in hormone therapy-resistant acquired cell lines and patient-derived xenografts (PDX), showed increased efficacy in CDK4/6 inhibitor treated acquired resistant cell lines versus their parental counterparts, and was efficacious in a PDX from a patient experiencing resistance to CDK4/6 inhibition. Bi-steric mTORC1-selective inhibition may be effective in overcoming multiple forms of therapy-resistance in ER+ breast cancers.