Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

mTOR Inhibition Enhances Delivery and Activity of Antisense Oligonucleotides in Uveal Melanoma Cells.

Uveal melanoma (UM) is the most common primary intraocular malignancy in adults. Owing to a lack of effective treatments, patients with metastatic disease have a median survival time of 6-12 months. We recently demonstrated that the Survival Associated Mitochondrial Melanoma Specific Oncogenic Non-coding RNA (SAMMSON) is essential for UM cell survival and that antisense oligonucleotide (ASO)-mediated silencing of SAMMSON impaired cell viability and tumor growth in vitro and in vivo. By screening a library of 2911 clinical stage compounds, we identified the mammalian target of rapamycin (mTOR) inhibitor GDC-0349 to synergize with SAMMSON inhibition in UM. Mechanistic studies revealed that mTOR inhibition enhanced uptake and reduced lysosomal accumulation of lipid complexed SAMMSON ASOs, improving SAMMSON knockdown and further decreasing UM cell viability. We found mTOR inhibition to also enhance target knockdown in other cancer cell lines as well as normal cells when combined with lipid nanoparticle complexed or encapsulated ASOs or small interfering RNAs (siRNAs). Our results are relevant to nucleic acid treatment in general and highlight the potential of mTOR inhibition to enhance ASO and siRNA-mediated target knockdown.

Longitudinal evaluation of serum microRNAs as biomarkers for neuroblastoma burden and therapeutic p53 reactivation.

Accurate assessment of treatment response and residual disease is indispensable for the evaluation of cancer treatment efficacy. However, performing tissue biopsies for longitudinal follow-up poses a major challenge in the management of solid tumours like neuroblastoma. In the present study, we evaluated whether circulating miRNAs are suitable to monitor neuroblastoma tumour burden and whether treatment-induced changes of miRNA abundance in the tumour are detectable in serum. We performed small RNA sequencing on longitudinally collected serum samples from mice carrying orthotopic neuroblastoma xenografts that were exposed to treatment with idasanutlin or temsirolimus. We identified 57 serum miRNAs to be differentially expressed upon xenograft tumour manifestation, out of which 21 were also found specifically expressed in the serum of human high-risk neuroblastoma patients. The murine serum levels of these 57 miRNAs correlated with tumour tissue expression and tumour volume, suggesting potential utility for monitoring tumour burden. In addition, we describe serum miRNAs that dynamically respond to p53 activation following treatment of engrafted mice with idasanutlin. We identified idasanutlin-induced serum miRNA expression changes upon one day and 11 days of treatment. By limiting to miRNAs with a tumour-related induction, we put forward hsa-miR-34a-5p as a potential pharmacodynamic biomarker of p53 activation in serum.

A 3'-end capture sequencing method for high-throughput targeted gene expression profiling.

Molecular phenotyping through shallow 3'-end RNA-sequencing workflows is increasingly applied in the context of large-scale chemical or genetic perturbation screens to study disease biology or support drug discovery. While these workflows enable accurate quantification of the most abundant genes, they are less effective for applications that require expression profiling of low abundant transcripts, like long noncoding RNAs (lncRNAs), or selected gene panels. To tackle these issues, we describe a workflow combining 3'-end library preparation with 3'-end hybrid capture probes and shallow RNA-sequencing for cost-effective, targeted quantification of subsets of (low abundant) genes across hundreds to thousands of samples. To assess the performance of the method, we designed a capture probe set for more than 100 mRNA and lncRNA target genes and applied the workflow to a cohort of 360 samples. When compared to standard 3'-end RNA-sequencing, 3'-end capture sequencing resulted in a more than 200-fold enrichment of target gene abundance while conserving relative intergene and intersample abundances. 3'-end RNA capture sequencing enables accurate targeted gene expression profiling at extremely shallow sequencing depth.

The long non-coding RNA SAMMSON is essential for uveal melanoma cell survival.

Long non-coding RNAs (lncRNAs) can exhibit cell-type and cancer-type specific expression profiles, making them highly attractive as therapeutic targets. Pan-cancer RNA sequencing data revealed broad expression of the SAMMSON lncRNA in uveal melanoma (UM), the most common primary intraocular malignancy in adults. Currently, there are no effective treatments for UM patients with metastatic disease, resulting in a median survival time of 6-12 months. We aimed to investigate the therapeutic potential of SAMMSON inhibition in UM. Antisense oligonucleotide (ASO)-mediated SAMMSON inhibition impaired the growth and viability of a genetically diverse panel of uveal melanoma cell lines. These effects were accompanied by an induction of apoptosis and were recapitulated in two uveal melanoma patient derived xenograft (PDX) models through subcutaneous ASO delivery. SAMMSON pulldown revealed several candidate interaction partners, including various proteins involved in mitochondrial translation. Consequently, inhibition of SAMMSON impaired global, mitochondrial and cytosolic protein translation levels and mitochondrial function in uveal melanoma cells. The present study demonstrates that SAMMSON expression is essential for uveal melanoma cell survival. ASO-mediated silencing of SAMMSON may provide an effective treatment strategy to treat primary and metastatic uveal melanoma patients.

MISpheroID: a knowledgebase and transparency tool for minimum information in spheroid identity.

Spheroids are three-dimensional cellular models with widespread basic and translational application across academia and industry. However, methodological transparency and guidelines for spheroid research have not yet been established. The MISpheroID Consortium developed a crowdsourcing knowledgebase that assembles the experimental parameters of 3,058 published spheroid-related experiments. Interrogation of this knowledgebase identified heterogeneity in the methodological setup of spheroids. Empirical evaluation and interlaboratory validation of selected variations in spheroid methodology revealed diverse impacts on spheroid metrics. To facilitate interpretation, stimulate transparency and increase awareness, the Consortium defines the MISpheroID string, a minimum set of experimental parameters required to report spheroid research. Thus, MISpheroID combines a valuable resource and a tool for three-dimensional cellular models to mine experimental parameters and to improve reproducibility.

Validation of Circular RNAs Using RT-qPCR After Effective Removal of Linear RNAs by Ribonuclease R.

Circular RNAs (circRNAs) are a class of endogenous noncoding RNAs that have been shown to play a role in normal development, homeostasis, and disease, including cancer. CircRNAs are formed through a process called back-splicing, which results in a covalently closed loop with a nonlinear back-spliced junction (BSJ). In general, circRNA BSJs are predicted in RNA sequencing data using one of numerous circRNA detection algorithms. Selected circRNAs are then typically validated using an orthogonal method such as reverse transcription quantitative PCR (RT-qPCR) with circRNA-specific primers. However, linear transcripts originating from endogenous trans-splicing can lead to false-positive signals both in RNA sequencing and in RT-qPCR experiments. Therefore, it is essential to perform the RT-qPCR validation step only after linear RNAs have been degraded using an exonuclease such as ribonuclease R (RNase R). Several RNase R protocols are available for circRNA detection using RNA sequencing or RT-qPCR. These protocols-which vary in enzyme concentration, RNA input amount, incubation times, and cleanup steps-typically lack a detailed validated standard protocol and fail to provide a range of conditions that deliver accurate results. As such, some protocols use RNase R concentrations that are too high, resulting in partial degradation of the target circRNAs. Here, we describe an optimized workflow for circRNA validation, combining RNase R treatment and RT-qPCR. First, we outline the steps for circRNA primer design and qPCR assay validation. Then, we describe RNase R treatment of total RNA and, importantly, a subsequent essential buffer cleanup step. Lastly, we outline the steps to perform the RT-qPCR and discuss the downstream data analyses. © 2021 Wiley Periodicals LLC. Basic Protocol 1: CircRNA primer design and qPCR assay validation Basic Protocol 2: RNase R treatment, cleanup, and RT-qPCR.

Charting Extracellular Transcriptomes in The Human Biofluid RNA Atlas.

Extracellular RNAs present in biofluids have emerged as potential biomarkers for disease. Where most studies focus on blood-derived fluids, other biofluids may be more informative. We present an atlas of messenger, circular, and small RNA transcriptomes of a comprehensive collection of 20 human biofluids. By means of synthetic spike-in controls, we compare RNA content across biofluids, revealing a 10,000-fold difference in concentration. The circular RNA fraction is increased in most biofluids compared to tissues. Each biofluid transcriptome is enriched for RNA molecules derived from specific tissues and cell types. Our atlas enables an informed selection of the most relevant biofluid to monitor particular diseases. To verify the biomarker potential in these biofluids, four validation cohorts representing a broad spectrum of diseases were profiled, revealing numerous differential RNAs between case and control subjects. Spike-normalized data are publicly available in the R2 web portal for further exploration.

Author Correction: Integrative analysis identifies lincRNAs up- and downstream of neuroblastoma driver genes.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

SMARTer single cell total RNA sequencing.

Single cell RNA sequencing methods have been increasingly used to understand cellular heterogeneity. Nevertheless, most of these methods suffer from one or more limitations, such as focusing only on polyadenylated RNA, sequencing of only the 3' end of the transcript, an exuberant fraction of reads mapping to ribosomal RNA, and the unstranded nature of the sequencing data. Here, we developed a novel single cell strand-specific total RNA library preparation method addressing all the aforementioned shortcomings. Our method was validated on a microfluidics system using three different cancer cell lines undergoing a chemical or genetic perturbation and on two other cancer cell lines sorted in microplates. We demonstrate that our total RNA-seq method detects an equal or higher number of genes compared to classic polyA[+] RNA-seq, including novel and non-polyadenylated genes. The obtained RNA expression patterns also recapitulate the expected biological signal. Inherent to total RNA-seq, our method is also able to detect circular RNAs. Taken together, SMARTer single cell total RNA sequencing is very well suited for any single cell sequencing experiment in which transcript level information is needed beyond polyadenylated genes.

Integrative analysis identifies lincRNAs up- and downstream of neuroblastoma driver genes.

Long intergenic non-coding RNAs (lincRNAs) are emerging as integral components of signaling pathways in various cancer types. In neuroblastoma, only a handful of lincRNAs are known as upstream regulators or downstream effectors of oncogenes. Here, we exploit RNA sequencing data of primary neuroblastoma tumors, neuroblast precursor cells, neuroblastoma cell lines and various cellular perturbation model systems to define the neuroblastoma lincRNome and map lincRNAs up- and downstream of neuroblastoma driver genes MYCN, ALK and PHOX2B. Each of these driver genes controls the expression of a particular subset of lincRNAs, several of which are associated with poor survival and are differentially expressed in neuroblastoma tumors compared to neuroblasts. By integrating RNA sequencing data from both primary tumor tissue and cancer cell lines, we demonstrate that several of these lincRNAs are expressed in stromal cells. Deconvolution of primary tumor gene expression data revealed a strong association between stromal cell composition and driver gene status, resulting in differential expression of these lincRNAs. We also explored lincRNAs that putatively act upstream of neuroblastoma driver genes, either as presumed modulators of driver gene activity, or as modulators of effectors regulating driver gene expression. This analysis revealed strong associations between the neuroblastoma lincRNAs MIAT and MEG3 and MYCN and PHOX2B activity or expression. Together, our results provide a comprehensive catalogue of the neuroblastoma lincRNome, highlighting lincRNAs up- and downstream of key neuroblastoma driver genes. This catalogue forms a solid basis for further functional validation of candidate neuroblastoma lincRNAs.

Prognostic and Therapeutic Implications of Circulating Androgen Receptor Gene Copy Number in Prostate Cancer Patients Using Droplet Digital Polymerase Chain Reaction.

BACKGROUND: Resistance mechanisms in the androgen receptor (AR) signaling pathway remain key drivers in the progression to castration-resistant prostate cancer (CRPC) and relapse under antihormonal therapy. MATERIALS AND METHODS: We evaluated the circulating AR gene copy number (CN) gain using droplet digital polymerase chain reaction in 21 control and 91 prostate cancer serum samples and its prognostic and therapeutic implications in prostate cancer. RESULTS: In CRPC, AR CN gain was associated with faster progression to CRPC (P = .026), a greater number of previous treatments (P = .045), and previous chemotherapy (P = .016). Comparing patients with and without CN gain, the median progression-free survival (PFS) in the abiraterone subgroup was 1.7 months versus not reached (P = .004), and the median overall survival (OS) was 7 months versus 20.9 months (P = .020). In the enzalutamide subgroup, PFS was 1.7 versus 10.8 months (P = .006), and OS was 6.1 versus 16.5 months (P = .042). In the taxane subgroup, PFS was 3.2 versus 6.5 months (P = .093), and OS was 3.9 months versus not reached (P = .026). The presence of more AR copies correlated with shorter androgen deprivation (P = .002), abiraterone (P = .022), enzalutamide (P = .008), and taxane (P = .039) therapy. CONCLUSION: Circulating AR CN gain predicts for a poor prognosis in CRPC. It is a promising biomarker predetermining rapid CRPC progression and predicting worse abiraterone and enzalutamide outcomes. Furthermore, it is associated with multiple previous treatments and previous chemotherapy.

Dual targeting of MDM2 and BCL2 as a therapeutic strategy in neuroblastoma.

Wild-type p53 tumor suppressor activity in neuroblastoma tumors is hampered by increased MDM2 activity, making selective MDM2 antagonists an attractive therapeutic strategy for this childhood malignancy. Since monotherapy in cancer is generally not providing long-lasting clinical responses, we here aimed to identify small molecule drugs that synergize with idasanutlin (RG7388). To this purpose we evaluated 15 targeted drugs in combination with idasanutlin in three p53 wild type neuroblastoma cell lines and identified the BCL2 inhibitor venetoclax (ABT-199) as a promising interaction partner. The venetoclax/idasanutlin combination was consistently found to be highly synergistic in a diverse panel of neuroblastoma cell lines, including cells with high MCL1 expression levels. A more pronounced induction of apoptosis was found to underlie the synergistic interaction, as evidenced by caspase-3/7 and cleaved PARP measurements. Mice carrying orthotopic xenografts of neuroblastoma cells treated with both idasanutlin and venetoclax had drastically lower tumor weights than mice treated with either treatment alone. In conclusion, these data strongly support the further evaluation of dual BCL2/MDM2 targeting as a therapeutic strategy in neuroblastoma.

Depletion of tRNA-halves enables effective small RNA sequencing of low-input murine serum samples.

The ongoing ascent of sequencing technologies has enabled researchers to gain unprecedented insights into the RNA content of biological samples. MiRNAs, a class of small non-coding RNAs, play a pivotal role in regulating gene expression. The discovery that miRNAs are stably present in circulation has spiked interest in their potential use as minimally-invasive biomarkers. However, sequencing of blood-derived samples (serum, plasma) is challenging due to the often low RNA concentration, poor RNA quality and the presence of highly abundant RNAs that dominate sequencing libraries. In murine serum for example, the high abundance of tRNA-derived small RNAs called 5' tRNA halves hampers the detection of other small RNAs, like miRNAs. We therefore evaluated two complementary approaches for targeted depletion of 5' tRNA halves in murine serum samples. Using a protocol based on biotinylated DNA probes and streptavidin coated magnetic beads we were able to selectively deplete 95% of the targeted 5' tRNA half molecules. This allowed an unbiased enrichment of the miRNA fraction resulting in a 6-fold increase of mapped miRNA reads and 60% more unique miRNAs detected. Moreover, when comparing miRNA levels in tumor-carrying versus tumor-free mice, we observed a three-fold increase in differentially expressed miRNAs.

Functional analysis of the p53 pathway in neuroblastoma cells using the small-molecule MDM2 antagonist nutlin-3.

Suppression of p53 activity is essential for proliferation and survival of tumor cells. A direct p53-activating compound, nutlin-3, was used in this study, together with p53 mutation analysis, to characterize p53 pathway defects in a set of 34 human neuroblastoma cell lines. We identified 9 cell lines (26%) with a p53 loss-of-function mutation, including 6 missense mutations, 1 nonsense mutation, 1 in-frame deletion, and 1 homozygous deletion of the 3' end of the p53 gene. Sensitivity to nutlin-3 was highly predictive of absence of p53 mutation. Signaling pathways downstream of p53 were functionally intact in 23 of 25 cell lines with wild-type p53. Knockdown and overexpression experiments revealed a potentiating effect of p14(ARF) expression on the response of neuroblastoma cells to nutlin-3. Our findings shed light on the spectrum of p53 pathway lesions in neuroblastoma cells, indicate that defects in effector molecules downstream of p53 are remarkably rare in neuroblastoma, and identify p14(ARF) as a determinant of the outcome of the response to MDM2 inhibition. These insights may prove useful for the clinical translation of evolving strategies aimed at p53 reactivation and for the development of new therapeutic approaches.

RNA pre-amplification enables large-scale RT-qPCR gene-expression studies on limiting sample amounts.

BACKGROUND: The quantitative polymerase chain reaction (qPCR) is a widely utilized method for gene-expression analysis. However, insufficient material often compromises large-scale gene-expression studies. The aim of this study is to evaluate an RNA pre-amplification method to produce micrograms of cDNA as input for qPCR. FINDINGS: The linear isothermal Ribo-SPIA pre-amplification method (WT-Ovation; NuGEN) was first evaluated by measuring the expression of 20 genes in RNA samples from six neuroblastoma cell lines and of 194 genes in two commercially available reference RNA samples before and after pre-amplification, and subsequently applied on a large panel of 738 RNA samples extracted from neuroblastoma tumours. All RNA samples were evaluated for RNA integrity and purity. Starting from 5 to 50 nanograms of total RNA the sample pre-amplification method was applied, generating approximately 5 microgams of cDNA, sufficient to measure more than 1000 target genes. The results obtained from this study show a constant yield of pre-amplified cDNA independent of the amount of input RNA; preservation of differential gene-expression after pre-amplification without introduction of substantial bias; no co-amplification of contaminating genomic DNA; no necessity to purify the pre-amplified material; and finally the importance of good RNA quality to enable pre-amplification. CONCLUSION: Application of this unbiased and easy to use sample pre-amplification technology offers great advantage to generate sufficient material for diagnostic and prognostic work-up and enables large-scale qPCR gene-expression studies using limited amounts of sample material.

Antitumor activity of the selective MDM2 antagonist nutlin-3 against chemoresistant neuroblastoma with wild-type p53.

BACKGROUND: Restoring p53 function by antagonizing its interaction with the negative regulator MDM2 is an appealing nongenotoxic approach to treating tumors with wild-type p53. Mutational inactivation of p53 is rare in neuroblastoma tumors at diagnosis and occurs in only a subset of multidrug-resistant neuroblastomas. METHODS: The antiproliferative and cytotoxic effect of nutlin-3, a small-molecule MDM2 antagonist, was examined in chemosensitive (UKF-NB-3) and matched chemoresistant neuroblastoma cells with wild-type p53 (UKF-NB-3(r)DOX20) or with mutant p53 (UKF-NB-3(r)VCR10). Activation of the p53 pathway was assessed by expression analysis of p53 target genes, flow cytometric cell cycle analysis, and apoptosis assays. Mice with established chemoresistant tumor xenografts were treated orally with nutlin-3 or vehicle control (n = 5-10 mice per group) and were used to evaluate effects on tumor growth, p53 pathway activity, and metastatic tumor burden. All statistical tests were two-sided. RESULTS: Nutlin-3 induced a similar activation of the p53 pathway in UKF-NB-3 and UKF-NB-3(r)DOX20 cells, as evidenced by increased expression of p53 target genes, G1 cell cycle arrest, and induction of apoptosis. No such response was observed in UKF-NB-3(r)VCR10 cells with mutant p53. Oral administration of nutlin-3 to UKF-NB-3(r)DOX20 xenograft-bearing mice led to inhibition of primary tumor growth (mean tumor volume after 3 weeks of treatment, nutlin-3- vs vehicle-treated mice: 772 vs 1661 mm3, difference = 890 mm3, 95% confidence interval = 469 to 1311 mm3, P < .001), p53 pathway activation, and reduction in the extent of metastatic disease. The growth of UKF-NB-3(r)VCR10 xenografts was unaffected by nutlin-3. CONCLUSIONS: Nutlin-3 activates the p53 pathway and suppresses tumor growth in this model system of chemoresistant neuroblastoma, provided that wild-type p53 is present.

Chromosome 3p microsatellite allelotyping in neuroblastoma: a report on the technical hurdles.

Pinpointing critical regions of recurrent loss may help localize tumor suppressor genes. To determine the regions of loss on chromosome 3p in neuroblastoma, we performed loss of heterozygosity analysis using 16 microsatellite markers in a series of 65 primary tumors and 29 neuroblastoma cell lines. In this study, we report the results and discuss the technical hurdles that we encountered during data generation and interpretation that are of relevance for current studies or tests employing microsatellites. To provide functional support for the implication of 3p tumor suppressor genes in this childhood malignancy, we performed a microcell-mediated chromosome 3 transfer in neuroblastoma cells.

Aberrant methylation of candidate tumor suppressor genes in neuroblastoma.

CpG island hypermethylation has been recognized as an alternative mechanism for tumor suppressor gene inactivation. In this study, we performed methylation-specific PCR (MSP) to investigate the methylation status of 10 selected tumor suppressor genes in neuroblastoma. Seven of the investigated genes (CD44, RASSF1A, CASP8, PTEN, ZMYND10, CDH1, PRDM2) showed high frequencies (> or =30%) of methylation in 33 neuroblastoma cell lines. In 42 primary neuroblastoma tumors, the frequencies of methylation were 69%, CD44; 71%, RASSF1A; 56%, CASP8; 25%, PTEN; 15%, ZMYND10; 8%, CDH1; and 0%, PRDM2. Furthermore, CASP8 and CDH1 hypermethylation was significantly associated with poor event-free survival. Meta-analysis of 115 neuroblastoma tumors demonstrated a significant correlation between CASP8 methylation and MYCN amplification. In addition, there was a correlation between ZMYND10 methylation and MYCN amplification. The MSP data, together with optimized mRNA re-expression experiments (in terms of concentration and time of treatment and use of proper reference genes) further strengthen the notion that epigenetic alterations could play a significant role in NB oncogenesis. This study thus warrants the need for a global profiling of gene promoter hypermethylation to identify genome-wide aberrantly methylated genes in order to further understand neuroblastoma pathogenesis and to identify prognostic methylation markers.