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.

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.

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.

Integrated analyses of translatome and proteome identify the rules of translation selectivity in RPS14-deficient cells.

In ribosomopathies, the Diamond-Blackfan anemia (DBA) or 5q- syndrome, ribosomal protein (RP) genes are affected by mutation or deletion, resulting in bone marrow erythroid hypoplasia. Unbalanced production of ribosomal subunits leading to a limited ribosome cellular content regulates translation at the expense of the master erythroid transcription factor GATA1. In RPS14-deficient cells mimicking 5q- syndrome erythroid defects, we show that the transcript length, codon bias of the coding sequence (CDS) and 3'UTR (untranslated region) structure are the key determinants of translation. In these cells, short transcripts with a structured 3'UTR and high codon adaptation index (CAI) showed a decreased translation efficiency. Quantitative analysis of the whole proteome confirmed that the post-transcriptional changes depended on the transcript characteristics that governed the translation efficiency in conditions of low ribosome availability. In addition, proteins involved in normal erythroid differentiation share most determinants of translation selectivity. Our findings thus indicate that impaired erythroid maturation due to 5q- syndrome may proceed from a translational selectivity at the expense of the erythroid differentiation program, and suggest that an interplay between the CDS and UTR may regulate mRNA translation.

Targeted genotyping of circulating tumor DNA for classical Hodgkin lymphoma monitoring: a prospective study.

The relevance of circulating tumor DNA (ctDNA) analysis as a liquid biopsy and minimal residual disease tool in the management of classical Hodgkin Lymphoma (cHL) patients was demonstrated in retrospective settings and remains to be confirmed in a prospective setting. We developed a targeted Next-Generation sequencing (NGS) panel for fast analysis (AmpliSeq technology) of nine commonly mutated genes in biopies and ctDNA of cHL patients. We then conducted a prospective trial to assess ctDNA follow up at diagnosis and after 2 cycles of chemotherapy (C2). Sixty cHL patients treated by first line conventional chemotherapy (BEACOPPescalated [21.3%], ABVD/ABVD-like [73.5%] and other regimens [5.2%, for elderly patients] were assessed in this non-interventional study. Median age of the patients was 33.5 years (range 20-86). Variants were identified in 42 (70%) patients. Mutations of NFKBIE, TNFAIP3, STAT6, PTPN1, B2M, XPO1, ITPKB, GNA13 and SOCS1 were found in 13.3%, 31.7%, 23.3%, 5%, 33.3%, 10%, 23.3%, 13.3% and 50% of patients, respectively. ctDNA concentration and genotype are correlated with clinical characteristics and presentation. Regarding early therapeutic response, 45 patients (83%, NA=6) had a negative positron emission tomography (PET) after C2 (Deauville Score 1-3). Mean of DeltaSUVmax after C2 was -78.8%. We analyzed ctDNA after C2 for 54 patients (90%). ctDNA became rapidly undetectable in all cases after C2. Variant detection in ctDNA is suitable to depict the genetic features of cHL at diagnosis and may help to assess early treatment response, in association with PET. Clinical Trial reference: NCT02815137.

Plasmacytoid dendritic cells proliferation associated with acute myeloid leukemia: phenotype profile and mutation landscape

Neoplasms involving plasmacytoid Dendritic Cells (pDCs) include Blastic pDC Neoplasms (BPDCN) and other pDC proliferations, where pDCs are associated with myeloid malignancies: most frequently Chronic MyeloMonocytic Leukemia (CMML) but also Acute Myeloid Leukemia (AML), hereafter named pDC-AML. We aimed to determine the reactive or neoplastic origin of pDCs in pDC-AML, and their link with the CD34+ blasts, monocytes or conventional DCs (cDCs) associated in the same sample, by phenotypic and molecular analyses (targeted NGS, 70 genes). We compared 15 pDC-AML at diagnosis with 21 BPDCN and 11 normal pDCs from healthy donors. CD45low CD34+ blasts were found in all cases (10-80% of medullar cells), associated with pDCs (4-36%), monocytes in 14 cases (1-10%) and cDCs (2 cases, 4.8-19%). pDCs in pDC-AML harbor a clearly different phenotype from BPDCN: CD4+ CD56- in 100% of cases, most frequently CD303+, CD304+ and CD34+; lower expression of cTCL1 and CD123 with isolated lymphoid markers (CD22/CD7/CD5) in some cases, suggesting a pre-pDC stage. In all cases, pDCs, monocytes and cDC are neoplastic since they harbor the same mutations as CD34+ blasts. RUNX1 is the most commonly mutated gene: detected in all AML with minimal differentiation (M0-AML) but not in the other cases. Despite low number of cases, the systematic association between M0-AML, RUNX1 mutations and an excess of pDC is puzzling. Further evaluation in a larger cohort is required to confirm RUNX1 mutations in pDC-AML with minimal differentiation and to investigate whether it represents a proliferation of blasts with macrophage and DC progenitor potential.

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 TFH-derived lymphomas divided into two subgroups according to a predominant TFH (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 TFH, 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.

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.

A lowered 26S proteasome activity correlates with mantle lymphoma cell lines resistance to genotoxic stress.

BACKGROUND: Mantle cell lymphoma (MCL) is a B-cell hemopathy characterized by the t(11;14) translocation and the aberrant overexpression of cyclin D1. This results in an unrestrained cell proliferation. Other genetic alterations are common in MCL cells such as SOX11 expression, mutations of ATM and/or TP53 genes, activation of the NF-κB signaling pathway and NOTCH receptors. These alterations lead to the deregulation of the apoptotic machinery and resistance to drugs. We observed that among a panel of MCL cell lines, REC1 cells were resistant towards genotoxic stress. We studied the molecular basis of this resistance. METHODS: We analyzed the cell response regarding apoptosis, senescence, cell cycle arrest, DNA damage response and finally the 26S proteasome activity following a genotoxic treatment that causes double strand DNA breaks. RESULTS: MCL cell lines displayed various sensitivity/resistance towards genotoxic stress and, in particular, REC1 cells did not enter apoptosis or senescence after an etoposide treatment. Moreover, the G2/M cell cycle checkpoint was deficient in REC1 cells. We observed that three main actors of apoptosis, senescence and cell cycle regulation (cyclin D1, MCL1 and CDC25A) failed to be degraded by the proteasome machinery in REC1 cells. We ruled out a default of the ßTrCP E3-ubiquitine ligase but detected a lowered 26S proteasome activity in REC1 cells compared to other cell lines. CONCLUSION: The resistance of MCL cells to genotoxic stress correlates with a low 26S proteasome activity. This could represent a relevant biomarker for a subtype of MCL patients with a poor response to therapies and a high risk of relapse.

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.

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

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

Detection 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.

Recurrent mutations of the exportin 1 gene (XPO1) and their impact on selective inhibitor of nuclear export compounds sensitivity in primary mediastinal B-cell lymphoma.

Primary mediastinal B-cell lymphoma (PMBL) is an entity of B-cell lymphoma distinct from the other molecular subtypes of diffuse large B-cell lymphoma (DLBCL). We investigated the prevalence, specificity, and clinical relevance of mutations of XPO1, which encodes a member of the karyopherin-ß nuclear transporters, in a large cohort of PMBL. PMBL cases defined histologically or by gene expression profiling (GEP) were sequenced and the XPO1 mutational status was correlated to genetic and clinical characteristics. The XPO1 mutational status was also assessed in DLBCL, Hodgkin lymphoma (HL) and mediastinal gray-zone lymphoma (MGZL).The biological impact of the mutation on Selective Inhibitor of Nuclear Export (SINE) compounds (KPT-185/330) sensitivity was investigated in vitro. XPO1 mutations were present in 28/117 (24%) PMBL cases and in 5/19 (26%) HL cases but absent/rare in MGZL (0/20) or DLBCL (3/197). A higher prevalence (50%) of the recurrent codon 571 variant (p.E571K) was observed in GEP-defined PMBL and was associated with shorter PFS. Age, International Prognostic Index and bulky mass were similar in XPO1 mutant and wild-type cases. KPT-185 induced a dose-dependent decrease in cell proliferation and increased cell-death in PMBL cell lines harboring wild type or XPO1 E571K mutant alleles. Experiments in transfected U2OS cells further confirmed that the XPO1 E571K mutation does not have a drastic impact on KPT-330 binding. To conclude the XPO1 E571K mutation represents a genetic hallmark of the PMBL subtype and serves as a new relevant PMBL biomarker. SINE compounds appear active for both mutated and wild-type protein. Am. J. Hematol. 91:923-930, 2016. © 2016 Wiley Periodicals, Inc.

Real life evaluation of AlphaMissense predictions in hematological malignancies.

High-throughput sequencing plays a pivotal role in hematological malignancy diagnostics, but interpreting missense mutations remains challenging. In this study, we used the newly available AlphaMissense database to assess the efficacy of machine learning to predict missense mutation effects and its impact to improve our ability to interpret them. Based on the analysis of 2073 variants from 686 patients analyzed for clinical purpose, we confirmed the very high accuracy of AlphaMissense predictions in a large real-life data set of missense mutations (AUC of ROC curve 0.95), and provided a comprehensive analysis of the discrepancies between AlphaMissense predictions and state of the art clinical interpretation.

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.

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.

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.