B-Myb Induces APOBEC3B Expression Leading to Somatic Mutation in Multiple Cancers

Her2 amplification, Rel-A, and Bach1 can influence APOBEC3A expression in breast cancer cells

APOBEC-induced mutations occur in 50% of sequenced human tumors, with APOBEC3A (A3A) being a major contributor to mutagenesis in breast cancer cells. The mechanisms that cause A3A activation and mutagenesis in breast cancers are still unknown. Here, we describe factors that influence basal A3A mRNA transcript levels in breast cancer cells. We found that basal A3A mRNA correlates with A3A protein levels and predicts the amount of APOBEC signature mutations in a panel of breast cancer cell lines, indicating that increased basal transcription may be one mechanism leading to breast cancer mutagenesis. We also show that alteration of ERBB2 expression can drive A3A mRNA levels, suggesting the enrichment of the APOBEC mutation signature in Her2-enriched breast cancer could in part result from elevated A3A transcription. Hierarchical clustering of transcripts in primary breast cancers determined that A3A mRNA was co-expressed with other genes functioning in viral restriction and interferon responses. However, reduction of STAT signaling via inhibitors or shRNA in breast cancer cell lines had only minor impact on A3A abundance. Analysis of single cell RNA-seq from primary tumors indicated that A3A mRNA was highest in infiltrating immune cells within the tumor, indicating that correlations of A3A with STAT signaling in primary tumors may be result from higher immune infiltrates and are not reflective of STAT signaling controlling A3A expression in breast cancer cells. Analysis of ATAC-seq data in multiple breast cancer cell lines identified two transcription factor sites in the APOBEC3A promoter region that could promote A3A transcription. We determined that Rel-A, and Bach1, which have binding sites in these peaks, elevated basal A3A expression. Our findings highlight a complex and variable set of transcriptional activators for A3A in breast cancer cells.

Genetic insights into carbohydrate sulfotransferase 8 and its impact on the immunotherapy efficacy of cancer

Immune checkpoint blockade (ICB) is a promising therapy for solid tumors, but its effectiveness depends on biomarkers that are not precise. Here, we utilized genome-wide association study to investigate the association between genetic variants and tumor mutation burden to interpret ICB response. We identified 16 variants (p < 5 × 10-8) probed to 17 genes on 9 chromosomes. Subsequent analysis of one of the most significant loci in 19q13.11 suggested that the rs111308825 locus at the enhancer is causal, as its A allele impairs KLF2 binding, leading to lower carbohydrate sulfotransferase 8 (CHST8) expression. Breast cancer cells expressing CHST8 suppress T cell activation, and Chst8 loss attenuates tumor growth in a syngeneic mouse model. Further investigation revealed that programmed death-ligand 1 (PD-L1) and its homologs could be sulfated by CHST8, resulting in M2-like macrophage enrichment in the tumor microenvironment. Finally, we confirmed that low-CHST8 tumors have better ICB response, supporting the genetic effect and clinical value of rs111308825 for ICB efficacy prediction.

Gene mutation analysis of oral submucous fibrosis cancerization in Hainan Island

Objective

The sequencing panel composed of 61 target genes was used to explore the related mutation genes of oral squamous cell carcinoma (OSCC) and oral submucous fibrosis (OSF) cancerization, so as to provide a theoretical basis for the early diagnosis of oral submucous fibrosis cancerization, find the most important mutations in OSF cancerization, and more targeted prevention of OSF cancerization.

Methods

A total of 74 clinically diagnosed samples were included, including 36 cases of OSCC and 38 cases of OSF cancer patients. DNA was extracted, and targeted gene panel sequencing technology was used to analyze the gene frequency of pathogenic mutation sites in clinical samples.

Results

Gene panel sequencing analysis showed that there were 69 mutations in 18 genes in OSCC and OSF cancerous specimens. The results of gene panel sequencing were screened, and 18 mutant genes were finally screened out and their mutation frequencies in the samples were analyzed. According to the frequency of gene mutations from high to low, they were TP53, FLT4, PIK3CA, CDKN2A, FGFR4, HRAS, BRCA1, PTPN11, NF1, KMT2A, RB1, PTEN, MSH2, MLH1, KMT2D, FLCN, BRCA2, APC. The mutation frequency of FLT4 gene was significantly higher than that of OSCC group (P < 0.05).

Conclusion

FLT4 gene may be related to OSF cancerization and is expected to be an early diagnostic biomarker for OSF cancerization.

ILF2 enhances the DNA cytosine deaminase activity of tumor mutator APOBEC3B in multiple myeloma cells

DNA cytosine deaminase APOBEC3B (A3B) is an endogenous source of mutations in many human cancers, including multiple myeloma. A3B proteins form catalytically inactive high molecular mass (HMM) complexes in nuclei, however, the regulatory mechanisms of A3B deaminase activity in HMM complexes are still unclear. Here, we performed mass spectrometry analysis of A3B-interacting proteins from nuclear extracts of myeloma cell lines and identified 30 putative interacting proteins. These proteins are involved in RNA metabolism, including RNA binding, mRNA splicing, translation, and regulation of gene expression. Except for SAFB, these proteins interact with A3B in an RNA-dependent manner. Most of these interacting proteins are detected in A3B HMM complexes by density gradient sedimentation assays. We focused on two interacting proteins, ILF2 and SAFB. We found that overexpressed ILF2 enhanced the deaminase activity of A3B by 30%, while SAFB did not. Additionally, siRNA-mediated knockdown of ILF2 suppressed A3B deaminase activity by 30% in HEK293T cell lysates. Based on these findings, we conclude that ILF2 can interact with A3B and enhance its deaminase activity in HMM complexes.

CENPA is one of the potential key genes associated with the proliferation and prognosis of ovarian cancer based on integrated bioinformatics analysis and regulated by MYBL2

Background

Ovarian cancer (OV) is a highly lethal disease, and the fifth leading cause of all cancer-related deaths in women. The study aimed to identify potential key genes associated with the proliferation and prognosis of OV.

Methods

Differentially expressed genes (DEGs) between ovarian cancer and normal tissues were screened by the robust rank aggregation (RRA) method. The expression of CENPA and MYBL2 were examined in SKOV3 and A2780 ovarian cancer cell lines and tumor tissues by qRT-PCR and western blot. Small RNA interference assays, plasmid overexpression assays and EdU assays were used to validate the proliferative effect of the MYBL2-CENPA axis in ovarian cancer cell lines. The ChIP assay was used to verify the direct regulation of MYBL2 on CENPA.

Results

133 up-regulated genes and 158 down-regulated genes were identified, and the up-regulated genes mainly enrichment in cell cycle. The three up-regulated gene with DNA separation (CENPA, CENPF and CEP55) might be tightly correlated with proliferation and prognosis of OV. Knockdown CENPA expression inhibited the proliferation of A2780 and SKOV3 cells After the knockout of MYBL2, the expression of CENPA significantly decreased. MYBL2 directly binds to the promoter region of CENPA.

Conclusions

The MYBL2-CENPA pathway plays an important role in the proliferation of ovarian cancer cells, suggesting that this pathway may be a potential target for the treatment of ovarian cancer.

A common variant in 11q23.3 associated with hyperlipidemia is mediated by the binding and regulation of GATA4

Large-scale genome-wide associations comprising multiple studies have identified hundreds of genetic loci commonly associated with hyperlipidemia-related phenotypes. However, single large cohort remains necessary in aiming to investigate ethnicity-specific genetic risks and mechanical insights. A community-based cohort comprising 23,988 samples that included both genotype and biochemical information was assembled for the genome-wide association analysis (GWAS) of hyperlipidemia. The analysis identified fifty genetic variants (P < 5 × 10⁻⁸) on five different chromosomes, and a subsequent validation analysis confirmed the significance of the lead variants. Integrated analysis combined with cell-based experiments of the most statistically significant locus in 11q23.3 revealed rs651821 (P = 4.52 × 10⁻⁷⁶) as the functional variant. We showed transcription factor GATA4 preferentially binds the T allele of rs651821, the protective allele for hyperlipidemia, which promoted APOA5 expression in liver cells and individuals with the TT genotype of rs651821. As GATA4-APOA5 axis maintains triglyceride homeostasis, GATA4 activation by phenylephrine implies synergism for lowering triglyceride levels in hyperlipidemia patients. Our study demonstrates that rs651821 mediates APOA5 activation via allele-specific regulation by GATA4. We suggest elevating GATA4 activity could provide a therapeutic potential for treating the development of hyperlipidemia.

m6A demethylation of cytidine deaminase APOBEC3B mRNA orchestrates arsenic-induced mutagenesis

The cytidine deaminase APOBEC3B (A3B) is an endogenous inducer of somatic mutations and causes chromosomal instability by converting cytosine to uracil in single-stranded DNA. Therefore, identification of factors and mechanisms that mediate A3B expression will be helpful for developing therapeutic approaches to decrease DNA mutagenesis. Arsenic (As) is one well-known mutagen and carcinogen, but the mechanisms by which it induces mutations have not been fully elucidated. Herein, we show that A3B is upregulated and required for As-induced DNA damage and mutagenesis. We found that As treatment causes a decrease of N6-methyladenosine (m6A) modification near the stop codon of A3B, consequently increasing the stability of A3B mRNA. We further reveal that the demethylase FTO is responsible for As-reduced m6A modification of A3B, leading to increased A3B expression and DNA mutation rates in a manner dependent on the m6A reader YTHDF2. Our in vivo data also confirm that A3B is a downstream target of FTO in As-exposed lung tissues. In addition, FTO protein is highly expressed and positively correlates with the protein levels of A3B in tumor samples from human non-small cell lung cancer patients. These findings indicate a previously unrecognized role of A3B in As-triggered somatic mutation and might open new avenues to reduce DNA mutagenesis by targeting the FTO/m6A axis.

MEK1/2 regulates APOBEC3B and polymerase iota-induced mutagenesis in head and neck cancer cells

Resistance to chemotherapy provides a major challenge in treatment of metastatic cancer. Prolonged exposure to almost any drug regimen leads to the formation of resistant subclones in almost all advanced solid tumors. Tumor heterogeneity because of intrinsic genetic instability is seen as one of the major contributing factors. In this work, we present evidence that genetic instability measured by mutation frequency is induced by treatment with the EGFR inhibitor afatinib or cisplatin in head and neck squamous cancer cells. We find that APOBEC3B and polymerase iota are upregulated, and inhibition of MEK1/2 by U0126 leads to downregulation on the protein level. Costimulation of afatnib and cisplatin with U0126 leads to a significantly lower mutation frequency. These findings may represent a molecular mechanism for dynamically controlling genetic instability during chemotherapy in head and neck squamous cell carcinoma (HNSCC) cancer cells.

The von Hippel-Lindau Cullin-RING E3 ubiquitin ligase regulates APOBEC3 cytidine deaminases

The 7 members of the A3 family of cytidine deaminases (A3A to A3H) share a conserved catalytic activity that converts cytidines in single-stranded (ss) DNA into uridines, thereby inducing mutations. After their initial identification as cell-intrinsic defenses against HIV and other retroviruses, A3s were also found to impair many additional viruses. Moreover, some of the A3 proteins (A3A, A3B, and A3H haplotype I) are dysregulated in cancer cells, thereby causing chromosomal mutations that can be selected to fuel progression of malignancy. Viral mechanisms that increase transcription of A3 genes or induce proteasomal degradation of A3 proteins have been characterized. However, only a few underlying biological mechanisms regulating levels of A3s in uninfected cells have been described. Here, we characterize that the von Hippel-Lindau tumor suppressor (pVHL), via its CRL^pVHL, induces degradation of all 7 A3 proteins. Two independent lines of evidence supported the conclusion that the multiprotein CRL^pVHL complex is necessary for A3 degradation. CRL^pVHL more effectively induced degradation of nuclear, procancer A3 (A3B) than the cytoplasmic, antiretroviral A3 (A3G). These results identify specific cellular factors that regulate A3s post-translationally.

Comprehensive Mapping of Key Regulatory Networks that Drive Oncogene Expression

Gene expression is controlled by the collective binding of transcription factors to cis-regulatory regions. Deciphering gene-centered regulatory networks is vital to understanding and controlling gene misexpression in human disease; however, systematic approaches to uncovering regulatory networks have been lacking. Here we present high-throughput interrogation of gene-centered activation networks (HIGAN), a pipeline that employs a suite of multifaceted genomic approaches to connect upstream signaling inputs, trans-acting TFs, and cis-regulatory elements. We apply HIGAN to understand the aberrant activation of the cytidine deaminase APOBEC3B, an intrinsic source of cancer hypermutation. We reveal that nuclear factor κB (NF-κB) and AP-1 pathways are the most salient trans-acting inputs, with minor roles for other inflammatory pathways. We identify a cis-regulatory architecture dominated by a major intronic enhancer that requires coordinated NF-κB and AP-1 activity with secondary inputs from distal regulatory regions. Our data demonstrate how integration of cis and trans genomic screening platforms provides a paradigm for building gene-centered regulatory networks.

PIK3CA Mutation in the ShortHER Randomized Adjuvant Trial for Patients with Early HER2+ Breast Cancer: Association with Prognosis and Integration with PAM50 Subtype

PURPOSE

We explored the prognostic effect of PIK3CA mutation in HER2⁺ patients enrolled in the ShortHER trial.

PATIENTS AND METHODS

The ShortHER trial randomized 1,253 patients with HER2⁺ breast cancer to 9 weeks or 1 year of adjuvant trastuzumab combined with chemotherapy. PIK3CA hotspot mutations in exon 9 and 20 were analyzed by pyrosequencing. Expression of 60 genes, including PAM50 genes was measured using the nCounter platform.

RESULTS

A mutation of the PIK3CA gene was detected in 21.7% of the 803 genotyped tumors. At a median follow-up of 7.7 years, 5-year disease-free survival (DFS) rates were 90.6% for PIK3CA mutated and 86.2% for PIK3CA wild-type tumors [HR, 0.84; 95% confidence interval (CI), 0.56-1.27; P = 0.417]. PIK3CA mutation showed a favorable prognostic impact in the PAM50 HER2-enriched subtype (n = 232): 5-year DFS 91.8% versus 76.1% (log-rank P = 0.049; HR, 0.46; 95% CI, 0.21-1.02). HER2-enriched/PIK3CA mutated versus wild-type tumors showed numerically higher tumor-infiltrating lymphocytes (TIL) and significant upregulation of immune-related genes (including CD8A, CD274, PDCD1, and MYBL2, a proliferation gene involved in immune processes). High TILs as well as the upregulation of PDCD1 and MYBL2 were associated with a significant DFS improvement within the HER2-enriched subtype (HR, 0.82; 95% CI, 0.68-0.99; P = 0.039 for 10% TILs increment; HR, 0.81; 95% CI, 0.65-0.99; P = 0.049 for PDCD1 expression; HR, 0.72; 95% CI, 0.53-0.99; P = 0.042 for MYBL2 expression).

CONCLUSIONS

PIK3CA mutation showed no prognostic impact in the ShortHER trial. Within the HER2-enriched molecular subtype, patients with PIK3CA mutated tumors showed better DFS versus PIK3CA wild-type, which may be partly explained by upregulation of immune-related genes.

APOBEC3B High Expression in Gastroenteropancreatic Neuroendocrine Neoplasms and Association With Lymph Metastasis

PURPOSE

Apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3B (APOBEC3B) is known as a source of mutations in multiple cancers. Gastroenteropancreatic neuroendocrine neoplasms (GEP-NENs) are a group of heterogeneous tumors. However, the expression and significance of APOBEC3B in GEP-NENs remains unclear.

MATERIALS AND METHODS

A total of 158 cases of GEP-NENs, including 78 cases of biopsy or endoscopic submucosal dissection resection specimens and 83 cases of surgical resection specimens were collected in this study. The cases were grouped according to tumor classification grade, including 42 cases of neuroendocrine tumors G1 (NET G1), 36 cases of NET G2, 36 cases of NET G3, 44 cases of neuroendocrine carcinoma (NEC). All of the 158 tumors were immunohistochemically studied using a polyclonal antibody against APOBEC3B. We evaluated APOBEC3B expression in GEP-NENs and investigated the relationships among the immunoreactivity of APOBEC3B, clinical and pathologic features, such as age, sex, tumor site, Ki67 cell proliferation index, and lymph metastasis.

RESULTS

A total of 33 cases (78.6%) of NET G1 showed high expression of APOBEC3B. A total of 28 cases (77.8%) of NET G2 demonstrated high expression of APOBEC3B. In NET G3 and NEC cases, the positive rates were 52.8% and 2.3%, respectively. The expression of APOBEC3B in NETs was significantly higher than that in NECs, NET G1 and NET G2 were higher than NET G3, and the difference was statistically significant. APOBEC3B high expression cases have lower lymph node metastasis rate, lower Ki67 cell proliferation index.

CONCLUSIONS

In this study, APOBEC3B is highly expressed in GEP-NETs and is a predictor of lymph node metastasis in NET G3 and NEC cases. These findings might provide new insights into the biological mechanisms of GEP-NENs tumorigenesis and progression.

A functional variant near XCL1 gene improves breast cancer survival via promoting cancer immunity

Most genome-wide association studies (GWASs) identify genetic variants for breast cancer occurrence. In contrast, few are for recurrence and mortality. We conducted a GWAS on breast cancer survival after diagnosis in estrogen receptor-positive patients, including 953 Taiwanese patients with 159 events. Through Cox proportional hazard models estimation, we identified 24 risk SNPs with p < 1 × 10^-5 . Based on imputation and integrated analysis, one SNP, rs1024176 (located in 1q24.2, p = 2.43 × 10^-5 ) was found to be a functional variant associated with breast cancer survival and XCL1 gene expression. A series of experimental approaches, including cell-based analyses and CRISPR/Cas9 genome-editing system, were then used and identified the transcription factor MYBL2 was able to discriminately bind to the A allele of rs1024176, the protective variant for breast cancer survival, which promoted XCL1 expression, but not to the G allele of rs1024176. The chemokine XCL1 attracts type 1 dendritic cells (DC1s) to the tumor microenvironment. In breast cancer tissues, we applied a two-step Mendelian randomization analysis, using expression quantitative trait loci as instrumental variables, to confirm higher XCL1 expression was correlated with higher DC1 signatures and favorable disease progression, through the causal effect of rs1024176-A allele. Our study supports the genetic effect on preventing breast cancer survival through XCL1-induced DC1 recruitment in tumor microenvironment.

APOBEC3B reporter myeloma cell lines identify DNA damage response pathways leading to APOBEC3B expression

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) DNA cytosine deaminase 3B (A3B) is a DNA editing enzyme which induces genomic DNA mutations in multiple myeloma and in various other cancers. APOBEC family proteins are highly homologous so it is especially difficult to investigate the biology of specifically A3B in cancer cells. To easily and comprehensively investigate A3B function in myeloma cells, we used CRISPR/Cas9 to generate A3B reporter cells that contain 3×FLAG tag and IRES-EGFP sequences integrated at the end of the A3B gene. These reporter cells stably express 3xFLAG tagged A3B and the reporter EGFP and this expression is enhanced by known stimuli, such as PMA. Conversely, shRNA knockdown of A3B decreased EGFP fluorescence and 3xFLAG tagged A3B protein levels. We screened a series of anticancer treatments using these cell lines and identified that most conventional therapies, such as antimetabolites or radiation, exacerbated endogenous A3B expression, but recent molecular targeted therapeutics, including bortezomib, lenalidomide and elotuzumab, did not. Furthermore, chemical inhibition of ATM, ATR and DNA-PK suppressed EGFP expression upon treatment with antimetabolites. These results suggest that DNA damage triggers A3B expression through ATM, ATR and DNA-PK signaling.

Protein kinase A inhibits tumor mutator APOBEC3B through phosphorylation

APOBEC3B cytidine deaminase (A3B) catalyzes cytosine into uracil in single-strand DNA and induces C-to-T mutations in genomic DNA of various types of tumors. Accumulation of APOBEC signature mutations is correlated with a worse prognosis for patients with breast cancer or multiple myeloma, suggesting that A3B activity might be a cause of the unfavorable DNA mutations and clonal evolution in these tumors. Phosphorylation of conserved threonine residues of other cytidine deaminases, activation induced deaminase (AID) and APOBEC3G, inhibits their activity. Here we show that protein kinase A (PKA) physically binds to A3B and phosphorylates Thr214. In vitro deaminase assays and foreign DNA editing assays in cells confirm that phosphomimetic A3B mutants, T214D and T214E, completely lose deaminase activity. Molecular dynamics simulation of A3B phosphorylation reveals that Thr214 phosphorylation disrupts binding between the phospho-A3B catalytic core and ssDNA. These mutants still inhibit retroviral infectivity at least partially, and also retain full anti-retrotransposition activity. These results imply that PKA-mediated phosphorylation inhibits A3B mutagenic activity without destructing its innate immune functions. Therefore, PKA activation could reduce further accumulation of mutations in A3B overexpressing tumors.

Suboptimal T-cell Therapy Drives a Tumor Cell Mutator Phenotype That Promotes Escape from First-Line Treatment

Antitumor T-cell responses raised by first-line therapies such as chemotherapy, radiation, tumor cell vaccines, and viroimmunotherapy tend to be weak, both quantitatively (low frequency) and qualitatively (low affinity). We show here that T cells that recognize tumor-associated antigens can directly kill tumor cells if used at high effector-to-target ratios. However, when these tumor-reactive T cells were present at suboptimal ratios, direct T-cell-mediated tumor cell killing was reduced and the ability of tumor cells to evolve away from a coapplied therapy (oncolytic or suicide gene therapy) was promoted. This T-cell-mediated increase in therapeutic resistance was associated with C to T transition mutations that are characteristic of APOBEC3 cytosine deaminase activity and was induced through a TNFα and protein kinase C-dependent pathway. Short hairpin RNA inhibition of endogenous APOBEC3 reduced rates of tumor escape from oncolytic virus or suicide gene therapy to those seen in the absence of antitumor T-cell coculture. Conversely, overexpression of human APOBEC3B in tumor cells enhanced escape from suicide gene therapy and oncolytic virus therapy both in vitro and in vivo Our data suggest that weak affinity or low frequency T-cell responses against tumor antigens may contribute to the ability of tumor cells to evolve away from first-line therapies. We conclude that immunotherapies need to be optimized as early as possible so that, if they do not kill the tumor completely, they do not promote treatment resistance.

Cancer mutational burden is shaped by G4 DNA, replication stress and mitochondrial dysfunction

A hallmark of cancer is genomic instability, which can enable cancer cells to evade therapeutic strategies. Here we employed a computational approach to uncover mechanisms underlying cancer mutational burden by focusing upon relationships between 1) translocation breakpoints and the thousands of G4 DNA-forming sequences within retrotransposons impacting transcription and exemplifying probable non-B DNA structures and 2) transcriptome profiling and cancer mutations. We determined the location and number of G4 DNA-forming sequences in the Genome Reference Consortium Human Build 38 and found a total of 358,605 covering ∼13.4 million bases. By analyzing >97,000 unique translocation breakpoints from the Catalogue Of Somatic Mutations In Cancer (COSMIC), we found that breakpoints are overrepresented at G4 DNA-forming sequences within hominid-specific SVA retrotransposons, and generally occur in tumors with mutations in tumor suppressor genes, such as TP53. Furthermore, correlation analyses between mRNA levels and exome mutational loads from The Cancer Genome Atlas (TCGA) encompassing >450,000 gene-mutation regressions revealed strong positive and negative associations, which depended upon tissue of origin. The strongest positive correlations originated from genes not listed as cancer genes in COSMIC; yet, these show strong predictive power for survival in most tumor types by Kaplan-Meier estimation. Thus, correlation analyses of DNA structure and gene expression with mutation loads complement and extend more traditional approaches to elucidate processes shaping genomic instability in cancer. The combined results point to G4 DNA, activation of cell cycle/DNA repair pathways, and mitochondrial dysfunction as three major factors driving the accumulation of somatic mutations in cancer cells.

Spatial statistical tools for genome-wide mutation cluster detection under a microarray probe sampling system

Mutation cluster analysis is critical for understanding certain mutational mechanisms relevant to genetic disease, diversity, and evolution. Yet, whole genome sequencing for detection of mutation clusters is prohibitive with high cost for most organisms and population surveys. Single nucleotide polymorphism (SNP) genotyping arrays, like the Mouse Diversity Genotyping Array, offer an alternative low-cost, screening for mutations at hundreds of thousands of loci across the genome using experimental designs that permit capture of de novo mutations in any tissue. Formal statistical tools for genome-wide detection of mutation clusters under a microarray probe sampling system are yet to be established. A challenge in the development of statistical methods is that microarray detection of mutation clusters is constrained to select SNP loci captured by probes on the array. This paper develops a Monte Carlo framework for cluster testing and assesses test statistics for capturing potential deviations from spatial randomness which are motivated by, and incorporate, the array design. While null distributions of the test statistics are established under spatial randomness via the homogeneous Poisson process, power performance of the test statistics is evaluated under postulated types of Neyman-Scott clustering processes through Monte Carlo simulation. A new statistic is developed and recommended as a screening tool for mutation cluster detection. The statistic is demonstrated to be excellent in terms of its robustness and power performance, and useful for cluster analysis in settings of missing data. The test statistic can also be generalized to any one dimensional system where every site is observed, such as DNA sequencing data. The paper illustrates how the informal graphical tools for detecting clusters may be misleading. The statistic is used for finding clusters of putative SNP differences in a mixture of different mouse genetic backgrounds and clusters of de novo SNP differences arising between tissues with development and carcinogenesis.

Human APOBEC3B interacts with the heterogenous nuclear ribonucleoprotein A3 in cancer cells

Human APOBEC3B (A3B), like other APOBEC3 members, is a cytosine deaminase which causes hypermutation of single stranded genome. Recent studies have shown that A3B is predominantly elevated in multiple cancer tissues and cell lines such as the bladder, cervix, lung, head and neck, and breast. Upregulation and activation of A3B in developing tumors can cause an unexpected cluster of mutations which promote cancer development and progression. The cellular proteins which facilitate A3B function through direct or indirect interactions remain largely unknown. In this study, we performed LC-MS-based proteomics to identify cellular proteins which coimmunoprecipitated with A3B. Our results indicated a specific interaction of A3B with hnRNP A3 (heterogeneous nuclear ribonucleoprotein). This interaction was verified by co-immunoprecipitation and was found to be RNA-dependent. Furthermore, A3B and hnRNP A3 colocalized as evident from immunofluorescence analysis.

Enzyme cycling contributes to efficient induction of genome mutagenesis by the cytidine deaminase APOBEC3B

The single-stranded DNA cytidine deaminases APOBEC3B, APOBEC3H haplotype I, and APOBEC3A can contribute to cancer through deamination of cytosine to form promutagenic uracil in genomic DNA. The enzymes must access single-stranded DNA during the dynamic processes of DNA replication or transcription, but the enzymatic mechanisms enabling this activity are not known. To study this, we developed a method to purify full length APOBEC3B and characterized it in comparison to APOBEC3A and APOBEC3H on substrates relevant to cancer mutagenesis. We found that the ability of an APOBEC3 to cycle between DNA substrates determined whether it was able to efficiently deaminate single-stranded DNA produced by replication and single-stranded DNA bound by replication protein A (RPA). APOBEC3 deaminase activity during transcription had a size limitation that inhibited APOBEC3B tetramers, but not APOBEC3A monomers or APOBEC3H dimers. Altogether, the data support a model in which the availability of single-stranded DNA is necessary, but alone not sufficient for APOBEC3-induced mutagenesis in cells because there is also a dependence on the inherent biochemical properties of the enzymes. The biochemical properties identified in this study can be used to measure the mutagenic potential of other APOBEC enzymes in the genome.

p53 controls expression of the DNA deaminase APOBEC3B to limit its potential mutagenic activity in cancer cells

Cancer genome sequencing has implicated the cytosine deaminase activity of apolipoprotein B mRNA editing enzyme catalytic polypeptide-like (APOBEC) genes as an important source of mutations in diverse cancers, with APOBEC3B (A3B) expression especially correlated with such cancer mutations. To better understand the processes directing A3B over-expression in cancer, and possible therapeutic avenues for targeting A3B, we have investigated the regulation of A3B gene expression. Here, we show that A3B expression is inversely related to p53 status in different cancer types and demonstrate that this is due to a direct and pivotal role for p53 in repressing A3B expression. This occurs through the induction of p21 (CDKN1A) and the recruitment of the repressive DREAM complex to the A3B gene promoter, such that loss of p53 through mutation, or human papilloma virus-mediated inhibition, prevents recruitment of the complex, thereby causing elevated A3B expression and cytosine deaminase activity in cancer cells. As p53 is frequently mutated in cancer, our findings provide a mechanism by which p53 loss can promote cancer mutagenesis.