Association of germline variants in the APOBEC3 region with cancer risk and enrichment with APOBEC-signature mutations in tumors

Identifying somatic fingerprints of cancers defined by germline and environmental risk factors

Numerous studies over the past generation have identified germline variants that increase specific cancer risks. Simultaneously, a revolution in sequencing technology has permitted high-throughput annotations of somatic genomes characterizing individual tumors. However, examining the relationship between germline variants and somatic alteration patterns is hugely challenged by the large numbers of variants in a typical tumor, the rarity of most individual variants, and the heterogeneity of tumor somatic fingerprints. In this article, we propose statistical methodology that frames the investigation of germline-somatic relationships in an interpretable manner. The method uses meta-features embodying biological contexts of individual somatic alterations to implicitly group rare mutations. Our team has used this technique previously through a multilevel regression model to diagnose with high accuracy tumor site of origin. Herein, we further leverage topic models from computational linguistics to achieve interpretable lower-dimensional embeddings of the meta-features. We demonstrate how the method can identify distinctive somatic profiles linked to specific germline variants or environmental risk factors. We illustrate the method using The Cancer Genome Atlas whole-exome sequencing data to characterize somatic tumor fingerprints in breast cancer patients with germline BRCA1/2 mutations and in head and neck cancer patients exposed to human papillomavirus.

APOBEC3A drives metastasis of high-grade serous ovarian cancer by altering epithelial-to-mesenchymal transition

High-grade serous ovarian cancer (HGSOC) is the most prevalent and aggressive histological subtype of ovarian cancer, and often presents with metastatic disease. The drivers of metastasis in HGSOC remain enigmatic. APOBEC3A (A3A), an enzyme that generates mutations across various cancers, has been proposed as a mediator of tumor heterogeneity and disease progression. However, the role of A3A in HGSOC has not been explored. Through analysis of genome sequencing from primary HGSOC, we observed an association between high levels of APOBEC3 mutagenesis and poor overall survival. We experimentally addressed this correlation by modeling A3A activity in HGSOC cell lines and mouse models which resulted in increased metastatic behavior of HGSOC cells in culture and distant metastatic spread in vivo . A3A activity in both primary and cultured HGSOC cells yielded consistent alterations in expression of epithelial-mesenchymal-transition (EMT) genes resulting in hybrid EMT and mesenchymal signatures, and providing a mechanism for their increased metastatic potential. Our findings define the prevalence of A3A mutagenesis in HGSOC and implicate A3A as a driver of HGSOC metastasis via EMT, underscoring its clinical relevance as a potential prognostic biomarker. Our study lays the groundwork for the development of targeted therapies aimed at mitigating the deleterious impact of A3A-driven EMT in HGSOC.

Regulation, functional impact, and therapeutic targeting of APOBEC3A in cancer

Enzymes of the apolipoprotein B mRNA editing catalytic polypeptide like (APOBEC) family are cytosine deaminases that convert cytosine to uracil in DNA and RNA. Among these proteins, APOBEC3 sub-family members, APOBEC3A (A3A) and APOBEC3B (A3B), are prominent sources of mutagenesis in cancer cells. The aberrant expression of A3A and A3B in cancer cells leads to accumulation of mutations with specific single-base substitution (SBS) signatures, characterized by C→T and C→G changes, in a number of tumor types. In addition to fueling mutagenesis, A3A and A3B, particularly A3A, induce DNA replication stress, DNA damage, and chromosomal instability through their catalytic activities, triggering a range of cellular responses. Thus, A3A/B have emerged as key drivers of genome evolution during cancer development, contributing to tumorigenesis, tumor heterogeneity, and therapeutic resistance. Yet, the expression of A3A/B in cancer cells presents a cancer vulnerability that can be exploited therapeutically. In this review, we discuss the recent studies that shed light on the mechanisms regulating A3A expression and the impact of A3A in cancer. We also review recent advances in the development of A3A inhibitors and provide perspectives on the future directions of A3A research.

The role of distinct APOBEC/ADAR mRNA levels in mutational signatures linked to aging and ultraviolet radiation

The APOBEC/AID family is known for its mutator activity, and recent evidence also supports the potential impact of ADARs. Furthermore, the mutator impacts of APOBEC/ADAR mutations have not yet been investigated. Assessment of pancancer TCGA exomes identified enriched somatic variants among exomes with nonsynonymous APOBEC1, APOBEC3B, APOBEC3C, ADAR, and ADARB1 mutations, compared to exomes with synonymous ones. Principal component (PC) analysis reduced the number of potential players to eight in cancer exomes/genomes, and to five in cancer types. Multivariate regression analysis was used to assess the impact of the PCs on each COSMIC mutational signature among pancancer exomes/genomes and particular cancers, identifying several novel links, including SBS17b, SBS18, and ID7 mainly determined by APOBEC1 mRNA levels; SBS40, ID1, and ID2 by age; SBS3 and SBS16 by APOBEC3A/APOBEC3B mRNA levels; ID5 and DBS9 by DNA repair/replication (DRR) defects; and SBS7a-d, SBS38, ID4, ID8, ID13, and DBS1 by ultraviolet (UV) radiation/ADARB1 mRNA levels. APOBEC/ADAR mutations appeared to potentiate the impact of DRR defects on several mutational signatures, and some factors seemed to inversely affect certain signatures. These findings potentially implicate certain APOBEC/ADAR mutations/mRNA levels in distinct mutational signatures, particularly APOBEC1 mRNA levels in aging-related signatures and ADARB1 mRNA levels in UV radiation-related signatures.

The cytidine deaminase APOBEC3A regulates nucleolar function to promote cell growth and ribosome biogenesis

Cancer initiates as a consequence of genomic mutations and its subsequent progression relies in part on increased production of ribosomes to maintain high levels of protein synthesis for unchecked cell growth. Recently, cytidine deaminases have been uncovered as sources of mutagenesis in cancer. In an attempt to form a connection between these 2 cancer driving processes, we interrogated the cytidine deaminase family of proteins for potential roles in human ribosome biogenesis. We identified and validated APOBEC3A and APOBEC4 as novel ribosome biogenesis factors through our laboratory's established screening platform for the discovery of regulators of nucleolar function in MCF10A cells. Through siRNA depletion experiments, we highlight APOBEC3A's requirement in making ribosomes and specific role within the processing and maturation steps that form the large subunit 5.8S and 28S ribosomal (r)RNAs. We demonstrate that a subset of APOBEC3A resides within the nucleolus and associates with critical ribosome biogenesis factors. Mechanistic insight was revealed by transient overexpression of both wild-type and a catalytically dead mutated APOBEC3A, which both increase cell growth and protein synthesis. Through an innovative nuclear RNA sequencing methodology, we identify only modest predicted APOBEC3A C-to-U target sites on the pre-rRNA and pre-mRNAs. Our work reveals a potential direct role for APOBEC3A in ribosome biogenesis likely independent of its editing function. More broadly, we found an additional function of APOBEC3A in cancer pathology through its function in ribosome biogenesis, expanding its relevance as a target for cancer therapeutics.

Long-read sequencing and structural variant characterization in 1,019 samples from the 1000 Genomes Project

Structural variants (SVs) contribute significantly to human genetic diversity and disease 1-4 . Previously, SVs have remained incompletely resolved by population genomics, with short-read sequencing facing limitations in capturing the whole spectrum of SVs at nucleotide resolution 5-7 . Here we leveraged nanopore sequencing 8 to construct an intermediate coverage resource of 1,019 long-read genomes sampled within 26 human populations from the 1000 Genomes Project. By integrating linear and graph-based approaches for SV analysis via pangenome graph-augmentation, we uncover 167,291 sequence-resolved SVs in these samples, considerably advancing SV characterization compared to population-wide short-read sequencing studies 3,4 . Our analysis details diverse SV classes-deletions, duplications, insertions, and inversions-at population-scale. LINE-1 and SVA retrotransposition activities frequently mediate transductions 9,10 of unique sequences, with both mobile element classes transducing sequences at either the 3'- or 5'-end, depending on the source element locus. Furthermore, analyses of SV breakpoint junctions suggest a continuum of homology-mediated rearrangement processes are integral to SV formation, and highlight evidence for SV recurrence involving repeat sequences. Our open-access dataset underscores the transformative impact of long-read sequencing in advancing the characterisation of polymorphic genomic architectures, and provides a resource for guiding variant prioritisation in future long-read sequencing-based disease studies.

Development and validation of a novel prognostic lncRNA signature based on the APOBEC3 family genes in gastric cancer

Introduction

Gastric Cancer (GC) refers to a prevalent malignant cancer accompanied by a weak prognosis. The APOBEC3 family genes and lncRNAs are linked with cancer progression. Nevertheless, there is still a scarcity of data concerning the prognostic value of APOBEC3-related lncRNAs in GC.

Methods

We extracted the data from GC samples, including transcriptome as well as clinical data, obtained from the TCGA database. Then, we screened for lncRNAs that were correlated with the APOBEC3 family genes and constructed an APOBEC3-related lncRNA prognostic signature (LPS) by utilizing univariate Cox and lasso regression analysis. Furthermore, we validated our constructed signature and evaluated it thoroughly, including analysis of its function, immunity, mutations, and clinical applications via multiple methods, including Metascape, GSEA, and analyses including TIC and TME, immune checkpoints, CNV and SNPs, Kaplan-Meier survival curves, nomogram, decision tree and drug prediction analysis. Finally, we overexpressed LINC01094 to evaluate the impacts on the proliferation as well as migration with regards to KATO-2 cells.

Results

We selected eight lncRNAs for our APOBEC3-related LPS, which is demonstrated as a valuable tool in predicting the individual GC patients' prognosis. Subsequently, we segregated the samples into subgroups of high-as well as low-risk relying on the risk score with regards to APOBEC3-related LPS. By performing functional analysis, we have shown that immune-as well as tumor-related pathways were enriched in high- and low-risk GC patients. Furthermore, immune analysis revealed a robust correlation between the APOBEC3-related LPS and immunity. We found that immune checkpoints were significantly associated with the APOBEC3-related LPS and were greatly exhibited in GC tumor and high-risk samples. Mutational analysis suggested that the mutational rate was greater in low-risk samples. Furthermore, we predicted small molecular drugs displayed greater sensitivity in patients categorized as high-risk. Moreover, the immune response was also better in high-risk patients. Of these drugs, dasatinib was significant in both methods and might be considered a potential novel drug for treating high-risk GC patients. Finally, we found that LINC01094 has the potential to enhance the migration, proliferation as well as inhibit apoptosis of KATO-2 in GC cells. And Dasatinib has an inhibitory effect on the migration as well as proliferation in GC cells.

Conclusion

We created a novel APOBEC3-related LPS in predicting the prognosis with regards to individual GC patients. Importantly, this APOBEC3-related LPS was closely associated with immunity and might guide clinical treatment.

Progress in understanding of relationship between inflammation and tumors

Over the past decade, there has been clear evidence that inflammation plays a key role in tumorigenesis. Tumor extrinsic inflammation is caused by many factors, including bacterial and viral infections, autoimmune diseases, obesity, smoking, excessive alcohol consumption, etc., all of which can increase cancer risk and stimulate malignant progression. Conversely, inflammation inherent in cancer or caused by cancer can be triggered by cancer-initiating mutations and can promote malignant progression through recruitment and activation of inflammatory cells. Both exogenous and endogenous inflammation can lead to immunosuppression, thus providing a preferred opportunity for tumor development. Studies have confirmed that chronic inflammation is involved in various steps of tumorigenesis, including cell transformation, promotion, survival, prolifer-ation, invasion, angiogenesis, and metastasis. Recent research has shed new light on the molecular and cellular circuits between inflammation and cancer. Two pathways have been preliminarily identified: Intrinsic and extrinsic. In the intrinsic pathway, genetic events leading to tumors initiate the expression of inflammatory related programs and guide the construction of the inflammatory microenvironment. In the extrinsic pathway, inflammatory conditions promote the development of cancer. This article reviews the recent progress in the understanding of the relationship between inflammation and tumors.

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.

APOBEC Mutagenesis in Cancer Development and Susceptibility

APOBEC cytosine deaminases are prominent mutators in cancer, mediating mutations in over 50% of cancers. APOBEC mutagenesis has been linked to tumor heterogeneity, persistent cell evolution, and therapy responses. While emerging evidence supports the impact of APOBEC mutagenesis on cancer progression, the understanding of its contribution to cancer susceptibility and malignant transformation is limited. We examine the existing evidence for the role of APOBEC mutagenesis in carcinogenesis on the basis of the reported associations between germline polymorphisms in genes encoding APOBEC enzymes and cancer risk, insights into APOBEC activities from sequencing efforts of both malignant and non-malignant human tissues, and in vivo studies. We discuss key knowledge gaps and highlight possible ways to gain a deeper understanding of the contribution of APOBEC mutagenesis to cancer development.

Data-adaptive and pathway-based tests for association studies between somatic mutations and germline variations in human cancers

Cancer is a disease driven by a combination of inherited genetic variants and somatic mutations. Recently available large-scale sequencing data of cancer genomes have provided an unprecedented opportunity to study the interactions between them. However, previous studies on this topic have been limited by simple, low statistical power tests such as Fisher's exact test. In this paper, we design data-adaptive and pathway-based tests based on the score statistic for association studies between somatic mutations and germline variations. Previous research has shown that two single-nucleotide polymorphism (SNP)-set-based association tests, adaptive sum of powered score (aSPU) and data-adaptive pathway-based (aSPUpath) tests, increase the power in genome-wide association studies (GWASs) with a single disease trait in a case-control study. We extend aSPU and aSPUpath to multi-traits, that is, somatic mutations of multiple genes in a cohort study, allowing extensive information aggregation at both SNP and gene levels.p $p$ -values from different parameters assuming varying genetic architecture are combined to yield data-adaptive tests for somatic mutations and germline variations. Extensive simulations show that, in comparison with some commonly used methods, our data-adaptive somatic mutations/germline variations tests can be applied to multiple germline SNPs/genes/pathways, and generally have much higher statistical powers while maintaining the appropriate type I error. The proposed tests are applied to a large-scale real-world International Cancer Genome Consortium whole genome sequencing data set of 2583 subjects, detecting more significant and biologically relevant associations compared with the other existing methods on both gene and pathway levels. Our study has systematically identified the associations between various germline variations and somatic mutations across different cancer types, which potentially provides valuable utility for cancer risk prediction, prognosis, and therapeutics.

Mutational impact of APOBEC3A and APOBEC3B in a human cell line and comparisons to breast cancer

A prominent source of mutation in cancer is single-stranded DNA cytosine deamination by cellular APOBEC3 enzymes, which results in signature C-to-T and C-to-G mutations in TCA and TCT motifs. Although multiple enzymes have been implicated, reports conflict and it is unclear which protein(s) are responsible. Here we report the development of a selectable system to quantify genome mutation and demonstrate its utility by comparing the mutagenic activities of three leading candidates-APOBEC3A, APOBEC3B, and APOBEC3H. The human cell line, HAP1, is engineered to express the thymidine kinase (TK) gene of HSV-1, which confers sensitivity to ganciclovir. Expression of APOBEC3A and APOBEC3B, but not catalytic mutant controls or APOBEC3H, triggers increased frequencies of TK mutation and similar TC-biased cytosine mutation profiles in the selectable TK reporter gene. Whole genome sequences from independent clones enabled an analysis of thousands of single base substitution mutations and extraction of local sequence preferences with APOBEC3A preferring YTCW motifs 70% of the time and APOBEC3B 50% of the time (Y = C/T; W = A/T). Signature comparisons with breast tumor whole genome sequences indicate that most malignancies manifest intermediate percentages of APOBEC3 signature mutations in YTCW motifs, mostly between 50 and 70%, suggesting that both enzymes contribute in a combinatorial manner to the overall mutation landscape. Although the vast majority of APOBEC3A- and APOBEC3B-induced single base substitution mutations occur outside of predicted chromosomal DNA hairpin structures, whole genome sequence analyses and supporting biochemical studies also indicate that both enzymes are capable of deaminating the single-stranded loop regions of DNA hairpins at elevated rates. These studies combine to help resolve a long-standing etiologic debate on the source of APOBEC3 signature mutations in cancer and indicate that future diagnostic and therapeutic efforts should focus on both APOBEC3A and APOBEC3B.

Analysis of Several Common APOBEC-type Mutations in Bladder Tumors Suggests Links to Viral Infection

FGFR3 and PIK3CA are among the most frequently mutated genes in bladder tumors. We hypothesized that recurrent mutations in these genes might be caused by common carcinogenic exposures such as smoking and other factors. We analyzed 2,816 bladder tumors with available data on FGFR3 and/or PIK3CA mutations, focusing on the most recurrent mutations detected in ≥10% of tumors. Compared to tumors with other FGFR3/PIK3CA mutations, FGFR3-Y375C was more common in tumors from smokers than never-smokers (P = 0.009), while several APOBEC-type driver mutations were enriched in never-smokers: FGFR3-S249C (P = 0.013) and PIK3CA-E542K/PIK3CA-E545K (P = 0.009). To explore possible causes of these APOBEC-type mutations, we analyzed RNA sequencing (RNA-seq) data from 798 bladder tumors and detected several viruses, with BK polyomavirus (BKPyV) being the most common. We then performed IHC staining for polyomavirus (PyV) Large T-antigen (LTAg) in an independent set of 211 bladder tumors. Overall, by RNA-seq or IHC-LTAg, we detected PyV in 26 out of 1,010 bladder tumors with significantly higher detection (P = 4.4 × 10-5), 25 of 554 (4.5%) in non-muscle-invasive bladder cancers (NMIBC) versus 1 of 456 (0.2%) of muscle-invasive bladder cancers (MIBC). In the NMIBC subset, the FGFR3/PIK3CA APOBEC-type driver mutations were detected in 94.7% (18/19) of PyV-positive versus 68.3% (259/379) of PyV-negative tumors (P = 0.011). BKPyV tumor positivity in the NMIBC subset with FGFR3- or PIK3CA-mutated tumors was also associated with a higher risk of progression to MIBC (P = 0.019). In conclusion, our results support smoking and BKPyV infection as risk factors contributing to bladder tumorigenesis in the general patient population through distinct molecular mechanisms.

PREVENTION RELEVANCE

Tobacco smoking likely causes one of the most common mutations in bladder tumors (FGFR3-Y375C), while viral infections might contribute to three others (FGFR3-S249C, PIK3CA-E542K, and PIK3CA-E545K). Understanding the causes of these mutations may lead to new prevention and treatment strategies, such as viral screening and vaccination.

Historical Pathogen-Driven Selection May Contribute to Contemporary Ethnic Difference in Bladder Cancer Susceptibility

BACKGROUND

The rationale for ethnic differences in bladder cancer (BCa) susceptibility is an important open question. In this study, we raised the hypothesis that the APOBEC3-rs1014971 variant associated with BCa risk and APOBEC-mutagenesis probably contribute to ethnic differences.

METHODS

We calculated the ethnicity-stratified 5-year age-adjusted incidence rates of BCa using the US SEER database. We performed somatic mutational-signature analyses and compared the APOBEC-related mutational contribution across BCa tumors in patients of different ethnicities. We analyzed the allele frequency distribution of APOBEC3-related rs1014971 in contemporary populations of different ethnicities and in ancient human genomes. We also analyzed the natural selection profiles and ages of the investigated SNPs.

RESULTS

We validated the ethnic difference in BCa risk using US SEER data, revealing Caucasians to be at >2-fold greater risk than Asians / Pacific islanders. In contemporary populations, we observed a coherent ethnic distribution in terms not only of the allele frequency of APOBEC3-related rs1014971, but also the mutational contribution of APOBEC-mediated mutagenesis in BCa tumors. Population genetics and ancient genome analyses further suggested that the diverse ethnic distribution of rs1014971 could be rooted in human evolution.

CONCLUSIONS

It is possible that APOBEC3-related rs1014971 is involved in the different BCa incidence across ethnic groups, and this difference is potentially derived from human evolution. Our findings suggested an evolutionary link between contemporary population-level variations in malignancy susceptibility and pathogen-driven selection in the past, not unlike previously reported cases of certain autoimmune and metabolic disorders.

Proteogenomics of different urothelial bladder cancer stages reveals distinct molecular features for papillary cancer and carcinoma in situ

The progression of urothelial bladder cancer (UC) is a complicated multi-step process. We perform a comprehensive multi-omics analysis of 448 samples from 190 UC patients, covering the whole spectrum of disease stages and grades. Proteogenomic integration analysis indicates the mutations of HRAS regulated mTOR signaling to form urothelial papilloma rather than papillary urothelial cancer (PUC). DNA damage is a key signaling pathway in the progression of carcinoma in situ (CIS) and related to APOBEC signature. Glucolipid metabolism increase and lower immune cell infiltration are associated with PUC compared to CIS. Proteomic analysis distinguishes the origins of invasive tumors (PUC-derived and CIS-derived), related to distinct clinical prognosis and molecular features. Additionally, loss of RBPMS, associated with CIS-derived tumors, is validated to increase the activity of AP-1 and promote metastasis. This study reveals the characteristics of two distinct branches (PUC and CIS) of UC progression and may eventually benefit clinical practice.

Genetic inhibitors of APOBEC3B-induced mutagenesis

The cytidine deaminases APOBEC3A (A3A) and APOBEC3B (A3B) are prominent mutators of human cancer genomes. However, tumor-specific genetic modulators of APOBEC-induced mutagenesis are poorly defined. Here, we used a screen to identify 61 gene deletions that increase A3B-induced mutations in yeast. We also determined whether each deletion was epistatic with Ung1 loss, which indicated whether the encoded factors participate in the homologous recombination (HR)-dependent bypass of A3B/Ung1-dependent abasic sites or suppress A3B-catalyzed deamination by protecting against aberrant formation of single-stranded DNA (ssDNA). We found that the mutation spectra of A3B-induced mutations revealed genotype-specific patterns of strand-specific ssDNA formation and nucleotide incorporation across APOBEC-induced lesions. Combining these three metrics, we were able to establish a multifactorial signature of APOBEC-induced mutations specific to (1) failure to remove H3K56 acetylation, (2) defective CTF18-RFC complex function, and (3) defective HR-mediated bypass of APOBEC-induced lesions. We extended these results by analyzing mutation data for human tumors and found BRCA1/2-deficient breast cancers display three- to fourfold more APOBEC-induced mutations. Mirroring our results in yeast, Rev1-mediated C-to-G substitutions are mainly responsible for increased APOBEC-signature mutations in BRCA1/2-deficient tumors, and these mutations associate with lagging strand synthesis during replication. These results identify important factors that influence DNA replication dynamics and likely the abundance of APOBEC-induced mutation during tumor progression. They also highlight a novel role for BRCA1/2 during HR-dependent lesion bypass of APOBEC-induced lesions during cancer cell replication.

Profiling the Somatic Mutational Landscape of Breast Tumors from Hispanic/Latina Women Reveals Conserved and Unique Characteristics

Somatic mutational profiling is increasingly being used to identify potential targets for breast cancer. However, limited tumor-sequencing data from Hispanic/Latinas (H/L) are available to guide treatment. To address this gap, we performed whole-exome sequencing (WES) and RNA sequencing on 146 tumors and WES of matched germline DNA from 140 H/L women in California. Tumor intrinsic subtype, somatic mutations, copy-number alterations, and expression profiles of the tumors were characterized and compared with data from tumors of non-Hispanic White (White) women in The Cancer Genome Atlas (TCGA). Eight genes were significantly mutated in the H/L tumors including PIK3CA, TP53, GATA3, MAP3K1, CDH1, CBFB, PTEN, and RUNX1; the prevalence of mutations in these genes was similar to that observed in White women in TCGA. Four previously reported Catalogue of Somatic Mutations in Cancer (COSMIC) mutation signatures (1, 2, 3, 13) were found in the H/L dataset, along with signature 16 that has not been previously reported in other breast cancer datasets. Recurrent amplifications were observed in breast cancer drivers including MYC, FGFR1, CCND1, and ERBB2, as well as a recurrent amplification in 17q11.2 associated with high KIAA0100 gene expression that has been implicated in breast cancer aggressiveness. In conclusion, this study identified a higher prevalence of COSMIC signature 16 and a recurrent copy-number amplification affecting expression of KIAA0100 in breast tumors from H/L compared with White women. These results highlight the necessity of studying underrepresented populations.

SIGNIFICANCE

Comprehensive characterization of genomic and transcriptomic alterations in breast tumors from Hispanic/Latina patients reveals distinct genetic alterations and signatures, demonstrating the importance of inclusive studies to ensure equitable care for patients. See related commentary by Schmit et al., p. 2443.

Analysis of somatic mutations and key driving factors of cervical cancer progression

We investigated the somatic mutations and key driving factors of cervical cancer by whole exome sequencing . We found 22,183 somatic single nucleotide variations (SNVs) in 52 paired samples. Somatic SNVs in cervical cancer were significantly higher than those in high-grade intraepithelial lesion and low-grade squamous intraepithelial lesion groups (P < 0.05). C → T/G accounted for 44.12% of base substitution. Copy number variation (false discovery rate < 0.05) was found in 57 chromosome regions. The three regions with significant differences between cervical cancer and non-cervical cancer groups were 1q21.1, 3q26.33, and 13q33.1, covering genes related to tumor proliferation, differentiation, and apoptosis. The frequency of human papillomavirus (HPV) insertion/integration and the number of "tCw" mutations in the cervical cancer group were significantly higher than those in the non-cervical cancer group (P < 0.05). The total number of mutations was positively correlated with the number of "tCw" mutations (R 2 = 0.7967). HPV insertion/integration (OR = 2.302, CI = 1.523-3.589, P = 0.0005), APOBEC enrichment (OR = 17.875, CI = 2.117-150.937, P = 0.001), and HLA-B*39 in HLA-I (OR = 6.435, CI = 0.823-48.919, P = 0.0042) were risk factors for cervical cancer, while HLA-DQB1*05 in HLA-II was a protective factor (OR = 0.426, CI = 0.197-0.910, P = 0.032). Conclusively, HPV insertion/integration, APOBEC mutagenesis, and HLA polymorphisms are high-risk factors for cervical cancer and may be causes of genome instability and somatic mutations. This study provides experimental data for revealing the molecular mechanism of cervical cancer.

Mutational Signatures in Gastric Cancer and Their Clinical Implications

Gastric cancer is characterised by high inter- and intratumour heterogeneity. The majority of patients are older than 65 years and the global burden of this disease is increasing due to the aging of the population. The disease is usually diagnosed at advanced stages, which is a consequence of nonspecific symptoms. Few improvements have been made at the level of noninvasive molecular diagnosis of sporadic gastric cancer, and therefore the mortality rate remains high. A new field of mutational signatures has emerged in the past decade with advances in the genome sequencing technology. These distinct mutational patterns in the genome, caused by exogenous and endogenous mutational processes, can be associated with tumour aetiology and disease progression, and could provide novel perception on the treatment possibilities. This review assesses the mutational signatures found in gastric cancer and summarises their potential for use in clinical setting as diagnostic or prognostic biomarkers. Associated treatment options and biomarkers already implemented in clinical use are discussed, together with those that are still being explored or are in clinical studies.

Genome-wide Association Study of Bladder Cancer Reveals New Biological and Translational Insights

BACKGROUND

Genomic regions identified by genome-wide association studies (GWAS) for bladder cancer risk provide new insights into etiology.

OBJECTIVE

To identify new susceptibility variants for bladder cancer in a meta-analysis of new and existing genome-wide genotype data.

DESIGN, SETTING, AND PARTICIPANTS

Data from 32 studies that includes 13,790 bladder cancer cases and 343,502 controls of European ancestry were used for meta-analysis.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSES

Log-additive associations of genetic variants were assessed using logistic regression models. A fixed-effects model was used for meta-analysis of the results. Stratified analyses were conducted to evaluate effect modification by sex and smoking status. A polygenic risk score (PRS) was generated on the basis of known and novel susceptibility variants and tested for interaction with smoking.

RESULTS AND LIMITATIONS

Multiple novel bladder cancer susceptibility loci (6p.22.3, 7q36.3, 8q21.13, 9p21.3, 10q22.1, 19q13.33) as well as improved signals in three known regions (4p16.3, 5p15.33, 11p15.5) were identified, bringing the number of independent markers at genome-wide significance (p < 5 × 10-8) to 24. The 4p16.3 (FGFR3/TACC3) locus was associated with a stronger risk for women than for men (p-interaction = 0.002). Bladder cancer risk was increased by interactions between smoking status and genetic variants at 8p22 (NAT2; multiplicative p value for interaction [pM-I] = 0.004), 8q21.13 (PAG1; pM-I = 0.01), and 9p21.3 (LOC107987026/MTAP/CDKN2A; pM-I = 0.02). The PRS based on the 24 independent GWAS markers (odds ratio per standard deviation increase 1.49, 95% confidence interval 1.44-1.53), which also showed comparable results in two prospective cohorts (UK Biobank, PLCO trial), revealed an approximately fourfold difference in the lifetime risk of bladder cancer according to the PRS (e.g., 1st vs 10th decile) for both smokers and nonsmokers.

CONCLUSIONS

We report novel loci associated with risk of bladder cancer that provide clues to its biological underpinnings. Using 24 independent markers, we constructed a PRS to stratify lifetime risk. The PRS combined with smoking history, and other established risk factors, has the potential to inform future screening efforts for bladder cancer.

PATIENT SUMMARY

We identified new genetic markers that provide biological insights into the genetic causes of bladder cancer. These genetic risk factors combined with lifestyle risk factors, such as smoking, may inform future preventive and screening strategies for bladder cancer.

Similar deamination activities but different phenotypic outcomes induced by APOBEC3 enzymes in breast epithelial cells

APOBEC3 (A3) enzymes deaminate cytosine to uracil in viral single-stranded DNA as a mutagenic barrier for some viruses. A3-induced deaminations can also occur in human genomes resulting in an endogenous source of somatic mutations in multiple cancers. However, the roles of each A3 are unclear since few studies have assessed these enzymes in parallel. Thus, we developed stable cell lines expressing A3A, A3B, or A3H Hap I using non-tumorigenic MCF10A and tumorigenic MCF7 breast epithelial cells to assess their mutagenic potential and cancer phenotypes in breast cells. The activity of these enzymes was characterized by γH2AX foci formation and in vitro deamination. Cell migration and soft agar colony formation assays assessed cellular transformation potential. We found that all three A3 enzymes had similar γH2AX foci formation, despite different deamination activities in vitro. Notably, in nuclear lysates, the in vitro deaminase activity of A3A, A3B, and A3H did not require digestion of cellular RNA, in contrast to that of A3B and A3H in whole-cell lysates. Their similar activities in cells, nonetheless, resulted in distinct phenotypes where A3A decreased colony formation in soft agar, A3B decreased colony formation in soft agar after hydroxyurea treatment, and A3H Hap I promoted cell migration. Overall, we show that in vitro deamination data do not always reflect cell DNA damage, all three A3s induce DNA damage, and the impact of each is different.

Clinical Implications of APOBEC3-Mediated Mutagenesis in Breast Cancer

Over recent years, members of the APOBEC3 family of cytosine deaminases have been implicated in increased cancer genome mutagenesis, thereby contributing to intratumor and intertumor genomic heterogeneity and therapy resistance in, among others, breast cancer. Understanding the available methods for clinical detection of these enzymes, the conditions required for their (dysregulated) expression, the clinical impact they have, and the clinical implications they may offer is crucial in understanding the current impact of APOBEC3-mediated mutagenesis in breast cancer. Here, we provide a comprehensive review of recent developments in the detection of APOBEC3-mediated mutagenesis and responsible APOBEC3 enzymes, summarize the pathways that control their expression, and explore the clinical ramifications and opportunities they pose. We propose that APOBEC3-mediated mutagenesis can function as a helpful predictive biomarker in several standard-of-care breast cancer treatment plans and may be a novel target for treatment.

APOBEC3A/B Polymorphism Is Not Associated with Human Papillomavirus Infection and Cervical Carcinogenesis

The persistence of a high-risk Human papillomavirus (HPV-HR) infection of the cervix results in different manifestations of lesions depending on the immunologic capacity of the host. Variations in apolipoprotein B mRNA editing enzyme catalytic polypeptide (APOBEC)-like genes, such as the APOBEC3A/B deletion hybrid polymorphism (A3A/B), may contribute to cervical malignancy in the presence of HPV. The aim of this study was to investigate the association between the A3A/B polymorphism and HPV infection and the development of cervical intraepithelial lesions and cervical cancer in Brazilian women. The study enrolled 369 women, who were categorized according to the presence of infection and subdivided according to the degree of intraepithelial lesion and cervical cancer. APOBEC3A/B was genotyped by allele-specific polymerase chain reaction (PCR). As for the A3A/B polymorphism, the distribution of genotypes was similar between groups and among the analyzed subgroups. There were no significant differences in the presence of infection or development of lesions, even after exclusion of confounding factors. This is the first study to show that the A3A/B polymorphism is not associated with HPV infection and the development of intraepithelial lesions and cervical cancer in Brazilian women.

Germline Cancer Gene Expression Quantitative Trait Loci Are Associated with Local and Global Tumor Mutations

Somatic mutations drive cancer development and are relevant to patient responses to treatment. Emerging evidence shows that variations in the somatic genome can be influenced by the germline genetic background. However, the mechanisms underlying these germline-somatic associations remain largely obscure. We hypothesized that germline variants can influence somatic mutations in a nearby cancer gene ("local impact") or a set of recurrently mutated cancer genes across the genome ("global impact") through their regulatory effect on gene expression. To test this hypothesis, tumor targeted sequencing data from 12,413 patients across 11 cancer types in the Dana-Farber Profile cohort were integrated with germline cancer gene expression quantitative trait loci (eQTL) from the Genotype-Tissue Expression Project. Variants that upregulate ATM expression were associated with a decreased risk of somatic ATM mutations across 8 cancer types. GLI2, WRN, and CBFB eQTL were associated with global tumor mutational burden of cancer genes in ovarian cancer, glioma, and esophagogastric carcinoma, respectively. An EPHA5 eQTL was associated with mutations in cancer genes specific to colorectal cancer, and eQTL related to expression of APC, WRN, GLI1, FANCA, and TP53 were associated with mutations in genes specific to endometrial cancer. These findings provide evidence that germline-somatic associations are mediated through expression of specific cancer genes, opening new avenues for research on the underlying biological processes.

SIGNIFICANCE

Analysis of associations between the germline genetic background and somatic mutations in patients with cancer suggests that germline variants can influence local and global tumor mutations by altering expression of cancer-related genes. See related commentary by Kar, p. 1165.

Genetic modifiers of APOBEC-induced mutagenesis

The cytidine deaminases APOBEC3A and APOBEC3B (A3B) are prominent mutators of human cancer genomes. However, tumor-specific genetic modulators of APOBEC-induced mutagenesis are poorly defined. Here, we utilized a screen to identify 61 gene deletions that increase A3B-induced mutations in yeast. Also, we determined whether each deletion was epistatic with UNG1 loss, which indicated whether the encoded factors participate in the error-free bypass of A3B/Ung1-dependent abasic sites or suppress A3B-catalyzed deamination by protecting against aberrant formation of single stranded DNA (ssDNA). Additionally, we determined that the mutation spectra of A3B-induced mutations revealed genotype-specific patterns of strand-specific ssDNA formation and nucleotide incorporation across APOBEC-induced lesions. Combining these three metrics we were able to establish a multifactorial signature of APOBEC-induced mutations specific to (1) failure to remove H3K56 acetylation, which results in extremely high A3B-induced mutagenesis, (2) defective CTF18-RFC complex function, which results in high levels of A3B induced mutations specifically on the leading strand template that synergistically increase with loss of UNG1, and (3) defective HR-mediated bypass of APOBEC-induced lesions, which were epistatic with Ung1 loss and result from increased Rev1-mediated C-to-G substitutions. We extended these results by analyzing mutation data for human tumors and found BRCA1/2-deficient breast cancer tumors display 3- to 4-fold more APOBEC-induced mutations. Mirroring our results in yeast, for BRCA1/2 deficient tumors Rev1-mediated C-to-G substitutions are solely responsible for increased APOBEC-signature mutations and these mutations occur on the lagging strand during DNA replication. Together these results identify important factors that influence the dynamics of DNA replication and likely the abundance of APOBEC-induced mutation during tumor progression as well as a novel mechanistic role for BRCA1/2 during HR-dependent lesion bypass of APOBEC-induced lesions during cancer cell replication.

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

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

Involvement of APOBEC3A/B Deletion in Mouse Mammary Tumor Virus (MMTV)-like Positive Human Breast Cancer

The association between mouse mammary tumor virus (MMTV)-like sequences and human breast cancer (BC) is largely documented in the literature, but further research is needed to determine how they influence carcinogenesis. APOBEC3 cytidine deaminases are viral restriction factors that have been implicated in cancer mutagenesis, and a germline deletion that results in the fusion of the APOBEC3A coding region with the APOBEC3B 3'-UTR has been linked to increased mutagenic potential, enhanced risk of BC development, and poor prognosis. However, little is known about factors influencing APOBEC3 family activation in cancer. Thus, we hypothesized that MMTV infection and APOBEC3-mediated mutagenesis may be linked in the pathogenesis of BC. We investigated APOBEC3A/B genotyping, MMTV-like positivity, and clinicopathological parameters of 209 BC patients. We show evidence for active APOBEC3-mediated mutagenesis in human-derived MMTV sequences and comparatively investigate the impact of APOBEC3A/B germline deletion in MMTV-like env positive and negative BC in a Brazilian cohort. In MMTV-like negative samples, APOBEC3A/B deletion was negatively correlated with tumor stage while being positively correlated with estrogen receptor expression. Although APOBEC3A/B was not associated with MMTV-like positivity, samples carrying both MMTV-like positivity and APOBEC3A/B deletion had the lowest age-at-diagnosis of all study groups, with all patients being less than 50 years old. These results indicate that APOBEC3 mutagenesis is active against MMTV-like sequences, and that APOBEC3A/B deletion might act along with the MMTV-like presence to predispose people to early-onset BC.

Integrated multi-omics analysis identifies CD73 as a prognostic biomarker and immunotherapy response predictor in head and neck squamous cell carcinoma

BACKGROUND

Advances in tumor immunotherapy have been developed for patients with advanced recurrent or metastatic (R/M) HNSCC. However, the response of most HNSCC patients to immune checkpoint inhibitors (ICI) remains unsatisfactory. CD73 is a promising target for tumor immunotherapy, but its role in HNSCC remains insufficient. In this study, we aim to explore the function of CD73 in HNSCC.

METHODS

Transcriptomic and clinical data of TCGA-HNSC were downloaded from UCSC Xena for analysis of CD73 mRNA expression and prognosis. Immunohistochemical assay were performed to validate the expression of CD73 in tumor tissues and its relationship with CD8+ T cells. GSEA analysis was performed with the "clusterProfiler" R package. Immune infiltration analysis was calculated with ESTIMATE, CIBERSORT and MCP-counter algorithms. Single-cell transcriptomic data was originated from GSE103322. Cell clustering, annotation and CD73 expression were from the TISCH database. Correlation data between CD73 and tumor signatures were obtained from the CancerSEA database. Somatic mutation data were obtained from TCGA-HNSC and analyzed by "maftools" R package. Immune efficacy prediction was performed using TIDE algorithm and validated with the IMvigor210 cohort.

RESULTS

Compared with normal tissues, both mRNA and protein expressions of CD73 were elevated in tumor tissues (P = 9.7×10-10, P = 7.6×10-5, respectively). Kaplan-Meier analysis revealed that patients with high expression of CD73 had worse overall survival (log-rank P = 0.0094), and CD73 could be used as a diagnostic factor for HNSCC (AUC = 0.778). Both bulk RNA-seq and single-cell RNA-seq analysis showed that high CD73 expression can promote EMT and metastasis, samples with high CD73 expression had reduced CD8+ T cells. Furthermore, it was found that CD73-high group was more prone to have mutations in TP53, HRAS and CDKN2A, and were negatively correlated with TMB (P = 0.0055) and MSI (P = 0.00034). Mutational signature analysis found that CD73 was associated with APOBEC signature. Immunotherapy efficacy analysis showed that CD73-high group was less sensitive to immune efficacy.

CONCLUSIONS

Our results demonstrate that CD73 has an inhibitory effect on the tumor microenvironment, and is more likely to be unresponsive to ICI therapy. Collectively, targeting CD73 may provide new insights for tumor targeted therapy and/or immunotherapy.

APOBEC3A /B deletion polymorphism and endometrial cancer risk

BACKGROUND

A common 30 kb deletion affecting the APOBEC3A and APOBEC3B genes has been linked to increased APOBEC activity and APOBEC-related mutational signatures in human cancers. The role of this deletion as a cancer risk factor remains controversial.

MATERIALS AND METHODS

We genotyped the APOBEC3A/B deletion in a sample of 1,470 Norwegian endometrial cancer cases and compared to 1,918 healthy controls. For assessment across Caucasian populations, we mined genotypes of the SNP rs12628403, which is in strong linkage disequilibrium with the deletion, in a GWAS dataset of 4,274 cases and 18,125 healthy controls, through the ECAC consortium.

RESULTS

We found the APOBEC3A/B deletion variant to be significantly associated with reduced risk of endometrial cancer among Norwegian women (OR = 0.75; 95% CI = 0.62-0.91; p = 0.003; dominant model). Similar results were found in the subgroup of endometrioid endometrial cancer (OR = 0.64; 95% CI = 0.51-0.79; p = 3.6 × 10^-5 ; dominant model). The observed risk reduction was particularly strong among individuals in the range of 50-60 years of age (OR = 0.51; 95% CI = 0.33-0.78; p = 0.002; dominant model). In the different populations included in the ECAC dataset, the ORs varied from 0.85 to 1.05. Although five out of six populations revealed ORs <1.0, the overall estimate was nonsignificant and, as such, did not formally validate the findings in the Norwegian cohort.

CONCLUSION

The APOBEC3A/B deletion polymorphism is associated with a decreased risk of endometrial cancer in the Norwegian population.

Addressing the benefits of inhibiting APOBEC3-dependent mutagenesis in cancer

Mutational signatures associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC)3 cytosine deaminase activity have been found in over half of cancer types, including some therapy-resistant and metastatic tumors. Driver mutations can occur in APOBEC3-favored sequence contexts, suggesting that mutagenesis by APOBEC3 enzymes may drive cancer evolution. The APOBEC3-mediated signatures are often detected in subclonal branches of tumor phylogenies and are acquired in cancer cell lines over long periods of time, indicating that APOBEC3 mutagenesis can be ongoing in cancer. Collectively, these and other observations have led to the proposal that APOBEC3 mutagenesis represents a disease-modifying process that could be inhibited to limit tumor heterogeneity, metastasis and drug resistance. However, critical aspects of APOBEC3 biology in cancer and in healthy tissues have not been clearly defined, limiting well-grounded predictions regarding the benefits of inhibiting APOBEC3 mutagenesis in different settings in cancer. We discuss the relevant mechanistic gaps and strategies to address them to investigate whether inhibiting APOBEC3 mutagenesis may confer clinical benefits in cancer.

Prognostic value of an APOBEC3 deletion polymorphism for glioma patients in Taiwan

OBJECTIVE

The molecular pathogenesis of malignant gliomas, characterized by diverse tumor histology with differential prognosis, remains largely unelucidated. An APOBEC3 deletion polymorphism, with a deletion in APOBEC3B, has been correlated to risk and prognosis in several cancers, but its role in glioma is unclear. The authors aimed to examine the clinical relevance of the APOBEC3 deletion polymorphism to glioma risk and survival in a glioma patient cohort in Taiwan.

METHODS

The authors detected deletion genotypes in 403 glioma patients and 1365 healthy individuals in Taiwan and correlated the genotypes with glioma risk, clinicopathological factors, patient survival, and patient sex. APOBEC3 gene family expression was measured and correlated to the germline deletion. A nomogram model was constructed to predict patient survival in glioma.

RESULTS

The proportion of APOBEC3B-/- and APOBEC3B+/- genotypes was higher in glioblastoma (GBM) patients than healthy individuals and correlated with higher GBM risk in males. A higher percentage of cases with APOBEC3B- was observed in male than female glioma patients. The presence of APOBEC3B-/- was correlated with better overall survival (OS) in male astrocytic glioma patients. No significant correlation of the genotypes to glioma risk and survival was observed in the female patient cohort. Lower APOBEC3B expression was observed in astrocytic glioma patients with APOBEC3B-/- and was positively correlated with better OS. A 5-factor nomogram model was constructed based on male patients with astrocytic gliomas in the study cohort and worked efficiently for predicting patient OS.

CONCLUSIONS

The germline APOBEC3 deletion was associated with increased GBM risk and better OS in astrocytic glioma patients in the Taiwan male population. The APOBEC3B deletion homozygote was a potential independent prognostic factor predicting better survival in male astrocytic glioma patients.

Prognostic Value of HPV E6 and APOBEC3B in Upper Urinary Tract Urothelial Carcinoma

Background

APOBEC mutation signature is common in upper urinary tract urothelial carcinoma (UTUC). When virus infection occurs, upregulated APOBEC plays an antiviral role by deoxycytidine deaminase activity. However, the carcinogenic roles of HPV E6 protein and APOBEC mutation signature in UTUC have not been investigated.

Aims

This study explored the relationship among HPV E6, APOBEC, and clinicopathological characteristics in patients with UTUC and impacts of their expression on the prognosis.

Methods

The expression of HPV E6 and APOBEC3B of 78 patients with UTUC was detected by immunohistochemistry. Correlation of HPV E6 and APOBEC3B expression levels with clinicopathological characteristics was statistically analyzed. Univariate and multivariate Cox regression analyses were used to evaluate the prognosis of HPV E6 and APOBEC3B for disease-free survival (DFS); survival analysis was performed using Kaplan-Meier methods.

Results

The expression of APOBEC3B was positively correlated with the expression of HPV E6 (r = 0.383, P = 0.001). HPV E6 was significantly increased in patients with stage I (χ² = 4.938, P = 0.026) and low-grade urothelial carcinoma (χ² = 3.939, P = 0.047), as well as in patients without LVI (χ² = 4.064, P = 0.044). Meanwhile, APOBEC3B was highly expressed in patients with stage I (χ² = 4.057, P = 0.044) and low-grade urothelial carcinoma (χ² = 7.153, P = 0.007). Multivariate Cox regression analysis revealed the APOBEC3B expression was the independent prognostic factor for DFS, Kaplan-Meier survival analysis showed that low expression of APOBEC3B and HPV E6 was significantly associated with the poor prognosis of UTUC patients.

Conclusion

HPV E6 expression is positively associated with APOBEC3B expression, and the high expression of HPV E6 and APOBEC3B is associated with favorable prognosis of patients with UTUC.

Pan-Cancer Analysis for Immune Cell Infiltration and Mutational Signatures Using Non-Negative Canonical Correlation Analysis

Mutational signatures indicate the mutational processes and substitution patterns in cancer cell genomes. However, the functional consequences of mutational signatures remain unclear, and there have been no comprehensive systematic studies to examine the relationships between the mutational signatures and the immune cell infiltration. Here, the relationship between mutational signatures and immune cell infiltration using non-negative canonical correlation analysis based on 8927 patients across 25 tumor types was investigated. By inspecting mutational signatures with the maximal coefficients determined by the non-negative canonical correlation analysis, the study identified mutational signatures related to immune cell infiltration composed of tumor microenvironments. The analysis was validated by showing that the genes associated with the identified mutational signatures were linked to overall survival by a Kaplan–Meier curve and a log-rank test and were mainly related to immunity by gene set enrichment analysis. These results will help expand our knowledge of tumor biology and recognize the functional roles and associations of immune systems with mutational signatures.

The impact of rare germline variants on human somatic mutation processes

Somatic mutations are an inevitable component of ageing and the most important cause of cancer. The rates and types of somatic mutation vary across individuals, but relatively few inherited influences on mutation processes are known. We perform a gene-based rare variant association study with diverse mutational processes, using human cancer genomes from over 11,000 individuals of European ancestry. By combining burden and variance tests, we identify 207 associations involving 15 somatic mutational phenotypes and 42 genes that replicated in an independent data set at a false discovery rate of 1%. We associate rare inherited deleterious variants in genes such as MSH3, EXO1, SETD2, and MTOR with two phenotypically different forms of DNA mismatch repair deficiency, and variants in genes such as EXO1, PAXIP1, RIF1, and WRN with deficiency in homologous recombination repair. In addition, we identify associations with other mutational processes, such as APEX1 with APOBEC-signature mutagenesis. Many of the genes interact with each other and with known mutator genes within cellular sub-networks. Considered collectively, damaging variants in the identified genes are prevalent in the population. We suggest that rare germline variation in diverse genes commonly impacts mutational processes in somatic cells.

Association between germline variants and somatic mutations in colorectal cancer

Colorectal cancer (CRC) is a heterogeneous disease with evidence of distinct tumor types that develop through different somatically altered pathways. To better understand the impact of the host genome on somatically mutated genes and pathways, we assessed associations of germline variations with somatic events via two complementary approaches. We first analyzed the association between individual germline genetic variants and the presence of non-silent somatic mutations in genes in 1375 CRC cases with genome-wide SNPs data and a tumor sequencing panel targeting 205 genes. In the second analysis, we tested if germline variants located within previously identified regions of somatic allelic imbalance were associated with overall CRC risk using summary statistics from a recent large scale GWAS (n≃125 k CRC cases and controls). The first analysis revealed that a variant (rs78963230) located within a CNA region associated with TLR3 was also associated with a non-silent mutation within gene FBXW7. In the secondary analysis, the variant rs2302274 located in CDX1/PDGFRB frequently gained/lost in colorectal tumors was associated with overall CRC risk (OR = 0.96, p = 7.50e-7). In summary, we demonstrate that an integrative analysis of somatic and germline variation can lead to new insights about CRC.

Induction of APOBEC3-mediated genomic damage in urothelium implicates BK polyomavirus (BKPyV) as a hit-and-run driver for bladder cancer

Limited understanding of bladder cancer aetiopathology hampers progress in reducing incidence. Mutational signatures show the anti-viral apolipoprotein B mRNA editing enzyme catalytic polypeptide (APOBEC) enzymes are responsible for the preponderance of mutations in bladder tumour genomes, but no causative viral agent has been identified. BK polyomavirus (BKPyV) is a common childhood infection that remains latent in the adult kidney, where reactivation leads to viruria. This study provides missing mechanistic evidence linking reactivated BKPyV-infection to bladder cancer risk. We used a mitotically-quiescent, functionally-differentiated model of normal human urothelium to examine BKPyV-infection. BKPyV-infection led to significantly elevated APOBEC3A and APOBEC3B protein, increased deaminase activity and greater numbers of apurinic/apyrimidinic sites in the host urothelial genome. BKPyV Large T antigen (LT-Ag) stimulated re-entry from G0 into the cell cycle through inhibition of retinoblastoma protein and activation of EZH2, E2F1 and FOXM1, with cells arresting in G2. The single-stranded DNA displacement loops formed in urothelial cells during BKPyV-infection interacted with LT-Ag to provide a substrate for APOBEC3-activity. Addition of interferon gamma (IFNγ) to infected urothelium suppressed expression of the viral genome. These results support reactivated BKPyV infections in adults as a risk factor for bladder cancer in immune-insufficient populations.

Mutagenic Activity of AID/APOBEC Deaminases in Antiviral Defense and Carcinogenesis

Proteins of the AID/APOBEC family are capable of cytidine deamination in nucleic acids forming uracil. These enzymes are involved in mRNA editing, protection against viruses, the introduction of point mutations into DNA during somatic hypermutation, and antibody isotype switching. Since these deaminases, especially AID, are potent mutagens, their expression, activity, and specificity are regulated by several intracellular mechanisms. In this review, we discuss the mechanisms of impaired expression and activation of AID/APOBEC proteins in human tumors and their role in carcinogenesis and tumor progression. Also, the diagnostic and potential therapeutic value of increased expression of AID/APOBEC in different types of tumors is analyzed. We assume that in the case of solid tumors, increased expression of endogenous deaminases can serve as a marker of response to immunotherapy since multiple point mutations in host DNA could lead to amino acid substitutions in tumor proteins and thereby increase the frequency of neoepitopes.

Somatic mutational profiles and germline polygenic risk scores in human cancer

BACKGROUND

The mutational profile of cancer reflects the activity of the mutagenic processes which have been operative throughout the lineage of the cancer cell. These processes leave characteristic profiles of somatic mutations called mutational signatures. Mutational signatures, including single-base substitution (SBS) signatures, may reflect the effects of exogenous or endogenous exposures.

METHODS

We used polygenic risk scores (PRS) to summarize common germline variation associated with cancer risk and other cancer-related traits and examined the association between somatic mutational profiles and germline PRS in 12 cancer types from The Cancer Genome Atlas. Somatic mutational profiles were constructed from whole-exome sequencing data of primary tumors. PRS were calculated for the 12 selected cancer types and 9 non-cancer traits, including cancer risk determinants, hormonal factors, and immune-mediated inflammatory diseases, using germline genetic data and published summary statistics from genome-wide association studies.

RESULTS

We found 17 statistically significant associations between somatic mutational profiles and germline PRS after Bonferroni correction (p < 3.15 × 10^-5), including positive associations between germline inflammatory bowel disease PRS and number of somatic mutations attributed to signature SBS1 in prostate cancer and APOBEC-related signatures in breast cancer. Positive associations were also found between age at menarche PRS and mutation counts of SBS1 in overall and estrogen receptor-positive breast cancer. Consistent with prior studies that found an inverse association between the pubertal development PRS and risk of prostate cancer, likely reflecting hormone-related mechanisms, we found an inverse association between age at menarche PRS and mutation counts of SBS1 in prostate cancer. Inverse associations were also found between several cancer PRS and tumor mutation counts.

CONCLUSIONS

Our analysis suggests that there are robust associations between tumor somatic mutational profiles and germline PRS. These may reflect the mechanisms through hormone regulation and immune responses that contribute to cancer etiology and drive cancer progression.

Comparison of somatic mutation landscapes in Chinese versus European breast cancer patients

Recent genomic studies suggest that Asian breast cancer (BC) may have distinct somatic features; however, most comparisons of BC genomic features across populations did not account for differences in age, subtype, and sequencing methods. In this study, we analyzed whole-exome sequencing (WES) data to characterize somatic copy number alterations (SCNAs) and mutation profiles in 98 Hong Kong BC (HKBC) patients and compared with those from The Cancer Genome Atlas of European ancestry (TCGA-EA, N = 686), which had similar distributions of age at diagnosis and PAM50 subtypes as in HKBC. We developed a two-sample Poisson model to compare driver gene selection pressure, which reflects the effect sizes of cancer driver genes, while accounting for differences in sample size, sequencing platforms, depths, and mutation calling methods. We found that somatic mutation and SCNA profiles were overall very similar between HKBC and TCGA-EA. The selection pressure for small insertions and deletions (indels) in GATA3 (false discovery rate (FDR) corrected p < 0.01) and single-nucleotide variants (SNVs) in TP53 (nominal p = 0.02, FDR corrected p = 0.28) was lower in HKBC than in TCGA-EA. Among the 13 signatures of single-base substitutions (SBS) that are common in BC, we found a suggestively higher contribution of SBS18 and a lower contribution of SBS1 in HKBC than in TCGA-EA, while the two APOBEC-induced signatures showed similar prevalence. Our results suggest that the genomic landscape of BC was largely very similar between HKBC and TCGA-EA, despite suggestive differences in some driver genes and mutational signatures that warrant future investigations in large and diverse Asian populations.

Transcriptional repression and apoptosis influence the effect of APOBEC3A/3B functional polymorphisms on biliary tract cancer risk

APOBEC3-related somatic mutations are predominant in biliary tract cancers (BTCs). We aimed to elucidate the roles of APOBEC3A/3B functional polymorphisms and their influencing factors on the development of cholangiocarcinoma and gallbladder cancer (GBC). Polymorphisms at the promoter regions of APOBEC3A and APOBEC3B were genotyped in 3231 participants using quantitative PCR. Dual-luciferase reporter assay was applied to investigate the promoter activity. The difference in gene accessibility between cholangiocarcinoma cells and GBC cells was analyzed through single-cell transposase accessible chromatin sequencing. The effect of APOBEC3A on apoptosis was examined by cytometry. It's found that rs2267401-G at the APOBEC3B promoter decreases cholangiocarcinoma risk (age-, gender-adjusted odds ratio [AOR], 0.69; 95% confidence interval [CI], 0.51-0.94) but increases GBC risk (AOR, 2.04; 95% CI, 1.35-3.10). rs2267401-G confers a decreased APOBEC3B promoter activity in cholangiocarcinoma cells but an increased activity in GBC cells, possibly because the transcriptional repressor TFAP2A is over-expressed in cholangiocarcinoma. Tumor necrosis factor-α (TNF-α) increases the level of APOBEC3B via inhibiting TFAP2A expression rather than directly increasing the accessibility of APOBEC3B promoter. APOBEC3A promoter rs12157810-C decreased the risks of cholangiocarcinoma and GBC, with an AOR (95% CI) of 0.80 (0.66-0.97) and 0.75 (0.59-0.95), respectively. rs12157810-C upregulated the promoter activity in both cholangiocarcinoma and GBC cells. TNF-α upregulated the activity of the APOBEC3A promoter with rs12157810-C via increasing the accessibility of Ets-1 p68. APOBEC3A overexpression attenuates cancer evolution by causing apoptosis, in contrast to APOBEC3B. The heterogeneity in the transcriptional regulation of APOBEC3B affects the evolutionary potential of cancer cells in the inflammatory microenvironment. This article is protected by copyright. All rights reserved.

APOBEC mediated mutagenesis drives genomic heterogeneity in endometriosis

Endometriosis is a benign gynecologic condition, acting as a precursor of certain histological subtypes of ovarian cancers. The epithelial cells of endometriotic tissues and normal uterine endometrium accumulated somatic mutations in cancer-associated genes such as phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) and Kirsten rat sarcoma (KRAS) proto-oncogene. To determine the genomic characteristic of endometriotic epithelial cells and normal uterine endometrium and to identify the predominant mutational process acting on them, we studied the somatic mutation profiles obtained from whole exome sequencing of 14 endometriotic epithelium and 11 normal uterine endometrium tissues and classified them into mutational signatures. We observed that single base substitutions 2/13 (SBS), attributed to Apolipoprotein B mRNA Editing Enzyme Catalytic Subunit (APOBEC) induced mutagenesis, were significant in endometriotic tissues, but not in the normal uterine endometrium. Additionally, the larger number and wider allele frequency distribution of APOBEC signature mutations, compared to cancer-associated driver mutations in endometriotic epithelium suggested APOBEC mutagenesis as an important source of mutational burden and heterogeneity in endometriosis. Further, the relative risk of enriched APOBEC signature mutations was higher in endometriosis patients who were carriers of APOBEC3A/3B germline deletion, a common polymorphism in East Asians which involves the complete loss of APOBEC3B coding region. Our results illustrate the significance of APOBEC induced mutagenesis in driving the genomic heterogeneity of endometriosis.

Impact of the APOBEC3A/B deletion polymorphism on risk of ovarian cancer

A germline 29.5-kb deletion variant removes the 3’ end of the APOBEC3A gene and a large part of APOBEC3B, creating a hybrid gene that has been linked to increased APOBEC3 activity and DNA damage in human cancers. We genotyped the APOBEC3A/B deletion in hospital-based samples of 1398 Norwegian epithelial ovarian cancer patients without detected BRCA1/2 germline mutations and compared to 1,918 healthy female controls, to assess the potential cancer risk associated with the deletion. We observed an association between APOBEC3A/B status and reduced risk for ovarian cancer (OR = 0.75; CI = 0.61–0.91; p = 0.003) applying the dominant model. Similar results were found in other models. The association was observed both in non-serous and serous cases (dominant model: OR = 0.69; CI = 0.50–0.95; p = 0.018 and OR = 0.77; CI = 0.62–0.96; p = 0.019, respectively) as well as within high-grade serous cases (dominant model: OR = 0.79; CI = 0.59–1.05). For validation purposes, we mined an available large multinational GWAS-based data set of > 18,000 cases and > 26,000 controls for SNP rs12628403, known to be in linkage disequilibrium with the APOBEC3A/B deletion. We found a non-significant trend for SNP rs12628403 being linked to reduced risk of ovarian cancer in general and similar trends for all subtypes. For clear cell cancers, the risk reduction reached significance (OR = 0.85; CI = 0.69–1.00).

Cancer Evo-Dev: A Theory of Inflammation-Induced Oncogenesis

Chronic inflammation is a prerequisite for the development of cancers. Here, we present the framework of a novel theory termed as Cancer Evolution-Development (Cancer Evo-Dev) based on the current understanding of inflammation-related carcinogenesis, especially hepatocarcinogenesis induced by chronic infection with hepatitis B virus. The interaction between genetic predispositions and environmental exposures, such as viral infection, maintains chronic non-resolving inflammation. Pollution, metabolic syndrome, physical inactivity, ageing, and adverse psychosocial exposure also increase the risk of cancer via inducing chronic low-grade smoldering inflammation. Under the microenvironment of non-resolving inflammation, pro-inflammatory factors facilitate the generation of somatic mutations and viral mutations by inducing the imbalance between the mutagenic forces such as cytidine deaminases and mutation-correcting forces including uracil-DNA glycosylase. Most cells with somatic mutations and mutated viruses are eliminated in survival competition. Only a small percentage of mutated cells survive, adapt to the hostile environment, retro-differentiate, and function as cancer-initiating cells via altering signaling pathways. These cancer-initiating cells acquire stem-ness, reprogram metabolic patterns, and affect the microenvironment. The carcinogenic process follows the law of "mutation-selection-adaptation". Chronic physical activity reduces the levels of inflammation via upregulating the activity and numbers of NK cells and lymphocytes and lengthening leukocyte telomere; downregulating proinflammatory cytokines including interleukin-6 and senescent lymphocytes especially in aged population. Anti-inflammation medication reduces the occurrence and recurrence of cancers. Targeting cancer stemness signaling pathways might lead to cancer eradication. Cancer Evo-Dev not only helps understand the mechanisms by which inflammation promotes the development of cancers, but also lays the foundation for effective prophylaxis and targeted therapy of various cancers.

Mutational signatures in esophageal squamous cell carcinoma from eight countries with varying incidence

Esophageal squamous cell carcinoma (ESCC) shows remarkable variation in incidence that is not fully explained by known lifestyle and environmental risk factors. It has been speculated that an unknown exogenous exposure(s) could be responsible. Here we combine the fields of mutational signature analysis with cancer epidemiology to study 552 ESCC genomes from eight countries with varying incidence rates. Mutational profiles were similar across all countries studied. Associations between specific mutational signatures and ESCC risk factors were identified for tobacco, alcohol, opium and germline variants, with modest impacts on mutation burden. We find no evidence of a mutational signature indicative of an exogenous exposure capable of explaining differences in ESCC incidence. Apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like (APOBEC)-associated mutational signatures single-base substitution (SBS)2 and SBS13 were present in 88% and 91% of cases, respectively, and accounted for 25% of the mutation burden on average, indicating that APOBEC activation is a crucial step in ESCC tumor development.

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 CRLpVHL, induces degradation of all 7 A3 proteins. Two independent lines of evidence supported the conclusion that the multiprotein CRLpVHL complex is necessary for A3 degradation. CRLpVHL 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.

Mutational signatures: emerging concepts, caveats and clinical applications

Whole-genome sequencing has brought the cancer genomics community into new territory. Thanks to the sheer power provided by the thousands of mutations present in each patient's cancer, we have been able to discern generic patterns of mutations, termed 'mutational signatures', that arise during tumorigenesis. These mutational signatures provide new insights into the causes of individual cancers, revealing both endogenous and exogenous factors that have influenced cancer development. This Review brings readers up to date in a field that is expanding in computational, experimental and clinical directions. We focus on recent conceptual advances, underscoring some of the caveats associated with using the mutational signature frameworks and highlighting the latest experimental insights. We conclude by bringing attention to areas that are likely to see advancements in clinical applications.

Effect of APOBEC3A functional polymorphism on renal cell carcinoma is influenced by tumor necrosis factor-α and transcriptional repressor ETS1

Human apolipoprotein B mRNA editing enzyme, catalytic polypeptide (APOBEC) 3 cytidine deaminases are the prominent drivers of somatic mutations in cancers. However, the effect of APOBEC3s functional polymorphisms on the development of renal cell carcinoma (RCC) remains unknown. Five genetic polymorphisms affecting the expression of APOBEC3A (A3A), APOBEC3B, and APOBEC4 and uracil DNA glycosylase (UNG) were genotyped in 728 RCC patients and 1500 healthy controls. The effects of tumor necrosis factor-α (TNFα) and interleukin-6 on the activity of the A3A promoter with rs12157810-A or -C in four RCC cell lines (786-O, A498, Caki2, ACHN) and two colorectal cancer cell lines (HCT116, SW620) were evaluated using dual-luciferase assays. Transcriptional repressors to the A3A promoter were identified by chromatin immunoprecipitation-quantitative PCR. The proapoptotic effect of A3A on RCC cells was evaluated using cytometry. The prognostic values of A3A and ETS1 were evaluated by the Cox regression analysis. The expressions of A3A and ETS1 were evaluated in clear cell RCC (ccRCC) specimens with different polymorphic genotypes using quantitative RT-PCR and immunohistochemistry. Of those functional polymorphisms, CC genotype at rs12157810 in the A3A promoter was significantly associated with a decreased risk of ccRCC, compared to the AA genotype (odds ratio adjusted for age and gender, 0.41, 95% confidence interval [CI], 0.28-0.57). Other polymorphic genotypes were not associated with the risk of RCC. The activity of the A3A promoter with rs12157810-C was significantly higher than that with rs12157810-A in the four RCC cell lines and two colorectal cancer cell lines. The activity of the A3A promoter with rs12157810-C was greatly up-regulated by TNFα and predominantly inhibited by a transcriptional repressor ETS1. The binding of ETS1 to the A3A promoter with rs12157810-C was looser than that with rs12157810-A. Ectopic expression of A3A significantly promoted apoptosis in ccRCC cells, rather than in colorectal cancer cells. Higher ETS1 expression predicted a favorable prognosis in ccRCC, with a hazard ratio of 0.58 (95% CI, 0.43-0.78). Rs121567810-C up-regulates the A3A promoter activity, possibly due to higher response to TNFα and looser transcriptional repression by ETS1. Up-regulation of A3A increases apoptosis, thus decreasing ccRCC risk in those carrying rs121567810-C.

Genetic Determinants of Somatic Selection of Mutational Processes in 3,566 Human Cancers

The somatic landscape of the cancer genome results from different mutational processes represented by distinct "mutational signatures." Although several mutagenic mechanisms are known to cause specific mutational signatures in cell lines, the variation of somatic mutational activities in patients, which is mostly attributed to somatic selection, is still poorly explained. Here, we introduce a quantitative trait, mutational propensity (MP), and describe an integrated method to infer genetic determinants of variations in the mutational processes in 3,566 cancers with specific underlying mechanisms. As a result, we report 2,314 candidate determinants with both significant germline and somatic effects on somatic selection of mutational processes, of which, 485 act via cancer gene expression and 1,427 act through the tumor-immune microenvironment. These data demonstrate that the genetic determinants of MPs provide complementary information to known cancer driver genes, clonal evolution, and clinical biomarkers. SIGNIFICANCE: The genetic determinants of the somatic mutational processes in cancer elucidate the biology underlying somatic selection and evolution of cancers and demonstrate complementary predictive power across cancer types.