Citations

Abstract Activation of NFkB has been noted in late B-cell malignancies, but no genetic basis for this activation has been shown except for rare mutations of the NFkB pathway in lymphoma and none have been described in MM. Bortezomib, highly tumoricidal in a subset of MM patients, inhibits NFkB suggesting it’s role as a therapeutic target. Here we show for the first time the critical importance of mutations that dysregulate this potential survival pathway in human myeloma cell lines (HMCL) and new diagnosis MM patients. We performed high-density oligonucleotide array CGH (Agilent 44k) on CD138+ selected MM cells from 68 patients and 42 unique HMCL. Initially we focused on small loci, bi-allelically deleted, as markers for inactivation of potential tumor suppressor genes. We next developed FISH probes for the identified loci to confirm the bi-allelic deletions seen by aCGH. We next sequenced all samples with monoallelic deletions, and additional unselected patients, and identified additional inactivating mutations. Multiple novel tumor suppressor genes (TSG) were identified and validated, whose loss of function (both in patients and HMCL) leads to activation of the non-canonical NFkB pathway including: TRAF2 (9q34), TRAF3 (14q32), BIRC2/BIRC3 (11q22) and CYLD (16q12). Functional studies showed that HMCL with deletion or mutations in these genes have markedly increased NFkB2 p52/p100 ratios relative to other HMCL (and for BIRC2/BIRC3 relative to an isogenic line that retains one copy of BIRC2/BIRC3), indicating genetic activation of the non-canonical NFkB pathway. Reintroduction of deleted TSG into HMCL inhibited cell growth and resulted in cell death (see data on abstract by J. Keats et al). We next focused our analysis (of combined aCGH and GEP data) on positive regulators of this pathway. Rare IgH translocations of NIK have recently been identified in MM. By aCGH we identified amplification and rearrangements of NIK. We identified dysregulated expression of TNF family receptors that activate the non-canonical pathway: CD40 (translocation), TACI (amplification) and LTBR. All of the HMCL with these mutations also showed increased NFkB2 p52/p100 ratios, consistent with activation of the non-canonical NFkB. These mutations correlate with a gene expression profile of NFkB activation that is detected in one half of patients. In summary, we have determined for the first time the genetic basis for constitutive NFkB activation - present in 19/42 HMCL, and we estimate in at least one quarter of newly diagnosed patients (approximately one half of those with evidence of NFkB activation by GEP). We propose a mechanism in which positive (TACI, CD40, LTBR) and negative regulators (TRAF2, TRAF3, BIRC2/BIRC3, CYLD) of the pathway converge on NIK to mediate the processing of NFKB2 p100 to p52. The preponderance of genetic evidence suggests that drugs that target the non-canonical NFkB pathway, will be most effective in treating the 50% of MM patients with NFkB activation.

Abstract Introduction: The analysis of what goes awry in a malignant cell has focused primarily on mutations of protein-encoding genes and their regulatory sequences. However, recent work on microRNA (miRNA) has shed light on the possible involvement of miRNA genes in human disease. miRNA are small “noncoding” or “non-messenger” RNAs of about 21–25 nucleotides in length that function as regulators of gene expression essentially by pairing to the mRNA of protein-coding genes to initiate mRNA degradation or repression of translation. We are interested in characterizing miRNA expression profiles in human myeloma cell lines (HMCLs) and myeloma patient samples. Multiple myeloma (MM) is a plasma cell tumor characterized by frequent chromosomal translocations. Materials and methods: The expression of miRNA in HMCLs, primary MM cells (hyperdiploid samples) and normal plasma cells was determined using human microRNA chips. The GenoExporerTM Human microRNA chips were developed by GenoSensor (Tempe, AZ) and the microarray contains 226 human microRNA sequences in addition to control sequences. The DNA oligo probes are synthesized and immobilized on the chips (1”x 3” standard glass microslide). The probes are designed based on the active mature miRNA sequences and some of their flanking sequences. Our analysis included seven myeloma patients (hyperdiploid MM), eight different HMCLs (MM1, 8226 with t(14;16); SKMM2, U266, INA-6 with t(11;14) and H929, KMS11, UTMC-2 with t(4;14)), and eight normal human plasma cells. The normal human plasma cells were collected at Mayo Clinic Rochester from bone samples of patients undergoing orthopedic surgery. For the assay 5–10 μg of total RNA per sample were used at a concentration of 1 μg/μL. The RNA is directly labeled with biotin and used as a target for the on-chip hybridization assays. A streptavidin-Alexa dye is used to stain the hybridized targets and the fluorescent signals are captured and analyzed. The gene signal intensities were normalized to tRNA signal intensity. Results: The primary MM cells and HMCLs displayed a distinctive miRNA expression profile compared with normal plasma cells. Using the GeneSpring 7 (Agilent Technologies, Palo Alto, CA) for gene expression analysis we have identified miRNA genes with significant variation in expression levels between tumor and normal samples including miR-125b, miR-133a, miR1, and miR-124a (p<0.01). Recently defined algorithms were used to identify putative targets for the miRNA genes of interest such as leukemia inhibitory factor (LIF) as targets of miR-125b and Stat3 and angiopoietin-1 precursor as targets of miR-124a (Lewis et al. 2003). Of note, miR-15 and miR-16, previously identified to be downregulated in CLL, were expressed at low levels in some MM patients and HMCLs, but not in the normal plasma cells. Further analysis will be carried out to validate the data from miRNA profiling using northern blot analysis or real-time RT-PCR to measure expression levels of miRNA genes and that of their target genes. Conclusion: Analysis of miRNA expression pattern (along with aCGH and gene expression profiling studies) will undoubtedly refine our understanding of the various genotypic subtypes of multiple myeloma.

Abstract H-MM affects about 40–50% of patients, and a similar proportion of MGUS/SMM. The genetic evolution of the disease and the clinical implications of the co-existence of 13q- or t(14q32) are unknown. We conducted a global survey of genetic abnormalities in H-MM in order to describe genetic evolution pathways. H-MM cases from the Mitelman cancer genetic database were combined with cases from the Mayo Clinic MM genetic database. All abnormalities (regional/whole chromosome gains/losses, breakpoints and translocations) were plotted against a 344-chromosome band profile created in a matrix format, and the number of genetic abnormalities per tumor (NAPT) logged. The abnormalities present in more than 5% of tumors were considered recurrent and each karyotype was scored for the presence or absence of each of the recurrent abnormalities. To investigate the temporal pattern of acquisition of genetic abnormalities, all tumors with a given abnormality were selected and the distributions of NAPT plotted. The modes of these distributions were used as a value of ‘lateness’, referred to as time of occurrence (TO). The impact of co-existence of 13q- and t(14q32) on progression free survival (PFS) and overall survival (OS) in H-MM was determined for MM patients entered into an Eastern Cooperative Oncology Group chemotherapy trial. Four hundred and sixty-nine H-MM patients were included in the karyotypic analysis. The median chromosome number was 53 (range 48–74) and the median NAPT was 10 (range 1–25). The distribution of NAPT revealed a normal distribution, suggesting an obligatory acquisition of a given number of genetic abnormalities. The most common abnormalities were trisomies of chromosomes 3,5,7,9,11,15,19 and 21 affecting a third or more of the patients. Deletions and abnormalities of even-numbered chromosomes occur in 20% or less of patients. None of the more than 300 specific translocations detected were recurrent. However, 5 translocation breakpoints were recurrent including 14q32, 8q24 and 11q13. TO analysis suggest that for H-MM, the trisomies are early genetic events, translocations occur midway and deletions tend to be later events. Amongst the trisomies, correlation analysis does not suggest a particular order of acquisition. Clustering analysis using principal component analysis suggests that in H-MM, IgH translocations are secondary events associated with disease evolution rather than defining unique biological entities. Consistent with this, survival of H-MM patients is not affected by 13q- status but is negatively impacted by the presence of t(14q32) especially those with unknown partners (p<0.05). The data presented suggest that H-MM fit a proposed aneuploid model of cancer due to auto-catalytic chromosome reassortment (Duesberg et al Cell Cycle2003:2:202–210). The presence of tumor specific trisomies and the apparent similar phenotype induced by various combinations of trisomies are consistent with this model.

Abstract Waldenström’s macroglobulinemia (WM) is an incurable low-grade lymphoproliferative disorder characterized by bone marrow (BM) infiltration of a clonal population of small B-lymphocytes, plasmacytoid lymphocytes and plasma cells that secrete monoclonal IgM antibody. Because it is difficult to obtain tumor metaphases for karyotype studies, few recurrent chromosomal abnormalities have been reported and the genetic basis of the disease remains poorly defined. Therefore, we performed a comprehensive high-resolution study to identify genomic abnormalities involved in WM pathogenesis. Fifty-seven WM patients were analyzed using a combination of array-based comparative genomic hybridization (aCGH)(N=42); gene expression profiling (GEP)(N=22), DNA sequencing (N=24) and cIgM-FISH (N=57). Agilent 244A and Affymetrix U133 platforms were used in the aCGH and GEP experiments, respectively. An NF-κB index was calculated from GEP data, as a surrogate marker of NF-κB transcription activity. NFKB2 sub cellular localization was determined by immunofluorescence staining. Overall, 83% of samples have aCGH defined chromosomal abnormalities, with a median of 3 abnormalities per patient (range 0 to 27). We identified 16 recurrent regions of copy number change found in >5% (3 or more patients); 10 deleted and 6 amplified regions. The most common abnormality was the entire or partial deletion of 6q, identified in 40% of cases. Four nonoverlapped minimal deleted regions were identified on 6q (MDR1 to MDR4), each of them being present in at least 14 of 17 patients with the abnormality. Gain of 6p was the second most common abnormality (17%) and its presence was always concomitant with 6q loss. Other recurrent deletions were 13q14 (10%) and 7q22, 8p, 11q22-q23, 11q23-q24 and 17p11-p13 (7% each). Copy gains were identified as partial or entire gains of chromosomes 18 (17%), 4 (12%) and 3 (10%), 8q (10%) and Xq27.1-q28 (10%). Three patients had either homozygous deletions or LOH with mutations affecting both alleles of TRAF3 on 14q32. TRAF3 inactivation was correlated with an increased NF-kB transcriptional signature and it was associated with activation of the non-canonical NFkB pathway. Additionally, we identified the inactivation of TNFAIP3, another negative regulator of the NF-kB pathways, in one of 24 patients. Finally, interstitial deletions identified in 4 patients at 13q14 identified a 1.1 Mb MDR, including MIRN15A and MIRN16-1. These findings confirm that focal 13q14 deletion is not restrict to B-CLL and that the MIRN15A/MIRN16-1 abnormalities might be common to several indolent B-cell diseases. Overall, we identified TRAF3 and TNFAIP3 inactivation in 5.3% and 2.4%, respectively. To our knowledge, this is the first study showing the inactivation of bona fide tumor suppressor genes in WM patients. Furthermore, these genes are negative regulators of NF-kB signaling pathway, highlighting its biologic importance in sustaining and limiting growth in WM. Mutational activation of this pathway, which is normally activated by ligand-receptor interactions within the BM microenvironment, suggest a therapeutic role for inhibitors of NF-KB pathway activation in the treatment of WM.

Abstract We have successfully generated the first immuno-competent transgenic mouse model of multiple myeloma (MM), Vk*MYC, in which sporadic MYC activation in the germinal center is AID dependent and mediated by the somatic hypermutation process. All of Vk*MYC mice develop monoclonal PC expansion by 50 weeks of age, that is manifested by high levels of serum IgGs and major M-spikes by SPEP. The PCs in Vk*MYC mice are fully differentiated, somatically mutated, have a very low proliferation index and are found exclusively in the BM. Associated with the long lifespan (medium survival = 661d), Vk*MYC mice develop anemia and diffuse bone disease (osteoporosis with low bone mineral density), with sporadic occurrence of lytic bone lesions and hind limb paralysis. The phenotypic fidelity of the Vk*MYC mice translates in a therapeutic fidelity, as Vk*MYC mice show response to known clinically active anti-myeloma drugs (bortezomib, melphalane, dexamethasone), but not to drugs that have shown no clinical activity as single agents (vincristine, hydroxyurea, fludarabine). Furthermore, acquired bortezomib resistance could efficiently be induced in Vk*MYC mice after only four cycles of suboptimal treatment, and studies are in progress to identify genetic events associated with bortezomib resistance. Preliminary aCGH experiments on a high resolution Agilent 244K array, using genomic DNA from six Vk*MYC CD138+ selected PCs, showed that the genomic complexity of Vk*MYC tumors is quite low, with only few unbalanced rearrangements, such as trisomy of chromosome 15 detected in two samples or trisomy 9 present in another tumor. The most common abnormalities found are small focal monoallelic deletions, six of which are present in more than one tumor. A focal homozygous deletion of several exons of a gene was identified in two tumors samples. Remarkably, a comparative analysis of the syntenic chromosomal region in human MM cell lines and primary patient samples identified similar recurrent homozygous deletion in this putative novel tumor suppressor gene. Overall, these data suggest that the genomic complexity in this inbred strain of transgenic mice is significantly less then in patients with MM, simplifying the search for critical loci involved in MM tumorigenesis and drug resistance and further validate the Vk*MYC mice as a particularly faithful model of human MM.

Abstract Constitutive activation of NF-κB signaling pathways has been identified in a variety of malignancies, however the molecular basis has remained largely unknown. Recently, we have shown mutations in 10 genes encoding positive and negative regulators of NF-κB signaling pathways in 20% of multiple myeloma (MM) patients, with TRAF3 being the most commonly affected gene. TRAF3 is a negative regulator of NF-κB signal and its inactivation results in the constitutive activation of the alternative or noncanonical NF-κB signaling pathway. With the exception of two recent studies in MM, there is no available data regarding the mechanisms that dysregulate NF-κB in other B-cell malignancies. The aim of this study was to analyze the TRAF3 status in a cohort of Waldenström’s macroglobulinemia (WM) patients. Forty-nine WM were analyzed by using a combination of array-based comparative genomic hybridization (aCGH) and cIgM-FISH. Genome sequencing from 28 of these patients was performed on all coding exons of TRAF3 and the adjacent intron-exon junctions. A NF-κB index, as a measure of NF-κB transcription activity, was calculated from 21 patients with gene expression data. The index is the mean expression level of four probe sets differentially expressed between human myeloma cell lines with identified mutations compared to those without identified NF-κB pathway mutations (CD74, IL2RG and 2x TNFAIP3). By performing aCGH we identified a small monoallelic deletion (∼460 Kb), including the TRAF3 locus in one of 22 patients. Genome sequencing identified a missense mutation (D483N) in the remaining allele. This substitution affects a highly conserved amino acid throughout vertebrate evolution, suggesting its functional relevance. We used cIgM-FISH to screen 27 patients with no available aCGH data and we found one biallelic TRAF3 deletion and one chromosome 14 monosomy. In the later case, no DNA was available to sequence the remaining allele. Both samples with biallelic inactivation of TRAF3 possess the two highest NFKB index values of our WM cohort, suggesting NF-κB transcription activity. Moreover, the case with biallelic deletion showed the lowest level of TRAF3 gene expression. Additionally, the entire TRAF3 coding region was sequenced in 28 of these 49 patients and one homozygous nonsense mutation (K286X) was identified. No GEP was available in this sample to calculate the NF-κB index. TRAF3 was only recently implicated in a malignant process, when inactivating mutations and biallelic deletions were reported in approximately 13% of MM patients, making it the most frequently inactivated tumor suppressor in MM. Our study has identified biallelic TRAF3 inactivation in at least 3 of 49 WM patients, suggesting that the constitutive activation of NF-κB signaling pathways may be involved in the pathogenesis of WM. Recent findings showing an association between TRAF3 loss of function and good response to bortezomib in MM patients suggest that proteasome inhibitor-based therapy may be an attractive option in the treatment of the WM patients with constitutive activation of NF-κB signaling pathways. Additional mutational screening of TRAF3 and other regulators of these pathways are ongoing to determine the extent of constitutive activation present in WM and other B-cell malignancies.

Abstract The availability of human myeloma cell lines (HMCLs) has drastically benefited out understanding of multiple myeloma (MM). They represent essential resources to the myeloma research community at large who use them to identify and more importantly validate oncogenic events in MM. Furthermore, they represent the front line tools in efforts to identify novel therapeutic agents for MM patients. However, the myeloma research community is not immune to problems of cell line contamination (mixed populations) or systemic errors where HMCLs are labeled incorrectly. As the community evolves towards the development of tailor medicine it is essential that our basic tools are both pure and correct. Otherwise, a potentially beneficial therapeutic agent designed for a specific patient subset maybe lost simply because of a simple laboratory mix-up at some point over the years. To that end we have developed a simple PCR based method that can be used to validate the purity and authenticity of HMCLs. First two reference collections of over 50 HMCL maintained at the NCI and Mayo Clinic Arizona were compared using microsatellite fingerprinting. This process identified several previously unnoticed discrepancies. First, in both labs several cell lines that reportedly originated from different patients had the same fingerprints. These include the following unexpected pairs; ANBL-6 and DP–6, CAG and ARP–1, KMS–11 and KMS–20, KMS–11 and KMS–21PE/BM, JIM3 and KPMM2. Second, both KP–6 and KAS6/1 showed a level of cross-contamination with ANBL–6. All of the discrepant cell lines were requested from the original contributor and those that were available were retested. After this step we had identified unique fingerprints for ANBL–6, DP–6, KMS–11, and KMS–20. Furthermore, pure populations of KP–6 and KAS6/1 were confirmed. However, CAG and ARP–1 still had an identical fingerprint. Since the use of microsatellite fingerprinting is an expensive task and the identification of cross-contamination is not very robust we set out to develop better method. We decided to use a simple qualitative PCR method, as the binomial presence or absence of a band makes the interpretation of the fingerprint extremely easy. Furthermore, the high sensitivity of a qualitative PCR can pick up a small amount of contamination in the culture that can be lost in microsatellite analysis. This method involves a multiplex PCR screen that simultaneously interrogates 10 regions of the human genome. The screen is divided into two reactions, one with a panel of highly polymorphic copy number variations (CNV) that are often completely absent in some individuals, and a second one based on known biologic events present in the cell lines. The integration of the results from the CNV and Biological Fingerprints will identify each HMCL as a unique entity. We validated the methodology by confirming that this qualitative screen could pick up the cross–contamination identified in the original KP–6 and KAS6/1 samples. As our understanding of myeloma increases and we strive to develop new therapies for specific biological entities it will be essential that groups confirm both the authenticity and purity of their test cultures. We hope that widespread use of this simple screening methodology will significantly improve the quality and reproducibility of assays carried out in labs around the world.

Abstract Events mediating transformation from the pre-malignant MGUS to MM is largely unknown. Previous studies have highlighted RAS mutation as a progression event. However, this abnormality is present in only a small subset of patients. Global gene expression differences between MGUS and MM have been reported but provided little coherent biological insights. To define progression events, we examined the differences in gene expression profiles (GEP, Affymetrix U133A genechip) between 101 MM and 22 MGUS using a powerful functional genomic tool, Gene Set Enrichment Analysis (Subramanium et al PNAS 2005), that allow extraction of biological information based on enrichment of genes associated with a phenotype within gene-sets derived experimentally or extracted from database. We found that a predominant theme relating to MM progression is the over-expression of proliferation related genes in MM highlighted by the enrichment of cell cycle related gene-sets. Another common theme is the enrichment of gene-sets relating to MYC activation. Seven of eight gene-sets related to MYC activation in the database are significantly enriched (p-value <0.05 and false detection rate <0.05). Leading edge (enriched genes) analysis of enriched MYC-related and proliferation-related gene-sets revealed a core, overlapping signature of 8 genes specific to MYC and not proliferation. These genes are further verified as MYC transcriptional targets using Metacore, a network analysis tool. In addition, 7 of the 8 genes had been validated as MYC target genes by chromatin immunoprecipitation previously (www.myccancergene.org). We validated this signature in independent publicly available cell lines and mouse gene expression experiments where this signature is present in mammary epithelial cell lines engineered to over-express MYC but not E2F3, SRC and BCAT, and mouse embryonic fibroblast over-expressing MYC but not HER2 or with RB knocked-out. In addition, analysis of proliferation index in these datasets showed that there is dissociation of proliferation from this signature. This showed conclusively that our signature is MYC specific and not related to proliferation. Of interest, the signature is also present in cell lines with activating RAS mutation, which also have high expression of MYC, suggesting that MYC activation is a downstream effect of activated RAS signaling. Indeed, in our patient cohort, most patients with RAS mutation express this signature. When applied to 2 large GEP MM patient cohorts, this signature was not present in normal plasma cells, rarely present in MGUS, and present in a significant proportion of MM and almost all the human myeloma cell lines. When analyzed in the context of other B-cell malignancies, the signature is very strong in Burkitt’s lymphoma, present in most myeloma but absent in chronic lymphocytic leukemia, Waldenstrom Macroglobulinemia, and normal B-cells. Our analysis therefore suggests that MYC activation may be a unifying pathological event in MGUS to MM progression. As MYC translocations occur in only about 13% of newly diagnosed MM patients, activating RAS mutations and other trans-mechanisms may mediate MYC activation. Our recent demonstration that sporadic MYC activation in a mouse strain prone to monoclonal gammopathy led to MM development provides important functional proof.

8007 Background: We recently described mutations in multiple genes in MM patients that lead to the activation of the non- canonical NF-KB pathway. Among the mutated genes, TRAF3 is the most commonly affected, and is inactivated in about 20% of patients. In this study, we further defined the clinical significance of TRAF3 inactivation in MM. Methods: Deletion of TRAF3 was identified by array- based comparative genomic hybridization using Agilent Human Genome CGH 44k oligonucleotide array and confirmed by fluorescent in situ hybridization. Mutations were detected by sequencing. We use a TRAF3 gene expression cut-off (derived from 127 Mayo clinic patients) as a surrogate for TRAF3 inactivation to investigate its association with Translocations and Cyclin D genetic subtypes, response rates, and survival using 2 large gene expression dataset: UAMS dataset (559 newly diagnosed patients treated with Total therapy II and III) and Millenium datasets (213 relapsed patients entered into randomized study of bortezomib versus dexamethasone). Results: A normalized expression level of less than 0.6 identified all but three of the known TRAF3 inactivated cases in the Mayo dataset (p<0.0001). Using this cut-off, TRAF3 inactivated cases were enriched for t(4;14) (26%, p=0.00002) and rarely D1 (mostly hyperdiploid) (9%, p=9.6E-10) in all 3 datasets. TRAF3 inactivation did not affect survival of patients treated with transplant-based therapy (UAMS and Mayo dataset) but was strongly associated with bortezomib response (89% TRAF3 low versus 40% TRAF3 normal, p0.6 had similar PFS with the two treatment. Conclusion: Mutations resulting in TRAF3 inactivation are common, particularly in non-hyperdiploid MM, and predict for response to bortezomib and resistance to dexamethasone suggesting that inhibition of the non-canonical NF-KB pathway may be one of the principle mechanisms of action of proteasome inhibitors in MM. No significant financial relationships to disclose.

Abstract The most common genetic abnormality in multiple myeloma is the deletion of chromosome 13, seen in almost half of newly diagnosed patients. Unlike chronic lymphocytic leukemia, where a recurrent minimally deleted region including MIR15A/MIR16-1 has been mapped, the deletions in multiple myeloma predominantly involve the entire chromosome and no specific driver gene has been identified. Additional candidate loci include RB1 and DIS3, but while biallelic deletion of RB1 is associated with disease progression, DIS3 is a common essential gene and complete inactivation is not observed. The Vk*MYC transgenic mouse model of multiple myeloma spontaneously acquires del(14), syntenic to human chromosome 13, and Rb1 complete inactivation, but not Dis3 mutations. Taking advantage of this model, we explored the role in multiple myeloma initiation and progression of two candidate loci on chromosome 13: RB1 and MIR15A/MIR16-1. Monoallelic deletion of Mir15a/Mir16-1, but not Rb1, was sufficient to accelerate the development of monoclonal gammopathy in wild-type mice and the progression of multiple myeloma in Vk*MYC mice, resulting in increased expression of Mir15a/Mir16-1 target genes and plasma cell proliferation, which was similarly observed in patients with multiple myeloma. Significance: In the absence of a defined, minimally deleted region the significance of del(13) in multiple myeloma has remained controversial. Here we show that haploinsufficiency of Mir15a/Mir16-1, but not Rb1, upregulates the cell cycle–regulatory network, inducing monoclonal gammopathy in mice and promoting multiple myeloma progression in both mice and men. See related commentary by Walker, p. 16. This article is highlighted in the In This Issue feature, p. 5