Identification of new driver and passenger mutations within APOBEC-induced hotspot mutations in bladder cancer

Exploring the Role of Clustered Mutations in Carcinogenesis and Their Potential Clinical Implications in Cancer

Abnormal cell proliferation and growth leading to cancer primarily result from cumulative genome mutations. Single gene mutations alone do not fully explain cancer onset and progression; instead, clustered mutations-simultaneous occurrences of multiple mutations-are considered to be pivotal in cancer development and advancement. These mutations can affect different genes and pathways, resulting in cells undergoing malignant transformation with multiple functional abnormalities. Clustered mutations influence cancer growth rates, metastatic potential, and drug treatment sensitivity. This summary highlights the various types and characteristics of clustered mutations to understand their associations with carcinogenesis and discusses their potential clinical significance in cancer. As a unique mutation type, clustered mutations may involve genomic instability, DNA repair mechanism defects, and environmental exposures, potentially correlating with responsiveness to immunotherapy. Understanding the characteristics and underlying processes of clustered mutations enhances our comprehension of carcinogenesis and cancer progression, providing new diagnostic and therapeutic approaches for cancer.

Systematically Evaluating Cell-Free DNA Fragmentation Patterns for Cancer Diagnosis and Enhanced Cancer Detection via Integrating Multiple Fragmentation Patterns

Cell-free DNA (cfDNA) fragmentation patterns have immense potential for early cancer detection. However, the definition of fragmentation varies, ranging from the entire genome to specific genomic regions. These patterns have not been systematically compared, impeding broader research and practical implementation. Here, 1382 plasma cfDNA sequencing samples from 8 cancer types are collected. Considering that cfDNA within open chromatin regions is more susceptible to fragmentation, 10 fragmentation patterns within open chromatin regions as features and employed machine learning techniques to evaluate their performance are examined. All fragmentation patterns demonstrated discernible classification capabilities, with the end motif showing the highest diagnostic value for cross-validation. Combining cross and independent validation results revealed that fragmentation patterns that incorporated both fragment length and coverage information exhibited robust predictive capacities. Despite their diagnostic potential, the predictive power of these fragmentation patterns is unstable. To address this limitation, an ensemble classifier via integrating all fragmentation patterns is developed, which demonstrated notable improvements in cancer detection and tissue-of-origin determination. Further functional bioinformatics investigations on significant feature intervals in the model revealed its impressive ability to identify critical regulatory regions involved in cancer pathogenesis.

Single-cell analyses reveal evolution mimicry during the specification of breast cancer subtype

Background: The stem or progenitor antecedents confer developmental plasticity and unique cell identities to cancer cells via genetic and epigenetic programs. A comprehensive characterization and mapping of the cell-of-origin of breast cancer using novel technologies to unveil novel subtype-specific therapeutic targets is still absent. Methods: We integrated 195,144 high-quality cells from normal breast tissues and 406,501 high-quality cells from primary breast cancer samples to create a large-scale single-cell atlas of human normal and cancerous breasts. Potential heterogeneous origin of malignant cells was explored by contrasting cancer cells against reference normal epithelial cells. Multi-omics analyses and both in vitro and in vivo experiments were performed to screen and validate potential subtype-specific treatment targets. Novel biomarkers of identified immune and stromal cell subpopulations were validated by immunohistochemistry in our cohort. Results: Tumor stratification based on cancer cell-of-origin patterns correlated with clinical outcomes, genomic aberrations and diverse microenvironment constitutions. We found that the luminal progenitor (LP) subtype was robustly associated with poor prognosis, genomic instability and dysfunctional immune microenvironment. However, the LP subtype patients were sensitive to neoadjuvant chemotherapy (NAC), PARP inhibitors (PARPi) and immunotherapy. The LP subtype-specific target PLK1 was investigated by both in vitro and in vivo experiments. Besides, large-scale single-cell profiling of breast cancer inspired us to identify a range of clinically relevant immune and stromal cell subpopulations, including subsets of innate lymphoid cells (ILCs), macrophages and endothelial cells. Conclusion: The present single-cell study revealed the cellular repertoire and cell-of-origin patterns of breast cancer. Combining single-cell and bulk transcriptome data, we elucidated the evolution mimicry from normal to malignant subtypes and expounded the LP subtype with vital clinical implications. Novel immune and stromal cell subpopulations of breast cancer identified in our study could be potential therapeutic targets. Taken together, Our findings lay the foundation for the precise prognostic and therapeutic stratification of breast cancer.

The aryl hydrocarbon receptor as a tumor modulator: mechanisms to therapy

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that is widely recognized to play important, but complex, modulatory roles in a variety of tumor types. In this review, we comprehensively summarize the increasingly controversial role of AhR as a tumor regulator and the mechanisms by which it alters tumor progression based on the cancer cell type. Finally, we discuss new and emerging strategies to therapeutically modulate AhR, focusing on novel agents that hold promise in current human clinical trials as well as existing FDA-approved drugs that could potentially be repurposed for cancer therapy.

Targeting aryl hydrocarbon receptor to prevent cancer in barrier organs

The skin, lung, and gut are important barrier organs that control how the body reacts to environmental stressors such as ultraviolet (UV) radiation, air pollutants, dietary components, and microorganisms. The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that plays an important role in maintaining homeostasis of barrier organs. AhR was initially discovered as a receptor for environmental chemical carcinogens such as polycyclic aromatic hydrocarbons (PAHs). Activation of AhR pathways by PAHs leads to increased DNA damage and mutations which ultimately lead to carcinogenesis. Ongoing evidence reveals an ever-expanding role of AhR. Recently, AhR has been linked to immune systems by the interaction with the development of natural killer (NK) cells, regulatory T (Treg) cells, and T helper 17 (Th17) cells, as well as the production of immunosuppressive cytokines. However, the role of AhR in carcinogenesis is not as straightforward as we initially thought. Although AhR activation has been shown to promote carcinogenesis in some studies, others suggest that it may act as a tumor suppressor. In this review, we aim to explore the role of AhR in the development of cancer that originates from barrier organs. We also examined the preclinical efficacy data of AhR agonists and antagonists on carcinogenesis to determine whether AhR modulation can be a viable option for cancer chemoprevention.

Proteogenomic Characterization of Bladder Cancer Reveals Sensitivity to Apoptosis Induced by Tumor Necrosis Factor-related Apoptosis-inducing Ligand in FGFR3-mutated Tumors

BACKGROUND

Molecular understanding of muscle-invasive (MIBC) and non-muscle-invasive (NMIBC) bladder cancer is currently based primarily on transcriptomic and genomic analyses.

OBJECTIVE

To conduct proteogenomic analyses to gain insights into bladder cancer (BC) heterogeneity and identify underlying processes specific to tumor subgroups and therapeutic outcomes.

DESIGN, SETTING, AND PARTICIPANTS

Proteomic data were obtained for 40 MIBC and 23 NMIBC cases for which transcriptomic and genomic data were already available. Four BC-derived cell lines harboring FGFR3 alterations were tested with interventions.

INTERVENTION

Recombinant tumor necrosis factor-related apoptosis-inducing ligand (TRAIL), second mitochondrial-derived activator of caspases mimetic (birinapant), pan-FGFR inhibitor (erdafitinib), and FGFR3 knockdown.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Proteomic groups from unsupervised analyses (uPGs) were characterized using clinicopathological, proteomic, genomic, transcriptomic, and pathway enrichment analyses. Additional enrichment analyses were performed for FGFR3-mutated tumors. Treatment effects on cell viability for FGFR3-altered cell lines were evaluated. Synergistic treatment effects were evaluated using the zero interaction potency model.

RESULTS AND LIMITATIONS

Five uPGs, covering both NMIBC and MIBC, were identified and bore coarse-grained similarity to transcriptomic subtypes underlying common features of these different entities; uPG-E was associated with the Ta pathway and enriched in FGFR3 mutations. Our analyses also highlighted enrichment of proteins involved in apoptosis in FGFR3-mutated tumors, not captured through transcriptomics. Genetic and pharmacological inhibition demonstrated that FGFR3 activation regulates TRAIL receptor expression and sensitizes cells to TRAIL-mediated apoptosis, further increased by combination with birinapant.

CONCLUSIONS

This proteogenomic study provides a comprehensive resource for investigating NMIBC and MIBC heterogeneity and highlights the potential of TRAIL-induced apoptosis as a treatment option for FGFR3-mutated bladder tumors, warranting a clinical investigation.

PATIENT SUMMARY

We integrated proteomics, genomics, and transcriptomics to refine molecular classification of bladder cancer, which, combined with clinical and pathological classification, should lead to more appropriate management of patients. Moreover, we identified new biological processes altered in FGFR3-mutated tumors and showed that inducing apoptosis represents a new potential therapeutic option.

Whole-genome mapping of APOBEC mutagenesis in metastatic urothelial carcinoma identifies driver hotspot mutations and a novel mutational signature

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) enzymes mutate specific DNA sequences and hairpin-loop structures, challenging the distinction between passenger and driver hotspot mutations. Here, we characterized 115 whole genomes of metastatic urothelial carcinoma (mUC) to identify APOBEC mutagenic hotspot drivers. APOBEC-associated mutations were detected in 92% of mUCs and were equally distributed across the genome, while APOBEC hotspot mutations (ApoHMs) were enriched in open chromatin. Hairpin loops were frequent targets of didymi (twins in Greek), two hotspot mutations characterized by the APOBEC SBS2 signature, in conjunction with an uncharacterized mutational context (Ap[C>T]). Next, we developed a statistical framework that identified ApoHMs as drivers in coding and non-coding genomic regions of mUCs. Our results and statistical framework were validated in independent cohorts of 23 non-metastatic UCs and 3,744 samples of 17 metastatic cancers, identifying cancer-type-specific drivers. Our study highlights the role of APOBEC in cancer development and may contribute to developing novel targeted therapy options for APOBEC-driven cancers.

Identification and segregation of genes with improved recurrent neural network trained with optimal gene level and mutation level features

Even though many different approaches have been employed to address the complex mutational heterogeneity of cancer, finding driver genes is still problematic since other genomic factors cannot be fully integrated for combined analyses. This research paper presents a novel gene identification and segregation model with five key processes (a) pre-processing, (b) treatment of class imbalances, (c) feature extraction, (d) feature selection, and (e) gene classification. To increase the data quality, the gathered initial information is first pre-processed utilizing data cleaning and data normalization. This turns the raw data into something that is both useful and effective. In actuality, the sample is skewed against drivers because passenger mutation markers appear in proportionally less instances than drivers do. To address the Class Imbalance Problem, improved K-Means + SMOTE are applied to the preprocessed data. The most crucial characteristics, including those at the gene and mutation levels, are then extracted from the balanced dataset. To lessen the computational load in terms of time, the best features from the retrieved features are selected using Forensic interpretation tailored hunger food search optimization (FIHFSO). The ideal features are used to train the deep learning classifier that conducts the separation procedure. In this research, an Improved Recurrent Neural Network (I-RNN) is used to make a final decision about genes. At 90% of learning percentage, the accuracy of the proposed method achieves 0.98% of 0.83, 0.81, 0.65, 0.80, 0.92 and 0.63% which is compared to the other methods like HGS, FBIO, AOA, AO, GOA and PRO respectively.

FGFR inhibition augments anti-PD-1 efficacy in murine FGFR3-mutant bladder cancer by abrogating immunosuppression

The combination of targeted therapy with immune checkpoint inhibition (ICI) is an area of intense interest. We studied the interaction of fibroblast growth factor receptor (FGFR) inhibition with ICI in urothelial carcinoma (UC) of the bladder, in which FGFR3 is altered in 50% of cases. Using an FGFR3-driven, Trp53-mutant genetically engineered murine model (UPFL), we demonstrate that UPFL tumors recapitulate the histology and molecular subtype of their FGFR3-altered human counterparts. Additionally, UPFL1 allografts exhibit hyperprogression to ICI associated with an expansion of T regulatory cells (Tregs). Erdafitinib blocked Treg proliferation in vitro, while in vivo ICI-induced Treg expansion was fully abrogated by FGFR inhibition. Combined erdafitinib and ICI resulted in high therapeutic efficacy. In aggregate, our work establishes that, in mice, co-alteration of FGFR3 and Trp53 results in high-grade, non-muscle-invasive UC and presents a previously underappreciated role for FGFR inhibition in blocking ICI-induced Treg expansion.

Hotspot propensity across mutational processes

The sparsity of mutations observed across tumours hinders our ability to study mutation rate variability at nucleotide resolution. To circumvent this, here we investigated the propensity of mutational processes to form mutational hotspots as a readout of their mutation rate variability at single base resolution. Mutational signatures 1 and 17 have the highest hotspot propensity (5-78 times higher than other processes). After accounting for trinucleotide mutational probabilities, sequence composition and mutational heterogeneity at 10 Kbp, most (94-95%) signature 17 hotspots remain unexplained, suggesting a significant role of local genomic features. For signature 1, the inclusion of genome-wide distribution of methylated CpG sites into models can explain most (80-100%) of the hotspot propensity. There is an increased hotspot propensity of signature 1 in normal tissues and de novo germline mutations. We demonstrate that hotspot propensity is a useful readout to assess the accuracy of mutation rate models at nucleotide resolution. This new approach and the findings derived from it open up new avenues for a range of somatic and germline studies investigating and modelling mutagenesis.

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.

The impact of mutational clonality in predicting the response to immune checkpoint inhibitors in advanced urothelial cancer

Immune checkpoint inhibitors (ICI) have revolutionized cancer treatment and can result in complete remissions even at advanced stages of the disease. However, only a small fraction of patients respond to the treatment. To better understand which factors drive clinical benefit, we have generated whole exome and RNA sequencing data from 27 advanced urothelial carcinoma patients treated with anti-PD-(L)1 monoclonal antibodies. We assessed the influence on the response of non-synonymous mutations (tumor mutational burden or TMB), clonal and subclonal mutations, neoantigen load and various gene expression markers. We found that although TMB is significantly associated with response, this effect can be mostly explained by clonal mutations, present in all cancer cells. This trend was validated in an additional cohort. Additionally, we found that responders with few clonal mutations had abnormally high levels of T and B cell immune markers, suggesting that a high immune cell infiltration signature could be a better predictive biomarker for this subset of patients. Our results support the idea that highly clonal cancers are more likely to respond to ICI and suggest that non-additive effects of different signatures should be considered for predictive models.

Mutation Hotspots Found in Bladder Cancer Aid Prediction of Carcinogenic Risk in Normal Urothelium

More than 80,000 new cases of bladder cancer are estimated to be diagnosed in 2023. However, the 5-year survival rate for bladder cancer has not changed in decades, highlighting the need for prevention. Numerous cancer-causing mutations are present in the urothelium long before signs of cancer arise. Mutation hotspots in cancer-driving genes were identified in non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC) tumor samples. Mutation burden within the hotspot regions was measured in normal urothelium with a low and high risk of cancer. A significant correlation was found between the mutation burden in normal urothelium and bladder cancer tissue within the hotspot regions. A combination of measured hotspot burden and personal risk factors was used to fit machine learning classification models. The efficacy of each model to differentiate between adjacent benign urothelium from bladder cancer patients and normal urothelium from healthy donors was measured. A random forest model using a combination of personal risk factors and mutations within MIBC hotspots yielded the highest AUC of 0.9286 for the prediction of high- vs. low-risk normal urothelium. Currently, there are no effective biomarkers to assess subclinical field disease and early carcinogenic progression in the bladder. Our findings demonstrate novel differences in mutation hotspots in NMIBC and MIBC and provide the first evidence for mutation hotspots to aid in the assessment of cancer risk in the normal urothelium. Early risk assessment and identification of patients at high risk of bladder cancer before the clinical presentation of the disease can pave the way for targeted personalized preventative therapy.

Epigenomic mapping identifies an enhancer repertoire that regulates cell identity in bladder cancer through distinct transcription factor networks

Muscle-invasive bladder cancer (BLCA) is an aggressive disease. Consensus BLCA transcriptomic subtypes have been proposed, with two major Luminal and Basal subgroups, presenting distinct molecular and clinical characteristics. However, how these distinct subtypes are regulated remains unclear. We hypothesized that epigenetic activation of distinct super-enhancers could drive the transcriptional programs of BLCA subtypes. Through integrated RNA-sequencing and epigenomic profiling of histone marks in primary tumours, cancer cell lines, and normal human urothelia, we established the first integrated epigenetic map of BLCA and demonstrated the link between subtype and epigenetic control. We identified the repertoire of activated super-enhancers and highlighted Basal, Luminal and Normal-associated SEs. We revealed super-enhancer-regulated networks of candidate master transcription factors for Luminal and Basal subgroups including FOXA1 and ZBED2, respectively. FOXA1 CRISPR-Cas9 mutation triggered a shift from Luminal to Basal phenotype, confirming its role in Luminal identity regulation and induced ZBED2 overexpression. In parallel, we showed that both FOXA1 and ZBED2 play concordant roles in preventing inflammatory response in cancer cells through STAT2 inhibition. Our study furthers the understanding of epigenetic regulation of muscle-invasive BLCA and identifies a co-regulated network of super-enhancers and associated transcription factors providing potential targets for the treatment of this aggressive disease.

OKseqHMM: a genome-wide replication fork directionality analysis toolkit

During each cell division, tens of thousands of DNA replication origins are co-ordinately activated to ensure the complete duplication of the human genome. However, replication fork progression can be challenged by many factors, including co-directional and head-on transcription-replication conflicts (TRC). Head-on TRCs are more dangerous for genome integrity. To study the direction of replication fork movement and TRCs, we developed a bioinformatics toolkit called OKseqHMM (https://github.com/CL-CHEN-Lab/OK-Seq, https://doi.org/10.5281/zenodo.7428883). Then, we used OKseqHMM to analyse a large number of datasets obtained by Okazaki fragment sequencing to directly measure the genome-wide replication fork directionality (RFD) and to accurately predict replication initiation and termination at a fine resolution in organisms including yeast, mouse and human. We also successfully applied our analysis to other genome-wide sequencing techniques that also contain RFD information (e.g. eSPAN, TrAEL-seq). Our toolkit can be used to predict replication initiation and fork progression direction genome-wide in a wide range of cell models and growth conditions. Comparing the replication and transcription directions allows identifying loci at risk of TRCs, particularly head-on TRCs, and investigating their role in genome instability by checking DNA damage data, which is of prime importance for human health.

Properties of non-coding mutation hotspots as urinary biomarkers for bladder cancer detection

Mutations at specific hotspots in non-coding regions of ADGRG6, PLEKHS1, WDR74, TBC1D12 and LEPROTL1 frequently occur in bladder cancer (BC). These mutations could function as biomarkers for the non-invasive detection of BC but this remains largely unexplored. Massively-parallel sequencing of non-coding hotspots was applied to 884 urine cell pellet DNAs: 591 from haematuria clinic patients (165 BCs, 426 non-BCs) and 293 from non-muscle invasive BC surveillance patients (29 with recurrence). Urine samples from 142 non-BC haematuria clinic patients were used to optimise variant calling. Non-coding mutations are readily detectable in the urine of BC patients and undetectable, or present at much lower frequencies, in the absence of BC. The mutations can be used to detect incident BC with 66% sensitivity (95% CI 58-75) at 92% specificity (95% CI 88-95) and recurrent disease with 55% sensitivity (95% CI 36-74) at 85% specificity (95% CI 80-89%) using a 2% variant allele frequency threshold. In the NMIBC surveillance setting, the detection of non-coding mutations in urine in the absence of clinically detectable disease was associated with an increased relative risk of future recurrence (RR = 4.62 (95% CI 3.75-5.48)). As urinary biomarkers, non-coding hotspot mutations behave similarly to driver mutations in BC-associated genes and could be included in biomarker panels for BC detection.

In-depth analysis of the interactions of various aryl hydrocarbon receptor ligands from a computational perspective

Aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that acts as a machinery that controls the expression of many genes, including cytochrome P450 CYP1A1, CYP1A2 and CYP1B1. It plays a principal role in numerous biological and toxicological functions, making it a promising target for developing therapeutic agents. Several novel small molecules targeting the AhR signaling pathway are currently under investigation as antitumor agents. Some have already advanced into clinical trials in patients with various tumors. Activation of AhR by diverse chemicals either endogenous or exogenous is initiated by the binding of these ligands to the PAS-B domain, which modulates AhR functions. There is, however, limited information about how various ligands interact with the PAS-B domain for activating or inhibiting the AhR. To better understand the mode of action of AhR agonists/antagonists. The current work proposes a combination of several computational tools to build dynamical models for the PAS-B domain bound to different ligands in mouse and human. Our findings reveal the essential roles of specific PAS-B residues (e.g., S365, V381& Q383), which mediate the AhR ligand-binding process. Our results also explain how these residues regulate the promiscuity of AhR in accommodating various chemicals in its binding PAS-B ligand-binding pocket.

APOBEC3: Friend or Foe in Human Papillomavirus Infection and Oncogenesis?

Human papillomavirus (HPV) infection is a causative agent of multiple human cancers, including cervical and head and neck cancers. In these HPV-positive tumors, somatic mutations are caused by aberrant activation of DNA mutators such as members of the apolipoprotein B messenger RNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) family of cytidine deaminases. APOBEC3 proteins are most notable for their restriction of various viruses, including anti-HPV activity. However, the potential role of APOBEC3 proteins in HPV-induced cancer progression has recently garnered significant attention. Ongoing research stems from the observations that elevated APOBEC3 expression is driven by HPV oncogene expression and that APOBEC3 activity is likely a significant contributor to somatic mutagenesis in HPV-positive cancers. This review focuses on recent advances in the study of APOBEC3 proteins and their roles in HPV infection and HPV-driven oncogenesis. Further, we discuss critical gaps and unanswered questions in our understanding of APOBEC3 in virus-associated cancers.

Identification of a signature of evolutionarily conserved stress-induced mutagenesis in cancer

The clustering of mutations observed in cancer cells is reminiscent of the stress-induced mutagenesis (SIM) response in bacteria. Bacteria deploy SIM when faced with DNA double-strand breaks in the presence of conditions that elicit an SOS response. SIM employs DinB, the evolutionary precursor to human trans-lesion synthesis (TLS) error-prone polymerases, and results in mutations concentrated around DNA double-strand breaks with an abundance that decays with distance. We performed a quantitative study on single nucleotide variant calls for whole-genome sequencing data from 1950 tumors, non-inherited mutations from 129 normal samples, and acquired mutations in 3 cell line models of stress-induced adaptive mutation. We introduce statistical methods to identify mutational clusters, quantify their shapes and tease out the potential mechanism that produced them. Our results show that mutations in both normal and cancer samples are indeed clustered and have shapes indicative of SIM. Clusters in normal samples occur more often in the same genomic location across samples than in cancer suggesting loss of regulation over the mutational process during carcinogenesis. Additionally, the signatures of TLS contribute the most to mutational cluster formation in both patient samples as well as experimental models of SIM. Furthermore, a measure of cluster shape heterogeneity was associated with cancer patient survival with a hazard ratio of 5.744 (Cox Proportional Hazard Regression, 95% CI: 1.824–18.09). Our results support the conclusion that the ancient and evolutionary-conserved adaptive mutation response found in bacteria is a source of genomic instability in cancer. Biological adaptation through SIM might explain the ability of tumors to evolve in the face of strong selective pressures such as treatment and suggests that the conventional ‘hit it hard’ approaches to therapy could prove themselves counterproductive.

Transcriptional and genetic alterations of cuproptosis-related genes correlated to malignancy and immune-infiltrate of esophageal carcinoma

Esophageal carcinoma (ESCA) is a common type of cancer with high mortality. Cuproptosis is a new type of cell death and is characterized by the dependence on mitochondrial respiration and protein lipoylation. However, the potential roles of cuproptosis-related genes (CRGs) in ESCA remain elusive. Here, we systematically assessed the transcriptional and genetic alterations of CRGs in ESCA. We identified a CRGs signature for ESCA patients. A 6-CRGs signature was constructed by the least absolute shrinkage and selection operator (LASSO) regression analysis along with the univariate cox regression analysis and differential genes analysis. The CRGs score could significantly stratify ESCA patients’ survival and a high CRGs score was significantly correlated with worse overall survival. Moreover, higher CRGs score indicated higher pathology grades and aberrant cell adhesion, possibly via the PI3K-AKT pathway, which could also underly their increased sensitivity to PI3K-AKT pathway inhibitors. In addition, patients with high CRGs tend to hold more mutation load and abnormal APOBEC mutation. Notably, a higher CRGs score was anomalously associated with more immune infiltration, which could explain its malignancy by increased PD-L1 stability and a higher proportion of bystander T cells. In conclusion, our report revealed the significance of cuproptosis in ESCA and may have therapeutic potential in activating the bystander T cells.

Recurrent exon-deleting activating mutations in AHR act as drivers of urinary tract cancer

Bladder cancer has a high recurrence rate and low survival of advanced stage patients. Few genetic drivers of bladder cancer have thus far been identified. We performed in-depth structural variant analysis on whole-genome sequencing data of 206 metastasized urinary tract cancers. In ~ 10% of the patients, we identified recurrent in-frame deletions of exons 8 and 9 in the aryl hydrocarbon receptor gene (AHR^Δe8-9), which codes for a ligand-activated transcription factor. Pan-cancer analyses show that AHR^Δe8-9 is highly specific to urinary tract cancer and mutually exclusive with other bladder cancer drivers. The ligand-binding domain of the AHR^Δe8-9 protein is disrupted and we show that this results in ligand-independent AHR-pathway activation. In bladder organoids, AHR^Δe8-9 induces a transformed phenotype that is characterized by upregulation of AHR target genes, downregulation of differentiation markers and upregulation of genes associated with stemness and urothelial cancer. Furthermore, AHR^Δe8-9 expression results in anchorage independent growth of bladder organoids, indicating tumorigenic potential. DNA-binding deficient AHR^Δe8-9 fails to induce transformation, suggesting a role for AHR target genes in the acquisition of the oncogenic phenotype. In conclusion, we show that AHR^Δe8-9 is a novel driver of urinary tract cancer and that the AHR pathway could be an interesting therapeutic target.

Combined exome and transcriptome sequencing of non-muscle-invasive bladder cancer: associations between genomic changes, expression subtypes, and clinical outcomes

BACKGROUND

Three-quarters of bladder cancer patients present with early-stage disease (non-muscle-invasive bladder cancer, NMIBC, UICC TNM stages Ta, T1 and Tis); however, most next-generation sequencing studies to date have concentrated on later-stage disease (muscle-invasive BC, stages T2+). We used exome and transcriptome sequencing to comprehensively characterise NMIBCs of all grades and stages to identify prognostic genes and pathways that could facilitate treatment decisions. Tumour grading is based upon microscopy and cellular appearances (grade 1 BCs are less aggressive, and grade 3 BCs are most aggressive), and we chose to also focus on the most clinically complex NMIBC subgroup, those patients with grade 3 pathological stage T1 (G3 pT1) disease.

METHODS

Whole-exome and RNA sequencing were performed in total on 96 primary NMIBCs including 22 G1 pTa, 14 G3 pTa and 53 G3 pT1s, with both exome and RNA sequencing data generated from 75 of these individual samples. Associations between genomic alterations, expression profiles and progression-free survival (PFS) were investigated.

RESULTS

NMIBCs clustered into 3 expression subtypes with different somatic alteration characteristics. Amplifications of ARNT and ERBB2 were significant indicators of worse PFS across all NMIBCs. High APOBEC mutagenesis and high tumour mutation burden were both potential indicators of better PFS in G3pT1 NMIBCs. The expression of individual genes was not prognostic in BCG-treated G3pT1 NMIBCs; however, downregulated interferon-alpha and gamma response pathways were significantly associated with worse PFS (adjusted p-value < 0.005).

CONCLUSIONS

Multi-omic data may facilitate better prognostication and selection of therapeutic interventions in patients with G3pT1 NMIBC. These findings demonstrate the potential for improving the management of high-risk NMIBC patients and warrant further prospective validation.

EpisomiR, a New Family of miRNAs, and Its Possible Roles in Human Diseases

MicroRNAs (miRNAs) are synthesized through a canonical pathway and play a role in human diseases, such as cancers and cardiovascular, neurodegenerative, psychiatric, and chronic inflammatory diseases. The development of sequencing technologies has enabled the identification of variations in noncoding miRNAs. These miRNA variants, called isomiRs, are generated through a non-canonical pathway, by several enzymes that alter the length and sequence of miRNAs. The isomiR family is, now, expanding further to include episomiRs, which are miRNAs with different modifications. Since recent findings have shown that isomiRs reflect the cell-specific biological function of miRNAs, knowledge about episomiRs and isomiRs can, possibly, contribute to the optimization of diagnosis and therapeutic technology for precision medicine.

Kronos scRT: a uniform framework for single-cell replication timing analysis

Mammalian genomes are replicated in a cell type-specific order and in coordination with transcription and chromatin organization. Currently, single-cell replication studies require individual processing of sorted cells, yielding a limited number (<100) of cells. Here, we develop Kronos scRT, a software for single-cell Replication Timing (scRT) analysis. Kronos scRT does not require a specific platform or cell sorting, which allows investigating large datasets obtained from asynchronous cells. By applying our tool to published data as well as droplet-based single-cell whole-genome sequencing data generated in this study, we exploit scRT from thousands of cells for different mouse and human cell lines. Our results demonstrate that although genomic regions are frequently replicated around their population average RT, replication can occur stochastically throughout S phase. Altogether, Kronos scRT allows fast and comprehensive investigations of the RT programme at the single-cell resolution for both homogeneous and heterogeneous cell populations.

A scalable approach to explore DNA replication in single cells reveals that although aneuploidy does not have a major impact on the pattern of replication, different cell types and sub-populations display distinguished replication paths.

Bioinformatic Analyzes of the Association Between Upregulated Expression of JUN Gene via APOBEC-Induced FLG Gene Mutation and Prognosis of Cervical Cancer

Globally, cervical cancer (CC) is the most common malignant tumor of the female reproductive system and its incidence is only second after breast cancer. Although screening and advanced treatment strategies have improved the rates of survival, some patients with CC still die due to metastasis and drug resistance. It is considered that cancer is driven by somatic mutations, such as single nucleotide, small insertions/deletions, copy number, and structural variations, as well as epigenetic changes. Previous studies have shown that cervical intraepithelial neoplasia is associated with copy number variants (CNVs) and/or mutations in cancer-related genes. Further, CC is also related to genetic mutations. The present study analyzed the data on somatic mutations of cervical squamous cell carcinoma (CESC) in the Cancer Genome Atlas database. It was evident that the Apolipoprotein B mRNA editing enzyme-catalyzed polypeptide-like (APOBEC)-related mutation of the FLG gene can upregulate the expression of the JUN gene and ultimately lead to poor prognosis for patients with CC. Therefore, the findings of the current study provide a new direction for future treatment of CC.

Prospectively defined patterns of APOBEC3A mutagenesis are prevalent in human cancers

Mutational signatures defined by single base substitution (SBS) patterns in cancer have elucidated potential mutagenic processes that contribute to malignancy. Two prevalent mutational patterns in human cancers are attributed to the APOBEC3 cytidine deaminase enzymes. Among the seven human APOBEC3 proteins, APOBEC3A is a potent deaminase and proposed driver of cancer mutagenesis. In this study, we prospectively examine genome-wide aberrations by expressing human APOBEC3A in avian DT40 cells. From whole-genome sequencing, we detect hundreds to thousands of base substitutions per genome. The APOBEC3A signature includes widespread cytidine mutations and a unique insertion-deletion (indel) signature consisting largely of cytidine deletions. This multi-dimensional APOBEC3A signature is prevalent in human cancer genomes. Our data further reveal replication-associated mutations, the rate of stem-loop and clustered mutations, and deamination of methylated cytidines. This comprehensive signature of APOBEC3A mutagenesis is a tool for future studies and a potential biomarker for APOBEC3 activity in cancer.

Apolipoprotein B mRNA-Editing Catalytic Polypeptide-Like-Induced Protein Changes in Estrogen Receptor-Positive, Human Epidermal Growth Factor Receptor 2-Negative Breast Cancer Throughout Disease Progression

PURPOSE

Apolipoprotein B mRNA-Editing Catalytic Polypeptide-like (APOBEC) enzymes are mutagenic factors contributing to tumor progression and therapy resistance. However, the effects of APOBEC-induced protein changes have not been systematically assessed. Here, we describe the effects of APOBEC on the coding sequence in primary and metastatic estrogen receptor–positive (ER+)/human epidermal growth factor receptor 2–negative (HER2–) breast cancer (BC).

METHODS

We determined the enrichment of amino acid (AA) changes resulting from APOBEC mutagenesis in 323 primary BC tumors and 424 metastatic breast cancer (mBC) lesions via comparison with a simulated mutational genomic landscape not under selection pressure. We subsequently explored genes with recurrent APOBEC-associated AA changes and investigated the clonality of individual APOBEC-associated mutations. Using public sequencing data from an independent primary BC and mBC cohort, we further confirm our findings by reporting genes having these enriched AA changes in an APOBEC context.

RESULTS

Our analysis demonstrated that several APOBEC-derived AA changes are significantly enriched compared with a simulated AA change distribution drawn at random. Among the enriched AA changes, Glutamate(E)>Lysine(K) and Glutamate(E)>Glutamine(Q) were mostly found at hotspots in oncogenes, whereas termination codons (Glutamine[Q]>STOP[X] and Serine[S]>STOP[X]) occurred in tumor suppressor genes and, mostly, not at hotspot locations. These mutations are found in genes contributing to BC initiation, eg, introduction of termination codons in TP53, MAP3K1, and CDH1 (in lobular BC) and two oncogenic hotspots in PIK3CA (p.E542K and p.E545K). In endocrine-resistant BC, we observed APOBEC-induced termination codons at ER-modulating genes, KMT2C, ARID1A, NF1, and ZFHX3. Finally, in mBC, compared with single PIK3CA mutations, dual PIK3CA mutations occurred more frequently in an APOBEC context (Fisher's exact P < .001).

CONCLUSION

Our results show that APOBEC mutagenesis recurrently targets various known drivers of BC initiation, progression, and endocrine resistance.

APOBEC drives early, metastatic and endocrine resistant disease in ER+/HER2- breast cancer

Primary Founder Mutations in the PRKDC Gene Increase Tumor Mutation Load in Colorectal Cancer

The clonal composition of a malignant tumor strongly depends on cellular dynamics influenced by the asynchronized loss of DNA repair mechanisms. Here, our aim was to identify founder mutations leading to subsequent boosts in mutation load. The overall mutation burden in 591 colorectal cancer tumors was analyzed, including the mutation status of DNA-repair genes. The number of mutations was first determined across all patients and the proportion of genes having mutation in each percentile was ranked. Early mutations in DNA repair genes preceding a mutational expansion were designated as founder mutations. Survival analysis for gene expression was performed using microarray data with available relapse-free survival. Of the 180 genes involved in DNA repair, the top five founder mutations were in PRKDC (n = 31), ATM (n = 26), POLE (n = 18), SRCAP (n = 18), and BRCA2 (n = 15). PRKDC expression was 6.4-fold higher in tumors compared to normal samples, and higher expression led to longer relapse-free survival in 1211 patients (HR = 0.72, p = 4.4 × 10^-3). In an experimental setting, the mutational load resulting from UV radiation combined with inhibition of PRKDC was analyzed. Upon treatments, the mutational load exposed a significant two-fold increase. Our results suggest PRKDC as a new key gene driving tumor heterogeneity.

Comprehensive bioinformatics analyses of APOBECs family and identification of APOBEC3D as the unfavorable prognostic biomarker in clear cell renal cell carcinoma

Purpose: At present, how early screening for ccRCC is still a thorny issue for urologists. Probing the mechanisms underlying the development of ccRCC and finding relevant prognostic biomarkers remains crucial. Therefore, we systematically analyzed the APOBEC family in this study and identified APOBEC3D as a prognostic biomarker. Methods: In this study, based on the TCGA database, we systematically assessed the expression and prognosis of the APOBEC family and analyzed potential bioinformatic pathways. We then constructed nomograms to predict the prognosis of ccRCC patients better. Afterward, we further focused on APOBEC3D in our data on ccRCC specimens. The APOBEC3D should be extensively studied in ccRCC in the future. Results: The results showed that the APOBEC family showed the most significant changes in expression in ccRCC. The pathway enrichment analysis showed that APOBEC3 family members mainly regulated cytidine and cytosine-related processes. Subsequently, the Cox regression was used to construct prognostic signature, and validated in ICGC and GEO databases. Next, a nomogram was created integrating clinical parameters showing good predictive performance. Finally, we screened for APOBEC3D and found in our clinical sample that patients with high expression of APOBEC3D had a worse prognosis. Conclusion: Based on these results, APOBEC family members play important roles in the development of ccRCC, and APOBEC3D could serve as the biomarker for predicting patient prognosis.

Review of Experimental Studies to Improve Radiotherapy Response in Bladder Cancer: Comments and Perspectives

Bladder cancer is among the top ten most common cancer types in the world. Around 25% of all cases are muscle-invasive bladder cancer, for which the gold standard treatment in the absence of metastasis is the cystectomy. In recent years, trimodality treatment associating maximal transurethral resection and radiotherapy combined with concurrent chemotherapy is increasingly used as an organ-preserving alternative. However, the use of this treatment is still limited by the lack of biomarkers predicting tumour response and by a lack of targeted radiosensitising drugs that can improve the therapeutic index, especially by limiting side effects such as bladder fibrosis. In order to improve the bladder-preserving treatment, experimental studies addressing these main issues ought to be considered (both in vitro and in vivo studies). Following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines for systematic reviews, we conducted a literature search in PubMed on experimental studies investigating how to improve bladder cancer radiotherapy with different radiosensitising agents using a comprehensive search string. We made comments on experimental model selection, experimental design and results, formulating the gaps of knowledge still existing: such as the lack of reliable predictive biomarkers of tumour response to chemoradiation according to the molecular tumour subtype and lack of efficient radiosensitising agents specifically targeting bladder tumour cells. We provided guidance to improve forthcoming studies, such as taking into account molecular characteristics of the preclinical models and highlighted the value of using patient-derived xenografts as well as syngeneic models. Finally, this review could be a useful tool to set up new radiation-based combined treatments with an improved therapeutic index that is needed for bladder preservation.