The topography of mutational processes in breast cancer genomes

Evaluating topography of mutational signatures with SigProfilerTopography

The mutations found in a cancer genome are shaped by diverse processes, each displaying a characteristic mutational signature that may be influenced by the genome's architecture. While prior analyses have evaluated the effect of topographical genomic features on mutational signatures, there has been no computational tool that can comprehensively examine this interplay. Here, we present SigProfilerTopography, a Python package that allows evaluating the effect of chromatin organization, histone modifications, transcription factor binding, DNA replication, and DNA transcription on the activities of different mutational processes. SigProfilerTopography elucidates the unique topographical characteristics of mutational signatures, unveiling their underlying biological and molecular mechanisms.

Replicative DNA polymerase epsilon and delta holoenzymes show wide-ranging inhibition at G-quadruplexes in the human genome

G-quadruplexes (G4s) are functional elements of the human genome, some of which inhibit DNA replication. We investigated replication of G4s within highly abundant microsatellite (GGGA, GGGT) and transposable element (L1 and SVA) sequences. We found that genome-wide, numerous motifs are located preferentially on the replication leading strand and the transcribed strand templates. We directly tested replicative polymerase ϵ and δ holoenzyme inhibition at these G4s, compared to low abundant motifs. For all G4s, DNA synthesis inhibition was higher on the G-rich than C-rich strand or control sequence. No single G4 was an absolute block for either holoenzyme; however, the inhibitory potential varied over an order of magnitude. Biophysical analyses showed the motifs form varying topologies, but replicative polymerase inhibition did not correlate with a specific G4 structure. Addition of the G4 stabilizer pyridostatin severely inhibited forward polymerase synthesis specifically on the G-rich strand, enhancing G/C strand asynchrony. Our results reveal that replicative polymerase inhibition at every G4 examined is distinct, causing complementary strand synthesis to become asynchronous, which could contribute to slowed fork elongation. Altogether, we provide critical information regarding how replicative eukaryotic holoenzymes navigate synthesis through G4s naturally occurring thousands of times in functional regions of the human genome.

Defining APOBEC-induced mutation signatures and modifying activities in yeast

APOBEC cytidine deaminases guard cells in a variety of organisms from invading viruses and foreign nucleic acids. Recently, several human APOBECs have been implicated in mutating evolving cancer genomes. Expression of APOBEC3A and APOBEC3B in yeast allowed experimental derivation of the substitution patterns they cause in dividing cells, which provided critical links to these enzymes in the etiology of the COSMIC single base substitution (SBS) signatures 2 and 13 in human tumors. Additionally, the ability to scale yeast experiments to high-throughput screens allows use of this system to also investigate cellular pathways impacting the frequency of APOBEC-induced mutation. Here, we present validated methods utilizing yeast to determine APOBEC mutation signatures, genetic interactors, and chromosomal substrate preferences. These methods can be employed to assess the potential of other human APOBECs and APOBEC orthologs in different species to contribute to cancer genome evolution as well as define the pathways that protect the nuclear genome from inadvertent APOBEC activity during viral restriction.

Differentiation signals induce APOBEC3A expression via GRHL3 in squamous epithelia and squamous cell carcinoma

Two APOBEC DNA cytosine deaminase enzymes, APOBEC3A and APOBEC3B, generate somatic mutations in cancer, thereby driving tumour development and drug resistance. Here, we used single-cell RNA sequencing to study APOBEC3A and APOBEC3B expression in healthy and malignant mucosal epithelia, validating key observations with immunohistochemistry, spatial transcriptomics and functional experiments. Whereas APOBEC3B is expressed in keratinocytes entering mitosis, we show that APOBEC3A expression is confined largely to terminally differentiating cells and requires grainyhead-like transcription factor 3 (GRHL3). Thus, in normal tissue, neither deaminase appears to be expressed at high levels during DNA replication, the cell-cycle stage associated with APOBEC-mediated mutagenesis. In contrast, in squamous cell carcinoma we find that, there is expansion of GRHL3expression and activity to a subset of cells undergoing DNA replication and concomitant extension of APOBEC3A expression to proliferating cells. These findings suggest that APOBEC3A may play a functional role during keratinocyte differentiation, and offer a mechanism for acquisition of APOBEC3A mutagenic activity in tumours.

Homologous recombination deficiency (HRD) diagnostics: underlying mechanisms and new perspectives

Homologous recombination deficiency (HRD) is considered a universal and effective sign of a tumor's sensitivity to poly(ADP-ribose) polymerase (PARP) inhibitors. HRD diagnostics have undergone several stages of transformations: from detection of point mutations in HR-related genes and large regions with loss of heterozygosity detected using single-nucleotide polymorphism arrays to whole-genome signatures of single-nucleotide variants, large genomic rearrangements (LGRs), and copy number alterations. All these methods have their own advantages and limitations. HRD tests, based on signatures of LGRs and copy number alterations, show in hindsight that some progenitor cells have possessed HRD status but not the current state of the genome. The aim of this review was to compare different methods of HRD detection and mechanisms of formation of HRD-specific LGRs. In the last several years, new data appeared implying a crucial role of proteins BRCA1 and BRCA2 in the resolution of stalled replication forks that may be associated with at least some of LGRs observed in HRD-positive tumors. Reviewing current knowledge on these mechanisms, distributions of different LGR types, and limitations of sequencing technologies and algorithms of data analysis, we offer some new perspectives on HRD diagnostics. We hope that this review will help to accelerate the development of new diagnostic approaches in this important field of molecular oncology.

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

The cytidine deaminase APOBEC3C has unique sequence and genome feature preferences

APOBEC proteins are cytidine deaminases that restrict the replication of viruses and transposable elements. Several members of the APOBEC3 family, APOBEC3A, APOBEC3B, and APOBEC3H-I, can access the nucleus and cause what is thought to be indiscriminate deamination of the genome, resulting in mutagenesis and genome instability. Although APOBEC3C is also present in the nucleus, the full scope of its deamination target preferences is unknown. By expressing human APOBEC3C in a yeast model system, I have defined the APOBEC3C mutation signature, as well as the preferred genome features of APOBEC3C targets. The APOBEC3C mutation signature is distinct from those of the known cancer genome mutators APOBEC3A and APOBEC3B. APOBEC3C produces DNA strand-coordinated mutation clusters, and APOBEC3C mutations are enriched near the transcription start sites of active genes. Surprisingly, APOBEC3C lacks the bias for the lagging strand of DNA replication that is seen for APOBEC3A and APOBEC3B. The unique preferences of APOBEC3C constitute a mutation profile that will be useful in defining sites of APOBEC3C mutagenesis in human genomes.

Clinical Pharmacist-Led Interventions for Improving Breast Cancer Management-A Scoping Review

Breast cancer is the leading cause of cancer-related death in women worldwide and the fifth most common cause of cancer death overall. Most women with breast cancer have a good prognosis if the cancer is detected at an early stage and the patients have access to the appropriate treatment and disease management. This study aims to evaluate the impact of pharmacist-led interventions on breast cancer management and health outcomes. A literature review was carried out through the scientific databases PubMed, Scopus, and Web of Science using predefined keywords. Only full-text original articles written in English that investigated the role of the pharmacist in the management of breast cancer were included in the final analysis. No publication date limits were set. A total of 1625 articles were retrieved from the electronic databases, of which 14 met the inclusion criteria. The current scoping review consists of different study types, including randomized controlled trials, cross-sectional studies, pre-post studies, retrospective cohort studies, quality improvement projects, case-control studies, and one pharmacoeconomic study. Pharmacists commonly provided the following interventions: consultations regarding chemotherapy treatment, risk assessment and patient education, adverse drug reactions and drug-drug interactions detection, and adherence assessment. This scoping review highlights the beneficial effects of the involvement of pharmacists in breast cancer management, such as better quality of life, reduced drug interaction risk, greater adherence rates, and improved patient knowledge. This confirms the importance of including the pharmacist in the oncology team caring for patients with breast cancer.

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.

Selective pressures of platinum compounds shape the evolution of therapy-related myeloid neoplasms

Therapy-related myeloid neoplasms (t-MN) arise as a complication of chemo- and/or radiotherapy. Although t-MN can occur both in adult and childhood cancer survivors, the mechanisms driving therapy-related leukemogenesis likely vary across different ages. Chemotherapy is thought to induce driver mutations in children, whereas in adults pre-existing mutant clones are selected by the exposure. However, selective pressures induced by chemotherapy early in life are less well studied. Here, we use single-cell whole genome sequencing and phylogenetic inference to show that the founding cell of t-MN in children starts expanding after cessation of platinum exposure. In patients with Li-Fraumeni syndrome, characterized by a germline TP53 mutation, we find that the t-MN already expands during treatment, suggesting that platinum-induced growth inhibition is TP53-dependent. Our results demonstrate that germline aberrations can interact with treatment exposures in inducing t-MN, which is important for the development of more targeted, patient-specific treatment regimens and follow-up.

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.

A Novel Germline Mutation of BRCA1 and Integrated Analysis With Somatic Mutation in a Chinese Multi-Cancer Family

The presence of mutations in the BRCA1 gene (MIM: 113705) is widely recognized as a significant genetic predisposition for ovarian cancer. This study investigated the genomic mutations in a Chinese family with a history of ovarian, breast, and rectal adenocarcinoma. A novel germline mutation (Phe1695Val) in BRCA1 was identified through whole-exome sequencing. Subsequently, we performed whole-genome sequencing to identify somatic mutations and analyze mutational signatures in individuals carrying the novel germline mutation. Our findings revealed a correlation between somatic mutational signatures and the BRCA1 germline mutation in the proband with ovarian cancer, while no such association was observed in the tumor tissue from the patient with breast cancer. Furthermore, distinct somatic driver mutations were identified, a truncated mutation in the TP53 gene in the ovarian tumor tissue, and a hotspot mutation in the PIK3CA gene in the breast cancer. According to our findings, the BRCA1 F1695V mutation is linked to ovarian cancer susceptibility in the family and causes specific somatic mutational profiles.

Strand-resolved mutagenicity of DNA damage and repair

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

Chromosome evolution screens recapitulate tissue-specific tumor aneuploidy patterns

Whole chromosome and arm-level copy number alterations occur at high frequencies in tumors, but their selective advantages, if any, are poorly understood. Here, utilizing unbiased whole chromosome genetic screens combined with in vitro evolution to generate arm- and subarm-level events, we iteratively selected the fittest karyotypes from aneuploidized human renal and mammary epithelial cells. Proliferation-based karyotype selection in these epithelial lines modeled tissue-specific tumor aneuploidy patterns in patient cohorts in the absence of driver mutations. Hi-C-based translocation mapping revealed that arm-level events usually emerged in multiples of two via centromeric translocations and occurred more frequently in tetraploids than diploids, contributing to the increased diversity in evolving tetraploid populations. Isogenic clonal lineages enabled elucidation of pro-tumorigenic mechanisms associated with common copy number alterations, revealing Notch signaling potentiation as a driver of 1q gain in breast cancer. We propose that intrinsic, tissue-specific proliferative effects underlie tumor copy number patterns in cancer.

Molecular analysis of cancer genomes in children with Lynch syndrome: Exploring causal associations

Lynch syndrome (LS) predisposes to cancer in adulthood and is caused by heterozygous germline variants in a mismatch repair (MMR) gene. Recent studies show an increased prevalence of LS among children with cancer, suggesting a causal relationship. For LS-spectrum (LSS) cancers, including high-grade gliomas and colorectal cancer, causality has been supported by typical MMR-related tumor characteristics, but for non-LSS cancers, causality is unclear. We characterized 20 malignant tumors of 18 children with LS, including 16 non-LSS tumors. We investigated second hits, tumor mutational load, mutational signatures and MMR protein expression. In all LSS tumors and three non-LSS tumors, we detected MMR deficiency caused by second hit somatic alterations. Furthermore, these MMR-deficient tumors carried driver variants that likely originated as a consequence of MMR deficiency. However, in 13 non-LSS tumors (81%), a second hit and MMR deficiency were absent, thus a causal link between LS and cancer development in these children is lacking. These findings demonstrate that causality of LS in children with cancer, which can be determined by molecular tumor characterization, seems to be restricted to specific tumor types. Large molecular and epidemiological studies are needed to further refine the tumor spectrum in children with LS.

Differentiation signals induce APOBEC3A expression via GRHL3 in squamous epithelia and squamous cell carcinoma

Two APOBEC (apolipoprotein-B mRNA editing enzyme catalytic polypeptide-like) DNA cytosine deaminase enzymes (APOBEC3A and APOBEC3B) generate somatic mutations in cancer, driving tumour development and drug resistance. Here we used single cell RNA sequencing to study APOBEC3A and APOBEC3B expression in healthy and malignant mucosal epithelia, validating key observations with immunohistochemistry, spatial transcriptomics and functional experiments. Whereas APOBEC3B is expressed in keratinocytes entering mitosis, we show that APOBEC3A expression is confined largely to terminally differentiating cells and requires Grainyhead-like transcription factor 3 (GRHL3). Thus, in normal tissue, neither deaminase appears to be expressed at high levels during DNA replication, the cell cycle stage associated with APOBEC-mediated mutagenesis. In contrast, we show that in squamous cell carcinoma tissues, there is expansion of GRHL3 expression and activity to a subset of cells undergoing DNA replication and concomitant extension of APOBEC3A expression to proliferating cells. These findings indicate a mechanism for acquisition of APOBEC3A mutagenic activity in tumours.

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.

Evaluating topography of mutational signatures with SigProfilerTopography

The mutations found in a cancer genome are shaped by diverse processes, each displaying a characteristic mutational signature that may be influenced by the genome's architecture. While prior analyses have evaluated the effect of topographical genomic features on mutational signatures, there has been no computational tool that can comprehensively examine this interplay. Here, we present SigProfilerTopography, a Python package that allows evaluating the effect of chromatin organization, histone modifications, transcription factor binding, DNA replication, and DNA transcription on the activities of different mutational processes. SigProfilerTopography elucidates the unique topographical characteristics of mutational signatures, unveiling their underlying biological and molecular mechanisms.

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.

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

DNA mismatch repair protects the genome from oxygen-induced replicative mutagenesis

DNA mismatch repair (MMR) corrects mismatched DNA bases arising from multiple sources including polymerase errors and base damage. By detecting spontaneous mutagenesis using whole genome sequencing of cultured MMR deficient human cell lines, we show that a primary role of MMR is the repair of oxygen-induced mismatches. We found an approximately twofold higher mutation rate in MSH6 deficient DLD-1 cells or MHL1 deficient HCT116 cells exposed to atmospheric conditions as opposed to mild hypoxia, which correlated with oxidant levels measured using electron paramagnetic resonance spectroscopy. The oxygen-induced mutations were dominated by T to C base substitutions and single T deletions found primarily on the lagging strand. A broad sequence context preference, dependence on replication timing and a lack of transcriptional strand bias further suggested that oxygen-induced mutations arise from polymerase errors rather than oxidative base damage. We defined separate low and high oxygen-specific MMR deficiency mutation signatures common to the two cell lines and showed that the effect of oxygen is observable in MMR deficient cancer genomes, where it best correlates with the contribution of mutation signature SBS21. Our results imply that MMR corrects oxygen-induced genomic mismatches introduced by a replicative process in proliferating cells.

Characterization of highly active mutational signatures in tumors from a large Chinese population

The majority of mutational signatures have been characterized in tumors from Western countries and the degree to which mutational signatures are similar or different in Eastern populations has not been fully explored. We leveraged a large-scale clinical sequencing cohort of tumors from a Chinese population containing 25 tumor types and found that the highly active mutational signatures were similar to those previously characterized1,2. The aristolochic acid signature SBS22 was observed in four soft tissue sarcomas and the POLE-associated signature SBS10 was observed in a gallbladder carcinoma. In lung adenocarcinoma, the polycyclic aromatic hydrocarbon (PAH) signature SBS4 was significantly higher in males compared to females but not associated with smoking status. The UV-associated signature SBS7 was significantly lower in cutaneous melanomas from the Chinese population compared to a similar American cohort. Overall, these results add to our understanding of the mutational processes that contribute to tumors from the Chinese population.

Clinical utility of genomic signatures for the management of early and metastatic triple-negative breast cancer

PURPOSE OF REVIEW

This comprehensive review aims to provide timely and relevant insights into the current therapeutic landscape for triple-negative breast cancer (TNBC) and the molecular features underlying this subtype. It emphasizes the need for more reliable biomarkers to refine prognostication and optimize therapy, considering the aggressive nature of TNBC and its limited targeted treatment options.

RECENT FINDINGS

The review explores the multidisciplinary management of early TNBC, which typically involves systemic chemotherapy, surgery, and radiotherapy. It highlights the emergence of immune checkpoint inhibitors (ICIs), poly(ADP-ribose) polymerase (PARP) inhibitors, and antibody-drug conjugates (ADCs) as promising therapeutic strategies for TNBC. Recent clinical trials investigating the use of ICIs in combination with chemotherapy and the approval of pembrolizumab and atezolizumab for PD-L1-positive metastatic TNBC are discussed. The efficacy of PARP inhibitors and ADCs in treating TNBC patients with specific genetic alterations is also highlighted.

SUMMARY

The findings discussed in this review have significant implications for clinical practice and research in TNBC. The identification of distinct molecular subtypes through gene expression profiling has enabled a better understanding of TNBC heterogeneity and its clinical implications. This knowledge has the potential to guide treatment decisions, as different subtypes display varying responses to neoadjuvant chemotherapy. Furthermore, the review emphasizes the importance of developing reliable genomic and transcriptomic signatures as biomarkers to refine patient prognostication and optimize therapy selection in TNBC. Integrating these signatures into clinical practice may lead to more personalized treatment approaches, improving outcomes for TNBC patients.

Intra-Patient Genomic Variations of Human Papillomavirus Type 31 in Cervical Cancer and Precancer

Human papillomavirus type 31 (HPV31) is detected less frequently in cervical cancer than two major causative types, HPV16 and HPV18. Here, we report a comprehensive analysis of HPV31 genome sequences in cervical lesions collected from Japanese women. Of 52 HPV31-positive cervical specimens analyzed by deep sequencing, 43 samples yielded complete genome sequences of around 7900 base pairs and 9 samples yielded partially deleted genome sequences. Phylogenetic analysis showed that HPV31 variant distribution was lineage A in 19 samples (36.5%), lineage B in 28 samples (53.8%), and lineage C in 5 samples (9.6%), indicating that lineage B variants are dominant among HPV31 infections in Japan. Deletions in the viral genome were found in the region from the E1 to L1 genes, but all the deleted genomes retained the E6/E7 genes. Among intra-patient nucleotide variations relative to a consensus genome sequence in each sample, C-to-T substitutions were most frequently detected, followed by T-to-C and C-to-A substitutions. High-frequency, intra-patient mutations (>10%) in cervical cancer samples were found in the E1, E2, and E7 genes, and all of them were nonsynonymous substitutions. The enrichment of high-frequency nonsynonymous substitutions strongly suggests that these intra-patient mutations are positively selected during the development of cervical cancer/precancer.

Predicting regional somatic mutation rates using DNA motifs

How the locus-specificity of epigenetic modifications is regulated remains an unanswered question. A contributing mechanism is that epigenetic enzymes are recruited to specific loci by DNA binding factors recognizing particular sequence motifs (referred to as epi-motifs). Using these motifs to predict biological outputs depending on local epigenetic state such as somatic mutation rates would confirm their functionality. Here, we used DNA motifs including known TF motifs and epi-motifs as a surrogate of epigenetic signals to predict somatic mutation rates in 13 cancers at an average 23kbp resolution. We implemented an interpretable neural network model, called contextual regression, to successfully learn the universal relationship between mutations and DNA motifs, and uncovered motifs that are most impactful on the regional mutation rates such as TP53 and epi-motifs associated with H3K9me3. Furthermore, we identified genomic regions with significantly higher mutation rates than the expected values in each individual tumor and demonstrated that such cancer-related regions can accurately predict cancer types. Interestingly, we found that the same mutation signatures often have different contributions to cancer-related and cancer-independent regions, and we also identified the motifs with the most contribution to each mutation signature.

APOBEC3B regulates R-loops and promotes transcription-associated mutagenesis in cancer

The single-stranded DNA cytosine-to-uracil deaminase APOBEC3B is an antiviral protein implicated in cancer. However, its substrates in cells are not fully delineated. Here APOBEC3B proteomics reveal interactions with a surprising number of R-loop factors. Biochemical experiments show APOBEC3B binding to R-loops in cells and in vitro. Genetic experiments demonstrate R-loop increases in cells lacking APOBEC3B and decreases in cells overexpressing APOBEC3B. Genome-wide analyses show major changes in the overall landscape of physiological and stimulus-induced R-loops with thousands of differentially altered regions, as well as binding of APOBEC3B to many of these sites. APOBEC3 mutagenesis impacts genes overexpressed in tumors and splice factor mutant tumors preferentially, and APOBEC3-attributed kataegis are enriched in RTCW motifs consistent with APOBEC3B deamination. Taken together with the fact that APOBEC3B binds single-stranded DNA and RNA and preferentially deaminates DNA, these results support a mechanism in which APOBEC3B regulates R-loops and contributes to R-loop mutagenesis in cancer.

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.

Topography of mutational signatures in human cancer

The somatic mutations found in a cancer genome are imprinted by different mutational processes. Each process exhibits a characteristic mutational signature, which can be affected by the genome architecture. However, the interplay between mutational signatures and topographical genomic features has not been extensively explored. Here, we integrate mutations from 5,120 whole-genome-sequenced tumors from 40 cancer types with 516 topographical features from ENCODE to evaluate the effect of nucleosome occupancy, histone modifications, CTCF binding, replication timing, and transcription/replication strand asymmetries on the cancer-specific accumulation of mutations from distinct mutagenic processes. Most mutational signatures are affected by topographical features, with signatures of related etiologies being similarly affected. Certain signatures exhibit periodic behaviors or cancer-type-specific enrichments/depletions near topographical features, revealing further information about the processes that imprinted them. Our findings, disseminated via the COSMIC (Catalog of Somatic Mutations in Cancer) signatures database, provide a comprehensive online resource for exploring the interactions between mutational signatures and topographical features across human cancer.

Spontaneous mutagenesis in human cells is controlled by REV1-Polymerase ζ and PRIMPOL

Translesion DNA synthesis (TLS) facilitates replication over damaged or difficult-to-replicate templates by employing specialized DNA polymerases. We investigate the effect on spontaneous mutagenesis of three main TLS control mechanisms: REV1 and PCNA ubiquitylation that recruit TLS polymerases and PRIMPOL that creates post-replicative gaps. Using whole-genome sequencing of cultured human RPE-1 cell clones, we find that REV1 and Polymerase ζ are wholly responsible for one component of base substitution mutagenesis that resembles homologous recombination deficiency, whereas the remaining component that approximates oxidative mutagenesis is reduced in PRIMPOL-/- cells. Small deletions in short repeats appear in REV1-/-PCNAK164R/K164R double mutants, revealing an alternative TLS mechanism. Also, 500-5,000 bp deletions appear in REV1-/- and REV3L-/- mutants, and chromosomal instability is detectable in REV1-/-PRIMPOL-/- cells. Our results indicate that TLS protects the genome from deletions and large rearrangements at the expense of being responsible for the majority of spontaneous base substitutions.

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.

Implications of Translesion DNA Synthesis Polymerases on Genomic Stability and Human Health

Replication fork arrest-induced DNA double strand breaks (DSBs) caused by lesions are effectively suppressed in cells due to the presence of a specialized mechanism, commonly referred to as DNA damage tolerance (DDT). In eukaryotic cells, DDT is facilitated through translesion DNA synthesis (TLS) carried out by a set of DNA polymerases known as TLS polymerases. Another parallel mechanism, referred to as homology-directed DDT, is error-free and involves either template switching or fork reversal. The significance of the DDT pathway is well established. Several diseases have been attributed to defects in the TLS pathway, caused either by mutations in the TLS polymerase genes or dysregulation. In the event of a replication fork encountering a DNA lesion, cells switch from high-fidelity replicative polymerases to low-fidelity TLS polymerases, which are associated with genomic instability linked with several human diseases including, cancer. The role of TLS polymerases in chemoresistance has been recognized in recent years. In addition to their roles in the DDT pathway, understanding noncanonical functions of TLS polymerases is also a key to unraveling their importance in maintaining genomic stability. Here we summarize the current understanding of TLS pathway in DDT and its implication for human health.

A mutation-level covariate model for mutational signatures

Mutational processes and their exposures in particular genomes are key to our understanding of how these genomes are shaped. However, current analyses assume that these processes are uniformly active across the genome without accounting for potential covariates such as strand or genomic region that could impact such activities. Here we suggest the first mutation-covariate models that explicitly model the effect of different covariates on the exposures of mutational processes. We apply these models to test the impact of replication strand on these processes and compare them to strand-oblivious models across a range of data sets. Our models capture replication strand specificity, point to signatures affected by it, and score better on held-out data compared to standard models that do not account for mutation-level covariate information.

Progesterone from ovulatory menstrual cycles is an important cause of breast cancer

Many factors, including reproductive hormones, have been linked to a woman's risk of developing breast cancer (BC). We reviewed the literature regarding the relationship between ovulatory menstrual cycles (MCs) and BC risk. Physiological variations in the frequency of MCs and interference with MCs through genetic variations, pathological conditions and or pharmaceutical interventions revealed a strong link between BC risk and the lifetime number of MCs. A substantial reduction in BC risk is observed in situations without MCs. In genetic or transgender situations with normal female breasts and estrogens, but no progesterone (P4), the incidence of BC is very low, suggesting an essential role of P4. During the MC, P4 has a strong proliferative effect on normal breast epithelium, whereas estradiol (E2) has only a minimal effect. The origin of BC has been strongly linked to proliferation associated DNA replication errors, and the repeated stimulation of the breast epithelium by P4 with each MC is likely to impact the epithelial mutational burden. Long-lived cells, such as stem cells, present in the breast epithelium, can carry mutations forward for an extended period of time, and studies show that breast tumors tend to take decades to develop before detection. We therefore postulate that P4 is an important factor in a woman's lifetime risk of developing BC, and that breast tumors arising during hormonal contraception or after menopause, with or without menopausal hormone therapy, are the consequence of the outgrowth of pre-existing neoplastic lesions, eventually stimulated by estrogens and some progestins.

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.

Clinical and genomic characterization of mutational signatures across human cancers

Mutational signatures, the generic patterns of mutations, are the footprints of both endogenous and exogenous factors that have influenced cancer development. To date, dozens of mutational signatures have been discerned through computational methods. However, the etiology, mutational properties, clonality, immunology and prognostic value of mutation signatures across cancer types are poorly understood. To address this, we extensively characterized mutational signatures across 8836 cancer samples spanning 42 cancer types. We confirmed and extended clinical and genomic features associated with mutation signatures. Mutation distribution analysis showed that most mutation processes were depleted in exons and APOBEC signatures (SBS2 and SBS13), the Pol-η related signature (SBS9) and SBS40 tended to contribute clustered mutations. We observed that age-related signatures (SBS1 and SBS5) and SBS40 tended to induce mutations affecting cancer genes and subclonal drivers posted by specific signatures (eg, mismatch repair deficiency-related signature SBS44) were unlikely subjected to positive selection. We also revealed early mutation signatures (eg, UV light exposure-related signature SBS7a) and signatures (eg, reactive oxygen species-related signature SBS18) predominated in the late stage of tumorigenesis. Comprehensive association analysis of mutation processes with microenvironment revealed that APOBEC- and mismatch repair deficiency-related signatures were positively associated with immune parameters, while age-related signatures showed negative correlations. In addition, prognostic association analysis showed that many signatures were favorable (eg, SBS9) or adverse factors (eg, SBS18) of patient survival. Our findings enhance appreciation of the role of mutational signatures in tumor evolution and underline their potential in immunotherapy guidance and prognostic prediction.

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.

Roles of trans-lesion synthesis (TLS) DNA polymerases in tumorigenesis and cancer therapy

DNA damage tolerance and mutagenesis are hallmarks and enabling characteristics of neoplastic cells that drive tumorigenesis and allow cancer cells to resist therapy. The 'Y-family' trans-lesion synthesis (TLS) DNA polymerases enable cells to replicate damaged genomes, thereby conferring DNA damage tolerance. Moreover, Y-family DNA polymerases are inherently error-prone and cause mutations. Therefore, TLS DNA polymerases are potential mediators of important tumorigenic phenotypes. The skin cancer-propensity syndrome xeroderma pigmentosum-variant (XPV) results from defects in the Y-family DNA Polymerase Pol eta (Polη) and compensatory deployment of alternative inappropriate DNA polymerases. However, the extent to which dysregulated TLS contributes to the underlying etiology of other human cancers is unclear. Here we consider the broad impact of TLS polymerases on tumorigenesis and cancer therapy. We survey the ways in which TLS DNA polymerases are pathologically altered in cancer. We summarize evidence that TLS polymerases shape cancer genomes, and review studies implicating dysregulated TLS as a driver of carcinogenesis. Because many cancer treatment regimens comprise DNA-damaging agents, pharmacological inhibition of TLS is an attractive strategy for sensitizing tumors to genotoxic therapies. Therefore, we discuss the pharmacological tractability of the TLS pathway and summarize recent progress on development of TLS inhibitors for therapeutic purposes.

Modeling the relationship between gene expression and mutational signature

Background

Mutational signatures computed from somatic mutations, allow an in-depth understanding of tumorigenesis and may illuminate early prevention strategies. Many studies have shown the regulation effects between somatic mutation and gene expression dysregulation.

Methods

We hypothesized that there are potential associations between mutational signature and gene expression. We capitalized upon RNA-seq data to model 49 established mutational signatures in 33 cancer types. Both accuracy and area under the curve were used as performance measures in five-fold cross-validation.

Results

A total of 475 models using unconstrained genes, and 112 models using protein-coding genes were selected for future inference purposes. An independent gene expression dataset on lung cancer smoking status was used for validation which achieved over 80% for both accuracy and area under the curve.

Conclusion

These results demonstrate that the associations between gene expression and somatic mutations can translate into the associations between gene expression and mutational signatures.

Clinical Utility of Genomic Tests Evaluating Homologous Recombination Repair Deficiency (HRD) for Treatment Decisions in Early and Metastatic Breast Cancer

Breast cancer is the most frequently occurring cancer worldwide. With its increasing incidence, it is a major public health problem, with many therapeutic challenges such as precision medicine for personalized treatment. Thanks to next-generation sequencing (NGS), progress in biomedical technologies, and the use of bioinformatics, it is now possible to identify specific molecular alterations in tumor cells-such as homologous recombination deficiencies (HRD)-enabling us to consider using DNA-damaging agents such as platinum salts or PARP inhibitors. Different approaches currently exist to analyze impairment of the homologous recombination pathway, e.g., the search for specific mutations in homologous recombination repair (HRR) genes, such as BRCA1/2; the use of genomic scars or mutational signatures; or the development of functional tests. Nevertheless, the role and value of these different tests in breast cancer treatment decisions remains to be clarified. In this review, we summarize current knowledge on the clinical utility of genomic tests, evaluating HRR deficiency for treatment decisions in early and metastatic breast cancer.

Proteogenomics decodes the evolution of human ipsilateral breast cancer

Ipsilateral breast tumor recurrence (IBTR) is a clinically important event, where an isolated in-breast recurrence is a potentially curable event but associated with an increased risk of distant metastasis and breast cancer death. It remains unclear if IBTRs are associated with molecular changes that can be explored as a resource for precision medicine strategies. Here, we employed proteogenomics to analyze a cohort of 27 primary breast cancers and their matched IBTRs to define proteogenomic determinants of molecular tumor evolution. Our analyses revealed a relationship between hormonal receptors status and proliferation levels resulting in the gain of somatic mutations and copy number. This in turn re-programmed the transcriptome and proteome towards a highly replicating and genomically unstable IBTRs, possibly enhanced by APOBEC3B. In order to investigate the origins of IBTRs, a second analysis that included primaries with no recurrence pinpointed proliferation and immune infiltration as predictive of IBTR. In conclusion, our study shows that breast tumors evolve into different IBTRs depending on hormonal status and proliferation and that immune cell infiltration and Ki-67 are significantly elevated in primary tumors that develop IBTR. These results can serve as a starting point to explore markers to predict IBTR formation and stratify patients for adjuvant therapy.

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.

Circulating miRNA Expression Profiling in Breast Cancer Molecular Subtypes: Applying Machine Learning Analysis in Bioinformatics

BACKGROUND/AIM

Breast cancer is a leading worldwide cause of female cancer-related morbidity and mortality. Since molecular characteristics increasingly guide disease management, demystifying breast tumor miRNA signature emerges as an essential step toward personalized care. This study aimed to investigate the variations in circulating miRNA expression profiles between breast cancer subtypes and healthy controls and to identify relevant target genes and molecular functions.

MATERIALS AND METHODS

MiRNA expression was tested by miScript™ miRNA PCR Array Human Cancer Pathway Finder kit, and subsequently, a machine learning approach was applied for miRNA profiling of the various breast cancer molecular subtypes.

RESULTS

Serum samples from patients with primary breast cancer (n=66) and healthy controls (n=16) were analyzed. MiR-21 was the single common molecule among all breast cancer subtypes. Furthermore, several miRNAs were found to be differentially expressed explicitly in the different subtypes; luminal A (miR-23b, miR-142, miR-29a, miR-181d, miR-16, miR-29b, miR-155, miR-181c), luminal B (miR-148a, let-7d, miR-92a, miR-34c, let-7b, miR-15a), HER2+ (miR-125b, miR-134, miR-98, miR-143, miR-138, miR-135b) and triple negative breast cancer (miR-17, miR-150, miR-210, miR-372, let-7f, miR-191, miR-133b, miR-146b, miR-7). Finally, miRNA-associated target genes and molecular functions were identified.

CONCLUSION

Applying a machine learning approach to delineate miRNA signatures of various breast cancer molecular subtypes allows further understanding of molecular disease characteristics that can prove clinically relevant.

Elevated Mutational Age in Blood of Children Treated for Cancer Contributes to Therapy-Related Myeloid Neoplasms

Childhood cancer survivors are confronted with various chronic health conditions like therapy-related malignancies. However, it is unclear how exposure to chemotherapy contributes to the mutation burden and clonal composition of healthy tissues early in life. Here, we studied mutation accumulation in hematopoietic stem and progenitor cells (HSPC) before and after cancer treatment of 24 children. Of these children, 19 developed therapy-related myeloid neoplasms (t-MN). Posttreatment HSPCs had an average mutation burden increase comparable to what treatment-naïve cells accumulate during 16 years of life, with excesses up to 80 years. In most children, these additional mutations were induced by clock-like processes, which are also active during healthy aging. Other patients harbored mutations that could be directly attributed to treatments like platinum-based drugs and thiopurines. Using phylogenetic inference, we demonstrate that most t-MN in children originate after the start of treatment and that leukemic clones become dominant during or directly after chemotherapy exposure.

SIGNIFICANCE

Our study shows that chemotherapy increases the mutation burden of normal blood cells in cancer survivors. Only few drugs damage the DNA directly, whereas in most patients, chemotherapy-induced mutations are caused by processes similar to those present during normal aging. This article is highlighted in the In This Issue feature, p. 1825.

Structural variants shape driver combinations and outcomes in pediatric high-grade glioma

We analyzed the contributions of structural variants (SVs) to gliomagenesis across 179 pediatric high-grade gliomas (pHGGs). The most recurrent SVs targeted MYC isoforms and receptor tyrosine kinases (RTKs), including an SV amplifying a MYC enhancer in 12% of diffuse midline gliomas (DMG), indicating an underappreciated role for MYC in pHGG. SV signature analysis revealed that tumors with simple signatures were TP53 wild type (TP53WT) but showed alterations in TP53 pathway members PPM1D and MDM4. Complex signatures were associated with direct aberrations in TP53, CDKN2A and RB1 early in tumor evolution and with later-occurring extrachromosomal amplicons. All pHGGs exhibited at least one simple-SV signature, but complex-SV signatures were primarily restricted to subsets of H3.3K27M DMGs and hemispheric pHGGs. Importantly, DMGs with complex-SV signatures were associated with shorter overall survival independent of histone mutation and TP53 status. These data provide insight into the impact of SVs on gliomagenesis and the mechanisms that shape them.

The shaping of cancer genomes with the regional impact of mutation processes

Mutation signature analysis has been used to infer the contributions of various DNA mutagenic-repair events in individual cancer genomes. Here, we build a statistical framework using a multinomial distribution to assign individual mutations to their cognate mutation signatures. We applied it to 47 million somatic mutations in 1925 publicly available cancer genomes to obtain a mutation signature map at the resolution of individual somatic mutations. Based on mutation signature-level genetic-epigenetic correlative analyses, mutations with transcriptional and replicative strand asymmetries show different enrichment patterns across genomes, and “transcribed” chromatin states and gene boundaries are particularly vulnerable to transcription-coupled repair activities. While causative processes of cancer-driving mutations can be diverse, as shown for converging effects of multiple mutational processes on TP53 mutations, the substantial fraction of recurrently mutated amino acids points to specific mutational processes, e.g., age-related C-to-T transition for KRAS p.G12 mutations. Our investigation of evolutionary trajectories with respect to mutation signatures further revealed that candidate pairs of early- vs. late-operative mutation processes in cancer genomes represent evolutionary dynamics of multiple mutational processes in the shaping of cancer genomes. We also observed that the local mutation clusters of kataegis often include mutations arising from multiple mutational processes, suggestive of a locally synchronous impact of multiple mutational processes on cancer genomes. Taken together, our examination of the genome-wide landscape of mutation signatures at the resolution of individual somatic mutations shows the spatially and temporally distinct mutagenesis-repair-replication histories of various mutational processes and their effects on shaping cancer genomes.

A statistical model that assigns non-hereditary DNA alterations known as somatic mutations to mutation “signatures” (groups of mutations arising from a specific biological process) on cancer genomes provides novel insights into disease evolution. Somatic mutations result from exposure to factors often linked to cancer development, such as tobacco or ultraviolet radiation. However, assigning a somatic mutation to a particular mutation “signature” remains challenging. The model created by Ruibin Xi (Peking University, China) and Tae-Min Kim (Catholic University of Korea, Seoul, South Korea) and co-workers grouped 47 million somatic mutations in 1925 cancer genomes into localized clusters before connecting them with mutation signatures. This strategy highlights the spatial and temporal patterns related to the origins of mutations, how the DNA strands are repaired and replicated, and how this influences the emerging cancer genome.

DNA nicks induce mutational signatures associated with BRCA1 deficiency

Analysis of human cancer genome sequences has revealed specific mutational signatures associated with BRCA1-deficient tumors, but the underlying mechanisms remain poorly understood. Here, we show that one-ended DNA double strand breaks (DSBs) converted from CRISPR/Cas9-induced nicks by DNA replication, not two-ended DSBs, cause more characteristic chromosomal aberrations and micronuclei in Brca1-deficient cells than in wild-type cells. BRCA1 is required for efficient homologous recombination of these nick-converted DSBs and suppresses bias towards long tract gene conversion and tandem duplication (TD) mediated by two-round strand invasion in a replication strand asymmetry. However, aberrant repair of these nick-converted one-ended DSBs, not that of two-ended DSBs in Brca1-deficient cells, generates mutational signatures such as small indels with microhomology (MH) at the junctions, translocations and small MH-mediated TDs, resembling those in BRCA1-deficient tumors. These results suggest a major contribution of DNA nicks to mutational signatures associated with BRCA1 deficiency in cancer and the underlying mechanisms.

KMT2C-deficient tumors have elevated APOBEC mutagenesis and genomic instability in multiple cancers

The histone methyltransferase KMT2C is among the most frequently mutated epigenetic modifier genes in cancer and plays an essential role in MRE11-dependent DNA replication fork restart. However, the effects of KMT2C deficiency on genomic instability during tumorigenesis are unclear. Analyzing 9,663 tumors from 30 cancer cohorts, we report that KMT2C mutant tumors have a significant excess of APOBEC mutational signatures in several cancer types. We show that KMT2C deficiency promotes APOBEC expression and deaminase activity, and compromises DNA replication speed and delays fork restart, facilitating APOBEC mutagenesis targeting single stranded DNA near stalled forks. APOBEC-mediated mutations primarily accumulate during early replication and tend to cluster along the genome and also in 3D nuclear domains. Excessive APOBEC mutational signatures in KMT2C mutant tumors correlate with elevated genome maintenance defects and signatures of homologous recombination deficiency. We propose that KMT2C deficiency is a likely promoter of APOBEC mutagenesis, which fosters further genomic instability during tumor progression in multiple cancer types.

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.

Biology and Model Predictions of the Dynamics and Heterogeneity of Chromatin-Nuclear Lamina Interactions

Associations of chromatin with the nuclear lamina, at the nuclear periphery, help shape the genome in 3 dimensions. The genomic landscape of lamina-associated domains (LADs) is well characterized, but much remains unknown on the physical and mechanistic properties of chromatin conformation at the nuclear lamina. Computational models of chromatin folding at, and interactions with, a surface representing the nuclear lamina are emerging in attempts to characterize these properties and predict chromatin behavior at the lamina in health and disease. Here, we highlight the heterogeneous nature of the nuclear lamina and LADs, outline the main 3-dimensional chromatin structural modeling methods, review applications of modeling chromatin-lamina interactions and discuss biological insights inferred from these models in normal and disease states. Lastly, we address perspectives on future developments in modeling chromatin interactions with the nuclear lamina.