APOBEC shapes tumor evolution and age at onset of lung cancer in smokers

APOBEC enzymes are part of the innate immunity and are responsible for restricting viruses and retroelements by deaminating cytosine residues1,2. Most solid tumors harbor different levels of somatic mutations attributed to the off-target activities of APOBEC3A (A3A) and/or APOBEC3B (A3B)3-6. However, how APOBEC3A/B enzymes shape the tumor evolution in the presence of exogenous mutagenic processes is largely unknown. Here, by combining deep whole-genome sequencing with multi-omics profiling of 309 lung cancers from smokers with detailed tobacco smoking information, we identify two subtypes defined by low (LAS) and high (HAS) APOBEC mutagenesis. LAS are enriched for A3B-like mutagenesis and KRAS mutations, whereas HAS for A3A-like mutagenesis and TP53 mutations. Unlike APOBEC3A, APOBEC3B expression is strongly associated with an upregulation of the base excision repair pathway. Hypermutation by unrepaired A3A and tobacco smoking mutagenesis combined with TP53-induced genomic instability can trigger senescence7, apoptosis8, and cell regeneration9, as indicated by high expression of pulmonary healing signaling pathway, stemness markers and distal cell-of-origin in HAS. The expected association of tobacco smoking variables (e.g., time to first cigarette) with genomic/epigenomic changes are not observed in HAS, a plausible consequence of frequent cell senescence or apoptosis. HAS have more neoantigens, slower clonal expansion, and older age at onset compared to LAS, particularly in heavy smokers, consistent with high proportions of newly generated, unmutated cells and frequent immuno-editing. These findings show how heterogeneity in mutational burden across co-occurring mutational processes and cell types contributes to tumor development, with important clinical implications.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

Changes in metabolic landscapes shape divergent but distinct mutational signatures and cytotoxic consequences of redox stress

Mutational signatures discerned in cancer genomes, in aging tissues and in cells exposed to toxic agents, reflect complex processes underlying transformation of cells from normal to dysfunctional. Due to its ubiquitous and chronic nature, redox stress contributions to cellular makeover remain equivocal. The deciphering of a new mutational signature of an environmentally-relevant oxidizing agent, potassium bromate, in yeast single strand DNA uncovered a surprising heterogeneity in the mutational signatures of oxidizing agents. NMR-based analysis of molecular outcomes of redox stress revealed profound dissimilarities in metabolic landscapes following exposure to hydrogen peroxide versus potassium bromate. The predominance of G to T substitutions in the mutational spectra distinguished potassium bromate from hydrogen peroxide and paraquat and mirrored the observed metabolic changes. We attributed these changes to the generation of uncommon oxidizing species in a reaction with thiol-containing antioxidants; a nearly total depletion of intracellular glutathione and a paradoxical augmentation of potassium bromate mutagenicity and toxicity by antioxidants. Our study provides the framework for understanding multidimensional processes triggered by agents collectively known as oxidants. Detection of increased mutational loads associated with potassium bromate-related mutational motifs in human tumors may be clinically relevant as a biomarker of this distinct type of redox stress.

Genomic and evolutionary classification of lung cancer in never smokers

Lung cancer in never smokers (LCINS) is a common cause of cancer mortality but its genomic landscape is poorly characterized. Here high-coverage whole-genome sequencing of 232 LCINS showed 3 subtypes defined by copy number aberrations. The dominant subtype (piano), which is rare in lung cancer in smokers, features somatic UBA1 mutations, germline AR variants and stem cell-like properties, including low mutational burden, high intratumor heterogeneity, long telomeres, frequent KRAS mutations and slow growth, as suggested by the occurrence of cancer drivers' progenitor cells many years before tumor diagnosis. The other subtypes are characterized by specific amplifications and EGFR mutations (mezzo-forte) and whole-genome doubling (forte). No strong tobacco smoking signatures were detected, even in cases with exposure to secondhand tobacco smoke. Genes within the receptor tyrosine kinase-Ras pathway had distinct impacts on survival; five genomic alterations independently doubled mortality. These findings create avenues for personalized treatment in LCINS.

Mixture Analyses of Air-sampled Pollen Extracts Can Accurately Differentiate Pollen Taxa

The daily pollen forecast provides crucial information for allergic patients to avoid exposure to specific pollen. Pollen counts are typically measured with air samplers and analyzed with microscopy by trained experts. In contrast, this study evaluated the effectiveness of identifying the component pollens using the metabolites extracted from an air-sampled pollen mixture. Ambient air-sampled pollen from Munich in 2016 and 2017 was visually identified from reference pollens and extracts were prepared. The extracts were lyophilized, rehydrated in optimal NMR buffers, and filtered to remove large proteins. NMR spectra were analyzed for pollen associated metabolites. Regression and decision-tree based algorithms using the concentration of metabolites, calculated from the NMR spectra outperformed algorithms using the NMR spectra themselves as input data for pollen identification. Categorical prediction algorithms trained for low, medium, high, and very high pollen count groups had accuracies of 74% for the tree, 82% for the grass, and 93% for the weed pollen count. Deep learning models using convolutional neural networks performed better than regression models using NMR spectral input, and were the overall best method in terms of relative error and classification accuracy (86% for tree, 89% for grass, and 93% for weed pollen count). This study demonstrates that NMR spectra of air-sampled pollen extracts can be used in an automated fashion to provide taxa and type-specific measures of the daily pollen count.

Somatic mutation distributions in cancer genomes vary with three-dimensional chromatin structure

Somatic mutations in driver genes may ultimately lead to the development of cancer. Understanding how somatic mutations accumulate in cancer genomes and the underlying factors that generate somatic mutations is therefore crucial for developing novel therapeutic strategies. To understand the interplay between spatial genome organization and specific mutational processes, we studied 3,000 tumor-normal-pair whole-genome datasets from 42 different human cancer types. Our analyses reveal that the change in somatic mutational load in cancer genomes is co-localized with topologically-associating-domain boundaries. Domain boundaries constitute a better proxy to track mutational load change than replication timing measurements. We show that different mutational processes lead to distinct somatic mutation distributions where certain processes generate mutations in active domains, and others generate mutations in inactive domains. Overall, the interplay between three-dimensional genome organization and active mutational processes has a substantial influence on the large-scale mutation-rate variations observed in human cancers.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella-another positive single stranded (ss) RNA virus. Each of the rubella virus isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella virus, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2-likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

APOBEC3-dependent kataegis and TREX1-driven chromothripsis during telomere crisis

Chromothripsis and kataegis are frequently observed in cancer and may arise from telomere crisis, a period of genome instability during tumorigenesis when depletion of the telomere reserve generates unstable dicentric chromosomes1-5. Here we examine the mechanism underlying chromothripsis and kataegis by using an in vitro telomere crisis model. We show that the cytoplasmic exonuclease TREX1, which promotes the resolution of dicentric chromosomes4, plays a prominent role in chromothriptic fragmentation. In the absence of TREX1, the genome alterations induced by telomere crisis primarily involve breakage-fusion-bridge cycles and simple genome rearrangements rather than chromothripsis. Furthermore, we show that the kataegis observed at chromothriptic breakpoints is the consequence of cytosine deamination by APOBEC3B. These data reveal that chromothripsis and kataegis arise from a combination of nucleolytic processing by TREX1 and cytosine editing by APOBEC3B.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella - another positive single stranded (ss) RNA virus. Each of the rubella isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2 - likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

Analysis of APOBEC-induced mutations in yeast strains with low levels of replicative DNA polymerases

Yeast strains with low levels of the replicative DNA polymerases (alpha, delta, and epsilon) have high levels of chromosome deletions, duplications, and translocations. By examining the patterns of mutations induced in strains with low levels of DNA polymerase by the human protein APOBEC3B (a protein that deaminates cytosine in single-stranded DNA), we show dramatically elevated amounts of single-stranded DNA relative to a wild-type strain. During DNA replication, one strand (defined as the leading strand) is replicated processively by DNA polymerase epsilon and the other (the lagging strand) is replicated as short fragments initiated by DNA polymerase alpha and extended by DNA polymerase delta. In the low DNA polymerase alpha and delta strains, the APOBEC-induced mutations are concentrated on the lagging-strand template, whereas in the low DNA polymerase epsilon strain, mutations occur on the leading- and lagging-strand templates with similar frequencies. In addition, for most genes, the transcribed strand is mutagenized more frequently than the nontranscribed strand. Lastly, some of the APOBEC-induced clusters in strains with low levels of DNA polymerase alpha or delta are greater than 10 kb in length.

Mutation signatures specific to DNA alkylating agents in yeast and cancers

Alkylation is one of the most ubiquitous forms of DNA lesions. However, the motif preferences and substrates for the activity of the major types of alkylating agents defined by their nucleophilic substitution reactions (SN1 and SN2) are still unclear. Utilizing yeast strains engineered for large-scale production of single-stranded DNA (ssDNA), we probed the substrate specificity, mutation spectra and signatures associated with DNA alkylating agents. We determined that SN1-type agents preferably mutagenize double-stranded DNA (dsDNA), and the mutation signature characteristic of the activity of SN1-type agents was conserved across yeast, mice and human cancers. Conversely, SN2-type agents preferably mutagenize ssDNA in yeast. Moreover, the spectra and signatures derived from yeast were detectable in lung cancers, head and neck cancers and tumors from patients exposed to SN2-type alkylating chemicals. The estimates of mutation loads associated with the SN2-type alkylation signature were higher in lung tumors from smokers than never-smokers, pointing toward the mutagenic activity of the SN2-type alkylating carcinogens in cigarettes. In summary, our analysis of mutations in yeast strains treated with alkylating agents, as well as in whole-exome and whole-genome-sequenced tumors identified signatures highly specific to alkylation mutagenesis and indicate the pervasive nature of alkylation-induced mutagenesis in cancers.

Comprehensive analysis of chromothripsis in 2,658 human cancers using whole-genome sequencing

Chromothripsis is a mutational phenomenon characterized by massive, clustered genomic rearrangements that occurs in cancer and other diseases. Recent studies in selected cancer types have suggested that chromothripsis may be more common than initially inferred from low-resolution copy-number data. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we analyze patterns of chromothripsis across 2,658 tumors from 38 cancer types using whole-genome sequencing data. We find that chromothripsis events are pervasive across cancers, with a frequency of more than 50% in several cancer types. Whereas canonical chromothripsis profiles display oscillations between two copy-number states, a considerable fraction of events involve multiple chromosomes and additional structural alterations. In addition to non-homologous end joining, we detect signatures of replication-associated processes and templated insertions. Chromothripsis contributes to oncogene amplification and to inactivation of genes such as mismatch-repair-related genes. These findings show that chromothripsis is a major process that drives genome evolution in human cancer.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

The repertoire of mutational signatures in human cancer

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature¹. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium² of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses^3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Infectious vaccine-derived rubella viruses emerge, persist, and evolve in cutaneous granulomas of children with primary immunodeficiencies

Rubella viruses (RV) have been found in an association with granulomas in children with primary immune deficiencies (PID). Here, we report the recovery and characterization of infectious immunodeficiency-related vaccine-derived rubella viruses (iVDRV) from diagnostic skin biopsies of four patients. Sequence evolution within PID hosts was studied by comparison of the complete genomic sequences of the iVDRVs with the genome of the vaccine virus RA27/3. The degree of divergence of each iVDRV correlated with the duration of persistence indicating continuous intrahost evolution. The evolution rates for synonymous and nonsynonymous substitutions were estimated to be 5.7 x 10-3 subs/site/year and 8.9 x 10-4 subs/site/year, respectively. Mutational spectra and signatures indicated a major role for APOBEC cytidine deaminases and a secondary role for ADAR adenosine deaminases in generating diversity of iVDRVs. The distributions of mutations across the genes and 3D hotspots for amino acid substitutions in the E1 glycoprotein identified regions that may be under positive selective pressure. Quasispecies diversity was higher in granulomas than in recovered infectious iVDRVs. Growth properties of iVDRVs were assessed in WI-38 fibroblast cultures. None of the iVDRV isolates showed complete reversion to wild type phenotype but the replicative and persistence characteristics of iVDRVs were different from those of the RA27/3 vaccine strain, making predictions of iVDRV transmissibility and teratogenicity difficult. However, detection of iVDRV RNA in nasopharyngeal specimen and poor neutralization of some iVDRV strains by sera from vaccinated persons suggests possible public health risks associated with iVDRV carriers. Detection of IgM antibody to RV in sera of two out of three patients may be a marker of virus persistence, potentially useful for identifying patients with iVDRV before development of lesions. Studies of the evolutionary dynamics of iVDRV during persistence will contribute to development of infection control strategies and antiviral therapies.

Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation

A single cancer genome can harbor thousands of clustered mutations. Mutation signature analyses have revealed that the origin of clusters are lesions in long tracts of single-stranded (ss) DNA damaged by apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) cytidine deaminases, raising questions about molecular mechanisms that generate long ssDNA vulnerable to hypermutation. Here, we show that ssDNA intermediates formed during the repair of gamma-induced bursts of double-strand breaks (DSBs) in the presence of APOBEC3A in yeast lead to multiple APOBEC-induced clusters similar to cancer. We identified three independent pathways enabling cluster formation associated with repairing bursts of DSBs: 5′ to 3′ bidirectional resection, unidirectional resection, and break-induced replication (BIR). Analysis of millions of mutations in APOBEC-hypermutated cancer genomes revealed that cancer tolerance to formation of hypermutable ssDNA is similar to yeast and that the predominant pattern of clustered mutagenesis is the same as in resection-defective yeast, suggesting that cluster formation in cancers is driven by a BIR-like mechanism. The phenomenon of genome-wide burst of clustered mutagenesis revealed by our study can play an important role in generating somatic hypermutation in cancers as well as in noncancerous cells.

This study uses yeast expressing a human cytidine deaminase to reveal simultaneous stretches of long single-strand DNA and multiple vast mutation clusters in a single eukaryotic cell repairing multiple double-strand breaks. This is reminiscent of the phenomenon of “kataegis” or hypermutation observed in cancer genomes, suggesting that a similar mechanism is involved.

Mutational signatures of redox stress in yeast single-strand DNA and of aging in human mitochondrial DNA share a common feature

Redox stress is a major hallmark of cancer. Analysis of thousands of sequenced cancer exomes and whole genomes revealed distinct mutational signatures that can be attributed to specific sources of DNA lesions. Clustered mutations discovered in several cancer genomes were linked to single-strand DNA (ssDNA) intermediates in various processes of DNA metabolism. Previously, only one clustered mutational signature had been clearly associated with a subclass of ssDNA-specific apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) cytidine deaminases. Others remain to be elucidated. We report here deciphering of the mutational spectra and mutational signature of redox stress in ssDNA of budding yeast and the signature of aging in human mitochondrial DNA. We found that the predominance of C to T substitutions is a common feature of both signatures. Measurements of the frequencies of hydrogen peroxide-induced mutations in proofreading-defective yeast mutants supported the conclusion that hydrogen peroxide-induced mutagenesis is not the result of increased DNA polymerase misincorporation errors but rather is caused by direct damage to DNA. Proteins involved in modulation of chromatin status play a significant role in prevention of redox stress-induced mutagenesis, possibly by facilitating protection through modification of chromatin structure. These findings provide an opportunity for the search and identification of the mutational signature of redox stress in cancers and in other pathological conditions and could potentially be used for informing therapeutic decisions. In addition, the discovery of such signatures that may be present in related organisms should also advance our understanding of evolution.

Molecular Evolution of Early-Onset Prostate Cancer Identifies Molecular Risk Markers and Clinical Trajectories

Identifying the earliest somatic changes in prostate cancer can give important insights into tumor evolution and aids in stratifying high- from low-risk disease. We integrated whole genome, transcriptome and methylome analysis of early-onset prostate cancers (diagnosis ≤55 years). Characterization across 292 prostate cancer genomes revealed age-related genomic alterations and a clock-like enzymatic-driven mutational process contributing to the earliest mutations in prostate cancer patients. Our integrative analysis identified four molecular subgroups, including a particularly aggressive subgroup with recurrent duplications associated with increased expression of ESRP1, which we validate in 12,000 tissue microarray tumors. Finally, we combined the patterns of molecular co-occurrence and risk-based subgroup information to deconvolve the molecular and clinical trajectories of prostate cancer from single patient samples.

Abstract LB-378: Comprehensive analysis of chromothripsis in 2,658 human cancers using whole-genome sequencing

Chromothripsis is a newly discovered mutational phenomenon involving massive, clustered genomic rearrangements that occurs in cancer and other diseases. Recent studies in cancer suggest that chromothripsis may be far more common than initially inferred from low resolution DNA copy number data. Here, we analyze the patterns of chromothripsis across 2,658 tumors spanning 39 cancer types using whole-genome sequencing data. We find that chromothripsis events are pervasive across cancers, with a frequency of >50% in several cancer types. Whereas canonical chromothripsis profiles display oscillations between two copy number states, a considerable fraction of the events involves multiple chromosomes as well as additional structural alterations. In addition to non-homologous end-joining, we detect signatures of replicative processes and templated insertions. Chromothripsis contributes to oncogene amplification as well as to inactivation of genes such as mismatch-repair related genes. These findings show that chromothripsis is a major process driving genome evolution in human cancer.

Citation Format: Isidro Cortés-Ciriano, June-Koo Lee, Ruibin Xi, Dhawal Jain, Youngsook L. Jung, Lixing Yang, Dmitry Gordenin, Leszek J. Klimczak, Cheng-Zhong Zhang, David S. Pellman, Peter J. Park. Comprehensive analysis of chromothripsis in 2,658 human cancers using whole-genome sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr LB-378.

Genomic Subtypes of Non-invasive Bladder Cancer with Distinct Metabolic Profile and Female Gender Bias in KDM6A Mutation Frequency

Bladder cancer incurs a higher lifetime treatment cost than other cancers due to frequent recurrence of non-invasive disease. Improved prognostic biomarkers and localized therapy are needed for this large patient group. We defined two major genomic subtypes of primary stage Ta tumors. One of these was characterized by loss of 9q including TSC1, increased KI67 labeling index, upregulated glycolysis, DNA repair, mTORC1 signaling, features of the unfolded protein response, and altered cholesterol homeostasis. Comparison with muscle-invasive bladder cancer mutation profiles revealed lower overall mutation rates and more frequent mutations in RHOB and chromatin modifier genes. More mutations in the histone lysine demethylase KDM6A were present in non-invasive tumors from females than males.

Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer

We report a comprehensive analysis of 412 muscle-invasive bladder cancers characterized by multiple TCGA analytical platforms. Fifty-eight genes were significantly mutated, and the overall mutational load was associated with APOBEC-signature mutagenesis. Clustering by mutation signature identified a high-mutation subset with 75% 5-year survival. mRNA expression clustering refined prior clustering analyses and identified a poor-survival "neuronal" subtype in which the majority of tumors lacked small cell or neuroendocrine histology. Clustering by mRNA, long non-coding RNA (lncRNA), and miRNA expression converged to identify subsets with differential epithelial-mesenchymal transition status, carcinoma in situ scores, histologic features, and survival. Our analyses identified 5 expression subtypes that may stratify response to different treatments.

The Impact of Environmental and Endogenous Damage on Somatic Mutation Load in Human Skin Fibroblasts

Accumulation of somatic changes, due to environmental and endogenous lesions, in the human genome is associated with aging and cancer. Understanding the impacts of these processes on mutagenesis is fundamental to understanding the etiology, and improving the prognosis and prevention of cancers and other genetic diseases. Previous methods relying on either the generation of induced pluripotent stem cells, or sequencing of single-cell genomes were inherently error-prone and did not allow independent validation of the mutations. In the current study we eliminated these potential sources of error by high coverage genome sequencing of single-cell derived clonal fibroblast lineages, obtained after minimal propagation in culture, prepared from skin biopsies of two healthy adult humans. We report here accurate measurement of genome-wide magnitude and spectra of mutations accrued in skin fibroblasts of healthy adult humans. We found that every cell contains at least one chromosomal rearrangement and 600–13,000 base substitutions. The spectra and correlation of base substitutions with epigenomic features resemble many cancers. Moreover, because biopsies were taken from body parts differing by sun exposure, we can delineate the precise contributions of environmental and endogenous factors to the accrual of genetic changes within the same individual. We show here that UV-induced and endogenous DNA damage can have a comparable impact on the somatic mutation loads in skin fibroblasts.

Somatic genomes are constantly accumulating changes caused by endogenous lesions, errors in DNA replication and repair, as well as environmental insults. Despite the importance of somatic genome instability in aging and age-related pathologies, including cancers, accurate measurements of mutation loads in healthy cells is still missing. In this study, we developed an experimental approach to accurately determine the somatic genome changes accrued in cell lineages over the lifetime of healthy humans. We show that the amounts and types of mutations in skin cells resemble many cancers, thus indicating that the mechanisms that lead to carcinogenesis are also functional in healthy cells. Moreover, sun-exposed skin cells have a higher mutation load attributable to ultraviolet radiation (UV) unlike cells from hips that were protected by clothing. Our work provides precise measurements of the mutation loads in single cells in human skin. Furthermore our data allowed defining the mutagenic impacts of environmental and endogenous processes within the same individual and led to conclusion that these processes have a comparable impact on the somatic mutation load.

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

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.

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.

Clustered mutations in yeast and in human cancers can arise from damaged long single-strand DNA regions

Mutations are typically perceived as random, independent events. We describe here nonrandom clustered mutations in yeast and in human cancers. Genome sequencing of yeast grown under chronic alkylation damage identified mutation clusters that extend up to 200 kb. A predominance of "strand-coordinated" changes of either cytosines or guanines in the same strand, mutation patterns, and genetic controls indicated that simultaneous mutations were generated by base alkylation in abnormally long single-strand DNA (ssDNA) formed at double-strand breaks (DSBs) and replication forks. Significantly, we found mutation clusters with analogous features in sequenced human cancers. Strand-coordinated clusters of mutated cytosines or guanines often resided near chromosome rearrangement breakpoints and were highly enriched with a motif targeted by APOBEC family cytosine-deaminases, which strongly prefer ssDNA. These data indicate that hypermutation via multiple simultaneous changes in randomly formed ssDNA is a general phenomenon that may be an important mechanism producing rapid genetic variation.

Detection of genome-wide polymorphisms in the AT-rich Plasmodium falciparum genome using a high-density microarray

Background

Genetic mapping is a powerful method to identify mutations that cause drug resistance and other phenotypic changes in the human malaria parasite Plasmodium falciparum. For efficient mapping of a target gene, it is often necessary to genotype a large number of polymorphic markers. Currently, a community effort is underway to collect single nucleotide polymorphisms (SNP) from the parasite genome. Here we evaluate polymorphism detection accuracy of a high-density 'tiling' microarray with 2.56 million probes by comparing single feature polymorphisms (SFP) calls from the microarray with known SNP among parasite isolates.

Results

We found that probe GC content, SNP position in a probe, probe coverage, and signal ratio cutoff values were important factors for accurate detection of SFP in the parasite genome. We established a set of SFP calling parameters that could predict mSFP (SFP called by multiple overlapping probes) with high accuracy (≥ 94%) and identified 121,087 mSFP genome-wide from five parasite isolates including 40,354 unique mSFP (excluding those from multi-gene families) and ~18,000 new mSFP, producing a genetic map with an average of one unique mSFP per 570 bp. Genomic copy number variation (CNV) among the parasites was also cataloged and compared.

Conclusion

A large number of mSFP were discovered from the P. falciparum genome using a high-density microarray, most of which were in clusters of highly polymorphic genes at chromosome ends. Our method for accurate mSFP detection and the mSFP identified will greatly facilitate large-scale studies of genome variation in the P. falciparum parasite and provide useful resources for mapping important parasite traits.

Novel directions in antipsychotic target identification using gene arrays

Schizophrenia is a major health problem that affects 2 million individuals in the United States. Antipsychotics offer considerable symptomatic relief and, although commonly discovered by screening with single biological targets, most interact with multiple receptors and signaling pathways. Considerable evidence from family and twin studies demonstrates genetic components and multiple chromosomal regions associated with schizophrenia. The polygenic nature of schizophrenia and multiple mechanisms for most effective agents indicate the need for broader approaches to target identification. Gene expression profiling of post-mortem human brain tissue simultaneously reveals the expression of many thousands of genes. A comparison of tissue from normals and patients provides a 'disease signature' of aberrantly expressed genes. 'Drug signatures' are the gene expression changes of cultured human or animal neurons treated with psychiatric drugs, and from animals chronically treated with these drugs. A selection of genes from disease and drug signatures can create a set of targets whose changes may better predict disease and its treatment by effective agents. This multi-parameter high throughput screening (MPHTS(SM)) approach evaluates the mRNA expression pattern of cultured cells exposed to candidate compounds. Compounds that normalize genes altered in schizophrenia may better address its underlying causes. Drugs that mimic gene expression changes that are consistently altered by effective antipsychotic agents provide a drug improvement strategy if efficacy is enhanced or side effects are attenuated.

Origin, targeting, and function of the apicomplexan plastid

The discovery of a plastid in Plasmodium, Toxoplasma and related protozoan parasites provides a satisfying resolution to several long-standing mysteries: the mechanism of action for various surprisingly effective antibiotics; the subcellular location of an enigmatic 35 kb episomal DNA; and the nature of an unusual intracellular structure containing multiple membranes. The apicomplexan plastid highlights the importance of lateral genetic transfer in evolution and provides an accessible system for the investigation of protein targeting to secondary endosymbiotic organelles. Combining molecular genetic identification of targeting signals with whole genome analysis promises to yield a complete picture of organellar metabolic pathways and new targets for drug design.

An enzyme similar to animal type II photolyases mediates photoreactivation in Arabidopsis

The important issue of photoreactivation DNA repair in plants has become even more interesting in recent years because a family of genes that are highly homologous to photoreactivating DNA repair enzymes but that function as blue light photoreceptors has been isolated. Here, we report the isolation of a novel photolyase-like sequence from Arabidopsis designated PHR1 (for photoreactivating enzyme). It shares little sequence similarity with either type I photolyases or the cryptochrome family of blue light photoreceptors. Instead, the PHR1 gene encodes an amino acid sequence with significant homology to the recently characterized type II photolyases identified in a number of prokaryotic and animal systems. PHR1 is a single-copy gene and is not expressed in dark-grown etiolated seedlings: the message is light inducible, which is similar to the expression profile for photoreactivation activity in plants. The PHR1 protein complements a photolyase-deficient mutant of Escherichia coli and thus confers photoreactivation activity. In addition, an Arabidopsis mutant that is entirely lacking in photolyase activity has been found to contain a lesion within this Arabidopsis type II photolyase sequence. We conclude that PHR1 represents a genuine plant photolyase gene and that the plant genes with homology to type I photolyases (the cryptochrome family of blue light photoreceptors) do not contribute to photoreactivation repair, at least in the case of Arabidopsis.

Multiple isoforms of Arabidopsis casein kinase I combine conserved catalytic domains with variable carboxyl-terminal extensions

Three cDNA clones encoding isoforms of casein kinase I (CKI) were isolated from Arabidopsis thaliana. One full-length clone, designated CKI1, contained an open reading frame of 1371 bp encoding a protein of 51,949 D with an isoelectric point of 9.7. In addition to the highly conserved catalytic domain (of about 300 amino acids), the Arabidopsis CKI isoforms contain 150 to 180 amino acid carboxyl-terminal extensions, which show among themselves a lower level of sequence conservation. These extensions do not show any sequence similarity to nonplant CKI isoforms, such as rat testis CKI delta, which is their closest isolated homolog, or to yeast CKI isoforms. Three additional isoforms of Arabidopsis CKI were found in the data bases of expressed sequence tags and/or were isolated serendipitously in nonspecific screening procedures by others. One of them also shows a carboxyl-terminal extension, but of only 80 amino acids. Casein kinase activity was detected in the soluble fraction of Escherichia coli strains expressing the CKI1 protein. This activity showed the crucial properties of CKI, including the ability to phosphorylate the D4 peptide, a specific substrate of CKI, and inhibition by N-(2-aminoethyl)-5-chloroisoquinoline-8-sulfonamide, a specific CKI inhibitor. Like several recombinant CKI isoforms from yeast, CKI1 was able to phosphorylate tyrosine-containing acidic polymers.

Reconstitution of Arabidopsis casein kinase II from recombinant subunits and phosphorylation of transcription factor GBF1

In contrast to the well-defined tetrameric structure of animal and yeast casein kinase II (CKII), plant CKII is found in two forms: a monomeric form and an oligomeric form whose subunit composition is not well defined. The Arabidopsis homologs of the catalytic subunit alpha (CKA1) and the regulatory subunit beta (CKB1) of CKII were expressed in Escherichia coli to examine their ability to form complexes, the effect of CKB1 on the catalytic activity, and the relationship of the recombinant enzymes to those isolated from plant material. Both subunits were found mainly in the inclusion body fraction in the bacterial expression strains, and they were solubilized and renatured with the recovery of catalytic (CKA1) and stimulatory (CKB1) activities. The combination of purified CKA1 and CKB1 proteins resulted in up to 100-fold stimulation of casein kinase activity compared with the CKA1 activity alone, showing that CKB1 has biochemical properties similar to those of the beta subunit from animals. CKA1 and CKB1 spontaneously assembled into a tetrameric complex, CKA1(2)CKB12, which had properties very similar to those of the oligomeric CKII form isolated from broccoli. However, the properties of the catalytic subunit CKA1 alone differed from those of the broccoli monomeric form of CKII-like activity. Phosphorylation of transcription factor GBF1 with the reconstituted CKA1(2)CKB1(2) enzyme resulted in stimulation of its DNA binding activity and retardation of the protein-DNA complex; these results are identical to those obtained previously with isolated nuclear CKII from broccoli.

Purification and characterization of casein kinase I from broccoli

Casein kinase I from broccoli was purified approximately 65,000-fold by chromatography on phosphocellulose, phenyl-Sepharose, CM-Sephacel, and affinity chromatography on N-(2-aminoethyl)-5-chloroisoquinolone-8-sulphonamide (CKI-7)-Sepharose. The catalytic subunit of casein kinase I was identified as a 36-38 kDa polypeptide doublet by using the technique of activity gel assay after SDS/PAGE with casein as a gel-incorporated substrate. A silver-stained polypeptide doublet of the same molecular mass constituted at least 95% of the protein in the final preparation, corresponding to a specific activity of approximately 1800 nmol/min per mg of protein. The enzyme was found to be a monomer by gel filtration and glycerol gradient sedimentation; the native molecular mass was calculated to be 34.2 kDa. These characteristics, as well as other essential features of plant casein kinase I activity, such as substrate specificity and sensitivity to inhibitors, were found to be similar to those established for animal casein kinase I. Broccoli casein kinase I showed weak immunological cross-reactivity with antibodies raised against bovine casein kinase I.

DNA binding activity of the Arabidopsis G-box binding factor GBF1 is stimulated by phosphorylation by casein kinase II from broccoli

To study the phosphorylation of one of the G-box binding factors from Arabidopsis (GBF1), we have obtained large amounts of this protein by expression in Escherichia coli. Bacterial GBF1 was shown to be phosphorylated very efficiently by nuclear extracts from broccoli. The phosphorylation activity was partially purified by chromatography on heparin-Sepharose and DEAE-cellulose and was characterized. It showed the essential features of casein kinase II activity: utilization of GTP in addition to ATP as a phosphate donor, strong inhibition by heparin, preference for acidic protein substrates, salt-induced binding to phosphocellulose, and salt-dependent deaggregation. The very low Km value for GBF1 (220 nM compared to approximately 10 microM for casein) was in the range observed for identified physiological substrates of casein kinase II. Phosphorylation of GBF1 resulted in stimulation of the G-box binding activity and formation of a slower migrating protein-DNA complex. The conditions of this stimulatory reaction fully corresponded to the properties of casein kinase II, in particular its dependence on the known phosphate donors. The DNA binding activity of the endogenous plant GBF was shown to be reduced by treatment with calf alkaline phosphatase; this reduction was diminished by addition of fluoride and phosphate or incubation in the presence of casein kinase II and ATP.

Biochemical Similarities between Soluble and Membrane-Bound Calcium-Dependent Protein Kinases of Barley

The soluble and membrane-bound forms of the calcium-dependent protein kinase from barley leaves (Hordeum vulgare L. cv. Borsoy) have been partially purified and compared. Both forms showed an active polypeptide of 37 kilodaltons on activity gels with incorporated histone as substrate. They eluted from chromatofocusing columns at an identical isoelectric point of pH 4.25 +/- 0.2, and also comigrated on several other chromatographic affinity media including Matrex Gel Blue A, histone-agarose, phenyl-Sepharose, and heparin-agarose. Both activities comigrated with chicken ovalbumin during gel filtration through Sephacryl S-200, indicating a native molecular mass of 45 kilodaltons. The activities share a number of enzymatic properties including salt and pH dependence, free calcium stimulation profile, substrate specificity, and Km values. The soluble activity was shown to bind to artificial lipid vesicles. These data suggest strongly that the soluble and membrane-bound calcium-dependent protein kinases from barley are very closely related or even identical.

Purification and characterization of DNA polymerase from the archaebacterium Sulfolobus acidocaldarius

DNA polymerase has been purified about 25,000-fold from the thermoacidophilic archaebacterium Sulfolobus acidocaldarius. On SDS-PAGE the enzyme was observed to have a molecular weight of 100 kDa and to be about 90% pure. The native molecular weight was 108 kDa indicating that the enzyme is composed of a single polypeptide. Activity gel analysis showed an active polypeptide of about 100 kDa. Under conditions promoting proteolysis this polypeptide was degraded to a slightly smaller form of 98 kDa. The enzyme has been characterized in respect to optimal assay conditions, template specificity, sensitivity to inhibitors and associated nuclease activities. The high temperature optimum of 65 degrees C should be emphasized. No substantial similarities have been found with other prokaryotic and eukaryotic DNA polymerases, although the enzyme bears certain resemblances to prokaryotic non-replicative polymerases.