TEQUILA-seq: a versatile and low-cost method for targeted long-read RNA sequencing.

Long-read RNA sequencing (RNA-seq) is a powerful technology for transcriptome analysis, but the relatively low throughput of current long-read sequencing platforms limits transcript coverage. One strategy for overcoming this bottleneck is targeted long-read RNA-seq for preselected gene panels. We present TEQUILA-seq, a versatile, easy-to-implement, and low-cost method for targeted long-read RNA-seq utilizing isothermally linear-amplified capture probes. When performed on the Oxford nanopore platform with multiple gene panels of varying sizes, TEQUILA-seq consistently and substantially enriches transcript coverage while preserving transcript quantification. We profile full-length transcript isoforms of 468 actionable cancer genes across 40 representative breast cancer cell lines. We identify transcript isoforms enriched in specific subtypes and discover novel transcript isoforms in extensively studied cancer genes such as TP53. Among cancer genes, tumor suppressor genes (TSGs) are significantly enriched for aberrant transcript isoforms targeted for degradation via mRNA nonsense-mediated decay, revealing a common RNA-associated mechanism for TSG inactivation. TEQUILA-seq reduces the per-reaction cost of targeted capture by 2-3 orders of magnitude, as compared to a standard commercial solution. TEQUILA-seq can be broadly used for targeted sequencing of full-length transcripts in diverse biomedical research settings.

The RNA-binding proteins hnRNP H and F regulate splicing of a MYC-dependent HRAS exon in prostate cancer cells.

The MYC proto-oncogene contributes to the pathogenesis of more than half of human cancers. Malignant transformation by MYC transcriptionally up-regulates the core pre-mRNA splicing machinery and causes misregulation of alternative splicing. However, our understanding of how splicing changes are directed by MYC is limited. We performed a signaling pathway-guided splicing analysis to identify MYC-dependent splicing events. These included an HRAS cassette exon repressed by MYC across multiple tumor types. To molecularly dissect the regulation of this HRAS exon, we used antisense oligonucleotide tiling to identify splicing enhancers and silencers in its flanking introns. RNA-binding motif prediction indicated multiple binding sites for hnRNP H and hnRNP F within these cis-regulatory elements. Using siRNA knockdown and cDNA expression, we found that both hnRNP H and F activate the HRAS cassette exon. Mutagenesis and targeted RNA immunoprecipitation implicate two downstream G-rich elements in this splicing activation. Analyses of ENCODE RNA-seq datasets confirmed hnRNP H regulation of HRAS splicing. Analyses of RNA-seq datasets across multiple cancers showed a negative correlation of HNRNPH gene expression with MYC hallmark enrichment, consistent with the effect of hnRNP H on HRAS splicing. Interestingly, HNRNPF expression showed a positive correlation with MYC hallmarks and thus was not consistent with the observed effects of hnRNP F. Loss of hnRNP H/F altered cell cycle progression and induced apoptosis in the PC3 prostate cancer cell line. Collectively, our results reveal mechanisms for MYC-dependent regulation of splicing and point to possible therapeutic targets in prostate cancers.

IRIS: Discovery of cancer immunotherapy targets arising from pre-mRNA alternative splicing.

Alternative splicing (AS) is prevalent in cancer, generating an extensive but largely unexplored repertoire of novel immunotherapy targets. We describe Isoform peptides from RNA splicing for Immunotherapy target Screening (IRIS), a computational platform capable of discovering AS-derived tumor antigens (TAs) for T cell receptor (TCR) and chimeric antigen receptor T cell (CAR-T) therapies. IRIS leverages large-scale tumor and normal transcriptome data and incorporates multiple screening approaches to discover AS-derived TAs with tumor-associated or tumor-specific expression. In a proof-of-concept analysis integrating transcriptomics and immunopeptidomics data, we showed that hundreds of IRIS-predicted TCR targets are presented by human leukocyte antigen (HLA) molecules. We applied IRIS to RNA-seq data of neuroendocrine prostate cancer (NEPC). From 2,939 NEPC-associated AS events, IRIS predicted 1,651 epitopes from 808 events as potential TCR targets for two common HLA types (A*02:01 and A*03:01). A more stringent screening test prioritized 48 epitopes from 20 events with "neoantigen-like" NEPC-specific expression. Predicted epitopes are often encoded by microexons of ≤30 nucleotides. To validate the immunogenicity and T cell recognition of IRIS-predicted TCR epitopes, we performed in vitro T cell priming in combination with single-cell TCR sequencing. Seven TCRs transduced into human peripheral blood mononuclear cells (PBMCs) showed high activity against individual IRIS-predicted epitopes, providing strong evidence of isolated TCRs reactive to AS-derived peptides. One selected TCR showed efficient cytotoxicity against target cells expressing the target peptide. Our study illustrates the contribution of AS to the TA repertoire of cancer cells and demonstrates the utility of IRIS for discovering AS-derived TAs and expanding cancer immunotherapies.