ERα-associated translocations underlie oncogene amplifications in breast cancer.

Focal copy-number amplification is an oncogenic event. Although recent studies have revealed the complex structure1-3 and the evolutionary trajectories4 of oncogene amplicons, their origin remains poorly understood. Here we show that focal amplifications in breast cancer frequently derive from a mechanism-which we term translocation-bridge amplification-involving inter-chromosomal translocations that lead to dicentric chromosome bridge formation and breakage. In 780 breast cancer genomes, we observe that focal amplifications are frequently connected to each other by inter-chromosomal translocations at their boundaries. Subsequent analysis indicates the following model: the oncogene neighbourhood is translocated in G1 creating a dicentric chromosome, the dicentric chromosome is replicated, and as dicentric sister chromosomes segregate during mitosis, a chromosome bridge is formed and then broken, with fragments often being circularized in extrachromosomal DNAs. This model explains the amplifications of key oncogenes, including ERBB2 and CCND1. Recurrent amplification boundaries and rearrangement hotspots correlate with oestrogen receptor binding in breast cancer cells. Experimentally, oestrogen treatment induces DNA double-strand breaks in the oestrogen receptor target regions that are repaired by translocations, suggesting a role of oestrogen in generating the initial translocations. A pan-cancer analysis reveals tissue-specific biases in mechanisms initiating focal amplifications, with the breakage-fusion-bridge cycle prevalent in some and the translocation-bridge amplification in others, probably owing to the different timing of DNA break repair. Our results identify a common mode of oncogene amplification and propose oestrogen as its mechanistic origin in breast cancer.

Clonal haematopoiesis as a risk factor for therapy-related myeloid neoplasms in patients with chronic lymphocytic leukaemia treated with chemo-(immuno)therapy.

Clonal haematopoiesis of indeterminate potential (CHIP) may predispose for the development of therapy-related myeloid neoplasms (t-MN). Using target next-generation sequencing (t-NGS) panels and digital droplet polymerase chain reactions (ddPCR), we studied the myeloid gene mutation profiles of patients with chronic lymphocytic leukaemia (CLL) who developed a t-MN after treatment with chemo-(immuno)therapy. Using NGS, we detected a total of 30 pathogenic/likely pathogenic (P/LP) variants in 10 of 13 patients with a t-MN (77%, median number of variants for patient: 2, range 0-6). The prevalence of CHIP was then backtracked in paired samples taken at CLL diagnosis in eight of these patients. Six of them carried at least one CHIP-variant at the time of t-MN (median: 2, range: 1-5), and the same variants were present in the CLL sample in five cases. CHIP variants were present in 34 of 285 patients from a population-based CLL cohort, which translates into a significantly higher prevalence of CHIP in patients with a CLL who developed a t-MN, compared to the population-based cohort (5/8, 62.5% vs. 34/285, 12%, p = 0.0001). Our data show that CHIP may be considered as a novel parameter affecting treatment algorithms in patients with CLL, and highlight the potential of using chemo-free therapies in CHIP-positive cases.

Comparative analysis of targeted next-generation sequencing panels for the detection of gene mutations in chronic lymphocytic leukemia: an ERIC multi-center study.

Next-generation sequencing (NGS) has transitioned from research to clinical routine, yet the comparability of different technologies for mutation profiling remains an open question. We performed a European multicenter (n=6) evaluation of three amplicon-based NGS assays targeting 11 genes recurrently mutated in chronic lymphocytic leukemia. Each assay was assessed by two centers using 48 pre-characterized chronic lymphocytic leukemia samples; libraries were sequenced on the Illumina MiSeq instrument and bioinformatics analyses were centralized. Across all centers the median percentage of target reads ≥100x ranged from 94.2- 99.8%. In order to rule out assay-specific technical variability, we first assessed variant calling at the individual assay level i.e., pairwise analysis of variants detected amongst partner centers. After filtering for variants present in the paired normal sample and removal of PCR/sequencing artefacts, the panels achieved 96.2% (Multiplicom), 97.7% (TruSeq) and 90% (HaloPlex) concordance at a variant allele frequency (VAF) >0.5%. Reproducibility was assessed by looking at the inter-laboratory variation in detecting mutations and 107 of 115 (93% concordance) mutations were detected by all six centers, while the remaining eight variants (7%) were undetected by a single center. Notably, 6 of 8 of these variants concerned minor subclonal mutations (VAF <5%). We sought to investigate low-frequency mutations further by using a high-sensitivity assay containing unique molecular identifiers, which confirmed the presence of several minor subclonal mutations. Thus, while amplicon-based approaches can be adopted for somatic mutation detection with VAF >5%, after rigorous validation, the use of unique molecular identifiers may be necessary to reach a higher sensitivity and ensure consistent and accurate detection of low-frequency variants.

Molecular characterization of a large unselected cohort of metastatic colorectal cancers in relation to primary tumor location, rare metastatic sites and prognosis.

Background: We have reported that BRAF V600E mutations and microsatellite instability-high (MSI-H) are more prevalent in a population-based cohort of metastatic colorectal cancer (mCRC) patients than has been reported from clinical trials or hospital-based patient groups. The aim was to explore if other mutations in mCRC differ in prevalence between these cohorts in relation to mismatch repair status and primary tumor location and if presence of bone or brain metastases is associated with any mutations.Material and methods: A population-based cohort of 798 mCRC patients from three regions in Scandinavia was used. Forty-four cancer related genes were investigated in a custom designed Ampliseq hotspot panel. Differences in survival were analyzed using the Kaplan-Meier estimator and the Cox regression analysis.Results: Determination of mutations was possible in 449/501 patients for 40/44 genes. Besides BRAF V600E, seen in 19% of the tumors, none of the other mutations appeared more prevalent than in trial cohorts. BRAF V600E and MSI-H, seen in 8%, were associated with poor prognosis as was right-sided primary tumor location (39%) when compared to left-sided and rectum together; however, in a multivariable regression, only the BRAF mutation retained its statistical significance. No other mutations were associated with poor prognosis. ERBB2 alterations were more common if bone metastases were present at diagnosis (17% vs. 4%, p = .011). No association was found for brain metastases. Fifty-two percent had an alteration that is treatable with an FDA-approved targeted therapy, chiefly by EGFR-inhibitor for RAS wild-type and a check-point inhibitor for MSI-H tumors.Conclusions: Right-sided tumor location, BRAF V600E mutations, but no other investigated mutation, and MSI-H are more commonly seen in an unselected cohort than is reported from clinical patient cohorts, likely because they indicate poor prognosis. Half of the patients have a tumor that is treatable with an already FDA-approved targeted drug for mCRC.

Whole-exome sequencing in relapsing chronic lymphocytic leukemia: clinical impact of recurrent RPS15 mutations.

Fludarabine, cyclophosphamide, and rituximab (FCR) is first-line treatment of medically fit chronic lymphocytic leukemia (CLL) patients; however, despite good response rates, many patients eventually relapse. Although recent high-throughput studies have identified novel recurrent genetic lesions in adverse prognostic CLL, the mechanisms leading to relapse after FCR therapy are not completely understood. To gain insight into this issue, we performed whole-exome sequencing of sequential samples from 41 CLL patients who were uniformly treated with FCR but relapsed after a median of 2 years. In addition to mutations with known adverse-prognostic impact (TP53, NOTCH1, ATM, SF3B1, NFKBIE, and BIRC3), a large proportion of cases (19.5%) harbored mutations in RPS15, a gene encoding a component of the 40S ribosomal subunit. Extended screening, totaling 1119 patients, supported a role for RPS15 mutations in aggressive CLL, with one-third of RPS15-mutant cases also carrying TP53 aberrations. In most cases, selection of dominant, relapse-specific subclones was observed over time. However, RPS15 mutations were clonal before treatment and remained stable at relapse. Notably, all RPS15 mutations represented somatic missense variants and resided within a 7 amino-acid, evolutionarily conserved region. We confirmed the recently postulated direct interaction between RPS15 and MDM2/MDMX and transient expression of mutant RPS15 revealed defective regulation of endogenous p53 compared with wild-type RPS15. In summary, we provide novel insights into the heterogeneous genetic landscape of CLL relapsing after FCR treatment and highlight a novel mechanism underlying clinical aggressiveness involving a mutated ribosomal protein, potentially representing an early genetic lesion in CLL pathobiology.

Frequent NFKBIE deletions are associated with poor outcome in primary mediastinal B-cell lymphoma.

We recently reported a truncating deletion in the NFKBIE gene, which encodes IκBε, a negative feedback regulator of NF-κB, in clinically aggressive chronic lymphocytic leukemia (CLL). Because preliminary data indicate enrichment of NFKBIE aberrations in other lymphoid malignancies, we screened a large patient cohort (n = 1460) diagnosed with different lymphoid neoplasms. While NFKBIE deletions were infrequent in follicular lymphoma, splenic marginal zone lymphoma, and T-cell acute lymphoblastic leukemia (<2%), slightly higher frequencies were seen in diffuse large B-cell lymphoma, mantle cell lymphoma, and primary central nervous system lymphoma (3% to 4%). In contrast, a remarkably high frequency of NFKBIE aberrations (46/203 cases [22.7%]) was observed in primary mediastinal B-cell lymphoma (PMBL) and Hodgkin lymphoma (3/11 cases [27.3%]). NFKBIE-deleted PMBL patients were more often therapy refractory (P = .022) and displayed inferior outcome compared with wild-type patients (5-year survival, 59% vs 78%; P = .034); however, they appeared to benefit from radiotherapy (P =022) and rituximab-containing regimens (P = .074). NFKBIE aberrations remained an independent factor in multivariate analysis (P = .003) and when restricting the analysis to immunochemotherapy-treated patients (P = .008). Whole-exome sequencing and gene expression profiling verified the importance of NF-κB deregulation in PMBL. In summary, we identify NFKBIE aberrations as a common genetic event across B-cell malignancies and highlight NFKBIE deletions as a novel poor-prognostic marker in PMBL.

Highly similar genomic landscapes in monoclonal B-cell lymphocytosis and ultra-stable chronic lymphocytic leukemia with low frequency of driver mutations.

Despite the recent discovery of recurrent driver mutations in chronic lymphocytic leukemia, the genetic factors involved in disease onset remain largely unknown. To address this issue, we performed whole-genome sequencing in 11 individuals with monoclonal B- cell lymphocytosis, both of the low-count and high-count subtypes, and 5 patients with ultra-stable chronic lymphocytic leukemia (>10 years without progression from initial diagnosis). All three entities were indistinguishable at the genomic level exhibiting low genomic complexity and similar types of somatic mutations. Exonic mutations were not frequently identified in putative chronic lymphocytic leukemia driver genes in all settings, including low-count monoclonal B-cell lymphocytosis. To corroborate these findings, we also performed deep sequencing in 11 known frequently mutated genes in an extended cohort of 28 monoclonal B-cell lymphocytosis/chronic lymphocytic leukemia cases. Interestingly, shared mutations were detected between clonal B cells and paired polymorphonuclear cells, strengthening the notion that at least a fraction of somatic mutations may occur before disease onset, likely at the hematopoietic stem cell level. Finally, we identified previously unreported non-coding variants targeting pathways relevant to B-cell and chronic lymphocytic leukemia development, likely associated with the acquisition of the characteristic neoplastic phenotype typical of both monoclonal B-cell lymphocytosis and chronic lymphocytic leukemia.

SF3B1 mutation identifies a distinct subset of myelodysplastic syndrome with ring sideroblasts.

Refractory anemia with ring sideroblasts (RARS) is a myelodysplastic syndrome (MDS) characterized by isolated erythroid dysplasia and 15% or more bone marrow ring sideroblasts. Ring sideroblasts are found also in other MDS subtypes, such as refractory cytopenia with multilineage dysplasia and ring sideroblasts (RCMD-RS). A high prevalence of somatic mutations of SF3B1 was reported in these conditions. To identify mutation patterns that affect disease phenotype and clinical outcome, we performed a comprehensive mutation analysis in 293 patients with myeloid neoplasm and 1% or more ring sideroblasts. SF3B1 mutations were detected in 129 of 159 cases (81%) of RARS or RCMD-RS. Among other patients with ring sideroblasts, lower prevalence of SF3B1 mutations and higher prevalence of mutations in other splicing factor genes were observed (P < .001). In multivariable analyses, patients with SF3B1 mutations showed significantly better overall survival (hazard ratio [HR], .37; P = .003) and lower cumulative incidence of disease progression (HR = 0.31; P = .018) compared with SF3B1-unmutated cases. The independent prognostic value of SF3B1 mutation was retained in MDS without excess blasts, as well as in sideroblastic categories (RARS and RCMD-RS). Among SF3B1-mutated patients, coexisting mutations in DNA methylation genes were associated with multilineage dysplasia (P = .015) but had no effect on clinical outcome. TP53 mutations were frequently detected in patients without SF3B1 mutation, and were associated with poor outcome. Thus, SF3B1 mutation identifies a distinct MDS subtype that is unlikely to develop detrimental subclonal mutations and is characterized by indolent clinical course and favorable outcome.

Tumor vessel up-regulation of INSR revealed by single-cell expression analysis of the tyrosine kinome and phosphatome in human cancers.

The tyrosine kinome and phosphatome harbor oncogenes and tumor suppressor genes and important regulators of angiogenesis and tumor stroma formation. To provide a better understanding of their potential roles in cancer, we analyzed the expression of 85 tyrosine kinases and 42 tyrosine phosphatases by in situ hybridization 48 human normal and 24 tumor tissue specimens. Nine-tenths of the assessed transcripts had tumor cell expression concordant with expression array databases. Further, pan-cancer expression of AATK, PTPRK, and PTPRU and expression of PTPRS in a subset of tumors were observed. To demonstrate tumor subcompartment resolution, we validated the predicted tumor stroma-specific markers HTRA1, HTRA3, MXRA5, MXRA8, and SERPING1 in situ. In addition to known vascular and stromal markers such as PDGFRB, we observed stromal expression of PTK6 and TNS1 and vascular expression of INSR, PTPRF, PTPRG, PTPRU, and TNS1, of which INSR emerged as a tumor-specific vessel marker. This study demonstrates the feasibility of large-scale analyses to chart the transcriptome in situ in human cancers and their ability to identify novel cancer biomarkers.

Targeted next-generation sequencing in chronic lymphocytic leukemia: a high-throughput yet tailored approach will facilitate implementation in a clinical setting.

Next-generation sequencing has revealed novel recurrent mutations in chronic lymphocytic leukemia, particularly in patients with aggressive disease. Here, we explored targeted re-sequencing as a novel strategy to assess the mutation status of genes with prognostic potential. To this end, we utilized HaloPlex targeted enrichment technology and designed a panel including nine genes: ATM, BIRC3, MYD88, NOTCH1, SF3B1 and TP53, which have been linked to the prognosis of chronic lymphocytic leukemia, and KLHL6, POT1 and XPO1, which are less characterized but were found to be recurrently mutated in various sequencing studies. A total of 188 chronic lymphocytic leukemia patients with poor prognostic features (unmutated IGHV, n=137; IGHV3-21 subset #2, n=51) were sequenced on the HiSeq 2000 and data were analyzed using well-established bioinformatics tools. Using a conservative cutoff of 10% for the mutant allele, we found that 114/180 (63%) patients carried at least one mutation, with mutations in ATM, BIRC3, NOTCH1, SF3B1 and TP53 accounting for 149/177 (84%) of all mutations. We selected 155 mutations for Sanger validation (variant allele frequency, 10-99%) and 93% (144/155) of mutations were confirmed; notably, all 11 discordant variants had a variant allele frequency between 11-27%, hence at the detection limit of conventional Sanger sequencing. Technical precision was assessed by repeating the entire HaloPlex procedure for 63 patients; concordance was found for 77/82 (94%) mutations. In summary, this study demonstrates that targeted next-generation sequencing is an accurate and reproducible technique potentially suitable for routine screening, eventually as a stand-alone test without the need for confirmation by Sanger sequencing.

Contribution of de novo retroelements to birth defects and childhood cancers.

Insertion of active retroelements-L1s, Alus, and SVAs-can disrupt proper genome function and lead to various disorders including cancer. However, the role of de novo retroelements (DNRTs) in birth defects and childhood cancers has not been well characterized due to the lack of adequate data and efficient computational tools. Here, we examine whole-genome sequencing data of 3,244 trios from 12 birth defect and childhood cancer cohorts in the Gabriella Miller Kids First Pediatric Research Program. Using an improved version of our tool xTea (x-Transposable element analyzer) that incorporates a deep-learning module, we identified 162 DNRTs, as well as 2 pseudogene insertions. Several variants are likely to be causal, such as a de novo Alu insertion that led to the ablation of a whole exon in the NF1 gene in a proband with brain tumor. We observe a high de novo SVA insertion burden in both high-intolerance loss-of-function genes and exons as well as more frequent de novo Alu insertions of paternal origin. We also identify potential mosaic DNRTs from embryonic stages. Our study reveals the important roles of DNRTs in causing birth defects and predisposition to childhood cancers.

Accurate and sensitive mutational signature analysis with MuSiCal.

Mutational signature analysis is a recent computational approach for interpreting somatic mutations in the genome. Its application to cancer data has enhanced our understanding of mutational forces driving tumorigenesis and demonstrated its potential to inform prognosis and treatment decisions. However, methodological challenges remain for discovering new signatures and assigning proper weights to existing signatures, thereby hindering broader clinical applications. Here we present Mutational Signature Calculator (MuSiCal), a rigorous analytical framework with algorithms that solve major problems in the standard workflow. Our simulation studies demonstrate that MuSiCal outperforms state-of-the-art algorithms for both signature discovery and assignment. By reanalyzing more than 2,700 cancer genomes, we provide an improved catalog of signatures and their assignments, discover nine indel signatures absent in the current catalog, resolve long-standing issues with the ambiguous 'flat' signatures and give insights into signatures with unknown etiologies. We expect MuSiCal and the improved catalog to be a step towards establishing best practices for mutational signature analysis.

Predicting response to immune checkpoint blockade therapy among mismatch repair-deficient patients using mutational signatures.

Despite the overall efficacy of immune checkpoint blockade (ICB) for mismatch repair deficiency (MMRD) across tumor types, a sizable fraction of patients with MMRD still do not respond to ICB. We performed mutational signature analysis of panel sequencing data (n = 95) from MMRD cases treated with ICB. We discover that T>C-rich single base substitution (SBS) signatures-SBS26 and SBS54 from the COSMIC Mutational Signatures catalog-identify MMRD patients with significantly shorter overall survival. Tumors with a high burden of SBS26 show over-expression and enriched mutations of genes involved in double-strand break repair and other DNA repair pathways. They also display chromosomal instability (CIN), likely related to replication fork instability, leading to copy number losses that trigger immune evasion. SBS54 is associated with transcriptional activity and not with CIN, defining a distinct subtype. Consistently, cancer cell lines with a high burden of SBS26 and SBS54 are sensitive to treatments targeting pathways related to their proposed etiology. Together, our analysis offers an explanation for the heterogeneous responses to ICB among MMRD patients and supports an SBS signature-based predictor as a prognostic biomarker for differential ICB response.

A pan-tissue survey of mosaic chromosomal alterations in 948 individuals.

Genetic mutations accumulate in an organism's body throughout its lifetime. While somatic single-nucleotide variants have been well characterized in the human body, the patterns and consequences of large chromosomal alterations in normal tissues remain largely unknown. Here, we present a pan-tissue survey of mosaic chromosomal alterations (mCAs) in 948 healthy individuals from the Genotype-Tissue Expression project, augmenting RNA-based allelic imbalance estimation with haplotype phasing. We found that approximately a quarter of the individuals carry a clonally-expanded mCA in at least one tissue, with incidence strongly correlated with age. The prevalence and genome-wide patterns of mCAs vary considerably across tissue types, suggesting tissue-specific mutagenic exposure and selection pressures. The mCA landscapes in normal adrenal and pituitary glands resemble those in tumors arising from these tissues, whereas the same is not true for the esophagus and skin. Together, our findings show a widespread age-dependent emergence of mCAs across normal human tissues with intricate connections to tumorigenesis.

BTK and PLCG2 remain unmutated in one-third of patients with CLL relapsing on ibrutinib.

Patients with chronic lymphocytic leukemia (CLL) progressing on ibrutinib constitute an unmet need. Though Bruton tyrosine kinase (BTK) and PLCG2 mutations are associated with ibrutinib resistance, their frequency and relevance to progression are not fully understood. In this multicenter retrospective observational study, we analyzed 98 patients with CLL on ibrutinib (49 relapsing after an initial response and 49 still responding after ≥1 year of continuous treatment) using a next-generation sequencing (NGS) panel (1% sensitivity) comprising 13 CLL-relevant genes including BTK and PLCG2. BTK hotspot mutations were validated by droplet digital polymerase chain reaction (ddPCR) (0.1% sensitivity). By integrating NGS and ddPCR results, 32 of 49 relapsing cases (65%) carried at least 1 hotspot BTK and/or PLCG2 mutation(s); in 6 of 32, BTK mutations were only detected by ddPCR (variant allele frequency [VAF] 0.1% to 1.2%). BTK/PLCG2 mutations were also identified in 6 of 49 responding patients (12%; 5/6 VAF <10%), of whom 2 progressed later. Among the relapsing patients, the BTK-mutated (BTKmut) group was enriched for EGR2 mutations, whereas BTK-wildtype (BTKwt) cases more frequently displayed BIRC3 and NFKBIE mutations. Using an extended capture-based panel, only BRAF and IKZF3 mutations showed a predominance in relapsing cases, who were enriched for del(8p) (n = 11; 3 BTKwt). Finally, no difference in TP53 mutation burden was observed between BTKmut and BTKwt relapsing cases, and ibrutinib treatment did not favor selection of TP53-aberrant clones. In conclusion, we show that BTK/PLCG2 mutations were absent in a substantial fraction (35%) of a real-world cohort failing ibrutinib, and propose additional mechanisms contributing to resistance.

Five Percent Variant Allele Frequency Is a Reliable Reporting Threshold for TP53 Variants Detected by Next Generation Sequencing in Chronic Lymphocytic Leukemia in the Clinical Setting.

The clinical significance of small TP53 clones detected with next generation sequencing (NGS) in chronic lymphocytic leukemia is an issue of active debate. According to the official guidelines, treatment decisions should be guided only by variants with variant allele frequency (VAF) ≥10%. We present data on 325 consecutive patients with chronic lymphocytic leukemia analyzed with NGS. In total 47 pathogenic/likely pathogenic (P/LP), TP53 variants were detected in 26 patients (8%). Eleven of these (23%) were in the 5% to 10% VAF range and reported according to our institutional policy. All TP53 variants in the 5% to 10% VAF range were confirmed (100% concordance) with a second NGS panel. Our results where further validated with the performance of Sanger sequencing and digital droplet PCR (ddPCR). In 12 patients with available fluorescence in situ hybridization data and TP53 mutations within 5% to 10% VAF, deletion of chromosome 17p (del(17p)) was detectable in only 1 patient. We propose a robust diagnostic algorithm, which allows the safe detection and reporting of TP53 variants with VAF down to 5% in the clinical setting. Our study provides evidence that NGS is equally potent to detect variants with VAF 5% to 10% compared to those with VAF 10% to 15%, highlighting the urgent need for harmonization of NGS methodologies across diagnostic laboratories.

Familial platelet disorder due to germline exonic deletions in RUNX1: a diagnostic challenge with distinct alterations of the transcript isoform equilibrium.

Germline pathogenic variants in RUNX1 are associated with familial platelet disorder with predisposition to myeloid malignancies (FPD/MM) with intragenic deletions in RUNX1 accounting for almost 7% of all reported variants. We present two new pedigrees with FPD/MM carrying two different germline RUNX1 intragenic deletions. The aforementioned deletions encompass exons 1-2 and 9-10 respectively, with the exon 9-10 deletion being previously unreported. RNA sequencing of patients carrying the exon 9-10 deletion revealed a fusion with LINC00160 resulting in a change in the 3' sequence of RUNX1. Expression analysis of the transcript isoform demonstrated altered RUNX1a/b/c ratios in carriers from both families compared to controls. Our data provide evidence on the impact of intragenic RUNX1 deletions on transcript isoform expression and highlight the importance of routinely performing copy number variant analysis in patients with suspected MM with germline predisposition.