Human germinal center transcriptional programs are de-synchronized in B cell lymphoma.

Most adult B cell lymphomas originate from germinal center (GC) B cells, but it is unclear to what extent B cells in overt lymphoma retain the functional dynamics of GC B cells or are blocked at a particular stage of the GC reaction. Here we used integrative single-cell analysis of phenotype, gene expression and variable-region sequence of the immunoglobulin heavy-chain locus to track the characteristic human GC B cell program in follicular lymphoma B cells. By modeling the cyclic continuum of GC B cell transitional states, we identified characteristic patterns of synchronously expressed gene clusters. GC-specific gene-expression synchrony was lost in single lymphoma B cells. However, distinct follicular lymphoma-specific cell states co-existed within single patient biopsies. Our data show that lymphoma B cells are not blocked in a GC B cell state but might adopt new dynamic modes of functional diversity, which opens the possibility of novel definitions of lymphoma identity.

The alternative RelB NF-κB subunit is a novel critical player in diffuse large B-cell lymphoma.

Diffuse large B-cell lymphoma (DLBCL) is the most frequent lymphoid malignancy affecting adults. The NF-κB transcription factor family is activated by 2 main pathways, the canonical and the alternative NF-κB activation pathway, with different functions. The alternative NF-κB pathway leads to activation of the transcriptionally active RelB NF-κB subunit. Alternative NF-κB activation status and its role in DLBCL pathogenesis remain undefined. Here, we reveal a frequent activation of RelB in a large cohort of DLBCL patients and cell lines, independently of their activated B-cell-like or germinal center B-cell-like subtype. RelB activity defines a new subset of patients with DLBCL and a peculiar gene expression profile and mutational pattern. Importantly, RelB activation does not correlate with the MCD genetic subtype, enriched for activated B-cell-like tumors carrying MYD88L265P and CD79B mutations that cooperatively activate canonical NF-κB, thus indicating that current genetic tools to evaluate NF-κB activity in DLBCL do not provide information on the alternative NF-κB activation. Furthermore, the newly defined RelB-positive subgroup of patients with DLBCL exhibits a dismal outcome after immunochemotherapy. Functional studies revealed that RelB confers DLBCL cell resistance to DNA damage-induced apoptosis in response to doxorubicin, a genotoxic agent used in the front-line treatment of DLBCL. We also show that RelB positivity is associated with high expression of cellular inhibitor of apoptosis protein 2 (cIAP2). Altogether, RelB activation can be used to refine the prognostic stratification of DLBCL and may contribute to subvert the therapeutic DNA damage response in a segment of patients with DLBCL.

Oncogenetic landscape of lymphomagenesis in coeliac disease.

OBJECTIVE: Enteropathy-associated T-cell lymphoma (EATL) is a rare but severe complication of coeliac disease (CeD), often preceded by low-grade clonal intraepithelial lymphoproliferation, referred to as type II refractory CeD (RCDII). Knowledge on underlying oncogenic mechanisms remains scarce. Here, we analysed and compared the mutational landscape of RCDII and EATL in order to identify genetic drivers of CeD-associated lymphomagenesis. DESIGN: Pure populations of RCDII-cells derived from intestinal biopsies (n=9) or sorted from blood (n=2) were analysed by whole exome sequencing, comparative genomic hybridisation and RNA sequencing. Biopsies from RCDII (n=50), EATL (n=19), type I refractory CeD (n=7) and uncomplicated CeD (n=18) were analysed by targeted next-generation sequencing. Moreover, functional in vitro studies and drug testing were performed in RCDII-derived cell lines. RESULTS: 80% of RCDII and 90% of EATL displayed somatic gain-of-functions mutations in the JAK1-STAT3 pathway, including a remarkable p.G1097 hotspot mutation in the JAK1 kinase domain in approximately 50% of cases. Other recurrent somatic events were deleterious mutations in nuclear factor kappa-light-chain-enhancer of activated B-cells (NF-κB) regulators TNFAIP3 and TNIP3 and potentially oncogenic mutations in TET2, KMT2D and DDX3X. JAK1 inhibitors, and the proteasome inhibitor bortezomib could block survival and proliferation of malignant RCDII-cell lines. CONCLUSION: Mutations activating the JAK1-STAT3 pathway appear to be the main drivers of CeD-associated lymphomagenesis. In concert with mutations in negative regulators of NF-κB, they may favour the clonal emergence of malignant lymphocytes in the cytokine-rich coeliac intestine. The identified mutations are attractive therapeutic targets to treat RCDII and block progression towards EATL.

Gene alterations in epigenetic modifiers and JAK-STAT signaling are frequent in breast implant-associated ALCL.

The oncogenic events involved in breast implant-associated anaplastic large cell lymphoma (BI-ALCL) remain elusive. To clarify this point, we have characterized the genomic landscape of 34 BI-ALCLs (15 tumor and 19 in situ subtypes) collected from 54 BI-ALCL patients diagnosed through the French Lymphopath network. Whole-exome sequencing (n = 22, with paired tumor/germline DNA) and/or targeted deep sequencing (n = 24) showed recurrent mutations of epigenetic modifiers in 74% of cases, involving notably KMT2C (26%), KMT2D (9%), CHD2 (15%), and CREBBP (15%). KMT2D and KMT2C mutations correlated with a loss of H3K4 mono- and trimethylation by immunohistochemistry. Twenty cases (59%) showed mutations in ≥1 member of the JAK/STAT pathway, including STAT3 (38%), JAK1 (18%), and STAT5B (3%), and in negative regulators, including SOCS3 (6%), SOCS1 (3%), and PTPN1 (3%). These mutations were more frequent in tumor-type samples than in situ samples (P = .038). All BI-ALCLs expressed pSTAT3, regardless of the mutational status of genes in the JAK/STAT pathway. Mutations in the EOMES gene (12%) involved in lymphocyte development, PI3K-AKT/mTOR (6%), and loss-of-function mutations in TP53 (12%) were also identified. Copy-number aberration (CNA) analysis identified recurrent alterations, including gains on chromosomes 2, 9p, 12p, and 21 and losses on 4q, 8p, 15, 16, and 20. Regions of CNA encompassed genes involved in the JAK/STAT pathway and epigenetic regulators. Our results show that the BI-ALCL genomic landscape is characterized by not only JAK/STAT activating mutations but also loss-of-function alterations of epigenetic modifiers.

Adverse outcome in follicular lymphoma is associated with MYC rearrangements but not MYC extra copies.

Follicular lymphomas (FLs) with MYC rearrangements (MYC-R) and extra copies of MYC (MYC-EC) are rare and the prognosis impact is uncertain. We conducted a retrospective study including 321 FL patients, among whom 259 (81%) had no 8q24 alterations and 62 (19%) were assigned to 8qAlt. Forty-five cases were classified as MYC-EC and six as MYC-R. MYC-R patients were significantly older (P = 0·008), had higher follicular lymphoma international prognostic index (FLIPI) stage (P = 0·05) and ß2-microglobulin (ß2m; P = 0·05). Among patients treated with immuno-chemotherapy, four presented a MYC-R and 25 a MYC-EC. Univariate analysis showed the absence of significant difference between MYC-EC and normal MYC (MYC-NL) regarding progression-free survival (PFS; HR1·3; 95% CI [0·4-1·6]) and specific overall survival (SOS; HR 1·6; 95% CI [0·4-5·7]). Those results were compared to data from the PRIMA trial. This confirmed that MYC-EC had no impact on PFS (P = 0·86) or SOS (P = 0·9). Conversely, MYC-R was associated with a trend to inferior outcome regarding PFS (HR : 6·1; 95% CI [2·2-17·1]; P = 0·00026), lymphoma-related death (SOS; HR 13·6; 95% CI [2·9-65]; P = 0·00014) and risk of transformation (transformation-free survival (TFS); HR 82·7; 95% CI [14·8-463·4]; P < 0·0001). In conclusion, MYC-EC has no prognostic impact in FL but MYC-R FL tended to be associated with an increased risk of transformation and poorer outcome.

A gene-expression profiling score for prediction of outcome in patients with follicular lymphoma: a retrospective training and validation analysis in three international cohorts.

BACKGROUND: Patients with follicular lymphoma have heterogeneous outcomes. Predictor models to distinguish, at diagnosis, between patients at high and low risk of progression are needed. The objective of this study was to use gene-expression profiling data to build and validate a predictive model of outcome for patients treated in the rituximab era. METHODS: A training set of fresh-frozen tumour biopsies was prospectively obtained from 160 untreated patients with high-tumour-burden follicular lymphoma enrolled in the phase 3 randomised PRIMA trial, in which rituximab maintenance was evaluated after rituximab plus chemotherapy induction (median follow-up 6·6 years [IQR 6·0-7·0]). RNA of sufficient quality was obtained for 149 of 160 cases, and Affymetrix U133 Plus 2.0 microarrays were used for gene-expression profiling. We did a multivariate Cox regression analysis to identify genes with expression levels associated with progression-free survival independently of maintenance treatment in a subgroup of 134 randomised patients. Expression levels from 95 curated genes were then determined by digital expression profiling (NanoString technology) in 53 formalin-fixed paraffin-embedded samples of the training set to compare the technical reproducibility of expression levels for each gene between technologies. Genes with high correlation (>0·75) were included in an L2-penalised Cox model adjusted on rituximab maintenance to build a predictive score for progression-free survival. The model was validated using NanoString technology to digitally quantify gene expression in 488 formalin-fixed, paraffin-embedded samples from three independent international patient cohorts from the PRIMA trial (n=178; distinct from the training cohort), the University of Iowa/Mayo Clinic Lymphoma SPORE project (n=201), and the Barcelona Hospital Clinic (n=109). All tissue samples consisted of pretreatment diagnostic biopsies and were confirmed as follicular lymphoma grade 1-3a. The patients were all treated with regimens containing rituximab and chemotherapy, possibly followed by either rituximab maintenance or ibritumomab-tiuxetan consolidation. We determined an optimum threshold on the score to predict patients at low risk and high risk of progression. The model, including the multigene score and the threshold, was initially evaluated in the three validation cohorts separately. The sensitivity and specificity of the score for the prediction of the risk of lymphoma progression at 2 years were assessed on the combined validation cohorts. FINDINGS: In the training cohort, the expression levels of 395 genes were associated with a risk of progression. 23 genes reflecting both B-cell biology and tumour microenvironment with correlation coefficients greater than 0·75 between the two technologies and sample types were retained to build a predictive model that identified a population at an increased risk of progression (p<0·0001). In a multivariate Cox model for progression-free survival adjusted on rituximab maintenance treatment and Follicular Lymphoma International Prognostic Index 1 (FLIPI-1) score, this predictor independently predicted progression (adjusted hazard ratio [aHR] of the high-risk group compared with the low-risk group 3·68, 95% CI 2·19-6·17 [p<0·0001]). The 5-year progression-free survival was 26% (95% CI 16-43) in the high-risk group and 73% (64-83) in the low-risk group. The predictor performances were confirmed in each of the individual validation cohorts (aHR comparing high-risk to low-risk groups 2·57 [95% CI 1·65-4·01] in cohort 1; 2·12 [1·32-3·39] in cohort 2; and 2·11 [1·01-4·41] in cohort 3). In the combined validation cohort, the median progression-free survival was 3·1 years (95% CI 2·4-4·8) in the high-risk group and 10·8 years (10·1-not reached) in the low-risk group (p<0·0001). The risk of lymphoma progression at 2 years was 38% (95% CI 29-46) in the high-risk group and 19% (15-24) in the low-risk group. In a multivariate analysis, the score predicted progression-free survival independently of anti-CD20 maintenance treatment and of the FLIPI score (aHR for the combined cohort 2·30, 95% CI 1·72-3·07). INTERPRETATION: We developed and validated a robust 23-gene expression-based predictor of progression-free survival that is applicable to routinely available formalin-fixed, paraffin-embedded tumour biopsies from patients with follicular lymphoma at time of diagnosis. Applying this score could allow individualised therapy for patients according to their risk category. FUNDING: Roche, SIRIC Lyric, LYSARC, National Institutes of Health, the Henry J Predolin Foundation, and the Spanish Plan Nacional de Investigacion.

BCL2 mutations do not confer adverse prognosis in follicular lymphoma patients treated with rituximab.

BCL2 mutations have been suggested to confer an adverse prognosis to follicular lymphoma (FL) patients, but their prognostic value has not been assessed in patients treated with a rituximab-containing regimen. Here we evaluated the prognostic value of BCL2 mutations in a large prospective cohort of 252 patients with FL treated with immunochemotherapy in the PRIMA randomized trial. Using a DNA-targeted sequencing approach, we detected amino acid altering mutations in 135 patients (54%) and showed that these mutations were probably mediated by the over-activation of AICDA (activation-induced cytidine deaminase) in the context of the t(14;18) translocation. The BCL2 variants identified in PRIMA patients affected the BH1, BH2, and BH3 functional motifs at a lower frequency than the N-terminus and flexible loop domain, with mostly conservative aminoacid changes. With a median follow-up of 6.7 years, we did not observe any impact of BCL2 mutations either on overall survival or progression-free survival.

Whole exome sequencing of relapsed/refractory patients expands the repertoire of somatic mutations in diffuse large B-cell lymphoma.

Despite the many efforts already spent to enumerate somatic mutations in diffuse large B-cell lymphoma (DLBCL), previous whole-genome and whole-exome studies conducted on patients of mixed outcomes failed at characterizing the 30% of patients who will relapse or resist current immunochemotherapies. To address this issue, we performed whole-exome sequencing of normal/tumoral DNA pairs in 14 relapsed/refractory (R/R) patients subclassified by full-transcriptome arrays (six activated B-cell like, three germinal center B-cell like, and five primary mediastinal B-cell lymphomas), from the LNH-03 LYSA clinical trial program. Aside from well-known DLBCL features, gene and pathway level recurrence analyses proposed several interesting leads including TBL1XR1 and activating mutations in IRF4 or in the insulin regulation pathway. Sequencing-based copy number analysis defined 23 short recurrently altered regions involving genes such as REL, CDKN2A, HYAL2, and TP53. Moreover, it highlighted mutations in genes such as GNA13, CARD11, MFHAS1, and PCLO as associated with secondary variant allele amplification events. The five primary mediastinal B-cell lymphomas (PMBL), while unexpected in a R/R cohort, showed a significantly higher mutation rate (P = 0.003) and provided many insights on this classical Hodgkin lymphoma related subtype. Novel genes such as XPO1, MFHAS1, and ITPKB were found particularly mutated, along with various cytokine-based signaling pathways. Among these analyses, somatic events in the NF-κB pathway were found preponderant in the three DLBCL subtypes, confirming its major implication in DLBCL aggressiveness and pinpointing several new candidate genes.

Detection of miRNA regulatory effect on triple negative breast cancer transcriptome.

Identifying key microRNAs (miRNAs) contributing to the genesis and development of a particular disease is a focus of many recent studies. We introduce here a rank-based algorithm to detect miRNA regulatory activity in cancer-derived tissue samples which combines measurements of gene and miRNA expression levels and sequence-based target predictions. The method is designed to detect modest but coordinated changes in the expression of sequence-based predicted target genes. We applied our algorithm to a cohort of 129 tumour and healthy breast tissues and showed its effectiveness in identifying functional miRNAs possibly involved in the disease. These observations have been validated using an independent publicly available breast cancer dataset from The Cancer Genome Atlas. We focused on the triple negative breast cancer subtype to highlight potentially relevant miRNAs in this tumour subtype. For those miRNAs identified as potential regulators, we characterize the function of affected target genes by enrichment analysis. In the two independent datasets, the affected targets are not necessarily the same, but display similar enriched categories, including breast cancer related processes like cell substrate adherens junction, regulation of cell migration, nuclear pore complex and integrin pathway. The R script implementing our method together with the datasets used in the study can be downloaded here (http://bioinfo-out.curie.fr/projects/targetrunningsum).