The Virtual Child.

SUMMARY: We are building the world's first Virtual Child-a computer model of normal and cancerous human development at the level of each individual cell. The Virtual Child will "develop cancer" that we will subject to unlimited virtual clinical trials that pinpoint, predict, and prioritize potential new treatments, bringing forward the day when no child dies of cancer, giving each one the opportunity to lead a full and healthy life.

Interleukin-15-armored GPC3-CAR T cells for patients with solid cancers.

Interleukin-15 (IL15) promotes the survival of T lymphocytes and enhances the antitumor properties of CAR T cells in preclinical models of solid neoplasms in which CAR T cells have limited efficacy1-4. Glypican-3 (GPC3) is expressed in a group of solid cancers5-10, and here we report the first evaluation in humans of the effects of IL15 co-expression on GPC3-CAR T cells. Cohort 1 patients (NCT02905188/NCT02932956) received GPC3-CAR T cells, which were safe but produced no objective antitumor responses and reached peak expansion at two weeks. Cohort 2 patients (NCT05103631/NCT04377932) received GPC3-CAR T cells that co-expressed IL15 (15.CAR), which mediated significantly increased cell expansion and induced a disease control rate of 66% and antitumor response rate of 33%. Infusion of 15.CAR T cells was associated with increased incidence of cytokine release syndrome, which was rapidly ameliorated by activation of the inducible caspase 9 safety switch. Compared to non-responders, tumor-infiltrating 15.CAR T cells from responders showed repression of SWI/SNF epigenetic regulators and upregulation of FOS and JUN family members as well as genes related to type I interferon signaling. Collectively, these results demonstrate that IL15 increases the expansion, intratumoral survival, and antitumor activity of GPC3-CAR T cells in patients.

ALK upregulates POSTN and WNT signaling to drive neuroblastoma.

Neuroblastoma is the most common extracranial solid tumor of childhood. While MYCN and mutant anaplastic lymphoma kinase (ALKF1174L) cooperate in tumorigenesis, how ALK contributes to tumor formation remains unclear. Here, we used a human stem cell-based model of neuroblastoma. Mis-expression of ALKF1174L and MYCN resulted in shorter latency compared to MYCN alone. MYCN tumors resembled adrenergic, while ALK/MYCN tumors resembled mesenchymal, neuroblastoma. Transcriptomic analysis revealed enrichment in focal adhesion signaling, particularly the extracellular matrix genes POSTN and FN1 in ALK/MYCN tumors. Patients with ALK-mutant tumors similarly demonstrated elevated levels of POSTN and FN1. Knockdown of POSTN, but not FN1, delayed adhesion and suppressed proliferation of ALK/MYCN tumors. Furthermore, loss of POSTN reduced ALK-dependent activation of WNT signaling. Reciprocally, inhibition of the WNT pathway reduced expression of POSTN and growth of ALK/MYCN tumor cells. Thus, ALK drives neuroblastoma in part through a feedforward loop between POSTN and WNT signaling.

BPTF cooperates with MYCN and MYC to link neuroblastoma cell cycle control to epigenetic cellular states.

The nucleosome remodeling factor BPTF is required for the deployment of the MYC-driven transcriptional program. Deletion of one Bptf allele delays tumor progression in mouse models of pancreatic cancer and lymphoma. In neuroblastoma, MYCN cooperates with the transcriptional core regulatory circuitry (CRC). High BPTF levels are associated with high-risk features and decreased survival. BPTF depletion results in a dramatic decrease of cell proliferation. Bulk RNA-seq, single-cell sequencing, and tissue microarrays reveal a positive correlation of BPTF and CRC transcription factor expression. Immunoprecipitation/mass spectrometry shows that BPTF interacts with MYCN and the CRC proteins. Genome-wide distribution analysis of BPTF and CRC in neuroblastoma reveals a dual role for BPTF: 1) it co-localizes with MYCN/MYC at the promoter of genes involved in cell cycle and 2) it co-localizes with the CRC at super-enhancers to regulate cell identity. The critical role of BPTF across neuroblastoma subtypes supports its relevance as a therapeutic target.

Percutaneous biopsy for the diagnosis, risk stratification, and molecular profiling of neuroblastoma: A single-center retrospective study.

PURPOSE: To determine whether percutaneous core needle biopsy (PCNB) is adequate for the diagnosis and full molecular characterization of newly diagnosed neuroblastoma. MATERIALS AND METHODS: Patients with newly diagnosed neuroblastoma who underwent PCNB in interventional radiology at a single center over a 5-year period were included. Pre-procedure imaging and procedure details were reviewed. Rates of diagnostic success and sufficiency for International Neuroblastoma Pathology Classification (INPC), risk stratification, and evaluation of genomic markers utilized in the Children's Oncology Group risk stratification, and status of the anaplastic lymphoma kinase (ALK) gene were assessed. RESULTS: Thirty-five patients (13 females, median age 2.4 years [interquartile range, IQR: 0.9-4.4] and median weight 12.4 kg [IQR: 9.6-18]) were included. Most had International Neuroblastoma Risk Group Stage M disease (n = 22, 63%). Median longest axis of tumor target was 8.8 cm [IQR: 6.1-12]. A 16-gauge biopsy instrument was most often used (n = 20, 57%), with a median of 20 cores [IQR: 13-23] obtained. Twenty-five specimens were assessed for adequacy, and 14 procedures utilized contrast-enhanced ultrasound guidance. There were two post-procedure bleeds (5.7%). Thirty-four of 35 procedures (97%) were sufficient for histopathologic diagnosis and risk stratification, 94% (n = 32) were sufficient for INPC, and 85% (n = 29) were sufficient for complete molecular characterization, including ALK testing. Biologic information was otherwise obtained from bone marrow (4/34, 12%) or surgery (1/34, 2.9%). The number of cores did not differ between patients with sufficient versus insufficient biopsies. CONCLUSION: In this study, obtaining multiple cores with PCNB resulted in a high rate of diagnosis and successful molecular profiling for neuroblastoma.

BARD1 germline variants induce haploinsufficiency and DNA repair defects in neuroblastoma.

BACKGROUND: High-risk neuroblastoma is a complex genetic disease that is lethal in more than 50% of patients despite intense multimodal therapy. Through genome-wide association studies (GWAS) and next-generation sequencing, we have identified common single nucleotide polymorphisms and rare, pathogenic or likely pathogenic germline loss-of-function variants in BARD1 enriched in neuroblastoma patients. The functional implications of these findings remain poorly understood. METHODS: We correlated BARD1 genotype with expression in normal tissues and neuroblastomas, along with the burden of DNA damage in tumors. To validate the functional consequences of germline pathogenic or likely pathogenic BARD1 variants, we used CRISPR-Cas9 to generate isogenic neuroblastoma (IMR-5) and control (RPE1) cellular models harboring heterozygous BARD1 loss-of-function variants (R112*, R150*, E287fs, and Q564*) and quantified genomic instability in these cells via next-generation sequencing and with functional assays measuring the efficiency of DNA repair. RESULTS: Both common and rare neuroblastoma-associated BARD1 germline variants were associated with lower levels of BARD1 mRNA and an increased burden of DNA damage. Using isogenic heterozygous BARD1 loss-of-function variant cellular models, we functionally validated this association with inefficient DNA repair. BARD1 loss-of-function variant isogenic cells exhibited reduced efficiency in repairing Cas9-induced DNA damage, ineffective RAD51 focus formation at DNA double-strand break sites, and enhanced sensitivity to cisplatin and poly (ADP-ribose) polymerase (PARP) inhibition both in vitro and in vivo. CONCLUSIONS: Taken together, we demonstrate that germline BARD1 variants disrupt DNA repair fidelity. This is a fundamental molecular mechanism contributing to neuroblastoma initiation that may have important therapeutic implications.

Germline pathogenic variants in neuroblastoma patients are enriched in BARD1 and predict worse survival.

BACKGROUND: Neuroblastoma is an embryonal cancer of the developing sympathetic nervous system. The genetic contribution of rare pathogenic or likely pathogenic germline variants in patients without a family history remains unclear. METHODS: Germline DNA sequencing was performed on 786 neuroblastoma patients. The frequency of rare cancer predisposition gene pathogenic or likely pathogenic variants in patients was compared with 2 cancer-free control cohorts. Matched tumor DNA sequencing was evaluated for second hits, and germline DNA array data from 5585 neuroblastoma patients and 23 505 cancer-free control children were analyzed to identify rare germline copy number variants. Patients with germline pathogenic or likely pathogenic variants were compared with those without to test for association with clinical characteristics, tumor features, and survival. RESULTS: We observed 116 pathogenic or likely pathogenic variants involving 13.9% (109 of 786) of neuroblastoma patients, representing a statistically significant excess burden compared with cancer-free participants (odds ratio [OR] = 1.60, 95% confidence interval [CI] = 1.27 to 2.00). BARD1 harbored the most statistically significant enrichment of pathogenic or likely pathogenic variants (OR = 32.30, 95% CI = 6.44 to 310.35). Rare germline copy number variants disrupting BARD1 were identified in patients but absent in cancer-free participants (OR = 29.47, 95% CI = 1.52 to 570.70). Patients harboring a germline pathogenic or likely pathogenic variant had a worse overall survival compared with those without (P = 8.6 x 10-3). CONCLUSIONS: BARD1 is an important neuroblastoma predisposition gene harboring both common and rare germline pathogenic or likely pathogenic variations. The presence of any germline pathogenic or likely pathogenic variant in a cancer predisposition gene was independently predictive of worse overall survival. As centers move toward paired tumor-normal sequencing at diagnosis, efforts should be made to centralize data and provide an infrastructure to support cooperative longitudinal prospective studies of germline pathogenic variation.

Surface and Global Proteome Analyses Identify ENPP1 and Other Surface Proteins as Actionable Immunotherapeutic Targets in Ewing Sarcoma.

PURPOSE: Ewing sarcoma is the second most common bone sarcoma in children, with 1 case per 1.5 million in the United States. Although the survival rate of patients diagnosed with localized disease is approximately 70%, this decreases to approximately 30% for patients with metastatic disease and only approximately 10% for treatment-refractory disease, which have not changed for decades. Therefore, new therapeutic strategies are urgently needed for metastatic and refractory Ewing sarcoma. EXPERIMENTAL DESIGN: This study analyzed 19 unique Ewing sarcoma patient- or cell line-derived xenografts (from 14 primary and 5 metastatic specimens) using proteomics to identify surface proteins for potential immunotherapeutic targeting. Plasma membranes were enriched using density gradient ultracentrifugation and compared with a reference standard of 12 immortalized non-Ewing sarcoma cell lines prepared in a similar manner. In parallel, global proteome analysis was carried out on each model to complement the surfaceome data. All models were analyzed by Tandem Mass Tags-based mass spectrometry to quantify identified proteins. RESULTS: The surfaceome and global proteome analyses identified 1,131 and 1,030 annotated surface proteins, respectively. Among surface proteins identified, both approaches identified known Ewing sarcoma-associated proteins, including IL1RAP, CD99, STEAP1, and ADGRG2, and many new cell surface targets, including ENPP1 and CDH11. Robust staining of ENPP1 was demonstrated in Ewing sarcoma tumors compared with other childhood sarcomas and normal tissues. CONCLUSIONS: Our comprehensive proteomic characterization of the Ewing sarcoma surfaceome provides a rich resource of surface-expressed proteins in Ewing sarcoma. This dataset provides the preclinical justification for exploration of targets such as ENPP1 for potential immunotherapeutic application in Ewing sarcoma. See related commentary by Bailey, p. 934.

A proteogenomic surfaceome study identifies DLK1 as an immunotherapeutic target in neuroblastoma.

Cancer immunotherapies have produced remarkable results in B-cell malignancies; however, optimal cell surface targets for many solid cancers remain elusive. Here, we present an integrative proteomic, transcriptomic, and epigenomic analysis of tumor specimens along with normal tissues to identify biologically relevant cell surface proteins that can serve as immunotherapeutic targets for neuroblastoma, an often-fatal childhood cancer of the developing nervous system. We apply this approach to human-derived cell lines (N=9) and cell/patient-derived xenograft (N=12) models of neuroblastoma. Plasma membrane-enriched mass spectrometry identified 1,461 cell surface proteins in cell lines and 1,401 in xenograft models, respectively. Additional proteogenomic analyses revealed 60 high-confidence candidate immunotherapeutic targets and we prioritized Delta-like canonical notch ligand 1 (DLK1) for further study. High expression of DLK1 directly correlated with the presence of a super-enhancer spanning the DLK1 locus. Robust cell surface expression of DLK1 was validated by immunofluorescence, flow cytometry, and immunohistochemistry. Short hairpin RNA mediated silencing of DLK1 in neuroblastoma cells resulted in increased cellular differentiation. ADCT-701, a DLK1-targeting antibody-drug conjugate (ADC), showed potent and specific cytotoxicity in DLK1-expressing neuroblastoma xenograft models. Moreover, DLK1 is highly expressed in several adult cancer types, including adrenocortical carcinoma (ACC), pheochromocytoma/paraganglioma (PCPG), hepatoblastoma, and small cell lung cancer (SCLC), suggesting potential clinical benefit beyond neuroblastoma. Taken together, our study demonstrates the utility of comprehensive cancer surfaceome characterization and credentials DLK1 as an immunotherapeutic target. Highlights: Plasma membrane enriched proteomics defines surfaceome of neuroblastomaMulti-omic data integration prioritizes DLK1 as a candidate immunotherapeutic target in neuroblastoma and other cancersDLK1 expression is driven by a super-enhancer DLK1 silencing in neuroblastoma cells results in cellular differentiation ADCT-701, a DLK1-targeting antibody-drug conjugate, shows potent and specific cytotoxicity in DLK1-expressing neuroblastoma preclinical models.

Norepinephrine transporter and vesicular monoamine transporter 2 tumor expression as a predictor of response to 131 I-MIBG in patients with relapsed/refractory neuroblastoma.

BACKGROUND: Prior studies suggest that norepinephrine transporter (NET) and vesicular monoamine transporter 2 (VMAT2) mediate meta-iodobenzylguanidine (MIBG) uptake and retention in neuroblastoma tumors. We evaluated the relationship between NET and VMAT2 tumor expression and clinical response to 131 I-MIBG therapy in patients with neuroblastoma. METHODS: Immunohistochemistry (IHC) was used to evaluate NET and VMAT2 protein expression levels on archival tumor samples (obtained at diagnosis or relapse) from patients with relapsed or refractory neuroblastoma treated with 131 I-MIBG. A composite protein expression H-score was determined by multiplying a semi-quantitative intensity value (0-3+) by the percentage of tumor cells expressing the protein. RESULTS: Tumor samples and clinical data were available for 106 patients, of whom 28.3% had partial response (PR) or higher. NET H-score was not significantly associated with response (≥PR), though the percentage of tumor cells expressing NET was lower among responders (median 80% for ≥PR vs. 90% for

Structural principles of peptide-centric chimeric antigen receptor recognition guide therapeutic expansion.

Peptide-centric chimeric antigen receptors (PC-CARs) recognize oncoprotein epitopes displayed by cell-surface human leukocyte antigens (HLAs) and offer a promising strategy for targeted cancer therapy. We have previously developed a PC-CAR targeting a neuroblastoma-associated PHOX2B peptide, leading to robust tumor cell lysis restricted by two common HLA allotypes. Here, we determine the 2.1-angstrom crystal structure of the PC-CAR-PHOX2B-HLA-A*24:02-ß2m complex, which reveals the basis for antigen-specific recognition through interactions with CAR complementarity-determining regions (CDRs). This PC-CAR adopts a diagonal docking mode, where interactions with both conserved and polymorphic HLA framework residues permit recognition of multiple HLA allotypes from the A9 serological cross-reactive group, covering a combined global population frequency of up to 46.7%. Biochemical binding assays, molecular dynamics simulations, and structural and functional analyses demonstrate that high-affinity PC-CAR recognition of cross-reactive pHLAs necessitates the presentation of a specific peptide backbone, where subtle structural adaptations of the peptide are critical for high-affinity complex formation, and CAR T cell killing. Our results provide a molecular blueprint for engineering CARs with optimal recognition of tumor-associated antigens in the context of different HLAs, while minimizing cross-reactivity with self-epitopes.

Advancing childhood cancer research through young investigator and advocate collaboration.

Cancer advocates and researchers share the same goal of driving science forward to create new therapies to cure more patients. The power of combining cancer researchers and advocates has become of increased importance due to their complementary expertise. Therefore, advocacy is a critical component of grant structures and has become embedded into the Stand Up 2 Cancer (SU2C) applications. To date, the optimal way to combine these skillsets and experiences to benefit the cancer community is currently unknown. The Saint Baldrick's Foundation (SBF)-SU2C now called St. Baldrick's Empowering Pediatric Immunotherapies for Childhood Cancer (EPICC) Team is comprised of a collaborative network across nine institutions in the United States and Canada. Since SU2C encourages incorporating advocacy into the team structure, we have assembled a diverse team of advocates and scientists by nominating a young investigator (YI) and advocate from each site. In order to further bridge this interaction beyond virtual monthly and yearly in person meetings, we have developed a questionnaire and conducted interviews. The questionnaire is focused on understanding each member's experience at the intersection between science/advocacy, comparing to previous experiences, providing advice on incorporating advocacy into team science and discussing how we can build on our work. Through creating a YI and advocate infrastructure, we have cultivated a supportive environment for meaningful conversation that impacts the entire research team. We see this as a model for team science by combining expertise to drive innovation forward and positively impact pediatric cancer patients, and perhaps those with adult malignancies. Significance: Questionnaire results show both advocates and YI's see this structure to be valuable and beneficial. YI's communicated their research to a non-scientific audience and learned advocate's experience. This was their first advocacy experience for most YIs. Advocates learned more about the research being conducted to provide hope. They can also aid with fundraising, publicity and lobbying. This collaboration improves science communication, designing patient-friendly clinical trials and sharing experience across institutions.

Targeting of intracellular oncoproteins with peptide-centric CARs.

The majority of oncogenic drivers are intracellular proteins, constraining their immunotherapeutic targeting to mutated peptides (neoantigens) presented by individual human leukocyte antigen (HLA) allotypes1. However, most cancers have a modest mutational burden that is insufficient for generating responses using neoantigen-based therapies2,3. Neuroblastoma is a paediatric cancer that harbours few mutations and is instead driven by epigenetically deregulated transcriptional networks4. Here we show that the neuroblastoma immunopeptidome is enriched with peptides derived from proteins essential for tumorigenesis. We focused on targeting the unmutated peptide QYNPIRTTF discovered on HLA-A*24:02, which is derived from the neuroblastoma-dependency gene and master transcriptional regulator PHOX2B. To target QYNPIRTTF, we developed peptide-centric chimeric antigen receptors (PC-CARs) through a counter panning strategy using predicted potentially cross-reactive peptides. We further proposed that PC-CARs can recognize peptides on additional HLA allotypes when presenting a similar overall molecular surface. Informed by our computational modelling results, we show that PHOX2B PC-CARs also recognize QYNPIRTTF presented by HLA-A*23:01, the most common non-A2 allele in people with African ancestry. Finally, we demonstrate potent and specific killing of neuroblastoma cells expressing these HLAs in vitro and complete tumour regression in mice. These data suggest that PC-CARs have the potential to expand the pool of immunotherapeutic targets to include non-immunogenic intracellular oncoproteins and allow targeting through additional HLA allotypes in a clinical setting.

Preclinical development of a chimeric antigen receptor T cell therapy targeting FGFR4 in rhabdomyosarcoma.

Pediatric patients with relapsed or refractory rhabdomyosarcoma (RMS) have dismal cure rates, and effective therapy is urgently needed. The oncogenic receptor tyrosine kinase fibroblast growth factor receptor 4 (FGFR4) is highly expressed in RMS and lowly expressed in healthy tissues. Here, we describe a second-generation FGFR4-targeting chimeric antigen receptor (CAR), based on an anti-human FGFR4-specific murine monoclonal antibody 3A11, as an adoptive T cell treatment for RMS. The 3A11 CAR T cells induced robust cytokine production and cytotoxicity against RMS cell lines in vitro. In contrast, a panel of healthy human primary cells failed to activate 3A11 CAR T cells, confirming the selectivity of 3A11 CAR T cells against tumors with high FGFR4 expression. Finally, we demonstrate that 3A11 CAR T cells are persistent in vivo and can effectively eliminate RMS tumors in two metastatic and two orthotopic models. Therefore, our study credentials CAR T cell therapy targeting FGFR4 to treat patients with RMS.

A MYCN-independent mechanism mediating secretome reprogramming and metastasis in MYCN-amplified neuroblastoma.

MYCN amplification (MNA) is a defining feature of high-risk neuroblastoma (NB) and predicts poor prognosis. However, whether genes within or in close proximity to the MYCN amplicon also contribute to MNA+ NB remains poorly understood. Here, we identify that GREB1, a transcription factor encoding gene neighboring the MYCN locus, is frequently coexpressed with MYCN and promotes cell survival in MNA+ NB. GREB1 controls gene expression independently of MYCN, among which we uncover myosin 1B (MYO1B) as being highly expressed in MNA+ NB and, using a chick chorioallantoic membrane (CAM) model, as a crucial regulator of invasion and metastasis. Global secretome and proteome profiling further delineates MYO1B in regulating secretome reprogramming in MNA+ NB cells, and the cytokine MIF as an important pro-invasive and pro-metastatic mediator of MYO1B activity. Together, we have identified a putative GREB1-MYO1B-MIF axis as an unconventional mechanism promoting aggressive behavior in MNA+ NB and independently of MYCN.

Polyclonal lymphoid expansion drives paraneoplastic autoimmunity in neuroblastoma.

Neuroblastoma is a lethal childhood solid tumor of developing peripheral nerves. Two percent of children with neuroblastoma develop opsoclonus myoclonus ataxia syndrome (OMAS), a paraneoplastic disease characterized by cerebellar and brainstem-directed autoimmunity but typically with outstanding cancer-related outcomes. We compared tumor transcriptomes and tumor-infiltrating T and B cell repertoires from 38 OMAS subjects with neuroblastoma to 26 non-OMAS-associated neuroblastomas. We found greater B and T cell infiltration in OMAS-associated tumors compared to controls and showed that both were polyclonal expansions. Tertiary lymphoid structures (TLSs) were enriched in OMAS-associated tumors. We identified significant enrichment of the major histocompatibility complex (MHC) class II allele HLA-DOB∗01:01 in OMAS patients. OMAS severity scores were associated with the expression of several candidate autoimmune genes. We propose a model in which polyclonal auto-reactive B lymphocytes act as antigen-presenting cells and drive TLS formation, thereby supporting both sustained polyclonal T cell-mediated anti-tumor immunity and paraneoplastic OMAS neuropathology.

Integrative Genomic Analyses Identify LncRNA Regulatory Networks across Pediatric Leukemias and Solid Tumors.

Long noncoding RNAs (lncRNA) play an important role in gene regulation and contribute to tumorigenesis. While pan-cancer studies of lncRNA expression have been performed for adult malignancies, the lncRNA landscape across pediatric cancers remains largely uncharted. Here, we curated RNA sequencing data for 1,044 pediatric leukemia and extracranial solid tumors and integrated paired tumor whole genome sequencing and epigenetic data in relevant cell line models to explore lncRNA expression, regulation, and association with cancer. A total of 2,657 lncRNAs were robustly expressed across six pediatric cancers, including 1,142 exhibiting histotype-elevated expression. DNA copy number alterations contributed to lncRNA dysregulation at a proportion comparable to protein coding genes. Application of a multidimensional framework to identify and prioritize lncRNAs impacting gene networks revealed that lncRNAs dysregulated in pediatric cancer are associated with proliferation, metabolism, and DNA damage hallmarks. Analysis of upstream regulation via cell type-specific transcription factors further implicated distinct histotype-elevated and developmental lncRNAs. Integration of these analyses prioritized lncRNAs for experimental validation, and silencing of TBX2-AS1, the top-prioritized neuroblastoma-specific lncRNA, resulted in significant growth inhibition of neuroblastoma cells, confirming the computational predictions. Taken together, these data provide a comprehensive characterization of lncRNA regulation and function in pediatric cancers and pave the way for future mechanistic studies. SIGNIFICANCE: Comprehensive characterization of lncRNAs in pediatric cancer leads to the identification of highly expressed lncRNAs across childhood cancers, annotation of lncRNAs showing histotype-specific elevated expression, and prediction of lncRNA gene regulatory networks.

ROGUE: an R Shiny app for RNA sequencing analysis and biomarker discovery.

BACKGROUND: The growing power and ever decreasing cost of RNA sequencing (RNA-Seq) technologies have resulted in an explosion of RNA-Seq data production. Comparing gene expression values within RNA-Seq datasets is relatively easy for many interdisciplinary biomedical researchers; however, user-friendly software applications increase the ability of biologists to efficiently explore available datasets. RESULTS: Here, we describe ROGUE (RNA-Seq Ontology Graphic User Environment, https://marisshiny. RESEARCH: chop.edu/ROGUE/ ), a user-friendly R Shiny application that allows a biologist to perform differentially expressed gene analysis, gene ontology and pathway enrichment analysis, potential biomarker identification, and advanced statistical analyses. We use ROGUE to identify potential biomarkers and show unique enriched pathways between various immune cells. CONCLUSIONS: User-friendly tools for the analysis of next generation sequencing data, such as ROGUE, will allow biologists to efficiently explore their datasets, discover expression patterns, and advance their research by allowing them to develop and test hypotheses.