Analysis of immunoglobulin/T-cell receptor repertoires by high-throughput RNA sequencing reveals a continuous dynamic of positive clonal selection in follicular lymphoma

Follicular lymphoma (FL) course is highly variable, making its clinical management challenging. In this incurable and recurring pathology, the interval between relapses tends to decrease while aggressiveness increases, sometimes resulting in the transformation to higher-grade lymphoma. These evolutions are particularly difficult to anticipate, resulting from complex clonal evolutions where multiple subclones compete and thrive due to their capacity to proliferate and resist therapies. Here, to apprehend further these processes, we used a high-throughput RNA sequencing approach to address simultaneously the B-cell immunoglobulin repertoires and T-cell immunoglobulin repertoires repertoires of lymphoma cells and their lymphoid microenvironment in a large cohort of 131 FL1/2-3A patients. Our data confirm the existence of a high degree of intra-clonal heterogeneity in this pathology, resulting from ongoing somatic hyper-mutation and class switch recombination. Through the evaluation of the Simpson ecological-diversity index, we show that the contribution of the cancerous cells increases during the course of the disease to the detriment of the reactive compartment, a phenomenon accompanied by a concomitant decrease in the diversity of the tumoral population. Clonal evolution in FL thus contrasts with many tumors, where clonal heterogeneity steadily increases over time and participates in treatment evasion. In this pathology, the selection of lymphoma subclones with proliferative advantages progressively outweighs clonal diversification, ultimately leading in extreme cases to transformation to high-grade lymphoma resulting from the rapid emergence of homogeneous subpopulations.

IGH 3'RR recombination uncovers a non-germinal center imprint and c-MYC-dependent IGH rearrangement in unmutated chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) is an incurable indolent non-Hodgkin lymphoma characterized by tumor B cells that weakly express a B-cell receptor. The mutational status of the variable region (IGHV) within the immunoglobulin heavy chain (IGH) locus is an important prognosis indicator and raises the question of the CLL cell of origin. Mutated IGHV gene CLL are genetically imprinted by activation-induced cytidine deaminase (AID). AID is also required for IGH rearrangements: class switch recombination and recombination between switch Mu (Sμ) and the 3' regulatory region (3'RR) (Sμ-3'RRrec). The great majority of CLL B cells being unswitched led us to examine IGH rearrangement blockade in CLL. Our results separated CLL into two groups on the basis of Sμ-3'RRrec counts per sample: Sμ-3'RRrecHigh cases (mostly unmutated CLL) and Sμ-3'RRrecLow cases (mostly mutated CLL), but not based on the class switch recombination junction counts. Sμ-3'RRrec appeared to be ongoing in Sμ-3'RRrecHigh CLL cells and comparison of Sμ-3'RRrec junction structural features pointed to different B-cell origins for both groups. In accordance with IGHV mutational status and PIM1 mutation rate, Sμ-3'RRrecHigh CLL harbor a non-germinal center experienced B-cell imprint while Sμ-3'RRrecLow CLL are from AID-experienced B cells from a secondary lymphoid organ. In addition to the proposals already made concerning the CLL cell of origin, our study highlights that analysis of IGH recombinatory activity can identify CLL cases from different origins. Finally, on-going Sμ-3'RRrec in Sμ-3'RRrecHigh cells appeared to presumably be the consequence of high c-MYC expression, as c-MYC overexpression potentiated IGH rearrangements and Sμ-3'RRrec, even in the absence of AID for the latter.

Diagnosis of chronic B-cell lymphoproliferative disease in peripheral blood = how machine learning may help to the interpretation of flow cytometry data

Flow cytometry (FCM) has become a method of choice for immunologic characterization of chronic lymphoproliferative disease (CLPD). To reduce the potential subjectivities of FCM data interpretation, we developed a machine learning random forest algorithm (RF) allowing unsupervised analysis. This assay relies on 16 parameters obtained from our FCM screening panel, routinely used in the exploration of peripheral blood (PB) samples (mean fluorescence intensity values (MFI) of CD19, CD45, CD5, CD20, CD200, CD23, HLA-DR, CD10 in CD19-gated B cells, ratio of kappa/Lambda, and different ratios of MFI B-cells/T-cells [CD20, CD200, CD23]). The RF algorithm was trained and validated on a large cohort of more than 300 annotated different CLPD cases (chronic B-cell leukemia, mantle cell lymphoma, marginal zone lymphoma, follicular lymphoma, splenic red pulp lymphoma, hairy cell leukemia) and non-tumoral selected from PB samples. The RF algorithm was able to differentiate tumoral from non-tumoral B-cells in all cases and to propose a correct CLPD classification in more than 90% of cases. In conclusion the RF algorithm could be proposed as an interesting help to FCM data interpretation allowing a first B-cells CLPD diagnostic hypothesis and/or to guide the management of complementary analysis (additional immunologic markers and genetic).

High PDL1/PDL2 gene expression correlates with worse outcome in primary mediastinal large B-cell lymphoma

Primary mediastinal B-cell lymphoma (PMBL) is an uncommon entity of aggressive B-cell lymphoma with an unusually good prognosis, except for 10-15% of chemotherapy-refractory cases. To identify earlier these higher risk patients, we performed molecular characterization of a retrospective multicenter cohort of patients treated with firstline immunochemotherapy. The traits of the patients with gene-expression profiling data (n = 120) were as follows: median age of 34 years (range, 18-67 years); female sex, 58.3%; elevated lactate dehydrogenase, 82.5%; Eastern Cooperative Oncology Group performance status score of 0 to 1, 85.7%; Ann Arbor stage I/II, 55%; International Prognostic Index score of 1 to 2, 64.4%; and median metabolic tumor volume, 290.4 cm3 (range, 15.7-1147.5 cm3). Among all 137 markers tested for correlation with survival data, only programmed death-ligand (PDL) 1 and PDL2 expression showed a prognostic impact. Overall, both PDL1 and PDL2 genes were highly expressed in 37 patients (30.8%; PDL1high/PDL2high). The baseline clinical characteristics of patients with PDL1high/PDL2high were similar to those of other patients. In univariate analysis, PDL1high/PDL2high status was associated with poor progression-free survival (PFS) (hazard ratio [HR], 4.292) and overall survival (OS; HR, 8.24). In multivariate analysis, PDL1high/PDL2high status was an independent prognostic factor of adverse outcomes (PFS: HR, 5.22; OS: HR, 10.368). We validated these results in an independent cohort of 40 patients and confirmed the significant association between PDL1high/PDL2high status and inferior PFS (HR, 6.11). High PDL1/PDL2 gene expression defines a population with strong immune privilege and poorer outcomes from standard chemotherapy who might benefit from firstline checkpoint inhibitor therapy.

Detection of salivary gland and sinonasal fusions by a next-generation sequencing based, ligation-dependent, multiplex RT-PCR assay

AIMS

The discovery of tumour type-specific gene fusion oncogenes in benign and malignant salivary gland and sinonasal (SGSN) tumours has significantly increased our knowledge about their molecular pathology and classification.

METHODS AND RESULTS

We developed a new targeted multiplexed next-generation sequencing (NGS)-based method that utilizes ligation dependent reverse-transcriptase polymerase chain reaction (LD-RT-PCR) to detect oncogenic fusion transcripts involving 116 genes, leading to 96 gene fusions known to be recurrently rearranged in these tumours. In all, 180 SGSN tumours (formalin-fixed, paraffin-embedded samples, 141 specimens and 39 core needle biopsies) from the REFCORpath (French network for rare head and neck cancers) with previously identified fusion genes by fluorescent in situ hybridisation (FISH), RT-PCR, or molecular immunohistochemistry were selected to test its specificity and sensitivity and validate its diagnostic use. Tested tumours encompassed 14 major tumours types, including secretory carcinoma, mucoepidermoid carcinoma, adenoid cystic carcinoma, salivary gland intraductal carcinoma, clear cell carcinoma, pleomorphic adenoma, adamantinoma-like Ewing Sarcoma, EWSR1::COLCA2 sinonasal sarcoma, DEK::AFF2 sinonasal carcinoma, and biphenotypic sinonasal sarcoma. In-frame fusion transcripts were detected in 97.8% of cases (176/180). Gene fusion assay results correlated with conventional techniques (immunohistochemistry [IHC], FISH, and RT-PCR) in 176/180 tumours (97.8%).

CONCLUSION

This targeted multiplexed NGS-based LD-RT-PCR method is a robust, highly sensitive method for the detection of recurrent gene fusions from routine clinical SGSN tumours. It can be easily customized to cover new fusions. These results are promising for implementing an integrated NGS system to rapidly detect genetic aberrations, facilitating accurate, genomics-based diagnoses, and accelerate time to precision therapies in SGSN tumours.

Complete hematologic response after belinostat treatment and allogeneic stem cell transplantation for multiple relapsed/refractory angioimmunoblastic T-cell lymphoma: A case report

Key Clinical Message

This case report highlights the potential of belinostat for the treatment of relapsed/refractory peripheral T-cell lymphomas, for which effective therapies are still scarce.

Abstract

Peripheral T-cell lymphomas have an aggressive disease course associated with poor outcomes. We report a young patient with highly pretreated relapsed/refractory nodal follicular helper T-cell lymphoma (angioimmunoblastic-type [nTFHL-AI]), who successfully received an allogeneic stem cell transplantation following belinostat therapy. The complete hematologic response achieved has lasted more than 2 years.

Disequilibrium between BRCA1 and BRCA2 Circular and Messenger RNAs Plays a Role in Breast Cancer

Breast cancer is a frequent disease for which the discovery of markers that enable early detection or prognostic assessment remains challenging. Circular RNAs (circRNAs) are single-stranded structures in closed loops that are produced by backsplicing. CircRNA and messenger RNA (mRNA) are generated co-transcriptionally, and backsplicing and linear splicing compete against each other. As mRNAs are key players in tumorigenesis, we hypothesize that a disruption of the balance between circRNAs and mRNAs could promote breast cancer. Hence, we developed an assay for a simultaneous study of circRNAs and mRNAs, which we have called splice and expression analyses by exon ligation and high-throughput sequencing (SEALigHTS). Following SEALigHTS validation for BRCA1 and BRCA2, our hypothesis was tested using an independent research set of 95 pairs from tumor and adjacent normal breast tissues. In this research set, ratios of BRCA1 and BRCA2 circRNAs/mRNAs were significantly lower in the tumor breast tissue compared to normal tissue (p = 1.6 × 10^-9 and p = 4.4 × 10^-5 for BRCA1 and BRCA2, respectively). Overall, we developed an innovative method to study linear splicing and backsplicing, described the repertoire of BRCA1 and BRCA2 circRNAs, including 15 novel ones, and showed for the first time that a disequilibrium between BRCA1 and BRCA2 circRNAs and mRNAs plays a role in breast cancer.

A novel SMARCA2-CREM fusion expending the molecular spectrum of salivary gland hyalinazing clear cell carcinoma beyond the FET genes

Hyalinizing clear cell carcinoma (HCCC) is a rare salivary gland carcinoma with a generally indolent behavior, characterized by recurrent chromosomal translocation involving EWSR1 (22q12.2) leading to two fusion genes EWSR1::ATF1 or EWSR1::CREM. We report one case of HCCC with a novel SMARCA2::CREM fusion, identified by targeted RNA next generation sequencing by LD-RT-PCR, which has until now never been described in salivary glands. The exon 4 of SMARCA2 is fused to exon 5 of CREM. This fusion has been described previously in only one tumor, a central nervous system tumor (intracranial mesenchymal tumor) but not in other FET::CREB fused tumors. This fusion was confirmed by CREM break-apart FISH and reverse transcriptase polymerase chain reaction (RT-PCR). The tumor cells showed retained expression of INI1, SMARCA2, and SMARCA4 by immunohistochemistry. We compare its clinical, histopathological, immunophenotypic, genetic features with those previously described in HCCC, FET::CREB fusion-positive. Our results added data suggesting that different histomolecular tumor subtypes seem to be included within the terminology "HCCC, FET::CREB fusion-positive," and that further series of cases are needed to better characterize them.

A Simple Molecular Tool for the Assessment of Kidney Transplant Biopsies

BACKGROUND

The Banff Classification for Allograft Pathology recommendations for the diagnosis of kidney transplant rejection includes molecular assessment of the transplant biopsy. However, implementation of molecular tools in clinical practice is still limited, partly due to the required expertise and financial investment. The reverse transcriptase multiplex ligation-dependent probe amplification (RT-MLPA) assay is a simple, rapid, and inexpensive assay that permits simultaneous evaluation of a restricted gene panel using paraffin-embedded tissue blocks. The aim of this study was to develop and validate a RT-MLPA assay for diagnosis and classification of rejection.

METHODS

A retrospective cohort of 220 kidney transplant biopsies from two centers, which included 52 antibody-mediated rejection, 51 T-cell-mediated rejection, and 117 no-rejection controls, was assessed. A 17-gene panel was identified on the basis of relevant pathophysiological pathways. A support vector machine classifier was developed. A subset of 109 biopsies was also assessed using the Nanostring Banff Human Organ Transplant panel to compare the two assays.

RESULTS

The support vector machine classifier train and test accuracy scores were 0.84 and 0.83, respectively. In the test cohort, the F1 score for antibody-mediated rejection, T-cell-mediated rejection, and control were 0.88, 0.86, and 0.69, respectively. Using receiver-operating characteristic curves, the area under the curve for class predictions was 0.96, 0.89, and 0.91, respectively, with a weighted average at 0.94. Classifiers' performances were highest for antibody-mediated rejection diagnosis with 94% correct predictions, compared with 88% correct predictions for control biopsies and 60% for T-cell-mediated rejection biopsies. Gene expression levels assessed by RT-MLPA and Nanostring were correlated: r = 0.68, P < 0.001. Equivalent gene expression profiles were obtained with both assays in 81% of the samples.

CONCLUSIONS

The 17-gene panel RT-MLPA assay, developed here for formalin-fixed paraffin-embedded kidney transplant biopsies, classified kidney transplant rejection with an overall accurate prediction ratio of 0.83.

PODCAST

This article contains a podcast at https://dts.podtrac.com/redirect.mp3/www.asn-online.org/media/podcast/CJASN/2023_04_10_CJN10100822.mp3.

Detection of sarcoma fusions by a next-generation sequencing based-ligation-dependent multiplex RT-PCR assay

Morphological, immunohistochemical, and molecular methods often need to be combined for accurate diagnosis and optimal clinical management of sarcomas. Here, we have developed, a new molecular diagnostic assay, for the detection of gene fusions in sarcomas. This targeted multiplexed next-generation sequencing (NGS)-based method utilizes ligation dependent reverse-transcriptase polymerase chain reaction (LD-RT-PCR-NGS) to detect oncogenic fusion transcripts involving 137 genes, leading to 139 gene fusions known to be recurrently rearranged in soft-tissue and bone tumors. 158 bone and soft-tissue tumors with previously identified fusion genes by fluorescent in situ hybridization (FISH) or RT-PCR were selected to test the specificity and the sensitivity of this assay. RNA were extracted from formalin-fixed paraffin-embedded (n = 143) or frozen (n = 15) material (specimen; n = 42 or core needle biopsies; n = 116). Tested tumors encompassed 23 major translocation-related sarcomas types, including Ewing and Ewing-like sarcomas, rhabdomyosarcomas, desmoplastic small round-cell tumors, clear-cell sarcomas, infantile fibrosarcomas, endometrial stromal sarcomas, epithelioid hemangioendotheliomas, alveolar soft-part sarcomas, biphenotypic sinonasal sarcomas, extraskeletal myxoid chondrosarcomas, myxoid/round-cell liposarcomas, dermatofibrosarcomas protuberans and solitary fibrous tumors. In-frame fusion transcripts were detected in 98.1% of cases (155/158). Gene fusion assay results correlated with conventional techniques (FISH and RT-PCR) in 155/158 tumors (98.1%). These data demonstrate that this assay is a rapid, robust, highly sensitive, and multiplexed targeted RNA sequencing assay for the detection of recurrent gene fusions on RNA extracted from routine clinical specimens of sarcomas (formalin-fixed paraffin-embedded or frozen). It facilitates the precise diagnosis and identification of tumors with potential targetable fusions. In addition, this assay can be easily customized to cover new fusions.

Clinico-biological features of T-cell acute lymphoblastic leukemia with fusion proteins

T-cell acute lymphoblastic leukemias (T-ALL) represent 15% of pediatric and 25% of adult ALL. Since they have a particularly poor outcome in relapsed/refractory cases, identifying prognosis factors at diagnosis is crucial to adapting treatment for high-risk patients. Unlike acute myeloid leukemia and BCP ALL, chromosomal rearrangements leading to chimeric fusion-proteins with strong prognosis impact are sparsely reported in T-ALL. To address this issue an RT-MPLA assay was applied to a consecutive series of 522 adult and pediatric T-ALLs and identified a fusion transcript in 20% of cases. PICALM-MLLT10 (4%, n = 23), NUP214-ABL1 (3%, n = 19) and SET-NUP214 (3%, n = 18) were the most frequent. The clinico-biological characteristics linked to fusion transcripts in a subset of 235 patients (138 adults in the GRAALL2003/05 trials and 97 children from the FRALLE2000 trial) were analyzed to identify their prognosis impact. Patients with HOXA trans-deregulated T-ALLs with MLLT10, KMT2A and SET fusion transcripts (17%, 39/235) had a worse prognosis with a 5-year EFS of 35.7% vs 63.7% (HR = 1.63; p = 0.04) and a trend for a higher cumulative incidence of relapse (5-year CIR = 45.7% vs 25.2%, HR = 1.6; p = 0.11). Fusion transcripts status in T-ALL can be robustly identified by RT-MLPA, facilitating risk adapted treatment strategies for high-risk patients.

UMI-Varcal: A Low-Frequency Variant Caller for UMI-Tagged Paired-End Sequencing Data

The rapid transition from traditional sequencing methods to Next-Generation Sequencing (NGS) has allowed for a faster and more accurate detection of somatic variants (Single-Nucleotide Variant (SNV) and Copy Number Variation (CNV)) in tumor cells. NGS technologies require a succession of steps during which false variants can be silently added at low frequencies. Filtering these artifacts can be a rather difficult task especially when the experiments are designed to look for very low frequency variants. Recently, adding unique molecular barcodes called UMI (Unique Molecular Identifier) to the DNA fragments appears to be a very effective strategy to specifically filter out false variants from the variant calling results (Kukita et al. DNA Res 22(4):269-277, 2015; Newman et al. Nat Biotechnol 34(5):547-555, 2016; Schmitt et al. Proc Natl Acad Sci U S A 109(36):14508-14513). Here, we describe UMI-VarCal (Sater et al. Bioinformatics 36:2718-2724, 2020), which can use the UMI information from UMI-tagged reads to offer a faster and more accurate variant calling analysis.

Pan-TRK Immunohistochemistry Is Highly Correlated With NTRK3 Gene Rearrangements in Salivary Gland Tumors

AIMS

Secretory carcinoma (SC) is characterized by ETV6 rearrangements, most often ETV6-NTRK3 fusion. Given its histologic overlap with other salivary gland tumors (SGTs), SCs can be difficult to diagnose without genetic confirmation. A recently developed pan-TRK (tropomyosin receptor kinase) antibody shows promise for identifying tumors with NTRK (neurotrophic tyrosine kinase receptor 3) fusions. The aim of this study was to evaluate the utility of pan-TRK immunohistochemistry in distinguishing SCs from mimics and selecting patients eligible for TRK inhibitor clinical trials. We examined whole-tissue sections from 111 SGTs with molecular characterization, including 26 SCs (23 with ETV6-NTRK3 fusion and 3 with ETV6-RET fusion detected by ligation-dependent reverse transcription-polymerase chain reaction, next-generation sequencing and 85 non-SC SGTs (no ETV6-NTRK3 fusion). Immunohistochemistry was performed with a pan-TRK rabbit monoclonal antibody. When any pan-TRK staining (nuclear or cytoplasmic with any staining intensity) was considered to indicate positivity, 22 of 23 SCs with ETV6-NTRK3 fusion (95.7%) and 33 of 85 non-SC (38.8%) salivary neoplasms were positive, mainly basal cell adenoma, pleomorphic adenomas, adenoid cystic carcinomas, and epithelial-myoepithelial carcinomas. All SCs with ETV6-RET fusion were entirely negative. When only nuclear pan-TRK staining with any staining intensity was considered positive, 18 of 23 SCs with ETV6-NTRK3 fusion (78.3%) were positive, 11 among them with diffuse staining (>30% of cells). All non-SCs and SCs with ETV6-RET fusion were entirely negative. In comparison to molecular analysis (ligation-dependent reverse transcription-polymerase chain reaction, next-generation sequencing), nuclear pan-TRK IHC has a sensitivity of 78.3% and a specificity of 100% for diagnosing SCs with ETV6-NTRK3 fusion, 69% and 100% for SCs (all fusions). Pan-TRK is a reasonable screening test for diagnosing SCs among SGTs when taking only nuclear staining into account. Although pan-TRK expression is not entirely sensitive for SCs, nuclear staining is highly specific for SCs with ETV6-NTRK3 fusion. The lack of pan-TRK immunoreactivity in a subset of SCs is suggestive of atypical exons 4 to 14 or exons 5 to 14 ETV6-NTRK3 fusion or non-NTRK alternative fusion partners such as ETV6-RET. Pan-TRK staining can serve as a strong diagnostic marker to distinguish SC from it mimics and to select patients eligible for TRK inhibitor clinical trials.

Concomitant occurrence of genetically distinct Hodgkin lymphoma and primary mediastinal lymphoma

Synchronous Hodgkin Lymphoma and Primary Mediastinal B-cell Lymphoma is possible, with molecular analyses proving the absence of clonal filiation between both entities. This suggests a common etiology but the existence of two divergent clones.

Detection of Gene Fusion Transcripts in Peripheral T-Cell Lymphoma Using a Multiplexed Targeted Sequencing Assay

The genetic basis of peripheral T-cell lymphoma (PTCL) is complex and encompasses several recurrent fusion transcripts discovered over the past years by means of massive parallel sequencing. However, there is currently no affordable and rapid technology for their simultaneous detection in clinical samples. Herein, we developed a multiplex ligation-dependent RT-PCR-based assay, followed by high-throughput sequencing, to detect 33 known PTCL-associated fusion transcripts. Anaplastic lymphoma kinase (ALK) fusion transcripts were detected in 15 of 16 ALK-positive anaplastic large-cell lymphomas. The latter case was further characterized by a novel SATB1_ALK fusion transcript. Among 239 other PTCLs, representative of nine entities, non-ALK fusion transcripts were detected in 24 samples, mostly of follicular helper T-cell (TFH) derivation. The most frequent non-ALK fusion transcript was ICOS_CD28 in nine TFH-PTCLs, one PTCL not otherwise specified, and one adult T-cell leukemia/lymphoma, followed by VAV1 rearrangements with multiple partners (STAP2, THAP4, MYO1F, and CD28) in five samples (three PTCL not otherwise specified and two TFH-PTCLs). The other rearrangements were CTLA4_CD28 (one TFH-PTCL), ITK_SYK (two TFH-PTCLs), ITK_FER (one TFH-PTCL), IKZF2_ERBB4 (one TFH-PTCL and one adult T-cell leukemia/lymphoma), and TP63_TBL1XR1 (one ALK-negative anaplastic large-cell lymphoma). All fusions detected by our assay were validated by conventional RT-PCR and Sanger sequencing in 30 samples with adequate material. The simplicity and robustness of this targeted multiplex assay make it an attractive tool for the characterization of these heterogeneous diseases.

An improved assay for detection of theranostic gene translocations and MET exon 14 skipping in thoracic oncology

Theranostic translocations may be difficult to detect by routine techniques, especially when specimens are exiguous. We recently demonstrated in a series of translocated lung adenocarcinomas that LD-RT-PCR (ligation-dependent reverse transcription polymerase chain reaction) assay could identify ALK, ROS1 and RET rearrangements with 64% sensitivity and 100% specificity. Here, we report an upgraded version of this assay used in a routine prospective cohort of lung carcinomas. Newly diagnosed lung carcinomas referred to the Rouen molecular platform between 15/05/2018 and 15/05/2019 for ALK and ROS1 IHC, genotyping (SNaPshot© +/- high-throughput genotyping) and sometimes FISH (standard routine process) were tested prospectively in parallel with the LD-RT-PCR assay designed to detect at one go ALK, ROS1 and RET translocations and MET exon 14 skipping. 413 tumors from 396 patients were included. LD-RT-PCR had a global sensitivity of 91.43% (standard routine process: 80%), with a specificity of 100%. It detected 15/18 ALK and 4/4 ROS1 translocated tumors, but also 6/6 tumors with MET exon 14 skipping retrieved by genotyping. In addition, it retrieved 7 alterations missed by the routine process, then confirmed by other means: 5 MET exon 14 skipping and 2 RET translocated tumors. Finally, it allowed to deny an effect on MET exon 14 skipping for 8 mutations detected by routine genotyping. We successfully implemented LD-RT-PCR in routine analysis. This technique is cheap, fast, sensitive, specific, and easily upgradable (e.g., NTRK translocations), but still requires IHC to be performed in parallel. Owing to its advantages, we recommend considering it, in parallel with IHC and genotyping, as an excellent cost-effective alternative, for the systematic testing of lung adenocarcinoma, to FISH and to more expensive and complex assays such as RNA-seq.

Improving high-resolution copy number variation analysis from next generation sequencing using unique molecular identifiers

BACKGROUND

Recently, copy number variations (CNV) impacting genes involved in oncogenic pathways have attracted an increasing attention to manage disease susceptibility. CNV is one of the most important somatic aberrations in the genome of tumor cells. Oncogene activation and tumor suppressor gene inactivation are often attributed to copy number gain/amplification or deletion, respectively, in many cancer types and stages. Recent advances in next generation sequencing protocols allow for the addition of unique molecular identifiers (UMI) to each read. Each targeted DNA fragment is labeled with a unique random nucleotide sequence added to sequencing primers. UMI are especially useful for CNV detection by making each DNA molecule in a population of reads distinct.

RESULTS

Here, we present molecular Copy Number Alteration (mCNA), a new methodology allowing the detection of copy number changes using UMI. The algorithm is composed of four main steps: the construction of UMI count matrices, the use of control samples to construct a pseudo-reference, the computation of log-ratios, the segmentation and finally the statistical inference of abnormal segmented breaks. We demonstrate the success of mCNA on a dataset of patients suffering from Diffuse Large B-cell Lymphoma and we highlight that mCNA results have a strong correlation with comparative genomic hybridization.

CONCLUSION

We provide mCNA, a new approach for CNV detection, freely available at https://gitlab.com/pierrejulien.viailly/mcna/ under MIT license. mCNA can significantly improve detection accuracy of CNV changes by using UMI.

Subcutaneous Rituximab-MiniCHOP Compared With Subcutaneous Rituximab-MiniCHOP Plus Lenalidomide in Diffuse Large B-Cell Lymphoma for Patients Age 80 Years or Older

PURPOSE

The prognosis of elderly patients with diffuse large B-cell lymphoma (DLBCL) is worse than that of young patients. An attenuated dose of chemotherapy-cyclophosphamide, doxorubicin, vincristine, and prednisone plus rituximab (R-miniCHOP)-is a good compromise between efficacy and safety in very elderly patients. In combination with R-CHOP (R2-CHOP), lenalidomide has an acceptable level of toxicity and may mitigate the negative prognosis of the non-germinal center B-cell-like phenotype. The Lymphoma Study association conducted a multicentric, phase III, open-label, randomized trial to compare R-miniCHOP and R2-miniCHOP.

PATIENTS AND METHODS

Patients of age 80 years or older with untreated DLBCL were randomly assigned into the R-miniCHOP21 group or the R2-miniCHOP21 group for six cycles and stratified according to CD10 expression and age. The first cycle of rituximab was delivered by IV on D1 after a prephase and then delivered subcutaneously on D1 of cycles 2-6. Lenalidomide was delivered at a dose of 10 mg once daily on D1-D14 of each cycle. The primary end point was overall survival (OS).

RESULTS

A total of 249 patients with new DLBCL were randomly assigned (127 R-miniCHOP and 122 R2-miniCHOP). The median age was 83 years (range, 80-96), and 55% of the patients were classified as non-GCB. The delivered dose for each R-miniCHOP compound was similar in both arms. Over a median follow-up of 25.1 months, the intention-to-treat analysis revealed that R2-miniCHOP did not improve OS (2-year OS 66% in R-miniCHOP and 65.7% in R2-miniCHOP arm, P = .98) in the overall population or in the non-GCB population. Grade 3-4 adverse events occurred in 53% of patients with R-miniCHOP and in 81% of patients with R2-miniCHOP.

CONCLUSION

The addition of lenalidomide to R-miniCHOP does not improve OS. Rituximab delivered subcutaneously was safe in this population.

UMI-VarCal: a new UMI-based variant caller that efficiently improves low-frequency variant detection in paired-end sequencing NGS libraries

MOTIVATION

Next-generation sequencing has become the go-to standard method for the detection of single-nucleotide variants in tumor cells. The use of such technologies requires a PCR amplification step and a sequencing step, steps in which artifacts are introduced at very low frequencies. These artifacts are often confused with true low-frequency variants that can be found in tumor cells and cell-free DNA. The recent use of unique molecular identifiers (UMI) in targeted sequencing protocols has offered a trustworthy approach to filter out artefactual variants and accurately call low-frequency variants. However, the integration of UMI analysis in the variant calling process led to developing tools that are significantly slower and more memory consuming than raw-reads-based variant callers.

RESULTS

We present UMI-VarCal, a UMI-based variant caller for targeted sequencing data with better sensitivity compared to other variant callers. Being developed with performance in mind, UMI-VarCal stands out from the crowd by being one of the few variant callers that do not rely on SAMtools to do their pileup. Instead, at its core runs an innovative homemade pileup algorithm specifically designed to treat the UMI tags in the reads. After the pileup, a Poisson statistical test is applied at every position to determine if the frequency of the variant is significantly higher than the background error noise. Finally, an analysis of UMI tags is performed, a strand bias and a homopolymer length filter are applied to achieve better accuracy. We illustrate the results obtained using UMI-VarCal through the sequencing of tumor samples and we show how UMI-VarCal is both faster and more sensitive than other publicly available solutions.

AVAILABILITY AND IMPLEMENTATION

The entire pipeline is available at https://gitlab.com/vincent-sater/umi-varcal-master under MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Correlations between baseline 18F-FDG PET tumour parameters and circulating DNA in diffuse large B cell lymphoma and Hodgkin lymphoma

Background

¹⁸F-FDG PET/CT is a standard for many B cell malignancies, while blood DNA measurements are emerging tools. Our objective was to evaluate the correlations between baseline PET parameters and circulating DNA in diffuse large B cell lymphoma (DLBCL) and classical Hodgkin lymphoma (cHL).

Methods

Twenty-seven DLBCL and forty-eight cHL were prospectively included. Twelve PET parameters were analysed. Spearman’s correlations were used to compare PET parameters each other and to circulating cell-free DNA ([cfDNA]) and circulating tumour DNA ([ctDNA]). p values were controlled by Benjamini–Hochberg correction.

Results

Among the PET parameters, three different clusters for tumour burden, fragmentation/massiveness and dispersion parameters were observed. Some PET parameters were significantly correlated with blood DNA parameters, including the total metabolic tumour surface (TMTS) describing the tumour–host interface (e.g. ρ = 0.81 p < 0.001 for [ctDNA] of DLBLC), the tumour median distance between the periphery and the centroid (medPCD) describing the tumour’s massiveness (e.g. ρ = 0.81 p < 0.001 for [ctDNA] of DLBLC) and the volume of the bounding box including tumours (TumBB) describing the disease’s dispersion (e.g. ρ = 0.83 p < 0.001 for [ctDNA] of DLBLC).

Conclusions

Some PET parameters describing tumour burden, fragmentation/massiveness and dispersion are significantly correlated with circulating DNA parameters of DLBCL and cHL patients. These results could help to understand the pathophysiology of B cell malignancies.

UMI-Gen: A UMI-based read simulator for variant calling evaluation in paired-end sequencing NGS libraries

Motivation

With Next Generation Sequencing becoming more affordable every year, NGS technologies asserted themselves as the fastest and most reliable way to detect Single Nucleotide Variants (SNV) and Copy Number Variations (CNV) in cancer patients. These technologies can be used to sequence DNA at very high depths thus allowing to detect abnormalities in tumor cells with very low frequencies. Multiple variant callers are publicly available and are usually efficient at calling out variants. However, when frequencies begin to drop under 1%, the specificity of these tools suffers greatly as true variants at very low frequencies can be easily confused with sequencing or PCR artifacts. The recent use of Unique Molecular Identifiers (UMI) in NGS experiments has offered a way to accurately separate true variants from artifacts. UMI-based variant callers are slowly replacing raw-read based variant callers as the standard method for an accurate detection of variants at very low frequencies. However, benchmarking done in the tools publication are usually realized on real biological data in which real variants are not known, making it difficult to assess their accuracy.

Results

We present UMI-Gen, a UMI-based read simulator for targeted sequencing paired-end data. UMI-Gen generates reference reads covering the targeted regions at a user customizable depth. After that, using a number of control files, it estimates the background error rate at each position and then modifies the generated reads to mimic real biological data. Finally, it will insert real variants in the reads from a list provided by the user.

Availability

The entire pipeline is available at https://gitlab.com/vincent-sater/umigen under MIT license.

Simultaneous detection of EGFR amplification and EGFRvIII variant using digital PCR-based method in glioblastoma

Epidermal growth factor receptor (EGFR) amplification and EGFR variant III (EGFRvIII, deletion of exons 2-7) are of clinical interest for glioblastoma. The aim was to develop a digital PCR (dPCR)-based method using locked nucleic acid (LNA)-based hydrolysis probes, allowing the simultaneous detection of the EGFR amplification and EGFRvIII variant. Sixty-two patients were included. An exploratory cohort (n = 19) was used to develop the dPCR assay using three selected amplicons within the EGFR gene, targeting intron 1 (EGFR1), junction of exon 3 and intron 3 (EGFR2) and intron 22 (EGFR3). The copy number of EGFR was estimated by the relative quantification of EGFR1, EGFR2 and EGFR3 amplicon droplets compared to the droplets of a reference gene. EGFRvIII was identified by comparing the copy number of the EGFR2 amplicon to either the EGFR1 or EGFR3 amplicon. dPCR results were compared to fluorescence in situ hybridization (FISH) and next-generation sequencing for amplification; and to RT-PCR-based method for EGFRvIII. The dPCR assay was then tested in a validation cohort (n = 43). A total of 8/19 EGFR-amplified and 5/19 EGFRvIII-positive tumors were identified in the exploratory cohort. Compared to FISH, the EGFR3 dPCR assay detected all EGFR-amplified tumors (8/8, 100%) and had the highest concordance with the copy number estimation by NGS. The concordance between RT-PCR and dPCR was also 100% for detecting EGFRvIII using an absolute difference of 10.8 for the copy number between EGFR2 and EGFR3 probes. In the validation cohort, the sensitivity and specificity of dPCR using EGFR3 probes were 100% for the EGFR amplification detection compared to FISH (19/19). EGFRvIII was detected by dPCR in 8 EGFR-amplified patients and confirmed by RT-PCR. Compared to FISH, the EGFR2/EGFR3 dPCR assay was estimated with a one-half cost value. These results highlight that dPCR allowed the simultaneous detection of EGFR amplification and EGFRvIII for glioblastoma.

Refining diffuse large B-cell lymphoma subgroups using integrated analysis of molecular profiles

Background

Gene expression profiling (GEP), next-generation sequencing (NGS) and copy number variation (CNV) analysis have led to an increasingly detailed characterization of the genomic profiles of DLBCL. The aim of this study was to perform a fully integrated analysis of mutational, genomic, and expression profiles to refine DLBCL subtypes. A comparison of our model with two recently published integrative DLBCL classifiers was carried out, in order to best reflect the current state of genomic subtypes.

Methods

223 patients with de novo DLBCL from the prospective, multicenter and randomized LNH-03B LYSA clinical trials were included. GEP data was obtained using Affymetrix GeneChip arrays, mutational profiles were established by Lymphopanel NGS targeting 34 key genes, CNV analysis was obtained by array CGH, and FISH and IHC were performed. Unsupervised independent component analysis (ICA) was applied to GEP data and integrated analysis of multi-level molecular data associated with each component (gene signature) was performed.

Findings

ICA identified 38 components reflecting transcriptomic variability across our DLBCL cohort. Many of the components were closely related to well-known DLBCL features such as cell-of-origin, stromal and MYC signatures. A component linked to gain of 19q13 locus, among other genomic alterations, was significantly correlated with poor OS and PFS. Through this integrated analysis, a high degree of heterogeneity was highlighted among previously described DLBCL subtypes.

Interpretation

The results of this integrated analysis enable a global and multi-level view of DLBCL, as well as improve our understanding of DLBCL subgroups.

Defining signatures of peripheral T-cell lymphoma with a targeted 20-marker gene expression profiling assay

Peripheral T-cell lymphoma comprises a heterogeneous group of mature non-Hodgkin lymphomas. Their diagnosis is challenging, with up to 30% of cases remaining unclassifiable and referred to as “not otherwise specified”. We developed a reverse transcriptase-multiplex ligation-dependent probe amplification gene expression profiling assay to differentiate the main T-cell lymphoma entities and to study the heterogeneity of the “not specified” category. The test evaluates the expression of 20 genes, including 17 markers relevant to T-cell immunology and lymphoma biopathology, one Epstein-Barr virus-related transcript, and variants of RHOA (G17V) and IDH2 (R172K/T). By unsupervised hierarchical clustering, our assay accurately identified 21 of 21 ALK-positive anaplastic large cell lymphomas, 16 of 16 extranodal natural killer (NK)/T-cell lymphomas, 6 of 6 hepatosplenic T-cell lymphomas, and 13 of 13 adult T-cell leukemia/lymphomas. ALK-negative anaplastic lymphomas (n=34) segregated into one cytotoxic cluster (n=10) and one non-cytotoxic cluster expressing Th2 markers (n=24) and enriched in DUSP22-rearranged cases. The 63 T_FH-derived lymphomas divided into two subgroups according to a predominant T_FH (n=50) or an enrichment in Th2 (n=13) signatures. We next developed a support vector machine predictor which attributed a molecular class to 27 of 77 not specified T-cell lymphomas: 17 T_FH, five cytotoxic ALK-negative anaplastic and five NK/T-cell lymphomas. Among the remaining cases, we identified two cell-of-origin subgroups corresponding to cytotoxic/Th1 (n=19) and Th2 (n=24) signatures. A reproducibility test on 40 cases yielded a 90% concordance between three independent laboratories. This study demonstrates the applicability of a simple gene expression assay for the classification of peripheral T-cell lymphomas. Its applicability to routinely-fixed samples makes it an attractive adjunct in diagnostic practice.

A recurrent clonally distinct Burkitt lymphoma case highlights genetic key events contributing to oncogenesis

Burkitt lymphoma (BL) is characterized by a translocation of the MYC oncogene that leads to the upregulation of MYC expression, cell growth and proliferation. It is well-established that MYC translocation is not a sufficient genetic event to cause BL. Next-generation sequencing has recently provided a comprehensive analysis of the landscape of additional genetic events that contribute to BL lymphomagenesis. Refractory BL or relapsing BL are almost always incurable as a result of the selection of a highly chemoresistant clonally related cell population. Conversely, a few BL recurrence cases arising from clonally distinct tumors have been reported and were associated with a favorable outcome similar to that reported for first-line treatment. Here, we used an unusual case of recurrent but clonally distinct EBV+ BL to highlight the key genetic events that drive BL lymphomagenesis. By whole exome sequencing, we established that ID3 gene was targeted by distinct mutations in the two clonally unrelated diseases, highlighting the crucial role of this gene during lymphomagenesis. We also detected a heterozygous E1021K PIK3CD mutation, thus increasing the spectrum of somatic mutations altering the PI3K signaling pathway in BL. Interestingly, this mutation is known to be associated with activated phosphoinositide 3-kinase delta syndrome (APDS). Finally, we also identified an inherited heterozygous truncating c.5791CT FANCM mutation that may contribute to the unusual recurrence of BL.

Non-invasive monitoring of diffuse large B-cell lymphoma by cell-free DNA high-throughput targeted sequencing: analysis of a prospective cohort

From a liquid biopsy, cell-free DNA (cfDNA) can provide information regarding basal tumoral genetic patterns and changes upon treatment. In a prospective cohort of 30 diffuse large B-cell lymphomas (DLBCL), we determined the clinical relevance of cfDNA using targeted next-generation sequencing and its correlation with PET scan imaging at the time of diagnosis and during treatment. Using a dedicated DLBCL panel, mutations were identified at baseline for 19 cfDNAs and profiles were consistent with expected DLBCL patterns. Tumor burden-related clinical and PET scan features (LDH, IPI, and metabolic tumor volume) were significantly correlated with the quantity of tumoral cfDNA. Among the four patients presenting additional mutations in their cfDNAs, three had high metabolic tumor volumes, suggesting that cfDNA more accurately reflects tumor heterogeneity than tissues biopsy itself. Mid-treatment, four patients still had basal mutations in their cfDNAs, including three in partial response according to their Deauville scores. Our study highlights the major interests in liquid biopsy, in particular in the context of bulky tumors where cfDNA allows capturing the entire tumoral mutation profile. Therefore, cfDNA analysis in DLBCL represents a complementary approach to PET scan imaging.

Application of the cghRA framework to the genomic characterization of Diffuse Large B-Cell Lymphoma

Motivation

Although sequencing-based technologies are becoming the new reference in genome analysis, comparative genomic hybridization arrays (aCGH) still constitute a simple and reliable approach for copy number analysis. The most powerful algorithms to analyze such data have been freely provided by the scientific community for many years, but combining them is a complex scripting task.

Results

The cghRA framework combines a user-friendly graphical interface and a powerful object-oriented command-line interface to handle a full aCGH analysis, as is illustrated in an original series of 107 Diffuse Large B-Cell Lymphomas. New algorithms for copy-number calling, polymorphism detection and minimal common region prioritization were also developed and validated. While their performances will only be demonstrated with aCGH, these algorithms could actually prove useful to any copy-number analysis, whatever the technique used.

Availability and implementation

R package and source for Linux, MS Windows and MacOS are freely available at http://bioinformatics.ovsa.fr/cghRA.

Contact

mareschal@ovsa.fr or fabrice.jardin@chb.unicancer.fr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Non-invasive detection of somatic mutations using next-generation sequencing in primary central nervous system lymphoma

Purpose

Primary central nervous system lymphomas (PCNSL) have recurrent genomic alterations. The main objective of our study was to demonstrate that targeted sequencing of circulating cell-free DNA (cfDNA) released by PCNSL at the time of diagnosis could identify somatic mutations by next-generation sequencing (NGS).

Patients and Methods

PlasmacfDNA and matched tumor DNA (tDNA) from 25 PCNSL patients were sequenced using an Ion Torrent Personal Genome Machine (Life Technologies^®). First, patient-specific targeted sequencing of identified somatic mutations in tDNA was performed. Then, a second sequencing targeting MYD88 c.T778C was performed and compared to plasma samples from 25 age-matched control patients suffering from other types of cancer.

Results

According to the patient-specific targeted sequencing, eight patients (32% [95% CI 15-54%]) had detectable somatic mutations in cfDNA. Considering MYD88 sequencing, six patients had the specific c.T778C alteration detected in plasma. Using a control group, the sensitivity was 24% [9-45%] and the specificity was 100%. Tumor volume or deep brain structure involvement did not influence the detection of somatic mutations in plasma.

Conclusion

This pilot study provided evidence that somatic mutations can be detected by NGS in the cfDNA of a subset of patients suffering from PCNSL.

BCL2 expression but not MYC and BCL2 coexpression predicts survival in elderly patients with diffuse large B-cell lymphoma independently of cell of origin in the phase 3 LNH03-6B trial

Background

Our aim was to evaluate whether the cell of origin (COO) as defined by the Hans algorithm and MYC/BCL2 coexpression, which are the two main biological risk factors in elderly patients treated with rituximab, cyclophosphamide, doxorubicin hydrochloride, vincristine, and prednisolone (R-CHOP), maintain their prognostic value in a large prospective clinical trial.

Patients and methods

We evaluated 285 paraffin-embedded samples from patients (60-80 years of age) enrolled in the Lymphoma Study Association trial LNH03-6B who were treated with R-CHOP. We correlated the COO defined by the transcriptome according to the Wright algorithm with that defined by the Hans algorithm in a subset of 62 tumors with available frozen tissue samples.

Results

The non-germinal center B-cell-like phenotype according to the Hans algorithm and BCL2 expression (but not MYC and BCL2 coexpression) predicted worse progression-free survival [hazard ratio (HR)=1.78, P = 0.003 and HR = 1.79, P = 0.003, respectively] and overall survival (HR = 1.85, P = 0.005 and HR = 1.67, P = 0.02, respectively) independently of the International Prognostic Index. The correlation between the Hans algorithm and the Wright algorithm was 91%, with an almost perfect concordance according to a kappa test (0.81).

Conclusions

Our results suggest that immunohistochemically defined COO remains a useful tool for predicting prognosis in diffuse large B-cell lymphoma when performed under optimized standardized conditions and that BCL2 expression may help to identify elderly patients at risk for relapse and who could potentially respond to anti-BCL2 targeted agents. In this prospective phase III trial, the coexpression of MYC and BCL2 does not appear to predict worse survival.

Clinical trial Number

NCT00144755.

Ligation-dependent RT-PCR: a new specific and low-cost technique to detect ALK, ROS, and RET rearrangements in lung adenocarcinoma

Detection of anaplastic lymphoma kinase (ALK), ROS proto-oncogene 1 (ROS1), and rearranged during transfection (RET) gene rearrangements in lung adenocarcinoma is usually performed by immunohistochemistry (IHC) screening followed by fluorescence in situ hybridization (FISH), which is an expensive and difficult technique. Ligation-dependent reverse transcription polymerase chain reaction (RT-PCR) multiplex technique can detect gene rearrangements using probes specifically hybridized to either side of the break point. PCR products are then sequenced by pyrosequencing or high throughput sequencing in order to identify the two genes involved. The reagent cost is <15 dollars per patient and results are available in 2 days. We have developed a 47-probe LD-RT-PCR kit especially for lung adenocarcinomas. Thirty-nine lung adenocarcinomas were studied: 24 ALK+, 14 ROS1+, and 1 RET+. ALK+ and ROS1+ were IHC+ (D5F3 Ventana for ALK and D4D6 Cell Signaling Technology for ROS1) and all cases were FISH+ (Vysis ALK Breakapart Probe Abbott for ALK, Zytolight SPEC ROS1 Dualcolor Breakapart Probe for ROS1 and Zytolight SPEC RET Dual Color Breakapart for RET); 14 wild type samples were included as negative controls. Using LD-RT-PCR, 15 rearrangements (63%) were detected in the ALK cases (gene partner: EML4 in all cases), 9 rearrangements (64%) in the ROS1 cases (gene partners: CD74 in 8 cases and SLC34A2 in 1 case) and 1 (100%) in the single RET case (gene partner: KIF5B). No rearrangement was found in the 14 negative control cases. Negative cases using LD-RT-PCR could be explained by the fact that some partner genes were not included in our assay and therefore could not be detected. Because it is an affordable, fast, and very simple technique, we propose using LD-RT-PCR when ALK immunostaining is positive. For LD-RT-PCR-negative cases, samples should then be analyzed by FISH.

Combination of baseline FDG PET/CT total metabolic tumour volume and gene expression profile have a robust predictive value in patients with diffuse large B-cell lymphoma

PURPOSE

This study evaluated the predictive significance of total metabolic tumour volume (TMTV) measured on baseline FDG PET/CT and its value in addition to gene expression profiling using a new method of gene analysis (rapid reverse transcriptase multiplex ligation-dependent probe amplification assay, RT-MLPA) in patients with diffuse large B-cell lymphoma treated with R-CHOP or R-CHOP-like chemotherapies.

METHODS

The analysis included 114 patients. TMTV was measured using a 41% SUVmax threshold and tumours were classified into GCB or ABC subtypes according to the RT-MLPA assay.

RESULTS

The median follow-up was 40 months. the 5-year progression-free survival (PFS) was 54% and the 5-year overall survival (OS) was 62%. The optimal TMTV cut-off value was 261 cm³. In 59 patients with a high TMTV the 5-year PFS and OS were 37% and 39%, respectively, in comparison with 72% and 83%, respectively, in 55 patients with a low TMTV (p = 0.0002 for PFS, p < 0.0001 for OS). ABC status was significantly associated with a worse prognosis. TMTV combined with molecular data identified three groups with very different outcomes: (1) patients with a low TMTV whatever their phenotype (n = 55), (2) patients with a high TMTV and GCB phenotype (n = 33), and (3) patients with a high TMTV and ABC phenotype (n = 26). In the three groups, 5-year PFS rates were 72%, 51% and 17% (p < 0.0001), and 5-year OS rates were 83%, 55% and 17% (p < 0.0001), respectively. In multivariate analysis, TMTV, ABC/GCB phenotype and International Prognostic Index were independent predictive factors for both PFS and OS (p < 0.05 for both).

CONCLUSIONS

This integrated risk model could lead to more accurate selection of patients that would allow better individualization of therapy.

Determination of Molecular Subtypes of Diffuse Large B-Cell Lymphoma Using a Reverse Transcriptase Multiplex Ligation-Dependent Probe Amplification Classifier: A CALYM Study

Diffuse large B-cell lymphoma (DLBCL) is the most common non-Hodgkin lymphoma. It includes three major subtypes termed germinal center B-cell-like, activated B-cell-like, and primary mediastinal B-cell lymphoma. With the emergence of novel targeted therapies, accurate methods capable of interrogating this cell-of-origin classification should soon become essential in the clinics. To address this issue, we developed a novel gene expression profiling DLBCL classifier based on reverse transcriptase multiplex ligation-dependent probe amplification. This assay simultaneously evaluates the expression of 21 markers, to differentiate primary mediastinal B-cell lymphoma, activated B-cell-like, germinal center B-cell-like, and also Epstein-Barr virus-positive DLBCLs. It was trained using 70 paraffin-embedded biopsies and validated using >160 independent samples. Compared with a reference classification established from Affymetrix U133 + 2 data, reverse transcriptase multiplex ligation-dependent probe amplification classified 85.0% samples into the expected subtype, comparing favorably with current diagnostic methods. This assay also proved to be highly efficient in detecting the MYD88 L265P mutation, even in archival paraffin-embedded tissues. This reliable, rapid, and cost-effective method uses common instruments and reagents and could thus easily be implemented into routine diagnosis workflows, to improve the management of these aggressive tumors.

Identification of Somatic Mutations in Primary Cutaneous Diffuse Large B-Cell Lymphoma, Leg Type by Massive Parallel Sequencing

To determine whether the mutational profile of primary cutaneous diffuse large B-cell lymphoma, leg type (PCLBCL-LT) is unique by comparison with other diffuse large B-cell lymphoma subtypes, we analyzed a total cohort of 20 PCLBCL-LT patients by using next-generation sequencing with a lymphoma panel designed for diffuse large B-cell lymphoma. We also analyzed 12 pairs of tumor and control DNA samples by whole-exome sequencing, which led us to perform resequencing of three selected genes not included in the lymphoma panel: TBL1XR1, KLHL6, and IKZF3. Our study clearly identifies an original mutational landscape of PCLBCL-LT with a very restricted set of highly recurrent mutations (>40%) involving MYD88 (p.L265P variant), PIM1, and CD79B. Other genes involved in B-cell signaling, NF-κB activation, or DNA modeling were found altered, notably TBL1XR1 (33%), MYC (26%) CREBBP (26%), and IRF4 (21%) or HIST1H1E (41%). MYD88^L265P variant was associated with copy number variations or copy neutral loss of heterozygosity in 60% of patients. The most frequent genetic losses involved CDKN2A/2B, TNFAIP3/A20, PRDM1, TCF3, and CIITA. Together, these results show that PCLBCL-LT exhibits a unique mutational landscape, combining highly recurrent hotspot mutations in genes involved in NF-kB and B-cell signaling pathways, which provides a rationale for using selective inhibitors of the B-cell receptor.

Oncogenic events rather than antigen selection pressure may be the main driving forces for relapse in diffuse large B-cell lymphomas

Little is known on the phylogenetic relationship between diagnostic and relapse clones of diffuse large B-cell lymphoma (DLBCL). We applied high throughput sequencing (HTS) of the VDJ locus of Immunoglobulin heavy chain (IGHV) on 14 DLBCL patients with serial samples, including tumor biopsies and/or peripheral blood mononuclear cells (PBMC). Phylogenetic data were consolidated with targeted sequencing and cytogenetics. Phylogeny clearly showed that DLBCL relapse could occur according either an early or a late divergent mode. These two modes of divergence were independent from the elapsed time between diagnosis and relapse. We found no significant features for antigen selection pressure in complementary determining region both at diagnosis and relapse for 9/12 pairs and a conserved negative selection pressure for the three remaining cases. Targeted HTS and conventional cytogenetics revealed a branched vs. linear evolution for 5/5 IGHV early divergent cases, but unexpected such "oncogenetic" branched evolution could be found in at least 2/7 IGHV late divergent cases. Thus, if BCR signaling is mandatory for DLBCL emergence, oncogenetic events under chemotherapy selection pressure may be the main driving forces at relapse. Finally, circulating subclones with divergent IGHV somatic hypermutations patterns from initial biopsy could be detected in PBMC at diagnosis for 4/6 patients and, for two of them, at least one was similar to the ones found at relapse. This study highlights that oncogenetic intraclonal diversity of DLBCL should be evaluated beyond the scope a single biopsy and represents a rationale for future investigations using peripheral blood for lymphoid malignancies genotyping. Am. J. Hematol. 92:68-76, 2017. © 2016 Wiley Periodicals, Inc.

Biological and Clinical Relevance of Associated Genomic Alterations in MYD88 L265P and non-L265P-Mutated Diffuse Large B-Cell Lymphoma: Analysis of 361 Cases

Purpose:MYD88 mutations, notably the recurrent gain-of-function L265P variant, are a distinguishing feature of activated B-cell like (ABC) diffuse large B-cell lymphoma (DLBCL), leading to constitutive NFκB pathway activation. The aim of this study was to examine the distinct genomic profiles of MYD88-mutant DLBCL, notably according to the presence of the L265P or other non-L265P MYD88 variants.Experimental Design: A cohort of 361 DLBCL cases (94 MYD88 mutant and 267 MYD88 wild-type) was submitted to next-generation sequencing (NGS) focusing on 34 genes to analyze associated mutations and copy number variations, as well as gene expression profiling, and clinical and prognostic analyses.Results: Importantly, we highlighted different genomic profiles for MYD88 L265P and MYD88 non-L265P-mutant DLBCL, shedding light on their divergent backgrounds. Clustering analysis also segregated subgroups according to associated genetic alterations among patients with the same MYD88 mutation. We showed that associated CD79B and MYD88 L265P mutations act synergistically to increase NFκB pathway activation, although the majority of MYD88 L265P-mutant cases harbors downstream NFκB alterations, which can predict BTK inhibitor resistance. Finally, although the MYD88 L265P variant was not an independent prognostic factor in ABC DLBCL, associated CD79B mutations significantly improved the survival of MYD88 L265P-mutant ABC DLBCL in our cohort.Conclusions: This study highlights the relative heterogeneity of MYD88-mutant DLBCL, adding to the field's knowledge of the theranostic importance of MYD88 mutations, but also of associated alterations, emphasizing the usefulness of genomic profiling to best stratify patients for targeted therapy. Clin Cancer Res; 23(9); 2232-44. ©2016 AACR.

Digital PCR for quantification of recurrent and potentially actionable somatic mutations in circulating free DNA from patients with diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is an aggressive and heterogeneous malignancy harboring frequent targetable activating somatic mutations. Emerging evidence suggests that circulating cell-free DNA (cfDNA) can be used to detect somatic variants in DLBCL using Next-Generation Sequencing (NGS) experiments. In this proof-of-concept study, we chose to develop simple and valuable digital PCR (dPCR) assays for the detection of recurrent exportin-1 (XPO1) E571K, EZH2 Y641N, and MYD88 L265P mutations in DLBCL patients, thereby identifying patients most likely to potentially benefit from targeted therapies. We demonstrated that our dPCR assays were sufficiently sensitive to detect rare XPO1, EZH2, and MYD88 mutations in plasma cfDNA, with a sensitivity of 0.05%. cfDNA somatic mutation detection by dPCR seems to be a promising technique in the management of DLBCL, in addition to NGS experiments.

Detection and prognostic value of recurrent exportin 1 mutations in tumor and cell-free circulating DNA of patients with classical Hodgkin lymphoma

Classical Hodgkin lymphoma is one of the most common lymphomas and shares clinical and genetic features with primary mediastinal B-cell lymphoma. In this retrospective study, we analyzed the recurrent hotspot mutation of the exportin 1 (XPO1, p.E571K) gene, previously identified in primary mediastinal B-cell lymphoma, in biopsies and plasma circulating cell-free DNA from patients with classical Hodgkin lymphoma using a highly sensitive digital PCR technique. A total of 94 patients were included in the present study. This widely expressed XPO1 E571K mutation is present in one quarter of classical Hodgkin lymphoma patients (24.2%). Mutated and wild-type classical Hodgkin lymphomas were similar regarding the main clinical features. Patients with a detectable XPO1 mutation at the end of treatment displayed a tendency toward shorter progression-free survival, as compared to patients with undetectable mutation in plasma cell-free DNA (2-year progression-free survival: 57.1%, 95% confidence interval: 30.1-100% versus 2-year progression-free survival: 90.5%, 95% confidence interval: 78.8-100%, respectively, P=0.0601). To conclude, the detection of the XPO1 E571K mutation in biopsy and plasma cell-free DNA by digital PCR may be used as a novel biomarker in classical Hodgkin lymphoma for both diagnosis and minimal residual disease, and pinpoints a crucial role of XPO1 in classical Hodgkin lymphoma pathogenesis. The detection of somatic mutation in the plasma cell-free DNA of patients represents a major technological advance in the context of liquid biopsies and noninvasive management of classical Hodgkin lymphoma.

Next-Generation Sequencing in Diffuse Large B-Cell Lymphoma Highlights Molecular Divergence and Therapeutic Opportunities: a LYSA Study

PURPOSE

Next-generation sequencing (NGS) has detailed the genomic characterization of diffuse large B-cell lymphoma (DLBCL) by identifying recurrent somatic mutations. We set out to design a clinically feasible NGS panel focusing on genes whose mutations hold potential therapeutic impact. Furthermore, for the first time, we evaluated the prognostic value of these mutations in prospective clinical trials.

EXPERIMENTAL DESIGN

A Lymphopanel was designed to identify mutations in 34 genes, selected according to literature and a whole exome sequencing study of relapsed/refractory DLBCL patients. The tumor DNA of 215 patients with CD20(+)de novo DLBCL in the prospective, multicenter, and randomized LNH-03B LYSA clinical trials was sequenced to deep, uniform coverage with the Lymphopanel. Cell-of-origin molecular classification was obtained through gene expression profiling with HGU133+2.0 Affymetrix GeneChip arrays.

RESULTS

The Lymphopanel was informative for 96% of patients. A clear depiction of DLBCL subtype molecular heterogeneity was uncovered with the Lymphopanel, confirming that activated B-cell-like (ABC), germinal center B-cell like (GCB), and primary mediastinal B-cell lymphoma (PMBL) are frequently affected by mutations in NF-κB, epigenetic, and JAK-STAT pathways, respectively. Novel truncating immunity pathway, ITPKB, MFHAS1, and XPO1 mutations were identified as highly enriched in PMBL. Notably, TNFAIP3 and GNA13 mutations in ABC patients treated with R-CHOP were associated with significantly less favorable prognoses.

CONCLUSIONS

This study demonstrates the contribution of NGS with a consensus gene panel to personalized therapy in DLBCL, highlighting the molecular heterogeneity of subtypes and identifying somatic mutations with therapeutic and prognostic impact. Clin Cancer Res; 22(12); 2919-28. ©2016 AACRSee related commentary by Lim and Elenitoba-Johnson, p. 2829.

Immunohistochemical and genomic profiles of diffuse large B-cell lymphomas: implications for targeted EZH2 inhibitor therapy?

Enhancer of Zeste Homolog 2 (EZH2) plays an essential epigenetic role in Diffuse Large B Cell Lymphoma (DLBCL) development. Recurrent somatic heterozygous gain-of-function mutations of EZH2 have been identified in DLBCL, most notably affecting tyrosine 641 (Y641), inducing hyper-trimethylation of H3K27 (H3K27me3). Novel EZH2 inhibitors are being tested in phase 1 and 2 clinical trials but no study has examined which patients would most benefit from this treatment. We evaluated the immunohistochemical (IHC) methylation profiles of 82 patients with DLBCL, as well as the mutational profiles of 32 patients with DLBCL using NGS analysis of a panel of 34 genes involved in lymphomagenesis. A novel IHC score based on H3K27me2 and H3K27me3 expression was developed, capable of distinguishing patients with wild-type (WT) EZH2 and patients with EZH2 Y641 mutations (p = 10⁻⁵). NGS analysis revealed a subclonal EZH2 mutation pattern in EZH2 mutant patients with WT-like IHC methylation profiles, while associated mutations capable of upregulating EZH2 were detected in WT EZH2 patients with mutant-like IHC methylation profiles. IHC and mutational profiles highlight in vivo hyper-H3K27me3 and hypo-H3K27me2 status, pinpoint associated activating mutations and determine EZH2 mutation clonality, maximizing EZH2 inhibitor potential by identifying patients most likely to benefit from treatment.