Correction: Whole exome sequencing of high-risk neuroblastoma identifies novel non-synonymous variants.

[This corrects the article DOI: 10.1371/journal.pone.0273280.].

Gene expression analysis revealed downregulation of complement receptor 1 in clonal B cells in cold agglutinin disease.

Cold agglutinin disease (CAD) is a rare B-cell lymphoproliferative disorder of the bone marrow, manifested by autoimmune hemolytic anemia caused by binding of monoclonal IgM autoantibodies to the I antigen. Underlying genetic changes have previously been reported, but their impact on gene expression profile has been unknown. Here, we define differentially expressed genes in CAD B cells. To unravel downstream alteration in cellular pathways, gene expression by RNA sequencing was undertaken. Clonal B-cell samples from 12 CAD patients and IgM-expressing memory B cells from 4 healthy individuals were analyzed. Differential expression analysis and filtering resulted in 93 genes with significant differential expression. Top upregulated genes included SLC4A1, SPTA1, YBX3, TESC, HBD, AHSP, TRAF1, HBA2, RHAG, CA1, SPTB, IL10, UBASH3B, ALAS2, HBA1, CRYM, RGCC, KANK2, and IGHV4-34. They were upregulated at least 8-fold, while complement receptor 1 (CR1/CD35) was downregulated 11-fold in clonal CAD B cells compared to control B cells. Flow cytometry analyses further confirmed reduced CR1 (CD35) protein expression by clonal CAD IgM+ B cells compared to IgM+ memory B cells in controls. CR1 (CD35) is an important negative regulator of B-cell activation and differentiation. Therefore, reduced CR1 (CD35) expression may increase activation, proliferation, and antibody production in CAD-associated clonal B cells.

Whole exome sequencing of high-risk neuroblastoma identifies novel non-synonymous variants.

Neuroblastoma (NBL), one of the main death-causing cancers in children, is known for its remarkable genetic heterogeneity and varied patient outcome spanning from spontaneous regression to widespread disease. Specific copy number variations and single gene rearrangements have been proven to be associated with biological behavior and prognosis; however, there is still an unmet need to enlarge the existing armamentarium of prognostic and therapeutic targets. We performed whole exome sequencing (WES) of samples from 18 primary tumors and six relapse samples originating from 18 NBL patients. Our cohort consists of 16 high-risk, one intermediate, and one very low risk patient. The obtained results confirmed known mutational hotspots in ALK and revealed other non-synonymous variants of NBL-related genes (TP53, DMD, ROS, LMO3, PRUNE2, ERBB3, and PHOX2B) and of genes cardinal for other cancers (KRAS, PIK3CA, and FLT3). Beyond, GOSeq analysis determined genes involved in biological adhesion, neurological cell-cell adhesion, JNK cascade, and immune response of cell surface signaling pathways. We were able to identify novel coding variants present in more than one patient in nine biologically relevant genes for NBL, including TMEM14B, TTN, FLG, RHBG, SHROOM3, UTRN, HLA-DRB1, OR6C68, and XIRP2. Our results may provide novel information about genes and signaling pathways relevant for the pathogenesis and clinical course in high-risk NBL.

Integrating whole genome sequencing, methylation, gene expression, topological associated domain information in regulatory mutation prediction: A study of follicular lymphoma.

A major challenge in human genetics is of the analysis of the interplay between genetic and epigenetic factors in a multifactorial disease like cancer. Here, a novel methodology is proposed to investigate genome-wide regulatory mechanisms in cancer, as studied with the example of follicular Lymphoma (FL). In a first phase, a new machine-learning method is designed to identify Differentially Methylated Regions (DMRs) by computing six attributes. In a second phase, an integrative data analysis method is developed to study regulatory mutations in FL, by considering differential methylation information together with DNA sequence variation, differential gene expression, 3D organization of genome (e.g., topologically associated domains), and enriched biological pathways. Resulting mutation block-gene pairs are further ranked to find out the significant ones. By this approach, BCL2 and BCL6 were identified as top-ranking FL-related genes with several mutation blocks and DMRs acting on their regulatory regions. Two additional genes, CDCA4 and CTSO, were also found in top rank with significant DNA sequence variation and differential methylation in neighboring areas, pointing towards their potential use as biomarkers for FL. This work combines both genomic and epigenomic information to investigate genome-wide gene regulatory mechanisms in cancer and contribute to devising novel treatment strategies.

11q Deletion or ALK Activity Curbs DLG2 Expression to Maintain an Undifferentiated State in Neuroblastoma.

High-risk neuroblastomas typically display an undifferentiated or poorly differentiated morphology. It is therefore vital to understand molecular mechanisms that block the differentiation process. We identify an important role for oncogenic ALK-ERK1/2-SP1 signaling in the maintenance of undifferentiated neural crest-derived progenitors through the repression of DLG2, a candidate tumor suppressor gene in neuroblastoma. DLG2 is expressed in the murine "bridge signature" that represents the transcriptional transition state when neural crest cells or Schwann cell precursors differentiate to chromaffin cells of the adrenal gland. We show that the restoration of DLG2 expression spontaneously drives neuroblastoma cell differentiation, highlighting the importance of DLG2 in this process. These findings are supported by genetic analyses of high-risk 11q deletion neuroblastomas, which identified genetic lesions in the DLG2 gene. Our data also suggest that further exploration of other bridge genes may help elucidate the mechanisms underlying the differentiation of NC-derived progenitors and their contribution to neuroblastomas.

MX 2 is a novel regulator of cell cycle in melanoma cells.

MX2 protein is a dynamin-like GTPase2 that has recently been identified as an interferon-induced restriction factor of HIV-1 and other primate lentiviruses. A single nucleotide polymorphism (SNP), rs45430, in an intron of the MX2 gene, was previously reported as a novel melanoma susceptibility locus in genome-wide association studies. Functionally, however, it is still unclear whether and how MX2 contributes to melanoma susceptibility and tumorigenesis. Here, we show that MX2 is differentially expressed in melanoma tumors and cell lines, with most metastatic cell lines showing lower MX2 expression than primary melanoma cell lines and melanocytes. Furthermore, high expression of MX2 RNA in primary melanoma tumors is associated with better patient survival. Overexpression of MX2 reduces in vivo proliferation partially through inhibition of AKT activation, suggesting that it can act as a tumor suppressor in melanoma. However, we have also identified a subset of melanoma cell lines with high endogenous MX2 expression where downregulation of MX2 leads to reduced proliferation. In these cells, MX2 downregulation interfered with DNA replication and cell cycle processes. Collectively, our data for the first time show that MX2 is functionally involved in the regulation of melanoma proliferation but that its function is context-dependent.

The MCL35 gene expression proliferation assay predicts high-risk MCL patients in a Norwegian cohort of younger patients given intensive first line therapy.

Patients with mantle cell lymphoma (MCL) generally have a dismal prognosis. Intensified induction treatment with rituximab and high dose cytarabine (R_HDAC), and consolidation with high-dose therapy with autologous stem cell support has resulted in 10-year overall survival (OS) higher than 60%. However, the clinical course varies. Diagnostic tools capable of stratifying patients include the MCL International Prognostic Index (MIPI), gene expression-based proliferation signature, Ki-67 proliferation index or tumour cell morphology. Here, we tested the performance of a newly developed Nanostring-based RNA expression-based proliferation assay (MCL35) on formalin-fixed paraffin-embedded tumour tissue from younger patients recruited in or treated according to Nordic MCL protocols compared to the prognosticators listed above. Seventy-four patients were included and the assay performed well in all cases except four, which had inadequate RNA quality. The patients were evenly distributed in the MCL35 low-, intermediate- and high-risk categories. MCL35 low- and intermediate- risk groups had overlapping progression-free survival (PFS), while patients in the high-risk category had significantly inferior PFS. Combining MCL35 with MIPI or the MIPI-C (MIPI with the addition of binary Ki67 score +/-30%) showed a better discrimination than either assessment alone. In conclusion, the MCL35 assay alone or combined with MIPI or MIPI-C scores can identify patients who still have a dismal outcome despite intensified treatment.

Folliculotropic Mycosis Fungoides with Skewed T-cell Receptor CDR3 Motif: Suggestive of Lipid-antigen Selection?

Folliculotropic mycosis fungoides (FMF), a variant of mycosis fungoides (MF) with distinct clinical features, is characterized by infiltration of malignant T cells in hair follicles. This raises the hypothesis that antigens in the hair follicle may contribute to the pathogenesis of FMF. T-cell receptor ß gene (TRB) sequences as well as dendritic cell subsets in patients with FMF (n = 21) and control patients with MF (n = 20) were studied to explore this hypothesis. A recurrent usage of the TRB junctional genes TRBJ2-1 and TRBJ2-7 was found in patients with FMF compared with those with MF. These genes contribute to an amino acid motif in the complementarity-determining region 3 (CDR3) of the T-cell receptor. This motif was previously found in T cells stimulated by lipids bound to CD1 on antigen-presenting cells. Additional immunohistochemical analysis revealed abundant CD1c- and CD1a- expressing dendritic cells in FMF. The combined findings support a role for lipid-antigen stimulation in FMF.

Distinct patterns of B-cell receptor signaling in non-Hodgkin lymphomas identified by single-cell profiling.

Kinases downstream of B-cell antigen receptor (BCR) represent attractive targets for therapy in non-Hodgkin lymphoma (NHL). As clinical responses vary, improved knowledge regarding activation and regulation of BCR signaling in individual patients is needed. Here, using phosphospecific flow cytometry to obtain malignant B-cell signaling profiles from 95 patients representing 4 types of NHL revealed a striking contrast between chronic lymphocytic leukemia (CLL) and mantle cell lymphoma (MCL) tumors. Lymphoma cells from diffuse large B-cell lymphoma patients had high basal phosphorylation levels of most measured signaling nodes, whereas follicular lymphoma cells represented the opposite pattern with no or very low basal levels. MCL showed large interpatient variability in basal levels, and elevated levels for the phosphorylated forms of AKT, extracellular signal-regulated kinase, p38, STAT1, and STAT5 were associated with poor outcome. CLL tumors had elevated basal levels for the phosphorylated forms of BCR-signaling nodes (Src family tyrosine kinase, spleen tyrosine kinase [SYK], phospholipase Cγ), but had low α-BCR-induced signaling. This contrasted MCL tumors, where α-BCR-induced signaling was variable, but significantly potentiated as compared with the other types. Overexpression of CD79B, combined with a gating strategy whereby signaling output was directly quantified per cell as a function of CD79B levels, confirmed a direct relationship between surface CD79B, immunoglobulin M (IgM), and IgM-induced signaling levels. Furthermore, α-BCR-induced signaling strength was variable across patient samples and correlated with BCR subunit CD79B expression, but was inversely correlated with susceptibility to Bruton tyrosine kinase (BTK) and SYK inhibitors in MCL. These individual differences in BCR levels and signaling might relate to differences in therapy responses to BCR-pathway inhibitors.

Comprehensive genome-wide transcription factor analysis reveals that a combination of high affinity and low affinity DNA binding is needed for human gene regulation.

BACKGROUND: High-throughput in vivo protein-DNA interaction experiments are currently widely used in gene regulation studies. Hitherto, comprehensive data analysis remains a challenge and for that reason most computational methods only consider the top few hundred or thousand strongest protein binding sites whereas weak protein binding sites are completely ignored. RESULTS: A new biophysical model of protein-DNA interactions, BayesPI2+, was developed to address the above-mentioned challenges. BayesPI2+ can be run in either a serial computation model or a parallel ensemble learning framework. BayesPI2+ allowed us to analyze all binding sites of the transcription factors, including weak binding that cannot be analyzed by other models. It is evaluated in both synthetic and real in vivo protein-DNA binding experiments. Analysing ESR1 and SPIB in breast carcinoma and activated B cell-like diffuse large B-cell lymphoma cell lines, respectively, revealed that the concerted binding to high and low affinity sites correlates best with gene expression. CONCLUSIONS: BayesPI2+ allows us to analyze transcription factor binding on a larger scale than hitherto achieved. By this analysis, we were able to demonstrate that genes are regulated by concerted binding to high and low affinity binding sites. The program and output results are publicly available at: http://folk.uio.no/junbaiw/BayesPI2Plus.

Lymphoplasmacytic lymphoma and marginal zone lymphoma in the bone marrow: paratrabecular involvement as an important distinguishing feature.

OBJECTIVES: The differential diagnosis between bone marrow involvement by lymphoplasmacytic lymphoma (LPL) and marginal zone lymphoma (MZL) is challenging because histology and immunophenotype of both diseases overlap. We revisited the diagnostic pathology features of both diseases in the bone marrow. METHODS: We studied a series of bone marrow trephine biopsy specimens from 59 patients with Waldenström macroglobulinemia without extramedullary involvement and bone marrow biopsy specimens from 23 patients with well-characterized MZL who also had bone marrow involvement. H&E- and immunoperoxidase-stained sections of bone marrow trephine biopsy specimens as well as flow cytometry and classic cytogenetics performed on aspirations were reviewed. The study was complemented with MYD88 L265P mutation analysis of all samples. RESULTS: The most distinguishing features of LPL with respect to MZL were focal paratrabecular involvement (P < .001), the presence of lymphoplasmacytoid cells (P < .001) and Dutcher bodies (P < .001), increased numbers of mast cells (P < .001), and the MYD88 L265P mutation (P < .001). CONCLUSIONS: LPL can be reliably distinguished from MZL in the bone marrow by using a combination of pathology characteristics. Our findings stress the diagnostic importance of using the combination of the following parameters for a correct LPL diagnosis: paratrabecular infiltration, the presence of lymphoplasmacytoid cells and cells with Dutcher bodies, and an increased number of mast cells in addition to the presence of MYD88 mutation.

A gene panel, including LRP12, is frequently hypermethylated in major types of B-cell lymphoma.

Epigenetic modifications and DNA methylation in particular, have been recognized as important mechanisms to alter gene expression in malignant cells. Here, we identified candidate genes which were upregulated after an epigenetic treatment of B-cell lymphoma cell lines (Burkitt's lymphoma, BL; Follicular lymphoma, FL; Diffuse large B-cell lymphoma, DLBCL activated B-cell like, ABC; and germinal center like, GCB) and simultaneously expressed at low levels in samples from lymphoma patients. Qualitative methylation analysis of 24 candidate genes in cell lines revealed five methylated genes (BMP7, BMPER, CDH1, DUSP4 and LRP12), which were further subjected to quantitative methylation analysis in clinical samples from 59 lymphoma patients (BL, FL, DLBCL ABC and GCB; and primary mediastinal B-cell lymphoma, PMBL). The genes LRP12 and CDH1 showed the highest methylation frequencies (94% and 92%, respectively). BMPER (58%), DUSP4 (32%) and BMP7 (22%), were also frequently methylated in patient samples. Importantly, all gene promoters were unmethylated in various control samples (CD19+ peripheral blood B cells, peripheral blood mononuclear cells and tonsils) as well as in follicular hyperplasia samples, underscoring a high specificity. The combination of LRP12 and CDH1 methylation could successfully discriminate between the vast majority of the lymphoma and control samples, emphasized by receiver operating characteristic analysis with a c-statistic of 0.999. These two genes represent promising epigenetic markers which may be suitable for monitoring of B-cell lymphoma.

Multiple distinct T-cell clones in folliculotropic mycosis fungoides.

Multiple distinct T-cell clones have been demonstrated in a subset of mycosis fungoides (MF), but have so far not been documented in folliculotropic MF, a clinical and histological variant of MF. We analyzed T-cell receptor (TCR) gene rearrangements in 8 patients with folliculotropic MF with multiple biopsies (range, 2-5) taken during the course of disease. Two patients had disease stage IA-IIA, 5 stage IIB-IVA, whereas data were not available for 1 patient. TCR ß and γ gene rearrangements were analyzed according to the BIOMED-2 PCR protocol. Multiple clonal TCR gene rearrangements indicating more than 1 T-cell clone were found in 5 patients. Although the number of patients is small, the finding of multiple distinct T-cell clones in 5 out of 8 patients suggests that chronic T-cell stimulation contributes to the development of folliculotropic MF.

Whole-genome integrative analysis reveals expression signatures predicting transformation in follicular lymphoma.

Transformation of follicular lymphoma (FL) to a more aggressive disease is associated with rapid progression and death. Existing molecular markers for transformation are few and their clinical impact is limited. Here, we report on a whole-genome study of DNA copy numbers and gene expression profiles in serial FL biopsies. We identified 698 genes with high correlation between gene expression and copy number, and the molecular network most enriched for these cis-associated genes. This network includes 14 cis-associated genes directly related to the nuclear factor κB (NF-κB) pathway. For each of these 14 genes, the correlated NF-κB target genes were identified and corresponding expression scores were defined. The scores for 6 of the cis-associated NFκB pathway genes (BTK, IGBP1, IRAK1, ROCK1, TMED7-TICAM2, and TRIM37) were significantly associated with transformation. The results suggest that genes regulating B-cell survival and activation are involved in transformation of FL.