Purpose: To determine the recurring DNA copy number alterations (CNA) in classical Hodgkin lymphoma (HL) by microarray-based comparative genomic hybridization (aCGH) using laser capture microdissected CD30+ Hodgkin and Reed–Sternberg (HRS) cells.

Experimental Design: Archived tissues from 27 CD30+ HL plus control samples were analyzed by DNA microarrays. The HL molecular karyotypes were compared with the genomic profiles of germinal center B cells and treatment outcome (chemotherapy responsive vs. primary refractory disease).

Results: Gains and losses observed in more than 35% of HL samples were localized to 22 and 12 chromosomal regions, respectively. Frequent gains (>65%) were associated with growth and proliferation, NF-κB activation, cell-cycle control, apoptosis, and immune and lymphoid development. Frequent losses (>40%) observed encompassed tumor suppressor genes (SPRY1, NELL1, and ID4, inhibitor of DNA binding 4), transcriptional repressors (TXNIP, thioredoxin interacting protein), SKP2 (S-phase kinase-associated protein 2; ubiquitin ligase component), and an antagonist of NF-κB activation (PPARGC1A). In comparison to the germinal center profiles, the most frequent imbalances in HL were losses in 5p13 (AMACR, GDNF, and SKP2), and gains in 7q36 (SHH, sonic hedgehog homolog) and 9q34 (ABL1, CDK9, LCN2, and PTGES). Gains (>35%) in the HL chemoresponsive patients housed genes known to regulate T-cell trafficking or NF-κB activation (CCL22, CX3CL1, CCL17, DOK4, and IL10), whereas the refractory samples showed frequent loss of 4q27 (interleukin; IL21/IL2) and 17p12, and gain of 19q13.3 (BCL3/RELB).

Conclusion: We identified nonrandom CNAs in the molecular karyotypes of classical HL. Several recurring genetic lesions correlated with disease outcome. These findings may be useful prognostic markers in the counseling and management of patients and for the development of novel therapeutic approaches in primary refractory HL. Clin Cancer Res; 17(10); 3443–54. ©2011 AACR.

Translational Relevance

We show the potential of identifying copy number alterations (CNA) in HL by aCGH using DNA extracted from microdissected CD30+ Hodgkin and Reed–Sternberg (HRS) cells. The molecular data show a clear relationship between CNAs of genes associated with the NF-κB signaling pathway, the tumor microenvironment, cell-cycle regulation, and apoptosis. We provide a comparison of chemosensitive and primary refractory Hodgkin lymphoma (HL) samples and describe potentially pathogenetic and prognostic CNA differences at disease onset. In particular, HL samples that showed losses of IL21/IL2 at 4q27, 11p14.3/SLC17A6, and 17p12, with gains of BCL3/19q13.31 were associated with primary refractory HL, whereas HL samples showing 16q13 gains or few genetic aberrations were most frequently associated with favorable IPS (international prognostic scores) and chemosensitive phenotypes. Due to the limited number of well-characterized HL samples available for ongoing studies, especially for primary refractory HL, collaboration among lymphoma clinical investigators is imperative to refine, and build on these promising clinical predictors.

The annual incidence of Hodgkin lymphoma (HL) is estimated at 3 cases per 100,000 persons, making this malignancy one of the most common lymphomas in the Western world (1). The characteristic pathologic feature of classical HL is the presence of Hodgkin and Reed-Sternberg (HRS) cells, which usually comprise less than 3% of the affected mixed cellular lesion. There is compelling evidence suggesting that the pathognomonic HRS cells are an outgrowth of a malignant clone derived from a “reprogrammed” germinal center (GC) B cell that no longer expresses B-lineage–specific genes, such as POU2F2, POU2AF1, and PU.1, and may express genetic markers characteristic of other hematopoietic lineages like IL21 (interleukin), CCL17, CSF1, ID2, CD3, and CD4 (2–4). Recurrent genetic lesions in critical hematopoietic transcription factors have led to the discovery that constitutive activation of the NF-κB signaling pathway is essential for HRS cell survival and proliferation (4). In particular, gains of REL and deletions or inactivating mutations of TNFAIP3 (TNF-α–induced protein 3), gene involved in the NF-κB signaling pathway, have been detected in about 40% of classical HL (5–7). To gain further insight into the pathogenesis of HL and its lineage infidelity, comprehensive genomic characterization of HRS cells from primary HL tumor samples is necessary; however, the scarcity of HRS cells in the HL lesions remains a key limiting factor in unraveling the molecular consequences of this malignancy.

Treatment with doxorubicin-based combination chemotherapy regimens leads to complete remission in the majority (70%–90%) of HL patients (1). Unfortunately, HL patients who fail front-line or second-line therapies, including stem cell transplantation, are quite often young and have a poor median survival (<3 years; ref. 8). In primary refractory HL (PR HL), the poor response to initial therapy may be attributed to nonrandom genetic alterations that underlie disease aggressiveness at onset. Comparison of the genetic aberrations of HRS cells from patients who respond favorably to standard HL treatment protocols with those who fail treatment has the potential to provide critical insight into the clinical behavior of HL and guide future treatment (9). Moreover, to develop novel HL therapies directed against rational targets, we need focused investigations that lead to more thorough understanding of the genetic alterations of HRS, its lineage infidelity consequences, and drug resistance.

Until recently, the majority of genetic and molecular biology HL data was gathered from studies using HL continuous cell lines established from patients with relapsed HL (10–13). Distinguishing primary genetic alterations from secondary events in HL pathogenesis remains a challenge. Recently, high-resolution cytogenetic techniques that combine laser capture microdissection with microarray-based comparative genomic hybridization (aCGH) technology have provided new opportunities to investigate genome-wide DNA alterations in lesions of limited size (13–15). Genomic DNA–based aCGH experiments show less variability than mRNA expression arrays, and suspected “cryptic” DNA aberrations can be confirmed easily by FISH. Importantly, robust prognostic features for many hematologic malignancies have been revealed using a combination of clinical parameters and genetic data gathered from conventional cytogenetics and molecular techniques. Given the paucity of genomic analyses of HRS cells, we sought to characterize their DNA copy number alterations (CNA) with bacterial artificial chromosome (BAC)-based aCGH using DNA extracted from laser capture microdissected (LCM) CD30+ HRS cells. The objectives of this study were to determine the recurring CNAs in HL, compare the findings of the malignant HL lesions to those detected in “benign” GC B cells–the normal counterpart of HRS cells, and define the most common CNA differences between the chemotherapy responsive (CR) and PR HL samples.

Patient and control samples

On approval from the Institutional Review Board of the City of Hope, primary formalin-fixed paraffin-embedded (FFPE) diagnostic samples were obtained from 27 patients including 15 patients with CR HL and 12 patients with PR HL. Clinicopathologic characteristics of the patients are summarized in Table 1. The international prognostic scores (IPS) of our study population were calculated as described (16). In this study, patients with an IPS of 2 or less were given a favorable designation and patients with an IPS of 3 or above were assigned to the unfavorable group. Control samples included 9 FFPE benign lymph node samples (4 males and 5 females) and GC cells from 10 FFPE reactive follicular hyperplasia (RFH) samples (2 males and 8 females). Each sample was submitted for conventional histopathologic processing to confirm HL involvement in the test samples or no evidence of malignancy in the control samples.

Table 1.

Clinical and pathologic characteristics of the HL patients

No. of patientsPercentage
Agea (range), y 32 (5–64) 27 100 
Gender (M/F) 13/14 27 48/52 
REAL/WHO classification 
 Nodular sclerosis  23 85 
 Mixed cellularity  11 
 Classical Hodgkin's disease, NOS  
Stage 
 Early (IA–IIB)  16 59 
 Advanced (IIIA–IVB)  11 41 
A versus B symptoms 
 A  13 48 
 B  14 52 
Bulky disease  26 
Bone marrow involvement 
IPS scoreb 
 Favorable  30 
 Unfavorable  17 63 
 Indeterminate  
Initial treatment 
 ABVD  20 74 
 Stanford V  11 
 BEACOPPc  11 
 COPP/ABV  
Treatment response 
 Sensitive  15 56 
 Refractory  12 44 
Transplant 
 Autologous  11 41 
 Allogeneic  
No. of patientsPercentage
Agea (range), y 32 (5–64) 27 100 
Gender (M/F) 13/14 27 48/52 
REAL/WHO classification 
 Nodular sclerosis  23 85 
 Mixed cellularity  11 
 Classical Hodgkin's disease, NOS  
Stage 
 Early (IA–IIB)  16 59 
 Advanced (IIIA–IVB)  11 41 
A versus B symptoms 
 A  13 48 
 B  14 52 
Bulky disease  26 
Bone marrow involvement 
IPS scoreb 
 Favorable  30 
 Unfavorable  17 63 
 Indeterminate  
Initial treatment 
 ABVD  20 74 
 Stanford V  11 
 BEACOPPc  11 
 COPP/ABV  
Treatment response 
 Sensitive  15 56 
 Refractory  12 44 
Transplant 
 Autologous  11 41 
 Allogeneic  

Abbreviation: NOS, not otherwise specified.

aAge, median (y).

bIPS is based on 7 potential unfavorable features at diagnosis: serum albumin ≤4 g/dL (40 g/L), hemoglobin ≤10.5 g/dL (105 g/L), male gender ≥45 years, stage IV, WBC ≥15,000/μL, and lymphocyte count ≤600/μL, and/or ≤8% of the WBC count (16).

cPediatric protocol CCG-59704.

FFPE tissue processing and laser capture microdissection

Five-micrometer serial sections from the FFPE tissue blocks were fixed onto PALM membrane slides (PEN-membrane; Zeiss) and processed as previously described (13, 14). A series of experiments designed to assess the impact of DNA source (e.g., archival material including frozen and FFPE), quantity, and amplification on array CGH was carried out to establish the FFPE aCGH protocol (13). Briefly, the slides were pretreated as follows for immunostaining: 1 hour at 65°C in a dry oven, 1 minute in xylene at room temperature, 5 minutes in 100% ethanol (2×), 5 minutes in 3% H2O2, and rinsed in dH2O. Antigen retrieval was carried out at 98°C for 30 minutes using the decloaking chamber and decloaker universal heat retrieval buffer (both from Biocare Medical). Test sample slides were stained with monoclonal mouse anti-human CD30 (Dako Inc.) using the DakoCytomation Autostainer (Dako Inc.) per manufacturer's instructions. Control samples were stained with CD20. A total of 150 CD30+ HRS cells per test sample and 150 CD20+ GC cells per control sample were LCM into microfuge caps containing water.

Whole-genome amplification

Isolated LCM cells were randomly fragmented and amplified in the LCM collection tubes using the GenomePlex Single-Cell Whole-Genome Amplification Kit (WGA4; Sigma-Aldrich) according to product instructions. The quality and quantity of the isolated amplified DNA were assessed electrophoretically and spectrophotometrically. Samples were processed prior to aCGH using the GenElute PCR Clean-Up Kit (Sigma-Aldrich). A single sample was run in triplicate to evaluate assay reproducibility.

BAC-based aCGH labeling and hybridization

One microgram of genomic DNA was labeled using Bioarray Kit (Enzo Life Sciences) prior to subsequent hybridization to a RPCI 19K BAC array (Roswell Park) using a Gentac hybridization station (Genomic Solutions). Sex-mismatched pooled DNA from 20 healthy donors was used as reference DNA. The minimal tiling RPCI BAC array containing about 19,000 BAC clones was used (backbone consisted of ∼4,600 BAC clones with a BAC center-to-center distance of 165 kb; human genome build 36.1; ref. 13). Slides were scanned on the Genepix 4200A Scanner (Molecular Devices).

Control assays, selection criteria, and outcome predictor

Two control studies were conducted to establish experimental reproducibility and quality control associated with multiple features of our study. These control studies had the objective of establishing standards for processing FFPE tissue, WGA, or eliminate problematic microarray calls associated with cross-hybridizing BACs or poor labeling/hybridization efficiency. First, 9 FFPE benign lymph node control samples were processed in an identical fashion to the test samples. Second, 1.05 ng of peripheral blood DNA (DNA equivalent of 150 cells) obtained from 6 normal, healthy individuals (3 males and 3 females) was amplified and analyzed by aCGH. Deviant calls from either of these control studies found in 2 or more microarrays were masked before the HL test samples were processed. Selection criteria for the reportable CNA calls in the HL samples were more than 5 BACs (∼700 kb in size) observed in 35% or more in either the sensitive or resistant samples with a value of P < 0.15 and not present in less than 20% (no more than 1 patient) of the benign control samples.

To determine which genes are presumably critical in the HL malignant process, GC B cells from 10 lymph node samples showing RFH were processed identically to the test samples and compared with the HL CNAs. Genes associated with RFH had to meet the following selection criteria: large (>700 kb) CNAs observed 15% or more commonly in the RFH group (2 or more samples) with a value of P < 0.15, and observed in less than 20% of the benign control group. The CNAs associated with RFH were filtered prior to the comparative analysis of the chemosensitive and primary refractory disease CNAs. This final analysis generated a list of genes that could be used as a potential outcome predictor to discriminate between CR HL and PR HL at presentation. To test the outcome predictor in a preliminary manner, 4 additional primary HL lesions were analyzed in a blinded yet similar manner to the original 27 HL test samples and categorized using the proposed predictor model.

Data analysis

Data images were quantified by ImaGene (BioDiscovery, Inc.), and the chromosome regions of abnormal copy number across the genome were detected by Nexus (BioDiscovery, Inc). To explore gene function and interactions, the CNAs were uploaded into Ingenuity Pathway Analysis tools 3.0 (IPA; http://www.ingenuity.com) for alternative and canonical pathway mapping. IPA calculates a per network significance score by applying the hypergeometric distribution calculated via the Fisher's exact test for 2 × 2 contingency tables. The P value is calculated by comparing the number of user-specified genes of interest that participate in a given pathway, relative to the total number of occurrences of these genes in all pathway annotations stored in the Ingenuity Pathways Knowledge Base (IPKB).

FISH confirmation studies

Locus-specific FISH analyses were conducted to confirm the recurring CNAs detected by aCGH. The probes used for FISH analysis were carefully chosen and mapped within the genomic coordinates of the specific chromosome region showing gain or loss. In each case, a control DNA FISH probe from the opposite chromosome arm was included to confirm the CNAs in relation to the ploidy level of the HRS cell. To validate CNAs by FISH, CD30-stained IHC slides with AEC detection were scanned, localized, and recorded on the slide using the Bioview Duet Image Analyzer (BioView, Ltd.) prior to FISH analysis as described previously (17).

Statistical evaluation

We background corrected and normalized the quantified images with print-tip loess after removing bad spots. We evaluated all pairwise correlations between normalized replicates by scatter plots and Pearson's correlations on the log2 ratios.

Defining the threshold and reproducibility of CNAs

From each patient and control FFPE sample, 150 cells (CD30+ HRS cells from test samples and CD20+ cells from control samples) were isolated from 5-μm thick tissue sections by LCM, amplified using WGA, and hybridized to a 19K whole-genome BAC array. Two control studies were conducted to reduce deviant calls associated with poor hybridization efficiency or amplification bias associated with WGA or FFPE tissues. In the first experiment, 1.05 ng of peripheral blood DNA (DNA equivalent of 150 cells) obtained from 6 normal individuals (3 males and 3 females) were WGA amplified and analyzed by aCGH. In the second study, 9 formalin-fixed benign lymph node control samples (4 males and 5 females) were processed similar to the test samples. A total of 235 of the 16,952 (1.4%) autosomal BACs showed a false-positive (gain or loss) result in the normal control samples; these BACs were masked from the test data analyses. In addition, no chromosomal site was called positive if its respective BACs showed gain or loss in more than 2 arrays from normal control samples. On the basis of these control data, and to minimize the risk of reporting false-positive results, thresholds for calling gains and losses were set at more than 35% of test samples. Reproducibility of the FFPE was carried out by analyzing 150 HRS cells from a single HL patient that were isolated independently from the block and processed by the LCM-WGA-aCGH method in triplicate. Supplementary Figure S1 shows the high correlation between the normalized log2 ratio values of each possible pair of replicates (Pearson r = 0.90, P values uniformly <0.0001), confirming reproducibility of the assay and the array.

CNAs detected in 27 primary FFPE diagnostic HL samples

Primary tissue blocks were obtained from 15 CR and 12 PR patients. The clinical and pathologic characteristics of the patients are listed in Table 1. Recurring chromosome imbalances were observed among all 27 HL samples. The most frequent gains (>65%) mapped to 2q37.3, 7p21.1, 8q24.3, 16q23.3-q24.3, 9q34.13-q34.3,14q32.33, 19q13.33, and 20q13.33 (Supplementary Data; Table 1) and included genes associated with growth and proliferation (AHR, BAI1, BOP1, COMMD5, LY6E, PTP4A3, SLURP1, CBFA2T3, SLC7A5, NOTCH1, FOXF1, TRAF2, IRF8, S100B, and MYH14), cell cycle (AKT1, CDK10, and SUMO3), drug metabolism (CYP11B1, CYP11B2, and SLC19A1), angiogenesis and cell adhesion (COL18A1, CDH4, ITGB2, and FOXC2), apoptosis regulation (FOXF1 and GPR132), immune and lymphatic development (CBFA2T3, IL17C, IRF8, CLEC11A, RXRA, SPIB, and ICOSLG) and invasion, metastasis or cancer relatedness (VAV2, PSCA, PTP4A3, GINS2, FUT7, TUBB2C, KLK, POLD1, and TFF2). Losses observed in more than 40% of HL samples encompassed SPRY1, NELL1, SLC1A3, GDNF, IL7R, SKP2 (S-phase kinase-associated protein 2), ID4 (inhibitor of DNA binding 4), PPARGC1A, and TXNIP (thioredoxin interacting protein). Because ploidy differences are inherent to HRS cells, to accurately validate the CNAs observed by aCGH, HRS-specific FISH confirmation studies using locus-specific probes that mapped within the genetic coordinates of the recurring CNAs were paired with a DNA control probe that localized to the opposite chromosome arm. These chromosome region–specific FISH analyses confirmed the presence or absence of aCGH-defined CNAs in CD30+ HRS for 15 regions in 3 to 5 patient samples for each paired FISH probe set (Fig. 1).

Figure 1.

FISH confirmation study to validate CNAs detected by aCGH. Top, genome aggregate for 12 PR HL samples. Arrow, gain of 19q13.31 in 83% of refractory samples. Middle, chromosome 19 drill down for an individual sample. Each circle represents an individual BAC on the array. Green line, threshold for gain; red line, threshold for CNA loss; yellow line, moving average; black lines, cluster values of gain or loss defined by a circular binary segmentation algorithm. The highlighted BAC (pink) was selected for FISH validation studies. RP11-876a24 BAC contains the genes for PVR (poliovirus receptor) and CEACAM19. Bottom, representative FISH images for sample BK01. A, bright-field image of CD30+ cell. B, the corresponding rhodamine filter image. RP11-876a24 was labeled with digoxigenin and detected with rhodamine anti-dig. Nine copies of this clone were detected. C, FITC (fluorescein isothiocyanate) image of internal control probe, only 6 copies were detected (net gain of 3 additional copies of PVR in cell).

Figure 1.

FISH confirmation study to validate CNAs detected by aCGH. Top, genome aggregate for 12 PR HL samples. Arrow, gain of 19q13.31 in 83% of refractory samples. Middle, chromosome 19 drill down for an individual sample. Each circle represents an individual BAC on the array. Green line, threshold for gain; red line, threshold for CNA loss; yellow line, moving average; black lines, cluster values of gain or loss defined by a circular binary segmentation algorithm. The highlighted BAC (pink) was selected for FISH validation studies. RP11-876a24 BAC contains the genes for PVR (poliovirus receptor) and CEACAM19. Bottom, representative FISH images for sample BK01. A, bright-field image of CD30+ cell. B, the corresponding rhodamine filter image. RP11-876a24 was labeled with digoxigenin and detected with rhodamine anti-dig. Nine copies of this clone were detected. C, FITC (fluorescein isothiocyanate) image of internal control probe, only 6 copies were detected (net gain of 3 additional copies of PVR in cell).

Close modal

Comparison of aCGH CNAs in 15 chemosensitive and 12 PR HL samples

The genomic differences observed at more than 35% in the CR HL patients clustered within 2 chromosomal bands on chromosome 16 (16q13 and 16p11.2) and showed gains of genes known to regulate T-cell trafficking or NF-κB activation (CCL22, CX3CL1, CCL17, DOK4, and IL10; Table 2). In comparison to the CR samples, the 12 PR patients showed frequent loss of genes localized to 1q21.2, 4q27, 5p15, 11p14.3, and 17p12 with recurrent gains of chromosome bands 2q37.2, 13q31.2, 19q13.3, and 22q13.33. These data were also compared with the IPS of our study population (Table 1). Of interest, the 39 unique genes we identified in regions with CNAs comparing patients with favorable versus unfavorable IPS scores showed remarkable overlap with the 63 unique genes identified in chemosensitive and primary resistant genomic profiles (Table 2), with 36 unique genes appearing in both lists. The probability of this occurring by chance follows a hypergeometric distribution (18, 19), and the associated P value was less than 0.0001.

Table 2.

Consistent CNAs observed in CR HL versus PR HL, not observed in GC cells

BandRegionEventLength, bpFrequency in refractoryFrequency in sensitiveGenes in sensitive points
16q13 chr16:55,456,619-56,127,978 Gain 671,359 25% 60% ARL2BP, CCDC102A, CCL17, CCL22, CETP, CIAPIN1, COQ9, CPNE2, CX3CL1, DOK4, HERPUD1, NIP30, NLRC5, PLLP, POLR2C, RSPRY1, SLC12A3 
16p11.2 chr16:32,592,341-33,172,220 Loss 579,879 17% 53% LOC729355, TP53TG3 
1q21.2 chr1:148,603,613-148,818,760 Loss 215,147 42% 7% ADAMTSL4, ECM1, KIAA0460, MCL1, TARS2 
2q37.2 chr2:236,067,474-236,698,859 Gain 631,385 50% 27% CENTG2 
4q27 chr4:120,634,997-123,761,662 Loss 3,126,665 50% 20% ADAD1, ANXA5, BBS7, C4orf31, CCNA2, EXOSC9, GPR103, IL2, IL21, KIAA1109, MAD2L1, PDE5A, TMEM155, TNIP3, TRPC3 
5p15.31 chr5:7,449,342-7,883,194 Loss 433,852 42% 0% ADCY2 
5p15.2 chr5:9,088,137-9,683,463 Loss 595,326 42% 7% SEMA5A, TAS2R1 
11p14.3 chr11:22,171,297-22,257,975 Loss 86,678 58% 7% TMEM16E 
13q31.2 chr13:87,122,870-87,129,871 Gain 7,001 58% 27% SLITRK5 
17p12 chr17:11,864,859-15,526,918 Loss 3,662,059 58% 13% CDRT1, CDRT15, CDRT4, COX10, ELAC2, FAM18B2, FLJ45831, HS3ST3A1, HS3ST3B1, MAP2K4, MYOCD, PMP22, RICH2, TEKT3, TRIM16 
19q13.31 chr19:49,894,197-49,955,141 Gain 60,944 50% 20% BCL3, RELB, CEACAM19 
22q13.33 chr22:49,092,963-49,100,010 Gain 7,047 50% 20% FAM116B 
BandRegionEventLength, bpFrequency in refractoryFrequency in sensitiveGenes in sensitive points
16q13 chr16:55,456,619-56,127,978 Gain 671,359 25% 60% ARL2BP, CCDC102A, CCL17, CCL22, CETP, CIAPIN1, COQ9, CPNE2, CX3CL1, DOK4, HERPUD1, NIP30, NLRC5, PLLP, POLR2C, RSPRY1, SLC12A3 
16p11.2 chr16:32,592,341-33,172,220 Loss 579,879 17% 53% LOC729355, TP53TG3 
1q21.2 chr1:148,603,613-148,818,760 Loss 215,147 42% 7% ADAMTSL4, ECM1, KIAA0460, MCL1, TARS2 
2q37.2 chr2:236,067,474-236,698,859 Gain 631,385 50% 27% CENTG2 
4q27 chr4:120,634,997-123,761,662 Loss 3,126,665 50% 20% ADAD1, ANXA5, BBS7, C4orf31, CCNA2, EXOSC9, GPR103, IL2, IL21, KIAA1109, MAD2L1, PDE5A, TMEM155, TNIP3, TRPC3 
5p15.31 chr5:7,449,342-7,883,194 Loss 433,852 42% 0% ADCY2 
5p15.2 chr5:9,088,137-9,683,463 Loss 595,326 42% 7% SEMA5A, TAS2R1 
11p14.3 chr11:22,171,297-22,257,975 Loss 86,678 58% 7% TMEM16E 
13q31.2 chr13:87,122,870-87,129,871 Gain 7,001 58% 27% SLITRK5 
17p12 chr17:11,864,859-15,526,918 Loss 3,662,059 58% 13% CDRT1, CDRT15, CDRT4, COX10, ELAC2, FAM18B2, FLJ45831, HS3ST3A1, HS3ST3B1, MAP2K4, MYOCD, PMP22, RICH2, TEKT3, TRIM16 
19q13.31 chr19:49,894,197-49,955,141 Gain 60,944 50% 20% BCL3, RELB, CEACAM19 
22q13.33 chr22:49,092,963-49,100,010 Gain 7,047 50% 20% FAM116B 

Comparison of HL samples with benign RFH

To delineate which genes are more likely to play a critical role in HL oncogenesis, normal GC cells from 10 RFH samples were analyzed and compared with the HL results. As described in Materials and Methods, strict criteria were established to call recurring gains and losses associated with CR HL, PR HL, and benign RFH. Figure 2A shows the number of genes concordant and discordant among these 3 groups observed in more than 35% of the samples tested in each group. The first assessment of these data identified 151 genes frequently altered in both CR HL and PR HL samples but not detected in the benign RFH samples, 81 genes commonly detected in the CR HL and 312 genes more frequently observed in the PR HL tumors and not in the other 2 groups. The 10 most frequent genetic alterations detected among the 151 HL CNAs included losses of genes in 5p13 (AMACR, GDNF, and SKP2), gains of 7q36 (SHH, sonic hedgehog homolog), and gains of 9q34 (ABL1, CDK9, ENG, LCN2, PTGES, and TSC1). A second assessment of the CR HL and PR HL gene subsets was conducted requiring a frequency difference of more than 20% and a value of P = 0.15 between the sensitive and resistant samples (Fig. 2B). In this analysis, 19 genes were clearly more common in chemosensitive samples, whereas 44 genes were more prevalent in the primary refractory tumors (Table 2).

Figure 2.

Comparison of HL samples with benign follicular hyperplasia (BFH). A, concordant and discordant genes among RFH follicular and the CR HL–sensitive and the PR HL–resistant samples. A total of 151 genes were frequently detected in both CR HL and PR HL samples but not in the RFH samples, 81 genes were commonly detected in the CR HL, and 312 genes more frequently observed in the PR HL tumors and not in the other 2 groups. B, a second assessment of the CR and PR genes showed the differences and similarity between the sensitive (sen) and resistant (res) samples.

Figure 2.

Comparison of HL samples with benign follicular hyperplasia (BFH). A, concordant and discordant genes among RFH follicular and the CR HL–sensitive and the PR HL–resistant samples. A total of 151 genes were frequently detected in both CR HL and PR HL samples but not in the RFH samples, 81 genes were commonly detected in the CR HL, and 312 genes more frequently observed in the PR HL tumors and not in the other 2 groups. B, a second assessment of the CR and PR genes showed the differences and similarity between the sensitive (sen) and resistant (res) samples.

Close modal

CNAs interaction networks and predictive modeling

Possible genetic interactions for those genes listed in Table 2 were examined in the context of the curated list of published molecular interactions in IPKB. In addition to an association with the NF-κB signaling pathway, IPA revealed connections between genes in the calcium signaling and IL2/IL21 signaling pathway (Fig. 3). Key genes from these studies were used to construct a putative outcome model for predicting chemosensitive or primary refractory disease (Fig. 4). Of interest, all 7 patients showing loss of the SLC17A6 (solute carrier family 17, sodium-dependent inorganic phosphate cotransporter member 6) gene had PR disease; however, the 4 patients showing gains of the SLC12A3 (solute carrier family 12 sodium/chloride transporters member 3) gene were CR. Second, patients showing no gain of 16q13 and loss of 4q27 or 17p12 were commonly associated with a refractory phenotype (9 of 11 patients). Third, patients showing 16q13 gain with either few CNAs (10 of 12), or no loss of 4q27 or 17p12 (3 of 3) were usually CR.

Figure 3.

Calcium signaling and IL21/IL2 pathway interactions in HL. Top scoring network of interactions among the differential CNAs in CR HL versus PR HL. The interactions indicate differences in the Ca2+ and IL21/IL2 signaling pathways. Inset, legend functional classes of the genes noted. Solid line, direct interaction. Dotted line, indirect interaction.

Figure 3.

Calcium signaling and IL21/IL2 pathway interactions in HL. Top scoring network of interactions among the differential CNAs in CR HL versus PR HL. The interactions indicate differences in the Ca2+ and IL21/IL2 signaling pathways. Inset, legend functional classes of the genes noted. Solid line, direct interaction. Dotted line, indirect interaction.

Close modal
Figure 4.

Gene profile and prediction outcome model. R, primary refractory patient; S, chemoresponsive patient; s, initially chemoresponsive but relapsed 10 years after disease onset; X, the prediction did not match with the medical record.

Figure 4.

Gene profile and prediction outcome model. R, primary refractory patient; S, chemoresponsive patient; s, initially chemoresponsive but relapsed 10 years after disease onset; X, the prediction did not match with the medical record.

Close modal

On the basis of the relatively high frequency of IL21/IL2 loss at 4q27 and BCL3 overrepresentation in our PR HL cohort, we hypothesized samples without gains in 16q13 and losses of IL21/IL2 at 4q27,11p14.3/SLC17A6, or 17p12 with or without concurrent gains of BCL3 portend a highly likely drug-resistant HL phenotype. To test this CR versus PR outcome model in a preliminary manner, we first recategorized all 27 patients according to the proposed model and predicted the outcome correctly in 23 (85%) of 27 patients (Fig. 4). Four newly acquired HL samples were then tested in a blind fashion. One sample showed loss of 4q27 without BCL3 gains and was correctly assigned primary refractory; the remaining 3 cases were correctly designated “chemosensitive;” however, 1 “chemosensitive” patient (BK21) who initially responded to treatment suffered relapse of disease 10 years later.

The identification of recurring cytogenetic and molecular alterations in HL has been thwarted by the infrequent occurrence of HRS cells in affected lesions (4). To gain further insight into the molecular genetic alterations in HL, we used a combined LCM-WGA-aCGH approach to analyze the genomic composition of HRS cells from 27 diagnostic samples from patients with classical HL. As expected, the molecular karyotypes of the HRS cells revealed numerous CNAs with recurring gains more common than losses due to the ploidy variation associated with mononucleated Hodgkin and multinucleated Reed–Sternberg cells. To refine our analysis, we limited this study to include the most common CNAs, defined as gains and losses in more than 35% of the HL samples, and confirmed by chromosome region–specific FISH analyses in CD30+ HRS.

Previous studies have shown that multiple signaling pathways and transcription factors are deregulated in HRS cells including NF-κB, JAK (Janus-activated kinase)-STAT, and PI3K (phosphoinositide-3-kinase)-AKT (4, 7, 20). Interestingly, regions with the most frequent gains in this study housed NOTCH1 and TRAF2, transcription factors, which mediate NF-κB activation (21, 22). Gains at the AKT1 locus support a role for constitutive activation of PI3K in HL (20, 23) and further support data linking activated AKT1 to NF-κB signaling, p53 inhibition, and enhanced resistance to apoptosis (24). Overrepresentation of the FOXC2 gene, a known transcriptional regulator of CXCR4 and delta-like 4 (DLL4; ref. 25) implies the downregulation of the Notch inhibitor Deltex1 to promote constitutive activation of NOTCH1 activity in HRS cells (21). In this study, gains of the forkhead (FOX) transcription factor FOXF1 gene were consistently observed with gains of SHH, supporting the supposition that NF-κB regulates SHH expression to further advance NF-κB–mediated proliferation and apoptosis resistance (26, 27). Other critical gains detected included the CDK10, a cyclin-dependent kinase gene, which plays a role in cellular proliferation; GPR132, a member of the G-protein–coupled receptor (GPCR) family that slows cell proliferation and repairs damaged DNA; and SUMO3, which encodes a small ubiquitin-related modifier involved in nuclear transport, DNA replication, breakage and repair, mitosis, signal transduction, and the cellular SUMOylation system (28).

Prominent losses observed in HRS cells targeted genes with prominent roles in cell-cycle progression and tumor suppression. Both SKP2 (5p13.2) and ID4 (6p22.3) have been linked to the pathogenesis of lymphoma (29, 30). Inactivation by hypermethylation of ID4 in follicular and diffuse large cell lymphoma suggests a tumor suppressor (TSG) role in malignant lymphoma (30). Similar studies in other cancers also indicate that SPRY1 and NELL1, 2 receptor tyrosine kinase signaling genes (31, 32) also have tumor suppressor roles. SKP2 (F-box protein family) encodes a critical component of the ubiquitin ligase complex that regulates the G1/S transition of the cell cycle by mediating target proteins such as CDKN1B cyclin-dependent kinase inhibitor and RELB (33). Consistent with the phenotype of HRS cells, SKP2-knockout experiments in mice show a hypoproliferative phenotype with serious cellular defects including nuclear enlargement, polyploidy, an increased number of chromosomes, and chromosomal instability (34). Two other notable losses observed in HL were the TXNIP, a gene required for natural killer cell maturation, and PPARGC1A, a gene reported to antagonize NF-κB activation and the inflammatory gene expression of PTGS2 (prostaglandin endoperoxide synthase-2; also known as cyclooxygenase-2 or COX2; ref. 35).

A comparison between the HRS and normal GC molecular karyotypes showed a shared relationship with gains in 4 transcription factors associated with hematopoietic stem cells and early lymphocyte development. IRF8 is a key transcription factor that modulates the myeloid versus lymphoid lineage choice by hematopoietic stem cells via interactions with the ETS family proteins encoded by SPI1 or SPIB. SPI1 encodes an ETS domain transcription factor (PU.1) that activates gene expression during myeloid and B-lymphoid cell development, whereas SPIB acts as a lymphoid-specific enhancer to promote development of plasmacytoid dendritic cells or natural IFN-producing cells (36, 37). Specifically within the GC, IRF8 is part of the transcriptional network governing B-cell lineage differentiation including the transcriptional regulation of activation-induced cytidine deaminase and BCL6 (37). Finally, ICOSLG is a ligand for the T-cell–specific cell surface receptor ICOS, and the product of this gene acts as a costimulatory signal for lymphoid proliferation and cytokine secretion to mediate local tissue responses to inflammatory conditions. Of interest, ICOS is highly expressed on tonsillar T cells, which are closely associated with B cells in the apical light zone of GC (38). Taken together, the commonality of CNAs observed between GC and HRS cells substantiates a hematopoietic stem cell–derived progenitor cell relationship (4, 39).

A further comparison of the GC and aCGH HL results defined a subset of 151 CNAs consistently associated with HL. The 10 most frequent gene alterations in HL localized to 3 chromosomal regions and involved genes associated with inflammation, cell-cycle progression, and NF-κB/PI3K signaling: losses of AMACR, GDNF, SKP2 at 5p13, gain of SHH/7q36, and gains of ABL1, CDK9, ENG, LCN2, PTGES, and TSC1 at 9q34. As noted earlier, loss of SKP2 appears to be critical for the HRS hypoproliferative phenotype. The short-lived protein Cdk9 (the product of the CDK9 gene) is also regulated by the SCF (SKP2) ubiquitin ligase complex, suggesting loss of SKP2 and gains of CDK9 allows for a stronger association with TRAF2 and continued NF-κB activation (40). LCN2 (lipocalin 2, previously known as NGAL) mediates inflammatory responses and has binding sites for NF-κB, the hematopoietic transcription factors GATA-1 and SPI1 (PU.1), as well as the matrix metalloproteinases (MMP), which protect them from degradation and promote tumor cell invasion and metastasis (41). A link between MMP and LCN2 in this study is intriguing, as Steidl and colleagues recently reported that overexpression of MMP11 is consistently associated with a tumor-associated macrophage expression profile and shortened survival in HL (9). Finally, gains of PTGS2 (formerly known as COX2) are thought to enhance tumor angiogenesis, suppress antitumor immunity, and possibly play an important role for cytokine-driven PTGS2 NF-κB activation (42).

Our genome-wide discrimination analyses on CNAs revealed distinct clusters that correlated closely with the favorable and unfavorable IPS in our HL patients suggesting the possibility of different chemoresponsive genomic profiles at disease presentation. To test this hypothesis, we compared the 2 distinct HL subgroups, selected 20 genes from 10 variable CNAs, and uploaded the data into IPA tools 3.0 to explore gene function and interactions. The data suggested that a key difference between the CR and PR molecular karyotypes may be opposing calcium signaling and IL21/IL2 pathway interactions within the tumor microenvironment. The CR patients showed frequent gains of 16q13, a chromosome region known to house genes that regulate T-cell trafficking or NF-κB activation (CCL22, CX3CL1, CCL17, DOK4, and IL10). HRS cells are known to secrete CCL17 and CCL22 which attracts TH2 and regulatory T cells; report higher numbers of CD4+CD25 (high) FOXP3+ regulatory T cells are reported to predict for improved survival and prognosis in non–Hodgkin lymphoma (43). Further characterization of the tumor microenvironment may provide insight on the differential regulation of regulatory T cell function and survival in HL.

In contrast, PR HL showed frequent loss of genes at 4q27, 5p15, 11p14.3, 17p12 and gains at 13q31.2 and 19q13.31. Losses of 2 cytokine genes, IL2 and IL21 within 4q27, were evident in about 50% of the PR HRS cells. IL21 exhibits potent antitumor responses, is homologous to IL2, and shares the common cytokine receptor γ chain of the IL2 receptor family (44). Although IL21 expression is usually restricted to subsets of CD4+ T cells and natural killer cells (44), two recent studies describe aberrant expression of IL21 activating STAT3 and STAT5 in HL continuous cell lines and a subset of primary HRS cells (3, 45). Moreover, Scheeren and colleagues (45) have shown that constitutive activation of STAT5 by IL21 results in activated NF-κB activity in immortalized B cells that have lost their B-cell phenotype and therefore mimic HRS cells, supporting a link between STAT5 and NF-κB in a subset of HL.

Conversely, the loss of IL21 in PR HL, as seen in our study, may prompt additional mutations in HRS cells to maintain cell survival in a cytokine-independent manner. One such complementary mutation may be gain of BCL3/RELB at 19q13.3. Previous to this report, only rare primary cases of HL (∼15%) showed BCL3 gains or translocations in HRS cells (15, 46); however, gains of BCL3 appear to be more common in HL cell lines (12). In this study, we show that BCL3 gains are common in PR HL but infrequent findings in CR HL (∼50% vs. 15%, respectively). Of interest, BCL3 not only regulates the NF-κB signaling pathway by modulating the DNA-binding activity of NF-κB transcription factors involved in apoptosis and cell growth but also inhibits p53-induced apoptosis (47). Because IL21 is a potent inducer of BCL3 (48), the absence of IL21 may trigger BCL3 gene amplification as a mechanism for HRS cell survival. Interestingly, gains or overexpression of BCL3 have been associated with poor prognosis in other lymphoid malignancies including HTLV-1–induced ATLL and multiple myeloma (48, 49). Finally, losses within chromosome band 17p12 could possibly result in genomic rearrangements leading to overexpression of CD68, a gene housed within the neighboring chromosome band 17p13.1. Increased numbers of CD68+ macrophages in HL have been correlated with short disease-free survival and a high likelihood of relapse post–autologous stem cell transplantation (9). Based on our aCGH HL findings, we propose HL samples without gains in 16q13 and losses of IL21/IL2 at 4q27, 11p14.3/SLC17A6 and/or 17p12 and possibly concurrent gains of BCL3 portend a chemorefractory HL phenotype. A preliminary evaluation of the model appeared to support our hypothesis and the finding of 16p13 gains in the CR samples suggests that genes housed in this region may predict for a chemosensitive HL phenotype. More detailed explorations of these loci are needed to identify precise and reproducible prognostic genetic markers in HL.

Even though BAC aCGH technology is robust and provides high signal-to-noise ratios when working with challenging samples, such as HRS cells, captured by LCM from FFPE tissue, there are some limitations to its use for development of new clinical prognostic assays in HL. Because aCGH does not detect inactivating mutations and our selection criteria precluded the detection of small deletions, loss of the small approximately 16-kb TNFAIP3 (or A20) at 6q23.3 was not detected in this study. We also did not observe amplification of the 42-kb REL gene within 2p16.1 as reported by others in about 30% to 40% of classical HL (5, 46); however, gains within chromosome band 19q13.3 which houses both BCL3 and RELB were observed. Most likely, our strict inclusion criteria of investigating only gains and losses observed in more than 35% of samples tested and within a 5-BAC region (>700 kb) precluded detection of the small gene gains and losses, suggesting the necessity for higher resolution microarray platforms in future studies.

Clarification of the genetic alterations and molecular mechanisms that underlie the oncogenic HRS cell are likely to improve our current understanding of HL. In this study, we identified DNA CNAs in HL by aCGH using DNA extracted from microdissected CD30+ HRS cells, showing the potential utility of this approach for defining prognostic parameters. Despite the limited sample size, this HL investigation is one of the largest to date analyzing the molecular karyotypes of well-characterized, primary diagnostic HL samples with long-term follow-up. The molecular data showed a clear relationship between CNAs of key genes associated with the NF-κB signaling pathway, the tumor microenvironment, cell-cycle regulation and apoptosis; however, the functional consequences of many of these genetic lesions and their complex interactions in HL tumorigenesis remain to be elucidated. Finally, the genetic differences observed between the chemosensitive and PR HL molecular karyotypes, in this study and others (5–7, 9), suggest that important pathogenetic and prognostic differences are present at disease onset. Given the paucity of PR HL tumor samples available for study and the rarity of HRS cells within these lesions, it is imperative that we unite our efforts for the collection of retrospective HL samples with long-term follow-up data for ongoing studies to confirm, refine, and define relevant biomarkers in PR HL. Further studies exploiting these genetic lesions and defining their role in the pathogenesis of HL will be essential in developing new targeted therapeutic approaches.

Added note in support of our findings

During the manuscript review process, an interesting and similar aCGH in classical HL study using laser microdissected HRS cells was published by Steidl and colleagues (50). Both studies agree that gains and losses of genes involved in the NF-κB signaling pathway are recurrent findings in classical HL. Unlike our study, Steidl and colleagues (50) report on the significance of gains and overexpression of the MDR gene ABCC1 in their treatment failure cohort, which consisted of lymph node samples collected from both treatment-naive and treatment-exposed HL patients. We did not find a correlation between gains of ABCC1 and drug resistance; however, our population consisted only of newly diagnosed pretreatment samples and we used a slightly different method (FFPE not frozen samples). Although our study size is smaller, we only utilized samples collected at HL disease presentation (no relapsed samples) to eliminate bias associated with treatment for our comparison of chemosensitive versus PR HL. This more homogenous patient population in combination with stringent gain and loss threshold requirements makes our data equally compelling.

M.L. Slovak is currently an employee of Quest Diagnostics. The other authors declared no potential conflicts of interest.

The authors thank Aaron Theisen for his critical reading of the manuscript and helpful suggestions.

This work was supported in part by a Lymphoma SPORE Developmental Research Award, Cancer Center grants at Roswell Park Cancer Institute (NCI CA016056) and the City of Hope (NCI CA33572 and CA-30206), The Tim Nesvig Lymphoma Research Fund, and a private donation in memory of Pearl Ruttenberg.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Diehl
V
,
Thomas
RK
,
Re
D
. 
Part II: Hodgkin's lymphoma–diagnosis and treatment
.
Lancet Oncol
2004
;
5
:
19
26
.
2.
Schwering
I
,
Bräuninger
A
,
Klein
U
,
Jungnickel
B
,
Tinguely
M
,
Diehl
V
, et al
Loss of the B-lineage-specific gene expression program in Hodgkin and Reed-Sternberg cells of Hodgkin lymphoma
.
Blood
2003
;
101
:
1505
12
.
3.
Lamprecht
B
,
Kreher
S
,
Anagnostopoulos
I
,
Jöhrens
K
,
Monteleone
G
,
Jundt
F
, et al
Aberrant expression of the Th2 cytokine IL-21 in Hodgkin lymphoma cells regulates STAT3 signaling and attracts Treg cells via regulation of MIP-3alpha
.
Blood
2008
;
112
:
3339
47
.
4.
Kuppers
R
. 
The biology of Hodgkin's lymphoma
.
Nat Rev Cancer
2009
;
9
:
15
27
.
5.
Joos
S
,
Menz
CK
,
Wrobel
G
,
Siebert
R
,
Gesk
S
,
Ohl
S
, et al
Classical Hodgkin lymphoma is characterized by recurrent copy number gains of the short arm of chromosome 2
.
Blood
2002
;
99
:
1381
7
.
6.
Martín-Subero
JI
,
Gesk
S
,
Harder
L
,
Sonoki
T
,
Tucker
PW
,
Schlegelberger
B
, et al
Recurrent involvement of the REL and BCL11A loci in classical Hodgkin lymphoma
.
Blood
2002
;
99
:
1474
7
.
7.
Schmitz
R
,
Hansmann
ML
,
Bohle
V
,
Martin-Subero
JI
,
Hartmann
S
,
Mechtersheimer
G
, et al
TNFAIP3 (A20) is a tumor suppressor gene in Hodgkin lymphoma and primary mediastinal B cell lymphoma
.
J Exp Med
2009
;
206
:
981
9
.
8.
Horning
S
,
Fanale
M
,
deVos
S
,
Borchmann
P
,
Illidge
T
,
Engert
A
, et al
Defining a population of Hodgkin lymphoma patients for novel therapeutics: an international effort [abstract]
.
Ann Oncol
2008
;
20
:
118
.
9.
Steidl
C
,
Lee
T
,
Shah
SP
,
Farinha
P
,
Han
G
,
Nayar
T
, et al
Tumor-associated macrophages and survival in classic Hodgkin's lymphoma
.
N Engl J Med
2010
;
362
:
875
85
.
10.
Giefing
M
,
Arnemann
J
,
Martin-SubSteidl
C
,
Lee
T
,
Shah
SP
,
Farinha
P
, et al
Identification of candidate tumour suppressor gene loci for Hodgkin and Reed-Sternberg cells by characterisation of homozygous deletions in classical Hodgkin lymphoma cell lines
.
Br J Haematol
2008
;
142
:
916
24
.
11.
Feys
T
,
Poppe
B
,
De Preter
K
,
Van Roy
N
,
Verhasselt
B
,
De Paepe
P
, et al
A detailed inventory of DNA copy number alterations in four commonly used Hodgkin's lymphoma cell lines
.
Haematologica
2007
;
92
:
913
20
.
12.
Mathas
S
,
Jöhrens
K
,
Joos
S
,
Lietz
A
,
Hummel
F
,
Janz
M
, et al
Elevated NF-kappaB p50 complex formation and Bcl-3 expression in classical Hodgkin, anaplastic large-cell, and other peripheral T-cell lymphomas
.
Blood
2005
;
106
:
4287
93
.
13.
Nowak
NJ
,
Miecznikowski
J
,
Moore
SR
,
Gaile
D
,
Bobadilla
D
,
Smith
DD
, et al
Challenges in array comparative genomic hybridization for the analysis of cancer samples
.
Genet Med
2007
;
9
:
585
95
.
14.
Moore
SR
,
Persons
DL
,
Sosman
JA
,
Bobadilla
D
,
Bedell
V
,
Smith
DD
, et al
Detection of copy number alterations in metastatic melanoma by a DNA fluorescence in situ hybridization probe panel and array comparative genomic hybridization: a southwest oncology group study (S9431)
.
Clin Cancer Res
2008
;
14
:
2927
35
.
15.
Hartmann
S
,
Martin-Subero
JI
,
Gesk
S
,
Hüsken
J
,
Giefing
M
,
Nagel
I
, et al
Detection of genomic imbalances in microdissected Hodgkin and Reed-Sternberg cells of classical Hodgkin's lymphoma by array-based comparative genomic hybridization
.
Haematologica
2008
;
93
:
1318
26
.
16.
Hasenclever
D
,
Diehl
V
. 
A prognostic score for advanced Hodgkin's disease. International Prognostic Factors Project on Advanced Hodgkin's Disease
.
N Engl J Med
1998
;
339
:
1506
14
.
17.
Bedell
V
,
Forman
SJ
,
Gaal
K
,
Pullarkat
V
,
Weiss
LM
,
Slovak
ML
. 
Successful application of a direct detection slide-based sequential phenotype/genotype assay using archived bone marrow smears and paraffin embedded tissue sections
.
J Mol Diagn
2007
;
9
:
589
97
.
18.
Fury
W
,
Batliwalla
F
,
Gregersen
PK
,
Li
W
. 
Overlapping probabilities of top ranking gene lists, hypergeometric distribution, and stringency of gene selection criterion
.
Conf Proc IEEE Eng Med Biol Soc
2006
;
1
:
5531
4
.
19.
Bankhead
A
,
Sach
I
,
Ni
C
,
LeMeur
N
,
Kruger
M
,
Ferrer
M
, et al
Knowledge based identification of essential signaling from genome-scale siRNA experiments
.
BMC Syst Biol
2009
;
3
:
80
.
20.
Dutton
A
,
Reynolds
GM
,
Dawson
CW
,
Young
LS
,
Murray
PG
. 
Constitutive activation of phosphatidyl-inositide 3 kinase contributes to the survival of Hodgkin's lymphoma cells through a mechanism involving Akt kinase and mTOR
.
J Pathol
2005
;
205
:
498
506
.
21.
Jundt
F
,
Acikgöz
O
,
Kwon
SH
,
Schwarzer
R
,
Anagnostopoulos
I
,
Wiesner
B
, et al
Aberrant expression of Notch1 interferes with the B-lymphoid phenotype of neoplastic B cells in classical Hodgkin lymphoma
.
Leukemia
2008
;
22
:
1587
94
.
22.
Zhang
L
,
Blackwell
K
,
Shi
Z
,
Habelhah
H
. 
The RING domain of TRAF2 plays an essential role in the inhibition of TNFalpha-induced cell death but not in the activation of NF-kappaB
.
J Mol Biol
2010
;
396
:
528
39
.
23.
Georgakis
GV
,
Li
Y
,
Rassidakis
GZ
,
Medeiros
LJ
,
Mills
GB
,
Younes
A
. 
Inhibition of the phosphatidylinositol-3 kinase/Akt promotes G1 cell cycle arrest and apoptosis in Hodgkin lymphoma
.
Br J Haematol
2006
;
132
:
503
11
.
24.
Jeong
SJ
,
Pise-Masison
CA
,
Radonovich
MF
,
Park
HU
,
Brady
JN
. 
Activated AKT regulates NF-kappaB activation, p53 inhibition and cell survival in HTLV-1-transformed cells
.
Oncogene
2005
;
24
:
6719
28
.
25.
Hayashi
H
,
Kume
T
. 
Foxc transcription factors directly regulate Dll4 and Hey2 expression by interacting with the VEGF-Notch signaling pathways in endothelial cells
.
PLoS One
2008
;
3
:
e2401
.
26.
Nakashima
H
,
Nakamura
M
,
Yamaguchi
H
,
Yamanaka
N
,
Akiyoshi
T
,
Koga
K
, et al
Nuclear factor-kappaB contributes to hedgehog signaling pathway activation through sonic hedgehog induction in pancreatic cancer
.
Cancer Res
2006
;
66
:
7041
9
.
27.
Kasperczyk
H
,
Baumann
B
,
Debatin
KM
,
Fulda
S
. 
Characterization of sonic hedgehog as a novel NF-kappaB target gene that promotes NF-kappaB-mediated apoptosis resistance and tumor growth in vivo
.
FASEB J
2009
;
23
:
21
33
.
28.
Galanty
Y
,
Belotserkovskaya
R
,
Coates
J
,
Polo
S
,
Miller
KM
,
Jackson
SP
. 
Mammalian SUMO E3-ligases PIAS1 and PIAS4 promote responses to DNA double-strand breaks
.
Nature
2009
;
462
:
935
9
.
29.
Seki
R
,
Ohshima
K
,
Fujisaki
T
,
Uike
N
,
Kawano
F
,
Gondo
H
, et al
Prognostic significance of S-phase kinase-associated protein 2 and p27kip1 in patients with diffuse large B-cell lymphoma: effects of rituximab
.
Ann Oncol
2010
;
21
:
833
41
.
30.
Hagiwara
K
,
Nagai
H
,
Li
Y
,
Ohashi
H
,
Hotta
T
,
Saito
H
. 
Frequent DNA methylation but not mutation of the ID4 gene in malignant lymphoma
.
J Clin Exp Hematop
2007
;
47
:
15
8
.
31.
Jin
Z
,
Mori
Y
,
Yang
J
,
Sato
F
,
Ito
T
,
Cheng
Y
, et al
Hypermethylation of the nel-like 1 gene is a common and early event and is associated with poor prognosis in early-stage esophageal adenocarcinoma
.
Oncogene
2007
;
26
:
6332
40
.
32.
Kwabi-Addo
B
,
Wang
J
,
Erdem
H
,
Vaid
A
,
Castro
P
,
Ayala
G
, et al
The expression of Sprouty1, an inhibitor of fibroblast growth factor signal transduction, is decreased in human prostate cancer
.
Cancer Res
2004
;
64
:
4728
35
.
33.
Nakayama
KI
,
Nakayama
K
. 
Regulation of the cell cycle by SCF-type ubiquitin ligases
.
Semin Cell Dev Biol
2005
;
16
:
323
33
.
34.
Nakayama
K
,
Nagahama
H
,
Minamishima
YA
,
Matsumoto
M
,
Nakamichi
I
,
Kitagawa
K
, et al
Targeted disruption of Skp2 results in accumulation of cyclin E and p27(Kip1), polyploidy and centrosome overduplication
.
EMBO J
2000
;
19
:
2069
81
.
35.
Ackerman
WEt
,
Zhang
XL
,
Rovin
BH
,
Kniss
DA
. 
Modulation of cytokine-induced cyclooxygenase 2 expression by PPARG ligands through NFkappaB signal disruption in human WISH and amnion cells
.
Biol Reprod
2005
;
73
:
527
35
.
36.
Schmidlin
H
,
Diehl
SA
,
Nagasawa
M
,
Scheeren
FA
,
Schotte
R
,
Uittenbogaart
CH
, et al
Spi-B inhibits human plasma cell differentiation by repressing BLIMP1 and XBP-1 expression
.
Blood
2008
;
112
:
1804
12
.
37.
Wang
H
,
Lee
CH
,
Qi
C
,
Tailor
P
,
Feng
J
,
Abbasi
S
, et al
IRF8 regulates B-cell lineage specification, commitment, and differentiation
.
Blood
2008
;
112
:
4028
38
.
38.
Hutloff
A
,
Dittrich
AM
,
Beier
KC
,
Eljaschewitsch
B
,
Kraft
R
,
Anagnostopoulos
I
, et al
ICOS is an inducible T-cell co-stimulator structurally and functionally related to CD28
.
Nature
1999
;
397
:
263
6
.
39.
Klein
U
,
Dalla-Favera
R
. 
Germinal centres: role in B-cell physiology and malignancy
.
Nat Rev Immunol
2008
;
8
:
22
33
.
40.
MacLachlan
TK
,
Sang
N
,
De Luca
A
,
Puri
PL
,
Levrero
M
,
Giordano
A
. 
Binding of CDK9 to TRAF2
.
J Cell Biochem
1998
;
71
:
467
78
.
41.
Yan
L
,
Borregaard
N
,
Kjeldsen
L
,
Moses
MA
. 
The high molecular weight urinary matrix metalloproteinase (MMP) activity is a complex of gelatinase B/MMP-9 and neutrophil gelatinase-associated lipocalin (NGAL). Modulation of MMP-9 activity by NGAL
.
J Biol Chem
2001
;
276
:
37258
65
.
42.
Ackerman
WEt
,
Summerfield
TL
,
Vandre
DD
,
Robinson
JM
,
Kniss
DA
. 
Nuclear factor-kappa B regulates inducible prostaglandin E synthase expression in human amnion mesenchymal cells
.
Biol Reprod
2008
;
78
:
68
76
.
43.
Ke
X
,
Wang
J
,
Li
L
,
Chen
IH
,
Wang
H
,
Yang
XF
. 
Roles of CD4+CD25(high) FOXP3 +Tregs in lymphomas and tumors are complex
.
Front Biosci
2008
;
13
:
3986
4001
.
44.
Spolski
R
,
Leonard
WJ
. 
Interleukin-21: basic biology and implications for cancer and autoimmunity
.
Annu Rev Immunol
2008
;
26
:
57
79
.
45.
Scheeren
FA
,
Diehl
SA
,
Smit
LA
,
Beaumont
T
,
Naspetti
M
,
Bende
RJ
, et al
IL-21 is expressed in Hodgkin lymphoma and activates STAT5: evidence that activated STAT5 is required for Hodgkin lymphomagenesis
.
Blood
2008
;
111
:
4706
15
.
46.
Martin-Subero
JI
,
Wlodarska
I
,
Bastard
C
,
Picquenot
JM
,
Höppner
J
,
Giefing
M
, et al
Chromosomal rearrangements involving the BCL3 locus are recurrent in classical Hodgkin and peripheral T-cell lymphoma
.
Blood
2006
;
108
:
401
2
;
author reply 402–3
.
47.
Kashatus
D
,
Cogswell
P
,
Baldwin
AS
. 
Expression of the Bcl-3 proto-oncogene suppresses p53 activation
.
Genes Dev
2006
;
20
:
225
35
.
48.
Brenne
AT
,
Fagerli
UM
,
Shaughnessy
JD
 Jr
,
Våtsveen
TK
,
TB
,
Hella
H
, et al
High expression of BCL3 in human myeloma cells is associated with increased proliferation and inferior prognosis
.
Eur J Haematol
2009
;
82
:
354
63
.
49.
Kim
YM
,
Sharma
N
,
Nyborg
JK
. 
The proto-oncogene Bcl3, induced by Tax, represses Tax-mediated transcription via p300 displacement from the human T-cell leukemia virus type 1 promoter
.
J Virol
2008
;
82
:
11939
47
.
50.
Steidl
C
,
Telenius
A
,
Shah
SP
,
Farinha
P
,
Barclay
L
,
Boyle
M
, et al
Genome-wide copy number analysis of Hodgkin Reed-Sternberg cells identifies recurrent imbalances with correlations to treatment outcome
.
Blood
2010
;
116
:
418
27
.