Cellular taxonomy of the preleukemic bone marrow niche of acute myeloid leukemia.

Leukemias arise from recurrent clonal mutations in hematopoietic stem/progenitor cells (HSPCs) that cause profound changes in the bone marrow microenvironment (BMM) favoring leukemic stem cell (LSC) growth over normal HSPCs. Understanding the cross talk between preleukemic mutated HSPCs and the BMM is critical to develop novel therapeutic strategies to prevent leukemogenesis. We hypothesize that preleukemic-LSCs (pLSCs) induce BMM changes critical for leukemogenesis. Using our AML-murine model, we performed single-cell RNA-sequencing of preleukemic BMM (pBMM) cells. We found normal HSC (nHSC)-regulating LepR+ mesenchymal stem cells, and endothelial cells were decreased, along with increases in CD55+ fibroblasts and pericytes. Preleukemic CD55+ fibroblasts had higher proliferation rates and decreased collagen expression, suggesting extracellular matrix remodeling during leukemogenesis. Importantly, co-culture assays found preleukemic CD55+ fibroblasts expanded pLSCs significantly over nHSCs. In conclusion, we have identified a distinct pBMM and a novel CD55+ fibroblast population that is expanded in pBMM that promote fitness of pLSCs over nHSCs.

Defining the transcriptome of PIK3CA-altered cells in a human capillary malformation using single cell long-read sequencing.

PIK3CA-related overgrowth spectrum (PROS) disorders are caused by somatic mosaic variants that result in constitutive activation of the phosphatidylinositol-3-kinase/AKT/mTOR pathway. Promising responses to molecularly targeted therapy have been reported, although identification of an appropriate agent can be hampered by the mosaic nature and corresponding low variant allele frequency of the causal variant. Moreover, our understanding of the molecular consequences of these variants-for example how they affect gene expression profiles-remains limited. Here we describe in vitro expansion of a human capillary malformation followed by molecular characterization using exome sequencing, single cell gene expression, and targeted long-read single cell RNA-sequencing in a patient with clinical features consistent with Megalencephaly-Capillary Malformation Syndrome (MCAP, a PROS condition). These approaches identified a targetable PIK3CA variant with expression restricted to PAX3+ fibroblast and undifferentiated keratinocyte populations. This study highlights the innovative combination of next-generation single cell sequencing methods to better understand unique transcriptomic profiles and cell types associated with MCAP, revealing molecular intricacies of this genetic syndrome.

Molecular characterization of gliomas and glioneuronal tumors amid Noonan syndrome: cancer predisposition examined.

Introduction: In the setting of pediatric and adolescent young adult cancer, increased access to genomic profiling has enhanced the detection of genetic variation associated with cancer predisposition, including germline syndromic conditions. Noonan syndrome (NS) is associated with the germline RAS pathway activating alterations and increased risk of cancer. Herein, we describe our comprehensive molecular profiling approach, the association of NS with glioma and glioneuronal tumors, and the clinical and histopathologic characteristics associated with the disease. Methods: Within an institutional pediatric cancer cohort (n = 314), molecular profiling comprised of paired somatic disease-germline comparator exome analysis, RNA sequencing, and tumor classification by DNA methylation analysis was performed. Results: Through the implementation of paired analysis, this study identified 4 of 314 (1.3%) individuals who harbored a germline PTPN11 variant associated with NS, of which 3 individuals were diagnosed with a glioma or glioneuronal tumor. Furthermore, we extend this study through collaboration with a peer institution to identify two additional individuals with NS and a glioma or glioneuronal tumor. Notably, in three of five (60%) individuals, paired genomic profiling led to a previously unrecognized diagnosis of Noonan syndrome despite an average age of cancer diagnosis of 16.8 years. The study of the disease-involved tissue identified signaling pathway dysregulation through somatic alteration of genes involved in cellular proliferation, survival, and differentiation. Discussion: Comparative pathologic findings are presented to enable an in-depth examination of disease characteristics. This comprehensive analysis highlights the association of gliomas and glioneuronal tumors with RASopathies and the potential therapeutic challenges and importantly demonstrates the utility of genomic profiling for the identification of germline cancer predisposition.

Myelomodulatory treatments augment the therapeutic benefit of oncolytic viroimmunotherapy in murine models of malignant peripheral nerve sheath tumors.

Introduction: Malignant peripheral nerve sheath tumors (MPNST) pose a significant therapeutic challenge due to high recurrence rates after surgical resection and a largely ineffective response to traditional chemotherapy. An alternative treatment strategy is oncolytic viroimmunotherapy, which can elicit a durable and systemic antitumor immune response and is Food and Drug Administration (FDA)-approved for the treatment of melanoma. Unfortunately, only a subset of patients responds completely, underscoring the need to address barriers hindering viroimmunotherapy effectiveness. Methods: Here we investigated the therapeutic utility of targeting key components of the MPNST immunosuppressive microenvironment to enhance viroimmunotherapy's antitumor efficacy in three murine models, one of which showed more immunogenic characteristics than the others. Results: Myelomodulatory therapy with pexidartinib, a small molecule inhibitor of CSF1R tyrosine kinase, and the oncolytic herpes simplex virus T-VEC exhibited the most significant increase in median survival time in the highly immunogenic model. Additionally, targeting myeloid cells with the myelomodulatory therapy trabectedin, a small molecule activator of caspase-8 dependent apoptosis, augmented the survival benefit of T-VEC in a less immunogenic MPNST model. However, tumor regressions or shrinkages were not observed. Depletion experiments confirmed that the enhanced survival benefit relied on a T cell response. Furthermore, flow cytometry analysis following combination viroimmunotherapy revealed decreased M2 macrophages and myeloid-derived suppressor cells and increased tumor-specific gp70+ CD8 T cells within the tumor microenvironment. Discussion: In summary, our findings provide compelling evidence for the potential to leverage viroimmunotherapy with myeloid cell targeting against MPNST and warrant further investigation.

Combinatorial macrophage induced innate immunotherapy against Ewing sarcoma: Turning "Two Keys" simultaneously.

BACKGROUND: Macrophages play important roles in phagocytosing tumor cells. However, tumors escape macrophage phagocytosis in part through the expression of anti-phagocytic signals, most commonly CD47. In Ewing sarcoma (ES), we found that tumor cells utilize dual mechanisms to evade macrophage clearance by simultaneously over-expressing CD47 and down-regulating cell surface calreticulin (csCRT), the pro-phagocytic signal. Here, we investigate the combination of a CD47 blockade (magrolimab, MAG) to inhibit the anti-phagocytic signal and a chemotherapy regimen (doxorubicin, DOX) to enhance the pro-phagocytic signal to induce macrophage phagocytosis of ES cells in vitro and inhibit tumor growth and metastasis in vivo. METHODS: Macrophages were derived from human peripheral blood monocytes by granulocyte-macrophage colony-stimulating factor (GM-CSF) and macrophage colony-stimulating factor (M-CSF). Flow cytometry- and microscopy-based in-vitro phagocytosis assays were performed to evaluate macrophage phagocytosis of ES cells. Annexin-V assay was performed to evaluate apoptosis. CD47 was knocked out by CRISPR/Cas9 approach. ES cell-based and patient-derived-xenograft (PDX)-based mouse models were utilized to assess the effects of MAG and/or DOX on ES tumor development and animal survival. RNA-Seq combined with CIBERSORTx analysis was utilized to identify changes in tumor cell transcriptome and tumor infiltrating immune cell profiling in MAG and/or DOX treated xenograft tumors. RESULTS: We found that MAG significantly increased macrophage phagocytosis of ES cells in vitro (p < 0.01) and had significant effect on reducing tumor burden (p < 0.01) and increasing survival in NSG mouse model (p < 0.001). The csCRT level on ES cells was significantly enhanced by DOX in a dose- and time-dependent manner (p < 0.01). Importantly, DOX combined with MAG significantly enhanced macrophage phagocytosis of ES cells in vitro (p < 0.01) and significantly decreased tumor burden (p < 0.01) and lung metastasis (p < 0.0001) and extended animal survival in vivo in two different mouse models of ES (p < 0.0001). Furthermore, we identified CD38, CD209, CD163 and CD206 as potential markers for ES-phagocytic macrophages. Moreover, we found increased M2 macrophage infiltration and decreased expression of Cd209 in the tumor microenvironment of MAG and DOX combinatorial therapy treated tumors. CONCLUSIONS: By turning "two keys" simultaneously to reactivate macrophage phagocytic activity, our data demonstrated an effective and highly translatable alternative therapeutic approach utilizing innate (tumor associated macrophages) immunotherapy against high-risk metastatic ES.

Germline susceptibility from broad genomic profiling of pediatric brain cancers.

Background: Identifying germline predisposition in CNS malignancies is of increasing clinical importance, as it contributes to diagnosis and prognosis, and determines aspects of treatment. The inclusion of germline testing has historically been limited due to challenges surrounding access to genetic counseling, complexity in acquiring a germline comparator specimen, concerns about the impact of findings, or cost considerations. These limitations were further defined by the breadth and scope of clinical testing to precisely identify complex variants as well as concerns regarding the clinical interpretation of variants including those of uncertain significance. Methods: In the course of conducting an IRB-approved protocol that performed genomic, transcriptomic and methylation-based characterization of pediatric CNS malignancies, we cataloged germline predisposition to cancer based on paired exome capture sequencing, coupled with computational analyses to identify variants in known cancer predisposition genes and interpret them relative to established clinical guidelines. Results: In certain cases, these findings refined diagnosis or prognosis or provided important information for treatment planning. Conclusions: We outline our aggregate findings on cancer predisposition within this cohort which identified 16% of individuals (27 of 168) harboring a variant predicting cancer susceptibility and contextualize the impact of these results in terms of treatment-related aspects of precision oncology.

Corrigendum: Myelomodulatory treatments augment the therapeutic benefit of oncolytic viroimmunotherapy in murine models of malignant peripheral nerve sheath tumors.

[This corrects the article DOI: 10.3389/fimmu.2024.1384623.].

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 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: Outcomes for patients with glioblastoma (GBM) 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: In this study, we report the findings of four patients enrolled onto the NeoVax clinical trial (NCT0342209). RESULTS: Immune reactivity to NeoVax neoantigens was assessed in peripheral blood mononuclear cells pre- and post-NeoVax for patients 1 to 3 using IFNγ-ELISPOT assay. A statistically significant increase in IFNγ producing T cells at the post-NeoVax time point for several neoantigens was observed. Furthermore, a post-NeoVax tumor biopsy was obtained from patient 3 and, upon evaluation, revealed evidence of infiltrating, clonally expanded T cells. CONCLUSIONS: Collectively, our findings suggest that NeoVax stimulated the 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 patients with GBM.

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 .

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.

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.

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.

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.

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.