The case for using mapped exonic non-duplicate reads when reporting RNA-sequencing depth: examples from pediatric cancer datasets.

BACKGROUND: The reproducibility of gene expression measured by RNA sequencing (RNA-Seq) is dependent on the sequencing depth. While unmapped or non-exonic reads do not contribute to gene expression quantification, duplicate reads contribute to the quantification but are not informative for reproducibility. We show that mapped, exonic, non-duplicate (MEND) reads are a useful measure of reproducibility of RNA-Seq datasets used for gene expression analysis. FINDINGS: In bulk RNA-Seq datasets from 2,179 tumors in 48 cohorts, the fraction of reads that contribute to the reproducibility of gene expression analysis varies greatly. Unmapped reads constitute 1-77% of all reads (median [IQR], 3% [3-6%]); duplicate reads constitute 3-100% of mapped reads (median [IQR], 27% [13-43%]); and non-exonic reads constitute 4-97% of mapped, non-duplicate reads (median [IQR], 25% [16-37%]). MEND reads constitute 0-79% of total reads (median [IQR], 50% [30-61%]). CONCLUSIONS: Because not all reads in an RNA-Seq dataset are informative for reproducibility of gene expression measurements and the fraction of reads that are informative varies, we propose reporting a dataset's sequencing depth in MEND reads, which definitively inform the reproducibility of gene expression, rather than total, mapped, or exonic reads. We provide a Docker image containing (i) the existing required tools (RSeQC, sambamba, and samblaster) and (ii) a custom script to calculate MEND reads from RNA-Seq data files. We recommend that all RNA-Seq gene expression experiments, sensitivity studies, and depth recommendations use MEND units for sequencing depth.

Comparative RNA-seq analysis aids in diagnosis of a rare pediatric tumor.

Gliomatosis peritonei is a rare pathologic finding that is associated with ovarian teratomas and malignant mixed germ cell tumors. The occurrence of gliomatosis as a mature glial implant can impart an improved prognosis to patients with immature ovarian teratoma, making prompt and accurate diagnosis important. We describe a case of recurrent immature teratoma in a 10-yr-old female patient, in which comparative analysis of the RNA sequencing gene expression data from the patient's tumor was used effectively to aid in the diagnosis of gliomatosis peritonei.

Comparative Tumor RNA Sequencing Analysis for Difficult-to-Treat Pediatric and Young Adult Patients With Cancer.

Importance: Pediatric cancers are epigenetic diseases; therefore, considering tumor gene expression information is necessary for a complete understanding of the tumorigenic processes. Objective: To evaluate the feasibility and utility of incorporating comparative gene expression information into the precision medicine framework for difficult-to-treat pediatric and young adult patients with cancer. Design, Setting, and Participants: This cohort study was conducted as a consortium between the University of California, Santa Cruz (UCSC) Treehouse Childhood Cancer Initiative and clinical genomic trials. RNA sequencing (RNA-Seq) data were obtained from the following 4 clinical sites and analyzed at UCSC: British Columbia Children's Hospital (n = 31), Lucile Packard Children's Hospital at Stanford University (n = 80), CHOC Children's Hospital and Hyundai Cancer Institute (n = 46), and the Pacific Pediatric Neuro-Oncology Consortium (n = 24). The study dates were January 1, 2016, to March 22, 2017. Exposures: Participants underwent tumor RNA-Seq profiling as part of 4 separate clinical trials at partner hospitals. The UCSC either downloaded RNA-Seq data from a partner institution for analysis in the cloud or provided a Docker pipeline that performed the same analysis at a partner institution. The UCSC then compared each participant's tumor RNA-Seq profile with more than 11 000 uniformly analyzed tumor profiles from pediatric and young adult patients with cancer, downloaded from public data repositories. These comparisons were used to identify genes and pathways that are significantly overexpressed in each patient's tumor. Results of the UCSC analysis were presented to clinical partners. Main Outcomes and Measures: Feasibility of a third-party institution (UCSC Treehouse Childhood Cancer Initiative) to obtain tumor RNA-Seq data from patients, conduct comparative analysis, and present analysis results to clinicians; and proportion of patients for whom comparative tumor gene expression analysis provided useful clinical and biological information. Results: Among 144 samples from children and young adults (median age at diagnosis, 9 years; range, 0-26 years; 72 of 118 [61.0%] male [26 patients sex unknown]) with a relapsed, refractory, or rare cancer treated on precision medicine protocols, RNA-Seq-derived gene expression was potentially useful for 99 of 144 samples (68.8%) compared with DNA mutation information that was potentially useful for only 34 of 74 samples (45.9%). Conclusions and Relevance: This study's findings suggest that tumor RNA-Seq comparisons may be feasible and highlight the potential clinical utility of incorporating such comparisons into the clinical genomic interpretation framework for difficult-to-treat pediatric and young adult patients with cancer. The study also highlights for the first time to date the potential clinical utility of harmonized publicly available genomic data sets.

Suppression of autophagy impedes glioblastoma development and induces senescence.

The function of macroautophagy/autophagy during tumor initiation or in established tumors can be highly distinct and context-dependent. To investigate the role of autophagy in gliomagenesis, we utilized a KRAS-driven glioblastoma mouse model in which autophagy is specifically disrupted via RNAi against Atg7, Atg13 or Ulk1. Inhibition of autophagy strongly reduced glioblastoma development, demonstrating its critical role in promoting tumor formation. Further supporting this finding is the observation that tumors originating from Atg7-shRNA injections escaped the knockdown effect and thereby still underwent functional autophagy. In vitro, autophagy inhibition suppressed the capacity of KRAS-expressing glial cells to form oncogenic colonies or to survive low serum conditions. Molecular analyses revealed that autophagy-inhibited glial cells were unable to maintain active growth signaling under growth-restrictive conditions and were prone to undergo senescence. Overall, these results demonstrate that autophagy is crucial for glioma initiation and growth, and is a promising therapeutic target for glioblastoma treatment.

Role of autophagy in histone deacetylase inhibitor-induced apoptotic and nonapoptotic cell death.

Autophagy is a cellular catabolic pathway by which long-lived proteins and damaged organelles are targeted for degradation. Activation of autophagy enhances cellular tolerance to various stresses. Recent studies indicate that a class of anticancer agents, histone deacetylase (HDAC) inhibitors, can induce autophagy. One of the HDAC inhibitors, suberoylanilide hydroxamic acid (SAHA), is currently being used for treating cutaneous T-cell lymphoma and under clinical trials for multiple other cancer types, including glioblastoma. Here, we show that SAHA increases the expression of the autophagic factor LC3, and inhibits the nutrient-sensing kinase mammalian target of rapamycin (mTOR). The inactivation of mTOR results in the dephosphorylation, and thus activation, of the autophagic protein kinase ULK1, which is essential for autophagy activation during SAHA treatment. Furthermore, we show that the inhibition of autophagy by RNAi in glioblastoma cells results in an increase in SAHA-induced apoptosis. Importantly, when apoptosis is pharmacologically blocked, SAHA-induced nonapoptotic cell death can also be potentiated by autophagy inhibition. Overall, our findings indicate that SAHA activates autophagy via inhibiting mTOR and up-regulating LC3 expression; autophagy functions as a prosurvival mechanism to mitigate SAHA-induced apoptotic and nonapoptotic cell death, suggesting that targeting autophagy might improve the therapeutic effects of SAHA.

ULK1.ATG13.FIP200 complex mediates mTOR signaling and is essential for autophagy.

Autophagy is a degradative process that recycles long-lived and faulty cellular components. It is linked to many diseases and is required for normal development. ULK1, a mammalian serine/threonine protein kinase, plays a key role in the initial stages of autophagy, though the exact molecular mechanism is unknown. Here we report identification of a novel protein complex containing ULK1 and two additional protein factors, FIP200 and ATG13, all of which are essential for starvation-induced autophagy. Both FIP200 and ATG13 are critical for correct localization of ULK1 to the pre-autophagosome and stability of ULK1 protein. Additionally, we demonstrate by using both cellular experiments and a de novo in vitro reconstituted reaction that FIP200 and ATG13 can enhance ULK1 kinase activity individually but both are required for maximal stimulation. Further, we show that ATG13 and ULK1 are phosphorylated by the mTOR pathway in a nutrient starvation-regulated manner, indicating that the ULK1.ATG13.FIP200 complex acts as a node for integrating incoming autophagy signals into autophagosome biogenesis.