Overall clinicopathologic and genetic features of FL
Clinicopathologic characteristics of 253 patients with FL were summarized in Table 1. To further explore the genetic features of FL, we performed DNA sequencing on 109 patients (WGS, n = 13; WES, n = 87; and targeted sequencing, n=9). Among them, a median of 59 (37-120) and 58 (3-131) somatic mutations in the coding regions per tumor tissue were identified by WES and WGS, respectively (Additional file 2: Figure S1A-C). A preference for missense mutation (Figure 1A) and G/C>A/T alteration (Figure 1B) were observed, including somatic single nucleotide variations (SNVs) and small insertions/deletions (INDELs), in analogous to the somatic SNV spectrum in other cancers [22]. Fifty-six candidate genes with predicted functional alterations and related to FL pathogenesis were analyzed. As shown in Figure 1C, KMT2D and CREBBP were the most frequently mutated genes, followed by BCL2, TNFRSF14, STAT6, EZH2, CARD11, ARID1A, and EP300. Seventy-one patients were positive for BCL2 translocation. We noted significantly decreased BCL2 translocation (Additional file 2: Figure S1D, p < 0.001) and mutations in KMT2D (Additional file 2: Figure S1E, p = 0.010) and CREBBP (p = 0.042) when compared with the previous literature in Western population [23].
Clinicopathologic and genetic features of HBV-associated FL
Among 253 patients, 17.0% patients (43/253) were HBsAg+ (Table 1). HBV DNA was quantified in all HBsAg+ FL (Additional file 3: Table S1) and below 103 IU/ml before rituximab-containing treatment. As shown in Table 1, HBsAg+ FL presented significantly younger age (p = 0.044), more frequent spleen involvement (p = 0.006), more histological grade FLIIIA (p < 0.001), and higher proliferation index (Ki67) (p < 0.001). As for HBsAg- subgroups, no significant difference was observed between HBsAg-/HBcAb- and HBsAg-/HBcAb+ patients, except that HBsAg-/HBcAb+ patients were older than HBsAg-/HBcAb- patients (Additional file 3: Table S2), probably due to HBV vaccination for newborns since 1992 [24].
Genetic aberrations were compared in 33 HBsAg+ FL and 76 HBsAg- FL. Totally, four clusters of genetic aberrations were defined according to gene functions (Figure 2A). HBsAg+ FL patients were significantly associated with increased mutations in immune modulators including TNFAIP3, CD70, CXCR4, PIM1, KLF2, FAS, CD58, and CD83, and decreased mutations in epigenetic modifiers including KMT2D and CREBBP. In terms of functional clusters (Figure 2B), higher mutations in immune modulators (p < 0.001), but lower mutations in epigenetic modifiers (p < 0.001) were observed in HBsAg+ FL. As for HBsAg- subgroups, no significant difference was observed between HBsAg-/HBcAb- and HBsAg-/HBcAb+ patients in genetic aberrations (Additional file 3: Table S3) or functional clusters (Additional file 3: Table S4).
Besides, decreased BCL2 translocation (p < 0.001) and BCL2 expression (p = 0.003) were detected in HBsAg+ FL patients (Figure 2C). As for HBsAg- subgroups, no significant differences were observed between HBsAg-/HBcAb+ and HBsAg-/HBcAb- patients (Additional file 3: Table S5). Meanwhile, different from those in hepatocellular carcinoma (e.g., viral gene integration into the liver cell genome) [25]. integration of viral genes into the genome of tumor cells was not detected in HBsAg+ FL using WGS, and HBsAg was not detected by immunohistochemistry method (Additional file 2: Figure S2A-S2C).
Enriched immune-related pathways in HBV-associated FL
Using transcriptomic sequencing, gene expression patterns of 17 HBsAg+ FL and 67 HBsAg- FL were compared. Overall, 1514 differentially expressed genes (DEGs) were identified (Additional file 3: Table S6). Enriched Gene Ontology (GO) terms of DEGs revealed that HBsAg+ FL was closely related to immune-related pathways, including T-cell activation, leukocyte differentiation, positive regulation of innate immune response, negative regulation of immune system process, activation of innate immune response, lymphocyte differentiation, antigen processing and presentation, TCR, NIK/NF-κB, positive regulation of leukocyte activation, TNF-mediated, and B-cell activation signaling pathways (Figure 3A). Moreover, ingenuity pathway analysis (IPA) of upstream regulators observed a number of inflammatory and effector cytokine pathways were activated in HBsAg+ FL, including interferon (IFN)-γ, TNF, CD3 complex, IFN-α, NF-κB complex, TCR complex, and B cell receptor (BCR) complex, while PTEN, BCL6, and BACH2 were inhibited (Figure 3B). Gene set enrichment analysis (GSEA) analysis confirmed significantly upregulated TCR signaling pathway (Additional file 2: Figure S3A) and TNF signaling pathway via NF-κB (Additional file 2: Figure S3B) in HBsAg+ FL. Altogether, these findings suggested that HBsAg+ FL exhibited an immune-inflamed phenotype.
Enriched antigen processing and presentation and IFN signatures in HBV-associated FL
As reported, immune-inflamed tumors were usually characterized by enhanced antigen processing and presentation, activation of IFN signatures, and infiltration of CD8+T cells [26, 27]. Interestingly, antigen processing and presentation was indeed enriched in HBsAg+ FL according to GO terms (Figure 3A), and GSEA confirmed this (Additional file 2: Figure S3C). Since MHC molecules were well-known for their role in antigen presentation [28], we evaluated the association of MHC molecules with HBsAg+ FL. Significantly higher expression of MHC class I molecules were observed in HBsAg+ FL (Figure 3C), including HLA-A (p = 0.001) and HLA-C (p = 0.002). As for MHC class II molecules (Figure 3D), significantly higher expression of HLA-DPB1, HLA-DQA1, HLA-DQA2, HLA-DQB1, HLA-DQB2, HLA-DRB1, and HLA-DRB5 was observed in HBsAg+ FL.
As aforementioned in Figure 3B, IFN-γ and IFN-α were activated upstream regulators in HBsAg+ FL. Indeed, expression of IFN signatures were significantly upregulated in HBsAg+ FL (Figure 3E), including 10-gene signature, 28-gene signature, 6-gene signature, and 18-gene signature as reported by Terrill K. McClanahan [29]. The detailed genes of signatures were summarized in Additional file 3: Table S7.
Enriched infiltration of CD8+T and Th1 cells in HBV-associated FL
T-cell subsets were distinguished by the xCell method using transcriptomic sequencing data (Figure 3F) [30]. HBsAg+ FL presented significantly increased total CD8+T cells (p = 0.010) and increased CD8+T cell subsets, including CD8+naïve T cell (p < 0.001), CD8+central memory T cell (CD8+Tcm, p = 0.003), and CD8+effector memory T cell (CD8+Tem, p = 0.003). As for CD4+T cell, Th1 cells were significantly higher in HBsAg+ FL (p < 0.001). Significant upregulation of individual genes related to CD8+naïve T cells, CD8+Tcm, CD8+Tem, and Th1 cells was also found in HBsAg+ FL (Figure 3G). The detailed genes of T-cell subsets were listed in Additional file 3: Table S8 as described by the xCell method [30]. Therefore, we concluded that HBV-associated FL was characterized by T-cell inflamed phenotype of CD8+T and Th1 cells.
Increased class-switched memory B cell signatures in HBV-associated FL
Since B cells can be efficiently activated by antigens under the help of T cells [31, 32], B-cell activation signaling pathway was enriched in HBsAg+ FL according to GO terms in Figure 3A, and BCR signaling pathway was upregulated according to GSEA (Additional file 2: Figure S3D). Meanwhile, activated upstream regulators in HBsAg+ FL included BCR complex (Figure 3B). When comparing the difference of B-cell subsets by the xCell [30], we observed decreased naïve B cells (p < 0.001) and increased class-switched memory B cells (MBCs) (p = 0.020) in HBsAg+ FL (Figure 4A). Indeed, expression of individual genes related to naïve B cells was downregulated, while expression of individual genes related to class-switched MBCs was upregulated in HBsAg+ FL (Figure 4B). The detailed genes of B-cell subsets were listed in Additional file 3: Table S8 as described by the xCell method [30]. The finding of increased class-switched MBCs in HBsAg+ FL suggested increased ongoing immune response and post-GC B cells upon chronic HBV antigen stimulation [33, 34].
Increased post-GC signatures in HBV-associated FL
Expression of the major transcription factors involved in GC development were subsequently analyzed (Figure 4C, left panel) [35, 36]. Significantly downregulation of genes involved in early GC initiation, late GC initiation, and GC maintenance (POU2AF1, MEFC, IRF8, BCL6, BACH2, and EBF1), while upregulation of genes involved in post-GC (IRF4 and PRDM1) were observed in HBsAg+ FL. Indeed, BCL6 and BACH2 were inhibited upstream regulators in HBsAg+ FL as shown in Figure 3B. Expression of key surface markers on GC B and post-GC B cells were also screened (Figure 4C, right panel) [37-39]. In parallel with the expression of transcription factors, lower expression of surface markers on GC B cells (CD10), but higher expression of surface markers on post-GC B cells (CD44) were observed in HBsAg+ FL. Using immunohistochemistry, protein expression of BCL6, MUM1, and CD10 were confirmed (Figure 4D), with significantly decreased expression of BCL6 (p = 0.005) and CD10 (p = 0.002), and increased expression of MUM1 (p < 0.001) in HBsAg+ FL.
Moreover, cell-of-origin (COO) was examined by transcriptomic sequencing data using the method previously reported in DLBCL [40]. Interestingly, HBsAg+ FL exhibited increased activated B cell like (ABC-like) subtype (p < 0.001, Figure 4E), which was reported to derive from post-GC B cells [36]. Detailed unsupervised clustering heatmap was shown in Additional file 2: Figure S4A. These data indicated that HBV-associated FL tumor cells were featured by a post-GC phenotype.
Increased CD8+T and Th1 cells in PBMCs co-cultured with TNFAIP3-mutant cells
To explore the role of immune modulator mutations in HBV-associated FL, the function of TNFAIP3 mutation was investigated. TNFAIP3-wt and TNFAIP3-mutant (L147fs, R706X) plasmids were transfected into the SC-1 cell line, and the transfection efficiency was confirmed by sequencing of RT-PCR products (Additional file 2: Figure S5A). Enhanced cell proliferation was observed in TNFAIP3-mutant cells (Figure 5A). Upregulated transcriptional factors involved in NF-κB signaling pathway were also confirmed in TNFAIP3-mutant cells by western blot (Additional file 2: Figure S5B).
When co-cultured with PBMCs, a significant increase in CD8+T cells, and higher expression of TIGIT on CD8+T cells were observed in TNFAIP3-mutant cells (Figure 5B). We also found a significant increase in Th1 cells and higher expression of VISTA on Th1 cells in TNFAIP3-mutant cells (Figure 5C).
Evidence of lenalidomide sensitivity in TNFAIP3-mutant cells in vitro and in vivo
TNFAIP3-mutant cells showed increased sensitivity to lenalidomide, with lower half maximal inhibitory concentration [IC50] at 48h and 72h (Figure 5D). When co-cultured with PBMCs, TNFAIP3-mutant cell growth was significantly inhibited by lenalidomide treatment (Figure 5E). Meanwhile, CD8+T cells of the co-cultured system with TNFAIP3-mutant cells were significantly increased, while the expression of TIGIT on CD8+T cells were decreased (Figure 5F). Th1 cells of the co-cultured system with TNFAIP3-mutant cells were significantly increased, while the expression of VISTA on Th1 cells were decreased (Figure 5G).
In vivo, zebrafish xenograft models showed decreased survival when injected with TNFAIP3-mutant cells (Figure 5H, upper panel), and prolonged survival upon lenalidomide treatment (Figure 5H, lower panel).
Clinical evidence of rituximab and lenalidomide sensitivity in HBV-associated FL
In clinical settings, significantly increased CD8+T cells, Th1 cells, as well as expression of TIGIT, and VISTA were observed in the tumor tissues with HBV-associated immune modulator mutations (TNFAIP3/CD70/CXCR4/PIM1/KLF2/FAS/CD58/CD83 mutations, Figure 6A). These findings suggested that lenalidomide may be more effective in HBsAg+ FL. Survival analysis was performed on 253 patients treated with rituximab-based immunochemotherapy or R2 regimen, including 43 HBsAg+ FL and 210 HBsAg- FL patients. The median follow-up time was 18.8 (range 2.0-88.1) months. For rituximab-based immunochemotherapy, HBsAg+ FL presented inferior PFS (p = 0.006, Figure 6B), with the 2-year PFS rate of 82.2%. Of note, for R2 regimen, there was no significant difference in PFS (p = 0.217, Figure 6C) between HBsAg+ FL and HBsAg- FL, with the 2-year PFS rate of 100.0% and 82.5%, respectively. Therefore, rituximab and lenalidomide appeared to abrogate a negative effect of HBV infection on FL outcome upon rituximab-based immunochemotherapy.