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.

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.

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.

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.

Consistent B Cell Receptor Immunoglobulin Features Between Siblings in Familial Chronic Lymphocytic Leukemia.

Key processes in the onset and evolution of chronic lymphocytic leukemia (CLL) are thought to include chronic (antigenic) activation of mature B cells through the B cell receptor (BcR), signals from the microenvironment, and acquisition of genetic alterations. Here we describe three families in which two or more siblings were affected by CLL. We investigated whether there are immunogenetic similarities in the leukemia-specific immunoglobulin heavy (IGH) and light (IGL/IGK) chain gene rearrangements of the siblings in each family. Furthermore, we performed array analysis to study if similarities in CLL-associated chromosomal aberrations are present within each family and screened for somatic mutations using paired tumor/normal whole-genome sequencing (WGS). In two families a consistent IGHV gene mutational status (one IGHV-unmutated, one IGHV-mutated) was observed. Intriguingly, the third family with four affected siblings was characterized by usage of the lambda IGLV3-21 gene, with the hallmark R110 mutation of the recently described clinically aggressive IGLV3-21R110 subset. In this family, the CLL-specific rearrangements in two siblings could be assigned to either stereotyped subset #2 or the immunogenetically related subset #169, both of which belong to the broader IGLV3-21R110 subgroup. Consistent patterns of cytogenetic aberrations were encountered in all three families. Furthermore, the CLL clones carried somatic mutations previously associated with IGHV mutational status, cytogenetic aberrations and stereotyped subsets, respectively. From these findings, we conclude that similarities in immunogenetic characteristics in familial CLL, in combination with genetic aberrations acquired, point towards shared underlying mechanisms behind CLL development within each family.

Prognostic and Predictive Implications of Cytogenetics and Genomics.

Chronic lymphocytic leukemia (CLL) is characterized by extreme genomic heterogeneity. Numerous recurrent genetic abnormalities are associated with dismal clinical outcome in patients treated with chemo(immuno)therapy, with aberrations of the TP53 gene being the main genomic abnormalities that dictate treatment choice. In the era of novel agents the predictive significance of the genomic aberrations is highly challenged as the results of the clinical trials performed thus far question the previously established unfavorable impact of genomic aberrations, even that of the TP53 gene. The prognostic and predictive value of the most common genomic abnormalities is discussed in the present review.

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.

Targeted sequencing reveals the somatic mutation landscape in a Swedish breast cancer cohort.

Breast cancer (BC) is a genetically heterogeneous disease with high prevalence in Northern Europe. However, there has been no detailed investigation into the Scandinavian somatic landscape. Here, in a homogeneous Swedish cohort, we describe the somatic events underlying BC, leveraging a targeted next-generation sequencing approach. We designed a 20.5 Mb array targeting coding and regulatory regions of genes with a known role in BC (n = 765). The selected genes were either from human BC studies (n = 294) or from within canine mammary tumor associated regions (n = 471). A set of predominantly estrogen receptor positive tumors (ER + 85%) and their normal tissue counterparts (n = 61) were sequenced to ~ 140 × and 85 × mean target coverage, respectively. MuTect2 and VarScan2 were employed to detect single nucleotide variants (SNVs) and copy number aberrations (CNAs), while MutSigCV (SNVs) and GISTIC (CNAs) algorithms estimated the significance of recurrent somatic events. The significantly mutated genes (q ≤ 0.01) were PIK3CA (28% of patients), TP53 (21%) and CDH1 (11%). However, histone modifying genes contained the largest number of variants (KMT2C and ARID1A, together 28%). Mutations in KMT2C were mutually exclusive with PI3KCA mutations (p ≤ 0. 001) and half of these affect the formation of a functional PHD domain. The tumor suppressor CDK10 was deleted in 80% of the cohort while the oncogene MDM4 was amplified. Mutational signature analyses pointed towards APOBEC deaminase activity (COSMIC signature 2) and DNA mismatch repair (COSMIC signature 6). We noticed two significantly distinct patterns related to patient age; TP53 being more mutated in the younger group (29% vs 9% of patients) and CDH23 mutations were absent from the older group. The increased somatic mutation prevalence in the histone modifying genes KMT2C and ARID1A distinguishes the Swedish cohort from previous studies. KMT2C regulates enhancer activation and assists tumor proliferation in a hormone-rich environment, possibly pointing to a role in ER + BC, especially in older cases. Finally, age of onset appears to affect the mutational landscape suggesting that a larger age-diverse population incorporating more molecular subtypes should be studied to elucidate the underlying mechanisms.

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.

Defining eligible patients for allele-selective chemotherapies targeting NAT2 in colorectal cancer.

Therapies targeting somatic bystander genetic events represent a new avenue for cancer treatment. We recently identified a subset of colorectal cancer (CRC) patients who are heterozygous for a wild-type and a low activity allele (NAT2*6) but lack the wild-type allele in their tumors due to loss of heterozygosity (LOH) at 8p22. These tumors were sensitive to treatment with a cytotoxic substrate of NAT2 (6-(4-aminophenyl)-N-(3,4,5-trimethoxyphenyl)pyrazin-2-amine, APA), and pointed to NAT2 loss being a therapeutically exploitable vulnerability of CRC tumors. To better estimate the total number of treatable CRC patients, we here determined whether tumor cells retaining also other NAT2 low activity variants after LOH respond to APA treatment. The prevalent low activity alleles NAT2*5 and NAT2*14, but not NAT2*7, were found to be low metabolizers with high sensitivity to APA. By analysis of two different CRC patient cohorts, we detected heterozygosity for NAT2 alleles targetable by APA, along with allelic imbalances pointing to LOH, in ~ 24% of tumors. Finally, to haplotype the NAT2 locus in tumor and patient-matched normal samples in a clinical setting, we develop and demonstrate a long-read sequencing based assay. In total, > 79.000 CRC patients per year fulfil genetic criteria for high sensitivity to a NAT2 LOH therapy and their eligibility can be assessed by clinical sequencing.