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.

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.

5' Rapid amplification of cDNA ends (5'RACE): A simpler method to analyze immunoglobulin genes and discover the value of the light chain in chronic lymphocytic leukemia.

The mutational status of the variable region of the immunoglobulin heavy chain gene (IGHV) is a very important biomarker for chronic lymphocytic leukemia (CLL) patients. However, the routine detection of IGHV mutational status is time-consuming and costly. Therefore, we performed 5' Rapid amplification of cDNA ends (5' RACE) in 81 CLL patients who previously underwent detection using Biomed-2. The agreement rate of these two methods was 93.8 %. Regarding the discordant cases, 5' RACE was more sensitive to identify unproductive and multiple rearrangements. Furthermore, 5' RACE can also be used to simultaneously sequence light chains. In most CLL cases, the mutational statuses of heavy and light chains are concordant, except in IGLV3-21. Most IGLV3-21 (24/25) rearrangement shared a similar LCDR3 (QVWDSSSDHPWV) and harbored a single point mutation, namely, IGLV3-21R110. Compared to mutated-CLL non IGLV3-2R110, IGLV3-21R110-CLL exhibited a shorter overall survival (OS) and time to first treatment (TTFT) (p = 0.05, p < 0.0001, respectively) even though 75 % (18/24) of these patients expressed mutated heavy chains. Altogether, IGLV3-21R110 defines a CLL subgroup with specific biological features and an unfavorable prognosis independent of the IGHV mutational status and emphasizes the important value of the light chain. This study is the first to use 5' RACE to detect the mutational status of IGH in CLL. Here, 5' RACE was a reliable and effective method to test the mutational status of heavy and light chains. In addition, 5' RACE can be combined with other assays in the NGS workflow to obtain more detailed insight into subclonal architecture and intraclonal diversity.

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.

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.

Circulating tumor DNA in primary mediastinal large B-cell lymphoma versus classical Hodgkin lymphoma: a retrospective study.

Few data exist concerning circulating tumor DNA (ctDNA) relevance in primary mediastinal B-cell lymphoma (PMBL). To explore this topic, we applied a 9-gene next-generation sequencing pipeline to samples from forty-four PMBL patients (median age 36.5 years). The primary endpoint was a similarity between paired biopsy/plasma mutational profiles. We detected at least one variant in 32 plasma samples (80%). The similarity between the biopsy and ctDNA genetic profiles for the 30 patients with paired mutated biopsy/plasma samples was greater than or equal to 80% in 19 patients (63.3%). We then compared PMBL ctDNA features with those of a cohort of Hodgkin lymphoma patients (n = 60). The top three mutated genes were SOCS1, TNFAIP3, and B2M in both lymphoma types. PMBL displayed more alterations in TNFAIP3 (71.9% vs. 46.3%, p = 0.029) and GNA13 (46.9% vs. 17.1%, p = 0.013) than cHL. Our 9-gene set may delineate tumor genotypes using ctDNA samples from both lymphoma types.

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.

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.

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.

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.

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

BACKGROUND: 18F-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.

Combining gene expression profiling and machine learning to diagnose B-cell non-Hodgkin lymphoma.

Non-Hodgkin B-cell lymphomas (B-NHLs) are a highly heterogeneous group of mature B-cell malignancies. Their classification thus requires skillful evaluation by expert hematopathologists, but the risk of error remains higher in these tumors than in many other areas of pathology. To facilitate diagnosis, we have thus developed a gene expression assay able to discriminate the seven most frequent B-cell NHL categories. This assay relies on the combination of ligation-dependent RT-PCR and next-generation sequencing, and addresses the expression of more than 130 genetic markers. It was designed to retrieve the main gene expression signatures of B-NHL cells and their microenvironment. The classification is handled by a random forest algorithm which we trained and validated on a large cohort of more than 400 annotated cases of different histology. Its clinical relevance was verified through its capacity to prevent important misclassification in low grade lymphomas and to retrieve clinically important characteristics in high grade lymphomas including the cell-of-origin signatures and the MYC and BCL2 expression levels. This accurate pan-B-NHL predictor, which allows a systematic evaluation of numerous diagnostic and prognostic markers, could thus be proposed as a complement to conventional histology to guide the management of patients and facilitate their stratification into clinical trials.

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.

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.