APOBEC Mutation Research Articles

Abstract P5-06-03: A solution to the APOBEC mutation paradox in breast cancer

Abstract Breast cancer has a major genetic component, and individual tumors can have thousands to tens-of-thousands of mutations. Large proportions are cytosine mutations in TCA and TCT trinucleotide motifs. A variety of studies have attributed this mutation signature to the APOBEC family of DNA cytosine deaminases, with previous literature favoring APOBEC3B as the most likely enzyme to bee responsible (from a total of nine active family members). However, breast tumor genomic DNA sequences from patients with a naturally occurring germline deletion of the entire APOBEC3B gene still show an intact APOBEC mutation signature. To resolve this paradox, we tested the hypothesis that the only other functionally dimorphic APOBEC family member, APOBEC3H, may be responsible. First, we performed an unbiased genetic analysis of TCGA data sets and showed that APOBEC3B-null tumors with this mutational bias have at least one copy of the haplotype-I variant of APOBEC3H, despite weak genetic linkage between these two genes (n=14/14). Remarkably, breast tumors without APOBEC3B and APOBEC3H haplotype-I showed no evidence for an APOBEC mutation signature (n=3/3). The proportion of APOBEC signature mutations between the cohort with and the cohort without APOBEC3H haplotype-I was highly significant (p<0.00001 and lower), indicating that these two enzymes alone may be fully responsible (and casting serious doubt on the potential role of APOBEC3A). Second, although deemed inactive in prior studies, we discovered that APOBEC3H haplotype-I has robust activity in biochemical assays with mutation of TCA-containing single-stranded DNA and in cellular assays with hypermutation of HIV-1 cDNA and mutation of a genomic DNA eGFP reporter gene. Third, the subcellular localization APOBEC3H haplotype-I and other expressed variants of this enzyme differed with only the latter showing significant nuclear localization in a panel of cell lines. The genomic mutagenicity and localization of APOBEC3H haplotype-I is attributable to a glycine at position 105, whereas other variants of this enzyme that have an arginine at this position, suggesting a general molecular mechanism. These studies combine to offer a parsimonious solution to the APOBEC mutation paradox with APOBEC3B and APOBEC3H haplotype-I, together, accounting for the strong “APOBEC mutation signature” observed in many breast cancers. Citation Format: Harris RS, Starrett GJ, McCann JL, Carpenter MA. A solution to the APOBEC mutation paradox in breast cancer [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P5-06-03.

Dysregulated Aid/Apobec Family Proteins Promote Genomic Instability in Multiple Myeloma

The AID/APOBEC family of cytidine deaminase proteins includes AID (activity induced deaminase), and 10 related APOBEC enzymes (A1, A2, A3A, A3B, A3C, A3D, A3F, A3G, A3H and A4). AID has been well-studied for its role in somatic hyper mutation and class switch recombination of immunoglobulin genes whereas APOBECs (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) have been shown to have roles in mRNA editing and in antiviral immunity. Dysregulated activity of APOBECs causes C >T transitions or C>G, C>A transversions in DNA. We have recently shown APOBEC signature mutation pattern in multiple myeloma (MM) genomes (Bolli et al Nat. Comm. 2014), and interestingly, the APOBEC mutation signature correlates with sub clonal diversity in myeloma. A role for the AID/APOBECs in generation of somatic mutations has also been proposed in a variety of other cancers based on identification of APOBEC signature mutationsIn order to understand which APOBECs are dysregulated in myeloma, we performed RNA sequencing analysis of primary myeloma cells from 409 newly-diagnosed MM patients and myeloma cell lines. Our analysis showed elevated expression of several APOBEC family members; mainly A3A, A3B, A3C, and A3G. We then optimized a plasmid-based functional assay and found high cytidine deaminase activity in extracts from a number of myeloma cell lines and patient derived CD138+ cells compared to CD138+ cells from healthy donors, suggesting that APOBECs are dysregulated in myeloma. We then investigated the impact of elevated APOBEC expression/function on overall genome maintenance and acquisition of genomic changes (such as amplifications, deletions) overtime. We used shRNA-mediated knockdown of specific APOBEC proteins in myeloma cell lines and investigated the acquisition of genomic changes in control and knockdown cells during their growth in culture, using SNP (Single Nucleotide Polymorphism) arrays and WGS (whole genome sequencing) platforms. Our results with both approaches showed significant reduction in the accumulation of copy number changes (both amplifications and deletions) and overall mutation load after APOBEC knockdown. Evaluation with both the SNP and WGS showed that when control and APOBEC knockdown cells were cultured for three weeks, the acquisition of new copy number and mutational changes throughout genome were reduced by ~50%.We next investigated the relationship between APOBEC expression/activity in MM and other DNA repair pathways. Using an in vitro HR activity assay, we measured HR activity in extracts from control and APOBEC knockdown cells. Depletion of APOBEC proteins resulted in 50-80% reduction in in vitro HR activity of the extracts. We also evaluated correlation between HR activity and gene expression using RNA-seq data from myeloma cells derived from 100 patients at diagnosis and identified the genes whose expression correlated with HR activity. Elevated expression of APOBECs 3D, 3G and 3F significantly correlated with high HR activity (R=0.3; P≤0.02), suggesting their relevance to HR. Analyzing genomic copy number information for each patient we have also observed significant correlation between higher expression of A3G and increased genomic instability in this dataset (P=0.0045).In summary, our study shows that dysregulated APOBECs induce mutations and genomic instability, and inhibiting APOBEC activity could reduce the rate of accumulation of ongoing genomic changes. This data sheds light on biology of the disease as well as clonal evolution. DisclosuresMunshi:Amgen: Consultancy; Oncopep: Patents & Royalties; Celgene: Consultancy; Janssen: Consultancy; Takeda: Consultancy; Merck: Consultancy; Pfizer: Consultancy.

Elucidating Genomic Characteristics of Lung Cancer Progression from In Situ to Invasive Adenocarcinoma.

To examine the diversity of somatic alterations and clonal evolution according to aggressiveness of disease, nineteen tumor-blood pairs of ‘formerly bronchiolo-alveolar carcinoma (BAC)’ which had been reclassified into preinvasive lesion (adenocarcinoma in situ; AIS), focal invasive lesion (minimally invasive adenocarcinoma; MIA), and invasive lesion (lepidic predominant adenocarcinoma; LPA and non-lepidic predominant adenocarcinoma; non-LPA) according to IASLC/ATS/ERS 2011 classification were explored by whole exome sequencing. Several distinct somatic alterations were observed compare to the lung adenocarcinoma study from the Cancer Genome Atlas (TCGA). There were higher numbers of tumors with significant APOBEC mutation fold enrichment (73% vs. 58% TCGA). The frequency of KRAS mutations was lower in our study (5% vs. 32% TCGA), while a higher number of mutations of RNA-splicing genes, RBM10 and U2AF1, were found (37% vs. 11% TCGA). We found neither mutational pattern nor somatic copy number alterations that were specific to AIS/MIA. We demonstrated that clonal cell fraction was the only distinctive feature that discriminated LPA/non-LPA from AIS/MIA. The broad range of clonal frequency signified a more branched clonal evolution at the time of diagnosis. Assessment of tumor clonal cell fraction might provide critical information for individualized therapy as a prognostic factor, however this needs further study.

CSN1 Somatic Mutations in Penile Squamous Cell Carcinoma.

Other than an association with HPV infection, little is known about the genetic alterations determining the development of penile cancer. Although penile cancer is rare in the developed world, it presents a significant burden in developing countries. Here, we report the findings of whole-exome sequencing (WES) to determine the somatic mutational landscape of penile cancer. WES was performed on penile cancer and matched germline DNA from 27 patients undergoing surgical resection. Targeted resequencing of candidate genes was performed in an independent 70 patient cohort. Mutation data were also integrated with DNA methylation and copy-number information from the same patients. We identified an HPV-associated APOBEC mutation signature and an NpCpG signature in HPV-negative disease. We also identified recurrent mutations in the novel penile cancer tumor suppressor genes CSN1(GPS1) and FAT1 Expression of CSN1 mutants in cells resulted in colocalization with AGO2 in cytoplasmic P-bodies, ultimately leading to the loss of miRNA-mediated gene silencing, which may contribute to disease etiology. Our findings represent the first comprehensive analysis of somatic alterations in penile cancer, highlighting the complex landscape of alterations in this malignancy. Cancer Res; 76(16); 4720-7. ©2016 AACR.

Abstract LB-323: The genomic landscape of non-muscle-invasive bladder cancer: implications for molecular classification and treatment

Abstract Non-muscle-invasive bladder cancers (NMIBC) and muscle-invasive bladder cancers (MIBC) show distinct molecular and clinical features and are considered to follow different pathogenesis pathways. MIBC commonly metastasise (>50%) and have poor prognosis whereas NMIBC, particularly low-stage low-grade tumors, rarely progress to invade muscle (<5%). Patients diagnosed with NMIBC frequently suffer disease recurrence, demanding long-term disease monitoring and repeated resection of recurrences. Large MIBC cohorts have been extensively interrogated on a genome-wide scale but relatively few NMIBC have been studied at this level. To address this we determined the mutational landscape of 82 non-muscle-invasive tumors (stage Ta grade 2) using whole-exome and targeted deep sequencing. Whole-exome sequencing identified an average of 124±79 synonymous and non-synonymous somatic mutations (single nucleotide substitutions and indels) per sample, giving mean and median somatic mutation rates of 2.41 and 1.64 per megabase, respectively. Overall, 47% of nucleotide substitutions were C>T transitions, followed by C>G transversions (24.8%). Similar to MIBC, APOBEC mutagenesis was the strongest source of mutation in NMIBC with 75% of samples showing up to 5-fold enrichment of the APOBEC mutation signature. Significantly mutated genes included genes implicated in epigenetic regulation and several genes that have not previously been reported as significantly mutated in bladder cancer. Comparison of mutation frequencies in NMIBC with those found in MIBC revealed both known differences (absence of TP53 mutations and high frequency of FGFR3, PIK3CA and STAG2 mutations in NMIBC) and novel findings including a higher frequency of mutations in chromatin-state regulators in NMIBC. Notably CDKN1A, RB1, ERCC2, ERBB3 and FBXW7, which are mutated in >10% of MIBC were not significantly mutated in TaG2 tumors, further delineating the distinct pathogenesis pathways of NMIBC and MIBC. Whole-genome expression array profiling and immunohistochemistry using a panel of markers were carried out to further examine the molecular features of NMIBC samples and define potentially clinically relevant subtypes. These data provide a detailed view of the genomic landscape of NMIBC that highlights chromatin modification as a key area for consideration in the development of potential future therapeutic approaches to the treatment of patients with non-muscle-invasive disease. Citation Format: Carolyn D. Hurst, Olivia Alder, Fiona M. Platt, Alastair Droop, Lucy F. Stead, Julie E. Burns, George J. Burghel, Sunjay Jain, Leszek J. Klimczak, Helen Lindsay, Jo-An Roulson, Claire F. Taylor, Helene Thygesen, Angus J. Cameron, Anne J. Ridley, Helene R. Mott, Dmitry A. Gordenin, Margaret A. Knowles. The genomic landscape of non-muscle-invasive bladder cancer: implications for molecular classification and treatment. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-323.

Correlation of Apobec Mrna Expression with overall Survival and pd-l1 Expression in Urothelial Carcinoma.

Metastatic urothelial carcinoma (mUC) has a very high mutational rate and is associated with an APOBEC mutation signature. We examined the correlation of APOBEC expression with overall survival (OS) and PD-L1 expression in a cohort of 73 mUC patients. mRNA expression of APOBEC3 family of genes (A3A, A3B, A3C, A3F_a, A3F_b, A3G, A3H) was measured using Nanostring. PD-L1 expression, evaluated by immunohistochemistry, on tumor infiltrating mononuclear cells (TIMCs) and tumor cells was scored from 0 to 4, with 2–4 being positive. Wilcoxon’s non-parametric tests assessed the association of APOBEC and PD-L1. The Cox regression model assessed the association of APOBEC with OS. All APOBEC genes were expressed in mUC. Increased A3A, A3D, and A3H expression associates with PD-L1 positive TIMCs (p = 0.0009, 0.009, 0.06). Decreased A3B expression was marginally associated with PD-L1 positive TIMCs expression (p = 0.05). Increased A3F_a and A3F_b expression was associated with increased expression of PD-L1 on tumor cells (p = 0.05). Increased expression of A3D and A3H was associated with longer OS (p = 0.0009). Specific APOBEC genes have different effects on mUC in terms of survival and PD-L1 expression. A3D and A3H may have the most important role in mUC as they are associated with OS and PD-L1 TIMC expression.

APOBEC-Induced Cancer Mutations Are Uniquely Enriched in Early-Replicating, Gene-Dense, and Active Chromatin Regions.

An antiviral component of the human innate immune system-the APOBEC cytidine deaminases-was recently identified as a prominent source of mutations in cancers. Here, we investigated the distribution of APOBEC-induced mutations across the genomes of 119 breast and 24 lung cancer samples. While the rate of most mutations is known to be elevated in late-replicating regions that are characterized by reduced chromatin accessibility and low gene density, we observed a marked enrichment of APOBEC mutations in early-replicating regions. This unusual mutagenesis profile may be associated with a higher propensity to form single-strand DNA substrates for APOBEC enzymes in early-replicating regions and should be accounted for in statistical analyses of cancer genome mutation catalogs aimed at understanding the mechanisms of carcinogenesis as well as highlighting genes that are significantly mutated in cancer.

Impact of APOBEC Mutations on CD8+ T Cell Recognition of HIV Epitopes Varies Depending on the Restricting HLA.

We previously showed that APOBEC-mediated mutations in HIV CD8 T-cell epitopes generally reduce recognition by CD8 T cells. Here, we examined this effect in the context of histocompatibility-linked leukocyte antigen (HLA) alleles differentially associated with disease progression rates. For HLA-B57-restricted epitopes, APOBEC mutations generally diminished CD8 T cell recognition. Conversely, recognition of HLA-B35-restricted epitopes was consistently enhanced. For epitopes that can be presented by either HLA-A2 or A3, the same APOBEC mutation had differential effects on CD8 T cell recognition, depending on the individual's HLA genotype. The pattern of HLA dependence provides additional evidence that APOBEC action is channeled toward cytotoxic CD8 T-cell escape.

Apobec Family Mutational Signatures Are Associated with Poor Prognosis Translocations in Multiple Myeloma

Aberrant chromosomal translocations are seen in ~40% of presenting patients and predominantly involve the IGH locus at 14q32. The five main translocations involving the IGH locus are t(4;14), t(6;14), t(11;14), t(14;16) and t(14;20), which result in over-expression of MMSET/FGFR3, CCND3, CCND1, MAF and MAFB, respectively. In previous clinical trials we have shown that the t(4;14), t(14;16) and t(14;20) are associated with a poor prognosis. In initial sequencing studies of myeloma it has been noted that the spectrum of mutations fall into two groups, one of which is characterised by an APOBEC signature. This signature comprises of C>T, C>G and C>A mutations in a TpC context and comprises only a subset of samples, with the rest having a rather generic mutation signature representing an intrinsic mutational process occurring as a result of the spontaneous deamination of methylated cytosine to thymine.Whole exome sequencing was performed on 463 presentation patients enrolled into the UK Myeloma XI trial. DNA was extracted from germline DNA and CD138+ plasma cells and whole exome sequencing was performed using SureSelect (Agilent). In addition to capturing the exome, extra baits were added covering the IGH, IGK, IGL and MYCloci in order to determine the breakpoints associated with translocations in these genes. Tumor and germline DNA were sequenced to a median of 60x and data processed to generate copy number, acquired variants and translocation breakpoints in the tumor. Progression-free and overall survival was measured from initial randomization and median follow up for this analysis was 25 months. These combined data allow us to examine the effect of translocations on the mutational spectra in myeloma and determine any associations with progression-free or overall survival.Translocations were detected in 232 (50.1%) patients of which 59 patients (12.7%) had a t(4;14), 86 patients (18.6%) a t(11;14), 17 patients (3.7%) a t(14;16), 5 patients (1%) a t(6;14) and 4 patients (0.9%) a t(14;20). MYC translocations were found in 85 patients (18.4%). Using the tiled regions we were able to detect a mutational signature, kataegis, where regional clustering of mutations can be indicative of somatic genomic rearrangements. We found the hallmarks of kataegis in 15 samples (3.2%), where there was enrichment for TpCpH mutations with an inter-mutational distance <1 kb. Where we detected kataegis surrounding MYC there was also an inter-chromosomal translocation involving either IGK or IGL. Interestingly, the partner chromosomes also showed signs of kataegis e.g. in the t(2;8) kataegis was found at IGK and MYC and in the t(8;22) kataegis was found at MYC and IGL. APOBECs are thought to be involved in the generation of kataegis and as such this co-localisation is indicative of APOBEC involvement in the generation of MYCbreakpoints.We found mutation of translocation partner oncogenes, in particular CCND1 was mutated in 10 patients with the t(11;14). There was an association of mutated CCND1 with a poor prognosis when compared with non-mutated t(11;14) patients (OS median of 20.2 months vs. not reached, p=0.005). Mutations were also seen in FGFR3, MAF and MAFB but only in the samples with the respective translocations. The mutations are likely due to somatic hypermutation mediated by AID, an APOBEC family member.We found that t(14;16) and t(14;20) samples have a significantly higher number of mutations compare to the other translocation groups (p=1.65x10-5). These mutations were enriched for those with an APOBEC signature (T(C>T)A, p=9.1x10-5; T(C>T)T, p=0.0014; T(C>G)A, p=0.001; T(C>G)T, p=0.0064), indicating that the ‘maf’ translocation groups are characterized by APOBEC signature mutations, specifically APOBEC3B. When samples are assigned to either an APOBEC or non-APOBEC group the ‘maf’ translocations account for 66.6% of samples in the APOBEC group but only 1.3% of the non-APOBEC group.Here we show three different mutational signatures mediated by the APOBEC family: translocation partner mutation by AID, APOBEC signature mediated by APOBEC3B, and kataegis mediated by an unknown APOBEC family member. We also show for the first time a clinical impact of APOBEC mutations and their association with a poor prognosis. The poor prognosis of this mutational signature is inextricably linked to a high mutation load and the adverse t(14;16) and t(14;20) translocation subgroups. DisclosuresWalker:Onyx Pharmaceuticals: Consultancy, Honoraria.

APOBEC3B gene overexpression in non-small-cell lung cancer.

Recent study results have demonstrated that a subclass of apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like (APOBEC) cytidine deaminase may induce mutation clusters in various types of cancer. From the Cancer Genome Altas, an APOBEC mutation pattern was identified in bladder, cervical, breast, head and neck and lung cancers. In the present study, APOBEC3B mRNA expression was investigated using quantitative reverse transcription-polymerase chain reaction (RT-qPCR) assay using LightCycler in surgically treated non-small-cell lung cancer (NSCLC) cases. Additionally, 88 surgically removed Japanese NSCLC cases were analyzed for mRNA level. The results showed that APOBEC3B/β-actin mRNA levels were significantly higher in lung cancer (1,598.481±6,465.781) when compared to adjacent normal lung tissues (2,116.639±8,337.331, P=0.5453). The tumor/normal (T/N) ratio of APOBEC3B/β-actin mRNA levels was not different within the gender, age, smoking status and pathological stages. The T/N ratio of APOBEC3B/β-actin mRNA levels was not significantly different in epidermal growth factor receptor (EGFR) or Kras mutation-positive cases as compared to the wild-type cases.

An APOBEC cytidine deaminase mutagenesis pattern is widespread in human cancers

Recent studies indicate that a subclass of APOBEC cytidine deaminases, which convert cytosine to uracil during RNA editing and retrovirus or retrotransposon restriction, may induce mutation clusters in human tumors. We show here that throughout cancer genomes APOBEC-mediated mutagenesis is pervasive and correlates with APOBEC mRNA levels. Mutation clusters in whole-genome and exome data sets conformed to the stringent criteria indicative of an APOBEC mutation pattern. Applying these criteria to 954,247 mutations in 2,680 exomes from 14 cancer types, mostly from The Cancer Genome Atlas (TCGA), showed a significant presence of the APOBEC mutation pattern in bladder, cervical, breast, head and neck, and lung cancers, reaching 68% of all mutations in some samples. Within breast cancer, the HER2-enriched subtype was clearly enriched for tumors with the APOBEC mutation pattern, suggesting that this type of mutagenesis is functionally linked with cancer development. The APOBEC mutation pattern also extended to cancer-associated genes, implying that ubiquitous APOBEC-mediated mutagenesis is carcinogenic.