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.

Transcriptomic signatures of human single skeletal muscle fibers in response to high-intensity interval exercise.

The heterogeneous fiber type composition of skeletal muscle makes it challenging to decipher the molecular signaling events driving the health- and performance benefits of exercise. We developed an optimized workflow for transcriptional profiling of individual human muscle fibers before, immediately after, and after 3 h of recovery from high-intensity interval cycling exercise. From a transcriptional point-of-view, we observe that there is no dichotomy in fiber activation, which could refer to a fiber being recruited or nonrecruited. Rather, the activation pattern displays a continuum with a more uniform response within fast versus slow fibers during the recovery from exercise. The transcriptome-wide response immediately after exercise is characterized by some distinct signatures for slow versus fast fibers, although the most exercise-responsive genes are common between the two fiber types. The temporal transcriptional waves further converge the gene signatures of both fiber types toward a more similar profile during the recovery from exercise. Furthermore, a large heterogeneity among all resting and exercised fibers was observed, with the principal drivers being independent of a slow/fast typology. This profound heterogeneity extends to distinct exercise responses of fibers beyond a classification based on myosin heavy chains. Collectively, our single-fiber methodological approach points to a substantial between-fiber diversity in muscle fiber responses to high-intensity interval exercise.NEW & NOTEWORTHY By development of a single-fiber transcriptomics technology, we assessed the transcriptional events in individual human skeletal muscle fibers upon high-intensity exercise. We demonstrate a large variability in transcriptional activation of fibers, with shared and distinct gene signatures for slow and fast fibers. The heterogeneous fiber-specific exercise response extends beyond this traditional slow/fast categorization. These findings expand on our understanding of exercise responses and uncover a profound between-fiber diversity in muscle fiber activation and transcriptional perturbations.

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.

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.

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.

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.

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.

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.

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.

EVI1 overexpression in t(3;17) positive myeloid malignancies results from juxtaposition of EVI1 to the MSI2 locus at 17q22.

Chromosomal translocations involving the EVI1 locus are a recurrent finding in myeloid leukemia and are associated with poor prognosis. In this study, we performed a detailed molecular characterization of the recurrent translocation t(3;17)(q26;q22) in 13 hematologic malignancies. The EVI1 gene locus was rearranged in all 13 patients and was associated with EVI1 overexpression. In 9 out of 13 patients, the 17q breakpoints clustered in a 250 kb region on band 17q22 encompassing the MSI2 (musashi homologue 2) gene. Expression analyses failed to demonstrate ectopic MSI2 expression or the presence of an MSI2/EVI1 fusion gene. In conclusion, we show for the first time that the t(3;17) is indeed a recurrent chromosomal aberration in myeloid malignancies. In keeping with findings in other recurrent 3q26 rearrangements, overexpression of the EVI1 gene appears to be the major contributor to leukemogenesis in patients with a t(3;17).

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.

Positional and functional mapping of a neuroblastoma differentiation gene on chromosome 11.

BACKGROUND: Loss of chromosome 11q defines a subset of high-stage aggressive neuroblastomas. Deletions are typically large and mapping efforts have thus far not lead to a well defined consensus region, which hampers the identification of positional candidate tumour suppressor genes. In a previous study, functional evidence for a neuroblastoma suppressor gene on chromosome 11 was obtained through microcell mediated chromosome transfer, indicated by differentiation of neuroblastoma cells with loss of distal 11q upon introduction of chromosome 11. Interestingly, some of these microcell hybrid clones were shown to harbour deletions in the transferred chromosome 11. We decided to further exploit this model system as a means to identify candidate tumour suppressor or differentiation genes located on chromosome 11. RESULTS: In a first step, we performed high-resolution array CGH DNA copy-number analysis in order to evaluate the chromosome 11 status in the hybrids. Several deletions in both parental and transferred chromosomes in the investigated microcell hybrids were observed. Subsequent correlation of these deletion events with the observed morphological changes lead to the delineation of three putative regions on chromosome 11: 11q25, 11p13-->11p15.1 and 11p15.3, that may harbour the responsible differentiation gene. CONCLUSION: Using an available model system, we were able to put forward some candidate regions that may be involved in neuroblastoma. Additional studies will be required to clarify the putative role of the genes located in these chromosomal segments in the observed differentiation phenotype specifically or in neuroblastoma pathogenesis in general.

Novel cryptic chromosomal rearrangements in childhood acute lymphoblastic leukemia detected by multiple color fluorescent in situ hybridization.

BACKGROUND AND OBJECTIVES: It is often difficult to obtain good karyotypes of cells from children with acute lymphoblastic leukemia (ALL) because of poor morphology and spreading. Detailed karyotyping can be further hampered by the presence of multiple rearrangements. Our objective was to search for cryptic rearrangements in childhood ALL. DESIGN AND METHODS: A series of eight cases of childhood ALL with at least two structural defects were selected and studied by multiple color fluorescent in situ hybridization (M-FISH). RESULTS: Four previously not reported translocations were detected: a t(14;20) (q32;q11.2) in a 3-year old girl with T-ALL, a cryptic t(7;11)(q35;q24) in association with a t(1;14)(p32;q32) in a patient with T-ALL and two translocations possibly involving the same 6q26 region on the distal end of the long arm of chromosome 6. Further FISH analysis on the t(7;11) indicated rearrangement of the TCRB locus at 7q35 suggesting that this t(7;11) leads to overexpression of an as yet unidentified gene at 11q24. This observation also triggered further screening for TCRB rearrangements in T-ALL. FISH analysis of the t(14;20) with an IGH locus-specific probe provided evidence for an unusual rearrangement of the IGH gene, in the variable gene segment region. Finally, we also observed cryptic insertions of AF4 and ETV6 in combination with complex rearrangements, leading to MLL/AF4 and ETV6/RUNX1 gene fusions. INTERPRETATION AND CONCLUSIONS: This study underscores the importance and power of M-FISH analysis in unraveling complex karyotypes and identifying cryptic chromosomal rearrangements. It also sheds some light on the implication of cryptic TCRB rearrangements in T-ALL.