Whole Exome Sequencing reveals NOTCH1 mutations in anaplastic large cell lymphoma and points to Notch both as a key pathway and a potential therapeutic target

Abstract

Patients diagnosed with Anaplastic Large Cell Lymphoma (ALCL) are still treated with toxic multi-agent chemotherapy and as many as 25-50% of patients relapse. To understand disease pathology and to uncover novel targets for therapy, Whole-Exome Sequencing (WES) of Anaplastic Lymphoma Kinase (ALK)+ ALCL was performed as well as Gene-Set Enrichment Analysis. This revealed that the T-cell receptor (TCR) and Notch pathways were the most enriched in mutations. In particular, variant T349P of NOTCH1, which confers a growth advantage to cells in which it is expressed, was detected in 12% of ALK+ and ALK- ALCL patient samples. Furthermore, we demonstrate that NPM-ALK promotes NOTCH1 expression through binding of STAT3 upstream of NOTCH1. Moreover, inhibition of NOTCH1 with γ-secretase inhibitors (GSIs) or silencing by shRNA leads to apoptosis; co-treatment in vitro with the ALK inhibitor Crizotinib led to additive/synergistic anti-tumour activity suggesting this may be an appropriate combination therapy for future use in the circumvention of ALK inhibitor resistance. Indeed, Crizotinib-resistant and sensitive ALCL were equally sensitive to GSIs. In conclusion, we show a variant in the extracellular domain of NOTCH1 that provides a growth advantage to cells and confirm the suitability of the Notch pathway as a second-line druggable target in ALK+ ALCL.

Introduction

Systemic anaplastic large cell lymphoma (ALCL) is a Tcell malignancy accounting for approximately 15% of all cases of pediatric lymphoma and 1-2% of adult lymphomas.¹ The majority of pediatric cases (>90%) express NPM-ALK, the result of the t(2;5)(p23;q35) translocation, creating a fusion between the nucleolar phosphoprotein gene nucleophosmin1 (NPM1) and anaplastic lymphoma kinase (ALK), leading to the ALK⁺ ALCL diagnostic entity, although other ALK fusion proteins also exist.^2,3 ALK fusion proteins induce the activation of several downstream signaling pathways involved in oncogenesis, including PI3K, ERK1/2 MAP kinase, and JAK-STAT.

Genetic studies conducted to date have not revealed consistent genetic abnormalities among ALK⁺ ALCL, although a higher frequency of genomic gains and losses has been associated with a poorer prognosis.^4-6 Genomic gains include the region encoding the NOTCH1 locus,⁴ which may explain why NOTCH1 is expressed in the majority of ALK⁺ ALCL and in some ALK^– ALCL.^7,8 In contrast, for ALK^– ALCL, mutations in the JAK1/STAT3 pathway have been described.⁹

ALK⁺ ALCL is largely a chemosensitive malignancy, although despite good initial responses, the relapse rate reaches 50% independently of the chemotherapy regimen used;^10-14 therefore, new therapies are needed for patients who do not respond to standard chemotherapy. There is also a need for less toxic treatment schedules for low-risk patients. Crizotinib, a small molecule ATP-competitive inhibitor of ALK/MET/ROS1 which is currently in clinical trials (NCT01606878, NCT02034981) for pediatric ALK⁺ lymphoid malignancies, given as monotherapy produces remission in more than 80% of relapsed patients, although rapid relapse on discontinuation of therapy has been reported.^11,15-17 As such, second-line treatment and combination therapies for relapsed patients are still required. We therefore conducted whole exome sequencing (WES) of 25 ALK⁺ ALCL tumors, validating data in a further 78 cases of ALCL to understand disease pathology and to uncover novel targets for therapy. Beside the T-cell receptor (TCR) pathway, the most commonly affected signaling axis is that involving NOTCH1, of which a T349P variant provides a growth advantage to cells. Furthermore, NOTCH1 is expressed as a consequence of NPM-ALK-driven STAT3 activity, a key signaling node in ALCL. Accordingly, the NOTCH1 pathway provides a target for second-line therapy, whereby γ-secretase inhibitors (GSI) show synergistic activity with inhibitors of ALK, and are efficacious as single agents in ALK-inhibitor resistant cell lines. Finally, we show evidence that NOTCH1 is a biomarker predictive of relapse risk.

Methods

Patients’ samples

Patients’ tumor tissues (at initial presentation) and matched peripheral blood DNA were obtained following patient/parental consent according to both the Declaration of Helsinki and local guidelines from the following institutions: Children’s Cancer and Leukaemia Group tissue bank, Newcastle, UK (Online Supplementary Table S1); Institut Universitaire du Cancer Toulouse, France; University Hospital Brno, Czech Republic; Biobank of the Medical University of Graz, Austria; Belarusian Center for Pediatric Oncology, Hematology and Immunology, Minsk, Belarus; Justus-Liebig University, Giessen, Germany; and Our Lady’s Children’s Hospital, Crumlin, Ireland (Online Supplementary Table S2). All tissues were obtained and processed with full ethical approval (NHS Research ethics committee reference numbers 07/Q0104/16, 06/MRE04/90 and 08/H0405/22+5).

The patient tissue microarray used here has been described elsewhere.¹⁸ Briefly, formalin-fixed paraffin-embedded (FFPE) tissue specimens from pediatric, patients with NPM-ALK⁺ ALCL treated in the Berlin-Frankfurt-Munster (BFM) group study NHL-BFM90, NHL-BFM95 or patients enrolled in the European intergroup trial ALCL99 between August 1998 and December 2008 were obtained from both male and female children (Online Supplementary Table S3) with informed consent and in accordance with the Declaration of Helsinki. Eligibility was confirmed by demonstration of NPMALK positivity of the tumor either by NPM-ALK polymerase chain reaction, two-color fluorescence in situ hybridization for the t(2;5), or nuclear and cytoplasmic staining for ALK. The inclusion criteria were fulfilled by 89 patients.

Whole exome sequencing

DNA was extracted from fresh-frozen tissue from patients (n=18; with a tumor content >90%) and matched peripheral blood from four patients using the QIAgen DNAeasy Blood and Tissue Kit (Qiagen, Hilden, Germany), following the manufacturer’s instructions. Library preparation was conducted using the Nextera Rapid Capture Exome Kit before samples were sequenced, at either the Washington State University Core (Illumina HiSeq2500) or the Molecular Diagnostics and Personalized Therapeutics Unit, University of Ha’il (Illumina MiSeq) (Online Supplementary Table S4). Sequencing data are available at the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra) under accession number PRJNA491296. The bioinformatic processing is detailed in the Online Supplementary Methods.

Immunohistochemistry, pathology and quantification

Immunohistochemistry was performed on FFPE sections with the conventional avidin–biotin–peroxidase method. Heat antigen retrieval was performed using citrate buffer, pH 6.1. Endogenous peroxidases were quenched by incubating sections in 3% H²O² in phosphate-buffered saline (PBS) for 10 min. An avidin/biotin and a protein block were subsequently performed. Primary antibodies (Online Supplementary Table S5) were added in 1% bovine serum albumin/PBS at 4°C overnight. Slides were incubated with biotinconjugated secondary antibodies and then with horse radish peroxidase (HRP) using the IDtect Super Stain System – HRP and developed under visual control using aminoethyl carbazole. Hemalaun counterstaining was performed and slides mounted with AquaTex. Sections were washed with PBS three times in between each step. Stained slides were assessed for cleaved NOTCH1 staining by an experienced pathologist (blinded with respect to clinicopathological parameters and patients’ outcome) using the histoscore system. Stained slides were scored qualitatively for the intensity of staining and classified as showing negative, weak, moderate or strong staining (to qualify for ‘moderate’ or ‘strong’ staining, at least 10% of cells had to stain positive). Analysis of event-free survival was performed as described previously,¹⁹ grouping negative and weak staining into ‘low cleaved NOTCH1 expression’ and moderate and strong into ‘high cleaved NOTCH1 expression’.

Statistical analyses

All experiments were executed in biological triplicates. The MTT, RealTimeGlo, apoptosis, cell cycle and quantitative polymerase chain reaction assays were additionally executed with technical triplicates. All plots are representative of the mean of the biological replicates, while the error bars represent the standard deviation. Two-tailed t-tests were used to calculate the P-value when comparing two samples (multiple comparisons were corrected using the Holm-Sidak method); when comparing more than two samples, two-way analysis of variance (ANOVA) was used (again, multiple comparisons were corrected using the Holm- Sidak method). Statistical tests were conducted using GraphPad PRISM 8 (Graphpad).

Results

The genomic profile of anaplastic large cell lymphoma

Eighteen ALK⁺ ALCL exomes sequenced in this study in addition to seven previously sequenced samples⁹ (Online Supplementary Table S1) were analyzed. The samples comprised 17 pediatric cases (≤18 years) and eight adult cases (Online Supplementary Table S4); a flowchart illustrating the cohorts of patients is shown in Online Supplementary Figure S1. All patient samples were collected at diagnosis. Data regarding variants found in at least a quarter of the patients are summarized in Figure 1, which shows that the most commonly mutated genes in both adult and pediatric cases were TYW1B, DEFB132 and KCNJ18 (the full list of variants can be found in Online Supplementary Table S6). None of the variants in these genes has been reported previously in hematologic malignancies and were not studied further here. Two patients presented with one mutation each in TP53 (COSMIC ID: COSM3958801 and COSM9969). We also studied copy number variations, but found no novel events larger than 100,000 bp present consistently in more than one sample at a sequencing depth of at least 50x (Online Supplementary Figure S2A, Online Supplementary Table S7), as previously observed.⁴ Among previously reported alterations in ALCL, a single copy gain on chromosome 7 was observed in three patient tumor samples (S3, S9 and S15)⁵ and a single copy loss on chromosome 17p was also seen in three patients (S9, S14 and S57).

The most predominant single nucleotide variants are non-synonymous and are present at higher levels in patients who subsequently relapsed

The majority of somatic variants detected in the 25 tumor samples were non-synonymous single nucleotide variants (39.4%), in keeping with a previous publication reporting that the ALK⁺ ALCL genome is largely stable.⁴ Single nucleotide variants were followed in frequency by frameshift and non-frameshift deletions and splice variants (24.1%, 10.8% and 10.6%, respectively), while the germline genome of ALCL patients points to an overwhelming presence of single nucleotide polymorphisms (89.3%) (Online Supplementary Figure S2B). The proportion of each type of variant detected differed between patient tumors (Figure 2B), although in general, pediatric patients known to have relapsed (n=9) had a significantly higher proportion of non-synonymous single nucleotide variants than patients who did not (n=9; P<0.0001) (Figure 2B), suggesting that a high percentage of non-synonymous single nucleotide variants at diagnosis may be indicative of relapse, although this requires validation in a larger dataset of patients treated with comparable treatment protocols.

Deficiency of DNA repair mechanisms and spontaneous deamination of 5-methyl cytosine are identified as signatures of anaplastic large cell lymphoma

Online Supplementary Figure S2C shows the prevalence, in representative patient S57, of the 96 variant types that were used to derive the mutational signatures (Online Supplementary Figure S2D). Examining the type of mutations in the patients for whom matched peripheral blood was available (n=11), showed an enrichment for mutational signatures 1, 3, 12 and 26²⁰ (Figure 2C). Interestingly, 1A is a signature based on the prevalence of C>T transitions at NpCpG trinucleotides and is associated with spontaneous deamination of 5-methyl-cytosine,²¹ whereas signature 3 has its roots in homologous recombination deficiency during DNA double-strand break repair.²⁰ The etiology of signature 12 has not yet been identified, although signature 26 is associated with a breakdown in DNA mismatch repair. The combination of signatures 3 and 26 may indicate, from an evolutionary perspective, how ALCL tumors accumulate mutations. Comparable patterns were found when comparing signatures to the COSMIC signature database²² (data not shown). There was no detectable difference between the mutational signature of pediatric (n=4) or adult (n=7) ALK⁺ ALCL patients (data not shown).

Gene set enrichment analysis confirms the importance of T-cell receptor signaling, but also of the Notch pathway in ALK+ anaplastic large cell lymphoma pathobiology

Gene set enrichment analysis (GSEA) of mutated genes showed that TCR signaling and Notch pathways are enriched across all five databases used (Figure 2D; Online Supplementary Table S8). Further analysis of the domains frequently found in the mutated genes revealed an enrichment in proteins with epidermal growth factor (EGF)-like or calcium-ion binding domains (Figure 2E), two features of the NOTCH1 protein, and indeed the locus of both of the NOTCH1 mutations identified in this study (see below). Twenty of the 25 patient tumors carry mutations in proteins of the Notch pathway with a range of one to four and a median of two mutations per patient (Online Supplementary Table S8). Furthermore, reactome network clustering analysis²³ showed TP53 as a key node, which is not unexpected as TP53 has been reported to play a key role in the pathogenesis of ALCL²⁴ (Online Supplementary Figure S2E). Given the importance of the Notch pathway in T-cell biology, particularly in the developing thymus, which we proposed tp be the origin of ALK⁺ ALCL,²⁵ and the previous implication of the NOTCH1 pathway in the pathogenesis of ALCL,⁸ the NOTCH1 mutations detected and the NOTCH1 pathway were explored for their role in the pathogenesis of ALK⁺ ALCL.

Figure 1.Summary of whole exome sequencing data. Individual results for each patient presented as a heatmap for genes found to be mutated in at least a quarter of patients. The number of different variants presenting in each gene is displayed for individual patients. The total number of variants identified for individual patients is also displayed above the graph. The panel on the right shows the percentage of patients presenting with at least one variant of the indicated gene. Pediatric and adult patients are separated. Patients for whom matched peripheral blood was sequenced are labeled with an ‘M’.

NOTCH1 variant T349P provides a growth advantage to cells

Of the 25 ALK⁺ ALCL tumor samples analyzed by WES, 24% presented with the NOTCH1 variant T349P, while 12% had the T311P variant. These data were validated by Sanger sequencing of a further 78 samples (including 18 of the samples previously analyzed by WES with a total of 55 ALK+ ALCL, and 23 ALK– ALCL) (Online Supplementary Table S2). In this validation cohort, the T349P variant was detected in 12% of patients (n=78; 15% of ALK^– patients and 9.3% of ALK⁺ patients) (Online Supplementary Figure S3A) and the T311P variant was found in 7.6% of patients (n=78; 10.2% of ALK⁺ patients, none in ALK^– patients) (Online Supplementary Figure S3B). In the majority of cases, tumors presented with a mutation at either T311P or T349P and therefore the overall incidence of patients with at least one mutation of the EGF-like domain was 18% (n=78), although one patient presented with both mutations. We detected two additional NOTCH1 mutations, H1190P and G1503S, although these were only detected in one patient each and so were not studied further. Of note, there was no significant difference in prognosis for patients presenting with NOTCH1 T349P and/or T311P over those with wild-type (WT) NOTCH1 when considering the whole cohort of patients (Online Supplementary Table S9, Online Supplementary Figure S3J). However, if adult patients were considered in isolation, there was a significant reduction in overall survival for those with NOTCH1 mutations (P<0.05) but these data are based on just three patients and so should be interpreted with caution.

Figure 2.Analysis of whole exome sequencing data yields insights into the ALK+ anaplastic large cell lymphoma genomic landscape. (A) Proportion of variant type for each patient. (B) Mean and standard deviation of the proportion of variant types for patients who did (n=9) or did not (n=9) relapse within 5 years after diagnosis (***P<0.001). (C) Mutational signatures²⁰ derived from the variant profiles of each of the 11 patients for whom we hold sequenced matched peripheral blood. (D, E) Scatter plots of the pathways (D) or domains (E) found to be enriched in our dataset, displaying the number of databases in which each hit was found to be enriched, along with the (-)log10 of the statistical enrichment P-value of the software in which each hit was found to be most enriched, and the corresponding number of genes involved. SNV: single nucleotide variant.

Variant T349P, at position 1045, is within the sixth of 34 exons of NOTCH1, which encodes one of the numerous EGF-like domains that make up the extracellular domain of NOTCH1 (EGF-like domain 9 of 36, which is a calcium-binding domain). NOTCH1 T349P was predicted to be a function-altering mutation by variant effect prediction software including SIFT (score=0.01) and PolyPhen (score=0.999), among others.^26-30 Furthermore, the COSMIC database shows that T311 and T349 are the two most frequently reported mutated amino acids at the presumed NOTCH1/JAG1 interface across a range of cancers (including chronic myelomonocytic leukemia,³¹ chronic lymphocytic leukemia,³² T-cell acute lymphoblastic leukemia,³³ rhabdomyoscaroma³⁴ and squamous cell carcinoma³⁵) (Figure 3A).

To determine the impact of the detected NOTCH1 mutants on cell proliferation, WT, T349P or T311P mutants of NOTCH1 were expressed in HEK293FT cells (Figure 3B and C). These cells were chosen because of their low levels of endogenous NOTCH1 expression. In comparison, ALCL cell lines express high levels of WT NOTCH1. A significant increase in ATP production, suggestive of enhanced cell proliferation (as determined by a RealTime-Glo assay), was observed at 72 h for cells expressing the NOTCH1 T349P mutant as compared to WT NOTCH1 in the absence of exogenously applied ligand (Figure 3D). The T311P mutant also led to an increase in proliferation, although this was not statistically significant when compared to WT NOTCH1. However, a significant increase in proliferation was detected when using an MTT assay (Online Supplementary Figure S3I). In addition, the transcriptional targets of NOTCH1 activity, HES1 and HEY1, were expressed at higher levels in the HEK293FT cells expressing NOTCH1 mutants compared to those transfected with either an empty vector or WT NOTCH1 (Online Supplementary Figure S3F and S3G). These data are suggestive of increased NOTCH1 activity as a consequence of the NOTCH1 T349P mutation compared to WT NOTCH1.

Figure 3.NOTCH1 T349P and T311P mutants confer a growth advantage to cells. (A) The most frequent variants (of the amino acids involved in the NOTCH1-JAG1 interface) of NOTCH1 reported on COSMIC, and the frequency at which they have been reported. (B) Quantitative polymerase chain reaction (qPCR) showing the foldincrease in transcript levels of NOTCH1 compared to that in cells transduced with an empty vector (n=3). (C) Western blot for full-length and cleaved intracellular NOTCH1 (ICN) and α-tubulin in HEK293FT cells 48 h after transfection to express wild-type (WT) or mutant (T349P or T311P) NOTCH1, or an empty vector (EV). Only the relevant sections of the whole blot are shown and the contrast of the whole image was modified in order to improve legibility. Data are representative of three biological repeats. (D) Proliferation of HEK293FT cells expressing the indicated proteins or empty vector control as measured by the Real Time Glo assay at 24, 48 and 72 h after transfection (***P<0.001; n=3). (E) Fold-change in expression of DLL1 over the WT NOTCH1 control as assessed by qPCR in HEK293FT cells expressing the indicated proteins and co-cultured with either WT or DLL1-expressing OP9 cells (*P<0.05; **P<0.01; n=3). All bar plots display the mean of biological replicates, and error bars represent the standard deviation.

To determine why the T349P NOTCH1 mutant led to enhanced cell proliferation, in silico modeling was performed to predict the effects of the mutation on protein conformation using a published crystal structure of NOTCH1 bound to one of its ligands (Jagged1; PDBID 5UK5) (Online Supplementary Figure S3C and S3D).³⁶ The model demonstrated that NOTCH1 residues 349 and 311 mediate binding to NOTCH1 ligands. Therefore, to determine whether increased cell proliferation (Figure 3D) was dependent on NOTCH1 ligand, NOTCH1 WT or mutant-expressing HEK293FT cells were co-cultured with OP9 cells expressing or not the NOTCH1 ligand DLL1. There was no discernible difference in proliferation between cells co-cultured with wild-type OP9, or OP9-DLL1 cells (Online Supplementary Figure S3E). However, increased transcription of endogenous DLL1 was observed on expression of the mutant forms of NOTCH1 compared to WT protein, suggesting that the mutant proteins might themselves lead to transcription of ligand in an autonomous fashion (Figure 3E). In an attempt to validate this, we silenced DLL1 using a specific siRNA (Online Supplementary Figure S3H), which reduced the proliferation advantage induced by the NOTCH1 T349P mutant (Online Supplementary Figure S3I).

NPM-ALK induces expression of NOTCH1 via STAT3 transcriptional activity in anaplastic large cell lymphoma

Given that ALCL cell lines express high levels of WT NOTCH1,⁸ the mechanism of NOTCH1 expression was investigated. NPM-ALK activity was inhibited by incubation with the ALK/ROS/MET inhibitor crizotinib (Figure 4A) or expression silenced with a specific shRNA (Figure 4D) in three and two cell lines, respectively. In both cases, a significant decrease in transcripts for NOTCH1 (Figure 4B and E) and its transcriptional target HES1 (Figure 4C and F), was observed, suggesting that NOTCH1 transcription and NOTCH1 activity are dependent on NPM-ALK. Given that STAT3 is a key nodal downstream target of NPM-ALK, STAT3 expression was inhibited by employing specific shRNA in ALCL cell lines (Figure 4G). As predicted, transcript levels of NOTCH1 (Figure 4H), HES1 (Figure 4I) and HEY1 (Figure 4J) were all significantly downregulated as a result of STAT3 silencing in both ALK⁺ and ALK^– ALCL cell lines (Figure 4H-J, Online Supplementary Figure S4A). Analysis of published chromatin immunoprecipitation sequencing data³⁷ of STAT3 binding sites in the ALCL cell lines SU-DHL1 and JB6, treated with either crizotinib or a vehicle control, showed a significant decrease in binding of STAT3 at the NOTCH1 gene in crizotinib-treated cells (Figure 4K, Online Supplementary Figure S4B). These data were validated by chromatin immunoprecipitation–quantitative polymerase chain reaction analysis whereby a significant decrease in binding of STAT3 at the NOTCH1 gene upon crizotinib treatment in the SUP-M2 cell lines was observed (Figure 4L).

Figure 4.STAT3-mediated regulation of NOTCH1 by NPM-ALK in ALK⁺ anaplastic large cell lymphoma. (A) Representative western blot for ALK, phospho-ALK and α- tubulin in ALK+ anaplastic large cell lymphoma (ALCL) cell lines when treated with 300 nM crizotinib or a vehicle control (dimethylsulfoxide, DMSO) for 6 h. Only the relevant sections of the whole blot are shown and the contrast of the whole image was modified in order to to improve legibility. Data are representative of three biological replicates. Densitometry is included, as fold-change over the vehicle control and loading control. (B, C) Fold-change expression of NOTCH1 (B) and HES1 (C) over vehicle control in the indicated ALK⁺ ALCL cell lines 48 h after treatment with crizotinib, as determined by quantitative polymerase chain reaction (qPCR) (*P<0.05; **P<0.01; ***P<0.001; n=3). (D-F) Fold-change expression of NPM-ALK (D), NOTCH1 (E) and HES1 (F) over non-targeting control in ALK⁺ ALCL cell lines 48 h after transduction with control non-targeting (NT) shRNA, or a shRNA targeting NPM-ALK, as determined by qPCR (**P<0.01; ***P<0.001; n=3). (G-J) Foldchange expression of STAT3 (G), NOTCH1 (H), HES1 (I) and HEY1 (J) over non-targeting NT control in ALK⁺ ALCL cell lines 48 h after transduction with control NT shRNA, or one of two shRNA targeting STAT3, as determined by qPCR (**P<0.01; ***P<0.001; n=3). (K) Binding of STAT3 to the promoter regions of NOTCH1 in SUDHL1 cells treated with a vehicle control (middle track) or crizotinib (lower track); the upper track is the input for two separate cell lines, data were obtained by analyzing previously published data.³⁷ (L) Chromatin immunoprecipitation-qPCR binding of STAT3 and GFP at the NOTCH1 promoter region, or at a negative control intergenic region, in SUP-M2 cells treated with either a vehicle control or crizotinib (300 nM) for 6 h, as determined by qPCR (***P<0.0001; n=3), expressed as the percentage of the total input. All measures of expression by qPCR were normalized to GAPDH and PPIA. All bar plots display the mean of biological replicates, and error bars represent standard deviations.

Figure 5.Silencing NOTCH1 expression in anaplastic large cell lymphoma cell lines inhibits cell proliferation and induces cell death. (A) Fold change expression of NOTCH1 over non-targeting (NT) control (normalized to GAPDH and PPIA) in the indicated ALK⁺ anaplastic large cell lymphoma (ALCL) cell lines 48 h after transduction with control non-targeting (NT) shRNA, or one of two shRNA targeting NOTCH1 as determined by quantitative polymerase chain reaction (**P<0.01; ***P<0.001; n=3). (B) Representative western blot for cleaved intracellular NOTCH1 (ICN) and α-tubulin in ALK⁺ ALCL cell lines 48 h after transduction with NT control shRNA, or one of three shRNA targeting NOTCH1. Only the relevant sections of the whole blot are shown and the contrast of the whole image was modified in order to improve legibility. Data are representative of three biological repeats. Densitometry is included as fold-change over the vehicle control and loading control. (C) Proliferation of ALK⁺ ALCL cell lines over the NT control shRNA, determined using an MTT assay 48 h after transduction with NT control shRNA, or one of two shRNA targeting NOTCH1 (*P<0.05; **P<0.01; ***P<0.001; n=3). (D) Quantification of the percentage of cells positive for annexin V and/or propidium iodide 48 h after transduction with NT control shRNA (left panel), or two shRNA targeting NOTCH1 (middle and right panels) (*P<0.05; ***P<0.001; n=3). (E-H) Fold-change of expression of HES1 (E), HEY1 (F), DTX1 (G) or MYC (H) over NT control shRNA transduced cells (normalized to GAPDH and PPIA) in the indicated ALK⁺ ALCL cell lines 48 h after transduction with NT control shRNA, or one of two shRNA targeting NOTCH1, as determined by quantitative polymerase chain reaction (*P<0.05; **P<0.01; ***P<0.001; n=3). All bar plots display the mean of biological replicates, and error bars represent standard deviations; the bar plots are color-coded as indicated in the Figure.

NOTCH1 is a therapeutic target in ALK– and ALK+ anaplastic large cell lymphoma

Given the expression and activity of NOTCH1 in ALCL, its utility as a therapeutic target was investigated. Expression of NOTCH1 was inhibited at both transcript (Figure 5A, Online Supplementary Figure S5A) and protein levels (Figure 5B) by shRNA in both ALK⁺ and ALK^– ALCL cell lines leading to a significant decrease in cell proliferation at 48 h, as measured by the MTT assay, in all four cell lines (Figure 5C, Online Supplementary Figure S5B) concomitant with a significant increase in the percentage of cells staining positive for either annexin-V or propidium iodide, suggestive of cell death by apoptosis (Figure 5D). As expected, NOTCH1’s transcriptional targets HES1 and HEY1 (Figures 5E and 5F, Online Supplementary Figure S5A) were also downregulated upon NOTCH1 silencing by shRNA. To further determine the signaling cascade downstream of NOTCH1, we made use of published microarray data examining the effect of GSI treatment on gene expression in T-cell acute lymphoblastic leukemia,^38-40 and found genes such as MYC and DTX1 to be potential targets of NOTCH1 signaling (Online Supplementary Figure S5C). In keeping with these data, we show that silencing NOTCH1 by shRNA in the ALK⁺ ALCL cell lines DEL, SU-DHL1 and SUP-M2 (and the ALK^– ALCL cell line FEPD) leads to significant decreases in both MYC and DTX1 transcript levels (Figure 5G and H, Online Supplementary Figure S5A) and protein levels (Online Supplementary Figure S5D), suggesting that NOTCH1 in ALCL signals through a number of pathways beyond HES1 and HEY1. Indeed, analysis of published microarray data suggests that the expression of MYC and NOTCH1, and DTX1 and NOTCH1 correlates in both ALK⁺ and ALK^– ALCL, but not in reactive lymph nodes (Online Supplementary Figure S5E-J).

Figure 6.ALK inhibitor sensitive and resistant anaplastic large cell lymphoma cell lines are sensitive to γ-secretase inhibitors. (A) Western blot for cleaved intracellular NOTCH1 (ICN) and α-tubulin in ALK⁺ anaplastic large cell lymphoma (ALCL) cell lines when treated with 1 mM GSI-I for 48 h. Only the relevant sections of the whole blot are shown and the contrast of the whole image was modified in order to improve legibility. Data are representative of three biological repeats. (B) Proliferation of a panel of ALK⁺ and ALK^– ALCL cell lines treated with 1 mM of GSI-I for 48 h, compared to vehicle control, as measured by the MTT assay (***P<0.0001; n=3). (C) Quantification of the percentage of cells positive for annexin V (AV) and/or propidium iodide (PI) when treated with either vehicle control or 1 mM GSI-I for 48 h (***P<0.001; n=3). (D) BLISS matrix showing the combination index on treating the indicated ALK⁺ ALCL cell lines with crizotinib and PF- 03084014 for 72 h (using a range of concentrations from 25 to 100 nM for crizotinib, and from 100 nM to 10 mM for PF-03084014). A combination index of <1 indicates synergy between drugs, 1 indicates additive effects, >1 indicates antagonistic effects (n=3). (E) Quantification of the percentage of cells positive for AV and/or PI when treated with either vehicle control, 50 nM crizotinib, 2 mM PF-03084014 or a combination of PF-03084014 and crizotinib for 48 h (NS: not significant; *P<0.05; ***P<0.00; n=3). (F) Proliferation over vehicle control of wild-type or crizotinib-resistant Karpas-299 cells when treated with increasing concentrations of GSI-I, as measured by the MTT assay (n=3). (G) Ten-year event-free survival of patients with ALK⁺ anaplastic large cell lymphoma showing either little or no (n=56), or strong (n=33) NOTCH1 expression. All bar plots display the mean of biological replicates and error bars represent standard deviations.

Given the reduced cell proliferation and increased cell death observed on shRNA-mediated knockdown of NOTCH1 expression, ALCL cell lines were incubated with two different GSI: GSI-I (Z-LLNle-CHO) and PF- 03084014.⁴¹ GSI-I inhibited NOTCH1 cleavage, as shown by a decrease in expression of the intracellular domain of NOTCH1 (ICN) 48 h after drug treatment (Figure 6A), concomitant with a decrease in cell proliferation, as shown by the MTT assay (Figure 6B), and an increase in apoptosis, as determined by positive cell surface staining for annexin V and/or intracellular propidium iodide (Figure 6C, Online Supplementary Figure S6A). In contrast, PF-03084014 did not significantly affect cell proliferation or apoptosis on its own at any of the concentrations assessed, ranging from 10 nM to 10 mM for up to 72 h of incubation (data not shown). None of the ALCL cell lines used in the research described here carries either the NOTCH1 T311P or T349P mutation.

Gamma secretase inhibitors synergize with ALK inhibitors to induce cell death

Co-incubation of three of four ALCL cell lines with either PF-03084014 or GSI-1 with crizotinib led to additive to synergistic activity in reducing cell proliferation, as indicated by a Bliss Independence Index of less than one across several concentrations (Figure 6D, Online Supplementary Figure S6B). Karpas 299 cells have different genetic defects, which may explain the antagonistic results in this cell line.⁴² Indeed, a significant increase in apoptotic cells was observed following 48 h of treatment with a combination of 50 nM crizotinib and 2 mM PF- 03084014 (Figure 6E, Online Supplementary Figure S6C).

Ideally, a single-agent ALK inhibitor would provide a less toxic frontline treatment approach in the future, although resistance would be expected to develop. Therefore, crizotinib-resistant ALCL cell lines were assessed for their sensitivity to GSI which led to a significant decrease in cell proliferation (Figure 6F, Online Supplementary Figure S6D-F). These data suggest the potential use of GSI as either second-line treatment for ALK inhibitor-resistant disease, or frontline therapy in combination with ALK tyrosine kinase inhibitors.

NOTCH1 is a potential biomarker predictive of relapse in ALK+ anaplastic large cell lymphoma

A clinically annotated FFPE tissue microarray of 89 ALK⁺ ALCL patient samples, biopsied at the time of the patients’ initial presentation, was analyzed for cleaved NOTCH1 protein expression (Online Supplementary Figure S6G). Of the 89 patient tumors assessed, 88.8% showed high cleaved NOTCH1 staining (moderate and strong staining categories), in keeping with previously published results.^7,8,43 Interestingly, patients with low cleaved NOTCH1 staining (negative and weak staining categories) are more likely to relapse (based on the 10-year event-free survival) and have a significantly worse prognosis (P<0.05) (Figure 6G).

Discussion

The genetics underlying ALK⁺ ALCL at the level of somatic mutations remains largely unknown. Being a relatively rare cancer with a common, well-characterized driving oncogenic event, more could be done to uncover other genetic alterations. We therefore sequenced 18 ALK⁺ ALCL tumors and analyzed their genome together with seven previously reported ALK⁺ cases⁹ in order to uncover pathogenic mechanisms and novel therapeutic targets. Mutational signature analysis showed signatures 1, 3 and 26 in all our patient samples. The latter two signatures have their roots in homologous recombination DNA double- strand-break repair and mismatch DNA repair deficiency.²⁰ This suggests that DNA damage repair mechanisms might be impaired in these patients, predisposing them to ALCL, perhaps through germline mutations in DNA repair proteins (such as a BRCA2 variant K3326X, COSM4984873, found in germline sample 67B⁴⁴). In contrast, signature 1, accounting for the majority of the total contribution, is associated with a prevalence of C>T transitions at NpCpG trinucleotides and spontaneous deamination of 5-methyl-cytosine, considered an age-related phenomenon due to endogenous mutational processes.²¹ In general, 5-methyl cytosine residues are unstable within DNA and are prone to mutation, representing hotspots for this activity.⁴⁵ This is paradoxical with the young age range of ALK⁺ ALCL patients, with the majority being young adults, and suggests that the mechanism(s) leading to such mutations may be ‘speeded up’. Indeed, CpG transitions may be a consequence of secondary factors that promote deamination, such as exogenous mutagens and carcinogens, for example polycyclic aromatic hydrocarbons.⁴⁵

To identify pathways that are key to ALCL biology, GSEA was employed and revealed a number of pathways commonly affected by mutations in ALK⁺ ALCL. Interestingly, the TCR signaling pathway was prominent in our analysis. It has previously been shown that NPMALK can substitute for key TCR-induced distal signaling pathways and silencing of proximal proteins has been shown in ALCL.⁴⁶ Another key pathway identified was NOTCH1; gain-of-function mutations in NOTCH1 have previously been identified in a number of other cancers, most notably in approximately 50 to 60% of T-cell acute lymphoblastic leukemia.^47,48 However, most of these mutations are in the intracellular domains of the protein, with few reported in the extracellular domains.^49-51 In contrast, novel mutations in the EGF-like domain of extracellular NOTCH1 were detected in 9.3% (T349P) and 10.2% (T311P) of ALK⁺ ALCL patients analyzed in this study. Predicted to be deleterious, the functional significance of these mutations was investigated using bioinformatics analysis. It has previously been shown that EGF-like domains 8 to 12 are important for NOTCH1 binding to its ligands.^36,52 Specifically, threonine bases, lost in T349P and T311P, within the EGF-like domains are post-translationally modified by O-linked glycosylation which is necessary for ligand-engaged NOTCH1 signaling.^53-55 Mutations to proline, a rigid and bulky amino acid, result in a change in tertiary structure, often forcing a change in b-sheet conformation (EGF-like domains are made of b-sheets among others). We demonstrated the positive impact of the T349P mutant on NOTCH1 activity, as shown by enhanced cell proliferation when expressed in HEK293 cells. We theorize that NOTCH1 T349P could modulate ligand binding (either directly or through modulation of calcium binding, particularly as calcium signaling is thought to be dysregulated in ALCL and calcium ions play an important role in NOTCH1 ligand binding^56,57).

Regardless of the presence of NOTCH1 mutations, NOTCH1 constitutes a therapeutic target in ALCL independently of ALK status; suppression of NOTCH1 expression or activity led to an increase in apoptosis, in keeping with previous reports.^7,8,58 Intriguingly, NPM-ALK has previously been shown to be sufficient to induce NOTCH1 expression.²⁵ Not only did we confirm this by silencing NPM-ALK in ALCL, but we also showed that NPM-ALK acts through STAT3, which binds to the NOTCH1 promoter. This could explain why we observed synergistic effects between crizotinib and GSI in inducing cell death despite our evidence that NPM-ALK promotes NOTCH1 expression. Indeed, studies have shown that crizotinib synergizes with brentuximab vedotin, which targets CD30-expressing cells, despite NPM-ALK having been shown to drive CD30 expression.⁵⁹ ALK inhibitors are now being added to frontline therapy⁶⁰ (e.g., in trial NCT01979536) although this combination has led to some unexpected toxicity.¹¹ NOTCH1 inhibition may therefore serve as a second-line treatment. Indeed, although GSI have suffered from gastrointestinal-related toxicity, isoform-specific GSI or antibody-based treatments that target NOTCH1 directly (NCT03422679) have shown more promise.⁶¹ Published studies variously describe that GSI-1 (both as a single treatment and in combination) and PF-03084014 are relatively well tolerated.^62-64

Ultimately, our data show that GSI and ALK inhibitors act additively/synergistically and induce apoptosis of ALCL cell lines, and furthermore that ALK inhibitor-resistant cell lines remain sensitive to NOTCH1 inhibition. Hence, inhibition of NOTCH1 via GSI might represent a therapeutic option for both treatment-naïve and ALK inhibitor-resistant ALCL – while NOTCH1 expression may be a viable biomarker predictive of relapse.

Footnotes

• Received September 19, 2019
• Accepted April 9, 2020

Correspondence

SUZANNE D. TURNER sdt36@cam.ac.uk

Disclosures

No conflicts of interest to disclose.

Contributions

HL: conceptualization, methodology, investigations, writing of the original draft, and visualization. NP: conceptualization, investigations, writing, review and editing. JDM: software, visualization, formal analysis, writing, review and editing. MS, SH, EN, MFM, AE, MC, LM, CGP, GH, CL, SP, AJ, AF, AM, OW, GI, LL, OM, WK, SM, IA, MW, GAAB and LK: resources. AFA: investigations, software. SPD: investigation, writing, review and editing. WW and CDW: resources, writing, review and editing. MZ: software and formal analysis. TLB: conceptualization, methodology, software, writing, review and editing. SDT: conceptualization, resources, writing, review, editing, supervision, project administration, and acquisition of funding.