The ability of a transcription factor to regulate its targets is modulated by a variety of genetic and epigenetic mechanisms. Alternative splicing can modulate gene function by adding or removing certain protein domains, and therefore affect the activity of protein. Reverse engineering of gene regulatory networks using gene expression profiles has proven valuable in dissecting the logical relationships among multiple proteins during the transcriptional regulation. However, it is unclear whether alternative splicing of certain proteins affects the activity of other transcription factors. In order to investigate the roles of alternative splicing during transcriptional regulation, we constructed a statistical model to infer whether the alternative splicing events of modulator proteins can affect the ability of key transcription factors in regulating the expression levels of their transcriptional targets. We tested our strategy in KIRC (Kidney Renal Clear Cell Carcinoma) using the RNA-seq data downloaded from TCGA (the Cancer Genomic Atlas). We identified 828of modulation relationships between the splicing levels of modulator proteins and activity levels of transcription factors. For instance, we found that the activity levels of GR (glucocorticoid receptor) protein, a key transcription factor in kidney, can be influenced by the splicing status of multiple proteins, including TP53, MDM2 (mouse double minute 2 homolog), RBM14 (RNA-binding protein 14) and SLK (STE20 like kinase). The influenced GR-targets are enriched by key cancer-related pathways, including p53 signaling pathway, TR/RXR activation, CAR/RXR activation, G1/S checkpoint regulation pathway, and G2/M DNA damage checkpoint regulation pathway. Our analysis suggests, for the first time, that exon inclusion levels of certain regulatory proteins can affect the activities of many transcription factors. Such analysis can potentially unravel a novel mechanism of how splicing variation influences the cellular function and provide important insights for how dysregulation of splicing outcome can lead to various diseases.

Neoantigen peptides arising from genetic alterations may serve as targets for personalized cancer vaccines and as positive predictors of response to immune checkpoint therapy. Mutations in genes regulating RNA splicing are common in hematological malignancies leading to dysregulated splicing and intron retention (IR). In this study, we investigated IR as a potential source of tumor neoantigens in multiple myeloma (MM) patients and the relationship of IR-induced neoantigens (IR-neoAg) with clinical outcomes. MM-specific IR events were identified in RNA-sequencing data from the Multiple Myeloma Research Foundation CoMMpass study after removing IR events that also occurred in normal plasma cells. We quantified the IR-neoAg load by assessing IR-induced novel peptides that were predicted to bind to major histocompatibility complex (MHC) molecules. We found that high IR-neoAg load was associated with poor overall survival in both newly diagnosed and relapsed MM patients. Further analyses revealed that poor outcome in MM patients with high IR-neoAg load was associated with high expression levels of T-cell co-inhibitory molecules and elevated interferon signaling activity. We also found that MM cells exhibiting high IR levels had lower MHC-II protein abundance and treatment of MM cells with a spliceosome inhibitor resulted in increased MHC-I protein abundance. Our findings suggest that IR-neoAg may represent a novel biomarker of MM patient clinical outcome and further that targeting RNA splicing may serve as a potential therapeutic strategy to prevent MM immune escape and promote response to checkpoint blockade.

While synonymous single-nucleotide variants (sSNVs) have largely been unstudied, since they do not alter protein sequence, mounting evidence suggests that they may affect RNA conformation, splicing, and the stability of nascent-mRNAs to promote various diseases. Accurately prioritizing deleterious sSNVs from a pool of neutral ones can significantly improve our ability of selecting functional genetic variants identified from various genome-sequencing projects, and, therefore, advance our understanding of disease etiology. In this study, we develop a computational algorithm to prioritize sSNVs based on their impact on mRNA splicing and protein function. In addition to genomic features that potentially affect splicing regulation, our proposed algorithm also includes dozens structural features that characterize the functions of alternatively spliced exons on protein function. Our systematical evaluation on thousands of sSNVs suggests that several structural features, including intrinsic disorder protein scores, solvent accessible surface areas, protein secondary structures, and known and predicted protein family domains, show significant differences between disease-causing and neutral sSNVs. Our result suggests that the protein structure features offer an added dimension of information while distinguishing disease-causing and neutral synonymous variants. The inclusion of structural features increases the predictive accuracy for functional sSNV prioritization.

The persistence of cancer stem cells (CSC) is believed to contribute to resistance to platinum-based chemotherapy and disease relapse in ovarian cancer, the fifth leading cause of cancer-related death among US women. HOXC transcript antisense RNA (HOTAIR) is a long, noncoding RNA (lncRNA) overexpressed in high-grade serous ovarian cancer and linked to chemoresistance. However, HOTAIR impacts chromatin dynamics in ovarian CSCs. Oncogenic lncRNA’s contributions to drug-resistant disease are incompletely understood. Here, we generated HOTAIR knockout (KO) high-grade serous ovarian cancer cell lines using paired CRISPR guide RNA design to investigate the function of HOTAIR. We show the loss of HOTAIR function resensitized ovarian cancer cells to platinum treatment and decreased the population of ovarian CSCs. Furthermore, HOTAIR KO inhibited the development of stemness-related phenotypes, including spheroid formation ability and expression of key stemness-associated genes ALDH1A1, NOTCH3, SOX9, and PROM1. HOTAIR KO altered the cellular transcriptome and chromatin accessibility landscape of multiple oncogenic-associated genes and pathways, including the NF-kB pathway. HOTAIR functions as an oncogene by recruiting enhancer of zeste homolog 2 (EZH2) to catalyze H3K27 trimethylation to suppress downstream tumor suppressor genes, and it was of interest to inhibit both HOTAIR and EZH2. In vivo, combining a HOTAIR inhibitor with an EZH2 inhibitor and platinum chemotherapy decreased tumor formation and increased survival. These results suggest a key role for HOTAIR in ovarian CSCs and malignant potential. Targeting HOTAIR in combination with epigenetic therapies may represent a therapeutic strategy to ameliorate ovarian cancer progression and resistance to platinum-based chemotherapy.

Study of genomic aberrations leading to immortalization of epithelial cells has been technically challenging due to the lack of isogenic models. To address this, we used healthy primary breast luminal epithelial cells of different genetic ancestry and their hTERT-immortalized counterparts to identify transcriptomic changes associated with immortalization. Elevated expression of TONSL (Tonsoku-like, DNA repair protein) was identified as one of the earliest events during immortalization. TONSL, which is located on chromosome 8q24.3, was found to be amplified in approximately 20% of breast cancers. TONSL alone immortalized primary breast epithelial cells and increased telomerase activity, but overexpression was insufficient for neoplastic transformation. However, TONSL-immortalized primary cells overexpressing defined oncogenes generated estrogen receptor–positive adenocarcinomas in mice. Analysis of a breast tumor microarray with approximately 600 tumors revealed poor overall and progression-free survival of patients with TONSL-overexpressing tumors. TONSL increased chromatin accessibility to pro-oncogenic transcription factors, including NF-κB and limited access to the tumor-suppressor p53. TONSL overexpression resulted in significant changes in the expression of genes associated with DNA repair hubs, including upregulation of several genes in the homologous recombination (HR) and Fanconi anemia pathways. Consistent with these results, TONSL-overexpressing primary cells exhibited upregulated DNA repair via HR. Moreover, TONSL was essential for growth of TONSL-amplified breast cancer cell lines in vivo, and these cells were sensitive to TONSL–FACT complex inhibitor CBL0137. Together, these findings identify TONSL as a regulator of epithelial cell immortalization to facilitate cancer initiation and as a target for breast cancer therapy. Significance: The chr.8q24.3 amplicon-resident gene TONSL is upregulated during the initial steps of tumorigenesis to support neoplastic transformation by increasing DNA repair and represents a potential therapeutic target for treating breast cancer.