Abstract
The pathogenesis of hepatosplenic T-cell lymphoma (HSTL), a rare entity mostly derived from γδ T cells and usually with a fatal outcome, remains largely unknown. In this study, HSTL samples (7γδ and 2αβ) and the DERL2 HSTL cell line were subjected to combined gene-expression profiling and array-based comparative genomic hybridization. Compared with other T-cell lymphomas, HSTL had a distinct molecular signature irrespective of TCR cell lineage. Compared with peripheral T-cell lymphoma, not otherwise specified and normal γδ T cells, HSTL overexpressed genes encoding NK-cell–associated molecules, oncogenes (FOS and VAV3), the sphingosine-1-phosphatase receptor 5 involved in cell trafficking, and the tyrosine kinase SYK, whereas the tumor-suppressor gene AIM1 (absent in melanoma 1) was among the most down-expressed. We found highly methylated CpG islands of AIM1 in DERL2 cells, and decitabine treatment induced a significant increase in AIM1 transcripts. Syk was present in HSTL cells and DERL2 cells contained phosphorylated Syk and were sensitive to a Syk inhibitor in vitro. Genomic profiles confirmed recurrent isochromosome 7q (n = 6/9) without alterations at the SYK and AIM1 loci. Our results identify a distinct molecular signature for HSTL and highlight oncogenic pathways that offer rationale for exploring new therapeutic options such as Syk inhibitors and demethylating agents.
Introduction
Hepatosplenic T-cell lymphoma (HSTL), originally described as hepatosplenic γδ T-cell lymphoma, is a rare lymphoma entity with a peculiar clinical presentation (hepatosplenomegaly without significant lymphadenopathy) and pathologic features (intrasinusal/sinusoidal infiltration by neoplastic Τ cells in the BM, spleen, and liver).1-3 The disease occurs predominantly in young adults, in association with a setting of long-term immunosuppression in solid organ transplantation recipients or with prolonged antigenic stimulation.4 Cases have also been reported in children treated by azathioprine and infliximab for Crohn disease.5 Whereas most HSTL cases are derived from the γδ subset, a few similar cases with an αβ phenotype have also been described,6,7 and the simplified designation “hepatosplenic T-cell lymphoma” was favored in the latest World Health Organization (WHO) classification.8 HSTL is associated with a recurrent isochromosome 7q and, less often, trisomy 8,9 but its pathogenesis remains largely unknown.
Despite relatively innocuous cytology, the disease is highly aggressive, with an almost certain fatal outcome and a median overall survival barely exceeding 1 year.4 Occasional long-time survivors have been reported, and few patients respond to cytarabine or deoxycoformycin.4,10 Therapeutic strategies curative in a significant proportion of other aggressive subtypes of lymphoma have proven ineffective in HSTL, and efficient treatment modalities remain to be defined.
Over the past years, genome-wide molecular profiling studies have contributed significant insights into the pathobiology of several T-cell lymphoma entities11-14 and have provided informative data on the multiple molecular subgroups in peripheral T-cell lymphoma, not otherwise specified (PTCL-NOS).15,16 Data on HSTL remain scarce.13,17 In the present study, we analyzed a series of HSTL samples in relation to normal γδ cells, PTCL-NOS, and extranodal NK/T-cell lymphoma, nasal-type (NKTCL), another entity derived from cytotoxic lymphocytes of the innate immune system. The aims of the present study were: (1) to characterize the molecular signature of HSTL, (2) to identify potential candidate pathways relevant to pathogenesis, and (3) to search for biomarkers useful for diagnostic purposes or in the future targeted therapies.
Methods
Patient characteristics and tumor samples
Nine HSTL patients with high-quality RNA and/or DNA extracted from frozen tumor samples were selected for this study. All patients had spleen, liver, and BM involvement without lymphadenopathies. Three patients had been included in previous studies.4,9 The main clinical, phenotypical, and molecular characteristics are summarized in Table 1. The tumor samples comprised 6 splenic tissue samples and 3 cell suspensions (from spleen, BM, and blood), 2 of which were enriched in tumor cells (samples HSTL_01 and HSTL_09). All cases were reviewed by 3 hematopathologists (L.d.L, Y.H., and P.G.) and diagnosed according to WHO criteria.8 The tumor cells had a CD3+, CD2+, CD5−, TiA1+, GzmB− immunophenotype and were negative for EBV. TCR lineage was determined by immunohistochemistry and/or flow cytometry for TCRβ and TCRδ chain expression and by GC-clamp multiplex PCR for TCRγ and/or δ chain rearrangements (the PCRδ–denaturing gradient gel electrophoresis procedure).18 In total, 7 patients with a δTCR1+, βF1− immunophenotype and/or a biallelic rearrangement of the TCRδ chain18,19 were classified as γδ HSTL and 2 patients with a δTCR1−, βF1+ phenotype as αβ HSTL. Four of 7 investigated patients had isochromosome 7q. Twelve additional HSTL cases were selected for validations (10 formalin-fixed tissues for immunohistochemistry and 2 frozen samples for RT-PCR analyses).
EBI ID . | Age, y . | Sex . | Diagnosis . | Organ . | Type of sample . | CD3 . | CD2 . | CD5 . | CD7 . | CD4 . | CD8 . | βF1 . | δTCR . | CD56 . | TiA1 . | GzmB . | EBV* . | T clone . | Cytogenetics . |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
HSTL_02 | 41 | F | HSTLγδ | Spleen | Tissue | + | + | - | - | - | - | - | +† | + | + | - | - | + | iso7q (+) |
HSTL_03 | 73 | F | HSTLγδ | Spleen | Tissue | + | + | - | - | - | + | - | + | + | + | + | - | + | iso7q (−) |
HSTL_04 | 38 | M | HSTLγδ | Spleen | Tissue | + | ND | - | - | - | - | + | +‡ | - | + | - | - | + | iso7q (+) trisomy 8 |
HSTL_05 | 56 | M | HSTL αβ | Spleen | Tissue | + | ND | - | ND | - | - | + | - | - | - | - | - | + | iso7q (+) |
HSTL_06 | 53 | F | HSTL αβ | Spleen | Tissue | + | - | - | + | - | - | + | - | + | + | - | - | + | iso7q (−) |
HSTL_08 | 24 | F | HSTLγδ | Blood | Cell | + | ND | - | + | - | - | - | + | + | + | ND | ND | + | |
HSTL_09 | 36 | M | HSTLγδ | Blood | Cell (purified CD56+, 77%) | + | + | - | + | - | - | - | + | + | + | - | ND | + | iso7q (+) |
HSTL_01 | 86 | F | HSTLγδ | Spleen | Cell (purified CD56+, 84%) | + | + | -/+ | ND | - | - | - | + | + | + | - | - | + | iso7q (−) |
HSTL_10 | 67 | M | HSTLγδ | Spleen | Tissue | + | + | - | + | - | - | - | +‡ | - | + | - | - | + | Iso7q (−) |
HSTL_07 DERL2 | 30 | M | HSTLγδ | Blood | Cell line | + | + | - | + | - | + | + | + | - | Iso7q (+) |
EBI ID . | Age, y . | Sex . | Diagnosis . | Organ . | Type of sample . | CD3 . | CD2 . | CD5 . | CD7 . | CD4 . | CD8 . | βF1 . | δTCR . | CD56 . | TiA1 . | GzmB . | EBV* . | T clone . | Cytogenetics . |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
HSTL_02 | 41 | F | HSTLγδ | Spleen | Tissue | + | + | - | - | - | - | - | +† | + | + | - | - | + | iso7q (+) |
HSTL_03 | 73 | F | HSTLγδ | Spleen | Tissue | + | + | - | - | - | + | - | + | + | + | + | - | + | iso7q (−) |
HSTL_04 | 38 | M | HSTLγδ | Spleen | Tissue | + | ND | - | - | - | - | + | +‡ | - | + | - | - | + | iso7q (+) trisomy 8 |
HSTL_05 | 56 | M | HSTL αβ | Spleen | Tissue | + | ND | - | ND | - | - | + | - | - | - | - | - | + | iso7q (+) |
HSTL_06 | 53 | F | HSTL αβ | Spleen | Tissue | + | - | - | + | - | - | + | - | + | + | - | - | + | iso7q (−) |
HSTL_08 | 24 | F | HSTLγδ | Blood | Cell | + | ND | - | + | - | - | - | + | + | + | ND | ND | + | |
HSTL_09 | 36 | M | HSTLγδ | Blood | Cell (purified CD56+, 77%) | + | + | - | + | - | - | - | + | + | + | - | ND | + | iso7q (+) |
HSTL_01 | 86 | F | HSTLγδ | Spleen | Cell (purified CD56+, 84%) | + | + | -/+ | ND | - | - | - | + | + | + | - | - | + | iso7q (−) |
HSTL_10 | 67 | M | HSTLγδ | Spleen | Tissue | + | + | - | + | - | - | - | +‡ | - | + | - | - | + | Iso7q (−) |
HSTL_07 DERL2 | 30 | M | HSTLγδ | Blood | Cell line | + | + | - | + | - | + | + | + | - | Iso7q (+) |
EBV investigation by in situ hybridization with EBER probes and/or immunohistochemistry with LMP-1 Ab.
HSTL_02 with a γδ HSTL phenotype at diagnosis disclosed a TCR-silent phenotype (βF1−δTCR1−) at the time of analysis.
δ chain rearrangement determination using a GC-clamp multiplex PCR–δ-DGGE procedure.
ND indicates not determined.
The study was approved by the institutional review board Comité de Protection des Personnes Ile de France IX, Créteil, France (number 08-009; June 5, 2008).
Cell lines and normal γδ T cells
Eight samples of normal peripheral blood activated γδ T cells (6 from Correia20 and 2 from Zhang21 ) and 3 samples of normal resting γδ T cells (from Correia20 ) were used for comparison. Normal γδ T cells sorted using a magnetic anti-TCRγ/δ microbeads kit (Miltenyi Biotec) from 3 spleens removed for benign conditions were also obtained for quantitative RT-PCR (qRT-PCR) validation. DERL2, a γδ cell line derived from HSTL,22 was cultured in IMDM supplemented with 2mM l-glutamine and 20% heat-inactivated human serum (Invitrogen) in the presence of recombinant human IL-2 (100 U/mL; Laboratoires Chiron). Centrifuged pellets of DERL2 cells were fixed in ethanol to construct paraffin-embedded blocks.14
Microarray procedures
GEP
For gene-expression profiling (GEP), HG-U133 Plus Version 2.0 Affymetrix array data were obtained from 9 HSTL patients and the DERL-2 cell line (Table 1) and deposited in ArrayExpress under accession number E-MTAB-638. HG-U133 Plus Version 2.0 profiles of 42 PTCL-NOS, 23 angioimmunoblastic T-cell lymphoma (AITL), 23 NKTCL, and 3 HSTL patients described previously by our group (E-TABM-783 and E-TABM-702, n = 40)11,14 and others (GSE19067 and GSE6338, n = 51)12,13 were used for analyses. Control samples (8 activated γδ T cells, 3 resting γδ T cells, 18 B cells, 6 spleens, and 9 lymph nodes) were collected by us (n = 6) or were from public sources (E-MEX-1601, GSE13906, GSE12195, GSE15271, and GSE7307).20,21,23,24 Affymetrix raw data of all samples were normalized in batch using the RMA algorithm. The clustering analysis of the 50 tumor samples from our series was performed as described previously.14 We used moderate t tests to identify genes differentially expressed between 2 groups of samples and area under the curve criteria to identify discriminant genes. KEGG and Biocarta pathways (and related genes) were obtained from ftp://ftp.genome.ad.jp/pub/kegg/pathways/hsa and http://www.biocarta.com, respectively. Pathways were ranked according to the scores of 3 different algorithms, including GlobalTest, SAM-GS, and GSA. Details are provided in supplemental Methods (available on the Blood Web site; see the Supplemental Materials link at the top of the online article).
aCGH
DNAs extracted from 7 HSTL patients and the DERL2 cell line were hybridized on Agilent Sureprint CGH array 4 × 180 k covering more than 170 000 coding and noncoding human sequences. Array-based comparative genomic hybridization (aCGH) raw Cy3/Cy5 intensities were quantile normalized independently for each sample; smoothing was performed using the tiling Array R package. Copy number status (gain/loss) determination in the smoothed segments was done as described previously.14 Details are given in supplemental Methods.
qRT-PCR analysis of candidate genes
The expression level of candidate genes (FOS, FOSB, AIM1, RHOB, ABCB1, VAV3, KIR3DS1, and S1PR5) identified in the GEP analyses was determined by TaqMan qRT-PCR (Applied Biosystems) in 11 primary HSTL tumors (including the 9 cases submitted to GEP analysis), the DERL2 cell line, and 3 normal splenic γδ T-cell samples as a control. All primers and probes were purchased from Applied Biosystems and the gene expression was measured using Mastercycler ep realplex2S system (Eppendorf). Quantifications were done in duplicate and mean values and SDs were calculated for each transcript, as described previously.25
Immunohistochemical analysis of selected candidate genes
Immunohistochemistry was performed on deparaffinized tissue sections using a standard indirect avidin-biotin immunoperoxidase method. After appropriate antigen retrieval, sections were stained for FosB (Cell Signaling Technology); CD56, CD163, and GSTP1 (Novocastra-Leica); β-catenin (BD Biosciences); Bcl-10 (Zymed); ICAM-1 (Atlas Antibodies); VCAM-1 and Syk (Santa Cruz Biotechnology); and Blimp-1 (BioLegend). For Granzyme H (4G5), a tyramide signal amplification system was applied, as described previously.14 Adequate controls were included.
Methylation analysis of AIM1
Two micrograms of genomic DNA was treated with bisulfite using the Epitect Bisulfite Kit (QIAGEN) according to the manufacturer's protocol. Nested bisulfite-specific PCR for the 2 CpG islands of AIM1 (absent in melanoma 1) was performed on DERL2 cells DNA (details, primers, and PCR programs used are given in supplemental Methods). PCR products were cloned in plasmid using the TOPO TA cloning kit for sequencing (Invitrogen). Ten clones for each CpG island were screened and sequenced on the ABI 3130 × 1 genetic analyzer (Applied Biosystems).
5-Aza-2′-deoxycytidine treatment of DERL2 cells
DERL2 cells were treated with the demethylating agent 5-Aza-2′-deoxycytidine (decitabine; Merck Chemicals) at 10 or 40μM for 96 hours (decitabine was added every 24 hours) alone or in combination with 500nM histone deacetylase inhibitor trichostatin A (TSA; Sigma-Aldrich). In the combined treatment (decitabine and TSA), cells were treated every 24 hours with 10 or 40μM decitabine for 96 hours, followed by 500nM TSA for the last 24 hours. Total RNA was extracted and qRT-PCR was performed to evaluate the reexpression of the demethylated gene AIM1 after decitabine treatment. In addition, the number of apoptotic cells was determined by annexin V+ 7-amino-actinomycin D–positive (7AAD+) cells using a cyan flow cytometer (Beckman Coulter).
Immunoblot analysis
Total cell proteins extracted from normal blood T cells purified after CD2 magnetic beads selection and activated or not by CD3 (1 μg/mL), CD28 (1 μg/mL), and IL-2 (100 U/mL); overnight serum deprived DERL2 cells with or without CD3, CD28, and IL-2 activation and overnight serum deprived SUDHL4 cells with or without BCR stimulation (anti-Fab′2, 4 μg/mL) were analyzed for the expression of Syk, phospho-Syk (Tyr525/526; C87C1; Cell Signaling technology), and β-actin (Sigma-Aldrich). Immunoblotting was as described previously.25
Apoptosis assay with a Syk inhibitor
Normal activated γδ T cells were obtained as described previously,26 and DERL2 cells were plated at 1 000 000/well and incubated for 48 hours in the absence or presence of Syk inhibitor II (Merck Chemicals) at 2 different concentrations (40 and 60μM). The number of apoptotic cells was determined by the number of annexin V+7AAD+ cells. For the normal activated γδ T cells, TCRVγ9 apoptotic cells were determined by measuring the percentage of TCRVγ9+ and 7AAD+ cells.
Results
HSTL as a molecularly distinct entity irrespective of TCR γδ or αβ lineage
Unsupervised consensus clustering of the expression profiles of 9 HSTL tumors and 41 additional T-cell lymphoma samples, including PTCL-NOS, AITL, and NKTCL, demonstrated a clear separation between HSTL and other T-cell lymphomas entities. The HSTL cluster comprised 2 subgroups consisting of HSTL tissue samples and HSTL cells (sorted cells and cell line; Figure 1). Interestingly, αβ and γδ HSTL clustered together. One nonhepatosplenic γδ T-cell lymphoma case (PTCL_20) with nodal presentation clustered within the NKTCL branch.
On supervised analysis, αβ and γδ HSTL appeared to have a highly similar gene signature. As illustrated in supplemental Figure 1, a similar level of transcripts of NK-cell–associated molecules, including killer immunoglobulin-like receptors (KIRs), was found in HSTL, irrespective of αβ or γδ lineage. In addition to the genes coding for the γδ and αβ TCR, which, as expected, were the most differentially expressed in both groups, only a small proportion of genes were overexpressed (eg, LCK, SYT11, BCL2L11, and KLRC1) or down-expressed (eg, CR1, TMPRSS3, and ITGA9) in γδ HSTL compared with αβ HSTL (supplemental Table 1). Comparison of HSTL GEPs with those of either PTCL-NOS or NKTCL revealed similar proportions of differentially expressed genes (see “H1 Proportion” in supplemental Methods) of 37% and 39%, respectively. However, in unsupervised clustering analysis, HSTL profiles appeared to be closer to that of NKTCL, another extranodal entity derived from cytotoxic cells, than to that of PTCL-NOS and AITL.
Compared with both PTCL-NOS and NKTCL, several T cell–associated transcripts encoding the γδ TCR subunit were among the top genes overexpressed in HSTL. Two of the genes found to be up-regulated in HSTL were S1PR5, a gene known to be involved in the homing of NK cells into the spleen,27 and ABCB1, a gene encoding for the p-glycoprotein multidrug transporter (MDR1). These findings were confirmed by qRT-PCR demonstrating high S1PR5 and ABCB1 transcripts in HSTL compared with PTCL-NOS and NKTCL (Figure 2). Genes encoding KIRs and other NK-cell–associated molecules (eg, NCAM1 and CD244) were specifically overexpressed compared with PTCL-NOS (Tables 2 and 3).
Genes overexpressed in HSTL compared with PTCL-NOS . | |
---|---|
Gene class and specific genes . | Fold change . |
NK-cell–associated molecules | |
KLRC3 | 73.95 |
KIR2DS2 | 33.03 |
KLRD1 | 25.39 |
CD16 | 11.83 |
KIR3DL2 | 21.58 |
KLRB1 | 17.52 |
KIR3DS1 | 16.88 |
KIR2DL2 | 20.35 |
KIR3DL1 | 19.85 |
KIR2DL3 | 16.35 |
KIR3DL3 | 18.95 |
KIR2DL1 | 17.42 |
KLRC4 | 14.41 |
NCAM | 18.04 |
CD244 | 10.65 |
KIR2DS5 | 11.59 |
KIR2DL5A | 12.02 |
KIR2DS1 | 9.47 |
KIR2DS3 | 9.05 |
KLRG1 | 9.39 |
KLRC1/KLRC2 | 10.37 |
KIR2DL3/KIR2DL4/KIR2DL5A | 6.04 |
CTSW (Cathepsin w) | 6.76 |
KLRK1 | 6.45 |
Microenvironment | |
Hemoglobin gamma | 17.62 |
Hemoglobin beta | 15.13 |
Hemoglobin alpha | 14.39 |
DEFA1/DEFA1B/DEFA3 | 13.64 |
CCL3 (MIP1) | 4.44 |
Chemokines | |
CXCL7 (PPBP) | 21.06 |
CXCL6 | 8.96 |
Multidrug resistance | |
ABCB1 (MDR1) | 11.72 |
ABCB1/ABCB4 | 6.74 |
Oncogenes | |
MYBL1 | 18.81 |
VAV3 | 3.86 |
Cell adhesion/cell-cell interaction | |
CD11d | 48.35 |
PGM5 | 6.90 |
TCR | |
TRA@/TRD@ | 13.22 |
TARP/TRGC2 | 10.94 |
Transcription factors | |
TCF21 | 12.87 |
PRDM16 | 7.86 |
WT1 | 3.70 |
Growth factors | |
IGFBP2 | 7.52 |
AREG (amphiregulin) | 8.96 |
PDGFD | 7.83 |
ECM interaction | |
ADAMTS17 | 12.93 |
Others | |
S1PR5 | 49.99 |
Genes overexpressed in HSTL compared with PTCL-NOS . | |
---|---|
Gene class and specific genes . | Fold change . |
NK-cell–associated molecules | |
KLRC3 | 73.95 |
KIR2DS2 | 33.03 |
KLRD1 | 25.39 |
CD16 | 11.83 |
KIR3DL2 | 21.58 |
KLRB1 | 17.52 |
KIR3DS1 | 16.88 |
KIR2DL2 | 20.35 |
KIR3DL1 | 19.85 |
KIR2DL3 | 16.35 |
KIR3DL3 | 18.95 |
KIR2DL1 | 17.42 |
KLRC4 | 14.41 |
NCAM | 18.04 |
CD244 | 10.65 |
KIR2DS5 | 11.59 |
KIR2DL5A | 12.02 |
KIR2DS1 | 9.47 |
KIR2DS3 | 9.05 |
KLRG1 | 9.39 |
KLRC1/KLRC2 | 10.37 |
KIR2DL3/KIR2DL4/KIR2DL5A | 6.04 |
CTSW (Cathepsin w) | 6.76 |
KLRK1 | 6.45 |
Microenvironment | |
Hemoglobin gamma | 17.62 |
Hemoglobin beta | 15.13 |
Hemoglobin alpha | 14.39 |
DEFA1/DEFA1B/DEFA3 | 13.64 |
CCL3 (MIP1) | 4.44 |
Chemokines | |
CXCL7 (PPBP) | 21.06 |
CXCL6 | 8.96 |
Multidrug resistance | |
ABCB1 (MDR1) | 11.72 |
ABCB1/ABCB4 | 6.74 |
Oncogenes | |
MYBL1 | 18.81 |
VAV3 | 3.86 |
Cell adhesion/cell-cell interaction | |
CD11d | 48.35 |
PGM5 | 6.90 |
TCR | |
TRA@/TRD@ | 13.22 |
TARP/TRGC2 | 10.94 |
Transcription factors | |
TCF21 | 12.87 |
PRDM16 | 7.86 |
WT1 | 3.70 |
Growth factors | |
IGFBP2 | 7.52 |
AREG (amphiregulin) | 8.96 |
PDGFD | 7.83 |
ECM interaction | |
ADAMTS17 | 12.93 |
Others | |
S1PR5 | 49.99 |
Genes down-regulated in HSTL compared with PTCL-NOS . | |
---|---|
Chemokines | |
CCL19 | 0.011 |
CXCL9 | 0.016 |
CCL18 | 0.041 |
CXCL10 | 0.044 |
Genes down-regulated in HSTL compared with PTCL-NOS . | |
---|---|
Chemokines | |
CCL19 | 0.011 |
CXCL9 | 0.016 |
CCL18 | 0.041 |
CXCL10 | 0.044 |
Gene class and specific genes . | Fold change . |
---|---|
TFH signature | |
CXCL13 | 0.057 |
CD200 | 0.165 |
ICOS | 0.158 |
CXCR5 | 0.246 |
Immunomodulation | |
IDO1 | 0.074 |
IL4I1 | 0.160 |
Tumor suppressor | |
AIM1 | 0.195 |
Others | |
TIAM1 | 0.097 |
CD5 | 0.209 |
Gene class and specific genes . | Fold change . |
---|---|
TFH signature | |
CXCL13 | 0.057 |
CD200 | 0.165 |
ICOS | 0.158 |
CXCR5 | 0.246 |
Immunomodulation | |
IDO1 | 0.074 |
IL4I1 | 0.160 |
Tumor suppressor | |
AIM1 | 0.195 |
Others | |
TIAM1 | 0.097 |
CD5 | 0.209 |
A complete list of genes is given in the supplemental Table 1.
Genes overexpressed in HSTL compared with NKTCL . | |
---|---|
Gene class and specific genes . | Fold change . |
NK-cell–associated molecules | |
KLRB1 | 17.26 |
KLRG1 | 6.81 |
Microenvironment | |
Hemoglobin delta | 32.87 |
Hemoglobin alpha | 22.79 |
Hemoglobin beta | 19.46 |
Hemoglobin gamma | 15.18 |
CD36 | 14.63 |
TCR | |
TARP//TRGC2 | 11.12 |
CD3D | 7.34 |
CD3G | 6.89 |
CD247 (CD3Z) | 3.82 |
Chemokines | |
CXCL7 | 23.27 |
CXCL6 | 10.73 |
Transcription factors | |
NR4A2 | 17.77 |
TCF21 | 11.72 |
GATA6 | 6.16 |
PRDM16 | 5.99 |
PRDM1 | 2.81 |
Signal transduction | |
SPRY2 | 9.87 |
PLCB1 | 9.66 |
MAP4K3 | 7.81 |
RHOB | 7.61 |
SPRY1 | 5.59 |
MAP3K5 | 2.35 |
Wnt signaling | |
TCF7L2 | 6.35 |
TLE1 | 5.55 |
Oncogenes | |
VAV3 | 5.25 |
FYN | 4.29 |
MAF | 3.46 |
Multidrug resistance | |
ABCB1 | 2.54 |
Growth factor | |
PDGFD | 9.26 |
Cell adhesion/ECM | |
CD11d | 28.67 |
ADAMTS17 | 9.53 |
Apoptosis | |
BCL2L11 | 4.59 |
Cell cycle | |
CCNL1 | 3.34 |
CCND3 | 2.42 |
Others | |
S1PR5 | 39.64 |
TMEM178 | 32.26 |
TANC1 | 32.24 |
FCRL3 | 20.27 |
CD5L | 14.89 |
CXXC5 | 14.83 |
Genes overexpressed in HSTL compared with NKTCL . | |
---|---|
Gene class and specific genes . | Fold change . |
NK-cell–associated molecules | |
KLRB1 | 17.26 |
KLRG1 | 6.81 |
Microenvironment | |
Hemoglobin delta | 32.87 |
Hemoglobin alpha | 22.79 |
Hemoglobin beta | 19.46 |
Hemoglobin gamma | 15.18 |
CD36 | 14.63 |
TCR | |
TARP//TRGC2 | 11.12 |
CD3D | 7.34 |
CD3G | 6.89 |
CD247 (CD3Z) | 3.82 |
Chemokines | |
CXCL7 | 23.27 |
CXCL6 | 10.73 |
Transcription factors | |
NR4A2 | 17.77 |
TCF21 | 11.72 |
GATA6 | 6.16 |
PRDM16 | 5.99 |
PRDM1 | 2.81 |
Signal transduction | |
SPRY2 | 9.87 |
PLCB1 | 9.66 |
MAP4K3 | 7.81 |
RHOB | 7.61 |
SPRY1 | 5.59 |
MAP3K5 | 2.35 |
Wnt signaling | |
TCF7L2 | 6.35 |
TLE1 | 5.55 |
Oncogenes | |
VAV3 | 5.25 |
FYN | 4.29 |
MAF | 3.46 |
Multidrug resistance | |
ABCB1 | 2.54 |
Growth factor | |
PDGFD | 9.26 |
Cell adhesion/ECM | |
CD11d | 28.67 |
ADAMTS17 | 9.53 |
Apoptosis | |
BCL2L11 | 4.59 |
Cell cycle | |
CCNL1 | 3.34 |
CCND3 | 2.42 |
Others | |
S1PR5 | 39.64 |
TMEM178 | 32.26 |
TANC1 | 32.24 |
FCRL3 | 20.27 |
CD5L | 14.89 |
CXXC5 | 14.83 |
Genes down-regulated in HSTL compared with NKTCL . | |
---|---|
Cytotoxic molecules | |
Granzyme B | 0.15 |
Granzyme K | 0.21 |
Chemokines | |
CXCL10 | 0.05 |
Apoptosis | |
BIRC5 | 0.26 |
Immunomodulation | |
IL4I1 | 0.26 |
Tumor suppressor | |
AIM1 | 0.26 |
Genes down-regulated in HSTL compared with NKTCL . | |
---|---|
Cytotoxic molecules | |
Granzyme B | 0.15 |
Granzyme K | 0.21 |
Chemokines | |
CXCL10 | 0.05 |
Apoptosis | |
BIRC5 | 0.26 |
Immunomodulation | |
IL4I1 | 0.26 |
Tumor suppressor | |
AIM1 | 0.26 |
A complete list of genes is given in the supplemental Table 1.
AIM1, reported as a tumor-suppressor gene,28 was underexpressed in HSTL compared with both PTCL-NOS and NKTCL. The genes underexpressed in HSTL compared with PTCL-NOS included genes encoding TFH-associated molecules (eg, CXCL13, ICOS, CD200, and CXCR5), others involved in immunomodulation (eg, IDO1 and IL4I1), and CD5, a T-cell surface antigen known to be absent in HSTL cells.
Candidate genes relevant to HSTL pathogenesis
To search for candidate genes relevant to the pathogenesis of HSTL, we compared the HSTL signature with that of normal γδ T cells. A selection of significantly overexpressed and underexpressed genes distinguishing HSTL from normal activated or resting γδ T cells is summarized in Table 4 (P < .01). Among the most overexpressed genes in HSTL were those related to NK-cell–associated molecules, such as KIRs, KLR, CD244, and NCAM1. The other overexpressed genes were related to oncogenes (eg, FOS, VAV3, MAF, and BRAF), cell adhesion (eg, VCAM1, CD11d, and ICAM1), tyrosine kinases (eg, SYK), signal transduction (eg, SPRY2, RHOB, MAP4K3, and SPRY1), the sonic hedgehog pathway (eg, GLI3, PRKAR2B, PRKACB, and PRKAR1A), the WNT pathway (eg, FRZB, TCF7L2, BAMBI, TLE1, CTNNB1, APC, and FZD5), and S1PR5. Interestingly, many of these genes also appeared to be overexpressed in HSTL-sorted cells compared with normal activated or resting cells γδ T cells (supplemental Table 1). This was the case for genes encoding FOS, VAV3, RHOB, S1PR5, and NK-cell–associated molecules, indicating that their overexpression could be attributed to the neoplastic cells. These findings were confirmed by qRT-PCR analysis (Figure 2) demonstrating overexpression of FOS, RHOB, VAV3, S1PR5, and KIR3DS1 mRNAs in HSTL-sorted cells, and to a lesser extent in HSTL tissues, compared with normal γδ T cells. Furthermore, another oncogene, FOSB (belonging to the AP1 family), was overexpressed in HSTL cells compared with activated γδ T cells. This finding was corroborated by qRT-PCR analysis (Figure 2D) and by immunohistochemistry showing the presence of FosB-positive cells coexpressing CD2 or TIA1 in HSTL (Figure 3A and Table 5). In addition, membrane CD56 (9/17) and cytoplasmic Syk (8/9) were demonstrated in the neoplastic cells infiltrating the sinuses in the spleen or BM (Figure 3C-D and Table 5).
Genes overexpressed in HSTL compared with normal γδ T cells (P < .01) . | ||
---|---|---|
Gene class and specific genes . | Fold change . | |
Activated . | Resting . | |
Oncogenes | ||
FOS* | 6.50 | |
VAV3† | 3.52 | 3.67 |
MAF | 4.42 | |
NK-cell–associated molecules | ||
KIR2DS2*† | 13.98 | 16.63 |
KLRC3*† | 11.33 | 14.87 |
KIR3DL1//KIR3DL2*† | 9.56 | 11.16 |
KIR2DL1*† | 8.90 | 7.57 |
KIR2DL2*† | 8.60 | 8.85 |
KIR2DS5*† | 7.26 | 5.19 |
KIR3DS1*† | 7.03 | 6.34 |
KIR2DL3*† | 6.83 | 9.06 |
KIR2DL5A*† | 6.40 | 6.09 |
KIR3DL3*† | 5.75 | 6.51 |
KLRC4*† | 5.70 | 5.17 |
KIR2DS3*† | 5.69 | 6.34 |
KIR2DS1*† | 5.00 | 5.17 |
KIR2DS4*† | 3.23 | 3.67 |
CD244* | 7.23 | 3.18 |
NCAM1* | 6.27 | |
Microenvironment | ||
Hemoglobin beta | 46.31 | 34.14 |
Hemoglobin alpha | 34.93 | 24.06 |
CXCL12 | 31.57 | 27.39 |
Immunoglobulin | 26.29 | 21.08 |
CD163 | 8.86 | |
Signal transduction | ||
SPRY2*† | 13.77 | 33.98 |
RHOB*† | 7.03 | |
MAP4K3*† | 5.96 | 9.18 |
SPRY1*† | 6.50 | 7.97 |
Sonic hedgehog | ||
GLI3* | 8.86 | 5.94 |
PRKAR2B* | 8.25 | 3.15 |
PRKACB† | 3.25 | 3.75 |
PRKAR1A† | 2.80 | 7.63 |
Tyrosine kinase | ||
SYK* | 10.59 | 2.68 |
Cell adhesion | ||
CD11d*† | 37.57 | 33.11 |
VCAM1 | 15.06 | 71.58 |
ICAM1 | 5.71 | |
Wnt pathways | ||
FRZB | 14.01 | 14.19 |
TCF7L2* | 10.96 | |
BAMBI | 8.58 | 7.52 |
TLE1* | 7.16 | 4.48 |
CTNNB1* | 3.99 | 11.24 |
APC | 3.60 | 2.33 |
FZD5 | 2.53 | 3.43 |
Others | ||
S1PR5* | 8.16 | 3.51 |
Genes overexpressed in HSTL compared with normal γδ T cells (P < .01) . | ||
---|---|---|
Gene class and specific genes . | Fold change . | |
Activated . | Resting . | |
Oncogenes | ||
FOS* | 6.50 | |
VAV3† | 3.52 | 3.67 |
MAF | 4.42 | |
NK-cell–associated molecules | ||
KIR2DS2*† | 13.98 | 16.63 |
KLRC3*† | 11.33 | 14.87 |
KIR3DL1//KIR3DL2*† | 9.56 | 11.16 |
KIR2DL1*† | 8.90 | 7.57 |
KIR2DL2*† | 8.60 | 8.85 |
KIR2DS5*† | 7.26 | 5.19 |
KIR3DS1*† | 7.03 | 6.34 |
KIR2DL3*† | 6.83 | 9.06 |
KIR2DL5A*† | 6.40 | 6.09 |
KIR3DL3*† | 5.75 | 6.51 |
KLRC4*† | 5.70 | 5.17 |
KIR2DS3*† | 5.69 | 6.34 |
KIR2DS1*† | 5.00 | 5.17 |
KIR2DS4*† | 3.23 | 3.67 |
CD244* | 7.23 | 3.18 |
NCAM1* | 6.27 | |
Microenvironment | ||
Hemoglobin beta | 46.31 | 34.14 |
Hemoglobin alpha | 34.93 | 24.06 |
CXCL12 | 31.57 | 27.39 |
Immunoglobulin | 26.29 | 21.08 |
CD163 | 8.86 | |
Signal transduction | ||
SPRY2*† | 13.77 | 33.98 |
RHOB*† | 7.03 | |
MAP4K3*† | 5.96 | 9.18 |
SPRY1*† | 6.50 | 7.97 |
Sonic hedgehog | ||
GLI3* | 8.86 | 5.94 |
PRKAR2B* | 8.25 | 3.15 |
PRKACB† | 3.25 | 3.75 |
PRKAR1A† | 2.80 | 7.63 |
Tyrosine kinase | ||
SYK* | 10.59 | 2.68 |
Cell adhesion | ||
CD11d*† | 37.57 | 33.11 |
VCAM1 | 15.06 | 71.58 |
ICAM1 | 5.71 | |
Wnt pathways | ||
FRZB | 14.01 | 14.19 |
TCF7L2* | 10.96 | |
BAMBI | 8.58 | 7.52 |
TLE1* | 7.16 | 4.48 |
CTNNB1* | 3.99 | 11.24 |
APC | 3.60 | 2.33 |
FZD5 | 2.53 | 3.43 |
Others | ||
S1PR5* | 8.16 | 3.51 |
Genes down-regulated in HSTL compared with normal γδ T cells (P < .01) . | ||
---|---|---|
Cytotoxic molecules | ||
Granulysin*† | 0.03 | 0.03 |
Granzyme H*† | 0.03 | 0.05 |
Granzyme B* | 0.04 | 0.17 |
Granzyme K*† | 0.26 | 0.10 |
Chemokines | ||
LTA* | 0.02 | |
TNF | 0.03 | |
IFNG* | 0.03 | |
Tumor suppressor | ||
AIM1*† | 0.34 | 0.24 |
Others | ||
CD5*† | 0.24 | 0.20 |
Genes down-regulated in HSTL compared with normal γδ T cells (P < .01) . | ||
---|---|---|
Cytotoxic molecules | ||
Granulysin*† | 0.03 | 0.03 |
Granzyme H*† | 0.03 | 0.05 |
Granzyme B* | 0.04 | 0.17 |
Granzyme K*† | 0.26 | 0.10 |
Chemokines | ||
LTA* | 0.02 | |
TNF | 0.03 | |
IFNG* | 0.03 | |
Tumor suppressor | ||
AIM1*† | 0.34 | 0.24 |
Others | ||
CD5*† | 0.24 | 0.20 |
A complete list of genes is given in the supplemental Table 1.
Genes also differentially expressed in HSTL cells compared with normal activated γδ T cells.
Genes also differentially expressed in HSTL cells compared with normal resting γδ T cells.
. | HSTL (n = 20) . | DERL2 . |
---|---|---|
FosB | 5*/14 | |
CD163 | 0/3 | ND |
ICAM1 | 0/7 | ND |
VCAM1 | 4†/12 | ND |
Granzyme H | 0/7 | |
Beta-catenin | 0/5 | ND |
Syk | 8/9 | + |
Blimp-1 | 7/7 | + |
Bcl10 | 5/5 | ND |
GSTP1 | 4/4 | ND |
CD56 | 9/17 | + |
. | HSTL (n = 20) . | DERL2 . |
---|---|---|
FosB | 5*/14 | |
CD163 | 0/3 | ND |
ICAM1 | 0/7 | ND |
VCAM1 | 4†/12 | ND |
Granzyme H | 0/7 | |
Beta-catenin | 0/5 | ND |
Syk | 8/9 | + |
Blimp-1 | 7/7 | + |
Bcl10 | 5/5 | ND |
GSTP1 | 4/4 | ND |
CD56 | 9/17 | + |
ND indicates not determined.
Double immunostaining showing presence of FosB-positive lymphocytes coexpressing CD2 or TIA1 in 3 investigated cases.
With weak (equivocal) staining in tumor cells.
Genes associated with cytotoxicity (eg, Granulysin, Granzyme H, Granzyme K, and Granzyme B), cytokines (eg, LTA, TNF, and IFNG), the tumor suppressor AIM1, and CD5 were among those genes significantly underexpressed in HSTL compared with their normal cell counterparts. The absence of expression of cytotoxic molecules such as Granzyme B (Table 1) and Granzyme H was shown by immunostaining (Figure 3H and Table 5), whereas the dramatic reduction of AIM1 transcripts was demonstrated by qRT-PCR in both HSTL tissues and sorted cells compared with normal γδ cells.
Biomarkers discriminating HSTL from other lymphomas and normal controls
Expression of several genes in the molecular signature of HSTL appeared most likely to be attributed to the normal spleen signature. This was the case for the ICAM-1, VCAM-1, and CD163 gene products, the expression of which was shown by immunohistochemistry to be restricted to nonneoplastic cells; that is, endothelial cells (VCAM1 and ICAM1) and macrophages (CD163). To build a biomarker list differentiating HSTL from other T-cell lymphoma entities and eliminating the observed microenvironmental signature, supervised tests (see supplemental Methods for details) were performed between HSTL and other lymphomas (ie, AITL, PTCL-NOS, and NKTCL), between HSTL and control samples (ie, normal lymph nodes, spleen, γδ T cells, and B cells), and between HSTL obtained from tissues and cell suspensions. This yielded a list of 54 genes distinguishing HSTL from other lymphoma entities and normal controls (supplemental Figure 2 and supplemental Tables 2-4). Among these 54 genes were those encoding KIR molecules (KIR3DS1), CD244, DTNBP1, and S1PR5. The application of this list of biomarkers to an independent series of lymphoma samples (including 4 HSTL, 15 AITL, and 11 PTCL-NOS samples) using a different platform17 allowed the discrimination of HSTL from other PTCLs (supplemental Figure 3).
Relevant pathways and genes potentially involved in HSTL resistance to therapy
Pathway analyses showed significant enrichment of several pathways that are involved in cancer development or that can be the target of innovative treatments, such as VEGF, MAPK, JAK-STAT, mTOR, Notch signaling pathways, cytokine-cytokine receptor interaction, and cell adhesion (Table 6 and supplemental Table 5). The Wnt signature was found to be deregulated, but activation of the pathway was not supported by immunohistochemistry because the expression of β-catenin was restricted to endothelial cells. Consistent with previous findings,17 pathways related to NK-cell–mediated cytotoxicity were among the most differentially expressed compared with other T-cell lymphoma entities and normal γδ T cells (Figure 4A). Furthermore, as illustrated in Figure 4B, many genes in the TCR signaling pathway (eg, VAV3 and BCL10) that were linked to the AP1 pathway (eg, FOS) were enriched in HSTL compared with PTCL-NOS and NKTCL. These finding were supported by the high levels of VAV3, FOS, and FOSB transcripts demonstrated by qRT-PCR and by immunohistochemical detection of Bcl-10 in the neoplastic cells of all investigated HSTL samples (5 of 5; Figure 3D and Table 5). Significant enrichment of genes in the pathways of Sprouty regulation of tyrosine kinase and the Sonic Hedgehog signaling, such as genes for the composition of cAMP-dependent kinase (eg, PRKAR2B, PRKAR1A, and PRKACB), the transcription factor Ci (eg, GLI3), the receptor patched (PTCH1), and GSK3B, were found in HSTL (Figure 4). Interestingly, several genes in the multidrug resistance signaling pathways (eg, ABCB1 and GSTP1) were enriched in HSTL (Figure 4E) compared with PTCL-NOS and NKTCL and PRDM1, another gene reported to be associated with resistance to therapy in PTCLs.29 High levels of ABCB1 transcripts were confirmed by qRT-PCR (Figure 2B), and expression of GSTP1 and Blimp-1 was demonstrated by immunohistochemistry in the HSTL cells of all evaluated cases, as well as in DERL2 cells (Figure 3). As illustrated in Figure 4F, HSTL cells were closer to activated γδ T cells than resting γδ T cells with respect to the expression of cell-cycle–related genes.
Pathway name . | GlobalTest P value . | |||
---|---|---|---|---|
HSTL vs PTCL-NOS . | HSTL vs NKTCL . | HSTL vs activated γδ T cell . | HSTL vs resting γδ T cell . | |
Biocarta 133: Ras-independent in NK-cell–mediated cytotoxicity | 3 × 10-6 | 3 × 10-2 | 3 × 10-4 | 3 × 10-3 |
KEGG 04650: NK-cell–mediated cytotoxicity | 5 × 10-6 | 7 × 10-3 | 7 × 10-5 | 3 × 10-3 |
KEGG 04370: VEGF signaling pathway | 6 × 10-6 | 1 × 10-3 | 2 × 10-4 | 3 × 10-3 |
KEGG 04060: cytokine-cytokine receptor interaction | 6 × 10-6 | 4 × 10-3 | 1 × 10-4 | 7 × 10-3 |
Biocarta 131: phospholipase C signaling pathway | 7 × 10-6 | 3 × 10-4 | 3 × 10-3 | 8 × 10-2 |
Biocarta 140: selective expression of chemokine receptors during T-cell polarization | 7 × 10-6 | 3 × 10-2 | 5 × 10-5 | 1 × 10-2 |
Biocarta 506: sprouty regulation of tyrosine kinase signals | 8 × 10-6 | 2 × 10-4 | 2 × 10-4 | 3 × 10-3 |
KEGG 04310: Wnt signaling pathway | 1 × 10-5 | 7 × 10-4 | 2 × 10-4 | 3 × 10-3 |
Biocarta 590: Sonic Hedgehog (Shh) pathway | 1 × 10-5 | 5 × 10-3 | 6 × 10-5 | 3 × 10-3 |
Biocarta 137: regulation of BAD phosphorylation | 1 × 10-5 | 7 × 10-3 | 9 × 10-5 | 3 × 10-3 |
KEGG 04010: MAPK signaling pathway | 1 × 10-5 | 3 × 10-3 | 2 × 10-4 | 3 × 10-3 |
KEGG 04340: Hedgehog signaling pathway | 3 × 10-5 | 1 × 10-3 | 5 × 10-4 | 2 × 10-2 |
Biocarta 586: multidrug resistance factors | 1 × 10-5 | 5 × 10-3 | 2 × 10-2 | NS |
Biocarta 465: role of MEF2D in T-cell apoptosis | 3 × 10-5 | 1 × 10-3 | 6 × 10-4 | 3 × 10-3 |
KEGG 04660: TCR signaling pathway | 3 × 10-5 | 3 × 10-3 | 8 × 10-5 | NS |
KEGG 04630: Jak-STAT signaling pathway | 4 × 10-5 | 9 × 10-3 | 1 × 10-4 | 3 × 10-3 |
KEGG 04620: TLR signaling pathway | 5 × 10-5 | 7 × 10-3 | 5 × 10-4 | 7 × 10-3 |
KEGG 04514: cell adhesion molecules | 6 × 10-5 | 3 × 10-3 | 3 × 10-4 | 7 × 10-3 |
KEGG 04210: apoptosis | 7 × 10-5 | 1 × 10-2 | 1 × 10-4 | 3 × 10-3 |
KEGG 04150: mTOR signaling pathway | 1 × 10-4 | 3 × 10-3 | 2 × 10-4 | 3 × 10-3 |
KEGG 04330: Notch signaling pathway | 2 × 10-4 | 1 × 10-3 | 3 × 10-4 | 2 × 10-2 |
KEGG 04670: leukocyte transendothelial migration | 3 × 10-4 | 1 × 10-2 | 1 × 10-3 | 3 × 10-3 |
KEGG 05200: pathways in cancer | 3 × 10-4 | 2 × 10-2 | 5 × 10-4 | 3 × 10-3 |
KEGG 04110: cell cycle | 1 × 10-2 | 5 × 10-3 | 3 × 10-3 | 3 × 10-3 |
KEGG 03050: proteasome | 5 × 10-3 | 3 × 10-3 | 2 × 10-4 | 1 × 10-2 |
Pathway name . | GlobalTest P value . | |||
---|---|---|---|---|
HSTL vs PTCL-NOS . | HSTL vs NKTCL . | HSTL vs activated γδ T cell . | HSTL vs resting γδ T cell . | |
Biocarta 133: Ras-independent in NK-cell–mediated cytotoxicity | 3 × 10-6 | 3 × 10-2 | 3 × 10-4 | 3 × 10-3 |
KEGG 04650: NK-cell–mediated cytotoxicity | 5 × 10-6 | 7 × 10-3 | 7 × 10-5 | 3 × 10-3 |
KEGG 04370: VEGF signaling pathway | 6 × 10-6 | 1 × 10-3 | 2 × 10-4 | 3 × 10-3 |
KEGG 04060: cytokine-cytokine receptor interaction | 6 × 10-6 | 4 × 10-3 | 1 × 10-4 | 7 × 10-3 |
Biocarta 131: phospholipase C signaling pathway | 7 × 10-6 | 3 × 10-4 | 3 × 10-3 | 8 × 10-2 |
Biocarta 140: selective expression of chemokine receptors during T-cell polarization | 7 × 10-6 | 3 × 10-2 | 5 × 10-5 | 1 × 10-2 |
Biocarta 506: sprouty regulation of tyrosine kinase signals | 8 × 10-6 | 2 × 10-4 | 2 × 10-4 | 3 × 10-3 |
KEGG 04310: Wnt signaling pathway | 1 × 10-5 | 7 × 10-4 | 2 × 10-4 | 3 × 10-3 |
Biocarta 590: Sonic Hedgehog (Shh) pathway | 1 × 10-5 | 5 × 10-3 | 6 × 10-5 | 3 × 10-3 |
Biocarta 137: regulation of BAD phosphorylation | 1 × 10-5 | 7 × 10-3 | 9 × 10-5 | 3 × 10-3 |
KEGG 04010: MAPK signaling pathway | 1 × 10-5 | 3 × 10-3 | 2 × 10-4 | 3 × 10-3 |
KEGG 04340: Hedgehog signaling pathway | 3 × 10-5 | 1 × 10-3 | 5 × 10-4 | 2 × 10-2 |
Biocarta 586: multidrug resistance factors | 1 × 10-5 | 5 × 10-3 | 2 × 10-2 | NS |
Biocarta 465: role of MEF2D in T-cell apoptosis | 3 × 10-5 | 1 × 10-3 | 6 × 10-4 | 3 × 10-3 |
KEGG 04660: TCR signaling pathway | 3 × 10-5 | 3 × 10-3 | 8 × 10-5 | NS |
KEGG 04630: Jak-STAT signaling pathway | 4 × 10-5 | 9 × 10-3 | 1 × 10-4 | 3 × 10-3 |
KEGG 04620: TLR signaling pathway | 5 × 10-5 | 7 × 10-3 | 5 × 10-4 | 7 × 10-3 |
KEGG 04514: cell adhesion molecules | 6 × 10-5 | 3 × 10-3 | 3 × 10-4 | 7 × 10-3 |
KEGG 04210: apoptosis | 7 × 10-5 | 1 × 10-2 | 1 × 10-4 | 3 × 10-3 |
KEGG 04150: mTOR signaling pathway | 1 × 10-4 | 3 × 10-3 | 2 × 10-4 | 3 × 10-3 |
KEGG 04330: Notch signaling pathway | 2 × 10-4 | 1 × 10-3 | 3 × 10-4 | 2 × 10-2 |
KEGG 04670: leukocyte transendothelial migration | 3 × 10-4 | 1 × 10-2 | 1 × 10-3 | 3 × 10-3 |
KEGG 05200: pathways in cancer | 3 × 10-4 | 2 × 10-2 | 5 × 10-4 | 3 × 10-3 |
KEGG 04110: cell cycle | 1 × 10-2 | 5 × 10-3 | 3 × 10-3 | 3 × 10-3 |
KEGG 03050: proteasome | 5 × 10-3 | 3 × 10-3 | 2 × 10-4 | 1 × 10-2 |
NS indicates not significant.
Isochromosome 7q is the main recurrent genomic aberration in HSTL
The aCGH findings confirmed recurrent isochromosome 7q in 4 of the 7 tested cases (60%), with only few other alterations comprising trisomy 8, as described previously,9 loss of the 10p arm and more focal alterations (−14q,+15q, −18q −21q, and −22q) in 2-6 cases (30%-85%; supplemental Figure 4). To determine whether iso7q affects the molecular signature, we analyzed the mean expression of genes located on the 7p and 7q arms for each sample (supplemental Figure 5). This analysis clearly showed overexpression of 7q and down-expression of 7p genes in all cases with the iso7q aberration as well as in 2 additional cases (HSTL_08 and HSTL_09) for which the iso7q status was unknown (due to the absence of aCGHa and FISH data), one of which was shown subsequently to be iso7q by FISH (HSTL_09). To identify potential candidate genes involved in iso7q, we also compared the expression level of iso7q+ HSTL cases with normal controls (ie, activated γδ T cells). Among the most repressed genes on the 7p arm were CYCS and IKZF1, which are associated with regulation of apoptosis, and HUS1 and CBX3, which are involved in DNA repair, whereas the most overexpressed genes on the 7q arm included the putative oncogene PTPN12 (supplemental Table 1). Moreover, pathway analysis comparing iso7q+ HSTL with activated or resting γδ T cells showed a significant enrichment of apoptosis-related pathways (supplemental Table 5).
AIM1 is methylated in HSTL
AIM1, initially described as a tumor-suppressor gene deleted in melanoma, was subsequently found to be deregulated in NKTCL.28,30 In the present study, we show significant down-expression of AIM1 mRNA in HSTL and DERL2 cells. Low levels of AIM1 mRNAs were confirmed by qRT-PCR in HTSL primary tumors and in DERL2 cells compared with normal γδ T cells (Figure 4A). We also searched for genetic alterations underlying AIM1 down-regulation. In the absence of genomic imbalance in the AIM1 locus at 6q21 in the aCGH data, we wondered whether AIM1 down-regulation in HSTL could be the result of methylation, as was reported previously in other hematologic malignancies such as myeloma and NKTCL.30,31 The promoter methylation status of AIM1 was analyzed in the HSTL-derived cell line DERL-2. Sequencing of AIM1 on bisulfite-treated DNA revealed methylation in both CpG islands (Figure 5B). Furthermore, treatment of DERL2 cells containing the AIM1-methylated promoter with the DNA demethylating agent decitabine with or without TSA resulted in a 2- to 3-fold increase in the AIM1 mRNA level compared with untreated cells, supporting the idea that AIM1 down-expression in HSTL could be reversed by demethylation (Figure 5C). Interestingly, the number of DERL2 apoptotic cells was increased after decitabine treatment (Figure 5D).
Syk as a candidate target for pharmacologic inhibition in HSTL
Whereas normal mature T lymphocytes lack Syk expression, the Syk protein has been recently reported to be aberrantly expressed in many PTCLs.32 In the present study, GEP analysis revealed overexpression of SYK mRNA in HSTL tissues, in HSTL-sorted cells, and in the DERL2 cell line compared with normal γδ T cells (Figure 6A). We also demonstrated Syk expression by neoplastic HSTL cells and DERL2 cells using immunohistochemistry and Western blot (Figure 3B, Table 5, and Figure 6B-C). Because the expression of Syk does not provide information on its activity, we analyzed the phosphorylation status of Syk in DERL2 cells using a phosphospecific anti-Syk (Tyr525/526) Ab. Syk was constitutively activated in DERL2 cells (Figure 6C). However, the mechanism leading to Syk activation is unclear because neither SYK amplification nor ITK-SYK fusion was detected by FISH analysis in DERL2 cells (data not shown). To determine whether Syk could contribute to the growth and survival of neoplastic HSTL cells, a chemical anti-Syk agent was used to inhibit Syk activity in DERL2 cells. Syk inhibition induced an significant increase in apoptosis in a dose-dependent manner in the DERL2 cell line compared with normal activated γδ T cells (Figure 6D).
Discussion
HSTL is a rare and distinct PTCL entity of unknown pathogenesis. The disease is characterized by a peculiar distribution of the neoplastic cells, most often of γδ T-cell derivation, within the sinuses of the spleen, liver, and BM and by a clinically fatal outcome. Previous scarce reports showed that HSTL has a distinct gene signature.13,17 Taking advantage of the exceptional availability of samples enriched in tumor cells, we showed herein by consensus clustering that HSTL expression profiles clustered in a single branch independently of the spleen, BM, or blood origin of the diagnostic samples, indicating that this signature is independent of the microenvironment. The distinct signature was further confirmed when a gene list was constructed by eliminating the signature of the normal spleen and lymph nodes. HSTL was initially reported as being of γδ T-cell derivation, but similar cases with a αβ phenotype have since been reported.
In the present study, we show that 2 cases of HSTL with a αβ phenotype had a gene signature and a genomic profile that were very similar to the classic γδ HSTL variant, providing strong molecular arguments for grouping γδ and αβ HSTL as a single entity, as was proposed in the current WHO classification.8 The similar αβ and γδ HSTL molecular signature was largely contributed by genes implicated in cytotoxic T- and NK-cell function and development, suggesting that both αβ and γδ variants are derived from subsets of cells of the innate immune system.8,33,34 Indeed, a subset of αβ T cells, referred as NKT cells, share characteristics with NK cells and γδ T cells (including expression of CD161 (KLRB1), CD16, CD56, and KIR molecules) and participate with NK cells and γδ T cells in the innate immune system.35 γδ T cells can accumulate in immunosuppressed patients and in the setting of chronic antigen stimulation. Expansion of γδ T cells has been reported particularly in patients after kidney transplantation, systemic lupus, Hodgkin lymphoma, or malaria4 and, more recently, also after treatment with anti-TNF agents.36
The elective distribution of neoplastic cells within the sinuses and sinusoids of the BM, spleen, and liver is remarkable and accounts for the peculiar clinical presentation of HSTL with hepatosplenomegaly and cytopenia but without lymphadenopathy and leukemia. Interestingly, one of the genes with the highest change of expression compared with PTCL-NOS and NKTCL and, to a lesser extent, normal γδ T cells, was sphingosine-1-phosphatase receptor 5 (S1PR5). S1PR5 encodes a member of the family of S1P receptors known to be involved in T- and B-cell exit from lymphoid organs.37 S1PR5 is preferentially expressed by mature CD56dim human NK cells38 and is involved in NK-cell trafficking in steady-state and inflammatory situations. It has been reported that S1PR5-deficient mice have defective homing of NK cells to the blood and spleen.27,38 Although the role of S1PR5 in γδ T cells is uncharacterized, the high level of S1PR5 mRNA in HSTL might explain the peculiar distribution and accumulation of neoplastic γδ cells in the spleen and BM, preventing tumor cells from exiting sinusoids. Whether cell-cell interactions involving the cell adhesion molecule LFA1 (CD11a), which is expressed by HSTL cells, and VCAM1, which is expressed by splenic sinusoidal cells, could also play a role in the clinicopathologic picture of the disease remains to be determined.
Although our GEP findings were based on the comparison with normal blood γδ T cells and not with splenic γδ T cells, the down- or overexpression of selected genes compared with normal splenic γδ T cells was validated by qRT-PCR. From the results of the present study, it appears that several candidate genes involved in distinct pathways (ie, TCR signaling, cell cycle, sonic hedgehog, sprouty regulation of tyrosine kinase pathways, and multidrug resistance) could potentially be associated with the pathophysiology of HSTL. One of the genes we found to be overexpressed in HSTL is MDR-1, the expression of which was significant for this lymphoma compared with normal γδ T cells, PTCL-NOS, and NKTCL. The P-glycoprotein (P-gp), encoded by MDR-1, is able to export drugs out of the cytoplasm, thus reducing the intracellular accumulation of cytotoxic drugs. P-gp has been reported to be overexpressed and associated with poor outcome in various hematologic malignancies including NKTCL and diffuse large B-cell lymphoma (DLBCL).39-41 Interestingly, sensitivity to aracytine is not influenced by P-gp overexpression in multidrug resistant cell lines.42 Although HSTL is reported as a fatal disease, recent studies have described long-time survivors after cytarabine-based chemotherapy.4 The present findings might provide a rationale for the use of cytarabine-containing regimens in HSTL. In addition to P-gp expression, we have also demonstrated that GST, another molecule involved in resistance to alkylating agents and anthracyclines,43 was also expressed by the neoplastic HSTL cells (as shown by immunohistochemistry). One could hypothesize that the combination of P-gp and GST-π expression might contribute to the resistance of HSTL to anthracycline-based chemotherapy regimens and thus to its poor outcome. High levels of GST-π are present in the neoplastic cells of mantle cell lymphoma and DLBCL44,45 and are correlated with poor prognosis in DLBCL.44 The role of the transcription factor Blimp-1 in PTCL is unclear. In the present study, we show that PRDM1 is expressed at the mRNA and protein levels in HSTL cells. Whereas inactivation of PRDM1 through epigenetic mechanisms has been reported recently in NKTCL, suggesting a role as tumor-suppressor gene,30,46 its expression has been associated with chemoresistance in T-cell lymphoma.29
Tumor development is a multistep process involving deregulation of oncogenes and tumor suppressor genes (TSGs). To date, no oncogene or TSG has been identified in HSTL. AIM1 is a gene with expression that has been shown to be altered in association with tumor suppression in a model of human melanoma.28 It is frequently deleted or methylated in several solid tumors and in myeloma and, more interestingly, in NK cell lines.14,30,31 These previous data and our present findings that AIM1 is reduced dramatically in HSTL, most likely because of promoter methylation, suggest that AIM1 might play a role as a TSG in HSTL oncogenesis and provide the rationale for testing demethylating agents in this disease. Our expression data also suggest that the recurrent iso7q aberration in HSTL could target putative TSG and oncogenes, including IKZF1 and PTPN12.
With respect to potential targets for novel therapies in HSTL, we also demonstrated Syk expression by neoplastic HSTL cells with constitutive activation in DERL2 cells compared with normal T cells. Syk is a protein tyrosine kinase involved in BCR signaling, and its activation is associated with cell growth and survival in B-cell lymphomas.47,48 Normal mature T lymphocytes lack Syk expression.49 However, the Syk protein has been reported recently as a feature common to most PTCLs; this may represent a novel therapeutic target, because, in addition to its aberrant expression in many PTCLs, inhibition of Syk induces apoptosis and blocks proliferation in T-cell lymphoma cell lines.32,49 An orally available Syk inhibitor is currently being tested for B-cell lymphomas in a phase 1/2 clinical trial, with encouraging results.50 In view of the present findings that Syk inhibition induced apoptosis in DERL2 cells in a dose-dependent manner, Syk could also be a potential candidate target for pharmacologic inhibition in HSTL.
In conclusion, the results of the present study show that the distinctive molecular signature of HSTL partly reflects the ontogeny and might explain several clinicopathologic features of the disease, including homing to the sinusoids of the spleen and BM. Our study highlights emerging oncogenic pathways that may be implicated in the pathophysiology of the disease and may offer a rationale for exploring new therapeutic options such as tyrosine kinase inhibitors and demethylating agents in these patients. Finally, this study also provides additional molecular arguments for grouping αβ and γδ HSTL as a single distinct entity.
The online version of this article contains a data supplement.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
The authors thank Christelle Thibault from the Institut de Génétique et de Biologie Moléculaire et Cellulaire platform (Affymetrix: Philippe Kastner); Joe A. Trapani for providing granzyme H Ab; L. Del Vecchio for providing the DERL2 cell line; Jean Jacques Fournié for providing BrHPP; Francois Radvanyi for critical opinions on data interpretation; Jérôme Cuvillier and Nadine Vailhen from Groupe d'Etude des Lymphomes de l'Adulte-Pathologie; and Virginie Fataccioli and Maryse Baia for technical assistance.
This work is part of the Carte d'Identité des Tumeurs program (http://cit.ligue-cancer.net/index.php/en) from the Ligue Nationale Contre le Cancer and of the Tenomic project supported by a Programme Hospitalier de Recherche Clinique and the Institut National du Cancer. This work was supported by Inserm, the Institut National du Cancer, the Plan Cancer of the Belgian government, and the Association pour la Recherche Thérapeutique, Génétique et Immunologique dans les Lymphomes. M.T. was supported by the Fondation pour la Recherche Médicale (DEQ 2010/0318253).
Authorship
Contribution: M.T., Y.H., N.M.-G., J. Bosq, T.M., P.G., and A.d.R. performed the research; M.T., A.d.R., and P.G. analyzed and interpreted the data; M.T., L.d.L., and P.G. designed the research; L.d.L., J. Bosq, J. Brière, M.-H.D.-L., F.B., J.S., E.M., and P.G. collected the data and/or tumor samples; A.d.R. and P.G. and supervised the bioinformatics analyses; T.M. contributed vital analytical tools; and M.T., L.d.L., A.d.R., and P.G. wrote the manuscript.
Conflict-of-interest disclosure: The authors declare no competing financial interests.
Correspondence: Philippe Gaulard, MD, Département de Pathologie and Inserm U955, Hôpital Henri Mondor, 94010 Créteil, France; e-mail: [email protected].
References
Author notes
A.d.R. and P.G. contributed equally to this work.
This feature is available to Subscribers Only
Sign In or Create an Account Close Modal