Analysis of mutational signatures in exomes from B-cell lymphoma cell lines suggest APOBEC3 family members to be involved in the pathogenesis of primary effusion lymphoma

Update: The molecular spectrum of virus-associated high-grade B-cell non-Hodgkin lymphomas

The vast spectrum of aggressive B-cell non-Hodgkin neoplasms (B-NHL) encompasses several infrequent entities occurring in association with viral infections, posing diagnostic challenges for practitioners. In the emerging era of precision oncology, the molecular characterization of malignancies has acquired paramount significance. The pathophysiological comprehension of specific entities and the identification of targeted therapeutic options have seen rapid development. However, owing to their rarity, not all entities have undergone exhaustive molecular characterization. Considerable heterogeneity exists in the extant body of work, both in terms of employed methodologies and the scale of cases studied. Presently, therapeutic strategies are predominantly derived from observations in diffuse large B-cell lymphoma (DLBCL), the most prevalent subset of aggressive B-NHL. Ongoing investigations into the molecular profiles of these uncommon virus-associated entities are progressively facilitating a clearer distinction from DLBCL, ultimately paving the way towards individualized therapeutic approaches. This review consolidates the current molecular insights into aggressive and virus-associated B-NHL, taking into consideration the recently updated 5th edition of the WHO classification of hematolymphoid tumors (WHO-5HAEM) and the International Consensus Classification (ICC). Additionally, potential therapeutically targetable susceptibilities are highlighted, offering a comprehensive overview of the present scientific landscape in the field.

Human APOBEC3B promotes tumor development in vivo including signature mutations and metastases

The antiviral DNA cytosine deaminase APOBEC3B has been implicated as a source of mutation in many cancers. However, despite years of work, a causal relationship has yet to be established in vivo. Here, we report a murine model that expresses tumor-like levels of human APOBEC3B. Animals expressing full-body APOBEC3B appear to develop normally. However, adult males manifest infertility, and older animals of both sexes show accelerated rates of carcinogenesis, visual and molecular tumor heterogeneity, and metastasis. Both primary and metastatic tumors exhibit increased frequencies of C-to-T mutations in TC dinucleotide motifs consistent with the established biochemical activity of APOBEC3B. Enrichment for APOBEC3B-attributable single base substitution mutations also associates with elevated levels of insertion-deletion mutations and structural variations. APOBEC3B catalytic activity is required for all of these phenotypes. Together, these studies provide a cause-and-effect demonstration that human APOBEC3B is capable of driving both tumor initiation and evolution in vivo.

APOBEC3-mediated mutagenesis in cancer: causes, clinical significance and therapeutic potential

Apolipoprotein B mRNA-editing enzyme, catalytic polypeptides (APOBECs) are cytosine deaminases involved in innate and adaptive immunity. However, some APOBEC family members can also deaminate host genomes to generate oncogenic mutations. The resulting mutations, primarily signatures 2 and 13, occur in many tumor types and are among the most common mutational signatures in cancer. This review summarizes the current evidence implicating APOBEC3s as major mutators and outlines the exogenous and endogenous triggers of APOBEC3 expression and mutational activity. The review also discusses how APOBEC3-mediated mutagenesis impacts tumor evolution through both mutagenic and non-mutagenic pathways, including by inducing driver mutations and modulating the tumor immune microenvironment. Moving from molecular biology to clinical outcomes, the review concludes by summarizing the divergent prognostic significance of APOBEC3s across cancer types and their therapeutic potential in the current and future clinical landscapes.

Human APOBEC3B promotes tumor heterogeneity in vivo including signature mutations and metastases

The antiviral DNA cytosine deaminase APOBEC3B has been implicated as a source of mutation in many different cancers. Despite over 10 years of work, a causal relationship has yet to be established between APOBEC3B and any stage of carcinogenesis. Here we report a murine model that expresses tumor-like levels of human APOBEC3B after Cre-mediated recombination. Animals appear to develop normally with full-body expression of APOBEC3B. However, adult males manifest infertility and older animals of both sexes show accelerated rates of tumorigenesis (mostly lymphomas or hepatocellular carcinomas). Interestingly, primary tumors also show overt heterogeneity, and a subset spreads to secondary sites. Both primary and metastatic tumors exhibit increased frequencies of C-to-T mutations in TC dinucleotide motifs consistent with the established biochemical activity of APOBEC3B. Elevated levels of structural variation and insertion-deletion mutations also accumulate in these tumors. Together, these studies provide the first cause-and-effect demonstration that human APOBEC3B is an oncoprotein capable of causing a wide range of genetic changes and driving tumor formation in vivo .

Whole-genome sequencing reveals complex genomic features underlying anti-CD19 CAR T-cell treatment failures in lymphoma

CD19-directed chimeric antigen receptor (CAR-19) T cells are groundbreaking immunotherapies approved for use against large B-cell lymphomas. Although host inflammatory and tumor microenvironmental markers associate with efficacy and resistance, the tumor-intrinsic alterations underlying these phenomena remain undefined. CD19 mutations associate with resistance but are uncommon, and most patients with relapsed disease retain expression of the wild-type receptor, implicating other genomic mechanisms. We therefore leveraged the comprehensive resolution of whole-genome sequencing to assess 51 tumor samples from 49 patients with CAR-19-treated large B-cell lymphoma. We found that the pretreatment presence of complex structural variants, APOBEC mutational signatures, and genomic damage from reactive oxygen species predict CAR-19 resistance. In addition, the recurrent 3p21.31 chromosomal deletion containing the RHOA tumor suppressor was strongly enriched in patients for whom CAR T-cell therapy failed. Pretreatment reduced expression or monoallelic loss of CD19 did not affect responses, suggesting CAR-19 therapy success and resistance are related to multiple mechanisms. Our study showed that tumor-intrinsic genomic alterations are key among the complex interplay of factors that underlie CAR-19 efficacy and resistance for large B-cell lymphomas.

A practical framework RNMF for exploring the association between mutational signatures and genes using gene cumulative contribution abundance

Background

Mutational signatures are somatic mutation patterns enriching operational mutational processes, which can provide abundant information about the mechanism of cancer. However, understanding of the pathogenic biological processes is still limited, such as the association between mutational signatures and genes.

Methods

We developed a simple and practical R package called RNMF (https://github.com/zhenzhang‐li/RNMF) for mutational signature analysis, including a key model of cumulative contribution abundance (CCA), which was designed to highlight the association between mutational signatures and genes and then applying it to a meta‐analysis of 1073 individuals with esophageal squamous cell carcinoma (ESCC).

Results

We revealed a number of known and previously undescribed SBS or ID signatures, and we found that APOBEC signatures (SBS2* and SBS13*) were closely associated with PIK3CA mutation, especially the E545k mutation. Furthermore, we found that age signature is closely related to the frequent mutation of TP53, of which R342* is highlighted due to strongly linked to age signature. In addition, the CCA matrix image data of genes in the signatures New, SBS3*, and SBS17b* were helpful for the preliminary evaluation of shortened survival outcome. These results can be extended to estimate the distribution of mutations or features, and study the potential impact of clinical factors.

Conclusions

In a word, RNMF can successfully achieve the correlation analysis of mutational signatures and genes, proving a strong theoretical basis for the study of mutational processes during tumor development.

KSHV/HHV8-mediated hematologic diseases

Kaposi sarcoma (KS) herpesvirus (KSHV), also known as human herpesvirus 8, is the causal agent of KS but is also pathogenetically related to several lymphoproliferative disorders, including primary effusion lymphoma (PEL)/extracavitary (EC) PEL, KSHV-associated multicentric Castleman disease (MCD), KSHV+ diffuse large B-cell lymphoma, and germinotropic lymphoproliferative disorder. These different KSHV-associated diseases may co-occur and may have overlapping features. KSHV, similar to Epstein-Barr virus (EBV), is a lymphotropic gammaherpesvirus that is preferentially present in abnormal lymphoid proliferations occurring in immunecompromised individuals. Notably, both KSHV and EBV can infect and transform the same B cell, which is frequently seen in KSHV+ EBV+ PEL/EC-PEL. The mechanisms by which KSHV leads to lymphoproliferative disorders is thought to be related to the expression of a few transforming viral genes that can affect cellular proliferation and survival. There are critical differences between KSHV-MCD and PEL/EC-PEL, the 2 most common KSHV-associated lymphoid proliferations, including viral associations, patterns of viral gene expression, and cellular differentiation stage reflected by the phenotype and genotype of the infected abnormal B cells. Advances in treatment have improved outcomes, but mortality rates remain high. Our deepening understanding of KSHV biology, clinical features of KSHV-associated diseases, and newer clinical interventions should lead to improved and increasingly targeted therapeutic interventions.

KRAS and RAS-MAPK Pathway Deregulation in Mature B Cell Lymphoproliferative Disorders

KRAS mutations account for the most frequent mutations in human cancers, and are generally correlated with disease aggressiveness, poor prognosis, and poor response to therapies. KRAS is required for adult hematopoiesis and plays a key role in B cell development and mature B cell proliferation and survival, proved to be critical for B cell receptor-induced ERK pathway activation. In mature B cell neoplasms, commonly seen in adults, KRAS and RAS-MAPK pathway aberrations occur in a relevant fraction of patients, reaching high recurrence in some specific subtypes like multiple myeloma and hairy cell leukemia. As inhibitors targeting the RAS-MAPK pathway are being developed and improved, it is of outmost importance to precisely identify all subgroups of patients that could potentially benefit from their use. Herein, we review the role of KRAS and RAS-MAPK signaling in malignant hematopoiesis, focusing on mature B cell lymphoproliferative disorders. We discuss KRAS and RAS-MAPK pathway aberrations describing type, incidence, mutual exclusion with other genetic abnormalities, and association with prognosis. We review the current therapeutic strategies applied in mature B cell neoplasms to counteract RAS-MAPK signaling in pre-clinical and clinical studies, including most promising combination therapies. We finally present an overview of genetically engineered mouse models bearing KRAS and RAS-MAPK pathway aberrations in the hematopoietic compartment, which are valuable tools in the understanding of cancer biology and etiology.

Long interspersed element-1 ribonucleoprotein particles protect telomeric ends in alternative lengthening of telomeres dependent cells

Malignant cells ensure telomere maintenance by the alternative lengthening of telomeres (ALT) in the absence of telomerase activity (TA). The retrotransposons "long interspersed nuclear element-1" (LINE-1, L1) are expressed in malignant cells and are primarily known to contribute to complex karyotypes. Here we demonstrate that LINE-1 ribonucleoprotein particles (L1-RNPs) expression is significantly higher in ALT⁺- versus in TA⁺-human glioma. Analyzing a role of L1-RNP in ALT, we show that L1-RNPs bind to telomeric repeat containing RNA (TERRA), which is critical for telomere stabilization and which is overexpressed in ALT⁺ cells. In turn, L1-RNP knockdown (KD) abrogated the nuclear retention of TERRA, resulted in increased telomeric DNA damage, decreased cell growth and reduced expression of ALT characteristics such as c-circles and PML-bodies. L1-RNP KD also decreased the expression of Shelterin- and the ALT-regulating protein Topoisomerase IIIα (TopoIIIα) indicating a more general role of L1-RNPs in supporting telomeric integrity in ALT. Our findings suggest an impact of L1-RNP on telomere stability in ALT⁺ dependent tumor cells. As L1-RNPs are rarely expressed in normal adult human tissue those elements might serve as a novel target for tumor ablative therapy.

Discovering novel mutation signatures by latent Dirichlet allocation with variational Bayes inference

MOTIVATION

A cancer genome includes many mutations derived from various mutagens and mutational processes, leading to specific mutation patterns. It is known that each mutational process leads to characteristic mutations, and when a mutational process has preferences for mutations, this situation is called a 'mutation signature.' Identification of mutation signatures is an important task for elucidation of carcinogenic mechanisms. In previous studies, analyses with statistical approaches (e.g. non-negative matrix factorization and latent Dirichlet allocation) revealed a number of mutation signatures. Nonetheless, strictly speaking, these existing approaches employ an ad hoc method or incorrect approximation to estimate the number of mutation signatures, and the whole picture of mutation signatures is unclear.

RESULTS

In this study, we present a novel method for estimating the number of mutation signatures-latent Dirichlet allocation with variational Bayes inference (VB-LDA)-where variational lower bounds are utilized for finding a plausible number of mutation patterns. In addition, we performed cluster analyses for estimated mutation signatures to extract novel mutation signatures that appear in multiple primary lesions. In a simulation with artificial data, we confirmed that our method estimated the correct number of mutation signatures. Furthermore, applying our method in combination with clustering procedures for real mutation data revealed many interesting mutation signatures that have not been previously reported.

AVAILABILITY AND IMPLEMENTATION

All the predicted mutation signatures with clustering results are freely available at http://www.f.waseda.jp/mhamada/MS/index.html. All the C++ source code and python scripts utilized in this study can be downloaded on the Internet (https://github.com/qkirikigaku/MS_LDA).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Determinants of Oligonucleotide Selectivity of APOBEC3B

APOBEC3B (A3B) is a prominent source of mutation in many cancers. To date, it has been difficult to capture the native protein-DNA interactions that confer A3B's substrate specificity by crystallography due to the highly dynamic nature of wild-type A3B active site. We use computational tools to restore a recent crystal structure of a DNA-bound A3B C-terminal domain mutant construct to its wild type sequence, and run molecular dynamics simulations to study its substrate recognition mechanisms. Analysis of these simulations reveal dynamics of the native A3Bctd-oligonucleotide interactions, including the experimentally inaccessible loop 1-oligonucleotide interactions. A second series of simulations in which the target cytosine nucleotide was computationally mutated from a deoxyribose to a ribose show a change in sugar ring pucker, leading to a rearrangement of the binding site and revealing a potential intermediate in the binding pathway. Finally, apo simulations of A3B, starting from the DNA-bound open state, experience a rapid and consistent closure of the binding site, reaching conformations incompatible with substrate binding. This study reveals a more realistic and dynamic view of the wild type A3B binding site and provides novel insights for structure-guided design efforts for A3B.

Biology and management of primary effusion lymphoma

Primary effusion lymphoma (PEL) is a rare B-cell malignancy that most often occurs in immunocompromised patients, such as HIV-infected individuals and patients receiving organ transplantation. The main characteristic of PEL is neoplastic effusions in body cavities without detectable tumor masses. The onset of the disease is associated with latent infection of human herpes virus 8/Kaposi sarcoma–associated herpes virus, and the normal counterpart of tumor cells is B cells with plasmablastic differentiation. A condition of immunodeficiency and a usual absence of CD20 expression lead to the expectation of the lack of efficacy of anti-CD20 monoclonal antibody; clinical outcomes of the disease remain extremely poor, with an overall survival at 1 year of ∼30%. Although recent progress in antiretroviral therapy has improved outcomes of HIV-infected patients, its benefit is still limited in patients with PEL. Furthermore, the usual high expression of programmed death ligand 1 in tumor cells, one of the most important immune-checkpoint molecules, results in the immune escape of tumor cells from the host immune defense, which could be the underlying mechanism of poor treatment efficacy. Molecular-targeted therapies for the activating pathways in PEL, including NF-κB, JAK/STAT, and phosphatidylinositol 3-kinase/AKT, have emerged to treat this intractable disease. A combination of immunological recovery from immune deficiency, overcoming the immune escape, and the development of more effective drugs will be vital for improving the outcomes of PEL patients in the future.

Genetic heterogeneity and mutational signature in Chinese Epstein-Barr virus-positive diffuse large B-cell lymphoma

Epstein-Barr virus (EBV)-positive diffuse large B-cell lymphoma (EBV+ DLBCL) is typically an aggressive tumor in elderly patients. However, in a subset of young patients, EBV+ DLBCL follows a relatively indolent clinical course and exhibits a good response to chemotherapy. This lymphoma comprises polymorphous lymphoma and large cell lymphomas subtypes, with the latter subtype showing a significantly poorer prognosis. It is unknown whether the genetic background differs between age groups and histopathological subtypes. To investigate the genetic basis, heterogeneity, and recurrently mutated genes in EBV+ DLBCL, we performed whole-exome sequencing of DNA from 11 tissue samples of this lymphoma. Sequencing revealed that the most common substitution was the transition C>T/G>A. Genetic features—including the numbers of mutated genes in exonic region, single-nucleotide variants (SNV), and indels—did not significantly differ between age groups or histological subtypes. Matching with the COSMIC database revealed that the main mutational signature was signature 3, which is associated with failure of DNA double-strand break-repair by homologous recombination. Mutant-Allele Tumor Heterogeneity (MATH) scores showed that EBV+ DLBCL exhibited broad intratumor heterogeneity, and were positively correlated with Ann Arbor Stage and ≥2 extranodal lesion sites. We identified 57 selected recurrently mutated genes. The most commonly mutated five genes—LNP1 (11/11), PRSS3 (10/11), MUC3A (9/11), FADS6 (9/11), and TRAK1 (8/11)—were validated by Sanger sequencing. These mutated genes have not previously been identified. Overall, our present results demonstrate the tremendous genetic heterogeneity underlying EBV+ DLBCLs, and highlight the need for personalized therapeutic approaches to treating these patients.

The effects of mutational processes and selection on driver mutations across cancer types

Epidemiological evidence has long associated environmental mutagens with increased cancer risk. However, links between specific mutation-causing processes and the acquisition of individual driver mutations have remained obscure. Here we have used public cancer sequencing data from 11,336 cancers of various types to infer the independent effects of mutation and selection on the set of driver mutations in a cancer type. First, we detect associations between a range of mutational processes, including those linked to smoking, ageing, APOBEC and DNA mismatch repair (MMR) and the presence of key driver mutations across cancer types. Second, we quantify differential selection between well-known alternative driver mutations, including differences in selection between distinct mutant residues in the same gene. These results show that while mutational processes have a large role in determining which driver mutations are present in a cancer, the role of selection frequently dominates.

The ubiquitous 'cancer mutational signature' 5 occurs specifically in cancers with deleted FHIT alleles

The FHIT gene is located at the fragile FRA3B locus where activation by carcinogen-induced and endogenous replication stress causes FHIT deletions even in normal cells over a lifetime. Our lab has shown that loss of FHIT expression causes genome instability and provides single-strand DNA substrates for APOBEC3B hypermutation, in line with evidence that FHIT locus deletions occur in many cancers. Based on these biological features, we hypothesized that FHIT loss drives development of COSMIC mutational signature 5 and here provide evidence, including data mining of >6,500 TCGA samples, that FHIT is the cancer-associated gene with copy number alterations correlating most significantly with signature 5 mutation rate. In addition, tissues of Fhit-deficient mice exhibit a mutational signature strongly resembling signature 5 (cosine similarity value = 0.89). We conclude that FHIT loss is a molecular determinant for signature 5 mutations, which occur in all cancer types early in cancer development, are clock-like, and accelerated by carcinogen exposure. Loss of FHIT caretaker function may be a predictive and preventive marker for cancer development.

Mutational signatures associated with tobacco smoking in human cancer

Tobacco smoking increases the risk of at least 17 classes of human cancer. We analyzed somatic mutations and DNA methylation in 5243 cancers of types for which tobacco smoking confers an elevated risk. Smoking is associated with increased mutation burdens of multiple distinct mutational signatures, which contribute to different extents in different cancers. One of these signatures, mainly found in cancers derived from tissues directly exposed to tobacco smoke, is attributable to misreplication of DNA damage caused by tobacco carcinogens. Others likely reflect indirect activation of DNA editing by APOBEC cytidine deaminases and of an endogenous clocklike mutational process. Smoking is associated with limited differences in methylation. The results are consistent with the proposition that smoking increases cancer risk by increasing the somatic mutation load, although direct evidence for this mechanism is lacking in some smoking-related cancer types.

Clock-like mutational processes in human somatic cells

During the course of a lifetime, somatic cells acquire mutations. Different mutational processes may contribute to the mutations accumulated in a cell, with each imprinting a mutational signature on the cell's genome. Some processes generate mutations throughout life at a constant rate in all individuals, and the number of mutations in a cell attributable to these processes will be proportional to the chronological age of the person. Using mutations from 10,250 cancer genomes across 36 cancer types, we investigated clock-like mutational processes that have been operating in normal human cells. Two mutational signatures show clock-like properties. Both exhibit different mutation rates in different tissues. However, their mutation rates are not correlated, indicating that the underlying processes are subject to different biological influences. For one signature, the rate of cell division may influence its mutation rate. This study provides the first survey of clock-like mutational processes operating in human somatic cells.