AutoGVP: a dockerized workflow integrating ClinVar and InterVar germline sequence variant classification.

SUMMARY: With the increasing rates of exome and whole genome sequencing, the ability to classify large sets of germline sequencing variants using up-to-date American College of Medical Genetics-Association for Molecular Pathology (ACMG-AMP) criteria is crucial. Here, we present Automated Germline Variant Pathogenicity (AutoGVP), a tool that integrates germline variant pathogenicity annotations from ClinVar and sequence variant classifications from a modified version of InterVar (PVS1 strength adjustments, removal of PP5/BP6). This tool facilitates large-scale, clinically focused classification of germline sequence variants in a research setting. AVAILABILITY AND IMPLEMENTATION: AutoGVP is an open source dockerized workflow implemented in R and freely available on GitHub at https://github.com/diskin-lab-chop/AutoGVP.

A transcriptome-based classifier to determine molecular subtypes in medulloblastoma.

Medulloblastoma is a highly heterogeneous pediatric brain tumor with five molecular subtypes, Sonic Hedgehog TP53-mutant, Sonic Hedgehog TP53-wildtype, WNT, Group 3, and Group 4, defined by the World Health Organization. The current mechanism for classification into these molecular subtypes is through the use of immunostaining, methylation, and/or genetics. We surveyed the literature and identified a number of RNA-Seq and microarray datasets in order to develop, train, test, and validate a robust classifier to identify medulloblastoma molecular subtypes through the use of transcriptomic profiling data. We have developed a GPL-3 licensed R package and a Shiny Application to enable users to quickly and robustly classify medulloblastoma samples using transcriptomic data. The classifier utilizes a large composite microarray dataset (15 individual datasets), an individual microarray study, and an RNA-Seq dataset, using gene ratios instead of gene expression measures as features for the model. Discriminating features were identified using the limma R package and samples were classified using an unweighted mean of normalized scores. We utilized two training datasets and applied the classifier in 15 separate datasets. We observed a minimum accuracy of 85.71% in the smallest dataset and a maximum of 100% accuracy in four datasets with an overall median accuracy of 97.8% across the 15 datasets, with the majority of misclassification occurring between the heterogeneous Group 3 and Group 4 subtypes. We anticipate this medulloblastoma transcriptomic subtype classifier will be broadly applicable to the cancer research and clinical communities.

annoFuse: an R Package to annotate, prioritize, and interactively explore putative oncogenic RNA fusions.

BACKGROUND: Gene fusion events are significant sources of somatic variation across adult and pediatric cancers and are some of the most clinically-effective therapeutic targets, yet low consensus of RNA-Seq fusion prediction algorithms makes therapeutic prioritization difficult. In addition, events such as polymerase read-throughs, mis-mapping due to gene homology, and fusions occurring in healthy normal tissue require informed filtering, making it difficult for researchers and clinicians to rapidly discern gene fusions that might be true underlying oncogenic drivers of a tumor and in some cases, appropriate targets for therapy. RESULTS: We developed annoFuse, an R package, and shinyFuse, a companion web application, to annotate, prioritize, and explore biologically-relevant expressed gene fusions, downstream of fusion calling. We validated annoFuse using a random cohort of TCGA RNA-Seq samples (N = 160) and achieved a 96% sensitivity for retention of high-confidence fusions (N = 603). annoFuse uses FusionAnnotator annotations to filter non-oncogenic and/or artifactual fusions. Then, fusions are prioritized if previously reported in TCGA and/or fusions containing gene partners that are known oncogenes, tumor suppressor genes, COSMIC genes, and/or transcription factors. We applied annoFuse to fusion calls from pediatric brain tumor RNA-Seq samples (N = 1028) provided as part of the Open Pediatric Brain Tumor Atlas (OpenPBTA) Project to determine recurrent fusions and recurrently-fused genes within different brain tumor histologies. annoFuse annotates protein domains using the PFAM database, assesses reciprocality, and annotates gene partners for kinase domain retention. As a standard function, reportFuse enables generation of a reproducible R Markdown report to summarize filtered fusions, visualize breakpoints and protein domains by transcript, and plot recurrent fusions within cohorts. Finally, we created shinyFuse for algorithm-agnostic interactive exploration and plotting of gene fusions. CONCLUSIONS: annoFuse provides standardized filtering and annotation for gene fusion calls from STAR-Fusion and Arriba by merging, filtering, and prioritizing putative oncogenic fusions across large cancer datasets, as demonstrated here with data from the OpenPBTA project. We are expanding the package to be widely-applicable to other fusion algorithms and expect annoFuse to provide researchers a method for rapidly evaluating, prioritizing, and translating fusion findings in patient tumors.

Aberrant splicing in B-cell acute lymphoblastic leukemia.

Aberrant splicing is a hallmark of leukemias with mutations in splicing factor (SF)-encoding genes. Here we investigated its prevalence in pediatric B-cell acute lymphoblastic leukemias (B-ALL), where SFs are not mutated. By comparing these samples to normal pro-B cells, we found thousands of aberrant local splice variations (LSVs) per sample, with 279 LSVs in 241 genes present in every comparison. These genes were enriched in RNA processing pathways and encoded ∼100 SFs, e.g. hnRNPA1. HNRNPA1 3'UTR was most pervasively mis-spliced, yielding the transcript subject to nonsense-mediated decay. To mimic this event, we knocked it down in B-lymphoblastoid cells and identified 213 hnRNPA1-regulated exon usage events comprising the hnRNPA1 splicing signature in pediatric leukemia. Some of its elements were LSVs in DICER1 and NT5C2, known cancer drivers. We searched for LSVs in other leukemia and lymphoma drivers and discovered 81 LSVs in 41 additional genes. Seventy-seven LSVs out of 81 were confirmed using two large independent B-ALL RNA-seq datasets, and the twenty most common B-ALL drivers, including NT5C2, showed higher prevalence of aberrant splicing than of somatic mutations. Thus, post-transcriptional deregulation of SF can drive widespread changes in B-ALL splicing and likely contributes to disease pathogenesis.

FBXW7ß isoform drives transcriptional activation of the proinflammatory TNF cluster in human pro-B cells.

Noncanonical exon usage plays many important roles in cellular phenotypes, but its contribution to human B-cell development remains sketchily understood. To fill this gap, we collected various B-cell fractions from bone marrow (BM) and tonsil donors, performed RNA sequencing, and examined transcript variants. We identified 150 genes that harbor local splicing variations in all pairwise comparisons. One of them encodes FBXW7, an E3 ubiquitin ligase implicated as a driver in several blood cancers. Surprisingly, we discovered that in normal human pro-B cells, the predominant transcript used an alternative first exon to produce the poorly characterized FBXW7ß isoform, previously thought to be restricted to neural tissues. The FBXW7ß transcript was also abundant in cell lines and primary samples of pediatric B-cell acute lymphoblastic leukemia (B-ALL), which originates in the BM. When overexpressed in a heterologous cell system, this transcript yielded the expected protein product, as judged by anti-FLAG immunoblotting and mass spectrometry. Furthermore, in REH B-ALL cells, FBXW7ß mRNA was the only FBXW7 isoform enriched in the polyribosome fraction. To shed light on possible functions of FBXW7ß, we used gain- and loss-of-function approaches and identified an FBXW7-dependent inflammatory gene signature, apparent in a subset of B-ALL with high FBXW7ß expression. This signature contained several members of the tumor necrosis factor superfamily, including those comprising the HLA Class III cluster (LTB, LST1, NCR3, LTA, and NFKBIL1). Our findings suggest that FBXW7ß expression drives proinflammatory responses, which could contribute to normal B-cell development, leukemogenesis, and responses to anticancer therapies.

AutoGVP: a dockerized workflow integrating ClinVar and InterVar germline sequence variant classification.

With the increasing rates of exome and whole genome sequencing, the ability to classify large sets of germline sequencing variants using up-to-date American College of Medical Genetics - Association for Molecular Pathology (ACMG-AMP) criteria is crucial. Here, we present Automated Germline Variant Pathogenicity (AutoGVP), a tool that integrates germline variant pathogenicity annotations from ClinVar and sequence variant classifications from a modified version of InterVar (PVS1 strength adjustments, removal of PP5/BP6). This tool facilitates large-scale, clinically-focused classification of germline sequence variants in a research setting.

Modulation of CD22 Protein Expression in Childhood Leukemia by Pervasive Splicing Aberrations: Implications for CD22-Directed Immunotherapies.

Downregulation of surface epitopes causes postimmunotherapy relapses in B-lymphoblastic leukemia (B-ALL). Here we demonstrate that mRNA encoding CD22 undergoes aberrant splicing in B-ALL. We describe the plasma membrane-bound CD22 Δex5-6 splice isoform, which is resistant to chimeric antigen receptor (CAR) T cells targeting the third immunoglobulin-like domain of CD22. We also describe splice variants skipping the AUG-containing exon 2 and failing to produce any identifiable protein, thereby defining an event that is rate limiting for epitope presentation. Indeed, forcing exon 2 skipping with morpholino oligonucleotides reduced CD22 protein expression and conferred resistance to the CD22-directed antibody-drug conjugate inotuzumab ozogamicin in vitro. Furthermore, among inotuzumab-treated pediatric patients with B-ALL, we identified one nonresponder in whose leukemic blasts Δex2 isoforms comprised the majority of CD22 transcripts. In a second patient, a sharp reduction in CD22 protein levels during relapse was driven entirely by increased CD22 exon 2 skipping. Thus, dysregulated CD22 splicing is a major mechanism of epitope downregulation and ensuing resistance to immunotherapy. SIGNIFICANCE: The mechanism(s) underlying downregulation of surface CD22 following CD22-directed immunotherapy remains underexplored. Our biochemical and correlative studies demonstrate that in B-ALL, CD22 expression levels are controlled by inclusion/skipping of CD22 exon 2. Thus, aberrant splicing of CD22 is an important driver/biomarker of de novo and acquired resistance to CD22-directed immunotherapies. See related commentary by Bourcier and Abdel-Wahab, p. 87. This article is highlighted in the In This Issue feature, p. 85.

Targeting CD123 in blastic plasmacytoid dendritic cell neoplasm using allogeneic anti-CD123 CAR T cells.

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a rare hematologic malignancy with poor outcomes with conventional therapy. Nearly 100% of BPDCNs overexpress interleukin 3 receptor subunit alpha (CD123). Given that CD123 is differentially expressed on the surface of BPDCN cells, it has emerged as an attractive therapeutic target. UCART123 is an investigational product consisting of allogeneic T cells expressing an anti-CD123 chimeric antigen receptor (CAR), edited with TALEN® nucleases. In this study, we examine the antitumor activity of UCART123 in preclinical models of BPDCN. We report that UCART123 have selective antitumor activity against CD123-positive primary BPDCN samples (while sparing normal hematopoietic progenitor cells) in the in vitro cytotoxicity and T cell degranulation assays; supported by the increased secretion of IFNγ by UCART123 cells when cultured in the presence of BPDCN cells. UCART123 eradicate BPDCN and result in long-term disease-free survival in a subset of primary patient-derived BPDCN xenograft mouse models. One potential challenge of CD123 targeting therapies is the loss of CD123 antigen through diverse genetic mechanisms, an event observed in one of three BPDCN PDX studied. In summary, these results provide a preclinical proof-of-principle that allogeneic UCART123 cells have potent anti-BPDCN activity.

Age-driven modulation of tRNA-derived fragments in Drosophila and their potential targets.

BACKGROUND: Development of sequencing technologies and supporting computation enable discovery of small RNA molecules that previously escaped detection or were ignored due to low count numbers. While the focus in the analysis of small RNA libraries has been primarily on microRNAs (miRNAs), recent studies have reported findings of fragments of transfer RNAs (tRFs) across a range of organisms. RESULTS: Here we describe Drosophila melanogaster tRFs, which appear to have a number of structural and functional features similar to those of miRNAs but are less abundant. As is the case with miRNAs, (i) tRFs seem to have distinct isoforms preferentially originating from 5' or 3' end of a precursor molecule (in this case, tRNA), (ii) ends of tRFs appear to contain short "seed" sequences matching conserved regions across 12 Drosophila genomes, preferentially in 3' UTRs but also in introns and exons; (iii) tRFs display specific isoform loading into Ago1 and Ago2 and thus likely function in RISC complexes; (iii) levels of loading in Ago1 and Ago2 differ considerably; and (iv) both tRF expression and loading appear to be age-dependent, indicating potential regulatory changes from young to adult organisms. CONCLUSIONS: We found that Drosophila tRF reads mapped to both nuclear and mitochondrial tRNA genes for all 20 amino acids, while previous studies have usually reported fragments from only a few tRNAs. These tRFs show a number of similarities with miRNAs, including seed sequences. Based on complementarity with conserved Drosophila regions we identified such seed sequences and their possible targets with matches in the 3'UTR regions. Strikingly, the potential target genes of the most abundant tRFs show significant Gene Ontology enrichment in development and neuronal function. The latter suggests that involvement of tRFs in the RNA interfering pathway may play a role in brain activity or brain changes with age.

The Msi Family of RNA-Binding Proteins Function Redundantly as Intestinal Oncoproteins.

Members of the Msi family of RNA-binding proteins have recently emerged as potent oncoproteins in a range of malignancies. MSI2 is highly expressed in hematopoietic cancers, where it is required for disease maintenance. In contrast to the hematopoietic system, colorectal cancers can express both Msi family members, MSI1 and MSI2. Here, we demonstrate that, in the intestinal epithelium, Msi1 and Msi2 have analogous oncogenic effects. Further, comparison of Msi1/2-induced gene expression programs and transcriptome-wide analyses of Msi1/2-RNA-binding targets reveal significant functional overlap, including induction of the PDK-Akt-mTORC1 axis. Ultimately, we demonstrate that concomitant loss of function of both MSI family members is sufficient to abrogate the growth of human colorectal cancer cells, and Msi gene deletion inhibits tumorigenesis in several mouse models of intestinal cancer. Our findings demonstrate that MSI1 and MSI2 act as functionally redundant oncoproteins required for the ontogeny of intestinal cancers.

The exoribonuclease Nibbler controls 3' end processing of microRNAs in Drosophila.

MicroRNAs (miRNAs) are endogenous noncoding small RNAs with important roles in many biological pathways; their generation and activity are under precise regulation [1-3]. Emerging evidence suggests that miRNA pathways are precisely modulated with controls at the level of transcription [4-8], processing [9-11], and stability [12, 13], with miRNA deregulation linked with diseases [14] and neurodegenerative disorders [15]. In the Drosophila miRNA biogenesis pathway, long primary miRNA transcripts undergo sequential cleavage [16-18] to release the embedded miRNAs. Mature miRNAs are then loaded into Argonaute1 (Ago1) within the RNA-induced silencing complex (RISC) [19, 20]. Intriguingly, we found that Drosophila miR-34 displays multiple isoforms that differ at the 3' end, suggesting a novel biogenesis mechanism involving 3' end processing. To define the cellular factors responsible, we performed an RNA interference (RNAi) screen and identified a putative 3'→5' exoribonuclease CG9247/nibbler essential for the generation of the smaller isoforms of miR-34. Nibbler (Nbr) interacts with Ago1 and processes miR-34 within RISC. Deep sequencing analysis revealed a larger set of multi-isoform miRNAs that are controlled by nibbler. These findings suggest that Nbr-mediated 3' end processing represents a critical step in miRNA maturation that impacts miRNA diversity.