Investigation of potential traces of pluripotency in germinal-center-derived B-cell lymphomas driven by MYC
Authors: R Wagener, M Lenz, B Schuldt, I Lenz, A Schuppert, R Siebert, F-J Müller
Journal: Blood Cancer Journal (2015)

Human pluripotent stem cells (hPSCs) are defined by their indefinite self-renewal in vitro and their potential to differentiate into cells representative of all lineages of the mammalian gastrulation embryo.1 Several similarities link these genomically normal cells with features salient to malignancies,2 for example, when injected into immunodeficient mice, hPSCs give rise to tumors closely resembling teratomas.1 hPSCs can be induced from somatic cells by the transient overexpression of four genes (POU5F1, KLF4, SOX2 and MYC). The induction of pluripotency in somatic cells shares again several molecular mechanisms with malignant transformation. Moreover, all of those four reprogramming factors possess oncogenic potential,2 with MYC being the most prominent oncogene activated in a wide range of malignancies. Among them, the prototype of a MYC-driven tumor is Burkitt lymphoma (BL).

BL is the most common B-cell lymphoma in childhood. Its biological hallmark are translocations juxtaposing the MYC oncogene next to one of the three immunoglobulin (IG) loci. Such an IG–MYC fusion resulting in a deregulated MYC expression is present in nearly all BL.3 In contrast to other B-cell lymphomas, in which IG–MYC or non-IG–MYC translocations frequently occur during clonal evolution and disease progression, it is widely accepted that the IG–MYC translocation is the initiating event in Burkitt lymphomagenesis.4

Studies in chicken bursa of Fabricius have shown that myc overexpression induces preneoplastic lesions, which subsequently lead to lymphoma in a cell population solely present in B-cell differentiation at the bursal stem cell stage.5, 6 Moreover, BL cells are exceptional among germinal-center-derived B-cell lymphomas, in that they usually show a dark zone gene expression profile and the dark zone is assumed to be the initiating part in the evolution of germinal centers after antigen stimulation.7 These findings as well as the incidence peak of BL between 4 and 7 years of age and the pluripotency-promoting capacity of MYC led us to the hypothesis that, BL might derive from or resemble a B-lymphoid cell with—potential—stem cell features.

To investigate whether BL cells share characteristics with hPSCs, we used a recently published bioinformatic assay to determine pluripotency in human cells.8 By this unbiased approach, we compared a highly discriminative, transcriptional model derived from a large database of hPSC gene expression profiles with genome-wide transcriptional profiles from various germinal center B-cell lymphoma entities including BL.

Gene expression data from BL, as well as other MYC-positive and MYC-negative mature aggressive B-cell lymphomas obtained by use of the Affymetrix (Santa Clara, CA, USA) U133A platform have been published previously.9 The raw dataset applied for the present analysis ({type:entrez-geo,attrs:{text:GSE4475,term_id:4475}}GSE4475) included a total of 221 mature aggressive B-cell lymphomas classified as molecular BL (mBL, n=44), non-mBLs (n=129) and lymphomas with a gene expression signature intermediate between mBL and non-mBL (intermediate lymphomas, n=48). Moreover, the cases were classified according to their MYC translocation status into those having a translocation juxtaposing MYC with an IG locus (IG–MYC positive), with a non-IG locus (non-IG–MYC positive) or lacking a MYC translocation (MYC negative). For contextualizing our findings, raw gene expression data of 173 primary germ cell tumors (GCTs) ({type:entrez-geo,attrs:{text:GSE3218,term_id:3218}}GSE3218, {type:entrez-geo,attrs:{text:GSE10783,term_id:10783}}GSE10783 and {type:entrez-geo,attrs:{text:GSE18155,term_id:18155}}GSE18155) and from a large expression atlas (E-MTAB-62 (ref. 10)) were downloaded. All raw data sets were preprocessed as described in the Supplementary Materials and Methods section.

Using the PluriTest algorithm, mBL could not be classified as pluripotent (Figure 1a). In comparison, most GCTs and teratomas showed strikingly higher Pluripotency Scores and lower Novelty Scores (Supplementary Figure S1), indicating a closer resemblance of most GCT transcription profiles to hPSC than any of the analyzed lymphomas to hPSCs.



Figure 1

PluriTest analysis. (a) PluriTest results in germ cell tumors, somatic tissues and B-cell lymphoma. The red and blue background encodes an empirical density map indicating the location of pluripotent (red) and non-pluripotent (blue) cells in the reference ...



Interestingly, when we compared the subgroups of the 221 analyzed mature aggressive B-cell lymphomas, we found that mBL scored slightly higher on the Pluripotency Score and lower on the Novelty Score than MYC-negative non-mBL (Figure 1b), effectively separating those two tumor entities in two groups (P-value<10−17) based on the resemblance to transcriptional patterns characteristic for hPSCs. To elucidate whether the translocation and subsequent overexpression of MYC leads to the increased Pluripotency Score and overall similarity to hPSCs as determined by the Novelty Score, non-mBL and intermediate B-cell lymphoma were color-coded based on their MYC translocation status (IG–MYC, non-IG–MYC or MYC negative). Figure 1c shows that the intermediates and non-mBL being MYC-positive mostly do not have a higher Pluripotency Score than the MYC-negative but are intermingled between mBL and MYC-negative non-mBL, suggesting that MYC status alone cannot explain the higher Pluripotency Score and lower Novelty Score of mBL compared with non-mBL.

Nevertheless, we find that in mBL 44.8% of Pluripotency Score signature genes are expressed at comparable levels to embryonic stem cells (Figure 1d). This is in line with recent findings that tumor cells share an overall gene expression pattern with PSC,11 which might be caused by the shared molecular and cellular features as, for example, certain genes drive a rapid proliferation rate both in mBL and hPSC.1 Along these lines, Kim et al.12 have shown that the similarities of hPSCs and some tumor-associated signatures are contributed by a distinct MYC regulatory network configuration. This agrees with our subsequent analyses, where we find that most mBL barely express genes mechanistically linked to pluripotency such as NANOG or POU5F1 (refer to Figure 1d), but express genes that are related to increased cell cycle and proliferative activity as BCL11A or PLCG2 at relatively high levels. Therefore, it is possible that both Pluripotency Score and Novelty Score are shifted toward hPSCs by a MYC-driven proliferation signature.

To summarize, using the PluriTest algorithm we could show that BL cells resemble more hPSCs than non-mBLs and intermediate lymphomas. In spite of this clear and distinguishing shift in the unbiased PluriTest Scores, we were not able to verify our hypothesis that BL might derive from or resemble cells with PSC features.