- Research article
- Open Access
- Open Peer Review
Identification of ADGRE5 as discriminating MYC target between Burkitt lymphoma and diffuse large B-cell lymphoma
BMC Cancer volume 19, Article number: 322 (2019)
MYC is a heterogeneously expressed transcription factor that plays a multifunctional role in many biological processes such as cell proliferation and differentiation. It is also associated with many types of cancer including the malignant lymphomas. There are two types of aggressive B-cell lymphoma, namely Burkitt lymphoma (BL) and a subgroup of diffuse large cell lymphoma (DLBCL), which both carry MYC translocations and overexpress MYC but both differ significantly in their clinical outcome. In DLBCL, MYC translocations are associated with an aggressive behavior and poor outcome, whereas MYC-positive BL show a superior outcome.
To shed light on this phenomenon, we investigated the different modes of actions of MYC in aggressive B-cell lymphoma cell lines subdivided into three groups: (i) MYC-positive BL, (ii) DLBCL with MYC translocation (DLBCLpos) and (iii) DLBCL without MYC translocation (DLBCLneg) for control. In order to identify genome-wide MYC-DNA binding sites a chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq) was performed. In addition, ChIP-Seq for H3K4me3 was used for determination of genomic regions accessible for transcriptional activity. These data were supplemented with gene expression data derived from RNA-Seq.
Bioinformatics integration of all data sets revealed different MYC-binding patterns and transcriptional profiles in MYC-positive BL and DLBCL cell lines indicating different functional roles of MYC for gene regulation in aggressive B-cell lymphomas. Based on this multi-omics analysis we identified ADGRE5 (alias CD97) - a member of the EGF-TM7 subfamily of adhesion G protein-coupled receptors - as a MYC target gene, which is specifically expressed in BL but not in DLBCL regardless of MYC translocation.
Our study describes a diverse genome-wide MYC-DNA binding pattern in BL and DLBCL cell lines with and without MYC translocations. Furthermore, we identified ADREG5 as a MYC target gene able to discriminate between BL and DLBCL irrespectively of the presence of MYC breaks in DLBCL. Since ADGRE5 plays an important role in tumor cell formation, metastasis and invasion, it might also be instrumental to better understand the different pathobiology of BL and DLBCL and help to explain discrepant clinical characteristics of BL and DLBCL.
The transcription factor MYC plays a multifunctional role in many cellular processes such as cell cycle progression, apoptosis and cellular transformation. Over-expression of MYC leads to an increased replication activity and is associated with different types of cancer. This holds also true for tumors of the immune system especially aggressive B-cell non-Hodgkin lymphomas (B-NHL) such as Burkitt lymphoma (BL) and diffuse large B-cell lymphoma (DLBCL). BL is an extremely fast growing tumor that carries immunoglobulin/MYC translocations in almost all cases. The tumor is predominantly found in male children but may also occur in adults especially with a compromised immune system. Treatment of BL is mainly based on high dose chemotherapy with usually favorable clinical outcome . In contrast, DLBCL rarely carries MYC rearrangements, which may be associated with both immunoglobulin and non-immunoglobulin genes. Whereas DLBCL without MYC translocation reveals long-term survival of 60–70% of the patients treated with combined immune-chemotherapy, DLBCL with MYC translocation – regardless of its translocation partner – shows a very poor clinical outcome [2,3,4,5,6,7,8]. It is currently unclear why BL and DLBCL with MYC translocations display this very different clinical course. In addition, molecular features for a precise stratification of patients into BL and DLBCL with MYC translocation are lacking despite the need for different treatment modalities. To determine the potentially different role of MYC in BL and DLBCL, we aimed at identifying their molecular features by means of chromatin immunoprecipitation combined with high-throughput sequencing (ChIP-Seq) and whole transcriptome shotgun sequencing (RNA-Seq) employing B-cell lymphoma cell lines. Validation of the results was performed with primary lymphoma tissue samples.
Three MYC break positive BL cell lines (Blue-1 / ACC-594; BL-2 / ACC-625 and BL-41 / ACC-160), two MYC break positive (Carnaval / ACC-724; U2932-R2 / ACC-633) and two MYC break negative (Karpas-422 ACC-32, U2932-R1 / ACC-633) DLBCL cell lines (overview Fig. 1a) were obtained in 2012 from the German Collection of Microorganisms and Cell Cultures (DSMZ). The sub-clones U2932-R1 and U2932-R2 were kindly provided by Dr. Quentmeier (DSMZ, Braunschweig, Germany) . All cell lines were negatively tested for mycoplasma contamination prior to use and are currently not listed as cross-contaminated or misidentified cell lines according the International Cell Line Authentication Committee (ICLAC). All cell lines were cultivated in RPMI 1640 medium supplemented with GlutaMAX™-I (Gibco, Thermo Fisher Scientific) and containing 20% of heat inactivated fetal bovine serum (PAN Biotech, Aidenbach, Germany) under a humidified atmosphere with 5% CO2 at 37 °C. Cells were thawed and continuously split 3 times per week for a maximum period of three weeks. Cell counting was performed on a BD Accuri C6 Flow Cytometer (BD Biosciences, New Jersey, United States) and cell viability was determined by propidium iodide (PI) – staining (BD Bioscience, Heidelberg, Germany) according to the manufacturer’s recommendations. Only cells, which exhibited more than 90% vitality, were used for further investigation.
1.5 × 106 vital cells were washed three times with PBS and lysed with protease inhibitors containing RIPA buffer supported by sonication. After measuring protein concentration using the BCA protein assay kit (Pierce, Thermo Fisher Scientific), protein lysates were separated under denaturing conditions via gels electrophoresis using 16% sodium dodecyl sulfate polyacrylamide gels (Invitrogen, California, United States) and transferred to Hybond-ECL nitrocellulose membranes (Amersham Biosciences, New Jersey, United States) by electroblotting. Membranes were blocked with a PBST 5% dry milk solution for 1 h followed by incubation with the respective primary antibody solution at 4 °C overnight. Subsequently, membranes were washed three times with PBST and incubated for 1 h with a secondary antibody conjugated with horseradish peroxidase (information on primary and secondary antibodies is available in Additional file 1: Table S1). Chemiluminescence was detected using HRP substrate (Luminata Forte, Merck Chemicals GmbH, Darmstadt, Germany) and FusionCapt Advance analysis Software (Fusion device, Vilber Lurmat GmBH, Eberhardzell, Germany).
Quantitative real-time PCR analysis
Total RNA was isolated from 1 × 106 vital cells after washing with PBS employing NucleoSpin RNA Kit (MACHEREY-NAGEL GmbH & Co. KG, Düren, Germany). RNA fluorometric quantification was performed by means of the Qubit RNA quantification assay (Thermo Fisher Scientific). Total RNA was reverse transcribed into complementary DNA (cDNA) using TaqMan reverse transcription reagents. Real-time PCR analysis was realized using TaqMan Real-Time PCR Master Mix on a Step One Plus Real-Time PCR System (Thermo Fisher Scientific). All procedures were performed according to the manufacturer’s recommendations. RT-PCR Taq-Man probes are listed in Additional file 1: Table S1. Relative RNA expression was calculated according to the comparative Ct method  using the average expression based on triplicates of two biological replicates of each cell line. For endogenous control b2-microglobulin (B2M) or succinate dehydrogenase complex, subunit A (SDHA) were used.
Chromatin immunoprecipitation (ChIP) was done according to published protocols [11, 12] with few modifications. Briefly, 2 × 107 vital cells were fixed for 10 min at 4 °C in medium containing 1% formaldehyde. After blocking with 0.1 M glycine and washing four times with PBS, the cells were snap frozen and stored at − 80 °C. After thawing on ice each cell pellet was resuspended in 5 mL cold LB1 lysis buffer, incubated for 10 min at 4 °C and for further 10 min ambient temperature in 5 mL LB2 lysis buffer before being finally dissolved in 3 mL LB3 buffer. Sonication was performed for 45 min [three cycles of 15 min each at high power in pulsed mode (30 s on and 30 s off)] using titanium rods combined with a Bioruptor Sonicator (Diagenode, Seraing, Belgium). After addition of 300 μL 10% (vol/vol) Triton X-100 and centrifugation the supernatant was removed, 50 μL of which were stored as input DNA sample. 1.5 mL of the supernatant was incubated with 10 μg MYC antibody or 5 μg H3K4me3 antibody at 4 °C overnight. For ChIP antibody information, refer to Additional file 1: Table S1. For precipitation of DNA indirectly bound to the respective antibody, 30 μg Dynabeads coupled with Protein G (Thermo Fisher Scientific) were added for each μg antibody and incubated for 3 h at 4 °C. Subsequently, the beads were washed and the immunoprecipitated (IP) DNA was eluted. Finally, the eluate (input DNA and IP DNA) was reverse cross-linked overnight at 65 °C followed by digestion with RNase A and Proteinase K. The resulting DNA was phenol/chloroform extracted, precipitated and the DNA was resuspended in 30 μL 10 mM Tris·HCl, pH 8.0. DNA was subjected to fluorometric quantification by the Qubit DNA quantification assay (Thermo Fisher Scientific). Ten ng of chromatin-immunoprecipitated DNA sample were processed with NEBNext ChIP-Seq Library Prep Master Mix Set for Illumina library generation according to the manufacturer’s recommendations. All amplified libraries were analyzed with the DNA 1000 Kit on the 2100 Bioanalyzer (Agilent, California, United States). Single-read NGS was done on an Illumina HiSeq 1500 system (50 cycles). Illumina adapters were trimmed from the raw sequence data and low quality bases and reads were removed with trimmomatic (LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36) . Sequence data was aligned to the main chromosomes of the human reference genome (GRCh38) with bowtie version 0.12.7 (−e 70 -k 1 -m 1 -n 2 –best) . H3K4 and MYC peaks were called with MACS2  with a q-value cut-off of 0.1 and the peaks from the two replicate ChIP samples were summarized with IDR , keeping all peaks with an IDR < 0.1. Final peaks were annotated to the nearest transcription start site (TSS) using gene annotation from Ensembl release 77. Only peaks with a maximum distance of 2000 bp to a TSS were kept. Artificial peaks were removed using the ENCODE blacklist (https://sites.google.com/site/anshulkundaje/projects/blacklists). Differential peaks between DLBCL with and without MYC break and BL were estimated using DiffBind  tool.
Total RNA was isolated from 1 × 106 lymphoma cells, which were previously spiked in with 1 × 105 insect cells (Schneider cells) for data calibration . The quality of the RNA was determined with an Agilent 4200 TapeStation and Software A.01.05 (Agilent, California, United States). 500 ng RNA per sample were processed using the Illumina TruSeq Stranded mRNA LT Sample Prep Kit following the manufacturer’s instructions to generate libraries for RNA sequencing. Samples were sequenced on a Hi-Seq 4000 (single read mode; length 150 bp) using the Illumina HiSeq 3000/4000 SBS 150 cycle kit. Sequence reads were aligned to a concatenated genome that consisted of the human (GRCh38) and the Drosophila melanogaster (BDGP5) reference genome, using STAR alignment tool  with default parameters. Gene annotation from Ensembl release 77 and feature Counts  with default parameters were used to assign read counts to human and Drosophila genes. Before differential gene expression analysis, we calculated DESeq2 sample sizeFactors  on the Drosophila gene counts and applied them to the human sample data. This way, gene expression levels of the cell lines were calibrated to the number of sample cells. Then, gene expression levels were modeled with a generalized linear model assuming negative binomial distributed data and categorical variables for the lymphoma type (BL or DLBCL) and MYC status (MYC break positive or negative). Gene expression changes were tested for significance with the Wald test and fold changes with an associated False Discovery Rate (FDR) below 0.05 were considered significant differentially expressed.
The SWATH-MS-based quantification of the proteins ADGRE5, BYSL and NPM1 was obtained from previously published data . SWATH-MS measurements were carried out on a TripleTOF 5600+ (Sciex, Darmstadt, Germany) coupled to an Ultimate 3000 nano-HPLC-system (Dionex, Idstein, Germany) using an 88 min-binary gradient. The PeakView 2.1 software (Sciex, Darmstadt, Germany) was employed for quantification of the peptides based on an in-house library. Only peptides with FDR < 1% and confidence > 95% were considered for quantification. Peptide intensities were summed up and normalized to total protein intensity. Statistical tests were conducted using heteroskedastic 1-way ANOVA.
Immunohistochemical staining was performed using sections derived from formalin-fixed paraffin-embedded cell line blocks (n = 12) and primary tissue samples (n = 38). The use of human primary tissue samples was approved by the Institutional Review Board of the Charité – Berlin (EA4/104/11). The immunostaining carried out using the Leica Bond-maX autostainer (Leica Biosystems, Illinois, United States) according to the manufacturer’s protocol. After heat-induced epitope retrieval, the sections were incubated with anti-c-myc and anti-CD97 (ADGRE5) rabbit antibodies, respectively (dilution 1:200). Horseradish peroxidase-labeled Anti-rabbit-IgG using the Bond Polymer Refine Detection Kit (Leica Biosystems, Illinois, United States) was employed to convert the chromogen substrate. Staining was performed with appropriate positive and negative controls.
First, we determined MYC mRNA and MYC protein expression by qRT-PCR, Western blotting and immunohistochemistry, respectively, in cell lines derived from BL, DLBCLpos and DLBCLneg patients (Fig. 1 b-d). With the exception of BL-41, all MYC break positive cell lines showed high expression of MYC mRNA. The level of MYC protein expression corresponded without exception with the presence of MYC breaks. The discrepant results between MYC RNA and MYC protein expression in BL-41 might reflect a longer half-life time of the MYC protein in BL-41 as compared to the other cell lines with MYC breaks [23,24,25,26,27]. Thus, less RNA is required to generate high amounts of MYC protein.
To investigate the MYC DNA-binding capabilities in BL and DLBCL, we performed MYC ChIP-Seq experiments to determine genome-wide MYC DNA-binding sites. Additional ChIP-Seq experiments for trimethylation of histone H3 at lysine 4 (H3K4me3) were carried out in order to locate genomic areas with open chromatin as indicators for potential transcriptional activity of nearby genes [28, 29]. To bioinformatically identify differential MYC DNA-binding sites the DiffBind package  was employed using a pairwise comparison of the cell groups (BL vs. DLBCL; BL vs. DLBCLpos; BL vs. DLBCLneg and DLBCLpos vs. DLBCLneg). Similar differential binding analysis was performed with H3K4me3 ChIP-Seq data to ascertain genome wide differential histone patterns and potential active transcriptional sites.
Detailed results of the bioinformatics analyses are available in Additional file 2: File S2(ChIP-Seq data), while Figs. 2 and 3 depict aggregated data. The overall number of MYC DNA-binding sites was higher (approx. 2-fold) in MYC break positive (BL, DLBCLpos) than MYC break negative (DLBCLneg) cells (Fig. 2a). Next, we explored whether genes associated with MYC-binding differed between the three groups of cell lines. Our data clearly indicate that there is not only a difference in the number of genes but in addition, that also different genes are targeted by MYC and/or H3K4 (Fig. 2b). To identify differential MYC-binding genes we performed a differential peak analysis comparing four combinations: BL vs. DLBCL, BL vs. DLBCLneg, DLBCLpos vs. DLBCLneg and BL vs. DLBCLpos. Cell lines carrying MYC breaks have more genes located in the vicinity of MYC-binding sites which leads to a higher number of differential MYC-binding peaks in relation to MYC break negative cell lines (Fig. 2c). Figure 3 highlights a list of twenty target genes selected that yielded the highest fold changes. The analysis of the MYC-binding motifs of MYC target genes showed an interesting distribution (Fig. 2d) with a preference for non-canonical E-Box motives (approx. 45%), while only 4% carried exclusively the classical canonic E-Box motif (CACGTG) and 19% both motifs. Strikingly, 32% of identified MYC targets genes displayed no known MYC-binding motifs. Non-canonical and/or canonical E-box was present in approx. 68% of MYC target genes, thus corroborating previous studies of global mapping of MYC-binding sites . However, the presence of E-box motives in the binding loci did not correlate with the regulation of associated genes [31, 32].
MYC-binding is not the sole factor for activating of gene expression and associated functional consequences. In order to gain a deeper insight into the transcriptional landscape and the impact of the various MYC-binding patterns, we performed RNA-Seq and correlated the results with the presence of MYC breaks and with the MYC and H3K4 binding patterns. In Fig. 4a the number of genes differentially expressed among the three cell line groups is given. The highest number of differentially expressed genes was found between BL and DLBCL in general, whereby the expression difference between BL and DLBCLpos was the lowest. This demonstrates that MYC has a major impact through activation of the same gene set, which constitutes a significant proportion of the entire transcriptome. In harmony with this notion, the comparison of the RNA-Seq data between BL and DLBCLneg, and DLBCLpos and DLBCLneg revealed very similar numbers of differentially expressed genes. This reinforces the similarity in the gene expression profiles of both types of MYC break positive cell lines. Lists of differentially expressed genes are given in Additional file 3: File S3(RNA-Seq data). In Fig. 4b some differentially expressed genes (from Additional file 3: File S3) are functionally grouped into clustering of differentiation (CD) molecules (B1), integrin molecules (B2) or MYC-related molecules (B3) and visualized as heat maps. Most of the identified CD molecules seem to be upregulated in MYC break positive (BL, DLBCLpos) cell lines compared to MYC break negative (DLBCLneg) cell lines.
To validate the data derived from genome-wide DNA-binding and gene expression, we performed additional gene-specific ChIP (MYC and H3K4me3) and RT-PCR experiments (Fig. 4c). The selection criteria for the target genes were MYC DNA-binding according to ChIP-Seq and differential expression according to RNA-Seq. Among the identified genes, ZAP-70, ADGRE5, CDK20, GPAM, SMAD1 and TERT were the most interesting. Genes lacking differentially expression such as LARS, FARSA and already described as MYC target genes like BYSL  and NMP1  were selected as positive control. By independent validation assays we were able to confirm the data derived from our ChIP-Seq and RNA-Seq analyses and demonstrate that a set of genes is able to reliably differentiate between lymphoma entities.
Interestingly ADGRE5 (previously designated as CD97 and marked by an arrow in Fig. 4, B1) shows strong selective expression in BL cell lines. On the other hand, some integrin molecules known to be binding partners of ADGRE5 are downregulated in BL as compared to DLBCL.
To validated this interesting outcome we quantify the proteins of ADGRE5 and already known homogeneous expressed MYC targets like BYSL and NPM1, obtained from previously published proteomic data  (Fig. 5 a) and western blot analysis (Fig. 5 c). Finally, we demonstrated the discriminating character of ADGRE5 between BL and DLBCL in additional immunostainings of cell lines and FFPE tissue samples (Fig. 5 d).
Our RNA-Seq and ChIP-Seq data showed a significant overrepresentation of ADGRE5 in BL as compared to DLBCL regardless of the MYC break status of the latter. This finding was confirmed by independent additional target-specific ChIP experiments, RT-PCR and re-analysis of published proteomic data  (see Figs. 4 C1-2 and 5 a). For further exploration, we selected ADGRE5 as an interesting candidate gene. ADGRE5 previously designated as CD97  belongs to the adhesion G protein-coupled receptor (GPCR) subfamily E and was the first receptor of this gene family found to be associated with cancer . ADGRE5 is a member of the EGF-TM7 (seven-span transmembrane protein containing epidermal growth factor domains; Fig. 5 b) protein family and is constitutively expressed in granulocytes, monocytes as well as in subsets of T- and B-cells [37,38,39,40,41,42]. An increased ADGRE5 expression is found in some types of leukemia [43,44,45,46]. Interestingly, differential expression of ADGRE5 has also been described for several solid cancers such as lung, thyroid and colorectal carcinomas, indicating a tumor and/or tissue-specific expression pattern [47,48,49]. Upregulation of ADGRE5 is often observed at the invading tumor front as well as in advanced tumor stages. Furthermore, ADGRE5 presents an unfavorable prognostic factor [50,51,52,53,54]. Depending on the cell type and tumor grade, ADGRE5 protein exists in three isoforms resulting from alternative splicing . ADGRE5 protein is cleaved by self-catalytic proteolysis into a large extracellular subunit, which contains three (EGF1,2,5), four (EGF1,2,3,5) or five (EGF1,2,3,4,5) extracellular N-terminal epidermal growth factor (EGF)-like domains, that are coupled to the seven-span transmembrane subunit (TM7) via an extended spacer region [56,57,58]. As a surface receptor molecule, ADGRE5 has the ability to bind ligands of the cellular and extracellular matrix, enhances proteolytic activity of matrix metalloproteinases (MMPs) and thus triggers secretion of chemokines . Finally, ADGRE5 interacts with CD55 [60,61,62,63], the glycosaminoglycan chondroitin sulfate [64, 65], integrin  or CD90  in an isoform-specific manner. Initial functional studies suggest that ADGRE5 is relevant for cell adhesion, migration and invasion [53, 59, 68].
To determine ADGRE5 isotype distribution in aggressive lymphoma, we performed Western blot analyses (Fig. 5c) and found that the short isoform (EGF1,2,5) is homogeneously expressed in all cell lines, while the largest isoform (EGF1,2,3,4,5) was preferentially present in BL cell lines. This is a very striking finding, as the EGF-like repeat 4, which has been reported to interact with chondroitin sulfate, is only found in the largest ADGRE5 isoform [64,65,66]. The interaction of chondroitin sulfate and ADGRE5 (EGF4) mediates cell adhesion and angiogenesis and plays an important role in the interaction of activated T-cells, dendritic cells and macrophages. This observation fits very nicely to features of BL, especially angiogenesis and macrophage attraction.
In order to determine the expression profile of ADGRE5 in primary patient specimens, we performed immunohistochemical staining (IHC). Figure 5d illustrates the higher expression of ADGRE5 on the cell surface of BL cell lines and primary BL patient specimens. In contrast, primary tissue specimens obtained from DLBCL patients and cell lines largely lacks ADGRE5 protein expression irrespectively of the presence of MYC breaks. (Additional IHC staining results are shown in Additional file 4: Figure S4). Table 5E summarizes the ADGRE5 IHC results obtained for 38 patients suffering from BL and DLBCL, respectively. Overall, ADGRE5 is significantly more frequently positive in BL patients (88%) as compared to DLBCLs patients that were mainly negative (80%). Thus, the data derived from our cell line experiments are nicely reflected in primary patient specimens.
Here we describe the impact of MYC in three types of aggressive B-cell non-Hodgkin lymphomas: BL and DLBCL with and without MYC break (DLBCLneg and DLBCLpos, respectively). More MYC-binding sites were found by MYC ChIP-Seq in BL and DLBCLpos as compared to DLBCLneg. Interestingly, MYC was found to be bound to different target genes in BL and DLBCLpos, which is also reflected by their gene expression differences. By combined analyses, ADGRE5 (CD97) was identified as an interesting differentially expressed MYC target gene, an observation confirmed by immunohistochemistry of primary FFPE patient samples. Its expression, in particular that of the largest ADGRE5 isoform (EGF1,2,3,4,5), was significantly higher in BL than DLBCL. Based on the reported function of the EGF4 repeat as a receptor for chondroitin sulfate, we hypothesize that this might contribute to some peculiar features of BL, namely macrophage attraction and angiogenesis, and potentially to the very good responsiveness to treatment. In addition, we suggest ADGRE5 as a marker to discriminate between BL and DLBCL (regardless of the presence of MYC breaks) in patient stratification.
B-cell non-Hodgkin lymphomas
Diffuse large B-cell lymphoma
DLBCL without MYC break
DLBCL with MYC break
seven-span transmembrane protein containing epidermal growth factor domains
False Discovery Rate
G protein-coupled receptor
Irreproducible Discovery Rate
Origin of Replication
Transcription start site
Sweetenham JW, Pearce R, Taghipour G, Blaise D, Gisselbrecht C, Goldstone AH. Adult Burkitt's and Burkitt-like non-Hodgkin's lymphoma - outcome for patients treated with high-dose therapy and autologous stem-cell transplantation in first remission or at relapse: results from the European Group for Blood and Marrow Transplantation. J Clin Oncol. 1996;14(9):2465–72.
Aquino G, Marra L, Cantile M, De Chiara A, Liguori G, Curcio MP, Sabatino R, Pannone G, Pinto A, Botti G, et al. MYC chromosomal aberration in differential diagnosis between Burkitt and other aggressive lymphomas. Infect Agents Cancer. 2013;8.
Chan WC, Armitage JO, Gascoyne R, Connors J, Close P, Jacobs P, Norton A, Lister TA, Pedrinis E, Cavalli F, et al. A clinical evaluation of the international lymphoma study group classification of non-Hodgkin's lymphoma. Blood. 1997;89(11):3909–18.
Copie-Bergman C, Cuilliere-Dartigues P, Baia M, Briere J, Delarue R, Canioni D, Salles G, Parrens M, Belhadj K, Fabiani B, et al. MYC-IG rearrangements are negative predictors of survival in DLBCL patients treated with immunochemotherapy: a GELA/LYSA study. Blood. 2015;126(22):2466–74.
Hummel M, Bentink S, Berger H, Klapper W, Wessendorf S, Barth TF, Bernd HW, Cogliatti SB, Dierlamm J, Feller AC, et al. A biologic definition of Burkitt's lymphoma from transcriptional and genomic profiling. N Engl J Med. 2006;354(23):2419–30.
Savage KJ, Johnson NA, Ben-Neriah S, Connors JM, Sehn LH, Farinha P, Horsman DE, Gascoyne RD. MYC gene rearrangements are associated with a poor prognosis in diffuse large B-cell lymphoma patients treated with R-CHOP chemotherapy. Blood. 2009;114(17):3533–7.
Stein H, Hummel M. Burkitt's and burkitt-like lymphoma. Molecular definition and value of the World Health Organisation's diagnostic criteria. Pathologe. 2007;28(1):41–5.
Visco C, Tzankov A, Xu-Monette ZY, Miranda RN, Tai YC, Li Y, Liu WM, d'Amore ESG, Li Y, Montes-Moreno S, et al. Patients with diffuse large B-cell lymphoma of germinal center origin with BCL2 translocations have poor outcome, irrespective of MYC status: a report from an international DLBCL rituximab-CHOP consortium program study. Haematologica. 2013;98(2):255–63.
Quentmeier H, Amini RM, Berglund M, Dirks WG, Ehrentraut S, Geffers R, Macleod RA, Nagel S, Romani J, Scherr M, et al. U-2932: two clones in one cell line, a tool for the study of clonal evolution. Leukemia. 2013;27(5):1155–64.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.
Lee TI, Johnstone SE, Young RA. Chromatin immunoprecipitation and microarray-based analysis of protein location. Nat Protoc. 2006;1(2):729–48.
Seitz V, Butzhammer P, Hirsch B, Hecht J, Gutgemann I, Ehlers A, Lenze D, Oker E, Sommerfeld A, von der Wall E, et al. Deep sequencing of MYC DNA-binding sites in Burkitt lymphoma. PLoS One. 2011;6(11):e26837.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.
Li QH, Brown JB, Huang HY, Bickel PJ. Measuring reproducibility of high-throughput experiments. Ann Appl Stat. 2011;5(3):1752–79.
Ross-Innes CS, Stark R, Teschendorff AE, Holmes KA, Ali HR, Dunning MJ, Brown GD, Gojis O, Ellis IO, Green AR, et al. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature. 2012;481(7381):389–U177.
Taruttis F, Feist M, Schwarzfischer P, Gronwald W, Kube D, Spang R, Engelmann JC. External calibration with Drosophila whole-cell spike-ins delivers absolute mRNA fold changes from human RNA-Seq and qPCR data. Biotechniques. 2017;62(2):53–61.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Schwarzfischer P, Reinders J, Dettmer K, Kleo K, Dimitrova L, Hummel M, Feist M, Kube D, Szczepanowski M, Klapper W, et al. Comprehensive metaboproteomics of Burkitt's and diffuse large B-cell lymphoma cell lines and primary tumor tissues reveals distinct differences in pyruvate content and metabolism. J Proteome Res. 2017.
Ciechanover A, DiGiuseppe JA, Schwartz AL, Brodeur GM. Degradation of MYCN oncoprotein by the ubiquitin system. Prog Clin Biol Res. 1991;366:37–43.
Dani C, Blanchard JM, Piechaczyk M, El Sabouty S, Marty L, Jeanteur P. Extreme instability of myc mRNA in normal and transformed human cells. Proc Natl Acad Sci U S A. 1984;81(22):7046–50.
Eick D, Piechaczyk M, Henglein B, Blanchard JM, Traub B, Kofler E, Wiest S, Lenoir GM, Bornkamm GW. Aberrant c-myc RNAs of Burkitt's lymphoma cells have longer half-lives. EMBO J. 1985;4(13B):3717–25.
Gregory MA, Hann SR. C-Myc proteolysis by the ubiquitin-proteasome pathway: stabilization of c-Myc in Burkitt's lymphoma cells. Mol Cell Biol. 2000;20(7):2423–35.
Hann SR. Role of post-translational modifications in regulating c-Myc proteolysis, transcriptional activity and biological function. Semin Cancer Biol. 2006;16(4):288–302.
Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K. High-resolution profiling of histone methylations in the human genome. Cell. 2007;129(4):823–37.
Yoo CB, Jones PA. Epigenetic therapy of cancer: past, present and future. Nat Rev Drug Discov. 2006;5(1):37–50.
Zeller KI, Zhao X, Lee CW, Chiu KP, Yao F, Yustein JT, Ooi HS, Orlov YL, Shahab A, Yong HC, et al. Global mapping of c-Myc binding sites and target gene networks in human B cells. Proc Natl Acad Sci U S A. 2006;103(47):17834–9.
Hulf T, Bellosta P, Furrer M, Steiger D, Svensson D, Barbour A, Gallant P. Whole-genome analysis reveals a strong positional bias of conserved dMyc-dependent E-boxes. Mol Cell Biol. 2005;25(9):3401–10.
Orian A, van Steensel B, Delrow J, Bussemaker HJ, Li L, Sawado T, Williams E, Loo LW, Cowley SM, Yost C, et al. Genomic binding by the Drosophila Myc, max, mad/Mnt transcription factor network. Genes Dev. 2003;17(9):1101–14.
Fukuda MN, Miyoshi M, Nadano D. The role of bystin in embryo implantation and in ribosomal biogenesis. Cell Mol Life Sci. 2008;65(1):92–9.
Zeller KI, Haggerty TJ, Barrett JF, Guo Q, Wonsey DR, Dang CV. Characterization of nucleophosmin (B23) as a Myc target by scanning chromatin immunoprecipitation. J Biol Chem. 2001;276(51):48285–91.
Hamann J, Aust G, Arac D, Engel FB, Formstone C, Fredriksson R, Hall RA, Harty BL, Kirchhoff C, Knapp B, et al. International Union of Basic and Clinical Pharmacology. XCIV. Adhesion G protein-coupled receptors. Pharmacol Rev. 2015;67(2):338–67.
Aust G, Eichler W, Laue S, Lehmann I, Heldin NE, Lotz O, Scherbaum WA, Dralle H, Hoang-Vu C. CD97: a dedifferentiation marker in human thyroid carcinomas. Cancer Res. 1997;57(9):1798–806.
Ancuta P, Liu KY, Misra V, Wacleche VS, Gosselin A, Zhou X, Gabuzda D. Transcriptional profiling reveals developmental relationship and distinct biological functions of CD16+ and CD16- monocyte subsets. BMC Genomics. 2009;10:403.
Eichler W, Hamann J, Aust G. Expression characteristics of the human CD97 antigen. Tissue Antigens. 1997;50(5):429–38.
Gasz B, Lenard L, Benko L, Borsiczky B, Szanto Z, Lantos J, Szabados S, Alotti N, Papp L, Roth E. Expression of CD97 and adhesion molecules on circulating leukocytes in patients undergoing coronary artery bypass surgery. Eur Surg Res. 2005;37(5):281–9.
Jaspars LH, Vos W, Aust G, Van Lier RA, Hamann J. Tissue distribution of the human CD97 EGF-TM7 receptor. Tissue Antigens. 2001;57(4):325–31.
Kop EN, Matmati M, Pouwels W, Leclercq G, Tak PP, Hamann J. Differential expression of CD97 on human lymphocyte subsets and limited effect of CD97 antibodies on allogeneic T-cell stimulation. Immunol Lett. 2009;123(2):160–8.
Veninga H, Becker S, Hoek RM, Wobus M, Wandel E, van der Kaa J, van der Valk M, de Vos AF, Haase H, Owens B, et al. Analysis of CD97 expression and manipulation: antibody treatment but not gene targeting curtails granulocyte migration. J Immunol. 2008;181(9):6574–83.
Bonardi F, Fusetti F, Deelen P, van Gosliga D, Vellenga E, Schuringa JJ. A proteomics and transcriptomics approach to identify leukemic stem cell (LSC) markers. Mol Cell Proteomics. 2013;12(3):626–37.
Coustan-Smith E, Song G, Clark C, Key L, Liu P, Mehrpooya M, Stow P, Su X, Shurtleff S, Pui CH, et al. New markers for minimal residual disease detection in acute lymphoblastic leukemia. Blood. 2011;117(23):6267–76.
Maiga A, Lemieux S, Pabst C, Lavallee VP, Bouvier M, Sauvageau G, Hebert J. Transcriptome analysis of G protein-coupled receptors in distinct genetic subgroups of acute myeloid leukemia: identification of potential disease-specific targets. Blood Cancer J. 2016;6:e431.
Mirkowska P, Hofmann A, Sedek L, Slamova L, Mejstrikova E, Szczepanski T, Schmitz M, Cario G, Stanulla M, Schrappe M, et al. Leukemia surfaceome analysis reveals new disease-associated features. Blood. 2013;121(25):E149–59.
Aust G, Zhu D, Van Meir EG, Xu L. Adhesion GPCRs in tumorigenesis. Handb Exp Pharmacol. 2016;234:369–96.
Boltze C, Schneider-Stock R, Aust G, Mawrin C, Dralle H, Roessner A, Cuong HV. CD97, CD95 and Fas-L clearly discriminate between chronic pancreatitis and pancreatic ductal adenocarcinoma in perioperative evaluation of cryocut sections. Pathol Int. 2002;52(2):83–8.
Liu Y, Chen L, Peng SY, Chen ZX, Hoang-Vu C. Role of CD97(stalk) and CD55 as molecular markers for prognosis and therapy of gastric carcinoma patients. J Zhejiang Univ Sci B. 2005;6(9):913–8.
Han SL, Xu C, Wu XL, Li JL, Liu Z, Zeng QQ. The impact of expressions of CD97 and its ligand CD55 at the invasion front on prognosis of rectal adenocarcinoma. Int J Color Dis. 2010;25(6):695–702.
He Z, Wu H, Jiao Y, Zheng J. Expression and prognostic value of CD97 and its ligand CD55 in pancreatic cancer. Oncol Lett. 2015;9(2):793–7.
Safaee M, Clark AJ, Oh MC, Ivan ME, Bloch O, Kaur G, Sun MZ, Kim JM, Oh T, Berger MS, et al. Overexpression of CD97 confers an invasive phenotype in glioblastoma cells and is associated with decreased survival of glioblastoma patients. PLoS One. 2013;8(4):e62765.
Steinert M, Wobus M, Boltze C, Schutz A, Wahlbuhl M, Hamann J, Aust G. Expression and regulation of CD97 in colorectal carcinoma cell lines and tumor tissues. Am J Pathol. 2002;161(5):1657–67.
Wu JS, Lei L, Wang SC, Gu DH, Zhang JH. Immunohistochemical expression and prognostic value of CD97 and its ligand CD55 in primary gallbladder carcinoma. J Biomed Biotechnol. 2012;2012:58767.
Bjarnadottir TK, Geirardsdottir K, Ingemansson M, Mirza MA, Fredriksson R, Schioth HB. Identification of novel splice variants of adhesion G protein-coupled receptors. Gene. 2007;387(1–2):38–48.
Gray JX, Haino M, Roth MJ, Maguire JE, Jensen PN, Yarme A, StetlerStevenson MA, Siebenlist U, Kelly K. CD97 is a processed, seven-transmembrane, heterodimeric receptor associated with inflammation. J Immunol. 1996;157(12):5438–47.
Hsiao CC, Chen HY, Chang GW, Lin HH. GPS autoproteolysis is required for CD97 to up-regulate the expression of N-cadherin that promotes homotypic cell-cell aggregation. FEBS Lett. 2011;585(2):313–8.
Hsiao CC, Cheng KF, Chen HY, Chou YH, Stacey M, Chang GW, Lin HH. Site-specific N-glycosylation regulates the GPS auto-proteolysis of CD97. FEBS Lett. 2009;583(19):3285–90.
Galle J, Sittig D, Hanisch I, Wobus M, Wandel E, Loeffler M, Aust G. Individual cell-based models of tumor-environment interactions: multiple effects of CD97 on tumor invasion. Am J Pathol. 2006;169(5):1802–11.
Abbott RJ, Spendlove I, Roversi P, Fitzgibbon H, Knott V, Teriete P, McDonnell JM, Handford PA, Lea SM. Structural and functional characterization of a novel T cell receptor co-regulatory protein complex, CD97-CD55. J Biol Chem. 2007;282(30):22023–32.
Chiu PL, Ng BH, Chang GW, Gordon S, Lin HH. Putative alternative trans-splicing of leukocyte adhesion-GPCR pre-mRNAs generates functional chimeric receptors. FEBS Lett. 2008;582(5):792–8.
Hamann J, Vogel B, vanSchijndel GMW, vanLier RAW. The seven-span transmembrane receptor CD97 has a cellular ligand (CD55, DAF). J Exp Med. 1996;184(3):1185–9.
Toomey CB, Cauvi DM, Pollard KM. The role of decay accelerating factor in environmentally induced and idiopathic systemic autoimmune disease. Autoimmune Dis. 2014;2014:452853.
Kwakkenbos MJ, Pouwels W, Matmati M, Stacey M, Lin HH, Gordon S, van Lier RA, Hamann J. Expression of the largest CD97 and EMR2 isoforms on leukocytes facilitates a specific interaction with chondroitin sulfate on B cells. J Leukoc Biol. 2005;77(1):112–9.
Stacey M, Chang GW, Davies JQ, Kwakkenbos MJ, Sanderson RD, Hamann J, Gordon S, Lin HH. The epidermal growth factor-like domains of the human EMR2 receptor mediate cell attachment through chondroitin sulfate glycosaminoglycans. Blood. 2003;102(8):2916–24.
Wang T, Ward Y, Tian L, Lake R, Guedez L, Stetler-Stevenson WG, Kelly K. CD97, an adhesion receptor on inflammatory cells, stimulates angiogenesis through binding integrin counterreceptors on endothelial cells. Blood. 2005;105(7):2836–44.
Wandel E, Saalbach A, Sittig D, Gebhardt C, Aust G. Thy-1 (CD90) is an interacting partner for CD97 on activated endothelial cells. J Immunol. 2012;188(3):1442–50.
Kobayashi T, Mitsuyama K, Yamasaki H, Masuda J, Takedatsu H, Kuwaki K, Yoshioka S, Nagayama K, Sata M. Microarray analyses of peripheral whole blood cells from ulcerative colitis patients: effects of leukocytapheresis. Int J Mol Med. 2013;31(4):789–96.
The authors would like to thank Edda von der Wall, Anke Sommerfeld and Hedwig Lammert for their excellent technical assistance and Dr. Hilmar Quentmeier for providing U-2932 sub clones. Special thanks go to Ulrike Krüger for Sequencing at the NGS core facility of the Berlin-Brandenburg Center for Regenerative Therapies.
This research was financially supported by the BMBF (Federal Ministry for Education and Research, Germany) in the context of the E:bio Consortium Molecular Mechanisms in Malignant Lymphoma with MYC-Deregulation (MMML-MYC-SYS), grant numbers 0316166H and 0316166G. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethics approval and consent to participate
The use of human primary tissue samples for the Immunohistostaining was approved by the Institutional Review Board of the Charité – Berlin (EA4/104/11).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Compilation of antibodies used for Chromatin Immunoprecipitation (ChIP), Western Blot (WB) or Immunohistochemistry (IHC) and list of TaqMan Assays used for RT-PCR analysis. *indicated endogenous control. (PDF 42 kb)
File S2. ChIP-Seq data of differentially bound MYC and H3K4 target genes derived from comparison of BL vs. DLBCL, BL vs. DLBCLneg, DLBCLneg vs. DLBCLpos and BL vs. DLBCLpos. (XLSX 1673 kb)
File S3. RNA-Seq data of differentially expressed genes derived from comparison of BL vs. DLBCL, BL vs. DLBCLneg, DLBCLneg vs. DLBCLpos and BL vs. DLBCLpos. (XLSX 13003 kb)
Figure S4. ADGRE5 IHC staining of cell line and patient tissue samples. (A) Burkitt Lymphoma (BL): Cell lines (1–5: Blue-1 Bl-41, BL-2, DG-75, CA-46) and primary tumor tissues (7–14) are manly ADGRE5 positive. (B) Diffuse large B cell lymphoma without MYC break (DLBCLneg): Cell lines (1–4: Karpass-422, U2932-R1, HT, WSU-DLCL2) and primary tumor tissues (5–19). (C) DLBCLpos: Cell lines (1–3: Carnaval, U2932-R2, SU-DHL-10) and primary tumor tissues (4–6). DLBCLpos and DLBCLneg are manly negative for ADGRE5. Strong positive staining in some tissue sections results from macrophages or T-cells. (PDF 4769 kb)
About this article
Cite this article
Kleo, K., Dimitrova, L., Oker, E. et al. Identification of ADGRE5 as discriminating MYC target between Burkitt lymphoma and diffuse large B-cell lymphoma. BMC Cancer 19, 322 (2019). https://doi.org/10.1186/s12885-019-5537-0