Genomic characterization of cervical lymph node metastases in papillary thyroid carcinoma following the Chornobyl accident.

Childhood radioactive iodine exposure from the Chornobyl accident increased papillary thyroid carcinoma (PTC) risk. While cervical lymph node metastases (cLNM) are well-recognized in pediatric PTC, the PTC metastatic process and potential radiation association are poorly understood. Here, we analyze cLNM occurrence among 428 PTC with genomic landscape analyses and known drivers (131I-exposed = 349, unexposed = 79; mean age = 27.9 years). We show that cLNM are more frequent in PTC with fusion (55%) versus mutation (30%) drivers, although the proportion varies by specific driver gene (RET-fusion = 71%, BRAF-mutation = 38%, RAS-mutation = 5%). cLNM frequency is not associated with other characteristics, including radiation dose. cLNM molecular profiling (N = 47) demonstrates 100% driver concordance with matched primary PTCs and highly concordant mutational spectra. Transcriptome analysis reveals 17 differentially expressed genes, particularly in the HOXC cluster and BRINP3; the strongest differentially expressed microRNA also is near HOXC10. Our findings underscore the critical role of driver alterations and provide promising candidates for elucidating the biological underpinnings of PTC cLNM.

Genomic profiling of mycosis fungoides identifies patients at high risk of disease progression.

ABSTRACT: Mycosis fungoides (MF) is the most prevalent primary cutaneous T-cell lymphoma, with an indolent or aggressive course and poor survival. The pathogenesis of MF remains unclear, and prognostic factors in the early stages are not well established. Here, we characterized the most recurrent genomic alterations using whole-exome sequencing of 67 samples from 48 patients from Lille University Hospital (France), including 18 sequential samples drawn across stages of the malignancy. Genomic data were analyzed on the Broad Institute's Terra bioinformatics platform. We found that gain7q, gain10p15.1 (IL2RA and IL15RA), del10p11.22 (ZEB1), or mutations in JUNB and TET2 are associated with high-risk disease stages. Furthermore, gain7q, gain10p15.1 (IL2RA and IL15RA), del10p11.22 (ZEB1), and del6q16.3 (TNFAIP3) are coupled with shorter survival. Del6q16.3 (TNFAIP3) was a risk factor for progression in patients at low risk. By analyzing the clonal heterogeneity and the clonal evolution of the cohort, we defined different phylogenetic pathways of the disease with acquisition of JUNB, gain10p15.1 (IL2RA and IL15RA), or del12p13.1 (CDKN1B) at progression. These results establish the genomics and clonality of MF and identify potential patients at risk of progression, independent of their clinical stage.

Clinical Outcomes and Evolution of Clonal Hematopoiesis in Patients with Newly Diagnosed Multiple Myeloma.

Clonal hematopoiesis (CH) at time of autologous stem cell transplant (ASCT) has been shown to be associated with decreased overall survival (OS) and progression-free survival (PFS) in patients with multiple myeloma not receiving immunomodulatory drugs (IMiD). However, the significance of CH in newly diagnosed patients, including transplant ineligible patients, and its effect on clonal evolution during multiple myeloma therapy in the era of novel agents, has not been well studied. Using our new algorithm to differentiate tumor and germline mutations from CH, we detected CH in approximately 10% of 986 patients with multiple myeloma from the Clinical Outcomes in MM to Personal Assessment of Genetic Profile (CoMMpass) cohort (40/529 transplanted and 59/457 non-transplanted patients). CH was associated with increased age, risk of recurrent bacterial infections and cardiovascular disease. CH at time of multiple myeloma diagnosis was not associated with inferior OS or PFS regardless of undergoing ASCT, and all patients benefited from IMiD-based therapies, irrespective of the presence of CH. Serial sampling of 52 patients revealed the emergence of CH over a median of 3 years of treatment, increasing its prevalence to 25%, mostly with DNMT3A mutations. SIGNIFICANCE: Using our algorithm to differentiate tumor and germline mutations from CH mutations, we detected CH in approximately 10% of patients with newly diagnosed myeloma, including both transplant eligible and ineligible patients. Receiving IMiDs improved outcomes irrespective of CH status, but the prevalence of CH significantly rose throughout myeloma-directed therapy.

Genomic signatures of past and present chromosomal instability in Barrett's esophagus and early esophageal adenocarcinoma.

The progression of precancerous lesions to malignancy is often accompanied by increasing complexity of chromosomal alterations but how these alterations arise is poorly understood. Here we perform haplotype-specific analysis of chromosomal copy-number evolution in the progression of Barrett's esophagus (BE) to esophageal adenocarcinoma (EAC) on multiregional whole-genome sequencing data of BE with dysplasia and microscopic EAC foci. We identify distinct patterns of copy-number evolution indicating multigenerational chromosomal instability that is initiated by cell division errors but propagated only after p53 loss. While abnormal mitosis, including whole-genome duplication, underlies chromosomal copy-number changes, segmental alterations display signatures of successive breakage-fusion-bridge cycles and chromothripsis of unstable dicentric chromosomes. Our analysis elucidates how multigenerational chromosomal instability generates copy-number variation in BE cells, precipitates complex alterations including DNA amplifications, and promotes their independent clonal expansion and transformation. In particular, we suggest sloping copy-number variation as a signature of ongoing chromosomal instability that precedes copy-number complexity. These findings suggest copy-number heterogeneity in advanced cancers originates from chromosomal instability in precancerous cells and such instability may be identified from the presence of sloping copy-number variation in bulk sequencing data.

Genomic and transcriptomic analysis of checkpoint blockade response in advanced non-small cell lung cancer.

Anti-PD-1/PD-L1 agents have transformed the treatment landscape of advanced non-small cell lung cancer (NSCLC). To expand our understanding of the molecular features underlying response to checkpoint inhibitors in NSCLC, we describe here the first joint analysis of the Stand Up To Cancer-Mark Foundation cohort, a resource of whole exome and/or RNA sequencing from 393 patients with NSCLC treated with anti-PD-(L)1 therapy, along with matched clinical response annotation. We identify a number of associations between molecular features and outcome, including (1) favorable (for example, ATM altered) and unfavorable (for example, TERT amplified) genomic subgroups, (2) a prominent association between expression of inducible components of the immunoproteasome and response and (3) a dedifferentiated tumor-intrinsic subtype with enhanced response to checkpoint blockade. Taken together, results from this cohort demonstrate the complexity of biological determinants underlying immunotherapy outcomes and reinforce the discovery potential of integrative analysis within large, well-curated, cancer-specific cohorts.

Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome.

Richter syndrome (RS) arising from chronic lymphocytic leukemia (CLL) exemplifies an aggressive malignancy that develops from an indolent neoplasm. To decipher the genetics underlying this transformation, we computationally deconvoluted admixtures of CLL and RS cells from 52 patients with RS, evaluating paired CLL-RS whole-exome sequencing data. We discovered RS-specific somatic driver mutations (including IRF2BP2, SRSF1, B2M, DNMT3A and CCND3), recurrent copy-number alterations beyond del(9p21)(CDKN2A/B), whole-genome duplication and chromothripsis, which were confirmed in 45 independent RS cases and in an external set of RS whole genomes. Through unsupervised clustering, clonally related RS was largely distinct from diffuse large B cell lymphoma. We distinguished pathways that were dysregulated in RS versus CLL, and detected clonal evolution of transformation at single-cell resolution, identifying intermediate cell states. Our study defines distinct molecular subtypes of RS and highlights cell-free DNA analysis as a potential tool for early diagnosis and monitoring.

MinimuMM-seq: Genome Sequencing of Circulating Tumor Cells for Minimally Invasive Molecular Characterization of Multiple Myeloma Pathology.

Multiple myeloma (MM) develops from well-defined precursor stages; however, invasive bone marrow (BM) biopsy limits screening and monitoring strategies for patients. We enumerated circulating tumor cells (CTC) from 261 patients (84 monoclonal gammopathy of undetermined significance, 155 smoldering multiple myeloma, and 22 MM), with neoplastic cells detected in 84%. We developed a novel approach, MinimuMM-seq, which enables the detection of translocations and copy-number abnormalities through whole-genome sequencing of highly pure CTCs. Application to CTCs in a cohort of 51 patients, 24 with paired BM, was able to detect 100% of clinically reported BM biopsy events and could replace molecular cytogenetics for diagnostic yield and risk classification. Longitudinal sampling of CTCs in 8 patients revealed major clones could be tracked in the blood, with clonal evolution and shifting dynamics of subclones over time. Our findings provide proof of concept that CTC detection and genomic profiling could be used clinically for monitoring and managing disease in MM. SIGNIFICANCE: In this study, we established an approach enabling the enumeration and sequencing of CTCs to replace standard molecular cytogenetics. CTCs harbored the same pathognomonic MM abnormalities as BM plasma cells. Longitudinal sampling of serial CTCs was able to track clonal dynamics over time and detect the emergence of high-risk genetic subclones. This article is highlighted in the In This Issue feature, p. 247.

Clinical Efficacy and Whole-Exome Sequencing of Liquid Biopsies in a Phase IB/II Study of Bazedoxifene and Palbociclib in Advanced Hormone Receptor-Positive Breast Cancer.

PURPOSE: Sensitivity to endocrine therapy (ET) is critical for the clinical benefit from the combination of palbociclib plus ET in hormone receptor-positive/HER2-negative (HR+/HER2-) advanced breast cancer. Bazedoxifene is a third-generation selective estrogen receptor (ER) modulator and selective ER degrader with activity in preclinical models of endocrine-resistant breast cancer, including models harboring ESR1 mutations. Clinical trials in healthy women showed that bazedoxifene is well tolerated. PATIENTS AND METHODS: We conducted a phase Ib/II study of bazedoxifene plus palbociclib in patients with HR+/HER2- advanced breast cancer who progressed on prior ET (N = 36; NCT02448771). RESULTS: The study met its primary endpoint, with a clinical benefit rate of 33.3%, and the safety profile was consistent with what has previously been seen with palbociclib monotherapy. The median progression-free survival (PFS) was 3.6 months [95% confidence interval (CI), 2.0-7.2]. An activating PIK3CA mutation at baseline was associated with a shorter PFS (HR = 4.4; 95% CI, 1.5-13; P = 0.0026), but activating ESR1 mutations did not impact the PFS. Longitudinal plasma circulating tumor DNA whole-exome sequencing (WES; N = 68 plasma samples) provided an overview of the tumor heterogeneity and the subclonal genetic evolution, and identified actionable mutations acquired during treatment. CONCLUSIONS: The combination of palbociclib and bazedoxifene has clinical efficacy and an acceptable safety profile in a heavily pretreated patient population with advanced HR+/HER2- breast cancer. These results merit continued investigation of bazedoxifene in breast cancer.

Molecular map of chronic lymphocytic leukemia and its impact on outcome.

Recent advances in cancer characterization have consistently revealed marked heterogeneity, impeding the completion of integrated molecular and clinical maps for each malignancy. Here, we focus on chronic lymphocytic leukemia (CLL), a B cell neoplasm with variable natural history that is conventionally categorized into two subtypes distinguished by extent of somatic mutations in the heavy-chain variable region of immunoglobulin genes (IGHV). To build the 'CLL map,' we integrated genomic, transcriptomic and epigenomic data from 1,148 patients. We identified 202 candidate genetic drivers of CLL (109 new) and refined the characterization of IGHV subtypes, which revealed distinct genomic landscapes and leukemogenic trajectories. Discovery of new gene expression subtypes further subcategorized this neoplasm and proved to be independent prognostic factors. Clinical outcomes were associated with a combination of genetic, epigenetic and gene expression features, further advancing our prognostic paradigm. Overall, this work reveals fresh insights into CLL oncogenesis and prognostication.

Genetic subtypes of smoldering multiple myeloma are associated with distinct pathogenic phenotypes and clinical outcomes.

Smoldering multiple myeloma (SMM) is a precursor condition of multiple myeloma (MM) with significant heterogeneity in disease progression. Existing clinical models of progression risk do not fully capture this heterogeneity. Here we integrate 42 genetic alterations from 214 SMM patients using unsupervised binary matrix factorization (BMF) clustering and identify six distinct genetic subtypes. These subtypes are differentially associated with established MM-related RNA signatures, oncogenic and immune transcriptional profiles, and evolving clinical biomarkers. Three genetic subtypes are associated with increased risk of progression to active MM in both the primary and validation cohorts, indicating they can be used to better predict high and low-risk patients within the currently used clinical risk stratification models.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma.

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

Genomic Correlates of Outcome in Tumor-Infiltrating Lymphocyte Therapy for Metastatic Melanoma.

PURPOSE: Adoptive cell therapy (ACT) of tumor-infiltrating lymphocytes (TIL) historically yields a 40%-50% response rate in metastatic melanoma. However, the determinants of outcome are largely unknown. EXPERIMENTAL DESIGN: We investigated tumor-based genomic correlates of overall survival (OS), progression-free survival (PFS), and response to therapy by interrogating tumor samples initially collected to generate TIL infusion products. RESULTS: Whole-exome sequencing (WES) data from 64 samples indicated a positive correlation between neoantigen load and OS, but not PFS or response to therapy. RNA sequencing analysis of 34 samples showed that expression of PDE1C, RTKN2, and NGFR was enriched in responders who had improved PFS and OS. In contrast, the expression of ELFN1 was enriched in patients with unfavorable response, poor PFS and OS, whereas enhanced methylation of ELFN1 was observed in patients with favorable outcomes. Expression of ELFN1, NGFR, and PDE1C was mainly found in cancer-associated fibroblasts and endothelial cells in tumor tissues across different cancer types in publicly available single-cell RNA sequencing datasets, suggesting a role for elements of the tumor microenvironment in defining the outcome of TIL therapy. CONCLUSIONS: Our findings suggest that transcriptional features of melanomas correlate with outcomes after TIL therapy and may provide candidates to guide patient selection.

ZBTB33 is mutated in clonal hematopoiesis and myelodysplastic syndromes and impacts RNA splicing.

Clonal hematopoiesis results from somatic mutations in cancer driver genes in hematopoietic stem cells. We sought to identify novel drivers of clonal expansion using an unbiased analysis of sequencing data from 84,683 persons and identified common mutations in the 5-methylcytosine reader, ZBTB33, as well as in YLPM1, SRCAP, and ZNF318. We also identified these mutations at low frequency in myelodysplastic syndrome patients. Zbtb33 edited mouse hematopoietic stem and progenitor cells exhibited a competitive advantage in vivo and increased genome-wide intron retention. ZBTB33 mutations potentially link DNA methylation and RNA splicing, the two most commonly mutated pathways in clonal hematopoiesis and MDS.

Lack of transgenerational effects of ionizing radiation exposure from the Chernobyl accident.

Effects of radiation exposure from the Chernobyl nuclear accident remain a topic of interest. We investigated germline de novo mutations (DNMs) in children born to parents employed as cleanup workers or exposed to occupational and environmental ionizing radiation after the accident. Whole-genome sequencing of 130 children (born 1987-2002) and their parents did not reveal an increase in the rates, distributions, or types of DNMs relative to the results of previous studies. We find no elevation in total DNMs, regardless of cumulative preconception gonadal paternal [mean = 365 milligrays (mGy), range = 0 to 4080 mGy] or maternal (mean = 19 mGy, range = 0 to 550 mGy) exposure to ionizing radiation. Thus, we conclude that, over this exposure range, evidence is lacking for a substantial effect on germline DNMs in humans, suggesting minimal impact from transgenerational genetic effects.

Radiation-related genomic profile of papillary thyroid carcinoma after the Chernobyl accident.

The 1986 Chernobyl nuclear power plant accident increased papillary thyroid carcinoma (PTC) incidence in surrounding regions, particularly for radioactive iodine (131I)-exposed children. We analyzed genomic, transcriptomic, and epigenomic characteristics of 440 PTCs from Ukraine (from 359 individuals with estimated childhood 131I exposure and 81 unexposed children born after 1986). PTCs displayed radiation dose-dependent enrichment of fusion drivers, nearly all in the mitogen-activated protein kinase pathway, and increases in small deletions and simple/balanced structural variants that were clonal and bore hallmarks of nonhomologous end-joining repair. Radiation-related genomic alterations were more pronounced for individuals who were younger at exposure. Transcriptomic and epigenomic features were strongly associated with driver events but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth after environmental radiation exposure.

Clonal hematopoiesis is associated with adverse outcomes in multiple myeloma patients undergoing transplant.

Multiple myeloma (MM) is a plasma-cell neoplasm that is treated with high-dose chemotherapy, autologous stem cell transplant (ASCT) and long-term immunomodulatory drug (IMiD) maintenance. The presence of somatic mutations in the peripheral blood is termed clonal hematopoiesis of indeterminate potential (CHIP) and is associated with adverse outcomes. Targeted sequencing of the stem cell product from 629 MM patients treated by ASCT at the Dana-Farber Cancer Institute (2003-2011) detects CHIP in 136/629 patients (21.6%). The most commonly mutated genes are DNMT3A, TET2, TP53, ASXL1 and PPM1D. Twenty-one from fifty-six patients (3.3%) receiving first-line IMiD maintenance develop a therapy-related myeloid neoplasm (TMN). However, regardless of CHIP status, the use of IMiD maintenance associates with improved PFS and OS. In those not receiving IMiD maintenance, CHIP is associated with decreased overall survival (OS) (HR:1.34, p = 0.02) and progression free survival (PFS) (HR:1.45, p < 0.001) due to an increase in MM progression.

Genomic Profiling of Smoldering Multiple Myeloma Identifies Patients at a High Risk of Disease Progression.

PURPOSE: Smoldering multiple myeloma (SMM) is a precursor condition of multiple myeloma (MM) with a 10% annual risk of progression. Various prognostic models exist for risk stratification; however, those are based on solely clinical metrics. The discovery of genomic alterations that underlie disease progression to MM could improve current risk models. METHODS: We used next-generation sequencing to study 214 patients with SMM. We performed whole-exome sequencing on 166 tumors, including 5 with serial samples, and deep targeted sequencing on 48 tumors. RESULTS: We observed that most of the genetic alterations necessary for progression have already been acquired by the diagnosis of SMM. Particularly, we found that alterations of the mitogen-activated protein kinase pathway (KRAS and NRAS single nucleotide variants [SNVs]), the DNA repair pathway (deletion 17p, TP53, and ATM SNVs), and MYC (translocations or copy number variations) were all independent risk factors of progression after accounting for clinical risk staging. We validated these findings in an external SMM cohort by showing that patients who have any of these three features have a higher risk of progressing to MM. Moreover, APOBEC associated mutations were enriched in patients who progressed and were associated with a shorter time to progression in our cohort. CONCLUSION: SMM is a genetically mature entity whereby most driver genetic alterations have already occurred, which suggests the existence of a right-skewed model of genetic evolution from monoclonal gammopathy of undetermined significance to MM. We identified and externally validated genomic predictors of progression that could distinguish patients at high risk of progression to MM and, thus, improve on the precision of current clinical models.

Genomic analyses of flow-sorted Hodgkin Reed-Sternberg cells reveal complementary mechanisms of immune evasion.

Classical Hodgkin lymphoma (cHL) is composed of rare malignant Hodgkin Reed-Sternberg (HRS) cells within an extensive, but ineffective, inflammatory/immune cell infiltrate. HRS cells exhibit near-universal somatic copy gains of chromosome 9p/9p24.1, which increase expression of the programmed cell death protein 1 (PD-1) ligands. To define genetic mechanisms of response and resistance to PD-1 blockade and identify complementary treatment targets, we performed whole-exome sequencing of flow cytometry-sorted HRS cells from 23 excisional biopsies of newly diagnosed cHLs, including 8 Epstein-Barr virus-positive (EBV+) tumors. We identified significantly mutated cancer candidate genes (CCGs) as well as somatic copy number alterations and structural variations and characterized their contribution to disease-defining immune evasion mechanisms and nuclear factor κB (NF-κB), JAK/STAT, and PI3K signaling pathways. EBV- cHLs had a higher prevalence of genetic alterations in the NF-κB and major histocompatibility complex class I antigen presentation pathways. In this young cHL cohort (median age, 26 years), we identified a predominant mutational signature of spontaneous deamination of cytosine- phosphate-guanines ("Aging"), in addition to apolipoprotein B mRNA editing catalytic polypeptide-like, activation-induced cytidine deaminase, and microsatellite instability (MSI)-associated hypermutation. In particular, the mutational burden in EBV- cHLs was among the highest reported, similar to that of carcinogen-induced tumors. Together, the overall high mutational burden, MSI-associated hypermutation, and newly identified genetic alterations represent additional potential bases for the efficacy of PD-1 blockade in cHL. Of note, recurrent cHL alterations, including B2M, TNFAIP3, STAT6, GNA13, and XPO1 mutations and 2p/2p15, 6p21.32, 6q23.3, and 9p/9p24.1 copy number alterations, were also identified in >20% of primary mediastinal B-cell lymphomas, highlighting shared pathogenetic mechanisms in these diseases.

Genomic analyses of PMBL reveal new drivers and mechanisms of sensitivity to PD-1 blockade.

Primary mediastinal large B-cell lymphomas (PMBLs) are aggressive tumors that typically present as large mediastinal masses in young women. PMBLs share clinical, transcriptional, and molecular features with classical Hodgkin lymphoma (cHL), including constitutive activation of nuclear factor κB (NF-κB), JAK/STAT signaling, and programmed cell death protein 1 (PD-1)-mediated immune evasion. The demonstrated efficacy of PD-1 blockade in relapsed/refractory PMBLs led to recent approval by the US Food and Drug Administration and underscored the importance of characterizing targetable genetic vulnerabilities in this disease. Here, we report a comprehensive analysis of recurrent genetic alterations -somatic mutations, somatic copy number alterations, and structural variants-in a cohort of 37 newly diagnosed PMBLs. We identified a median of 9 genetic drivers per PMBL, including known and newly identified components of the JAK/STAT and NF-κB signaling pathways and frequent B2M alterations that limit major histocompatibility complex class I expression, as in cHL. PMBL also exhibited frequent, newly identified driver mutations in ZNF217 and an additional epigenetic modifier, EZH2. The majority of these alterations were clonal, which supports their role as early drivers. In PMBL, we identified several previously uncharacterized molecular features that may increase sensitivity to PD-1 blockade, including high tumor mutational burden, microsatellite instability, and an apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) mutational signature. The shared genetic features between PMBL and cHL provide a framework for analyzing the mechanism of action of PD-1 blockade in these related lymphoid malignancies.

RNA sequence analysis reveals macroscopic somatic clonal expansion across normal tissues.

How somatic mutations accumulate in normal cells is poorly understood. A comprehensive analysis of RNA sequencing data from ~6700 samples across 29 normal tissues revealed multiple somatic variants, demonstrating that macroscopic clones can be found in many normal tissues. We found that sun-exposed skin, esophagus, and lung have a higher mutation burden than other tested tissues, which suggests that environmental factors can promote somatic mosaicism. Mutation burden was associated with both age and tissue-specific cell proliferation rate, highlighting that mutations accumulate over both time and number of cell divisions. Finally, normal tissues were found to harbor mutations in known cancer genes and hotspots. This study provides a broad view of macroscopic clonal expansion in human tissues, thus serving as a foundation for associating clonal expansion with environmental factors, aging, and risk of disease.