HTSplotter: An end-to-end data processing, analysis and visualisation tool for chemical and genetic in vitro perturbation screening.

In biomedical research, high-throughput screening is often applied as it comes with automatization, higher-efficiency, and more and faster results. High-throughput screening experiments encompass drug, drug combination, genetic perturbagen or a combination of genetic and chemical perturbagen screens. These experiments are conducted in real-time assays over time or in an endpoint assay. The data analysis consists of data cleaning and structuring, as well as further data processing and visualisation, which, due to the amount of data, can easily become laborious, time-consuming and error-prone. Therefore, several tools have been developed to aid researchers in this process, but these typically focus on specific experimental set-ups and are unable to process data of several time points and genetic-chemical perturbagen screens. To meet these needs, we developed HTSplotter, a web tool and Python module that performs automatic data analysis and visualization of visualization of eitherendpoint or real-time assays from different high-throughput screening experiments: drug, drug combination, genetic perturbagen and genetic-chemical perturbagen screens. HTSplotter implements an algorithm based on conditional statements to identify experiment types and controls. After appropriate data normalization, including growth rate normalization, HTSplotter executes downstream analyses such as dose-response relationship and drug synergism assessment by the Bliss independence (BI), Zero Interaction Potency (ZIP) and Highest Single Agent (HSA) methods. All results are exported as a text file and plots are saved in a PDF file. The main advantage of HTSplotter over other available tools is the automatic analysis of genetic-chemical perturbagen screens and real-time assays where growth rate and perturbagen effect results are plotted over time. In conclusion, HTSplotter allows for the automatic end-to-end data processing, analysis and visualisation of various high-throughput in vitro cell culture screens, offering major improvements in terms of versatility, efficiency and time over existing tools.

Exploring the cell-free total RNA transcriptome in diffuse large B-cell lymphoma and primary mediastinal B-cell lymphoma patients as biomarker source in blood plasma liquid biopsies.

Introduction: Diffuse large B-cell lymphoma (DLBCL) and primary mediastinal B-cell lymphoma (PMBCL) are aggressive histological subtypes of non-Hodgkin's lymphoma. Improved understanding of the underlying molecular pathogenesis has led to new classification and risk stratification tools, including the development of cell-free biomarkers through liquid biopsies. The goal of this study was to investigate cell-free RNA (cfRNA) biomarkers in DLBCL and PMBCL patients. Materials and methods: Blood plasma samples (n=168) and matched diagnostic formalin-fixed paraffin-embedded (FFPE) tissue samples (n=69) of DLBCL patients, PMBCL patients and healthy controls were collected between 2016-2021. Plasma samples were collected at diagnosis, at interim evaluation, after treatment, and in case of refractory or relapsed disease. RNA was extracted from 200 µl plasma using the miRNeasy serum/plasma kit and from FFPE tissue using the miRNeasy FFPE kit. RNA was subsequently sequenced on a NovaSeq 6000 instrument using the SMARTer Stranded Total RNA-seq pico v3 library preparation kit. Results: Higher cfRNA concentrations were demonstrated in lymphoma patients compared to healthy controls. A large number of differentially abundant genes were identified between the cell-free transcriptomes of DLBCL patients, PMBCL patients, and healthy controls. Overlap analyses with matched FFPE samples showed that blood plasma has a unique transcriptomic profile that significantly differs from that of the tumor tissue. As a good concordance between tissue-derived gene expression and the immunohistochemistry Hans algorithm for cell-of-origin (COO) classification was demonstrated in the FFPE samples, but not in the plasma samples, a 64-gene cfRNA classifier was developed that can accurately determine COO in plasma. High plasma levels of a 9-gene signature (BECN1, PRKCB, COPA, TSC22D3, MAP2K3, UQCRHL, PTMAP4, EHD1, NAP1L1 pseudogene) and a 5-gene signature (FTH1P7, PTMAP4, ATF4, FTH1P8, ARMC7) were significantly associated with inferior progression-free and overall survival in DLBCL patients, respectively, independent of the NCCN-IPI score. Conclusion: Total RNA sequencing of blood plasma samples allows the analysis of the cell-free transcriptome in DLBCL and PMBCL patients and demonstrates its unexplored potential in identifying diagnostic, cell-of-origin, and prognostic cfRNA biomarkers.

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.

Whole transcriptome profiling of liquid biopsies from tumour xenografted mouse models enables specific monitoring of tumour-derived extracellular RNA.

While cell-free DNA (cfDNA) is widely being investigated, free circulating RNA (extracellular RNA, exRNA) has the potential to improve cancer therapy response monitoring and detection due to its dynamic nature. However, it remains unclear in which blood subcompartment tumour-derived exRNAs primarily reside. We developed a host-xenograft deconvolution framework, exRNAxeno, with mapping strategies to either a combined human-mouse reference genome or both species genomes in parallel, applicable to exRNA sequencing data from liquid biopsies of human xenograft mouse models. The tool enables to distinguish (human) tumoural RNA from (murine) host RNA, to specifically analyse tumour-derived exRNA. We applied the combined pipeline to total exRNA sequencing data from 95 blood-derived liquid biopsy samples from 30 mice, xenografted with 11 different tumours. Tumoural exRNA concentrations are not determined by plasma platelet levels, while host exRNA concentrations increase with platelet content. Furthermore, a large variability in exRNA abundance and transcript content across individual mice is observed. The tumoural gene detectability in plasma is largely correlated with the RNA expression levels in the tumour tissue or cell line. These findings unravel new aspects of tumour-derived exRNA biology in xenograft models and open new avenues to further investigate the role of exRNA in cancer.

mRNA Capture Sequencing and RT-qPCR for the Detection of Pathognomonic, Novel, and Secondary Fusion Transcripts in FFPE Tissue: A Sarcoma Showcase.

We assess the performance of mRNA capture sequencing to identify fusion transcripts in FFPE tissue of different sarcoma types, followed by RT-qPCR confirmation. To validate our workflow, six positive control tumors with a specific chromosomal rearrangement were analyzed using the TruSight RNA Pan-Cancer Panel. Fusion transcript calling by FusionCatcher confirmed these aberrations and enabled the identification of both fusion gene partners and breakpoints. Next, whole-transcriptome TruSeq RNA Exome sequencing was applied to 17 fusion gene-negative alveolar rhabdomyosarcoma (ARMS) or undifferentiated round cell sarcoma (URCS) tumors, for whom fluorescence in situ hybridization (FISH) did not identify the classical pathognomonic rearrangements. For six patients, a pathognomonic fusion transcript was readily detected, i.e., PAX3-FOXO1 in two ARMS patients, and EWSR1-FLI1, EWSR1-ERG, or EWSR1-NFATC2 in four URCS patients. For the 11 remaining patients, 11 newly identified fusion transcripts were confirmed by RT-qPCR, including COPS3-TOM1L2, NCOA1-DTNB, WWTR1-LINC01986, PLAA-MOB3B, AP1B1-CHEK2, and BRD4-LEUTX fusion transcripts in ARMS patients. Additionally, recurrently detected secondary fusion transcripts in patients diagnosed with EWSR1-NFATC2-positive sarcoma were confirmed (COPS4-TBC1D9, PICALM-SYTL2, SMG6-VPS53, and UBE2F-ALS2). In conclusion, this study shows that mRNA capture sequencing enhances the detection rate of pathognomonic fusions and enables the identification of novel and secondary fusion transcripts in sarcomas.

IL-10 producing regulatory B cells are decreased in blood from smokers and COPD patients.

BACKGROUND: Two opposing B cell subsets have been defined based on their cytokine profile: IL-6 producing effector B cells (B-effs) versus IL-10 producing regulatory B cells (B-regs) that respectively positively or negatively regulate immune responses. B-regs are decreased and/or impaired in many autoimmune diseases and inflammatory conditions. Since there is increasing evidence that links B cells and B cell-rich lymphoid follicles to the pathogenesis of COPD, the aim of this study was to investigate the presence and function of B-regs in COPD. METHODS: First, presence of IL-10 producing regulatory B cells in human lung tissue was determined by immunohistochemistry. Secondly, quantification of IL-10 + B-regs and IL-6 + B-effs in peripheral blood mononuclear cells (PBMCs) from healthy controls, smokers without airflow limitation, and COPD patients (GOLD stage I-IV) was performed by flow cytometry. Thirdly, we exposed blood-derived B cells from COPD patients in vitro to cigarette smoke extract (CSE) and quantified IL-10 + B-regs and IL-6 + B-effs. Furthermore, we aimed at restoring the perturbed IL10 production by blocking BAFF. Fourthly, we determined mRNA expression of transcription factors involved in IL-10 production in FACS sorted memory- and naive B cells upon exposure to medium or CSE. RESULTS: The presence of IL-10 producing regulatory B cells in parenchyma and lymphoid follicles in lungs was confirmed by immunohistochemistry. The percentage of IL-10 + B-regs was significantly decreased in blood-derived memory B cell subsets from smokers without airflow limitation and patients with COPD, compared to never smokers. Furthermore, the capacity of B cells to produce IL-10 was reduced upon in vitro exposure to CSE and this could not be restored by BAFF-blockade. Finally, upon CSE exposure, mRNA levels of the transcription factors IRF4 and HIF-1α, were decreased in memory B cells. CONCLUSION: Decreased numbers and impaired function of B-regs in smokers and patients with COPD might contribute to the initiation and progression of the disease.

CiLiQuant: Quantification of RNA Junction Reads Based on Their Circular or Linear Transcript Origin.

Distinguishing circular RNA reads from reads derived from the linear host transcript is a challenging task because of sequence overlap. We developed a computational approach, CiLiQuant, that determines the relative circular and linear abundance of transcripts and gene loci using back-splice and unambiguous forward-splice junction reads generated by existing mapping and circular RNA discovery tools.

CIRCprimerXL: Convenient and High-Throughput PCR Primer Design for Circular RNA Quantification.

Circular RNA (circRNA) is a class of endogenous non-coding RNA characterized by a back-splice junction (BSJ). In general, large-scale circRNA BSJ detection is performed based on RNA sequencing data, followed by the selection and validation of circRNAs of interest using RT-qPCR with circRNA-specific PCR primers. Such a primer pair is convergent and functional on the circRNA template but divergent and non-functional on the linear host gene. Although a few circRNA primer design pipelines have been published, none of them offer large-scale, easy-to-use circRNA primer design. Other limitations are that these tools generally do not take into account assay specificity, secondary structures, and SNPs in the primer annealing regions. Furthermore, these tools are limited to circRNA primer design for humans (no other organisms possible), and no wet-lab validation is demonstrated. Here, we present CIRCprimerXL, a circRNA RT-qPCR assay design pipeline based on the primer design framework primerXL. CIRCprimerXL takes a circRNA BSJ position as input, and designs BSJ-spanning primers using Primer3. The user can choose to use the unspliced or spliced circRNA sequence as template. Prior to primer design, sequence regions with secondary structures and common SNPs are flagged. Next, the primers are filtered based on predicted specificity and the absence of secondary structures of the amplicon to select a suitable primer pair. Our tool is both available as a user-friendly web tool and as a stand-alone pipeline based on Docker and Nextflow, allowing users to run the pipeline on a wide range of computer infrastructures. The CIRCprimerXL Nextflow pipeline can be used to design circRNA primers for any species by providing the appropriate reference genome. The CIRCprimerXL web tool supports circRNA primer design for human, mouse, rat, zebrafish, Xenopus tropicalis, and C. elegans. The design process can easily be scaled up for the qPCR assay design of tens of thousands of circRNAs within a couple of hours. We show how CIRCprimerXL has been successfully used to design qPCR assays for over 15,000 human circRNAs of which 20 were empirically validated. CIRCprimerXL software, documentation, and test data can be found at: https://github.com/OncoRNALab/CIRCprimerXL. CIRCprimerXL is also implemented as a webtool at: https://circprimerxl.cmgg.be.

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.

Comprehensive RNA dataset of tissue and plasma from patients with esophageal cancer or precursor lesions.

In the past decades, the incidence of esophageal adenocarcinoma has increased dramatically in Western populations. Better understanding of disease etiology along with the identification of novel prognostic and predictive biomarkers are urgently needed to improve the dismal survival probabilities. Here, we performed comprehensive RNA (coding and non-coding) profiling in various samples from 17 patients diagnosed with esophageal adenocarcinoma, high-grade dysplastic or non-dysplastic Barrett's esophagus. Per patient, a blood plasma sample, and a healthy and disease esophageal tissue sample were included. In total, this comprehensive dataset consists of 102 sequenced libraries from 51 samples. Based on this data, 119 expression profiles are available for three biotypes, including miRNA (51), mRNA (51) and circRNA (17). This unique resource allows for discovery of novel biomarkers and disease mechanisms, comparison of tissue and liquid biopsy profiles, integration of coding and non-coding RNA patterns, and can serve as a validation dataset in other RNA landscaping studies. Moreover, structural RNA differences can be identified in this dataset, including protein coding mutations, fusion genes, and circular RNAs.

RNA biomarkers from proximal liquid biopsy for diagnosis of ovarian cancer.

BACKGROUND: Most ovarian cancer patients are diagnosed at an advanced stage and have a high mortality rate. Current screening strategies fail to improve prognosis because markers that are sensitive for early stage disease are lacking. This medical need justifies the search for novel approaches using utero-tubal lavage as a proximal liquid biopsy. METHODS: In this study, we explore the extracellular transcriptome of utero-tubal lavage fluid obtained from 26 ovarian cancer patients and 48 controls using messenger RNA (mRNA) capture and small RNA sequencing. RESULTS: We observed an enrichment of ovarian and fallopian tube specific messenger RNAs in utero-tubal lavage fluid compared to other human biofluids. Over 300 mRNAs and 41 miRNAs were upregulated in ovarian cancer samples compared with controls. Upregulated genes were enriched for genes involved in cell cycle activation and proliferation, hinting at a tumor-derived signal. CONCLUSION: This is a proof-of-principle that mRNA capture sequencing of utero-tubal lavage fluid is technically feasible, and that the extracellular transcriptome of utero-tubal lavage should be further explored in larger cohorts to assess the diagnostic value of the biomarkers identified in this study. IMPACT: Proximal liquid biopsy from the gynecologic tract is a promising source for mRNA and miRNA biomarkers for diagnosis of early-stage ovarian cancer.

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.

The MicroRNA Landscape of Acute Beta Cell Destruction in Type 1 Diabetic Recipients of Intraportal Islet Grafts.

Ongoing beta cell death in type 1 diabetes (T1D) can be detected using biomarkers selectively discharged by dying beta cells into plasma. microRNA-375 (miR-375) ranks among the top biomarkers based on studies in animal models and human islet transplantation. Our objective was to identify additional microRNAs that are co-released with miR-375 proportionate to the amount of beta cell destruction. RT-PCR profiling of 733 microRNAs in a discovery cohort of T1D patients 1 h before/after islet transplantation indicated increased plasma levels of 22 microRNAs. Sub-selection for beta cell selectivity resulted in 15 microRNAs that were subjected to double-blinded multicenter analysis. This led to the identification of eight microRNAs that were consistently increased during early graft destruction: besides miR-375, these included miR-132/204/410/200a/429/125b, microRNAs with known function and enrichment in beta cells. Their potential clinical translation was investigated in a third independent cohort of 46 transplant patients by correlating post-transplant microRNA levels to C-peptide levels 2 months later. Only miR-375 and miR-132 had prognostic potential for graft outcome, and none of the newly identified microRNAs outperformed miR-375 in multiple regression. In conclusion, this study reveals multiple beta cell-enriched microRNAs that are co-released with miR-375 and can be used as complementary biomarkers of beta cell death.

The RNA Atlas expands the catalog of human non-coding RNAs.

Existing compendia of non-coding RNA (ncRNA) are incomplete, in part because they are derived almost exclusively from small and polyadenylated RNAs. Here we present a more comprehensive atlas of the human transcriptome, which includes small and polyA RNA as well as total RNA from 300 human tissues and cell lines. We report thousands of previously uncharacterized RNAs, increasing the number of documented ncRNAs by approximately 8%. To infer functional regulation by known and newly characterized ncRNAs, we exploited pre-mRNA abundance estimates from total RNA sequencing, revealing 316 microRNAs and 3,310 long non-coding RNAs with multiple lines of evidence for roles in regulating protein-coding genes and pathways. Our study both refines and expands the current catalog of human ncRNAs and their regulatory interactions. All data, analyses and results are available for download and interrogation in the R2 web portal, serving as a basis for future exploration of RNA biology and function.

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.

Performance assessment of total RNA sequencing of human biofluids and extracellular vesicles.

RNA profiling has emerged as a powerful tool to investigate the biomarker potential of human biofluids. However, despite enormous interest in extracellular nucleic acids, RNA sequencing methods to quantify the total RNA content outside cells are rare. Here, we evaluate the performance of the SMARTer Stranded Total RNA-Seq method in human platelet-rich plasma, platelet-free plasma, urine, conditioned medium, and extracellular vesicles (EVs) from these biofluids. We found the method to be accurate, precise, compatible with low-input volumes and able to quantify a few thousand genes. We picked up distinct classes of RNA molecules, including mRNA, lncRNA, circRNA, miscRNA and pseudogenes. Notably, the read distribution and gene content drastically differ among biofluids. In conclusion, we are the first to show that the SMARTer method can be used for unbiased unraveling of the complete transcriptome of a wide range of biofluids and their extracellular vesicles.