Cancer Mutational Processes Vary in Their Association with Replication Timing and Chromatin Accessibility

Strategies for engineering oncolytic viruses to enhance cancer immunotherapy

Non-small cell lung cancer (NSCLC) is the predominant form of lung cancer and is characterized by rapid metastasis and high mortality, presenting a challenge for early-stage treatment modalities. The heterogeneity of NSCLC's tumor microenvironment (TME) significantly influences the efficacy of anti-PD-1 immune checkpoint inhibitors (ICIs) therapy, leading to varied patient responses. This review characterized different strains of oncolytic viruses in NSCLC and the different gene edits in pre-existing oncolytic viruses. This study also aimed to provide strategies to enhance anti-PD-1 therapy in NSCLC by engineering oncolytic viruses (OVs). This study offers insights into the genomic adaptations necessary for OVs targeting NSCLC, identify genetic determinants of anti-PD-1 response variability, and propose genomic edits to bolster therapy effectiveness. The primary goal of this study is to present a theoretically designed OV with a detailed genomic framework capable of enhancing the response to anti-PD-1 therapy, thereby advancing the field of cancer immunotherapy.

Replication timing alterations are associated with mutation acquisition during breast and lung cancer evolution

During each cell cycle, the process of DNA replication timing is tightly regulated to ensure the accurate duplication of the genome. The extent and significance of alterations in this process during malignant transformation have not been extensively explored. Here, we assess the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find that 6%-18% of the cancer genome exhibits ART, with regions with a change from early to late replication displaying an increased mutation rate and distinct mutational signatures. Whereas regions changing from late to early replication contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART occurs relatively early during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

The mutagenic consequences of defective DNA repair

Multiple separate repair mechanisms safeguard the genome against various types of DNA damage, and their failure can increase the rate of spontaneous mutagenesis. The malfunction of distinct repair mechanisms leads to genomic instability through different mutagenic processes. For example, defective mismatch repair causes high base substitution rates and microsatellite instability, whereas homologous recombination deficiency is characteristically associated with deletions and chromosome instability. This review presents a comprehensive collection of all mutagenic phenotypes associated with the loss of each DNA repair mechanism, drawing on data from a variety of model organisms and mutagenesis assays, and placing greatest emphasis on systematic analyses of human cancer datasets. We describe the latest theories on the mechanism of each mutagenic process, often explained by reliance on an alternative repair pathway or the error-prone replication of unrepaired, damaged DNA. Aided by the concept of mutational signatures, the genomic phenotypes can be used in cancer diagnosis to identify defective DNA repair pathways.

Cancer mutational signatures identification in clinical assays using neural embedding-based representations

While mutational signatures provide a plethora of prognostic and therapeutic insights, their application in clinical-setting, targeted gene panels is extremely limited. We develop a mutational representation model (which learns and embeds specific mutation signature connections) that enables prediction of dominant signatures with only a few mutations. We predict the dominant signatures across more than 60,000 tumors with gene panels, delineating their landscape across different cancers. Dominant signature predictions in gene panels are of clinical importance. These included UV, tobacco, and apolipoprotein B mRNA editing enzyme, catalytic polypeptide (APOBEC) signatures that are associated with better survival, independently from mutational burden. Further analyses reveal gene and mutation associations with signatures, such as SBS5 with TP53 and APOBEC with FGFR3S249C. In a clinical use case, APOBEC signature is a robust and specific predictor for resistance to epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs). Our model provides an easy-to-use way to detect signatures in clinical setting assays with many possible clinical implications for an unprecedented number of cancer patients.

Cell cycle gene alterations associate with a redistribution of mutation risk across chromosomal domains in human cancers

Mutations in human cells exhibit increased burden in heterochromatic, late DNA replication time (RT) chromosomal domains, with variation in mutation rates between tissues mirroring variation in heterochromatin and RT. We observed that regional mutation risk further varies between individual tumors in a manner independent of cell type, identifying three signatures of domain-scale mutagenesis in >4,000 tumor genomes. The major signature reflects remodeling of heterochromatin and of the RT program domains seen across tumors, tissues and cultured cells, and is robustly linked with higher expression of cell proliferation genes. Regional mutagenesis is associated with loss of activity of the tumor-suppressor genes RB1 and TP53, consistent with their roles in cell cycle control, with distinct mutational patterns generated by the two genes. Loss of regional heterogeneity in mutagenesis is associated with deficiencies in various DNA repair pathways. These mutation risk redistribution processes modify the mutation supply towards important genes, diverting the course of somatic evolution.

Single-cell DNA replication dynamics in genomically unstable cancers

Dysregulated DNA replication is both a cause and a consequence of aneuploidy, yet the dynamics of DNA replication in aneuploid cell populations remains understudied. We developed a new method, PERT, for inferring cell-specific DNA replication states from single-cell whole genome sequencing, and investigated clone-specific DNA replication dynamics in >50,000 cells obtained from a collection of aneuploid and clonally heterogeneous cell lines, xenografts and primary cancer tissues. Clone replication timing (RT) profiles correlated with future copy number changes in serially passaged cell lines. Cell type was the strongest determinant of RT heterogeneity, while whole genome doubling and mutational process were associated with accumulation of late S-phase cells and weaker RT associations. Copy number changes affecting chromosome X had striking impact on RT, with loss of the inactive X allele shifting replication earlier, and loss of inactive Xq resulting in reactivation of Xp. Finally, analysis of time series xenografts illustrate how cell cycle distributions approximate clone proliferation, recapitulating expected relationships between proliferation and fitness in treatment-naive and chemotherapeutic contexts.

Both cell autonomous and non-autonomous processes modulate the association between replication timing and mutation rate

Cancer somatic mutations are the product of multiple mutational and repair processes, some of which are tightly associated with DNA replication. Mutation rates (MR) are known to be higher in late replication timing (RT) regions, but different processes can affect this association. Systematic analysis of the mutational landscape of 2787 tumors from 32 tumor types revealed that approximately one third of the tumor samples show weak association between replication timing and mutation rate. Further analyses revealed that those samples have unique mutational signatures and are enriched with mutations in genes involved in DNA replication, DNA repair and chromatin structure. Surprisingly, analysis of differentially expressed genes between weak and strong RT-MR association groups revealed that tumors with weak association are enriched with genes associated with cell-cell communication and the immune system, suggesting a non-autonomous response to DNA damage.

The landscape of somatic mutations in lymphoblastoid cell lines

Somatic mutations have important biological ramifications while exerting substantial rate, type, and genomic location heterogeneity. Yet, their sporadic occurrence makes them difficult to study at scale and across individuals. Lymphoblastoid cell lines (LCLs), a model system for human population and functional genomics, harbor large numbers of somatic mutations and have been extensively genotyped. By comparing 1,662 LCLs, we report that the mutational landscape of the genome varies across individuals in terms of the number of mutations, their genomic locations, and their spectra; this variation may itself be modulated by somatic trans-acting mutations. Mutations attributed to the translesion DNA polymerase η follow two different modes of formation, with one mode accounting for the hypermutability of the inactive X chromosome. Nonetheless, the distribution of mutations along the inactive X chromosome appears to follow an epigenetic memory of the active form.

Cell-type specificity of the human mutation landscape with respect to DNA replication dynamics

The patterns of genomic mutations are associated with various genomic features, most notably late replication timing, yet it remains contested which mutation types and signatures relate to DNA replication dynamics and to what extent. Here, we perform high-resolution comparisons of mutational landscapes between lymphoblastoid cell lines, chronic lymphocytic leukemia tumors, and three colon adenocarcinoma cell lines, including two with mismatch repair deficiency. Using cell-type-matched replication timing profiles, we demonstrate that mutation rates exhibit heterogeneous replication timing associations among cell types. This cell-type heterogeneity extends to the underlying mutational pathways, as mutational signatures show inconsistent replication timing bias between cell types. Moreover, replicative strand asymmetries exhibit similar cell-type specificity, albeit with different relationships to replication timing than mutation rates. Overall, we reveal an underappreciated complexity and cell-type specificity of mutational pathways and their relationship to replication timing.

Mutational signatures association with replication timing in normal cells reveals similarities and differences with matched cancer tissues

Mutational signatures' association with replication timing (RT) has been studied in cancer samples, but the RT distribution of somatic mutations in non-cancerous cells was only minimally explored. Here, we performed comprehensive analyses of mutational signatures in 2.9 million somatic mutations across multiple non-cancerous tissues, stratified by early and late RT regions. We found that many mutational processes are active mainly or solely in early RT, such as SBS16 in hepatocytes and SBS88 in the colon, or in late RT, such as SBS4 in lung and hepatocytes, and SBS18 across many tissues. The two ubiquitous signatures, SBS1 and SBS5, showed late and early bias, respectively, across multiple tissues and in mutations representing germ cells. We also performed a direct comparison with cancer samples in 4 matched tissue-cancer types. Unexpectedly, while for most signatures the RT bias was consistent in normal tissue and in cancer, we found that SBS1's late RT bias is lost in cancer.

Regional mutational signature activities in cancer genomes

Cancer genomes harbor a catalog of somatic mutations. The type and genomic context of these mutations depend on their causes and allow their attribution to particular mutational signatures. Previous work has shown that mutational signature activities change over the course of tumor development, but investigations of genomic region variability in mutational signatures have been limited. Here, we expand upon this work by constructing regional profiles of mutational signature activities over 2,203 whole genomes across 25 tumor types, using data aggregated by the Pan-Cancer Analysis of Whole Genomes (PCAWG) consortium. We present GenomeTrackSig as an extension to the TrackSig R package to construct regional signature profiles using optimal segmentation and the expectation-maximization (EM) algorithm. We find that 426 genomes from 20 tumor types display at least one change in mutational signature activities (changepoint), and 306 genomes contain at least one of 54 recurrent changepoints shared by seven or more genomes of the same tumor type. Five recurrent changepoint locations are shared by multiple tumor types. Within these regions, the particular signature changes are often consistent across samples of the same type and some, but not all, are characterized by signatures associated with subclonal expansion. The changepoints we found cannot strictly be explained by gene density, mutation density, or cell-of-origin chromatin state. We hypothesize that they reflect a confluence of factors including evolutionary timing of mutational processes, regional differences in somatic mutation rate, large-scale changes in chromatin state that may be tissue type-specific, and changes in chromatin accessibility during subclonal expansion. These results provide insight into the regional effects of DNA damage and repair processes, and may help us localize genomic and epigenomic changes that occur during cancer development.

The comparison of cancer gene mutation frequencies in Chinese and U.S. patient populations

Knowing the mutation frequency of cancer genes in China is crucial for reducing the global health burden. We integrate the tumor epidemiological statistics with cancer gene mutation rates identified in 11,948 cancer patients to determine their weighted proportions within a Chinese cancer patient cohort. TP53 (51.4%), LRP1B (13.4%), PIK3CA (11.6%), KRAS (11.1%), EGFR (10.6%), and APC (10.5%) are identified as the top mutated cancer genes in China. Additionally, 18 common cancer types from both China and U.S. cohorts are analyzed and classified into three patterns principally based upon TP53 mutation rates: TP53-Top, TP53-Plus, and Non-TP53. Next, corresponding similarities and prominent differences are identified upon comparing the mutational profiles from both cohorts. Finally, the potential population-specific and environmental risk factors underlying the disparities in cancer gene mutation rates between the U.S. and China are analyzed. Here, we show and compare the mutation rates of cancer genes in Chinese and U.S. population cohorts, for a better understanding of the associated etiological and epidemiological factors, which are important for cancer prevention and therapy.

KIF26B in the Prognosis and Immune Biomarking of Various Cancers: A Pan-Cancer Study

KIF26B has been identified as an oncogene in several tumors; however, its utility as a prognostic indicator for various cancers has not yet been comprehensively evaluated. Here, we first examined how KIF26B intervenes in thirty-three cancers within the TCGA database, including potential immunological functions, and how it affects the prognosis. Based on the open databases TCGA, TIMER2, GEPIA2, GTEx, CPTAC, and HPA, we found that, when compared with normal tissues, KIF26B is overexpressed in 22 tumor tissues. Following a survival analysis, a relationship between the expression of KIF26B and the prognosis of various cancers was observed. Among the genetic alterations assessed, mutations were the most frequent. On the contrary, high phosphorylation levels of S977 were detected in breast cancer, KIRC, LUAD, and UCEC. We also found positive or negative correlations between KIF26B and the immune infiltration of endothelial cells and cancer-associated fibroblast infiltration. This could imply that patients may benefit from immunotherapy. Finally, KEGG pathways and GO enrichment analyses were implemented to identify the molecular mechanisms of KIF26B. This study illustrates the function of KIF26B from a pan-cancer perspective and offers a new horizon for cancer prognostic and immunotherapeutic investigations.