Search

In addition to morphologic analysis, molecular diagnostic work up of Spitz tumours is often of great value for their accurate diagnosis/classification. Nowadays, next-generation sequencing (NGS) is the predominant screening method in molecular diagnostics. Up to 80% of these melanocytic neoplasms comprise gene fusions as genetic anomalies for which the driver codes for a protein harbouring a kinase domain. However, because of the variety of fusion partners the use of PCR-based targeted enrichment NGS methods is not recommended. We describe a series of four Spitz tumour samples in which distinct gene fusions were detected by hybridisation-based capture NGS (TPM3::ALK, LIMA1::ROS1, LRRFIP2::ROS1 and MYO5A::RET). Two of these fusions are not previously described. All 4 fusions were confirmed by reverse transcription-PCR. These findings demonstrate the need for molecular analysis that can detect unknown fusions in Spitz neoplasms for optimal diagnosis., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

We performed long-read transcriptome and proteome profiling of pathogen-stimulated peripheral blood mononuclear cells (PBMCs) from healthy donors to discover new transcript and protein isoforms expressed during immune responses to diverse pathogens. Long-read transcriptome profiling reveals novel sequences and isoform switching induced upon pathogen stimulation, including transcripts that are difficult to detect using traditional short-read sequencing. Widespread loss of intron retention occurs as a common result of all pathogen stimulations. We highlight novel transcripts of NFKB1 and CASP1 that may indicate novel immunological mechanisms. RNA expression differences did not result in differences in the amounts of secreted proteins. Clustering analysis of secreted proteins revealed a correlation between chemokine (receptor) expression on the RNA and protein levels in C. albicans- and poly(I:C)-stimulated PBMCs. Isoform aware long-read sequencing of pathogen-stimulated immune cells highlights the potential of these methods to identify novel transcripts, revealing a more complex transcriptome landscape than previously appreciated., Competing Interests: The authors declare no competing interests., (© 2024 The Author(s).)

Objectives: Mutational signatures (MS) are gaining traction for deriving therapeutic insights for immune checkpoint inhibition (ICI). We asked if MS attributions from comprehensive targeted sequencing assays are reliable enough for predicting ICI efficacy in non-small cell lung cancer (NSCLC)., Methods: Somatic mutations of m = 126 patients were assayed using panel-based sequencing of 523 cancer-related genes. In silico simulations of MS attributions for various panels were performed on a separate dataset of m = 101 whole genome sequenced patients. Non-synonymous mutations were deconvoluted using COSMIC v3.3 signatures and used to test a previously published machine learning classifier., Results: The ICI efficacy predictor performed poorly with an accuracy of 0.51 -0.09 +0.09 , average precision of 0.52 -0.11 +0.11 , and an area under the receiver operating characteristic curve of 0.50 -0.09 +0.10 . Theoretical arguments, experimental data, and in silico simulations pointed to false negative rates (FNR) related to panel size. A secondary effect was observed, where deconvolution of small ensembles of point mutations lead to reconstruction errors and misattributions., Conclusion: MS attributions from current targeted panel sequencing are not reliable enough to predict ICI efficacy. We suggest that, for downstream classification tasks in NSCLC, signature attributions be based on whole exome or genome sequencing instead., Competing Interests: Declaration of Competing Interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: HCD: None to declare; KC: None to declare; GF: None to declare; HJMG: Consulting or Advisory Role: Novartis, Lilly. TJNH: Advisory/consultancy fees from AstraZeneca, Bristol-Myers-Squibb, Illumina, Merck Sharp Dohme, Roche, and research grants/funding from AstraZeneca, Hoffmann-La Roche. BM: None to declare; ES: Honoraria/speakers fee: Bio-Rad, Roche, Agena Bioscience, Illumina, Lilly; Consulting or Advisory Role: MSD/Merck, Astellas, Bayer, BMS, Agena Bioscience, Janssen Cilag (Johnson & Johnson), Novartis, Roche, AstraZeneca, Amgen, Lilly; Research Funding: Biocartis, Bio-Rad, Roche, Agena Bioscience, AstraZeneca, InVitae/Archer (all paid to UMCG); Travel, Accommodations, Expenses: Roche Molecular Diagnostics, Bio-Rad. PJV: None to declare; GAL: Research grant from Boehringer-Ingelheim. Owns shares in Genomics PLC.; BvE: None to declare., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

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., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Our incomplete knowledge of the human transcriptome impairs the detection of disease-causing variants, in particular if they affect transcripts only expressed under certain conditions. These transcripts are often lacking from reference transcript sets, such as Ensembl/GENCODE and RefSeq, and could be relevant for establishing genetic diagnoses. We present SUsPECT (Solving Unsolved Patient Exomes/gEnomes using Custom Transcriptomes), a pipeline based on the Ensembl Variant Effect Predictor (VEP) to predict variant impact on custom transcript sets, such as those generated by long-read RNA-sequencing, for downstream prioritization. Our pipeline predicts the functional consequence and likely deleteriousness scores for missense variants in the context of novel open reading frames predicted from any transcriptome. We demonstrate the utility of SUsPECT by uncovering potential mutational mechanisms of pathogenic variants in ClinVar that are not predicted to be pathogenic using the reference transcript annotation. In further support of SUsPECT's utility, we identified an enrichment of immune-related variants predicted to have a more severe molecular consequence when annotating with a newly generated transcriptome from stimulated immune cells instead of the reference transcriptome. Our pipeline outputs crucial information for further prioritization of potentially disease-causing variants for any disease and will become increasingly useful as more long-read RNA sequencing datasets become available., (© 2023. The Author(s).)

Introduction: Next-generation sequencing applications are becoming indispensable for clinical diagnostics. These experiments require numerous wet- and dry-laboratory steps, each one increasing the probability of a sample swap or contamination. Therefore, identity confirmation at the end of the process is recommended to ensure the right data are used for each patient., Methods: We tested three commercially available, single nucleotide polymorphism (SNP)-based sample tracking kits in a diagnostic workflow to evaluate their ease of use and performance. The coverage uniformity, on-target specificity, sample identification, and genotyping performance were determined to assess the reliability and cost effectiveness of each kit., Results and Discussion: Hands-on time and manual steps are almost identical for the kits from pxlence and Nimagen. The Swift kit has an extra purification step, making it the longest and most demanding protocol. Furthermore, the Swift kit failed to correctly genotype 26 of the 46 samples. The Nimagen kit identified all but one sample and the pxlence kit unambiguously identified all samples, making it the most reliable and robust kit of this evaluation. The Nimagen kit showed poor on-target mapping rates, resulting in deeper sequencing needs and higher sequencing costs compared with the other two kits., Conclusion: Our conclusion is that the Human Sample ID kit from pxlence is the most cost effective of the three tested tools for DNA sample tracking and identification., (© 2022. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

Maintaining high sensitivity while limiting false positives is a key challenge in peptide identification from mass spectrometry data. Here, we investigate the effects of integrating the machine learning-based postprocessor Percolator into our spectral library searching tool COSS (CompOmics Spectral library Searching tool). To evaluate the effects of this postprocessing, we have used 40 data sets from 2 different projects and have searched these against the NIST and MassIVE spectral libraries. The searching is carried out using 2 spectral library search tools, COSS and MSPepSearch with and without Percolator postprocessing, and using sequence database search engine MS-GF+ as a baseline comparator. The addition of the Percolator rescoring step to COSS is effective and results in a substantial improvement in sensitivity and specificity of the identifications. COSS is freely available as open source under the permissive Apache2 license, and binaries and source code are found at https://github.com/compomics/COSS.

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Vromman, Anckaert, Vandesompele and Volders.)

Accumulating evidence highlights the role of long non-coding RNAs (lncRNAs) in cellular homeostasis, and their dysregulation in disease settings. Most lncRNAs function by interacting with proteins or protein complexes. While several orthogonal methods have been developed to identify these proteins, each method has its inherent strengths and limitations. Here, we combine two RNA-centric methods ChIRP-MS and RNA-BioID to obtain a comprehensive list of proteins that interact with the well-known lncRNA HOTAIR. Overexpression of HOTAIR has been associated with a metastasis-promoting phenotype in various cancers. Although HOTAIR is known to bind with PRC2 and LSD1 protein complexes, only very limited unbiased comprehensive approaches to map its interactome have been performed. Both ChIRP-MS and RNA-BioID data sets show an association of HOTAIR with mitoribosomes, suggesting that HOTAIR has functions independent of its (post-)transcriptional mode-of-action., (© 2022. The Author(s).)

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., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Roughly half of all high-risk neuroblastoma patients present with MYCN amplification. The molecular consequences of MYCN overexpression in this aggressive pediatric tumor have been studied for decades, but thus far, our understanding of the early initiating steps of MYCN-driven tumor formation is still enigmatic. We performed a detailed transcriptome landscaping during murine TH-MYCN-driven neuroblastoma tumor formation at different time points. The neuroblastoma dependency factor MEIS2, together with ASCL1, was identified as a candidate tumor-initiating factor and shown to be a novel core regulatory circuit member in adrenergic neuroblastomas. Of further interest, we found a KEOPS complex member ( gm6890 ), implicated in homologous double-strand break repair and telomere maintenance, to be strongly upregulated during tumor formation, as well as the checkpoint adaptor Claspin ( CLSPN ) and three chromosome 17q loci CBX2 , GJC1 and LIMD2 . Finally, cross-species master regulator analysis identified FOXM1, together with additional hubs controlling transcriptome profiles of MYCN-driven neuroblastoma. In conclusion, time-resolved transcriptome analysis of early hyperplastic lesions and full-blown MYCN-driven neuroblastomas yielded novel components implicated in both tumor initiation and maintenance, providing putative novel drug targets for MYCN-driven neuroblastoma.

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., (© 2021 Wiley Periodicals LLC.)

Discovery of variant peptides such as a single amino acid variant (SAAV) in shotgun proteomics data is essential for personalized proteomics. Both the resolution of shotgun proteomics methods and the search engines have improved dramatically, allowing for confident identification of SAAV peptides. However, it is not yet known if these methods are truly successful in accurately identifying SAAV peptides without prior genomic information in the search database. We studied this in unprecedented detail by exploiting publicly available long-read RNA sequences and shotgun proteomics data from the gold standard reference cell line NA12878. Searching spectra from this cell line with the state-of-the-art open modification search engine ionbot against carefully curated search databases resulted in 96.7% false-positive SAAVs and an 85% lower true positive rate than searching with peptide search databases that incorporate prior genetic information. While adding genetic variants to the search database remains indispensable for correct peptide identification, inclusion of long-read RNA sequences in the search database contributes only 0.3% new peptide identifications. These findings reveal the differences in SAAV detection that result from various approaches, providing guidance to researchers studying SAAV peptides and developers of peptide spectrum identification tools.

The placenta is the interface between mother and fetus and inadequate function contributes to short and long-term ill-health. The placenta is absent from most large-scale RNA-Seq datasets. We therefore analyze long and small RNAs (~101 and 20 million reads per sample respectively) from 302 human placentas, including 94 cases of preeclampsia (PE) and 56 cases of fetal growth restriction (FGR). The placental transcriptome has the seventh lowest complexity of 50 human tissues: 271 genes account for 50% of all reads. We identify multiple circular RNAs and validate 6 of these by Sanger sequencing across the back-splice junction. Using large-scale mass spectrometry datasets, we find strong evidence of peptides produced by translation of two circular RNAs. We also identify novel piRNAs which are clustered on Chr1 and Chr14. PE and FGR are associated with multiple and overlapping differences in mRNA, lincRNA and circRNA but fewer consistent differences in small RNAs. Of the three protein coding genes differentially expressed in both PE and FGR, one encodes a secreted protein FSTL3 (follistatin-like 3). Elevated serum levels of FSTL3 in pregnant women are predictive of subsequent PE and FGR. To aid visualization of our placenta transcriptome data, we develop a web application ( https://www.obgyn.cam.ac.uk/placentome/ ).

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., (© The Author(s) 2020. Published by Oxford University Press.)

Spectral similarity searching to identify peptide-derived MS/MS spectra is a promising technique, and different spectrum similarity search tools have therefore been developed. Each of these tools, however, comes with some limitations, mainly because of low processing speed and issues with handling large databases. Furthermore, the number of spectral data formats supported is typically limited, which also creates a threshold to adoption. We have therefore developed COSS (CompOmics Spectral Searching), a new and user-friendly spectral library search tool supporting two scoring functions. COSS also includes decoy spectra generation for result validation. We have benchmarked COSS on three different spectral libraries and compared the results with established spectral searching tools and a sequence database search tool. Our comparison showed that COSS more reliably identifies spectra, is capable of handling large data sets and libraries, and is an easy to use tool that can run on low computer specifications. COSS binaries and source code can be freely downloaded from https://github.com/compomics/COSS.

γδ 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., (© 2020 The Authors.)

Background: To understand biology and differences among various tissues or cell types, one typically searches for molecular features that display characteristic abundance patterns. Several specificity metrics have been introduced to identify tissue-specific molecular features, but these either require an equal number of replicates per tissue or they can't handle replicates at all., Results: We describe a non-parametric specificity score that is compatible with unequal sample group sizes. To demonstrate its usefulness, the specificity score was calculated on all GTEx samples, detecting known and novel tissue-specific genes. A webtool was developed to browse these results for genes or tissues of interest. An example python implementation of SPECS is available at https://github.com/celineeveraert/SPECS. The precalculated SPECS results on the GTEx data are available through a user-friendly browser at specs.cmgg.be., Conclusions: SPECS is a non-parametric method that identifies known and novel specific-expressed genes. In addition, SPECS could be adopted for other features and applications.

In recent years, technological advances in transcriptome profiling revealed that the repertoire of human RNA molecules is more diverse and extended than originally thought. This diversity and complexity mainly derive from a large ensemble of noncoding RNAs. Because of their key roles in cellular processes important for normal development and physiology, disruption of noncoding RNA expression is intrinsically linked to human disease, including cancer. Therefore, studying the noncoding portion of the transcriptome offers the prospect of identifying novel therapeutic and diagnostic targets. Although evidence of the relevance of noncoding RNAs in cancer is accumulating, we still face many challenges when it comes to accurately profiling their expression levels. Some of these challenges are inherent to the technologies employed, whereas others are associated with characteristics of the noncoding RNAs themselves. In this review, we discuss the challenges related to long noncoding RNA expression profiling, highlight how cancer long noncoding RNAs provide new opportunities for cancer diagnosis and treatment, and reflect on future developments., (© 2018 Wiley Periodicals, Inc.)

Transcription initiates at both coding and noncoding genomic elements, including mRNA and long noncoding RNA (lncRNA) core promoters and enhancer RNAs (eRNAs). However, each class has a different expression profile with lncRNAs and eRNAs being the most tissue specific. How these complex differences in expression profiles and tissue specificities are encoded in a single DNA sequence remains unresolved. Here, we address this question using computational approaches and massively parallel reporter assays (MPRA) surveying hundreds of promoters and enhancers. We find that both divergent lncRNA and mRNA core promoters have higher capacities to drive transcription than nondivergent lncRNA and mRNA core promoters, respectively. Conversely, intergenic lncRNAs (lincRNAs) and eRNAs have lower capacities to drive transcription and are more tissue specific than divergent genes. This higher tissue specificity is strongly associated with having less complex transcription factor (TF) motif profiles at the core promoter. We experimentally validated these findings by testing both engineered single-nucleotide deletions and human single-nucleotide polymorphisms (SNPs) in MPRA. In both cases, we observe that single nucleotides associated with many motifs are important drivers of promoter activity. Thus, we suggest that high TF motif density serves as a robust mechanism to increase promoter activity at the expense of tissue specificity. Moreover, we find that 22% of common SNPs in core promoter regions have significant regulatory effects. Collectively, our findings show that high TF motif density provides redundancy and increases promoter activity at the expense of tissue specificity, suggesting that specificity of expression may be regulated by simplicity of motif usage., (© 2019 Mattioli et al.; Published by Cold Spring Harbor Laboratory Press.)

While long non-coding RNA (lncRNA) research in the past has primarily focused on the discovery of novel genes, today it has shifted towards functional annotation of this large class of genes. With thousands of lncRNA studies published every year, the current challenge lies in keeping track of which lncRNAs are functionally described. This is further complicated by the fact that lncRNA nomenclature is not straightforward and lncRNA annotation is scattered across different resources with their own quality metrics and definition of a lncRNA. To overcome this issue, large scale curation and annotation is needed. Here, we present the fifth release of the human lncRNA database LNCipedia (https://lncipedia.org). The most notable improvements include manual literature curation of 2482 lncRNA articles and the use of official gene symbols when available. In addition, an improved filtering pipeline results in a higher quality reference lncRNA gene set.

The 2017 Dagstuhl Seminar on Computational Proteomics provided an opportunity for a broad discussion on the current state and future directions of the generation and use of peptide tandem mass spectrometry spectral libraries. Their use in proteomics is growing slowly, but there are multiple challenges in the field that must be addressed to further increase the adoption of spectral libraries and related techniques. The primary bottlenecks are the paucity of high quality and comprehensive libraries and the general difficulty of adopting spectral library searching into existing workflows. There are several existing spectral library formats, but none captures a satisfactory level of metadata; therefore, a logical next improvement is to design a more advanced, Proteomics Standards Initiative-approved spectral library format that can encode all of the desired metadata. The group discussed a series of metadata requirements organized into three designations of completeness or quality, tentatively dubbed bronze, silver, and gold. The metadata can be organized at four different levels of granularity: at the collection (library) level, at the individual entry (peptide ion) level, at the peak (fragment ion) level, and at the peak annotation level. Strategies for encoding mass modifications in a consistent manner and the requirement for encoding high-quality and commonly seen but as-yet-unidentified spectra were discussed. The group also discussed related topics, including strategies for comparing two spectra, techniques for generating representative spectra for a library, approaches for selection of optimal signature ions for targeted workflows, and issues surrounding the merging of two or more libraries into one. We present here a review of this field and the challenges that the community must address in order to accelerate the adoption of spectral libraries in routine analysis of proteomics datasets.

The landscape of somatic copy-number alterations (SCNAs) affecting long non-coding RNAs (lncRNAs) in human cancers remains largely unexplored. While the majority of lncRNAs remain to be functionally characterized, several have been implicated in cancer development and metastasis. Considering the plethora of lncRNAs genes that have been currently reported, it is conceivable that many more lncRNAs might function as oncogenes or tumor suppressor genes. We devised a strategy to detect focal lncRNA SCNAs using a custom DNA microarray platform probing 10,519 lncRNA genes. By screening a panel of 80 cancer cell lines, we detected numerous focal aberrations targeting one or multiple lncRNAs without affecting neighboring protein-coding genes. These focal aberrations are highly suggestive for a tumor suppressive or oncogenic role of the targeted lncRNA gene. Although functional validation remains an essential step in the further characterization of the involved candidate cancer lncRNAs, our results provide a direct way of prioritizing candidate lncRNAs that are involved in cancer pathogenesis.

The combination of genome-edited human embryonic stem cells (hESCs) and subsequent neural differentiation is a powerful tool to study neurodevelopmental disorders. Since the naïve state of pluripotency has favourable characteristics for efficient genome-editing, we optimized a workflow for the CRISPR/Cas9 system in these naïve stem cells. Editing efficiencies of respectively 1.3-8.4% and 3.8-19% were generated with the Cas9 nuclease and the D10A Cas9 nickase mutant. Next to this, wildtype and genome-edited naïve hESCs were successfully differentiated to neural progenitor cells. As a proof-of-principle of our workflow, two monoclonal genome-edited naïve hESCs colonies were obtained for TUNA, a long non-coding RNA involved in pluripotency and neural differentiation. In these genome-edited hESCs, an effect was seen on expression of TUNA, although not on neural differentiation potential. In conclusion, we optimized a genome-editing workflow in naïve hESCs that can be used to study candidate genes involved in neural differentiation and/or functioning.

Over the past decade, long noncoding RNAs (lncRNAs) have emerged as novel functional entities of the eukaryotic genome. However, the scientific community remains divided over the amount of true noncoding transcripts among the large number of unannotated transcripts identified by recent large scale and deep RNA-sequencing efforts. Here, we systematically exclude possible technical reasons underlying the absence of lncRNA-encoded proteins in mass spectrometry data sets, strongly suggesting that the large majority of lncRNAs is indeed not translated.

Metastatic melanoma of the skin has a high mortality despite the recent introduction of targeted therapy and immunotherapy. Long non-coding RNAs (lncRNAs) are defined as transcripts of more than 200 nucleotides in length that lack protein-coding potential. There is growing evidence that lncRNAs play an important role in gene regulation, including oncogenesis. We present 13 lncRNA genes involved in the pathogenesis of cutaneous melanoma through a variety of pathways and molecular interactions. Some of these lncRNAs are possible biomarkers or therapeutic targets for malignant melanoma.

Background: Reconstructing transcript models from RNA-sequencing (RNA-seq) data and establishing these as independent transcriptional units can be a challenging task. Current state-of-the-art tools for long non-coding RNA (lncRNA) annotation are mainly based on evolutionary constraints, which may result in false negatives due to the overall limited conservation of lncRNAs., Results: To tackle this problem we have developed the Zipper plot, a novel visualization and analysis method that enables users to simultaneously interrogate thousands of human putative transcription start sites (TSSs) in relation to various features that are indicative for transcriptional activity. These include publicly available CAGE-sequencing, ChIP-sequencing and DNase-sequencing datasets. Our method only requires three tab-separated fields (chromosome, genomic coordinate of the TSS and strand) as input and generates a report that includes a detailed summary table, a Zipper plot and several statistics derived from this plot., Conclusion: Using the Zipper plot, we found evidence of transcription for a set of well-characterized lncRNAs and observed that fewer mono-exonic lncRNAs have CAGE peaks overlapping with their TSSs compared to multi-exonic lncRNAs. Using publicly available RNA-seq data, we found more than one hundred cases where junction reads connected protein-coding gene exons with a downstream mono-exonic lncRNA, revealing the need for a careful evaluation of lncRNA 5'-boundaries. Our method is implemented using the statistical programming language R and is freely available as a webtool.

In microRNA (miRNA) target prediction, typically two levels of information need to be modeled: the number of potential miRNA binding sites present in a target mRNA and the genomic context of each individual site. Single model structures insufficiently cope with this complex training data structure, consisting of feature vectors of unequal length as a consequence of the varying number of miRNA binding sites in different mRNAs. To circumvent this problem, we developed a two-layered, stacked model, in which the influence of binding site context is separately modeled. Using logistic regression and random forests, we applied the stacked model approach to a unique data set of 7990 probed miRNA-mRNA interactions, hereby including the largest number of miRNAs in model training to date. Compared to lower-complexity models, a particular stacked model, named miSTAR (miRNA stacked model target prediction; www.mi-star.org), displays a higher general performance and precision on top scoring predictions. More importantly, our model outperforms published and widely used miRNA target prediction algorithms. Finally, we highlight flaws in cross-validation schemes for evaluation of miRNA target prediction models and adopt a more fair and stringent approach., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

RNAcentral is a database of non-coding RNA (ncRNA) sequences that aggregates data from specialised ncRNA resources and provides a single entry point for accessing ncRNA sequences of all ncRNA types from all organisms. Since its launch in 2014, RNAcentral has integrated twelve new resources, taking the total number of collaborating database to 22, and began importing new types of data, such as modified nucleotides from MODOMICS and PDB. We created new species-specific identifiers that refer to unique RNA sequences within a context of single species. The website has been subject to continuous improvements focusing on text and sequence similarity searches as well as genome browsing functionality. All RNAcentral data is provided for free and is available for browsing, bulk downloads, and programmatic access at http://rnacentral.org/., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Database Url: http://www.decoderna.org., (© The Author(s) 2017. Published by Oxford University Press.)

Long non-coding RNAs (lncRNAs) form a new class of RNA molecules implicated in various aspects of protein coding gene expression regulation. To study lncRNAs in cancer, we generated expression profiles for 1707 human lncRNAs in the NCI60 cancer cell line panel using a high-throughput nanowell RT-qPCR platform. We describe how qPCR assays were designed and validated and provide processed and normalized expression data for further analysis. Data quality is demonstrated by matching the lncRNA expression profiles with phenotypic and genomic characteristics of the cancer cell lines. This data set can be integrated with publicly available omics and pharmacological data sets to uncover novel associations between lncRNA expression and mRNA expression, miRNA expression, DNA copy number, protein coding gene mutation status or drug response.

Recently, exome sequencing led to the identification of causal mutations in 16-31% of patients with intellectual disability (ID), leaving the underlying cause for many patients unidentified. In this context, the noncoding part of the human genome remains largely unexplored. For many long non-coding RNAs (lncRNAs) a crucial role in neurodevelopment and hence the human brain is anticipated. Here we aimed at identifying lncRNAs associated with neuronal development and ID. Therefore, we applied an integrated genomics approach, harnessing several public epigenetic datasets. We found that the presence of neuron-specific H3K4me3 confers the highest specificity for genes involved in neurodevelopment and ID. Based on the presence of this feature and GWAS hits for CNS disorders, we identified 53 candidate lncRNA genes. Extensive expression profiling on human brain samples and other tissues, followed by Gene Set Enrichment Analysis indicates that at least 24 of these lncRNAs are indeed implicated in processes such as synaptic transmission, nervous system development and neurogenesis. The bidirectional or antisense overlapping orientation relative to multiple coding genes involved in neuronal processes supports these results. In conclusion, we identified several lncRNA genes putatively involved in neurodevelopment and CNS disorders, providing a resource for functional studies.

Restoration of the antitumor activity of p53 could offer a promising approach for the treatment of neuroblastoma. MicroRNAs (miRNAs) are important mediators of p53 activity, but their role in the p53 response has not yet been comprehensively addressed in neuroblastoma. Therefore, we set out to characterize alterations in miRNA expression that are induced by p53 activation in neuroblastoma cells. Genome-wide miRNA expression analysis showed that miR-34a-5p, miR-182-5p, miR-203a, miR-222-3p, and miR-432-5p are upregulated following nutlin-3 treatment in a p53 dependent manner. The function of miR-182-5p, miR-203a, miR-222-3p, and miR-432-5p was analyzed by ectopic overexpression of miRNA mimics. We observed that these p53-regulated miRNAs inhibit the proliferation of neuroblastoma cells to varying degrees, with the most profound growth inhibition recorded for miR-182-5p. Overexpression of miR-182-5p promoted apoptosis in some neuroblastoma cell lines and induced neuronal differentiation of NGP cells. Using Chromatin Immunoprecipitation-qPCR (ChIP-qPCR), we did not observe direct binding of p53 to MIR182, MIR203, MIR222, and MIR432 in neuroblastoma cells. Taken together, our findings yield new insights in the network of p53-regulated miRNAs in neuroblastoma.

Structural genomic variations play an important role in human disease and phenotypic diversity. With the rise of high-throughput sequencing tools, mate-pair/paired-end/single-read sequencing has become an important technique for the detection and exploration of structural variation. Several analysis tools exist to handle different parts and aspects of such sequencing based structural variation analyses pipelines. A comprehensive analysis platform to handle all steps, from processing the sequencing data, to the discovery and visualization of structural variants, is missing. The ViVar platform is built to handle the discovery of structural variants, from Depth Of Coverage analysis, aberrant read pair clustering to split read analysis. ViVar provides you with powerful visualization options, enables easy reporting of results and better usability and data management. The platform facilitates the processing, analysis and visualization, of structural variation based on massive parallel sequencing data, enabling the rapid identification of disease loci or genes. ViVar allows you to scale your analysis with your work load over multiple (cloud) servers, has user access control to keep your data safe and is easy expandable as analysis techniques advance. URL: https://www.cmgg.be/vivar/

Genetic studies in T-cell acute lymphoblastic leukemia have uncovered a remarkable complexity of oncogenic and loss-of-function mutations. Amongst this plethora of genetic changes, NOTCH1 activating mutations stand out as the most frequently occurring genetic defect, identified in more than 50% of T-cell acute lymphoblastic leukemias, supporting a role as an essential driver for this gene in T-cell acute lymphoblastic leukemia oncogenesis. In this study, we aimed to establish a comprehensive compendium of the long non-coding RNA transcriptome under control of Notch signaling. For this purpose, we measured the transcriptional response of all protein coding genes and long non-coding RNAs upon pharmacological Notch inhibition in the human T-cell acute lymphoblastic leukemia cell line CUTLL1 using RNA-sequencing. Similar Notch dependent profiles were established for normal human CD34(+) thymic T-cell progenitors exposed to Notch signaling activity in vivo. In addition, we generated long non-coding RNA expression profiles (array data) from ex vivo isolated Notch active CD34(+) and Notch inactive CD4(+)CD8(+) thymocytes and from a primary cohort of 15 T-cell acute lymphoblastic leukemia patients with known NOTCH1 mutation status. Integration of these expression datasets with publicly available Notch1 ChIP-sequencing data resulted in the identification of long non-coding RNAs directly regulated by Notch activity in normal and malignant T cells. Given the central role of Notch in T-cell acute lymphoblastic leukemia oncogenesis, these data pave the way for the development of novel therapeutic strategies that target hyperactive Notch signaling in human T-cell acute lymphoblastic leukemia., (Copyright© Ferrata Storti Foundation.)

Since 2002, information on individual microRNAs (miRNAs), such as reference names and sequences, has been stored in miRBase, the reference database for miRNA annotation. As a result of progressive insights into the miRNome and its complexity, miRBase underwent addition and deletion of miRNA records, changes in annotated miRNA sequences and adoption of more complex naming schemes over time. Unfortunately, miRBase does not allow straightforward assessment of these ongoing miRNA annotation changes, which has resulted in substantial ambiguity regarding miRNA identity and sequence in public literature, in target prediction databases and in content on various commercially available analytical platforms. As a result, correct interpretation, comparison and integration of miRNA study results are compromised, which we demonstrate here by assessing the impact of ignoring sequence annotation changes. To address this problem, we developed miRBase Tracker (www.mirbasetracker.org), an easy-to-use online database that keeps track of all historical and current miRNA annotation present in the miRBase database. Three basic functionalities allow researchers to keep their miRNA annotation up-to-date, reannotate analytical miRNA platforms and link published results with outdated annotation to the latest miRBase release. We expect miRBase Tracker to increase the transparency and annotation accuracy in the field of miRNA research., Database Url: www.mirbasetracker.org., (© The Author(s) 2014. Published by Oxford University Press.)

Here, we present LNCipedia (http://www.lncipedia.org), a novel database for human long non-coding RNA (lncRNA) transcripts and genes. LncRNAs constitute a large and diverse class of non-coding RNA genes. Although several lncRNAs have been functionally annotated, the majority remains to be characterized. Different high-throughput methods to identify new lncRNAs (including RNA sequencing and annotation of chromatin-state maps) have been applied in various studies resulting in multiple unrelated lncRNA data sets. LNCipedia offers 21 488 annotated human lncRNA transcripts obtained from different sources. In addition to basic transcript information and gene structure, several statistics are determined for each entry in the database, such as secondary structure information, protein coding potential and microRNA binding sites. Our analyses suggest that, much like microRNAs, many lncRNAs have a significant secondary structure, in-line with their presumed association with proteins or protein complexes. Available literature on specific lncRNAs is linked, and users or authors can submit articles through a web interface. Protein coding potential is assessed by two different prediction algorithms: Coding Potential Calculator and HMMER. In addition, a novel strategy has been integrated for detecting potentially coding lncRNAs by automatically re-analysing the large body of publicly available mass spectrometry data in the PRIDE database. LNCipedia is publicly available and allows users to query and download lncRNA sequences and structures based on different search criteria. The database may serve as a resource to initiate small- and large-scale lncRNA studies. As an example, the LNCipedia content was used to develop a custom microarray for expression profiling of all available lncRNAs.