APRIL Drives a Coordinated but Diverse Response as a Foundation for Plasma Cell Longevity.

Ab-secreting cells survive in niche microenvironments, but cellular responses driven by particular niche signals are incompletely defined. The TNF superfamily member a proliferation-inducing ligand (APRIL) can support the maturation of transitory plasmablasts into long-lived plasma cells. In this study, we explore the biological programs established by APRIL in human plasmablasts. Under conditions allowing the maturation of ex vivo- or in vitro-generated plasmablasts, we find that APRIL drives activation of ERK, p38, and JNK, accompanied by a classical NF-κB response and activation of the AKT/FOXO1 pathway. Time-course gene expression data resolve coordinated transcriptional responses propagated via immediate early genes and NF-κB targets and converging onto modules of genes enriched for MYC targets and metabolism/cell growth-related pathways. This response is shared between APRIL and an alternate TNF superfamily member CD40L but is not a feature of alternative niche signals delivered by IFN-α or SDF1. However, APRIL and CD40L responses also diverge. CD40L drives expression of genes related to the activated B cell state whereas APRIL does not. Thus, APRIL establishes a broad foundation for plasma cell longevity with features of cellular refueling while being uncoupled from support of the B cell state.

A System for In Vitro Generation of Mature Murine Plasma Cells Uncovers Differential Blimp-1/Prdm1 Promoter Usage.

Upon encounter with Ag, B cells undergo a sequential process of differentiation to become Ab-secreting plasma cells. Although the key drivers of differentiation have been identified, research has been limited by the lack of in vitro models recapitulating the full process for murine B cells. In this study, we describe methodology using BCR or TLR ligation to obtain plasma cells that are phenotypically mature, have exited cell cycle and express a gene signature concordant with long-lived plasma cells. Dependent on the initial stimuli, the transcriptomes also show variation including the enhanced expression of matrisome components after BCR stimulation, suggestive of unique functional properties for the resultant plasma cells. Moreover, using the new culture conditions we demonstrate that alternative promoter choice regulating the expression of the master transcription factor Blimp-1/Prdm1 can be observed; when the canonical B cell promoter for Prdm1 is deleted, differentiating B cells exhibit flexibility in the choice of promoter, dictated by the initiating stimulus, with preferential maintenance of expression following exposure to TLR ligation. Thus our system provides a readily tractable model for furthering our understanding of plasma cell biology.

Comparative analysis of gene expression platforms for cell-of-origin classification of diffuse large B-cell lymphoma shows high concordance.

Cell-of-origin subclassification of diffuse large B cell lymphoma (DLBCL) into activated B cell-like (ABC), germinal centre B cell-like (GCB) and unclassified (UNC) or type III by gene expression profiling is recommended in the latest update of the World Health Organization's classification of lymphoid neoplasms. There is, however, no accepted gold standard method or dataset for this classification. Here, we compare classification results using gene expression data for 68 formalin-fixed paraffin-embedded DLBCL samples measured on four different gene expression platforms (Illumina wG-DASLTM arrays, Affymetrix PrimeView arrays, Illumina TrueSeq RNA sequencing and the HTG EdgeSeq DLBCL Cell of Origin Assay EU using an established platform agnostic classification algorithm (DAC) and the classifier native to the HTG platform, which is CE marked for in vitro diagnostic use (CE-IVD). Classification methods and platforms show a high level of concordance, with agreement in at least 80% of cases and rising to much higher levels for classifications of high confidence. Our results demonstrate that cell-of-origin classification by gene expression profiling on different platforms is robust, and that the use of the confidence value alongside the classification result is important in clinical applications.

Growth Factor-like Gene Regulation Is Separable from Survival and Maturation in Antibody-Secreting Cells.

Recurrent mutational activation of the MAP kinase pathway in plasma cell myeloma implicates growth factor-like signaling responses in the biology of Ab-secreting cells (ASCs). Physiological ASCs survive in niche microenvironments, but how niche signals are propagated and integrated is poorly understood. In this study, we dissect such a response in human ASCs using an in vitro model. Applying time course expression data and parsimonious gene correlation network analysis (PGCNA), a new approach established by our group, we map expression changes that occur during the maturation of proliferating plasmablast to quiescent plasma cell under survival conditions including the potential niche signal TGF-ß3. This analysis demonstrates a convergent pattern of differentiation, linking unfolded protein response/endoplasmic reticulum stress to secretory optimization, coordinated with cell cycle exit. TGF-ß3 supports ASC survival while having a limited effect on gene expression including upregulation of CXCR4. This is associated with a significant shift in response to SDF1 in ASCs with amplified ERK1/2 activation, growth factor-like immediate early gene regulation and EGR1 protein expression. Similarly, ASCs responding to survival conditions initially induce partially overlapping sets of immediate early genes without sustaining the response. Thus, in human ASCs growth factor-like gene regulation is transiently imposed by niche signals but is not sustained during subsequent survival and maturation.

Gene-expression profiling of bortezomib added to standard chemoimmunotherapy for diffuse large B-cell lymphoma (REMoDL-B): an open-label, randomised, phase 3 trial.

BACKGROUND: Biologically distinct subtypes of diffuse large B-cell lymphoma can be identified using gene-expression analysis to determine their cell of origin, corresponding to germinal centre or activated B cell. We aimed to investigate whether adding bortezomib to standard therapy could improve outcomes in patients with these subtypes. METHODS: In a randomised evaluation of molecular guided therapy for diffuse large B-cell lymphoma with bortezomib (REMoDL-B), an open-label, adaptive, randomised controlled, phase 3 superiority trial, participants were recruited from 107 cancer centres in the UK (n=94) and Switzerland (n=13). Eligible patients had previously untreated, histologically confirmed diffuse large B-cell lymphoma with sufficient diagnostic material from initial biopsies for gene-expression profiling and pathology review; were aged 18 years or older; had ECOG performance status of 2 or less; had bulky stage I or stage II-IV disease requiring full-course chemotherapy; had measurable disease; and had cardiac, lung, renal, and liver function sufficient to tolerate chemotherapy. Patients initially received one 21-day cycle of standard rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisolone (R-CHOP; rituximab 375 mg/m2, cyclophosphamide 750 mg/m2, doxorubicin 50 mg/m2, and vincristine 1·4 mg/m2 [to a maximum of 2 mg total dose] intravenously on day 1 of the cycle, and prednisolone 100 mg orally once daily on days 1-5). During this time, we did gene-expression profiling using whole genome cDNA-mediated annealing, selection, extension, and ligation assay of tissue from routine diagnostic biopsy samples to determine the cell-of-origin subtype of each participant (germinal centre B cell, activated B cell, or unclassified). Patients were then centrally randomly assigned (1:1) via a web-based system, with block randomisation stratified by international prognostic index score and cell-of-origin subtype, to continue R-CHOP alone (R-CHOP group; control), or with bortezomib (RB-CHOP group; experimental; 1·3 mg/m2 intravenously or 1·6 mg/m2 subcutaneously) on days 1 and 8 for cycles two to six. If RNA extracted from the diagnostic tissues was of insufficient quality or quantity, participants were given R-CHOP as per the control group. The primary endpoint was 30-month progression-free survival, for the germinal centre and activated B-cell population. The primary analysis was on the modified intention-to-treat population of activated and germinal centre B-cell population. Safety was assessed in all participants who were given at least one dose of study drug. We report the progression-free survival and safety outcomes for patients in the follow-up phase after the required number of events occurred. This study was registered at ClinicalTrials.gov, number NCT01324596, and recruitment and treatment has completed for all participants, with long-term follow-up ongoing. FINDINGS: Between June 2, 2011, and June 10, 2015, 1128 eligible patients were registered, of whom 918 (81%) were randomly assigned to receive treatment (n=459 to R-CHOP, n=459 to RB-CHOP), comprising 244 (26·6%) with activated B-cell disease, 475 (51·7%) with germinal centre B cell disease, and 199 (21·7%) with unclassified disease. At a median follow-up of 29·7 months (95% CI 29·0-32·0), we saw no evidence for a difference in progression-free survival in the combined germinal centre and activated B-cell population between R-CHOP and RB-CHOP (30-month progression-free survival 70·1%, 95% CI 65·0-74·7 vs 74·3%, 69·3-78·7; hazard ratio 0·86, 95% CI 0·65-1·13; p=0·28). The most common grade 3 or worse adverse event was haematological toxicity, reported in 178 (39·8%) of 447 patients given R-CHOP and 187 (42·1%) of 444 given RB-CHOP. However, RB-CHOP was not associated with increased haematological toxicity and 398 [87·1%] of 459 participants assigned to receive RB-CHOP completed six cycles of treatment. Grade 3 or worse neuropathy occurred in 17 (3·8%) patients given RB-CHOP versus eight (1·8%) given R-CHOP. Serious adverse events occurred in 190 (42·5%) patients given R-CHOP, including five treatment-related deaths, and 223 (50·2%) given RB-CHOP, including four treatment-related deaths. INTERPRETATION: This is the first large-scale study in diffuse large B-cell lymphoma to use real-time molecular characterisation for prospective stratification, randomisation, and subsequent analysis of biologically distinct subgroups of patients. The addition of bortezomib did not improve progression-free survival. FUNDING: Janssen-Cilag, Bloodwise, and Cancer Research UK.

Biallelic interferon regulatory factor 8 mutation: A complex immunodeficiency syndrome with dendritic cell deficiency, monocytopenia, and immune dysregulation.

BACKGROUND: The homozygous K108E mutation of interferon regulatory factor 8 (IRF8) is reported to cause dendritic cell (DC) and monocyte deficiency. However, more widespread immune dysfunction is predicted from the multiple roles ascribed to IRF8 in immune cell development and function. OBJECTIVE: We sought to describe the effect on hematopoiesis and immunity of the compound heterozygous R83C/R291Q mutation of IRF8, which is present in a patient with recurrent viral infection, granuloproliferation, and intracerebral calcification. METHODS: Variant IRF8 alleles were identified by means of exome sequencing, and their function was tested by using reporter assays. The cellular phenotype was studied in detail by using flow cytometry, functional immunologic assay transcriptional profiling, and antigen receptor profiling. RESULTS: Both mutations affected conserved residues, and R291Q is orthologous to R294, which is mutated in the BXH2 IRF8-deficient mouse. R83C showed reduced nuclear translocation, and neither mutant was able to regulate the Ets/IRF composite element or interferon-stimulated response element, whereas R291Q retained BATF/JUN interactions. DC deficiency and monocytopenia were observed in blood, dermis, and lung lavage fluid. Granulocytes were consistently increased, dysplastic, and hypofunctional. Natural killer cell development and maturation were arrested. TH1, TH17, and CD8+ memory T-cell differentiation was significantly reduced, and T cells did not express CXCR3. B-cell development was impaired, with fewer memory cells, reduced class-switching, and lower frequency and complexity of somatic hypermutation. Cell-specific gene expression was widely disturbed in interferon- and IRF8-regulated transcripts. CONCLUSIONS: This analysis defines the clinical features of human biallelic IRF8 deficiency, revealing a complex immunodeficiency syndrome caused by DC and monocyte deficiency combined with widespread immune dysregulation.

Network Analysis Identifies Proinflammatory Plasma Cell Polarization for Secretion of ISG15 in Human Autoimmunity.

Plasma cells (PCs) as effectors of humoral immunity produce Igs to match pathogenic insult. Emerging data suggest more diverse roles exist for PCs as regulators of immune and inflammatory responses via secretion of factors other than Igs. The extent to which such responses are preprogrammed in B-lineage cells or can be induced in PCs by the microenvironment is unknown. In this study, we dissect the impact of IFNs on the regulatory networks of human PCs. We show that core PC programs are unaffected, whereas PCs respond to IFNs with distinctive transcriptional responses. The IFN-stimulated gene 15 (ISG15) system emerges as a major transcriptional output induced in a sustained fashion by IFN-α in PCs and linked both to intracellular conjugation and ISG15 secretion. This leads to the identification of ISG15-secreting plasmablasts/PCs in patients with active systemic lupus erythematosus. Thus, ISG15-secreting PCs represent a distinct proinflammatory PC subset providing an Ig-independent mechanism of PC action in human autoimmunity.

SPIB and BATF provide alternate determinants of IRF4 occupancy in diffuse large B-cell lymphoma linked to disease heterogeneity.

Interferon regulatory factor 4 (IRF4) is central to the transcriptional network of activated B-cell-like diffuse large B-cell lymphoma (ABC-DLBCL), an aggressive lymphoma subgroup defined by gene expression profiling. Since cofactor association modifies transcriptional regulatory input by IRF4, we assessed genome occupancy by IRF4 and endogenous cofactors in ABC-DLBCL cell lines. IRF4 partners with SPIB, PU.1 and BATF genome-wide, but SPIB provides the dominant IRF4 partner in this context. Upon SPIB knockdown IRF4 occupancy is depleted and neither PU.1 nor BATF acutely compensates. Integration with ENCODE data from lymphoblastoid cell line GM12878, demonstrates that IRF4 adopts either SPIB- or BATF-centric genome-wide distributions in related states of post-germinal centre B-cell transformation. In primary DLBCL high-SPIB and low-BATF or the reciprocal low-SPIB and high-BATF mRNA expression links to differential gene expression profiles across nine data sets, identifying distinct associations with SPIB occupancy, signatures of B-cell differentiation stage and potential pathogenetic mechanisms. In a population-based patient cohort, SPIBhigh/BATFlow-ABC-DLBCL is enriched for mutation of MYD88, and SPIBhigh/BATFlow-ABC-DLBCL with MYD88-L265P mutation identifies a small subgroup of patients among this otherwise aggressive disease subgroup with distinct favourable outcome. We conclude that differential expression of IRF4 cofactors SPIB and BATF identifies biologically and clinically significant heterogeneity among ABC-DLBCL.

In vitro generation of long-lived human plasma cells.

Plasma cells (PCs), the terminal effectors of humoral immunity, are short-lived unless supported by niche environments in which they may persist for years. No model system has linked B cell activation with niche function to allow the in vitro generation of long-lived PCs. Thus, the full trajectory of B cell terminal differentiation has yet to be investigated in vitro. In this article, we describe a robust model for the generation of polyclonal long-lived human PCs from peripheral blood B cells. After a proliferative plasmablast phase, PCs persist in the absence of cell division, with viability limited only by elective culture termination. Conservative predictions for PC life expectancy are 300 d, but with the potential for significantly longer life spans for some cells. These long-lived PCs are preferentially derived from memory B cells, and acquire a CD138(high) phenotype analogous to that of human bone marrow PCs. Analysis of gene expression across the system defines clusters of genes with related dynamics and linked functional characteristics. Importantly, genes in these differentiation clusters demonstrate a similar overall pattern of expression for in vitro and ex vivo PCs. In vitro PCs are fully reprogrammed to a secretory state and are adapted to their secretory load, maintaining IgG secretion of 120 pg/cell/day in the absence of XBP1 mRNA splicing. By establishing a set of conditions sufficient to allow the development and persistence of mature human PCs in vitro, to our knowledge, we provide the first platform with which to sequentially explore and manipulate each stage of human PC differentiation.

IDHwt glioblastomas can be stratified by their transcriptional response to standard treatment, with implications for targeted therapy.

BACKGROUND: Glioblastoma (GBM) brain tumors lacking IDH1 mutations (IDHwt) have the worst prognosis of all brain neoplasms. Patients receive surgery and chemoradiotherapy but tumors almost always fatally recur. RESULTS: Using RNA sequencing data from 107 pairs of pre- and post-standard treatment locally recurrent IDHwt GBM tumors, we identify two responder subtypes based on longitudinal changes in gene expression. In two thirds of patients, a specific subset of genes is upregulated from primary to recurrence (Up responders), and in one third, the same genes are downregulated (Down responders), specifically in neoplastic cells. Characterization of the responder subtypes indicates subtype-specific adaptive treatment resistance mechanisms that are associated with distinct changes in the tumor microenvironment. In Up responders, recurrent tumors are enriched in quiescent proneural GBM stem cells and differentiated neoplastic cells, with increased interaction with the surrounding normal brain and neurotransmitter signaling, whereas Down responders commonly undergo mesenchymal transition. ChIP-sequencing data from longitudinal GBM tumors suggests that the observed transcriptional reprogramming could be driven by Polycomb-based chromatin remodeling rather than DNA methylation. CONCLUSIONS: We show that the responder subtype is cancer-cell intrinsic, recapitulated in in vitro GBM cell models, and influenced by the presence of the tumor microenvironment. Stratifying GBM tumors by responder subtype may lead to more effective treatment.

Spontaneous EBV-Reactivation during B Cell Differentiation as a Model for Polymorphic EBV-Driven Lymphoproliferation.

Epstein-Barr virus (EBV)-driven B cell neoplasms arise from the reactivation of latently infected B cells. In a subset of patients, EBV was seen to drive a polymorphous lymphoproliferative disorder (LPD) in which B cell differentiation was retained. In this work, spontaneous EBV reactivation following B cell mitogen stimulation was shown to provide a potential model of polymorphic EBV-driven LPD. Here, we developed an in vitro model of plasma cell (PC) differentiation from peripheral blood memory B cells. To assess the frequency and phenotypes of EBV-associated populations derived during differentiation, we analysed eight differentiations during the PC stage with a targeted single-cell gene expression panel. We identified subpopulations of EBV-gene expressing cells with PC and/or B cell expression features in differentiations from all tested donors. EBV-associated cells varied in frequency, ranging from 3-28% of cells. Most EBV-associated cells expressed PC genes such as XBP1 or MZB1, and in all samples these included a quiescent PC fraction that lacked cell a cycle gene expression. With increasing EBV-associated cells, populations with B cell features became prominent, co-expressing a germinal centre (GC) and activating B cell gene patterns. The presence of highly proliferative EBV-associated cells was linked to retained MS4A1/CD20 expression and IGHM and IGHD co-expression, while IGHM class-switched cells were enriched in quiescent PC fractions. Thus, patterns of gene expression in primary EBV reactivation were shown to include features related to GC B cells, which was also observed in EBV-transformed lymphoblastoid cell lines. This suggests a particular association between spontaneously developing EBV-expansions and IgM+ IgD+ non-switched B cells.

Whole genome expression profiling based on paraffin embedded tissue can be used to classify diffuse large B-cell lymphoma and predict clinical outcome.

This study tested the validity of whole-genome expression profiling (GEP) using RNA from formalin-fixed, paraffin-embedded (FFPE) tissue to sub-classify Diffuse Large B-cell Lymphoma (DLBCL), in a population based cohort of 172 patients. GEP was performed using Illumina Whole Genome cDNA-mediated Annealing, Selection, extension & Ligation, and tumours were classified into germinal centre (GCB), activated B-cell (ABC) and Type-III subtypes. The method was highly reproducible and reliably classified cell lines of known phenotype. GCB and ABC subtypes were each characterized by unique gene expression signatures consistent with previously published data. A significant relationship between subtype and survival was observed, with ABC having the worst clinical outcome and in a multivariate survival model only age and GEP class remained significant. This effect was not seen when tumours were classified by immunohistochemistry. There was a significant association between age and subtype (mean ages ABC - 72·8 years, GC - 68·4 years, Type-III - 64·5 years). Older patients with ABC subtype were also over-represented in patients who died soon after diagnosis. The relationship between prognosis and subtype improved when only patients assigned to the three categories with the highest level of confidence were analysed. This study demonstrates that GEP-based classification of DLBCL can be applied to RNA extracted from routine FFPE samples and has potential for use in stratified medicine trials and clinical practice.

An extended set of PRDM1/BLIMP1 target genes links binding motif type to dynamic repression.

The transcriptional repressor B lymphocyte-induced maturation protein-1 (BLIMP1) regulates gene expression and cell fate. The DNA motif bound by BLIMP1 in vitro overlaps with that of interferon regulatory factors (IRFs), which respond to inflammatory/immune signals. At such sites, BLIMP1 and IRFs can antagonistically regulate promoter activity. In vitro motif selection predicts that only a subset of BLIMP1 or IRF sites is subject to antagonistic regulation, but the extent to which antagonism occurs is unknown, since an unbiased assessment of BLIMP1 occupancy in vivo is lacking. To address this, we identified an extended set of promoters occupied by BLIMP1. Motif discovery and enrichment analysis demonstrate that multiple motif variants are required to capture BLIMP1 binding specificity. These are differentially associated with CpG content, leading to the observation that BLIMP1 DNA-binding is methylation sensitive. In occupied promoters, only a subset of BLIMP1 motifs overlap with IRF motifs. Conversely, a distinct subset of IRF motifs is not enriched amongst occupied promoters. Genes linked to occupied promoters containing overlapping BLIMP1/IRF motifs (e.g. AIM2, SP110, BTN3A3) are shown to constitute a dynamic target set which is preferentially activated by BLIMP1 knock-down. These data confirm and extend the competitive model of BLIMP1 and IRF interaction.

GO-At: in silico prediction of gene function in Arabidopsis thaliana by combining heterogeneous data.

Despite recent advances, accurate gene function prediction remains an elusive goal, with very few methods directly applicable to the plant Arabidopsis thaliana. In this study, we present GO-At (gene ontology prediction in A. thaliana), a method that combines five data types (co-expression, sequence, phylogenetic profile, interaction and gene neighbourhood) to predict gene function in Arabidopsis. Using a simple, yet powerful two-step approach, GO-At first generates a list of genes ranked in descending order of probability of functional association with the query gene. Next, a prediction score is automatically assigned to each function in this list based on the assumption that functions appearing most frequently at the top of the list are most likely to represent the function of the query gene. In this way, the second step provides an effective alternative to simply taking the 'best hit' from the first list, and achieves success rates of up to 79%. GO-At is applicable across all three GO categories: molecular function, biological process and cellular component, and can assign functions at multiple levels of annotation detail. Furthermore, we demonstrate GO-At's ability to predict functions of uncharacterized genes by identifying ten putative golgins/Golgi-associated proteins amongst 8219 genes of previously unknown cellular component and present independent evidence to support our predictions. A web-based implementation of GO-At (http://www.bioinformatics.leeds.ac.uk/goat) is available, providing a unique resource for plant researchers to make predictions for uncharacterized genes and predict novel functions in Arabidopsis.

A dichotomy of gene regulatory associations during the activated B-cell to plasmablast transition.

The activated B-cell (ABC) to plasmablast transition encompasses the cusp of antibody-secreting cell (ASC) differentiation. We explore this transition with integrated analysis in human cells, focusing on changes that follow removal from CD40-mediated signals. Within hours of input signal loss, cell growth programs shift toward enhanced proliferation, accompanied by ER-stress response, and up-regulation of ASC features. Clustering of genomic occupancy for IRF4, BLIMP1, XBP1, and CTCF with histone marks identifies a dichotomy: XBP1 and IRF4 link to induced but not repressed gene modules in plasmablasts, whereas BLIMP1 links to modules of ABC genes that are repressed, but not to activated genes. Between ABC and plasmablast states, IRF4 shifts away from AP1/IRF composite elements while maintaining occupancy at IRF and ETS/IRF elements. This parallels the loss of BATF expression, which is identified as a potential BLIMP1 target. In plasmablasts, IRF4 acquires an association with CTCF, a feature maintained in plasma cell myeloma lines. Thus, shifting occupancy links IRF4 to both ABC and ASC gene expression, whereas BLIMP1 occupancy links to repression of the activation state.

TLR-mediated activation of Waldenström macroglobulinemia B cells reveals an uncoupling from plasma cell differentiation.

Waldenström macroglobulinemia (WM) is a rare malignancy in which clonal B cells infiltrate the bone marrow and give rise to a smaller compartment of neoplastic plasma cells that secrete monoclonal immunoglobulin M paraprotein. Recent studies into underlying mutations in WM have enabled a much greater insight into the pathogenesis of this lymphoma. However, there is considerably less characterization of the way in which WM B cells differentiate and how they respond to immune stimuli. In this study, we assess WM B-cell differentiation using an established in vitro model system. Using T-cell-dependent conditions, we obtained CD138+ plasma cells from WM samples with a frequency similar to experiments performed with B cells from normal donors. Unexpectedly, a proportion of the WM B cells failed to upregulate CD38, a surface marker that is normally associated with plasmablast transition and maintained as the cells proceed with differentiation. In normal B cells, concomitant Toll-like receptor 7 (TLR7) activation and B-cell receptor cross-linking drives proliferation, followed by differentiation at similar efficiency to CD40-mediated stimulation. In contrast, we found that, upon stimulation with TLR7 agonist R848, WM B cells failed to execute the appropriate changes in transcriptional regulators, identifying an uncoupling of TLR signaling from the plasma cell differentiation program. Provision of CD40L was sufficient to overcome this defect. Thus, the limited clonotypic WM plasma cell differentiation observed in vivo may result from a strict requirement for integrated activation.

Distinct genetic changes reveal evolutionary history and heterogeneous molecular grade of DLBCL with MYC/BCL2 double-hit.

Using a Burkitt lymphoma-like gene expression signature, we recently defined a high-risk molecular high-grade (MHG) group mainly within germinal centre B-cell like diffuse large B-cell lymphomas (GCB-DLBCL), which was enriched for MYC/BCL2 double-hit (MYC/BCL2-DH). The genetic basis underlying MHG-DLBCL and their aggressive clinical behaviour remain unknown. We investigated 697 cases of DLBCL, particularly those with MYC/BCL2-DH (n = 62) by targeted sequencing and gene expression profiling. We showed that DLBCL with MYC/BCL2-DH, and those with BCL2 translocation, harbour the characteristic mutation signatures that are associated with follicular lymphoma and its high-grade transformation. We identified frequent MYC hotspot mutations that affect the phosphorylation site (T58) and its adjacent amino acids, which are important for MYC protein degradation. These MYC mutations were seen in a subset of cases with MYC translocation, but predominantly in those of MHG. The mutations were more frequent in double-hit lymphomas with IG as the MYC translocation partner, and were associated with higher MYC protein expression and poor patient survival. DLBCL with MYC/BCL2-DH and those with BCL2 translocation alone are most likely derived from follicular lymphoma or its precursor lesion, and acquisition of MYC pathogenic mutations may augment MYC function, resulting in aggressive clinical behaviour.

Parsimonious Gene Correlation Network Analysis (PGCNA): a tool to define modular gene co-expression for refined molecular stratification in cancer.

Cancers converge onto shared patterns that arise from constraints placed by the biology of the originating cell lineage and microenvironment on programs driven by oncogenic events. Here we define consistent expression modules reflecting this structure in colon and breast cancer by exploiting expression data resources and a new computationally efficient approach that we validate against other comparable methods. This approach, Parsimonious Gene Correlation Network Analysis (PGCNA), allows comparison of network structures between these cancer types identifying shared modules of gene co-expression reflecting: cancer hallmarks, functional and structural gene batteries, copy number variation and biology of originating lineage. These networks along with the mapping of outcome data at gene and module level provide an interactive resource that generates context for relationships between genes within and between such modules. Assigning module expression values (MEVs) provides a tool to summarize network level gene expression in individual cases illustrating potential utility in classification and allowing analysis of linkage between module expression and mutational state. Exploiting TCGA data thus defines both recurrent patterns of association between module expression and mutation at data-set level, and exemplifies the polarization of mutation patterns with the leading edge of module expression at individual case level. We illustrate the scalable nature of the approach within immune response related modules, which in the context of breast cancer demonstrates the selective association of immune subsets, in particular mast cells, with the underlying mutational pattern. Together our analyses provide evidence for a generalizable framework to enhance molecular stratification in cancer.

Molecular High-Grade B-Cell Lymphoma: Defining a Poor-Risk Group That Requires Different Approaches to Therapy.

PURPOSE: Biologic heterogeneity is a feature of diffuse large B-cell lymphoma (DLBCL), and the existence of a subgroup with poor prognosis and phenotypic proximity to Burkitt lymphoma is well known. Conventional cytogenetics identifies some patients with rearrangements of MYC and BCL2 and/or BCL6 (double-hit lymphomas) who are increasingly treated with more intensive chemotherapy, but a more biologically coherent and clinically useful definition of this group is required. PATIENTS AND METHODS: We defined a molecular high-grade (MHG) group by applying a gene expression-based classifier to 928 patients with DLBCL from a clinical trial that investigated the addition of bortezomib to standard rituximab plus cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP) therapy. The prognostic significance of MHG was compared with existing biomarkers. We performed targeted sequencing of 70 genes in 400 patients and explored molecular pathology using gene expression signature databases. Findings were validated in an independent data set. RESULTS: The MHG group comprised 83 patients (9%), with 75 in the cell-of-origin germinal center B-cell-like group. MYC rearranged and double-hit groups were strongly over-represented in MHG but comprised only one half of the total. Gene expression analysis revealed a proliferative phenotype with a relationship to centroblasts. Progression-free survival rate at 36 months after R-CHOP in the MHG group was 37% (95% CI, 24% to 55%) compared with 72% (95% CI, 68% to 77%) for others, and an analysis of treatment effects suggested a possible positive effect of bortezomib. Double-hit lymphomas lacking the MHG signature showed no evidence of worse outcome than other germinal center B-cell-like cases. CONCLUSION: MHG defines a biologically coherent high-grade B-cell lymphoma group with distinct molecular features and clinical outcomes that effectively doubles the size of the poor-prognosis, double-hit group. Patients with MHG may benefit from intensified chemotherapy or novel targeted therapies.

Deleterious SNP prediction: be mindful of your training data!

MOTIVATION: To predict which of the vast number of human single nucleotide polymorphisms (SNPs) are deleterious to gene function or likely to be disease associated is an important problem, and many methods have been reported in the literature. All methods require data sets of mutations classified as 'deleterious' or 'neutral' for training and/or validation. While different workers have used different data sets there has been no study of which is best. Here, the three most commonly used data sets are analysed. We examine their contents and relate this to classifiers, with the aims of revealing the strengths and pitfalls of each data set, and recommending a best approach for future studies. RESULTS: The data sets examined are shown to be substantially different in content, particularly with regard to amino acid substitutions, reflecting the different ways in which they are derived. This leads to differences in classifiers and reveals some serious pitfalls of some data sets, making them less than ideal for non-synonymous SNP prediction. AVAILABILITY: Software is available on request from the authors.