Applications of Immunogenomics to Cancer.

Cancer immunogenomics originally was framed by research supporting the hypothesis that cancer mutations generated novel peptides seen as "non-self" by the immune system. The search for these "neoantigens" has been facilitated by the combination of new sequencing technologies, specialized computational analyses, and HLA binding predictions that evaluate somatic alterations in a cancer genome and interpret their ability to produce an immune-stimulatory peptide. The resulting information can characterize a tumor's neoantigen load, its cadre of infiltrating immune cell types, the T or B cell receptor repertoire, and direct the design of a personalized therapeutic.

The next-generation sequencing revolution and its impact on genomics.

Genomics is a relatively new scientific discipline, having DNA sequencing as its core technology. As technology has improved the cost and scale of genome characterization over sequencing's 40-year history, the scope of inquiry has commensurately broadened. Massively parallel sequencing has proven revolutionary, shifting the paradigm of genomics to address biological questions at a genome-wide scale. Sequencing now empowers clinical diagnostics and other aspects of medical care, including disease risk, therapeutic identification, and prenatal testing. This Review explores the current state of genomics in the massively parallel sequencing era.

Genomic landscape of non-small cell lung cancer in smokers and never-smokers.

We report the results of whole-genome and transcriptome sequencing of tumor and adjacent normal tissue samples from 17 patients with non-small cell lung carcinoma (NSCLC). We identified 3,726 point mutations and more than 90 indels in the coding sequence, with an average mutation frequency more than 10-fold higher in smokers than in never-smokers. Novel alterations in genes involved in chromatin modification and DNA repair pathways were identified, along with DACH1, CFTR, RELN, ABCB5, and HGF. Deep digital sequencing revealed diverse clonality patterns in both never-smokers and smokers. All validated EFGR and KRAS mutations were present in the founder clones, suggesting possible roles in cancer initiation. Analysis revealed 14 fusions, including ROS1 and ALK, as well as novel metabolic enzymes. Cell-cycle and JAK-STAT pathways are significantly altered in lung cancer, along with perturbations in 54 genes that are potentially targetable with currently available drugs.

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.

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 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.

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.

EGFR internal tandem duplications in fusion-negative congenital and neonatal spindle cell tumors.

Next-generation sequencing (NGS) assays can sensitively detect somatic variation, and increasingly can enable the identification of complex structural rearrangements. A subset of infantile spindle cell sarcomas, particularly congenital mesoblastic nephromas with classic or mixed histology, have structural rearrangement in the form of internal tandem duplications (ITD) involving EGFR. We performed prospective analysis to identify EGFR ITD through clinical or research studies, as well as retrospective analysis to quantify the frequency of EGFR ITD in pediatric sarcomas. Within our institution, three tumors with EGFR ITD were prospectively identified, all occurring in patients less than 1 year of age at diagnosis, including two renal tumors and one mediastinal soft tissue tumor. These three cases exhibited both cellular and mixed cellular and classic histology. All patients had no evidence of disease progression off therapy, despite incomplete resection. To extend our analysis and quantify the frequency of EGFR ITD in pediatric sarcomas, we retrospectively analyzed a cohort of tumors (n = 90) that were previously negative for clinical RT-PCR-based fusion testing. We identified EGFR ITD in three analyzed cases, all in patients less than 1 year of age (n = 18; 3/18, 17%). Here we expand the spectrum of tumors with EGFR ITD to congenital soft tissue tumors and report an unusual example of an EGFR ITD in a tumor with cellular congenital mesoblastic nephroma histology. We also highlight the importance of appropriate test selection and bioinformatic analysis for identification of this genomic alteration that is unexpectedly common in congenital and infantile spindle cell tumors.

Infantile metastatic ependymoma with a novel molecular profile and favorable outcome to intensive chemotherapy without irradiation: Case-based review.

Ependymal tumors are the third most common brain tumor under 14 years old. Even though metastatic disease is a rare event, it affects mostly young children and carries an adverse prognosis. The factors associated with dissemination and the best treatment approach have not yet been established and there is limited published data on how to manage metastatic disease, especially in patients under 3 years of age. We provide a review of the literature on clinical characteristics and radiation-sparing treatments for metastatic ependymoma in children under 3 years of age treated. The majority (73%) of the identified cases were above 12 months old and had the PF as the primary site at diagnosis. Chemotherapy-based approaches, in different regimens, were used with radiation reserved for progression or relapse. The prognosis varied among the studies, with an average of 50%-58% overall survival. This study also describes the case of a 7-month-old boy with metastatic posterior fossa (PF) ependymoma, for whom we identified a novel SPECC1L-RAF1 gene fusion using a patient-centric comprehensive molecular profiling protocol. The patient was successfully treated with intensive induction chemotherapy followed by high-dose chemotherapy and autologous hematopoietic progenitor cell rescue (AuHSCR). Currently, the patient is in continuous remission 5 years after his diagnosis, without radiation therapy. The understanding of the available therapeutic approaches may assist physicians in their management of such patients. This report also opens the perspective of newly identified molecular alterations in metastatic ependymomas that might drive more chemo-sensitive tumors.

Manufacture and Characterization of Good Manufacturing Practice-Compliant SARS-COV-2 Cytotoxic T Lymphocytes.

BACKGROUND: Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). SARS-CoV-2 virus-specific cytotoxic T-cell lymphocytes (vCTLs) could provide a promising modality in COVID-19 treatment. We aimed to screen, manufacture, and characterize SARS-CoV-2-vCTLs generated from convalescent COVID-19 donors using the CliniMACS Cytokine Capture System (CCS). METHODS: Donor screening was done by stimulation of convalescent COVID-19 donor peripheral blood mononuclear cells with viral peptides and identification of interferonγ (IFN-γ)+ CD4 and CD8 T cells using flow cytometry. Clinical-grade SARS-CoV-2-vCTLs were manufactured using the CliniMACS CCS. The enriched SARS-CoV-2-vCTLs were characterized by T-cell receptor sequencing, mass cytometry, and transcriptome analysis. RESULTS: Of the convalescent donor blood samples, 93% passed the screening criteria for clinical manufacture. Three validation runs resulted in enriched T cells that were 79% (standard error of the mean 21%) IFN-γ+ T cells. SARS-CoV-2-vCTLs displayed a highly diverse T-cell receptor repertoire with enhancement of both memory CD8 and CD4 T cells, especially in CD8 TEM, CD4 TCM, and CD4 TEMRA cell subsets. SARS-CoV-2-vCTLs were polyfunctional with increased gene expression in T-cell function, interleukin, pathogen defense, and tumor necrosis factor superfamily pathways. CONCLUSIONS: Highly functional SARS-CoV-2-vCTLs can be rapidly generated by direct cytokine enrichment (12 hours) from convalescent donors. CLINICAL TRIALS REGISTRATION: NCT04896606.

A sequence motif within chromatin entry sites directs MSL establishment on the Drosophila X chromosome.

The Drosophila MSL complex associates with active genes specifically on the male X chromosome to acetylate histone H4 at lysine 16 and increase expression approximately 2-fold. To date, no DNA sequence has been discovered to explain the specificity of MSL binding. We hypothesized that sequence-specific targeting occurs at "chromatin entry sites," but the majority of sites are sequence independent. Here we characterize 150 potential entry sites by ChIP-chip and ChIP-seq and discover a GA-rich MSL recognition element (MRE). The motif is only slightly enriched on the X chromosome ( approximately 2-fold), but this is doubled when considering its preferential location within or 3' to active genes (>4-fold enrichment). When inserted on an autosome, a newly identified site can direct local MSL spreading to flanking active genes. These results provide strong evidence for both sequence-dependent and -independent steps in MSL targeting of dosage compensation to the male X chromosome.

Orthotopic patient-derived xenografts of paediatric solid tumours.

Paediatric solid tumours arise from endodermal, ectodermal, or mesodermal lineages. Although the overall survival of children with solid tumours is 75%, that of children with recurrent disease is below 30%. To capture the complexity and diversity of paediatric solid tumours and establish new models of recurrent disease, here we develop a protocol to produce orthotopic patient-derived xenografts at diagnosis, recurrence, and autopsy. Tumour specimens were received from 168 patients, and 67 orthotopic patient-derived xenografts were established for 12 types of cancer. The origins of the patient-derived xenograft tumours were reflected in their gene-expression profiles and epigenomes. Genomic profiling of the tumours, including detailed clonal analysis, was performed to determine whether the clonal population in the xenograft recapitulated the patient's tumour. We identified several drug vulnerabilities and showed that the combination of a WEE1 inhibitor (AZD1775), irinotecan, and vincristine can lead to complete response in multiple rhabdomyosarcoma orthotopic patient-derived xenografts tumours in vivo.

Detection of brain somatic variation in epilepsy-associated developmental lesions.

OBJECTIVE: Epilepsy-associated developmental lesions, including malformations of cortical development and low-grade developmental tumors, represent a major cause of drug-resistant seizures requiring surgical intervention in children. Brain-restricted somatic mosaicism has been implicated in the genetic etiology of these lesions; however, many contributory genes remain unidentified. METHODS: We enrolled 50 children who were undergoing epilepsy surgery into a translational research study. Resected tissue was divided for clinical neuropathologic evaluation and genomic analysis. We performed exome and RNA sequencing to identify somatic variation and we confirmed our findings using high-depth targeted DNA sequencing. RESULTS: We uncovered candidate disease-causing somatic variation affecting 28 patients (56%), as well as candidate germline variants affecting 4 patients (8%). In agreement with previous studies, we identified somatic variation affecting solute carrier family 35 member A2 (SLC35A2) and mechanistic target of rapamycin kinase (MTOR) pathway genes in patients with focal cortical dysplasia. Somatic gains of chromosome 1q were detected in 30% (3 of 10) of patients with Type I focal cortical dysplasia (FCD)s. Somatic variation in mitogen-activated protein kinase (MAPK) pathway genes (i.e., fibroblast growth factor receptor 1 [FGFR1], FGFR2, B-raf proto-oncogene, serine/threonine kinase [BRAF], and KRAS proto-oncogene, GTPase [KRAS]) was associated with low-grade epilepsy-associated developmental tumors. RNA sequencing enabled the detection of somatic structural variation that would have otherwise been missed, and which accounted for more than one-half of epilepsy-associated tumor diagnoses. Sampling across multiple anatomic regions revealed that somatic variant allele fractions vary widely within epileptogenic tissue. Finally, we identified putative disease-causing variants in genes not yet associated with focal cortical dysplasia. SIGNIFICANCE: These results further elucidate the genetic basis of structural brain abnormalities leading to focal epilepsy in children and point to new candidate disease genes.

PTEN somatic mutations contribute to spectrum of cerebral overgrowth.

Phosphatase and tensin homologue (PTEN) regulates cell growth and survival through inhibition of the mammalian target of rapamycin (MTOR) signalling pathway. Germline genetic variation of PTEN is associated with autism, macrocephaly and PTEN hamartoma tumour syndromes. The effect of developmental PTEN somatic mutations on nervous system phenotypes is not well understood, although brain somatic mosaicism of MTOR pathway genes is an emerging cause of cortical dysplasia and epilepsy in the paediatric population. Here we report two somatic variants of PTEN affecting a single patient presenting with intractable epilepsy and hemimegalencephaly that varied in clinical severity throughout the left cerebral hemisphere. High-throughput sequencing analysis of affected brain tissue identified two somatic variants in PTEN. The first variant was present in multiple cell lineages throughout the entire hemisphere and associated with mild cerebral overgrowth. The second variant was restricted to posterior brain regions and affected the opposite PTEN allele, resulting in a segmental region of more severe malformation, and the only neurons in which it was found by single-nuclei RNA-sequencing had a unique disease-related expression profile. This study reveals brain mosaicism of PTEN as a disease mechanism of hemimegalencephaly and furthermore demonstrates the varying effects of single- or bi-allelic disruption of PTEN on cortical phenotypes.

Immunological ignorance is an enabling feature of the oligo-clonal T cell response to melanoma neoantigens.

The impact of intratumoral heterogeneity (ITH) and the resultant neoantigen landscape on T cell immunity are poorly understood. ITH is a widely recognized feature of solid tumors and poses distinct challenges related to the development of effective therapeutic strategies, including cancer neoantigen vaccines. Here, we performed deep targeted DNA sequencing of multiple metastases from melanoma patients and observed ubiquitous sharing of clonal and subclonal single nucleotide variants (SNVs) encoding putative HLA class I-restricted neoantigen epitopes. However, spontaneous antitumor CD8+ T cell immunity in peripheral blood and tumors was restricted to a few clonal neoantigens featuring an oligo-/monoclonal T cell-receptor (TCR) repertoire. Moreover, in various tumors of the 4 patients examined, no neoantigen-specific TCR clonotypes were identified despite clonal neoantigen expression. Mature dendritic cell (mDC) vaccination with tumor-encoded amino acid-substituted (AAS) peptides revealed diverse neoantigen-specific CD8+ T responses, each composed of multiple TCR clonotypes. Isolation of T cell clones by limiting dilution from tumor-infiltrating lymphocytes (TILs) permitted functional validation regarding neoantigen specificity. Gene transfer of TCRαß heterodimers specific for clonal neoantigens confirmed correct TCR clonotype assignments based on high-throughput TCRBV CDR3 sequencing. Our findings implicate immunological ignorance of clonal neoantigens as the basis for ineffective T cell immunity to melanoma and support the concept that therapeutic vaccination, as an adjunct to checkpoint inhibitor treatment, is required to increase the breadth and diversity of neoantigen-specific CD8+ T cells.

Novel morphologic findings in PLAG1-rearranged soft tissue tumors.

Oncogenesis in PLAG1-rearranged tumors often results from PLAG1 transcription factor overexpression driven by promoter-swapping between constitutively expressed fusion partners. PLAG1-rearranged tumors demonstrate diverse morphologies. This study adds to this morphologic heterogeneity by introducing two tumors with PLAG1 rearrangements that display distinct histologic features. The first arose in the inguinal region of a 3-year-old, appeared well-circumscribed with a multinodular pattern, and harbored two fusions: ZFHX4-PLAG1 and CHCHD7-PLAG1. The second arose in the pelvic cavity of a 15-year-old girl, was extensively infiltrative and vascularized with an adipocytic component, and demonstrated a COL3A1-PLAG1 fusion. Both showed low-grade cytomorphology, scarce mitoses, no necrosis, and expression of CD34 and desmin. The ZFHX4-/CHCHD7-PLAG1-rearranged tumor showed no evidence of recurrence after 5 months. By contrast, the COL3A1-PLAG1-rearranged tumor quickly recurred following primary excision with positive margins; subsequent re-excision with adjuvant chemotherapy resulted in no evidence of recurrence after 2 years. While both tumors show overlap with benign and malignant fibroblastic and fibrovascular neoplasms, they also display divergent features. These cases highlight the importance of appropriate characterization in soft tissue tumors with unusual clinical and histologic characteristics.

Gastroblastoma with a novel EWSR1-CTBP1 fusion presenting in adolescence.

Gastroblastomas are rare tumors with a biphasic epithelioid/spindle cell morphology that typically present in early adulthood and have recurrent MALAT1-GLI1 fusions. We describe an adolescent patient with Wiskott-Aldrich syndrome who presented with a large submucosal gastric tumor with biphasic morphology. Despite histologic features consistent with gastroblastoma, a MALAT1-GLI1 fusion was not found in this patient's tumor; instead, comprehensive molecular profiling identified a novel EWSR1-CTBP1 fusion and no other significant genetic alterations. The tumor also overexpressed NOTCH and FGFR by RNA profiling. The novel fusion and expression profile suggest a role for epithelial-mesenchymal transition in this tumor, with potential implications for the pathogenesis of biphasic gastric tumors such as gastroblastoma.