Improvement of tumor neoantigen detection by high field asymmetric waveform ion mobility mass spectrometry.

Cancer neoantigens have been shown to elicit cancer-specific T-cell responses and have garnered much attention for their roles in both spontaneous and therapeutically induced antitumor responses. Mass spectrometry (MS) profiling of tumor immunopeptidomes has been used, in part, to identify MHC-bound mutant neoantigen ligands. However, under standard conditions, MS-based detection of such rare but clinically relevant neoantigens is relatively insensitive, requiring 300 million cells or more. Here, to quantitatively define the minimum detectable amounts of therapeutically relevant MHC-I and MHC-II neoantigen peptides, we analyzed different dilutions of immunopeptidomes isolated from the well-characterized T3 mouse methylcholanthrene (MCA)-induced cell line by MS. Using either data-dependent acquisition (DDA) or parallel reaction monitoring (PRM), we established the minimum amount of material required to detect the major T3 neoantigens in the presence or absence of high field asymmetric waveform ion mobility spectrometry (FAIMS). This analysis yielded a 14-fold enhancement of sensitivity in detecting the major T3 MHC-I neoantigen (mLama4) with FAIMS-PRM compared with PRM without FAIMS, allowing ex-vivo detection of this neoantigen from an individual 100 mg T3 tumor. These findings were then extended to two other independent MCA-sarcoma lines (1956 and F244). This study demonstrates that FAIMS substantially increases the sensitivity of MS-based characterization of validated neoantigens from tumors.

A novel IKZF1 variant in a family with autosomal dominant CVID: A case for expanding exon coverage in inborn errors of immunity.

Common variable immune deficiency (CVID) is a heterogenous group of disorders characterized by varying degrees of hypogammaglobulinemia, recurrent infections, and autoimmunity. Currently, pathogenic variants are identified in approximately 20-30% of CVID cases. Here we report a 3-generation family with autosomal dominant Common Variable Immunodeficiency (CVID) diagnosed in 9 affected individuals. Although primary immune deficiency panels and exome sequencing were non-diagnostic, whole genome sequencing revealed a novel, pathogenic c.499C > T: p.His167Tyr variant in IKZF1, a critical regulator of B cell development. Functional testing done through pericentromeric heterochromatin localization and light shift chemiluminescent electrophoretic mobility shift assay confirmed the variant's deleterious effect via a haploinsufficiency mechanism. Our findings expand the spectrum of known IKZF1 mutations and contribute to a more comprehensive understanding of CVID's genetic heterogeneity. Furthermore, this case underscores the importance of considering whole genome sequencing for comprehensive genetic diagnosis when concern for a monogenic inborn errors of immunity is high.

Integrating multisector molecular characterization into personalized peptide vaccine design for patients with newly diagnosed glioblastoma.

PURPOSE: Glioblastoma (GBM) patient outcomes remain poor despite multimodality treatment with surgery, radiation, and chemotherapy. There are few immunotherapy options due to the lack of tumor immunogenicity. Several clinical trials have reported promising results with cancer vaccines. To date, studies have used data from a single tumor site to identify targetable antigens, but this approach limits the antigen pool and is antithetical to the heterogeneity of GBM. We have implemented multisector sequencing to increase the pool of neoantigens across the GBM genomic landscape that can be incorporated into personalized peptide vaccines called NeoVax. PATIENTS AND METHODS: Here, we report the findings of four subjects enrolled onto the NeoVax clinical trial (NCT0342209). RESULTS: Immune reactivity to NeoVax neoantigens was assessed in peripheral blood mononuclear cells (PBMCs) pre- and post-NeoVax for subjects 1-3 using IFNg-ELISPOT assay. A statistically significant increase in IFNg producing T cells at the post-NeoVax time point for several neoantigens was observed. Furthermore, a post-NeoVax tumor biopsy was obtained from subject 3 and, upon evaluation, revealed evidence of infiltrating, clonally expanded T cells. CONCLUSIONS: Collectively, our findings suggest NeoVax did stimulate expansion of neoantigen-specific effector T cells and provide encouraging results to aid in the development of future neoantigen vaccine-based clinical trials in patients with GBM. Herein, we demonstrate the feasibility of incorporating multisector sampling in cancer vaccine design and provide information on the clinical applicability of clonality, distribution, and immunogenicity of the neoantigen landscape in GBM patients.

Trabectedin Enhances Oncolytic Virotherapy by Reducing Barriers to Virus Spread and Cytotoxic Immunity in Preclinical Pediatric Bone Sarcoma.

We previously reported that the DNA alkylator and transcriptional-blocking chemotherapeutic agent trabectedin enhances oncolytic herpes simplex viroimmunotherapy in human sarcoma xenograft models, though the mechanism remained to be elucidated. Here we report trabectedin disrupts the intrinsic cellular anti-viral response which increases viral transcript spread throughout the human tumor cells. We also extended our synergy findings to syngeneic murine sarcoma models, which are poorly susceptible to virus infection. In the absence of robust virus replication, we found trabectedin enhanced viroimmunotherapy efficacy by reducing immunosuppressive macrophages and stimulating granzyme expression in infiltrating T and NK cells to cause immune-mediated tumor regressions. Thus, trabectedin enhances both the direct virus-mediated killing of tumor cells and the viral-induced activation of cytotoxic effector lymphocytes to cause tumor regressions across models. Our data provide a strong rationale for clinical translation as both mechanisms should be simultaneously active in human patients.

SNVstory: inferring genetic ancestry from genome sequencing data.

BACKGROUND: Genetic ancestry, inferred from genomic data, is a quantifiable biological parameter. While much of the human genome is identical across populations, it is estimated that as much as 0.4% of the genome can differ due to ancestry. This variation is primarily characterized by single nucleotide variants (SNVs), which are often unique to specific genetic populations. Knowledge of a patient's genetic ancestry can inform clinical decisions, from genetic testing and health screenings to medication dosages, based on ancestral disease predispositions. Nevertheless, the current reliance on self-reported ancestry can introduce subjectivity and exacerbate health disparities. While genomic sequencing data enables objective determination of a patient's genetic ancestry, existing approaches are limited to ancestry inference at the continental level. RESULTS: To address this challenge, and create an objective, measurable metric of genetic ancestry we present SNVstory, a method built upon three independent machine learning models for accurately inferring the sub-continental ancestry of individuals. We also introduce a novel method for simulating individual samples from aggregate allele frequencies from known populations. SNVstory includes a feature-importance scheme, unique among open-source ancestral tools, which allows the user to track the ancestral signal broadcast by a given gene or locus. We successfully evaluated SNVstory using a clinical exome sequencing dataset, comparing self-reported ethnicity and race to our inferred genetic ancestry, and demonstrate the capability of the algorithm to estimate ancestry from 36 different populations with high accuracy. CONCLUSIONS: SNVstory represents a significant advance in methods to assign genetic ancestry, opening the door to ancestry-informed care. SNVstory, an open-source model, is packaged as a Docker container for enhanced reliability and interoperability. It can be accessed from https://github.com/nch-igm/snvstory .

A LINE-1 mediated deletion resulting in germline retinoblastoma predisposition.

Retinoblastoma is an ocular cancer associated with genomic variation in the RB1 gene. In individuals with bilateral retinoblastoma, a germline variant in RB1 is identified in virtually all cases. We describe herein an individual with bilateral retinoblastoma for whom multiple clinical lab assays performed by outside commercial laboratories failed to identify a germline RB1 variant. Paired tumor/normal exome sequencing, long-read whole genome sequencing, and long-read isoform sequencing was performed on a translational research basis ultimately identified a germline likely de novo Long Interspersed Nuclear Element (LINE)-1 mediated deletion resulting in a premature stop of translation of RB1 as the underlying genetic cause of retinoblastoma in this individual. Based on these research findings, the LINE-1 mediated deletion was confirmed via Sanger sequencing in our clinical laboratory, and results were reported in the patient's medical record to allow for appropriate genetic counseling.

Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity.

BACKGROUND: Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms. RESULTS: In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms. CONCLUSION: This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

Overview of modern genomic tools for diagnosis and precision therapy of childhood solid cancers.

PURPOSE OF REVIEW: The application of technology and computational analyses to generate new data types from pediatric solid cancers is transforming diagnostic accuracy. This review provides an overview of such new capabilities in the pursuit of improved treatment for essentially rare and underserved diseases that are the highest cause of mortality in children over one year of age. Sophisticated ways of identifying therapeutic vulnerabilities for highly personalized treatment are presented alongside cutting-edge disease response monitoring by liquid biopsy. RECENT FINDINGS: Precision molecular profiling data are now being combined with conventional pathology-based evaluation of pediatric cancer tissues. The resulting diagnostic information can be used to guide therapeutic decision-making, including the use of small molecule inhibitors and of immunotherapies. Integrating somatic and germline variant profiles constitutes a critical component of this emerging paradigm, as does tissue-of-origin derivation from methylation profiling, and rapid screening of potential therapies. These new approaches are poised for use in disease response and therapy resistance monitoring. SUMMARY: The integration of clinical molecular profiling data with pathology can provide a highly precise diagnosis, identify therapeutic vulnerabilities, and monitor patient responses, providing next steps toward precision oncology for improved outcomes, including reducing lifelong treatment-related sequelae.

Childhood cancer data initiative: Status report.

In March 2023, over 800 researchers, clinicians, patients, survivors, and advocates from the pediatric oncology community met to discuss the progress of the National Cancer Institute's Childhood Cancer Data Initiative. We present here the status of the initiative's efforts in building its data ecosystem and updates on key programs, especially the Molecular Characterization Initiative and the planned Coordinated National Initiative for Rare Cancers in Children and Young Adults. These activities aim to improve access to childhood cancer data, foster collaborations, facilitate integrative data analysis, and expand access to molecular characterization, ultimately leading to the development of innovative therapeutic approaches.

DGIdb 5.0: rebuilding the drug-gene interaction database for precision medicine and drug discovery platforms.

The Drug-Gene Interaction Database (DGIdb, https://dgidb.org) is a publicly accessible resource that aggregates genes or gene products, drugs and drug-gene interaction records to drive hypothesis generation and discovery for clinicians and researchers. DGIdb 5.0 is the latest release and includes substantial architectural and functional updates to support integration into clinical and drug discovery pipelines. The DGIdb service architecture has been split into separate client and server applications, enabling consistent data access for users of both the application programming interface (API) and web interface. The new interface was developed in ReactJS, and includes dynamic visualizations and consistency in the display of user interface elements. A GraphQL API has been added to support customizable queries for all drugs, genes, annotations and associated data. Updated documentation provides users with example queries and detailed usage instructions for these new features. In addition, six sources have been added and many existing sources have been updated. Newly added sources include ChemIDplus, HemOnc, NCIt (National Cancer Institute Thesaurus), Drugs@FDA, HGNC (HUGO Gene Nomenclature Committee) and RxNorm. These new sources have been incorporated into DGIdb to provide additional records and enhance annotations of regulatory approval status for therapeutics. Methods for grouping drugs and genes have been expanded upon and developed as independent modular normalizers during import. The updates to these sources and grouping methods have resulted in an improvement in FAIR (findability, accessibility, interoperability and reusability) data representation in DGIdb.

Novel inherited CDX2 variant segregating in a family with diverse congenital malformations of the genitourinary system.

Anorectal malformations (ARMs) constitute a group of congenital defects of the gastrointestinal and urogenital systems. They affect males and females, with an estimated worldwide prevalence of 1 in 5000 live births. These malformations are clinically heterogeneous and can be part of a syndromic presentation (syndromic ARM) or as a nonsyndromic entity (nonsyndromic ARM). Despite the well-recognized heritability of nonsyndromic ARM, the genetic etiology in most patients is unknown. In this study, we describe three siblings with diverse congenital anomalies of the genitourinary system, anemia, delayed milestones, and skeletal anomalies. Genome sequencing identified a novel, paternally inherited heterozygous Caudal type Homeobox 2 (CDX2) variant (c.722A > G (p.Glu241Gly)), that was present in all three affected siblings. The variant identified in this family is absent from population databases and predicted to be damaging by most in silico pathogenicity tools. So far, only two other reports implicate variants in CDX2 with ARMs. Remarkably, the individuals described in these studies had similar clinical phenotypes and genetic alterations in CDX2 CDX2 encodes a transcription factor and is considered the master regulator of gastrointestinal development. This variant maps to the homeobox domain of the encoded protein, which is critical for interaction with DNA targets. Our finding provides a potential molecular diagnosis for this family's condition and supports the role of CDX2 in anorectal anomalies. It also highlights the clinical heterogeneity and variable penetrance of ARM predisposition variants, another well-documented phenomenon. Finally, it underscores the diagnostic utility of genomic profiling of ARMs to identify the genetic etiology of these defects.

Post-zygotic rescue of meiotic errors causes brain mosaicism and focal epilepsy.

Somatic mosaicism is a known cause of neurological disorders, including developmental brain malformations and epilepsy. Brain mosaicism is traditionally attributed to post-zygotic genetic alterations arising in fetal development. Here we describe post-zygotic rescue of meiotic errors as an alternate origin of brain mosaicism in patients with focal epilepsy who have mosaic chromosome 1q copy number gains. Genomic analysis showed evidence of an extra parentally derived chromosome 1q allele in the resected brain tissue from five of six patients. This copy number gain is observed only in patient brain tissue, but not in blood or buccal cells, and is strongly enriched in astrocytes. Astrocytes carrying chromosome 1q gains exhibit distinct gene expression signatures and hyaline inclusions, supporting a novel genetic association for astrocytic inclusions in epilepsy. Further, these data demonstrate an alternate mechanism of brain chromosomal mosaicism, with parentally derived copy number gain isolated to brain, reflecting rescue in other tissues during development.

Expanding the Clinical Utility of Targeted RNA Sequencing Panels beyond Gene Fusions to Complex, Intragenic Structural Rearrangements.

Gene fusions are a form of structural rearrangement well established as driver events in pediatric and adult cancers. The identification of such events holds clinical significance in the refinement, prognostication, and provision of treatment in cancer. Structural rearrangements also extend beyond fusions to include intragenic rearrangements, such as internal tandem duplications (ITDs) or exon-level deletions. These intragenic events have been increasingly implicated as cancer-promoting events. However, the detection of intragenic rearrangements may be challenging to resolve bioinformatically with short-read sequencing technologies and therefore may not be routinely assessed in panel-based testing. Within an academic clinical laboratory, over three years, a total of 608 disease-involved samples (522 hematologic malignancy, 86 solid tumors) underwent clinical testing using Anchored Multiplex PCR (AMP)-based RNA sequencing. Hematologic malignancies were evaluated using a custom Pan-Heme 154 gene panel, while solid tumors were assessed using a custom Pan-Solid 115 gene panel. Gene fusions, ITDs, and intragenic deletions were assessed for diagnostic, prognostic, or therapeutic significance. When considering gene fusions alone, we report an overall diagnostic yield of 36% (37% hematologic malignancy, 41% solid tumors). When including intragenic structural rearrangements, the overall diagnostic yield increased to 48% (48% hematologic malignancy, 45% solid tumor). We demonstrate the clinical utility of reporting structural rearrangements, including gene fusions and intragenic structural rearrangements, using an AMP-based RNA sequencing panel.

Ten-Patient Trial: Remarkable Responses in Pediatric Cancers.

A clinical trial of nivolumab in 10 patients with pediatric cancer with high tumor mutational burden demonstrated complete responses in 50% of patients. This result recapitulates multiple clinical trial results in high mutation-burden adult cancers and may redefine best practice in the setting of germline DNA mismatch repair-based susceptibility. See related article by Das et al., p. 4770.

Oncolytic virus-driven immune remodeling revealed in mouse medulloblastomas at single cell resolution.

Oncolytic viruses, modified for tumor-restricted infection, are a promising cancer immunotherapeutic, yet much remains to be understood about factors driving their activity and outcome in the tumor microenvironment. Here, we report that oncolytic herpes simplex virus C134, previously found to exert T cell-dependent efficacy in mouse models of glioblastoma, exerts T cell-independent efficacy in mouse models of medulloblastoma, indicating this oncolytic virus uses different mechanisms in different tumors. We investigated C134's behavior in mouse medulloblastomas, using single cell RNA sequencing to map C134-induced gene expression changes across cell types, timepoints, and medulloblastoma subgroup models at whole-transcriptome resolution. Our work details substantial oncolytic virus-induced transcriptional remodeling of medulloblastoma-infiltrating immune cells, 10 subpopulations of monocytes and macrophages collectively demonstrating M1-like responses to C134, and suggests C134 be investigated as a potential new therapy for medulloblastoma.

Genetic characterization of primary and metastatic high-grade serous ovarian cancer tumors reveals distinct features associated with survival.

High-grade serous ovarian cancer (HGSC) is the most lethal histotype of ovarian cancer and the majority of cases present with metastasis and late-stage disease. Over the last few decades, the overall survival for patients has not significantly improved, and there are limited targeted treatment options. We aimed to better characterize the distinctions between primary and metastatic tumors based on short- or long-term survival. We characterized 39 matched primary and metastatic tumors by whole exome and RNA sequencing. Of these, 23 were short-term (ST) survivors (overall survival (OS) < 3.5 years) and 16 were long-term (LT) survivors (OS > 5 years). We compared somatic mutations, copy number alterations, mutational burden, differential gene expression, immune cell infiltration, and gene fusion predictions between the primary and metastatic tumors and between ST and LT survivor cohorts. There were few differences in RNA expression between paired primary and metastatic tumors, but significant differences between the transcriptomes of LT and ST survivors in both their primary and metastatic tumors. These findings will improve the understanding of the genetic variation in HGSC that exist between patients with different prognoses and better inform treatments by identifying new targets for drug development.

Genetically modified IL2 bone-marrow-derived myeloid cells reprogram the glioma immunosuppressive tumor microenvironment.

Gliomas are one of the leading causes of cancer-related death in the adolescent and young adult (AYA) population. Two-thirds of AYA glioma patients are affected by low-grade gliomas (LGGs), but there are no specific treatments. Malignant progression is supported by the immunosuppressive stromal component of the tumor microenvironment (TME) exacerbated by M2 macrophages and a paucity of cytotoxic T cells. A single intravenous dose of engineered bone-marrow-derived myeloid cells that release interleukin-2 (GEMys-IL2) was used to treat mice with LGGs. Our results demonstrate that GEMys-IL2 crossed the blood-brain barrier, infiltrated the TME, and reprogrammed the immune cell composition and transcriptome. Moreover, GEMys-IL2 extended survival in an LGG immunocompetent mouse model. Here, we report the efficacy of an in vivo approach that demonstrates the potential for a cell-mediated innate immunotherapy designed to enhance the recruitment of activated effector T and natural killer cells within the glioma TME.

A novel transcriptional signature identifies T-cell infiltration in high-risk paediatric cancer.

BACKGROUND: Molecular profiling of the tumour immune microenvironment (TIME) has enabled the rational choice of immunotherapies in some adult cancers. In contrast, the TIME of paediatric cancers is relatively unexplored. We speculated that a more refined appreciation of the TIME in childhood cancers, rather than a reliance on commonly used biomarkers such as tumour mutation burden (TMB), neoantigen load and PD-L1 expression, is an essential prerequisite for improved immunotherapies in childhood solid cancers. METHODS: We combined immunohistochemistry (IHC) with RNA sequencing and whole-genome sequencing across a diverse spectrum of high-risk paediatric cancers to develop an alternative, expression-based signature associated with CD8+ T-cell infiltration of the TIME. Furthermore, we explored transcriptional features of immune archetypes and T-cell receptor sequencing diversity, assessed the relationship between CD8+ and CD4+ abundance by IHC and deconvolution predictions and assessed the common adult biomarkers such as neoantigen load and TMB. RESULTS: A novel 15-gene immune signature, Immune Paediatric Signature Score (IPASS), was identified. Using this signature, we estimate up to 31% of high-risk cancers harbour infiltrating T-cells. In addition, we showed that PD-L1 protein expression is poorly correlated with PD-L1 RNA expression and TMB and neoantigen load are not predictive of T-cell infiltration in paediatrics. Furthermore, deconvolution algorithms are only weakly correlated with IHC measurements of T-cells. CONCLUSIONS: Our data provides new insights into the variable immune-suppressive mechanisms dampening responses in paediatric solid cancers. Effective immune-based interventions in high-risk paediatric cancer will require individualised analysis of the TIME.