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.

Benchmarking digital PCR partition classification methods with empirical and simulated duplex data.

Digital PCR (dPCR) is a highly accurate technique for the quantification of target nucleic acid(s). It has shown great potential in clinical applications, like tumor liquid biopsy and validation of biomarkers. Accurate classification of partitions based on end-point fluorescence intensities is crucial to avoid biased estimators of the concentration of the target molecules. We have evaluated many clustering methods, from general-purpose methods to specific methods for dPCR and flowcytometry, on both simulated and real-life data. Clustering method performance was evaluated by simulating various scenarios. Based on our extensive comparison of clustering methods, we describe the limits of these methods, and formulate guidelines for choosing an appropriate method. In addition, we have developed a novel method for simulating realistic dPCR data. The method is based on a mixture distribution of a Poisson point process and a skew-$t$ distribution, which enables the generation of irregularities of cluster shapes and randomness of partitions between clusters ('rain') as commonly observed in dPCR data. Users can fine-tune the model parameters and generate labeled datasets, using their own data as a template. Besides, the database of experimental dPCR data augmented with the labeled simulated data can serve as training and testing data for new clustering methods. The simulation method is available as an R Shiny app.

Artifacts and biases of the reverse transcription reaction in RNA sequencing.

RNA sequencing has spurred a significant number of research areas in recent years. Most protocols rely on synthesizing a more stable complementary DNA (cDNA) copy of the RNA molecule during the reverse transcription reaction. The resulting cDNA pool is often wrongfully assumed to be quantitatively and molecularly similar to the original RNA input. Sadly, biases and artifacts confound the resulting cDNA mixture. These issues are often overlooked or ignored in the literature by those that rely on the reverse transcription process. In this review, we confront the reader with intra- and intersample biases and artifacts caused by the reverse transcription reaction during RNA sequencing experiments. To fight the reader's despair, we also provide solutions to most issues and inform on good RNA sequencing practices. We hope the reader can use this review to their advantage, thereby contributing to scientifically sound RNA studies.

Nucleic acid degradation as barrier to gene delivery: a guide to understand and overcome nuclease activity.

Gene therapy is on its way to revolutionize the treatment of both inherited and acquired diseases, by transferring nucleic acids to correct a disease-causing gene in the target cells of patients. In the fight against infectious diseases, mRNA-based therapeutics have proven to be a viable strategy in the recent Covid-19 pandemic. Although a growing number of gene therapies have been approved, the success rate is limited when compared to the large number of preclinical and clinical trials that have been/are being performed. In this review, we highlight some of the hurdles which gene therapies encounter after administration into the human body, with a focus on nucleic acid degradation by nucleases that are extremely abundant in mammalian organs, biological fluids as well as in subcellular compartments. We overview the available strategies to reduce the biodegradation of gene therapeutics after administration, including chemical modifications of the nucleic acids, encapsulation into vectors and co-administration with nuclease inhibitors and discuss which strategies are applied for clinically approved nucleic acid therapeutics. In the final part, we discuss the currently available methods and techniques to qualify and quantify the integrity of nucleic acids, with their own strengths and limitations.

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.

Nationwide quality assurance of high-throughput diagnostic molecular testing during the SARS-CoV-2 pandemic: role of the Belgian National Reference Centre.

Since the onset of the coronavirus disease (COVID-19) pandemic in Belgium, UZ/KU Leuven has played a crucial role as the National Reference Centre (NRC) for respiratory pathogens, to be the first Belgian laboratory to develop and implement laboratory developed diagnostic assays for SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) and later to assess the quality of commercial kits. To meet the growing demand for decentralised testing, both clinical laboratories and government-supported high-throughput platforms were gradually deployed across Belgium. Consequently, the role of the NRC transitioned from a specialised testing laboratory to strengthening capacity and coordinating quality assurance. Here, we outline the measures taken by the NRC, the national public health institute Sciensano and the executing clinical laboratories to ensure effective quality management of molecular testing throughout the initial two years of the pandemic (March 2020 to March 2022).

RNA Extraction Method Impacts Quality Metrics and Sequencing Results in Formalin-Fixed, Paraffin-Embedded Tissue Samples.

Archived formalin-fixed, paraffin-embedded (FFPE) tissue samples are being increasingly used in molecular cancer research. Compared with fresh-frozen tissue, the nucleic acid analysis of FFPE tissue is technically more challenging. This study aimed to compare the impact of 3 different RNA extraction methods on yield, quality, and sequencing-based gene expression results in FFPE samples. RNA extraction was performed in 16 FFPE tumor specimens from patients with diffuse large B-cell lymphoma and in reference FFPE material from microsatellite-stable and microsatellite-instable cell lines (3 replicates each) using 2 silica-based procedures (A, miRNeasy FFPE; C, iCatcher FFPE Tissue RNA) and 1 isotachophoresis-based procedure (B, Ionic FFPE to Pure RNA). The RNA yield; RNA integrity, as reflected by the distribution value 200; and RNA purity, as reflected by the 260/280 and the 260/230 nm absorbance ratios, were determined. The RNA was sequenced on the NovaSeq 6000 instrument using the TruSeq RNA Exome and SMARTer Stranded Total RNA-Seq Pico v3 library preparations kits. Our results highlight the impact of RNA extraction methodology on both preanalytical and sequencing-based gene expression results. Overall, methods B and C outperformed method A because these showed significantly higher fractions of uniquely mapped reads, an increased number of detectable genes, a lower fraction of duplicated reads, and better representation of the B-cell receptor repertoire. Differences among the extraction methods were generally more explicit for the total RNA sequencing method than for the exome-capture sequencing method. Importantly, the predicative value of quality metrics varies among extraction kits, and caution should be applied when comparing and interpreting results obtained using different methods.

Blocking Abundant RNA Transcripts by High-Affinity Oligonucleotides during Transcriptome Library Preparation.

BACKGROUND: RNA sequencing has become the gold standard for transcriptome analysis but has an inherent limitation of challenging quantification of low-abundant transcripts. In contrast to microarray technology, RNA sequencing reads are proportionally divided in function of transcript abundance. Therefore, low-abundant RNAs compete against highly abundant - and sometimes non-informative - RNA species. RESULTS: We developed an easy-to-use strategy based on high-affinity RNA-binding oligonucleotides to block reverse transcription and PCR amplification of specific RNA transcripts, thereby substantially reducing their abundance in the final sequencing library. To demonstrate the broad application potential of our method, we applied it to different transcripts and library preparation strategies, including YRNAs in small RNA sequencing of human blood plasma, mitochondrial rRNAs in both 3' end sequencing and long-read sequencing, and MALAT1 in single-cell 3' end sequencing. We demonstrate that the blocking strategy is highly efficient, reproducible, specific, and generally results in better transcriptome coverage and complexity. CONCLUSION: Our method does not require modifications of the library preparation procedure apart from simply adding blocking oligonucleotides to the RT reaction and can thus be easily integrated into virtually any RNA sequencing library preparation protocol.

Closing the circle: current state and perspectives of circular RNA databases.

Circular RNAs (circRNAs) are covalently closed RNA molecules that have been linked to various diseases, including cancer. However, a precise function and working mechanism are lacking for the larger majority. Following many different experimental and computational approaches to identify circRNAs, multiple circRNA databases were developed as well. Unfortunately, there are several major issues with the current circRNA databases, which substantially hamper progression in the field. First, as the overlap in content is limited, a true reference set of circRNAs is lacking. This results from the low abundance and highly specific expression of circRNAs, and varying sequencing methods, data-analysis pipelines, and circRNA detection tools. A second major issue is the use of ambiguous nomenclature. Thus, redundant or even conflicting names for circRNAs across different databases contribute to the reproducibility crisis. Third, circRNA databases, in essence, rely on the position of the circRNA back-splice junction, whereas alternative splicing could result in circRNAs with different length and sequence. To uniquely identify a circRNA molecule, the full circular sequence is required. Fourth, circRNA databases annotate circRNAs' microRNA binding and protein-coding potential, but these annotations are generally based on presumed circRNA sequences. Finally, several databases are not regularly updated, contain incomplete data or suffer from connectivity issues. In this review, we present a comprehensive overview of the current circRNA databases and their content, features, and usability. In addition to discussing the current issues regarding circRNA databases, we come with important suggestions to streamline further research in this growing field.

Digital PCR Partition Classification.

BACKGROUND: Partition classification is a critical step in the digital PCR data analysis pipeline. A range of partition classification methods have been developed, many motivated by specific experimental setups. An overview of these partition classification methods is lacking and their comparative properties are often unclear, likely impacting the proper application of these methods. CONTENT: This review provides a summary of all available digital PCR partition classification approaches and the challenges they aim to overcome, serving as a guide for the digital PCR practitioner wishing to apply them. We additionally discuss strengths and weaknesses of these methods, which can further guide practitioners in vigilant application of these existing methods. This review provides method developers with ideas for improving methods or designing new ones. The latter is further stimulated by our identification and discussion of application gaps in the literature, for which there are currently no or few methods available. SUMMARY: This review provides an overview of digital PCR partition classification methods, their properties, and potential applications. Ideas for further advances are presented and may bolster method development.

Digital PCR-based evaluation of nucleic acid extraction kit performance for the co-purification of cell-free DNA and RNA.

BACKGROUND: Blood plasma, one of the most studied liquid biopsies, contains various molecules that have biomarker potential for cancer detection, including cell-free DNA (cfDNA) and cell-free RNA (cfRNA). As the vast majority of cell-free nucleic acids in circulation are non-cancerous, a laboratory workflow with a high detection sensitivity of tumor-derived nucleic acids is a prerequisite for precision oncology. One way to meet this requirement is by the combined analysis of cfDNA and cfRNA from the same liquid biopsy sample. So far, no study has systematically compared the performance of cfDNA and cfRNA co-purification to increase sensitivity. RESULTS: First, we set up a framework using digital PCR (dPCR) technology to quantify cfDNA and cfRNA from human blood plasma in order to compare cfDNA/cfRNA co-purification kit performance. To that end, we optimized two dPCR duplex assays, designed to quantify both cfDNA and cfRNA with the same assays, by ensuring that primers and probes are located within a highly abundant exon. Next, we applied our optimized workflow to evaluate the co-purification performance of two manual and two semi-automated methods over a range of plasma input volumes (0.06-4 mL). Some kits result in higher nucleic acid concentrations in the eluate, while consuming only half of the plasma volume. The combined nucleic acid quantification systematically results in higher nucleic acid concentrations as compared to a parallel quantification of cfDNA and cfRNA in the eluate. CONCLUSIONS: We provide a framework to evaluate the performance of cfDNA/cfRNA co-purification kits and have tested two manual and two semi-automated co-purification kits in function of the available plasma input amount and the intended use of the nucleic acid eluate. We demonstrate that the combined quantification of cfDNA and cfRNA has a benefit compared to separate quantification. We foresee that the results of this study are instrumental for clinical applications to help increase mutation detection sensitivity, allowing improved disease detection and monitoring.

Transient nuclear lamin A/C accretion aids in recovery from vapor nanobubble-induced permeabilisation of the plasma membrane.

Vapor nanobubble (VNB) photoporation is a physical method for intracellular delivery that has gained significant interest in the past decade. It has successfully been used to introduce molecular cargo of diverse nature into different cell types with high throughput and minimal cytotoxicity. For translational purposes, it is important to understand whether and how photoporation affects cell homeostasis. To obtain a comprehensive view on the transcriptional rewiring that takes place after VNB photoporation, we performed a longitudinal shotgun RNA-sequencing experiment. Six hours after photoporation, we found a marked upregulation of LMNA transcripts as well as their protein products, the A-type lamins. At the same time point, we observed a significant increase in several heterochromatin marks, suggesting a global stiffening of the nucleus. These molecular features vanished 24 h after photoporation. Since VNB-induced chromatin condensation was prolonged in LMNA knockout cells, A-type lamins may be required for restoring the nucleus to its original state. Selective depletion of A-type lamins reduced cell viability after VNB photoporation, while pharmacological stimulation of LMNA transcription increased the percentage of successfully transfected cells that survived after photoporation. Therefore, our results suggest that cells respond to VNB photoporation by temporary upregulation of A-type lamins to facilitate their recovery.

Distinct Notch1 and BCL11B requirements mediate human γδ/αß T cell development.

γδ and αß T cells have unique roles in immunity and both originate in the thymus from T-lineage committed precursors through distinct but unclear mechanisms. Here, we show that Notch1 activation is more stringently required for human γδ development compared to αß-lineage differentiation and performed paired mRNA and miRNA profiling across 11 discrete developmental stages of human T cell development in an effort to identify the potential Notch1 downstream mechanism. Our data suggest that the miR-17-92 cluster is a Notch1 target in immature thymocytes and that miR-17 can restrict BCL11B expression in these Notch-dependent T cell precursors. We show that enforced miR-17 expression promotes human γδ T cell development and, consistently, that BCL11B is absolutely required for αß but less for γδ T cell development. This study suggests that human γδ T cell development is mediated by a stage-specific Notch-driven negative feedback loop through which miR-17 temporally restricts BCL11B expression and provides functional insights into the developmental role of the disease-associated genes BCL11B and the miR-17-92 cluster in a human context.

Melanoma addiction to the long non-coding RNA SAMMSON.

Focal amplifications of chromosome 3p13-3p14 occur in about 10% of melanomas and are associated with a poor prognosis. The melanoma-specific oncogene MITF resides at the epicentre of this amplicon. However, whether other loci present in this amplicon also contribute to melanomagenesis is unknown. Here we show that the recently annotated long non-coding RNA (lncRNA) gene SAMMSON is consistently co-gained with MITF. In addition, SAMMSON is a target of the lineage-specific transcription factor SOX10 and its expression is detectable in more than 90% of human melanomas. Whereas exogenous SAMMSON increases the clonogenic potential in trans, SAMMSON knockdown drastically decreases the viability of melanoma cells irrespective of their transcriptional cell state and BRAF, NRAS or TP53 mutational status. Moreover, SAMMSON targeting sensitizes melanoma to MAPK-targeting therapeutics both in vitro and in patient-derived xenograft models. Mechanistically, SAMMSON interacts with p32, a master regulator of mitochondrial homeostasis and metabolism, to increase its mitochondrial targeting and pro-oncogenic function. Our results indicate that silencing of the lineage addiction oncogene SAMMSON disrupts vital mitochondrial functions in a cancer-cell-specific manner; this silencing is therefore expected to deliver highly effective and tissue-restricted anti-melanoma therapeutic responses.

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.

Recurrent chromosomal imbalances provide selective advantage to human embryonic stem cells under enhanced replicative stress conditions.

Human embryonic stem cells (hESCs) and embryonal tumors share a number of common features, including a compromised G1/S checkpoint. Consequently, these rapidly dividing hESCs and cancer cells undergo elevated levels of replicative stress, inducing genomic instability that drives chromosomal imbalances. In this context, it is of interest that long-term in vitro cultured hESCs exhibit a remarkable high incidence of segmental DNA copy number gains, some of which are also highly recurrent in certain malignancies such as 17q gain (17q+). The selective advantage of DNA copy number changes in these cells has been attributed to several underlying processes including enhanced proliferation. We hypothesized that these recurrent chromosomal imbalances become rapidly embedded in the cultured hESCs through a replicative stress driven Darwinian selection process. To this end, we compared the effect of hydroxyurea-induced replicative stress vs normal growth conditions in an equally mixed cell population of isogenic euploid and 17q + hESCs. We could show that 17q + hESCs rapidly overtook normal hESCs. Our data suggest that recurrent chromosomal segmental gains provide a proliferative advantage to hESCs under increased replicative stress, a process that may also explain the highly recurrent nature of certain imbalances in cancer.

SPsimSeq: semi-parametric simulation of bulk and single-cell RNA-sequencing data.

SUMMARY: SPsimSeq is a semi-parametric simulation method to generate bulk and single-cell RNA-sequencing data. It is designed to simulate gene expression data with maximal retention of the characteristics of real data. It is reasonably flexible to accommodate a wide range of experimental scenarios, including different sample sizes, biological signals (differential expression) and confounding batch effects. AVAILABILITY AND IMPLEMENTATION: The R package and associated documentation is available from https://github.com/CenterForStatistics-UGent/SPsimSeq. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The Dangers of Using Cq to Quantify Nucleic Acid in Biological Samples: A Lesson From COVID-19.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA quantities, measured by reverse transcription quantitative PCR (RT-qPCR), have been proposed to stratify clinical risk or determine analytical performance targets. We investigated reproducibility and how setting diagnostic cutoffs altered the clinical sensitivity of coronavirus disease 2019 (COVID-19) testing. METHODS: Quantitative SARS-CoV-2 RNA distributions [quantification cycle (Cq) and copies/mL] from more than 6000 patients from 3 clinical laboratories in United Kingdom, Belgium, and the Republic of Korea were analyzed. Impact of Cq cutoffs on clinical sensitivity was assessed. The June/July 2020 INSTAND external quality assessment scheme SARS-CoV-2 materials were used to estimate laboratory reported copies/mL and to estimate the variation in copies/mL for a given Cq. RESULTS: When the WHO-suggested Cq cutoff of 25 was applied, the clinical sensitivity dropped to about 16%. Clinical sensitivity also dropped to about 27% when a simulated limit of detection of 106 copies/mL was applied. The interlaboratory variation for a given Cq value was >1000 fold in copies/mL (99% CI). CONCLUSION: While RT-qPCR has been instrumental in the response to COVID-19, we recommend Cq (cycle threshold or crossing point) values not be used to set clinical cutoffs or diagnostic performance targets due to poor interlaboratory reproducibility; calibrated copy-based units (used elsewhere in virology) offer more reproducible alternatives. We also report a phenomenon where diagnostic performance may change relative to the effective reproduction number. Our findings indicate that the disparities between patient populations across time are an important consideration when evaluating or deploying diagnostic tests. This is especially relevant to the emergency situation of an evolving pandemic.

Candidate RNA biomarkers in biofluids for early diagnosis of ovarian cancer: A systematic review.

Ovarian cancer is often diagnosed in an advanced stage and is associated with a high mortality rate. It is assumed that early detection of ovarian cancer could improve patient outcomes. Unfortunately, effective screening methods for early diagnosis of ovarian cancer are still lacking. Extracellular RNAs circulating in human biofluids can reliably be measured and are emerging as potential biomarkers in cancer. In this systematic review, we present 75 RNA biomarkers detectable in human biofluids that have been studied for early diagnosis of ovarian cancer. The majority of these markers are microRNAs identified using RT-qPCR or microarrays in blood-based fluids. A handful of studies used RNA-sequencing and explored alternative fluids, such as urine and ascites. Candidate RNA biomarkers that were more abundant in biofluids of ovarian cancer patients compared to controls in at least two independent studies include miR-21, the miR-200 family, miR-205, miR-10a and miR-346. Amongst the markers confirmed to be lower in at least two studies are miR-122, miR-193a, miR-223, miR-126 and miR-106b. While these biomarkers show promising diagnostic potential, further validation is required before implementation in routine clinical care. Challenges related to biomarker validation and reflections on future perspectives to accelerate progress in this field are discussed.

Custom long non-coding RNA capture enhances detection sensitivity in different human sample types.

Long non-coding RNAs (lncRNAs) are a heterogeneous group of transcripts that lack protein coding potential and display regulatory functions in various cellular processes. As a result of their cell- and cancer-specific expression patterns, lncRNAs have emerged as potential diagnostic and therapeutic targets. The accurate characterization of lncRNAs in bulk transcriptome data remains challenging due to their low abundance compared to protein coding genes. To tackle this issue, we describe a unique short-read custom lncRNA capture sequencing approach that relies on a comprehensive set of 565,878 capture probes for 49,372 human lncRNA genes. This custom lncRNA capture approach was evaluated on various sample types ranging from artificial high-quality RNA mixtures to more challenging formalin-fixed paraffin-embedded tissue and biofluid material. The custom enrichment approach allows the detection of a more diverse repertoire of lncRNAs, with better reproducibility and higher coverage compared to classic total RNA-sequencing.