Skip to main content

Whole-exome sequencing analysis identifies distinct mutational profile and novel prognostic biomarkers in primary gastrointestinal diffuse large B-cell lymphoma



Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive non-Hodgkin lymphoma, and about 10% of DLBCL cases primarily occur in the gastrointestinal tract. Previous reports have revealed that primary gastrointestinal-DLBCL (pGI-DLBCL) harbors different genetic mutations from other nodal or extranodal DLBCL. However, the exonic mutation profile of pGI-DLBCL has not been fully addressed.


We performed whole-exome sequencing of matched tumor tissues and blood samples from 53 pGI-DLBCL patients. The exonic mutation profiles were screened, and the correlations between genetic mutations and clinicopathological characteristics were analyzed.


A total of 6,588 protein-altering events were found and the five most frequent mutated genes in our pGI-DLBCL cohort were IGLL5 (47%), TP53 (42%), BTG2 (28%), P2RY8 (26%) and PCLO (23%). Compared to the common DLBCL, significantly less or absence of MYD88 (0%), EZH2 (0%), BCL2 (2%) or CD79B (8%) mutations were identified in pGI-DLBCL. The recurrent potential driver genes were mainly enriched in pathways related to signal transduction, infectious disease and immune regulation. In addition, HBV infection had an impact on the mutational signature in pGI-DLBCL, as positive HBsAg was significantly associated with the TP53 and LRP1B mutations, two established tumor suppressor genes in many human cancers. Moreover, IGLL5 and LRP1B mutations were significantly correlated with patient overall survival and could serve as two novel prognostic biomarkers in pGI-DLBCL.


Our study provides a comprehensive view of the exonic mutation profile of the largest pGI-DLBCL cohort to date. The results could facilitate the clinical development of novel therapeutic and prognostic biomarkers for pGI-DLBCL.


The incident rate of non-Hodgkin lymphomas (NHLs) in most contries has considerably increased in recent years [1]. Diffuse large B-cell lymphoma (DLBCL) is the most common subtype of NHLs, accounting for nearly one-third of all lymphoid neoplasm in China annually [2, 3]. Though at least two DLBCL subtypes have been identified by RNA expression profiles, the germinal center B-cell-like (GCB) subtype and the activated B-cell-like (ABC) subtype, DLBCL still represents a clinical heterogenous disease due to its complex and diverse histological characteristics [4, 5]. DLBCL patients often present with an aggressive clinical behavior, but most of them can be cured by the standard regimen based on rituximab plus cyclophosphomide, doxorubicin, vincristine and prednisone (R-CHOP) [6]. The application of next-generation sequencing has helped reveal a deep degree of molecular and genetic heterogeneity in hematological diseases, and confirmed that genetic aberrations contribute to occurrence and progression of DLBCL [7, 8].

DLBCL arises from extranodal organs in about 30% of total cases, and one third of extranodal DLBCL cases occur in the gastrointestinal tract, making it the most common primary extranodal site [9, 10]. Patient prognosis and recurrence risk of extranodal DLBCL vary according to the primary site of origin, which may harbor different genetic mutations clarified by high through-put sequencing studies [11, 12]. Primary gastrointestinal DLBCL (pGI-DLBCL) has a significantly decreased level of MYD88 and CD79B mutations compared to nodal DLBCL and other extranodal DLBCL in immune-privileged sites, such as central nervous system and testis [13, 14]. Meanwhile, genetic mutations of MYC or BCL2 rearrangements could be related to the survival and prognosis of pGI-DLBCL patients [15, 16]. The genetic mutation profiles discovered by more in-depth analysis revealed that pGI-DLBCL may have different modes of pathogenesis and progression from non-gastrointestinal DLBCL. Recently, by analyzing a small group of patients using whole-exome sequencing (WES), a study by Li et al. has shed a light on the genetic mutations in pGI-DLBCL [17]. However, comprehensive research focusing on the mutational landscape of pGI-DLBCL, and the correlation between its genetic mutations and clinicopathological features are still rare.

In the present study, we aimed to derive the predictive mutational profile by performing capture-based targeted WES on 53 Chinese pGI-DLBCL patients. The association between clinical characteristics and genetic alterations was also explored. In addition, we tried to identify genetic mutations possibly affecting patient survival and their underlying mechanisms. Our study provided a deeper insight into the genetic features of pGI-DLBCL, which may be helpful to clarify the lymphomagenesis process and develop putative therapeutic and prognostic biomarkers for this disease.

Materials and methods

Patient Cohort

Fifty-three patients diagnosed with pGI-DLBCL according to the criteria defined by Lewin et al. [18] were recruited in this study. All patients underwent partial gastrectomy or enterectomy plus R-CHOP based therapy in our hospital spanning from January 1, 2011 to July 21, 2021. Forty-six surgical resection specimens, seven biopsy specimens and matched patient peripheral blood mononuclear cells (PBMCs) were used for sequencing study. All specimens were reviewed by two independent hematopathologists (Yan Huang and Hai-Ling Liu) according to the 2017 World Health Organization classification criteria [19]. The corresponding medical records of all patients were reviewed to obtain the clinicopathological information. The study was approved by the institutional review board at the Sixth Affiliated Hospital of Sun Yat-Sen University.


Tumor DNA was isolated from five 5-μm-thick sections of formalin-fixed paraffin-embedded tumor tissues with a minimum of 70% neoplastic cells using QIAamp FFPE DNA Tissue Kit (Qiagen, USA), and the paired normal control DNA of PBMCs was extracted with DNeasy Tissue and Blood Kit (Qiagen, USA) according to the manufacturer’s instructions. Degradation and contamination were monitored on a 1% agarose gel, and the concentration was measured by using a Qubit® DNA Assay Kit in a Qubit® 2.0 Fluorometer (Life Technologies, USA). Qualified genomic DNA from tumors and matched PBMCs from 53 pGI-DLBCL patients were fragmented by Covaris technology with resultant library fragments of 180–280 bp, and then adapters were ligated to both ends of the fragments. Extracted DNA was then amplified by ligation-mediated PCR (LM-PCR), purified, and hybridized to the Agilent SureSelect Human Exome V6 (Santa Clara, USA) for enrichment, and nonhybridized fragments were then washed out. Both uncaptured and captured LM-PCR products were subjected to real-time PCR to estimate the magnitude of enrichment. Each captured library was then loaded onto the Illumina HiSeq X platform (Hangzhou Jichenjunchuang Medical Laboratory Co., Ltd, Beijing, China). We performed high-throughput sequencing for each captured library independently. Tumor and normal DNA samples were sequenced to an average depth of > 100 × and > 40 × in targeted exonic regions, respectively.

Genomic analysis

After generating raw data through base calling, paired-end reads were trimmed to remove stretches of low-quality bases (< Q10) and adapters in the sequences. The clean reads were mapped to NCBI Build 37 (hg19) using BWA-0.7.12 mem with the default settings. SAMtools-1.2 was used to sort and index all the BAM files; PicardTools-1.119 was used to remove the duplicates; and GATK-3.3–0 was used for InDel realignment and base quality score recalibration. MuTect-1.1.4 and Strelka were used to call somatic SNVs and InDels in the paired normal and tumor samples. Variants identified in the 1,000 Genomes database ( with a frequency > 1% (unless they were in the Catalog of Somatic Mutations in Cancer (COSMIC) database) or in the Exome Aggregation Consortium ( with a frequency > 0.1% were discarded from the analysis. Variants with an alternate allele depth < 2 and a frequency < 5% were also excluded. In addition, SNVs and InDels were filtered to remove benign changes predicted by the following predictive software programs, including Polyphen2, MutationTaster, Mutation Assessor, FATHMM, Radial SVM, LR, SIFT, and LRT. ANNOVAR was used to annotate all the somatic mutations after filtering.

Pathway enrichment analysis

Gene clustering analysis of the driver mutations was performed by Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool ( as previously described [20]. Only the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis which evaluates the modules at the functional level of the selected genes was executed. Bonferroni P value < 0.05 was set as the cut-off criterion and regarded as statistically significant.

Statistical analysis

Statistical analysis was performed using R version 4.1.2 and GraphPad Prism version 7 (La Jolla, CA, USA). The Mann–Whitney U test and the Spearman rank correlation test were employed to analyze the relationship between the mutated genes and clinicopathological characteristics. Survival analysis was performed using Kaplan–Meier curves and compared using the log-rank test. Comparative test differences were considered significant if the 2-tailed P value was < 0.05 otherwise indicated.


Clinicopathological characteristics of the pGI-DLBCL patient cohort

The clinicopathological characteristics of the pGI-DLBCL patient cohort were summarized in Table 1 and Additional file 1: Table S1. Of note, we included 53 patients diagnosed with pGI-DLBCL in this study, which consisted of 40 males and 13 females, respectively. Tumors were primarily originated from the stomach of 11 patients, small intestine of 29 patients, or large intestine of 13 patients. Helicobacter pylori (Hp) or hepatitis B virus (HBV) infection was positive in 21 (39.6%) or 11 cases (20.8%), respectively. According to the Hans algorithm, 33 and 20 patients were classified as GCB (62.3%) and non-GCB (37.7%) DLBCL subtypes based on the immunohistochemical features. The cohort included 35 patients in clinical stage I or II, and 18 patients in clinical stage III or IV. By the end of the current study, the follow-up duration of the patients was as long as 128.4 months with 11 dead records.

Table 1 Clinicopathological characteristics of 53 pGI-DLBCL patients

Exonic mutational profile of pGI-DLBCL

We performed WES of patient-derived tumor tissue and matched blood DNA. Collectively, 6,588 protein-altering mutational events spanning 3,229 genes were identified from our patient cohort. Of these, 5,489 were missense variants, 171 were in frame insertions or deletions, 394 were frameshift variants, 187 were splice site mutations, 23 were start lost mutations, 13 were stop lost mutations, and 311 were stop gain mutations. The spectrum of the top 40 frequently mutated genes was presented in Fig. 1 and the mutational profile of the entire cohort was summarized in Additional file 2: Table S2. The gene with the highest mutation rate was IGLL5 (mutated in 47% pGI-DLBCL patients), which is also the top 1 mutated gene reported in HBV-related DLBCL [21]. Other most frequently mutated genes (≥ 15%) included TP53, BTG2, P2RY8, PCLO, HIST1H1E, IGHM, KMT2D, CSMD3, MUC16, RYR2, CCND3, DUSP2, FAT4, IGHJ6, CARD11, HIST1H1C, LRP1B, MYC, NBPF1, SI. The genome-wide mutational signatures were also characterized according to the 96 possible mutation types [22]. Three highly confident mutational signatures were extracted from our patient cohort. Of these 3 mutation signatures, signatures 1 and 3 were fitted with COSMIC signature 1 and 26, which have been linked to age and defective DNA mismatch repair in cancer, respectively. Meanwhile signature 2, which was mainly characterized by T to G mutations, was not correlated with any COSMIC signature (Fig. 2).

Fig. 1
figure 1

Top 40 mutated genes in 53 pGI-DLBCL patients. The bar graph on the top indicates the absolute number of exonic mutations in each patient. Top 40 frequently mutated genes constitute the individual rows and are arranged by their mutation rates displayed on the right. Each column represents a patient and each row represents a gene. The histogram on the right shows the number of mutations in each gene. The tracks at the bottom provide information on gender, the molecular subtype sorted by Hans algorithm, the primary tumor sites and the IPI that are color-coded as indicated in the legend. TMB: tumor mutational burden

Fig. 2
figure 2

Major mutational signatures were identified according to the alphabetical 96-substitution classifications from 53 pGI-DLBCL patients. The probability bars for the six types of substitutions are displayed in different colors. The mutation types are on the horizontal axes, whereas vertical axes differ between individual signatures for visualization of their patterns and indicate the percentage of mutations attributed to specific mutation types

Potential driver mutations in pGI-DLBCL

In order to identify potential driver mutations in pGI-DLBCL, we compared the mutation profile of our patient cohort with those pathogenic genes associated with human tumors, which have been published and indexed in the COSMIC, MDG125 [23], SMG127 [24], CDG291 datasets [25]. A total of 417 potential driver genes were identified (Table 2). Among these genes, 30 commonly mutated driver genes were found in at least 5 pGI-DLBCL patients, including TP53, P2RY8, KMT2D, MUC16, CSMD3, FAT4, CCND3, HIST1H1C, CARD11, MYC, LRP1B, B2M, TET2, FOXO1, EBF1, BTG1, SETD1B, BCR, COL3A1, DDX3X, AHNAK2, PIM1, ID3, DNM2, PTPN6, FAT1, ROBO2, NFKBIA, BCL7A, SGK1. Next, we used those potential driver genes shared by at least 2 pGI-DLBCL patients to perform gene clustering analysis with the aid of DAVID algorithm. The result revealed that these recurrent driver genes were mainly enriched in pathways related to human cancers, signal transduction, cell metabolism, infection disease and immune regulation. Important signal transduction pathways were substantially affected such as thyroid hormone signaling, central carbon metabolism, HBV infection, FoxO signaling and B cell receptor signaling (Fig. 3 and Additional file 3: Table S3). These results indicated that abnormal signal transduction cascades, altered cell metabolism and virus infection may jointly contribute to the pathogenesis of pGI-DLBCL.

Table 2 Potential driver mutations in pGI-DLBCL
Fig. 3
figure 3

Heatmap of potential oncogenic pathways affected by exonic mutations in 53 pGI-DLBCL patients. A Thyroid hormone signaling pathway. B Central carbon metabolism in cancer. C Hepatitis B. D FoxO signaling pathway. E B cell receptor signaling pathway. The mutation rate of each gene is displayed on the right of each row. The histogram on the right shows the number of mutations in each gene

Associations between clinicopathological characteristics and exonic mutations in pGI-DLCBL pateints

We analyzed the correlations between the status of top 30 mutated genes and the clinicopathological characteristics, such as age, gender, Hp or HBV infection, LDH level, Eastern Cooperative Oncology Group (ECOG) score, B symptoms, International Prognostic Index (IPI), tumor stage, etc. The result was displayed in Fig. 4, and the correlations with statistical significance were summarized in Additional file 4: Table S4. Interestingly, younger patients tended to have FAT4 and FOXO1 mutations, and patients with non-GCB tumors were correlated with CARD11 mutations. Hp infection showed no association with any parameter, however, HBV infection seemed to be related to certain mutations in pGI-DLBCL, as positive HBsAg was significantly associated with the mutations of TP53 and LRP1B, two important tumor suppressor genes (TSGs) reported in many human cancers (Fig. 5A, B). Moreover, HBsAg positive pGI-DLBCL patients have a significant shorter overall survival (OS), when compared to those without HBV infection (Fig. 5C). These results indicated that genetic mutations in pGI-DLBCL patients were associated with certain clinicopathological parameters, and HBV infection could possibly cause worse prognosis due to mutation in TSGs.

Fig. 4
figure 4

The Spearman correlation matrix between major clinicopathological parameters and the status of top 30 mutated genes across 53 pGI-DLBCL patients. The correlations were obtained by deriving Spearman's correlation coefficients. Red represents a positive correlation and blue represents a negative correlation. The cross mark in each box denotes that the correlation did not reach statistical significance

Fig. 5
figure 5

HBV infection was associated with certain mutations and patient OS in pGI-DLBCL. A, B The bar graph indicates the Spearman’s correlation between HBsAg and TP53 (A) or LRP1B (B) mutation. The stacked percentage for each group is shown and the number in the bar denotes patient number count for each group. C OS for pGI-DLBCL patients stratified by HBsAg status

Mutations correlated with patient survival in pGI-DLBCL

In order to find potential genetic mutations with predictive value for patient OS, we performed survival analysis with the top 30 mutated genes in our pGI-DLBCL patient cohort. Most of the mutated genes were not significantly associated with patient OS. However, we did observe that patients with IGLL5 mutations presented with a better OS, and LRP1B mutations led to a shorter OS (Fig. 6A). A large proportion of the mutations in IGLL5 were missense variants located at its N-terminus uncharacterized domains, and the LRP1B mutations were all missense variants evenly distributed across the entire protein structure (Fig. 6B and Additional file 5: Table S5). How these mutations affect individual gene function and the patient survival needs further exploration.

Fig. 6
figure 6

IGLL5 and LRP1B mutations were correlated with OS in pGI-DLBCL. A OS for pGI-DLBCL patients stratified by IGLL5 (upper panel) or LRP1B (lower panel) mutation. B Lollipop plots with the distribution of somatic mutations on the linear protein and domains of IGLL5 (upper panel) or LRP1B (lower panel) in pGI-DLBCL. Each lollipop denotes a unique mutation location, and its height represents the number of observed mutations. Colored bars indicate the individual protein domains. The type of the mutation is indicated in the legend


In the current study, we performed WES of the largest cohort of pGI-DLBCL to date and identified putative cancer driver mutations and their enriched signaling pathways. We also revealed that HBV infection had an impact on the exonic mutation profile pGI-DLBCL, and mutations of IGLL5 and LRP1B genes could predict patient survival, which to our knowledge, was previously unreported by others.

In accordance with the previous reports [17], our analysis of the pGI-DLBCL exome confirmed the high prevalence of mutations in the cell cycle and apoptosis regulatory pathway, with potential tumor driver mutations in TP53 (22/53), CCND3 (9/53) and MYC (8/53) in over 60% patients. TP53 mutations displayed a significantly increased frequency and MYD88 (0/53), NFKBIE (4/53) or CD79B (4/53) mutations were less or not found in our pGI-DLBCL cohort, suggesting that the pathogenesis of pGI-DLBCL were different from the nodal or other extranodal DLBCL, which relies on an activated NF-κB signaling pathway due to the common mutations in the above mentioned MYD88, NFKBIE, or CD79B genes [26]. Furthermore, mutation frequencies of MUC16 (10/53), CSMD3 (10/53), RYR2 (10/53), FAT4 (9/53), TET2 (7/53), EBF1 (7/53) and SETD1B (7/53), which functions at the transcriptional regulation, epigenetic modification or either cellular attachment, were also increased compared to those in common DLBCL according to COSMIC database. Third, we also identified a relatively large proportion of gene mutations, like P2RY8 (14/53), LRP1B (8/53), B2M (7/53), BCR (6/53), that seldom mentioned by other DLBCL sequencing studies but may probably become the oncogenic events by modulating the B cell migratory behavior and signaling activation [27, 28]. Therefore, we hypothesized that the mutation signature of pGI-DLBCL was different from other DLBCL subtypes, and the potential oncogenic driver mutations should be validated by further research.

Another important finding of our study was that HBV infection may affect the mutation spectrum of pGI-DLBCL. We showed that the oncogenic driver mutations were significantly enriched in the HBV regulatory pathway, and patients with positive HBsAg status had a relatively shorter OS and were more likely to carry TP53 and LRP1B mutations, both of which are supposed to function as TSGs during lymphomagenesis process. Previous studies have shown that HBV infection could cause an enhanced rate of mutagenesis and a distinct set of mutation targets in common DLBCL genome [21]. It is worth mentioning that the three genes, namely IGLL5, TP53 and BTG2, are among the top 5 most mutated genes among their and our WES data. Interestingly, LRP1B have been described as a common target gene for HBV integration in liver cancer [29]. In addition, meta-analysis also revealed that patients infected with HBV had a higher risk of developing DLBCL, and those HBsAg-positive DLBCL patients tended to be diagnosed at a younger age with a more advanced clinical stage and worse outcome [30, 31]. Our study presents the first genomic analysis reinforcing the relationship between HBV infection and the mutation signature of pGI-DLBCL. However, further investigations are needed to verify the interactive mechanism between HBV integration and pGI-DLBCL genome, and how the HBV-related mutations affect the pathogenesis and development processes of pGI-DLBCL disease.

Highlighting the clinical significance of our finding, we identified that two recurrent mutations, IGLL5 and LRP1B, could serve as prognostic biomarkers for pGI-DLBCL patients. Although the function of IGLL5 has not been clarified, pervious reports have shown that it was commonly mutated in DLBCL [32, 33] and is homologous to IGLL1, a gene which is critical for B-cell development [34]. In chronic lymphocytic leukaemia (CLL), IGLL5 mutations were associated with a trend towards decreased overall gene expression, and patients bearing IGLL5 mutations were suggestive for the low-risk of CLL [35], which to some extent, was consistent to our result showing that IGLL5 mutated pGI-DLBCL patients had a better OS. On the other hand, LRP1B is giant membrane molecule that is among the most altered genes in human malignancies [36]. Functional studies have confirmed that LRP1B expression in cancer cells could reduce in vitro cell proliferation and migration abilities, and also suppress in vivo tumorigenicity in mouse models [37,38,39,40]. Genetic alteration events, such as deletions, point mutations or frameshift mutations commonly led to the inactivation of this TSG [41,42,43]. Therefore, it is speculated that LRP1B mutations found in our pGI-DLBCL cohort was associated with the impairment of its gene function, which could cause inferior result on disease progression. Despite we first propose that mutations of IGLL5 and LRP1B were significantly related to the survival of pGI-DLBCL patients, there is still a lack of detailed information on how the mutations affect their expression and/or functional role. Some research suggested that Tumor mutation burden estimated by cancer gene panels (CGPs) could be a potential predictor for prognostic stratification of Chinese DLBCL patients [44]. However, IGLL5 and LRP1B discovered in our study as potential biomarkers for the therapeutics or prognosis of pGI-DLBCL remain to be fully elucidated.

In summary, we performed a comprehensive analysis of the exonic mutation profile of the largest pGI-DLBCL cohort to date, which was characterized by an increased mutation frequency in TP53 and MYC, and a decrease rate or absence of MYD88 or CD79B alteration. We also revealed that HBV infection was related to the mutational signature and patient prognosis of pGI-DLBCL. IGLL5 and LRP1B could serve as predictive biomarkers for patient survival. Our study provides a deeper understanding of the genomic information of pGI-DLBCL and could facilitate the clinical development of novel therapeutic and prognostic biomarkers for pGI-DLBCL.

Availability of data and materials

The data that support the findings of this study are available from the corresponding authors upon reasonable request.



Activated B-cell-like


Chronic lymphocytic leukaemia


Catalog of somatic mutations in cancer


Database for annotation, visualization and integrated discovery


Diffuse large B-cell lymphoma;


Eastern cooperative oncology group


Germinal center B-cell-like


Hepatitis B virus


Helicobacter pylori


Insertion or deletion


International prognostic index


Kyoto encyclopedia of genes and genomes


Ligation-mediated PCR


Non-Hodgkin lymphoma


Overall survival


Peripheral blood mononuclear cells


Primary gastrointestinal diffuse large B-cell lymphoma


Rituximab plus cyclophosphomide, doxorubicin, vincristine and prednisone


Single nucleotide variant


Tumor suppressor gene


Whole-exome sequencing


  1. Keykhaei M, Masinaei M, Mohammadi E, Azadnajafabad S, Rezaei N, Saeedi Moghaddam S, Rezaei N, Nasserinejad M, Abbasi-Kangevari M, Malekpour MR, Ghamari SH, Haghshenas R, Koliji K, Kompani F, Farzadfar F. A global, regional, and national survey on burden and quality of care index (QCI) of hematologic malignancies; global burden of disease systematic analysis 1990–2017. Exp Hematol Oncol. 2021;10(1):11.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Susanibar-Adaniya S, Barta SK. 2021 Update on diffuse large B cell lymphoma: a review of current data and potential applications on risk stratification and management. Am J Hematol. 2021;96(5):617–29.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Sun J, Yang Q, Lu Z, He M, Gao L, Zhu M, Sun L, Wei L, Li M, Liu C, Zheng J, Liu W, Li G, Chen J. Distribution of lymphoid neoplasms in China: analysis of 4,638 cases according to the world health organization classification. Am J Clin Pathol. 2012;138(3):429–34.

    Article  PubMed  Google Scholar 

  4. Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A, Boldrick JC, Sabet H, Tran T, Yu X, Powell JI, Yang L, Marti GE, Moore T, Hudson J Jr, Lu L, Lewis DB, Tibshirani R, Sherlock G, Chan WC, Greiner TC, Weisenburger DD, Armitage JO, Warnke R, Levy R, Wilson W, Grever MR, Byrd JC, Botstein D, Brown PO, Staudt LM. Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature. 2000;403(6769):503–11.

    Article  CAS  PubMed  Google Scholar 

  5. Schmitz R, Wright GW, Huang DW, Johnson CA, Phelan JD, Wang JQ, Roulland S, Kasbekar M, Young RM, Shaffer AL, Hodson DJ, Xiao W, Yu X, Yang Y, Zhao H, Xu W, Liu X, Zhou B, Du W, Chan WC, Jaffe ES, Gascoyne RD, Connors JM, Campo E, Lopez-Guillermo A, Rosenwald A, Ott G, Delabie J, Rimsza LM, Tay Kuang Wei K, Zelenetz AD, Leonard JP, Bartlett NL, Tran B, Shetty J, Zhao Y, Soppet DR, Pittaluga S, Wilson WH, Staudt LM. Genetics and Pathogenesis of Diffuse Large B-Cell Lymphoma. N Engl J Med. 2018;378(15):1396–407.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Coiffier B, Lepage E, Briere J, Herbrecht R, Tilly H, Bouabdallah R, Morel P, Van Den Neste E, Salles G, Gaulard P, Reyes F, Lederlin P, Gisselbrecht C. CHOP chemotherapy plus rituximab compared with CHOP alone in elderly patients with diffuse large-B-cell lymphoma. N Engl J Med. 2002;346(4):235–42.

    Article  CAS  PubMed  Google Scholar 

  7. Reddy A, Zhang J, Davis NS, Moffitt AB, Love CL, Waldrop A, Leppa S, Pasanen A, Meriranta L, Karjalainen-Lindsberg ML, Norgaard P, Pedersen M, Gang AO, Hogdall E, Heavican TB, Lone W, Iqbal J, Qin Q, Li G, Kim SY, Healy J, Richards KL, Fedoriw Y, Bernal-Mizrachi L, Koff JL, Staton AD, Flowers CR, Paltiel O, Goldschmidt N, Calaminici M, Clear A, Gribben J, Nguyen E, Czader MB, Ondrejka SL, Collie A, Hsi ED, Tse E, Au-Yeung RKH, Kwong YL, Srivastava G, Choi WWL, Evens AM, Pilichowska M, Sengar M, Reddy N, Li S, Chadburn A, Gordon LI, Jaffe ES, Levy S, Rempel R, Tzeng T, Happ LE, Dave T, Rajagopalan D, Datta J, Dunson DB, Dave SS. Genetic and functional drivers of diffuse large B cell lymphoma. Cell. 2017;171(2):481–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Ren Y, Xiao F, Cheng F, Huang X, Li J, Wang X, Lang W, Zhou X, Lan J, Tong H. Whole exome sequencing reveals a novel LRBA mutation and clonal hematopoiesis in a common variable immunodeficiency patient presented with hemophagocytic lymphohistiocytosis. Exp Hematol Oncol. 2021;10(1):38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Yin X, Xu A, Fan F, Huang Z, Cheng Q, Zhang L, Sun C, Hu Y. Incidence and mortality trends and risk prediction nomogram for extranodal diffuse large b-cell lymphoma: an analysis of the surveillance, epidemiology, and end results database. Front Oncol. 2019;9:1198.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Olszewska-Szopa M, Wrobel T. Gastrointestinal non-Hodgkin lymphomas. Adv Clin Exp Med. 2019;28(8):1119–24.

    Article  PubMed  Google Scholar 

  11. Ollila TA, Olszewski AJ. Extranodal diffuse large b cell lymphoma: molecular features, prognosis, and risk of central nervous system recurrence. Curr Treat Options Oncol. 2018;19(8):38.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Rovira J, Karube K, Valera A, Colomer D, Enjuanes A, Colomo L, Martinez-Trillos A, Gine E, Dlouhy I, Magnano L, Delgado J, Martinez A, Villamor N, Campo E, Lopez-Guillermo A. MYD88 L265P mutations, but no other variants, identify a subpopulation of DLBCL patients of activated B-cell origin, extranodal involvement, and poor outcome. Clin Cancer Res. 2016;22(11):2755–64.

    Article  CAS  PubMed  Google Scholar 

  13. Frick M, Bettstetter M, Bertz S, Schwarz-Furlan S, Hartmann A, Richter T, Lenze D, Hummel M, Dreyling M, Lenz G, Gaumann A. Mutational frequencies of CD79B and MYD88 vary greatly between primary testicular DLBCL and gastrointestinal DLBCL. Leuk Lymphoma. 2018;59(5):1260–3.

    Article  CAS  PubMed  Google Scholar 

  14. Kraan W, van Keimpema M, Horlings HM, Schilder-Tol EJ, Oud ME, Noorduyn LA, Kluin PM, Kersten MJ, Spaargaren M, Pals ST. High prevalence of oncogenic MYD88 and CD79B mutations in primary testicular diffuse large B-cell lymphoma. Leukemia. 2014;28(3):719–20.

    Article  CAS  PubMed  Google Scholar 

  15. Kawajiri A, Maruyama D, Maeshima AM, Nomoto J, Makita S, Kitahara H, Miyamoto KI, Fukuhara S, Suzuki T, Munakata W, Tajima K, Itami J, Taniguchi H, Kobayashi Y, Tobinai K. Impact of the double expression of MYC and BCL2 on outcomes of localized primary gastric diffuse large B-cell lymphoma patients in the rituximab era. Blood Cancer J. 2016;6(9): e477.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. 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.

    Article  CAS  PubMed  Google Scholar 

  17. Li P, Chai J, Chen Z, Liu Y, Wei J, Liu Y, Zhao D, Ma J, Wang K, Li X, Shao Y, Gong L, Zhang W, Guo S, Yan Q, Li M, Fan L, Wang Z. Genomic mutation profile of primary gastrointestinal diffuse large B-cell lymphoma. Front Oncol. 2021;11: 622648.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Lewin KJ, Ranchod M, Dorfman RF. Lymphomas of the gastrointestinal tract: a study of 117 cases presenting with gastrointestinal disease. Cancer. 1978;42(2):693–707.

    Article  CAS  PubMed  Google Scholar 

  19. Swerdlow SH, Campo E, Harris NL, Jaffe ES, Pileri SA, Stein H, Thiele J. 2017. WHO Classification of Tumours of Haematopoietic and Lymphoid Tissues. Revised 4th ed. Lyon: Inernational Agency for Research on Cancer. 2017.

  20. Li SS, Lian YF, Huang YL, Huang YH, Xiao J. Overexpressing PLOD family genes predict poor prognosis in gastric cancer. J Cancer. 2020;11(1):121–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Ren W, Ye X, Su H, Li W, Liu D, Pirmoradian M, Wang X, Zhang B, Zhang Q, Chen L, Nie M, Liu Y, Meng B, Huang H, Jiang W, Zeng Y, Li W, Wu K, Hou Y, Wiman KG, Li Z, Zhang H, Peng R, Zhu S, Pan-Hammarstrom Q. Genetic landscape of hepatitis B virus-associated diffuse large B-cell lymphoma. Blood. 2018;131(24):2670–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Alexandrov LB, Kim J, Haradhvala NJ, Huang MN, Tian Ng AW, Wu Y, Boot A, Covington KR, Gordenin DA, Bergstrom EN, Islam SMA, Lopez-Bigas N, Klimczak LJ, McPherson JR, Morganella S, Sabarinathan R, Wheeler DA, Mustonen V, Group PMSW, Getz G, Rozen SG, Stratton MR, Consortium P. The repertoire of mutational signatures in human cancer. Nature. 2020;578(7793):94–101.

    Article  Google Scholar 

  23. Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA Jr, Kinzler KW. Cancer genome landscapes. Science. 2013;339(6127):1546–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Kandoth C, McLellan MD, Vandin F, Ye K, Niu B, Lu C, Xie M, Zhang Q, McMichael JF, Wyczalkowski MA, Leiserson MDM, Miller CA, Welch JS, Walter MJ, Wendl MC, Ley TJ, Wilson RK, Raphael BJ, Ding L. Mutational landscape and significance across 12 major cancer types. Nature. 2013;502(7471):333–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Tamborero D, Gonzalez-Perez A, Perez-Llamas C, Deu-Pons J, Kandoth C, Reimand J, Lawrence MS, Getz G, Bader GD, Ding L, Lopez-Bigas N. Comprehensive identification of mutational cancer driver genes across 12 tumor types. Sci Rep. 2013;3:2650.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Visco C, Tanasi I, Quaglia FM, Ferrarini I, Fraenza C, Krampera M. Oncogenic mutations of MYD88 and CD79B in diffuse large B-Cell lymphoma and implications for clinical practice. Cancers. 2020;12(10):2913.

    Article  CAS  PubMed Central  Google Scholar 

  27. Young RM, Phelan JD, Wilson WH, Staudt LM. Pathogenic B-cell receptor signaling in lymphoid malignancies: new insights to improve treatment. Immunol Rev. 2019;291(1):190–213.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lu E, Wolfreys FD, Muppidi JR, Xu Y, Cyster JG. S-Geranylgeranyl-L-glutathione is a ligand for human B cell-confinement receptor P2RY8. Nature. 2019;567(7747):244–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Ding D, Lou X, Hua D, Yu W, Li L, Wang J, Gao F, Zhao N, Ren G, Li L, Lin B. Recurrent targeted genes of hepatitis B virus in the liver cancer genomes identified by a next-generation sequencing-based approach. PLoS Genet. 2012;8(12): e1003065.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Chen MH, Hsiao LT, Chiou TJ, Liu JH, Gau JP, Teng HW, Wang WS, Chao TC, Yen CC, Chen PM. High prevalence of occult hepatitis B virus infection in patients with B cell non-Hodgkin’s lymphoma. Ann Hematol. 2008;87(6):475–80.

    Article  PubMed  Google Scholar 

  31. Deng L, Song Y, Young KH, Hu S, Ding N, Song W, Li X, Shi Y, Huang H, Liu W, Zheng W, Wang X, Xie Y, Lin N, Tu M, Ping L, Ying Z, Zhang C, Sun Y, Zhu J. Hepatitis B virus-associated diffuse large B-cell lymphoma: unique clinical features, poor outcome, and hepatitis B surface antigen-driven origin. Oncotarget. 2015;6(28):25061–73.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Evrard SM, Pericart S, Grand D, Amara N, Escudie F, Gilhodes J, Bories P, Traverse-Glehen A, Dubois R, Brousset P, Parrens M, Laurent C. Targeted next generation sequencing reveals high mutation frequency of CREBBP, BCL2 and KMT2D in high-grade B-cell lymphoma with MYC and BCL2 and/or BCL6 rearrangements. Haematologica. 2019;104(4):e154–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. de Miranda NF, Georgiou K, Chen L, Wu C, Gao Z, Zaravinos A, Lisboa S, Enblad G, Teixeira MR, Zeng Y, Peng R, Pan-Hammarstrom Q. Exome sequencing reveals novel mutation targets in diffuse large B-cell lymphomas derived from Chinese patients. Blood. 2014;124(16):2544–53.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Thompson EC, Cobb BS, Sabbattini P, Meixlsperger S, Parelho V, Liberg D, Taylor B, Dillon N, Georgopoulos K, Jumaa H, Smale ST, Fisher AG, Merkenschlager M. Ikaros DNA-binding proteins as integral components of B cell developmental-stage-specific regulatory circuits. Immunity. 2007;26(3):335–44.

    Article  CAS  PubMed  Google Scholar 

  35. Kasar S, Kim J, Improgo R, Tiao G, Polak P, Haradhvala N, Lawrence MS, Kiezun A, Fernandes SM, Bahl S, Sougnez C, Gabriel S, Lander ES, Kim HT, Getz G, Brown JR. Whole-genome sequencing reveals activation-induced cytidine deaminase signatures during indolent chronic lymphocytic leukaemia evolution. Nat Commun. 2015;6:8866.

    Article  CAS  PubMed  Google Scholar 

  36. Principe C, de DionisioSousa IJ, Prazeres H, Soares P, Lima RT. LRP1B: a giant lost in cancer translation. Pharmaceuticals. 2021;14(9):836.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Prazeres H, Torres J, Rodrigues F, Pinto M, Pastoriza MC, Gomes D, Cameselle-Teijeiro J, Vidal A, Martins TC, Sobrinho-Simoes M, Soares P. Chromosomal, epigenetic and microRNA-mediated inactivation of LRP1B, a modulator of the extracellular environment of thyroid cancer cells. Oncogene. 2011;30(11):1302–17.

    Article  CAS  PubMed  Google Scholar 

  38. Lu YJ, Wu CS, Li HP, Liu HP, Lu CY, Leu YW, Wang CS, Chen LC, Lin KH, Chang YS. Aberrant methylation impairs low density lipoprotein receptor-related protein 1B tumor suppressor function in gastric cancer. Genes Chromosom Cancer. 2010;49(5):412–24.

    CAS  PubMed  Google Scholar 

  39. Sonoda I, Imoto I, Inoue J, Shibata T, Shimada Y, Chin K, Imamura M, Amagasa T, Gray JW, Hirohashi S, Inazawa J. Frequent silencing of low density lipoprotein receptor-related protein 1B (LRP1B) expression by genetic and epigenetic mechanisms in esophageal squamous cell carcinoma. Cancer Res. 2004;64(11):3741–7.

    Article  CAS  PubMed  Google Scholar 

  40. Ni S, Hu J, Duan Y, Shi S, Li R, Wu H, Qu Y, Li Y. Down expression of LRP1B promotes cell migration via RhoA/Cdc42 pathway and actin cytoskeleton remodeling in renal cell cancer. Cancer Sci. 2013;104(7):817–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Liu CX, Musco S, Lisitsina NM, Forgacs E, Minna JD, Lisitsyn NA. LRP-DIT, a putative endocytic receptor gene, is frequently inactivated in non-small cell lung cancer cell lines. Cancer Res. 2000;60(7):1961–7.

    CAS  PubMed  Google Scholar 

  42. Tabouret E, Labussiere M, Alentorn A, Schmitt Y, Marie Y, Sanson M. LRP1B deletion is associated with poor outcome for glioblastoma patients. J Neurol Sci. 2015;358(1–2):440–3.

    Article  CAS  PubMed  Google Scholar 

  43. Chen H, Chong W, Wu Q, Yao Y, Mao M, Wang X. Association of LRP1B mutation with tumor mutation burden and outcomes in melanoma and non-small cell lung cancer patients treated with immune check-point blockades. Front Immunol. 2019;10:1113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Chen C, Liu S, Jiang X, Huang L, Chen F, Wei X, Guo H, Shao Y, Li Y, Li W. Tumor mutation burden estimated by a 69-gene-panel is associated with overall survival in patients with diffuse large B-cell lymphoma. Exp Hematol Oncol. 2021;10(1):20.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


This study was supported by National Natural Science Foundation of China (grant number 81902397), Major Talent Training Project of the Third Affiliated Hospital of Sun Yat-sen University (granted to Yi-Fan Lian), General Planned Project of Guangzhou Science and Technology (grant number 202201010950 and 202102080283), Fundamental Research Funds for the Central Universities (grant number 22qntd3401), Key Project of Rural Science and Technology Commissioner of Guangdong Province (grant number KPT20190263), Research Fund of Wu Jieping Medical Foundation (grant number 320.6750.2020-01-32) and Bethune-Tuoyi Young and Middle-aged Doctors' Research Ability Development Project (grant number BQE-TY-SSP(C-3)-S-03).

Author information

Authors and Affiliations



LSS, LYF and XJ conceived and designed the study. LSS, ZXH, LTZ, CTY, CDM, XLX, GXQ, CK and HWJ collected samples and patient information. LSS, ZXH and LHL performed the experiment. LHL and HY reviewed and confirmed the specimens. LSS, HY, LYF and XJ analyzed the data. LSS and LYF wrote the manuscript. LYF and XJ supervised the project and provided funding. All authors contributed to the article and approved the submitted version. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yan Huang, Yi-Fan Lian or Jian Xiao.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Institute Research Ethics Committee at the Sixth Affiliated Hospital of Sun Yat-sen University. Written informed consent was obtained from each patient.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Clinicopathological information of 53 pGI-DLBCL patients.

Additional file 2: Table S2.

Exonic mutation profile of 53 pGI-DLBCL patients.

Additional file 3: Table S3.

KEGG enrichment results of recurrent driver genes in pGI-DLBCL.

Additional file 4: Table S4.

Summary of the statistically significant correlations in the matrix.

Additional file 5: Table S5.

Summary of IGLL5 and LRP1B mutations in pGI-DLBCL.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, SS., Zhai, XH., Liu, HL. et al. Whole-exome sequencing analysis identifies distinct mutational profile and novel prognostic biomarkers in primary gastrointestinal diffuse large B-cell lymphoma. Exp Hematol Oncol 11, 71 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: