Whole genome analysis reveals the genomic complexity in metastatic cutaneous squamous cell carcinoma

Metastatic cutaneous squamous cell carcinoma (cSCC) is associated with a high risk of recurrence and poor prognosis. There is limited published data exploring whole genome sequencing (WGS). The aim of this project was to provide the first comprehensive genomic understanding of the state of metastatic cSCC. In this study, we used WGS on matched tumor and blood DNA to detect somatic genetic alterations from 25 patients with regional metastases of head and neck cSCC. Our computational analyses interrogate clinical impacts of these genetic alterations on metastatic cSCC across the cohort for both the coding and non-coding genome. In the non-coding genome, 3UTR regions of EVC (48%), PPP1R1A (48%) and LUM (16%) were significantly functionally altered (Q-value < 0.05). Further, significant functional alterations are observed in the tumor suppressing lncRNA LINC01003 ( 68% of specimens, Q-value: 0.0158). In addition, significant recurrent copy number loss in tumor suppressor genes KANSL1 and PTPRD and gain in CALR, CCND1 and FGF3 was observed for coding regions. SNVs driver analyses predicted TP53, CDKN2A, as potential drivers of the metastasis cSCC (using 3 different tools). Indel signature analysis highlight dominance of ID signature 13 followed by ID8 & ID9. Interestingly, ID 9 has previously been shown to have no association with skin melanoma, unlike ID 13 and 8, suggesting some point of difference between these two skin-based diseases. The overall landscape of variation in metastatic cSCC is dominated by cell cycle and DNA repair disruption.

Automated prediction of the clinical impact of structural copy number variations

Copy number variants (CNVs) play an important role in many biological processes, including the development of genetic diseases, making them attractive targets for genetic analyses. The interpretation of the effect of these structural variants is a challenging problem due to highly variable numbers of gene, regulatory, or other genomic elements affected by the CNV. This led to the demand for the interpretation tools that would relieve researchers, laboratory diagnosticians, genetic counselors, and clinical geneticists from the laborious process of annotation and classification of CNVs. We designed and validated a prediction method (ISV; Interpretation of Structural Variants) that is based on boosted trees which takes into account annotations of CNVs from several publicly available databases. The presented approach achieved more than 98% prediction accuracy on both copy number loss and copy number gain variants while also allowing CNVs being assigned “uncertain” significance in predictions. We believe that ISV’s prediction capability and explainability have a great potential to guide users to more precise interpretations and classifications of CNVs.

Genomic alterations and evolution of cell clusters in metastatic invasive micropapillary carcinoma of the breast

Invasive micropapillary carcinoma (IMPC) has very high rates of lymphovascular invasion and lymph node metastasis and has been reported in several organs. However, the genomic mechanisms underlying its metastasis are unclear. Here, we perform whole-genome sequencing of tumor cell clusters from primary IMPC and paired axillary lymph node metastases. Cell clusters in multiple lymph node foci arise from a single subclone of the primary tumor. We find evidence that the monoclonal metastatic ancestor in primary IMPC shares high frequency copy-number loss of PRDM16 and IGSF9 and the copy number gain of ALDH2. Immunohistochemistry analysis further shows that low expression of IGSF9 and PRDM16 and high expression of ALDH2 are associated with lymph node metastasis and poor survival of patients with IMPC. We expect these genomic and evolutionary profiles to contribute to the accurate diagnosis of IMPC.

Pan-Genomic Sequencing Reveals Actionable CDKN2A/2B Deletions and Kataegis in Anaplastic Thyroid Carcinoma

Anaplastic thyroid carcinoma (ATC) is a lethal malignancy characterized by poor response to conventional therapies. Whole-genome sequencing (WGS) analyses of this tumor type are limited, and we therefore interrogated eight ATCs using WGS and RNA sequencing. Five out of eight cases (63%) displayed cyclin-dependent kinase inhibitor 2A (CDKN2A) abnormalities, either copy number loss (n = 4) or truncating mutations (n = 1). All four cases with loss of the CDKN2A locus (encoding p16 and p14arf) also exhibited loss of the neighboring CDKN2B gene (encoding p15ink4b), and displayed reduced CDKN2A/2B mRNA levels. Mutations in established ATC-related genes were observed, including TP53, BRAF, ARID1A, and RB1, and overrepresentation of mutations were also noted in 13 additional cancer genes. One of the more predominant mutational signatures was intimately coupled to the activity of Apolipoprotein B mRNA-editing enzyme, the catalytic polypeptide-like (APOBEC) family of cytidine deaminases implied in kataegis, a focal hypermutation phenotype, which was observed in 4/8 (50%) cases. We corroborate the roles of CDKN2A/2B in ATC development and identify kataegis as a recurrent phenomenon. Our findings pinpoint clinically relevant alterations, which may indicate response to CDK inhibitors, and focal hypermutational phenotypes that may be coupled to improved responses using immune checkpoint inhibitors.

CD8+ T Cell-Based Molecular Classification With Heterogeneous Immunogenomic Landscapes and Clinical Significance of Clear Cell Renal Cell Carcinoma

The tumor microenvironment (TME) exerts a high impact on tumor biology and immunotherapy. The heterogeneous phenotypes and the clinical significance of CD8+ T cells in TME have not been fully elucidated. Here, a comprehensive immunogenomic analysis based on multi-omics data was performed to investigate the clinical significance and tumor heterogeneity between CD8+ T cell-related molecular clusters. We identified two distinct molecular clusters of ccRCC (C1 and C2) in TCGA and validated in E-MTAB-1980 cohorts. The C1 cluster was characterized by unfavorable prognosis, increased expression levels of CD8+ T cell exhaustion markers, high immune infiltration levels as well as more immune escape mechanisms. The C2 cluster was featured by favorable prognosis, elevated expression levels of CD8+ T cell effector markers, low load of copy number loss and low frequency of 9p21.3 deletion. Moreover, the effect of molecular classifications on Nivolumab therapeutic efficacy in the CheckMate 025 cohort was examined, and the C2 cluster exhibited a better prognosis. Taken together, we determine two CD8+ T cell-related molecular clusters in ccRCC, and provide new insights for evaluating the functions of CD8+ T cells. Our molecular classification is a potential strategy for prognostic prediction and immunotherapeutic guidance for ccRCC patients.

Cold and heterogeneous T cell repertoire is associated with copy number aberrations and loss of immune genes in small-cell lung cancer

Small-cell lung cancer (SCLC) is speculated to harbor complex genomic intratumor heterogeneity (ITH) associated with high recurrence rate and suboptimal response to immunotherapy. Here, using multi-region whole exome/T cell receptor (TCR) sequencing as well as immunohistochemistry, we reveal a rather homogeneous mutational landscape but extremely cold and heterogeneous TCR repertoire in limited-stage SCLC tumors (LS-SCLCs). Compared to localized non-small cell lung cancers, LS-SCLCs have similar predicted neoantigen burden and genomic ITH, but significantly colder and more heterogeneous TCR repertoire associated with higher chromosomal copy number aberration (CNA) burden. Furthermore, copy number loss of IFN-γ pathway genes is frequently observed and positively correlates with CNA burden. Higher mutational burden, higher T cell infiltration and positive PD-L1 expression are associated with longer overall survival (OS), while higher CNA burden is associated with shorter OS in patients with LS-SCLC.

Genomic Analysis of Cutaneous CD30-Positive Lymphoproliferative Disorders

Primary cutaneous CD30+ T-cell lymphoproliferative disorders (CD30+ LPD) are the second most common cutaneous lymphomas. According to the World Health Organization (WHO), CD30+ LPD include primary cutaneous anaplastic large cell lymphoma (pcALCL) and lymphomatoid papulosis (LyP) as well as borderline lesions. pcALCL and LyP is thought to represent two ends of a spectrum of diseases that have different clinical presentations, clinical courses, and prognoses in their classic forms, but share the same histology of medium to large CD30+ atypical lymphoid cell infiltrates. Because the behavior of these entities is different clinically and prognostically, we aim to search for oncogenic genomic variants using whole exome sequencing (WES) that drive the development of LyP and pcALCL. Clinical information, pathology, immunohistochemistry, and T-cell rearrangements on six cases of LyP and five cases of pcALCL were reviewed to confirm the rendered diagnosis prior to WES of all specimens. All cases of CD30 LPD had recurrent mutations in at least one of the epigenetic modifying genes, with the most frequent mutations found in the mixed lineage methyltransferase family involved in the methylation of H3K4: CREBBP (27%), KMT2A (36%), KMT2D (36%), SETD2 (27%), and SMARCA4 (27%) (Figure). Within the JAK/STAT pathway, mutations of STAT3 were observed in 2 cases of pcALCL (n=2, 18%), STAT5B in one case of LyP (n=1, 9%) and JAK1 mutation in one case of LyP(n=1, 9%). Lastly, mutations were also identified in the T-cell signaling pathway. 3/5 cases of pcALCL demonstrated loss of function mutations in EOMES, a T-box transcription factor important in lymphocyte development. We had two cases of C-ALCL and one case of LyP with NOTCH1 mutations supporting the importance of this gene in T-cell lymphomas. TP53 mutations and copy number loss is more characteristic of aggressive lymphomas and these abnormalities were absent in pcALCL and LyP and may explain the indolent nature of CD30+ LPD. While the JAK/STAT pathway, T-signaling, and epigenetic alterations possibly play a role in the pathogenesis of these diseases, genes involved in cell-cycle control and apoptosis (ie. TP53) that are characteristic of more aggressive diseases were not identified. Extending this investigation to a larger number of samples will allow us to detect additional mutations that may help distinguish between CD30+ LPDs of LyP and pcALCL as well as other histologic mimics such as systemic ALCL and large cell transformation of MF. Figure 1 Figure 1. Disclosures Abdulla: Caris Life Licenses: Current Employment. Zain: Kiyoaw Kirin, Secura Bio, Seattle Genetics: Honoraria; Secura Bio, DaichiSankyo, Abbvie: Research Funding; Secura Bio, Ono , Legend, Kiyowa Kirin, Myeloid Therapeutics Verastem Daichi Sankyo: Consultancy.

PATH-24. COPY NUMBER ALTERATIONS IN NOTCH PATHWAY GENES ARE PROGNOSTIC IN A SUBSET OF DIFFUSE ASTROCYTIC GLIOMAS

Inactivating mutations in NOTCH1 occur in many cancer types and are frequently observed in IDH-mutant, 1p/19q-codeleted oligodendroglioma. Although the role of NOTCH1 as a tumor suppressor in diffuse glioma has become appreciated in human tissue and small animal models, the spectrum of inactivating mutations in Notch pathway genes in diffuse astrocytic gliomas has not been well described. To address this, we queried the TCGA lower-grade glioma and glioblastoma datasets to establish the extent of inactivation of Notch pathway genes, specifically by cataloging single nucleotide variants and those with copy number loss or deletion. Key alteration frequencies were found to be similar in two-independent glioma cohorts (Col, MSK). Notch pathway genes with inactivating alterations (overwhelmingly copy number loss) were present in 77% of TCGA diffuse gliomas. Across all diffuse gliomas, DLL3 loss was the most common alteration (TCGA 31%). For IDH-mutant diffuse astrocytic gliomas, JAG2 loss was the most common alteration (TCGA 23.0%, Col 35%, MSK 27%). DLL1 loss and MAML1 loss were mutually exclusive (p&lt; 0.001) in TCGA IDH-mutant astrocytomas with a combined frequency of 39% (Col 47%, MSK 56%). The presence of any alteration in the top 10 altered Notch pathway genes indicated a shorter progression-free survival (p = 0.028) for TCGA IDH-mutant diffuse astrocytomas. For IDH-wildtype diffuse astrocytic gliomas, EP300 loss was the most common inactivating alteration (TCGA 35.4%, Col 49%, MSK 38%). EP300 loss, DLL1 loss, DLL4 loss were mutually exclusive (p = 0.006) in TCGA IDH-wildtype diffuse astrocytic gliomas with a combined frequency of 61% (Col 72%, MSK 66%). The presence of alterations in any of these three genes indicated a decreased overall survival (p = 0.045) in TCGA IDH-wildtype diffuse astrocytic gliomas. Overall, loss of differential Notch pathway genes has prognostic implications in both IDH-wildtype and IDH-mutant diffuse astrocytic gliomas.

Integrated Multiomics Analyses Revealing Different Molecular Profiles Between Early- and Late-Stage Lung Adenocarcinoma

The molecular differences in genetic and epigenetic profiling between early-stage (ES) and late-stage (LS) lung adenocarcinoma (LUAD), which might help to understand cancer progression and biomarker guided precision treatment, need further be investigated. In this study, we performed comprehensive analysis using multi-omics next-generation sequencing (NGS) on tissue samples from 7 ES (stage I) and 10 LS (stage III/IV) LUAD patients to study molecular characteristics between the two groups. Characterization of the genomic and transcriptomic profiles showed stage-specific somatic mutations, copy number variations (CNVs) and differentially expressed genes (DEGs). LS samples tend to have more TP53, ERBB2 and CHD4 mutations. Gene copy number loss occurs in immune-related gene pathways in the late stage of LUAD. ATAC-seq analysis showed that LS samples harbored more open chromatin peaks around promoter regions and transcription start sites (TSS) than ES samples. We then identified the known transcription factor (TF) binding motifs for the differentially abundant ATAC-seq peaks between the ES and LS samples and found distinct regulatory mechanisms related to each stage. Furthermore, integrative analysis of ATAC-seq with WGS and RNA-seq data showed that the degree of chromatin accessibility is related to copy number changes, and the open chromatin regions could directly regulate the expression of some DEGs. In conclusion, we performed a comprehensive multi-omics analysis of the early and late stages of LUAD and highlighted some important molecular differences in regulatory mechanisms during cancer progression. Those findings help to further understand mechanism and biomarker related targeted therapy.