Tumor-specific molecules are needed across diverse areas of oncology for use in early detection, diagnosis, prognosis and therapy. Large and growing public databases of transcriptome sequencing data (RNA-seq) derived from tumors and normal tissues hold the potential of yielding tumor-specific molecules, but because the data are new they have not been fully explored for this purpose. We have developed custom bioinformatic algorithms and used them with 296 high-grade serous ovarian (HGS-OvCa) tumor and 1,839 normal RNA-seq datasets to identify mRNA isoforms with tumor-specific expression. We rank prioritized isoforms by likelihood of being expressed in HGS-OvCa tumors and not in normal tissues and analyzed 671 top-ranked isoforms by high-throughput RT-qPCR. Six of these isoforms were expressed in a majority of the 12 tumors examined but not in 18 normal tissues. An additional 11 were expressed in most tumors and only one normal tissue, which in most cases was fallopian or colon. Of the 671 isoforms, the topmost 5% (n = 33) ranked based on having tumor-specific or highly restricted normal tissue expression by RT-qPCR analysis are enriched for oncogenic, stem cell/cancer stem cell, and early development loci--including ETV4, FOXM1, LSR, CD9, RAB11FIP4, and FGFRL1. Many of the 33 isoforms are predicted to encode proteins with unique amino acid sequences, which would allow them to be specifically targeted for one or more therapeutic strategies--including monoclonal antibodies and T-cell-based vaccines. The systematic process described herein is readily and rapidly applicable to the more than 30 additional tumor types for which sufficient amounts of RNA-seq already exist.

Abstract Background Short-term laboratory evolution of bacteria followed by genomic sequencing provides insight into the mechanism of adaptive evolution, such as the number of mutations needed for adaptation, genotype-phenotype relationships, and the reproducibility of adaptive outcomes. Results In the present study, we describe the genome sequencing of 11 endpoints of Escherichia coli that underwent 60-day laboratory adaptive evolution under growth rate selection pressure in lactate minimal media. Two to eight mutations were identified per endpoint. Generally, each endpoint acquired mutations to different genes. The most notable exception was an 82 base-pair deletion in the rph-pyrE operon that appeared in 7 of the 11 adapted strains. This mutation conferred an approximately 15% increase to the growth rate when experimentally introduced to the wild-type background and resulted in an approximately 30% increase to growth rate when introduced to a background already harboring two adaptive mutations. Additionally, most endpoints had a mutation in a regulatory gene (crp or relA, for example) or the RNA polymerase. Conclusions The 82 base-pair deletion found in the rph-pyrE operon of many endpoints may function to relieve a pyrimidine biosynthesis defect present in MG1655. In contrast, a variety of regulators acquire mutations in the different endpoints, suggesting flexibility in overcoming regulatory challenges in the adaptation.

Background

Deregulation of RNA editing by adenosine deaminases acting on dsRNA (ADARs) has been implicated in the progression of diverse human cancers including hematopoietic malignancies such as chronic myeloid leukemia (CML). Inflammation-associated activation of ADAR1 occurs in leukemia stem cells specifically in the advanced, often drug-resistant stage of CML known as blast crisis. However, detection of cancer stem cell-associated RNA editing by RNA sequencing in these rare cell populations can be technically challenging, costly and requires PCR validation. The objectives of this study were to validate RNA editing of a subset of cancer stem cell-associated transcripts, and to develop a quantitative RNA editing fingerprint assay for rapid detection of aberrant RNA editing in human malignancies.

Methods

To facilitate quantification of cancer stem cell-associated RNA editing in exons and intronic or 3'UTR primate-specific Alu sequences using a sensitive, cost-effective method, we established an in vitro RNA editing model and developed a sensitive RNA editing fingerprint assay that employs a site-specific quantitative PCR (RESSq-PCR) strategy. This assay was validated in a stably-transduced human leukemia cell line, lentiviral-ADAR1 transduced primary hematopoietic stem and progenitor cells, and in primary human chronic myeloid leukemia stem cells.

Results

In lentiviral ADAR1-expressing cells, increased RNA editing of MDM2, APOBEC3D, GLI1 and AZIN1 transcripts was detected by RESSq-PCR with improved sensitivity over sequencing chromatogram analysis. This method accurately detected cancer stem cell-associated RNA editing in primary chronic myeloid leukemia samples, establishing a cancer stem cell-specific RNA editing fingerprint of leukemic transformation that will support clinical development of novel diagnostic tools to predict and prevent cancer progression.

Conclusions

RNA editing quantification enables rapid detection of malignant progenitors signifying cancer progression and therapeutic resistance, and will aid future RNA editing inhibitor development efforts.

Mutations in the splicing factor SF3B1 are found in several cancer types and have been associated with various splicing defects. Using transcriptome sequencing data from chronic lymphocytic leukemia, breast cancer and uveal melanoma tumor samples, we show that hundreds of cryptic 3' splice sites (3'SSs) are used in cancers with SF3B1 mutations. We define the necessary sequence context for the observed cryptic 3' SSs and propose that cryptic 3'SS selection is a result of SF3B1 mutations causing a shift in the sterically protected region downstream of the branch point. While most cryptic 3'SSs are present at low frequency (<10%) relative to nearby canonical 3'SSs, we identified ten genes that preferred out-of-frame cryptic 3'SSs. We show that cancers with mutations in the SF3B1 HEAT 5-9 repeats use cryptic 3'SSs downstream of the branch point and provide both a mechanistic model consistent with published experimental data and affected targets that will guide further research into the oncogenic effects of SF3B1 mutation.

Formative research suggests that a human embryonic stem cell-specific alternative splicing gene regulatory network, which is repressed by Muscleblind-like (MBNL) RNA binding proteins, is involved in cell reprogramming. In this study, RNA sequencing, splice isoform-specific quantitative RT-PCR, lentiviral transduction, and in vivo humanized mouse model studies demonstrated that malignant reprogramming of progenitors into self-renewing blast crisis chronic myeloid leukemia stem cells (BC LSCs) was partially driven by decreased MBNL3. Lentiviral knockdown of MBNL3 resulted in reversion to an embryonic alternative splice isoform program typified by overexpression of CD44 transcript variant 3, containing variant exons 8-10, and BC LSC proliferation. Although isoform-specific lentiviral CD44v3 overexpression enhanced chronic phase chronic myeloid leukemia (CML) progenitor replating capacity, lentiviral shRNA knockdown abrogated these effects. Combined treatment with a humanized pan-CD44 monoclonal antibody and a breakpoint cluster region - ABL proto-oncogene 1, nonreceptor tyrosine kinase (BCR-ABL1) antagonist inhibited LSC maintenance in a niche-dependent manner. In summary, MBNL3 down-regulation-related reversion to an embryonic alternative splicing program, typified by CD44v3 overexpression, represents a previously unidentified mechanism governing malignant progenitor reprogramming in malignant microenvironments and provides a pivotal opportunity for selective BC LSC detection and therapeutic elimination.

Post-transcriptional adenosine-to-inosine RNA editing mediated by adenosine deaminase acting on RNA1 (ADAR1) promotes cancer progression and therapeutic resistance. However, ADAR1 editase-dependent mechanisms governing leukemia stem cell (LSC) generation have not been elucidated. In blast crisis chronic myeloid leukemia (BC CML), we show that increased JAK2 signaling and BCR-ABL1 amplification activate ADAR1. In a humanized BC CML mouse model, combined JAK2 and BCR-ABL1 inhibition prevents LSC self-renewal commensurate with ADAR1 downregulation. Lentiviral ADAR1 wild-type, but not an editing-defective ADAR1(E912A) mutant, induces self-renewal gene expression and impairs biogenesis of stem cell regulatory let-7 microRNAs. Combined RNA sequencing, qRT-PCR, CLIP-ADAR1, and pri-let-7 mutagenesis data suggest that ADAR1 promotes LSC generation via let-7 pri-microRNA editing and LIN28B upregulation. A small-molecule tool compound antagonizes ADAR1's effect on LSC self-renewal in stromal co-cultures and restores let-7 biogenesis. Thus, ADAR1 activation represents a unique therapeutic vulnerability in LSCs with active JAK2 signaling.