Methods behind neoantigen prediction for personalized anticancer vaccines.

Next to conventional cancer therapies, immunotherapies such as immune checkpoint inhibitors have broadened the cancer treatment landscape over the past decades. Recent advances in next generation sequencing and bioinformatics technologies have made it possible to identify a patient's own immunogenic neoantigens. These cancer neoantigens serve as important targets for personalized immunotherapy which has the benefit of being more active and effective in targeting cancer cells. This paper is a step-by-step guide discussing the different analyses and challenges encountered during in-silico neoantigen prediction. The protocol describes all the tools and steps required for the identification of immunogenic neoantigens.

Challenges in neoantigen-directed therapeutics.

A fundamental prerequisite for the efficacy of cancer immunotherapy is the presence of functional, antigen-specific T cells within the tumor. Neoantigen-directed therapy is a promising strategy that aims at targeting the host's immune response against tumor-specific antigens, thereby eradicating cancer cells. Initial forays have been made in clinical environments utilizing vaccines and adoptive cell therapy; however, many challenges lie ahead. We provide an in-depth overview of the current state of the field with an emphasis on in silico neoantigen discovery and the clinical aspects that need to be addressed to unlock the full potential of this therapy.

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.

Performance Evaluation of Three DNA Sample Tracking Tools in a Whole Exome Sequencing Workflow.

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.

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.

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

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

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.

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.

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.

Vaccination with early ferroptotic cancer cells induces efficient antitumor immunity.

BACKGROUND: Immunotherapy represents the future of clinical cancer treatment. The type of cancer cell death determines the antitumor immune response and thereby contributes to the efficacy of anticancer therapy and long-term survival of patients. Induction of immunogenic apoptosis or necroptosis in cancer cells does activate antitumor immunity, but resistance to these cell death modalities is common. Therefore, it is of great importance to find other ways to kill tumor cells. Recently, ferroptosis has been identified as a novel, iron-dependent form of regulated cell death but whether ferroptotic cancer cells are immunogenic is unknown. METHODS: Ferroptotic cell death in murine fibrosarcoma MCA205 or glioma GL261 cells was induced by RAS-selective lethal 3 and ferroptosis was analyzed by flow cytometry, atomic force and confocal microscopy. ATP and high-mobility group box 1 (HMGB1) release were detected by luminescence and ELISA assays, respectively. Immunogenicity in vitro was analyzed by coculturing of ferroptotic cancer cells with bone-marrow derived dendritic cells (BMDCs) and rate of phagocytosis and activation/maturation of BMDCs (CD11c+CD86+, CD11c+CD40+, CD11c+MHCII+, IL-6, RNAseq analysis). The tumor prophylactic vaccination model in immune-competent and immune compromised (Rag-2-/-) mice was used to analyze ferroptosis immunogenicity. RESULTS: Ferroptosis can be induced in cancer cells by inhibition of glutathione peroxidase 4, as evidenced by confocal and atomic force microscopy and inhibitors' analysis. We demonstrate for the first time that ferroptosis is immunogenic in vitro and in vivo. Early, but not late, ferroptotic cells promote the phenotypic maturation of BMDCs and elicit a vaccination-like effect in immune-competent mice but not in Rag-2-/- mice, suggesting that the mechanism of immunogenicity is very tightly regulated by the adaptive immune system and is time dependent. Also, ATP and HMGB1, the best-characterized damage-associated molecular patterns involved in immunogenic cell death, have proven to be passively released along the timeline of ferroptosis and act as immunogenic signal associated with the immunogenicity of early ferroptotic cancer cells. CONCLUSIONS: These results pave the way for the development of new therapeutic strategies for cancers based on induction of ferroptosis, and thus broadens the current concept of immunogenic cell death and opens the door for the development of new strategies in cancer immunotherapy.

A Combined RNA Preservation and Extraction Protocol for Gene Expression Studies in Cacao Beans.

Despite the high economic importance of cacao beans, few RNA-based studies have been conducted on this plant material and hence no optimal RNA-extraction has been reported. Moreover, extraction of high-quality RNA from recalcitrant cacao bean tissue has shown many difficulties and requires optimization. Furthermore, cacao beans are mostly found at remote and under-resourced locations, which pressures the outsourcing of such analysis and thereby demands RNA-stable preservation and transportation of cacao beans. This study aims to select an appropriate RNA extraction and preservation/transportation method for cacao beans. For this purpose, three sample homogenization and five extraction protocols on cacao beans were compared. In addition, 13 preservation conditions-differing in tissue crushing degree, preservation method, duration, and temperature-were compared and evaluated. A comparative analysis revealed that CTAB-based homogenization and extraction outcompeted all tested commercial protocols in RNA yield and integrity, respectively. Preservation at -80°C affected RNA quality the least, whereas freeze-drying was most suitable for transportation at room temperature for maximum 1 week. The cacao bean RNA obtained from the selected methods were compatible for downstream applications. The results of this study will facilitate on-field sampling and transportation of genetically sensitive cacao material prior to cacao bean transcriptomic studies. In addition, valuable insights on sample homogenization, extraction, preservation, and transportation have been provided, which is of interest to every plant geneticist.

Long non-coding RNAs and latent HIV - A search for novel targets for latency reversal.

The latent cellular reservoir of HIV is recognized as the major barrier to cure from HIV infection. Long non-coding RNAs (lncRNAs) are more tissue and cell type-specific than protein coding genes, and may represent targets of choice for HIV latency reversal. Using two in vitro primary T-cell models, we identified lncRNAs dysregulated in latency. PVT1 and RP11-347C18.3 were up-regulated in common between the two models, and RP11-539L10.2 was down-regulated. The major component of the latent HIV reservoir, memory CD4+ T-cells, had higher expression of these lncRNAs, compared to naïve T-cells. Guilt-by-association analysis demonstrated that lncRNAs dysregulated in latency were associated with several cellular pathways implicated in HIV latency establishment and maintenance: proteasome, spliceosome, p53 signaling, and mammalian target of rapamycin (MTOR). PVT1, RP11-347C18.3, and RP11-539L10.2 were down-regulated by latency reversing agents, suberoylanilide hydroxamic acid and Romidepsin, suggesting that modulation of lncRNAs is a possible secondary mechanism of action of these compounds. These results will facilitate prioritization of lncRNAs for evaluation as targets for HIV latency reversal. Importantly, our study provides insights into regulatory function of lncRNA during latent HIV infection.

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

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

The development of a novel SNP genotyping assay to differentiate cacao clones.

In this study, a double-mismatch allele-specific (DMAS) qPCR SNP genotyping method has been designed, tested and validated specifically for cacao, using 65 well annotated international cacao reference accessions retrieved from the Center for Forestry Research and Technology Transfer (CEFORTT) and the International Cocoa Quarantine Centre (ICQC). In total, 42 DMAS-qPCR SNP genotyping assays have been validated, with a 98.05% overall efficiency in calling the correct genotype. In addition, the test allowed for the identification of 15.38% off-types and two duplicates, highlighting the problem of mislabeling in cacao collections and the need for conclusive genotyping assays. The developed method showed on average a high genetic diversity (He = 0.416) and information index (I = 0.601), making it applicable to assess intra-population variation. Furthermore, only the 13 most informative markers were needed to achieve maximum differentiation. This simple, effective method provides robust and accurate genotypic data which allows for more efficient resource management (e.g. tackling mislabeling, conserving valuable genetic material, parentage analysis, genetic diversity studies), thus contributing to an increased knowledge on the genetic background of cacao worldwide. Notably, the described method can easily be integrated in other laboratories for a wide range of objectives and organisms.

Antisense Oligonucleotide-Based Downregulation of the G56R Pathogenic Variant Causing NR2E3-Associated Autosomal Dominant Retinitis Pigmentosa.

The recurrent missense variant in Nuclear Receptor Subfamily 2 Group E Member 3 (NR2E3), c.166G>A, p.(Gly56Arg) or G56R, underlies 1%-2% of cases with autosomal dominant retinitis pigmentosa (adRP), a frequent, genetically heterogeneous inherited retinal disease (IRD). The mutant NR2E3 protein has a presumed dominant negative effect (DNE) by competition for dimer formation with Cone-Rod Homeobox (CRX) but with abolishment of DNA binding, acting as a repressor in trans. Both the frequency and DNE of G56R make it an interesting target for allele-specific knock-down of the mutant allele using antisense oligonucleotides (AONs), an emerging therapeutic strategy for IRD. Here, we designed gapmer AONs with or without a locked nucleic acid modification at the site of the mutation, which were analyzed for potential off-target effects. Next, we overexpressed wild type (WT) or mutant NR2E3 in RPE-1 cells, followed by AON treatment. Transcript and protein levels of WT and mutant NR2E3 were detected by reverse transcription quantitative polymerase chain reaction (RT-qPCR) and Western blot respectively. All AONs showed a general knock-down of mutant and WT NR2E3 on RNA and protein level, showing the accessibility of the region for AON-induced knockdown. Further modifications are needed however to increase allele-specificity. In conclusion, we propose the first proof-of-concept for AON-mediated silencing of a single nucleotide variation with a dominant negative effect as a therapeutic approach for NR2E3-associated adRP.

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

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

Cost-effective and robust genotyping using double-mismatch allele-specific quantitative PCR.

For a wide range of diseases, SNPs in the genome are the underlying mechanism of dysfunction. Therefore, targeted detection of these variations is of high importance for early diagnosis and (familial) screenings. While allele-specific PCR has been around for many years, its adoption for SNP genotyping or somatic mutation detection has been hampered by its low discriminating power and high costs. To tackle this, we developed a cost-effective qPCR based method, able to detect SNPs in a robust and specific manner. This study describes how to combine the basic principles of allele-specific PCR (the combination of a wild type and variant primer) with the straightforward readout of DNA-binding dye based qPCR technology. To enhance the robustness and discriminating power, an artificial mismatch in the allele-specific primer was introduced. The resulting method, called double-mismatch allele-specific qPCR (DMAS-qPCR), was successfully validated using 12 SNPs and 15 clinically relevant somatic mutations on 48 cancer cell lines. It is easy to use, does not require labeled probes and is characterized by high analytical sensitivity and specificity. DMAS-qPCR comes with a complimentary online assay design tool, available for the whole scientific community, enabling researchers to design custom assays and implement those as a diagnostic test.

Long noncoding RNA expression profiling in cancer: Challenges and opportunities.

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.

A mechanistic classification of clinical phenotypes in neuroblastoma.

Neuroblastoma is a pediatric tumor of the sympathetic nervous system. Its clinical course ranges from spontaneous tumor regression to fatal progression. To investigate the molecular features of the divergent tumor subtypes, we performed genome sequencing on 416 pretreatment neuroblastomas and assessed telomere maintenance mechanisms in 208 of these tumors. We found that patients whose tumors lacked telomere maintenance mechanisms had an excellent prognosis, whereas the prognosis of patients whose tumors harbored telomere maintenance mechanisms was substantially worse. Survival rates were lowest for neuroblastoma patients whose tumors harbored telomere maintenance mechanisms in combination with RAS and/or p53 pathway mutations. Spontaneous tumor regression occurred both in the presence and absence of these mutations in patients with telomere maintenance-negative tumors. On the basis of these data, we propose a mechanistic classification of neuroblastoma that may benefit the clinical management of patients.