Nuri A. Temiz

APOBEC3B is an enzymatic source of mutation in breast cancer

Several mutations are required for cancer development, and genome sequencing has revealed that many cancers, including breast cancer, have somatic mutation spectra dominated by C-to-T transitions. Most of these mutations occur at hydrolytically disfavoured non-methylated cytosines throughout the genome, and are sometimes clustered. Here we show that the DNA cytosine deaminase APOBEC3B is a probable source of these mutations. APOBEC3B messenger RNA is upregulated in most primary breast tumours and breast cancer cell lines. Tumours that express high levels of APOBEC3B have twice as many mutations as those that express low levels and are more likely to have mutations in TP53. Endogenous APOBEC3B protein is predominantly nuclear and the only detectable source of DNA C-to-U editing activity in breast cancer cell-line extracts. Knockdown experiments show that endogenous APOBEC3B correlates with increased levels of genomic uracil, increased mutation frequencies, and C-to-T transitions. Furthermore, induced APOBEC3B overexpression causes cell cycle deviations, cell death, DNA fragmentation, γ-H2AX accumulation and C-to-T mutations. Our data suggest a model in which APOBEC3B-catalysed deamination provides a chronic source of DNA damage in breast cancers that could select TP53 inactivation and explain how some tumours evolve rapidly and manifest heterogeneity.

Evidence for APOBEC3B mutagenesis in multiple human cancers

Thousands of somatic mutations accrue in most human cancers, and their causes are largely unknown. We recently showed that the DNA cytidine deaminase APOBEC3B accounts for up to half of the mutational load in breast carcinomas expressing this enzyme. Here we address whether APOBEC3B is broadly responsible for mutagenesis in multiple tumor types. We analyzed gene expression data and mutation patterns, distributions and loads for 19 different cancer types, with over 4,800 exomes and 1,000,000 somatic mutations. Notably, APOBEC3B is upregulated, and its preferred target sequence is frequently mutated and clustered in at least six distinct cancers: bladder, cervix, lung (adenocarcinoma and squamous cell carcinoma), head and neck, and breast. Interpreting these findings in the light of previous genetic, cellular and biochemical studies, the most parsimonious conclusion from these global analyses is that APOBEC3B-catalyzed genomic uracil lesions are responsible for a large proportion of both dispersed and clustered mutations in multiple distinct cancers.

A Sleeping Beauty forward genetic screen identifies new genes and pathways driving osteosarcoma development and metastasis

Osteosarcomas are sarcomas of the bone, derived from osteoblasts or their precursors, with a high propensity to metastasize. Osteosarcoma is associated with massive genomic instability, making it problematic to identify driver genes using human tumors or prototypical mouse models, many of which involve loss of Trp53 function. To identify the genes driving osteosarcoma development and metastasis, we performed a Sleeping Beauty (SB) transposon-based forward genetic screen in mice with and without somatic loss of Trp53. Common insertion site (CIS) analysis of 119 primary tumors and 134 metastatic nodules identified 232 sites associated with osteosarcoma development and 43 sites associated with metastasis, respectively. Analysis of CIS-associated genes identified numerous known and new osteosarcoma-associated genes enriched in the ErbB, PI3K-AKT-mTOR and MAPK signaling pathways. Lastly, we identified several oncogenes involved in axon guidance, including Sema4d and Sema6d, which we functionally validated as oncogenes in human osteosarcoma.

Human Papillomavirus E6 Triggers Upregulation of the Antiviral and Cancer Genomic DNA Deaminase APOBEC3B

ABSTRACT Several recent studies have converged upon the innate immune DNA cytosine deaminase APOBEC3B (A3B) as a significant source of genomic uracil lesions and mutagenesis in multiple human cancers, including those of the breast, head/neck, cervix, bladder, lung, ovary, and other tissues. A3B is upregulated in these tumor types relative to normal tissues, but the mechanism is unclear. Because A3B also has antiviral activity in multiple systems and is a member of the broader innate immune response, we tested the hypothesis that human papillomavirus (HPV) infection causes A3B upregulation. We found that A3B mRNA expression and enzymatic activity were upregulated following transfection of a high-risk HPV genome and that this effect was abrogated by inactivation of E6. Transduction experiments showed that the E6 oncoprotein alone was sufficient to cause A3B upregulation, and a panel of high-risk E6 proteins triggered higher A3B levels than did a panel of low-risk or noncancer E6 proteins. Knockdown experiments in HPV-positive cell lines showed that endogenous E6 is required for A3B upregulation. Analyses of publicly available head/neck cancer data further support this relationship, as A3B levels are higher in HPV-positive cancers than in HPV-negative cancers. Taken together with the established role for high-risk E6 in functional inactivation of TP53 and published positive correlations in breast cancer between A3B upregulation and genetic inactivation of TP53, our studies suggest a model in which high-risk HPV E6, possibly through functional inactivation of TP53, causes derepression of A3B gene transcription. This would lead to a mutator phenotype that explains the observed cytosine mutation biases in HPV-positive head/neck and cervical cancers. IMPORTANCE The innate immune DNA cytosine deaminase APOBEC3B (A3B) accounts for a large proportion of somatic mutations in cervical and head/neck cancers, but nothing is known about the mechanism responsible for its upregulation in these tumor types. Almost all cervical carcinomas and large proportions of head/neck tumors are caused by human papillomavirus (HPV) infection. Here, we establish a mechanistic link between HPV infection and A3B upregulation. The E6 oncoprotein of high-risk, but not low-risk, HPV types triggers A3B upregulation, supporting a model in which TP53 inactivation causes a derepression of A3B gene transcription and elevated A3B enzyme levels. This virus-induced mutator phenotype provides a mechanistic explanation for A3B signature mutations observed in HPV-positive head/neck and cervical carcinomas and may also help to account for the preferential cancer predisposition caused by high-risk HPV isolates. The innate immune DNA cytosine deaminase APOBEC3B (A3B) accounts for a large proportion of somatic mutations in cervical and head/neck cancers, but nothing is known about the mechanism responsible for its upregulation in these tumor types. Almost all cervical carcinomas and large proportions of head/neck tumors are caused by human papillomavirus (HPV) infection. Here, we establish a mechanistic link between HPV infection and A3B upregulation. The E6 oncoprotein of high-risk, but not low-risk, HPV types triggers A3B upregulation, supporting a model in which TP53 inactivation causes a derepression of A3B gene transcription and elevated A3B enzyme levels. This virus-induced mutator phenotype provides a mechanistic explanation for A3B signature mutations observed in HPV-positive head/neck and cervical carcinomas and may also help to account for the preferential cancer predisposition caused by high-risk HPV isolates.

Non-B DB v2.0: a database of predicted non-B DNA-forming motifs and its associated tools

The non-B DB, available at http://nonb.abcc.ncifcrf.gov, catalogs predicted non-B DNA-forming sequence motifs, including Z-DNA, G-quadruplex, A-phased repeats, inverted repeats, mirror repeats, direct repeats and their corresponding subsets: cruciforms, triplexes and slipped structures, in several genomes. Version 2.0 of the database revises and re-implements the motif discovery algorithms to better align with accepted definitions and thresholds for motifs, expands the non-B DNA-forming motifs coverage by including short tandem repeats and adds key visualization tools to compare motif locations relative to other genomic annotations. Non-B DB v2.0 extends the ability for comparative genomics by including re-annotation of the five organisms reported in non-B DB v1.0, human, chimpanzee, dog, macaque and mouse, and adds seven additional organisms: orangutan, rat, cow, pig, horse, platypus and Arabidopsis thaliana. Additionally, the non-B DB v2.0 provides an overall improved graphical user interface and faster query performance.

The role of methylation in the intrinsic dynamics of B- and Z-DNA.

Methylation of cytosine at the 5-carbon position (5mC) is observed in both prokaryotes and eukaryotes. In humans, DNA methylation at CpG sites plays an important role in gene regulation and has been implicated in development, gene silencing, and cancer. In addition, the CpG dinucleotide is a known hot spot for pathologic mutations genome-wide. CpG tracts may adopt left-handed Z-DNA conformations, which have also been implicated in gene regulation and genomic instability. Methylation facilitates this B-Z transition but the underlying mechanism remains unclear. Herein, four structural models of the dinucleotide d(GC)5 repeat sequence in B-, methylated B-, Z-, and methylated Z-DNA forms were constructed and an aggregate 100 nanoseconds of molecular dynamics simulations in explicit solvent under physiological conditions was performed for each model. Both unmethylated and methylated B-DNA were found to be more flexible than Z-DNA. However, methylation significantly destabilized the BII, relative to the BI, state through the Gp5mC steps. In addition, methylation decreased the free energy difference between B- and Z-DNA. Comparisons of α/γ backbone torsional angles showed that torsional states changed marginally upon methylation for B-DNA, and Z-DNA. Methylation-induced conformational changes and lower energy differences may contribute to the transition to Z-DNA by methylated, over unmethylated, B-DNA and may be a contributing factor to biological function.

The DNA cytosine deaminase APOBEC3B promotes tamoxifen resistance in ER-positive breast cancer

An antiviral enzyme promotes drug resistance in breast cancer. Breast tumors often display extreme genetic heterogeneity characterized by hundreds of gross chromosomal aberrations and tens of thousands of somatic mutations. Tumor evolution is thought to be ongoing and driven by multiple mutagenic processes. A major outstanding question is whether primary tumors have preexisting mutations for therapy resistance or whether additional DNA damage and mutagenesis are necessary. Drug resistance is a key measure of tumor evolvability. If a resistance mutation preexists at the time of primary tumor presentation, then the intended therapy is likely to fail. However, if resistance does not preexist, then ongoing mutational processes still have the potential to undermine therapeutic efficacy. The antiviral enzyme APOBEC3B (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3B) preferentially deaminates DNA C-to-U, which results in signature C-to-T and C-to-G mutations commonly observed in breast tumors. We use clinical data and xenograft experiments to ask whether APOBEC3B contributes to ongoing breast tumor evolution and resistance to the selective estrogen receptor modulator, tamoxifen. First, APOBEC3B levels in primary estrogen receptor–positive (ER+) breast tumors inversely correlate with the clinical benefit of tamoxifen in the treatment of metastatic ER+ disease. Second, APOBEC3B depletion in an ER+ breast cancer cell line results in prolonged tamoxifen responses in murine xenograft experiments. Third, APOBEC3B overexpression accelerates the development of tamoxifen resistance in murine xenograft experiments by a mechanism that requires the enzyme’s catalytic activity. These studies combine to indicate that APOBEC3B promotes drug resistance in breast cancer and that inhibiting APOBEC3B-dependent tumor evolvability may be an effective strategy to improve efficacies of targeted cancer therapies.

Guanine Holes Are Prominent Targets for Mutation in Cancer and Inherited Disease

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G•C bp in the context of all 64 5′-NGNN-3′ motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

RNA sequencing of Sleeping Beauty transposon-induced tumors detects transposon-RNA fusions in forward genetic cancer screens

Forward genetic screens using Sleeping Beauty (SB)-mobilized T2/Onc transposons have been used to identify common insertion sites (CISs) associated with tumor formation. Recurrent sites of transposon insertion are commonly identified using ligation-mediated PCR (LM-PCR). Here, we use RNA sequencing (RNA-seq) data to directly identify transcriptional events mediated by T2/Onc. Surprisingly, the majority (∼80%) of LM-PCR identified junction fragments do not lead to observable changes in RNA transcripts. However, in CIS regions, direct transcriptional effects of transposon insertions are observed. We developed an automated method to systematically identify T2/Onc-genome RNA fusion sequences in RNA-seq data. RNA fusion-based CISs were identified corresponding to both DNA-based CISs (Cdkn2a, Mycl1, Nf2, Pten, Sema6d, and Rere) and additional regions strongly associated with cancer that were not observed by LM-PCR (Myc, Akt1, Pth, Csf1r, Fgfr2, Wisp1, Map3k5, and Map4k3). In addition to calculating recurrent CISs, we also present complementary methods to identify potential driver events via determination of strongly supported fusions and fusions with large transcript level changes in the absence of multitumor recurrence. These methods independently identify CIS regions and also point to cancer-associated genes like Braf. We anticipate RNA-seq analyses of tumors from forward genetic screens will become an efficient tool to identify causal events.

Searching for Non‐B DNA‐Forming Motifs Using nBMST (Non‐B DNA Motif Search Tool)

This unit describes basic protocols on using the non‐B DNA Motif Search Tool (nBMST) to search for sequence motifs predicted to form alternative DNA conformations that differ from the canonical right‐handed Watson‐Crick double‐helix, collectively known as non‐B DNA, and on using the associated PolyBrowse, a GBrowse–based genomic browser. The nBMST is a Web‐based resource that allows users to submit one or more DNA sequences to search for inverted repeats (cruciform DNA), mirror repeats (triplex DNA), direct/tandem repeats (slipped/hairpin structures), G4 motifs (tetraplex, G‐quadruplex DNA), alternating purine‐pyrimidine tracts (left‐handed Z‐DNA), and A‐phased repeats (static bending). The nBMST is versatile, simple to use, does not require bioinformatics skills, and can be applied to any type of DNA sequences, including viral and bacterial genomes, up to an aggregate of 20 megabasepairs (Mbp). Curr. Protoc. Hum. Genet. 73:18.7.1‐18.7.22.