SEARCH · PubMed Health

Results for “mutation bias”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2025

Nonhypermutator Cancers Access Driver Mutations Through Reversals in Germline Mutational Bias.

Cancer is an evolutionary disease driven by mutations in asexually reproducing somatic cells. In asexual microbes, bias reversals in the mutation spectrum can speed adaptation by increasing access to previously undersampled beneficial mutations. By analyzing tumors from 20 tissues, along with normal tissue and the germline, we demonstrate this effect in cancer. Nonhypermutated tumors reverse the germline mutation bias and have consistent spectra across tissues. These spectra changes carry the signature of hypoxia, and they facilitate positive selection in cancer genes. Hypermutated and nonhypermutated tumors thus acquire driver mutations differently: hypermutated tumors by higher mutation rates and nonhypermutated tumors by changing the mutation spectrum to reverse the germline mutation bias.

Neoplasms↗

PubMed2024

Non-hypermutator cancers access driver mutations through reversals in germline mutational bias.

Cancer is an evolutionary disease driven by mutations in asexually-reproducing somatic cells. In asexual microbes, bias reversals in the mutation spectrum can speed adaptation by increasing access to previously undersampled beneficial mutations. By analyzing tumors from 20 tissues, along with normal tissue and the germline, we demonstrate this effect in cancer. Non-hypermutated tumors reverse the germline mutation bias and have consistent spectra across tissues. These spectra changes carry the signature of hypoxia, and they facilitate positive selection in cancer genes. Hypermutated and non-hypermutated tumors thus acquire driver mutations differently: hypermutated tumors by higher mutation rates and non-hypermutated tumors by changing the mutation spectrum to reverse the germline mutation bias.

Journal Article↗

PubMed2017

Single-sperm sequencing reveals the accelerated mitochondrial mutation rate in male Daphnia pulex (Crustacea, Cladocera).

Mutation rate in the nuclear genome differs between sexes, with males contributing more mutations than females to their offspring. The male-biased mutation rates in the nuclear genome is most likely to be driven by a higher number of cell divisions in spermatogenesis than in oogenesis, generating more opportunities for DNA replication errors. However, it remains unknown whether male-biased mutation rates are present in mitochondrial DNA (mtDNA). Although mtDNA is maternally inherited and male mtDNA mutation typically does not contribute to genetic variation in offspring, male mtDNA mutations are critical for male reproductive health. In this study, we measured male mtDNA mutation rate using publicly available whole-genome sequences of single sperm of the freshwater microcrustacean Daphnia pulex Using a stringent mutation detection pipeline, we found that the male mtDNA mutation rate is 3.32 × 10-6 per site per generation. All the detected mutations are heteroplasmic base substitutions, with 57% of mutations converting G/C to A/T nucleotides. Consistent with the male-biased mutation in the nuclear genome, the male mtDNA mutation rate in D. pulex is approximately 20 times higher than the female rate per generation. We propose that the elevated mutation rate per generation in male mtDNA is consistent with an increased number of cell divisions during male gametogenesis.

Animals↗

PubMed2026

Mutation rate heterogeneity biases variant effect prediction and reveals genuine mutational robustness.

Variant effect predictors (VEPs) are widely used to interpret the functional consequences of human genetic variation. Because most methods rely on sequence conservation, they implicitly treat conservation as evidence of functional constraint. However, substitution patterns across a phylogeny reflect not only selection but also differences in underlying mutation rates. Here, we show that this creates a systematic confounding: most VEPs capture mutation rate variation and misinterpret it as variation in functional importance. Widely used conservation metrics exhibit a related bias; in particular, phyloP scores correlate strongly with mutation rate even at putatively neutral sites. Consequently, variants at low-mutation-rate sites tend to be predicted as more damaging, and variants at highly mutable sites as more tolerated, than warranted by their true functional impact. We also identify a distinct biological signal in experimental measurements of mutational effects on protein stability: amino acid substitutions that are more likely to arise are, on average, less destabilizing than rarer substitutions. This provides empirical support for mutational robustness in the context of protein stability. However, this relationship is insufficient to explain the mutation-rate dependence observed in current VEP outputs. Together, our findings show that mutation rate heterogeneity systematically biases current variant effect prediction frameworks, highlight the need to model mutation probabilities explicitly in future VEPs, and reveal a genuine biological signal of mutational robustness.

conservation scores↗

PubMed2025

Bistable Mutation-Selection Equilibria and Violations of Fisher's Theorem in Tetraploids: Insights from Nonlinear Dynamics.

Polyploidy and whole genome duplication (WGD) are widespread biological phenomena with substantial cellular, meiotic, and genetic effects. Despite their prevalence and significance across the tree of life, population genetics theory for polyploids is not well developed. The lack of theoretical models limits our understanding of polyploid evolution and restricts our ability to harness polyploidy for crop improvement amidst increasing environmental stress. To address this gap, we developed and analyzed deterministic models of mutation-selection balance for tetraploids under polysomic (autotetraploid) and disomic (allotetraploid) inheritance patterns and arbitrary dominance relationships. We also introduced a new mathematical framework based on ordinary differential equations and nonlinear dynamics for analyzing the models. We find that autotetraploids approach Hardy-Weinberg Equilibrium 33% faster than allotetraploids, but the different tetraploid inheritance models show little differences in mutation load and allele frequency at mutation-selection balance. Our model also reveals two bistable points of mutation-selection balance for dominant alleles with biased mutation rates over a wide range of selection coefficients in the tetraploid models compared to bistability in only a narrow range for diploids. Finally, using discrete time simulations, we explore the temporal dynamics of allele frequency and fitness change and compare these dynamics to the predictions of Fisher's Fundamental Theorem of Natural Selection. While Fisher's predictions generally hold, we show that the bistable dynamics for dominant mutations fundamentally alter the associated temporal dynamics. Overall, this work develops foundational theoretical models that will facilitate the development of population genetic models and methodologies to study evolution in empirical tetraploid populations.

Fisher’s Fundamental Theorem↗

PubMed2025

Unequally Abundant Chromosomes and Unusual Collections of Transferred Sequences Characterize Mitochondrial Genomes of Gastrodia (Orchidaceae), One of the Largest Mycoheterotrophic Plant Genera.

The mystery of genomic alternations in heterotrophic plants is among the most intriguing in evolutionary biology. Compared to plastid genomes (plastomes) with parallel size reduction and gene loss, mitochondrial genome (mitogenome) variation in heterotrophic plants remains underexplored in many aspects. To further unravel the evolutionary outcomes of heterotrophy, we present a comparative mitogenomic study with 13 de novo assemblies of Gastrodia (Orchidaceae), one of the largest fully mycoheterotrophic plant genera, and its relatives. Analyzed Gastrodia mitogenomes range from 0.56 to 2.1 Mb, each consisting of numerous, unequally abundant chromosomes or contigs. Size variation might have evolved through chromosome rearrangements followed by stochastic loss of "dispensable" chromosomes, with deletion-biased mutations. The discovery of a hyper-abundant (∼15 times intragenomic average) chromosome in two assemblies represents the hitherto most extreme copy number variation in any mitogenomes, with similar architectures discovered in two metazoan lineages. Transferred sequence contents highlight asymmetric evolutionary consequences of heterotrophy: despite drastically reduced intracellular plastome transfers convergent across heterotrophic plants, their rarity of horizontally acquired sequences sharply contrasts parasitic plants, where massive transfers from their hosts prevail. Rates of sequence evolution are markedly elevated but not explained by copy number variation, extending prior findings of accelerated molecular evolution from parasitic to heterotrophic plants. Putative evolutionary scenarios for these mitogenomic convergence and divergence fit well with the common (e.g. plastome contraction) and specific (e.g. host identity) aspects of the two heterotrophic types. These idiosyncratic mycoheterotrophs expand known architectural variability of plant mitogenomes and provide mechanistic insights into their content and size variation.

Genome, Mitochondrial↗

PubMed2025

Mammalian antiviral proteins ZAP and KHNYN can independently restrict CpG-enriched avian viruses.

Zoonotic viruses are an omnipresent threat to global health. Influenza A virus (IAV) transmits between birds, livestock, and humans. Proviral host factors involved in the cross-species interface are well known. Less is known about antiviral mechanisms that suppress IAV zoonoses. We observed CpG dinucleotide depletion in human IAV relative to avian IAV. Notably, human ZAP selectively depletes CpG-enriched viral RNAs with its cofactor KHNYN. ZAP is conserved in tetrapods, but we uncovered that avian species lack KHNYN. We found that chicken ZAP may not affect IAV (PR8) or CpG-enriched IAV (PR8CG). Human ZAP or KHNYN independently restricted CpG-enriched IAV PR8CG by overexpression in chicken cells and by combined knockout in human cells. Additionally, mammalian ZAP-L and KHNYN also independently restricted an avian retrovirus (ROSV). Curiously, platypus KHNYN, the most divergent from eutherian mammals, was also capable of independent restriction of multiple diverse viruses. We suggest that some mammalian KHNYN can act as a bona fide restriction factor with cell-autonomous activity. Furthermore, we speculate that through repeated contact between avian viruses and mammalian hosts, protein changes may accompany CpG-biased mutations or reassortment to evade mammalian ZAP and KHNYN.

Animals↗

PubMed2026

Natural Selection Drives Codon Usage Bias in the Mitochondrial Genome of Ligula intestinalis (Linnaeus, 1758) Gmelin, 1790 (Cestoda: Diphyllobothriidea): Insights from Comparative Genomics and Optimal Codon Identification.

Codon usage bias (CUB) is a useful indicator of evolutionary forces shaping mitochondrial genomes. Codon usage bias in mitochondrial genomes of Diphyllobothriidae and especially in Ligula intestinalis was characterized. The roles of natural selection and mutation pressure in framing this bias were evaluated on the basis of 12 protein-coding genes in Diphyllobothriidae. The complete mitogenome (13,725 bp) of L. intestinalis comprises 12 protein-coding genes (PCGs), 22 tRNAs, and two rRNAs, all positioned on the heavy strand, and contains an overall AT content of 66.15%. The mean CAI (0.176), CBI (-0.105), and ENC (45.33) and an evident preference for U-ending codons observed in all examined genes indicate weak CUB. Neutrality, ENC, and PR2 plots consistently demonstrate that natural selection is the predominant force driving CUB and contributes approximately 56% in L. intestinalis and 83% in other Diphyllobothriidea species, with mutation pressure playing a secondary role. Phylogenetic reconstruction supported the monophyly of Diphyllobothriidea, confirmed the paraphyly of Diphyllobothrium as traditionally defined, and placed Ligula and Digramma as sister taxa. These findings clarify the evolutionary constraints governing codon usage in cestode mitogenomes and provide practical resources for codon optimization in heterologous gene expression and genetic studies of this economically important parasite.

Diphyllobothriidea↗

PubMed2025

SURROGATE SELECTION OVERSAMPLES EXPANDED T CELL CLONOTYPES.

Surrogate selection is an experimental design that without sequencing any DNA can restrict a sample of cells to those carrying certain genomic mutations. In immunological disease studies, this design may provide a relatively easy approach to enrich a lymphocyte sample with cells relevant to the disease response because the emergence of neutral mutations associates with the proliferation history of clonal subpopulations. A statistical analysis of clonotype sizes provides a structured, quantitative perspective on this useful property of surrogate selection. Our model specification couples within-clonotype birth-death processes with an exchangeable model across clonotypes. Beyond enrichment questions about the surrogate selection design, our framework enables a study of sampling properties of elementary sample diversity statistics; it also points to new statistics that may usefully measure the burden of somatic genomic alterations associated with clonal expansion. We examine statistical properties of immunological samples governed by the coupled model specification, and we illustrate calculations in surrogate selection studies of melanoma and in single-cell genomic studies of T cell repertoires.

Bayes’s rule↗

PubMed2025

Alternative splice acceptor site in MSH4 gene is responsible for male sterility conferred by ms5 in soybean.

In soybean breeding, using the recessive male-sterile ms5 gene, derived from fast neutron mutagenesis, for recurrent selection is advantageous because of the d2 locus, which controls cotyledon color in mature seeds and can be used as a phenotypic selection marker for ms5 male sterility. However, occasional self-fertilization occurs because of the elimination of d2 linkage and instability of male sterility. Elucidating the mechanism and the gene responsible for ms5 male sterility may resolve these problems. Using fine mapping with 15 simple sequence repeat (SSR) markers, we narrowed down the candidate ms5 locus to a 54-kbp region. Bulked-DNA analysis using next-generation sequencing revealed a deletion as a candidate variation in the region. This 15-bp deletion and a nucleotide substitution were identified in intron 1 of MutS homolog (GmMSH4), which modulates chromosomal recombination in meiosis. The ms5 transcript contained a novel exon with a premature termination codon. This exon originated from an alternative splice acceptor site caused by the deletion and nucleotide substitution, disrupting gene function. Co-segregation of male sterility with five independent mutations in GmMSH4 was confirmed using progeny of mutant lines. Mutations in GmMSH4 led to biased DNA partitioning during meiosis, resulting in collapsed or enlarged pollen and suggesting that ms5 male sterility is caused by the failure of pollen formation during meiosis due to the loss of function of GmMSH4. These findings could help explain the mechanism of instability of ms5 male sterility and improve the efficiency of recurrent selection using DNA markers in soybean breeding.

Glycine max↗

PubMed2025

Sex differences in the developing human cortex intersect with genetic risk of neurodevelopmental disorders.

Autism is highly heritable and diagnosed more frequently in males than females. To identify neurodevelopmental processes that might present sex-biased vulnerability, we generated transcriptomic and epigenomic profiles of cell types present in the prenatally developing human cerebral cortex of 27 males and 21 females. By intersecting sex-biased molecular signatures and genes with de novo mutations in male and female autistic probands, we reveal two points of vulnerability contributing to the sex-biased penetrance in neurodevelopmental disorders (NDDs). First, we show that NDD risk genes are biased towards higher expression in females, identifying the NDD gene MEF2C as a critical transcription factor for female-biased expression. Second, we identify a significant contribution of X chromosome genes to NDD pathobiology. We construct a gene regulatory map of X-linked risk genes to enable functional studies of genetic variants that likely disrupt gene expression in the developing brains of autistic males. Together, these results point towards an outsized contribution of the X-chromosome to both the origin of sex differences in the developing human cortex and NDD vulnerability. We propose a model where female-biased vulnerability is driven by coding variation within genes while male-biased vulnerability is driven by noncoding variation in regulatory elements that affect gene expression.

Sex differences↗

PubMed2026

Type 1 interferon perturbates clonal competition by reshaping human blood development.

Inflammation accelerates evolutionary dynamics of hematopoietic stem cells (HSCs) in clonal hematopoiesis and myeloid neoplasms. We studied HSCs, progenitors and immune cells from patients with myeloproliferative neoplasms at baseline and following interferon-α (IFNα) treatment, the only therapy to deplete mutated stem cells. We deployed single-cell multiomics methods that distinguish the IFNα effects on mutated stem cells from the admixed wild-type HSCs, with respect to their differentiation, transcriptomes, immunophenotypes and chromatin accessibility. IFNα simultaneously activated HSCs into two polarized states: a lymphoid progenitor expansion associated with an anti-inflammatory state and an inflammatory myeloid progenitor state derived from HSCs. The augmented lymphoid differentiation balanced the typical myeloproliferative-neoplasm-induced myeloid bias, associated with normalized blood counts. Somatic mutations modified the effects of IFNα on HSC differentiation and cell cycle entry rates. Clonal fitness upon IFNα exposure was due to resistance of CALR- or JAK2-mutated stem cells to differentiate into inflammatory myeloid progenitors.

Journal Article↗

PubMed2026

Codon bias variation in Staphylococcus aureus.

BACKGROUND: Staphylococcus aureus causes a multiplicity of human diseases acquired in community and healthcare settings alike around the globe. While most studies focus on coding changes to assess genome evolution and study genetic adaptation, interrogation of silent mutations in the form of synonymous codon usage bias is less well-studied. As such, understanding of patterns in codon bias at the gene and genome levels, and how codon bias impacts protein expression in S. aureus remains incomplete. METHODS: The codon bias of 2,565 protein encoding genes from NCTC 8325 was queried against all publicly available closed S. aureus genomes. Using public BioSample data, genomes were sorted by disease state, submitting institution, and collection site. Codon bias was assessed at the level of gene and genome using the codon adaptation index (CAI), calculated using 30S and 50S ribosomal genes. Gene set enrichment analysis was applied to determine associations between physiological functions, CAI gene scores, and interquartile ranges. CAI scores were also compared to an in vitro S. aureus proteomics database to correlate codon bias and protein expression. RESULTS: CAI scores varied within and between isolates at the gene and genome levels. Genes with ribosome-associated functions were most enriched among high CAI genes, and had low CAI interquartile ranges (IQR), suggesting selective pressure to maintain high expression of these genes across all S. aureus isolates. Genome sequences submitted by Aga Khan University Hospital, Nairobi, Kenya were most different from others. For the LAC USA 300 strain, CAI and protein expression were moderately positively correlated (cor = 0.534, p < 2.2e-16). CONCLUSIONS: Codon bias in S. aureus was shown to vary between gene, and to be a source of genetic variation between isolates; CAI and in vitro protein expression were positively correlated.

Staphylococcus aureus↗

PubMed2025

The mutation landscape of Daphnia obtusa reveals evolutionary forces shaping genome stability.

Spontaneous mutations are the primary source of genetic variation and play a central role in shaping evolutionary processes. To investigate mutational dynamics in Daphnia obtusa, we generated a chromosome-level genome assembly spanning 129.4 Mb across 12 chromosomes, encompassing 15,321 predicted protein-coding genes. Leveraging whole-genome sequencing of eight mutation accumulation (MA) lines propagated for an average of 482 generations (spanning over 20 years), we estimated a spontaneous single nucleotide mutation (SNM) rate of 2.23 × 10-9 and an indel mutation rate of 2.75 × 10-10 per site per generation. The SNM spectrum was strongly biased toward C:G > T:A transitions. Comparative analyses with natural population data revealed that exonic mutations observed in the MA lines were significantly less likely to be present in standing variation than intronic or intergenic mutations, suggesting that purifying selection in natural populations acts to remove deleterious alleles. We also identified 48 de novo loss-of-heterozygosity (LOH) events, comprising 8 heterozygous deletions and 40 gene conversion events. The genome-wide gene conversion rate was estimated at 2.62 × 10-5 per heterozygous site per generation. These findings provide a comprehensive view of the mutation spectrum, selective pressures, and mechanisms underlying genome stability in D. obtusa.

Daphnia obtusa↗

PubMed2025

Suppression of LKB1-mutant lung adenocarcinoma by natural killer cells from females.

BACKGROUND: This study addressed the enigma of sex differences in smoking-related lung cancer, particularly focusing on the low LKB1 mutation frequency in female patients with lung adenocarcinoma. METHODS: Sex bias was studied with a genetically engineered mouse model and various tail-vein injection models. Immune cells were analyzed by antibody-depletion study, flow cytometry, and immunofluorescence. The relevance of our findings to human disease was validated by evaluating various lung adenocarcinoma datasets. All statistical tests are 2-sided. RESULTS: A statistically significant percentage of females are resistant to LKB1-mutant tumor formation in our models, reflecting this sex difference in humans. Natural killer (NK) cells were identified as a critical factor in this sex-biased response. This sex difference was observed primarily in LKB1-mutant lung adenocarcinoma, probably due to their low major histocompatibility complex class I level, making them the ideal target for NK cells through the missing-self recognition. Although females resistant to LKB1-mutant lung adenocarcinoma formation did not have enhancement of any specific NK subpopulation, our immunofluorescence analysis revealed high numbers of NKs in female lungs even with the presence of LKB1-mutant lung adenocarcinoma. Our gene set enrichment analysis of The Cancer Genome Atlas-lung adenocarcinoma dataset also showed that female LKB1-mutant lung adenocarcinoma patients have a stronger NK-mediated response after adjusting for other male-female differences using the LKB1 wild-type lung adenocarcinoma dataset. CONCLUSION: Females have a stronger NK-mediated response against LKB1-mutant lung adenocarcinoma, which was present in our mouse model and the human lung adenocarcinoma dataset. This study revealed a novel role of NK cells in suppressing LKB1-mutant lung adenocarcinoma in females, which should be assessed in the clinical setting in the future.

Killer Cells, Natural↗

PubMed2026

First complete mitochondrial genome of Uzelothrips scabrosus (Thysanoptera: Uzelothripidae) provides insights into gene rearrangements and phylogenetic position within Terebrantia.

The family Uzelothripidae is represented by a single genus Uzelothrips and can be distinguished from others by the presence of whip-like antennae, a circular ventral sensorium on antennal segment III, a well-developed tentorium, and a membranous ovipositor. Here, we generated the first complete mitochondrial genome of Uzelothrips scabrosus (15,674 bp) using next-generation sequencing to explore the gene rearrangements and phylogenetic relationships. It consists of 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and two putative control regions. The genome exhibits strong AT bias (71.35%) with negative AT and GC skew. Codon usage analyses indicate a strong bias towards A/U-ending codons and influenced by both natural selection and mutation pressure. All PCGs were under purifying selection, with cox1 being the most conserved and nad4L the most variable. The gene order of the family Uzelothripidae is highly rearranged compared to the ancestral insect gene order. Comparative analysis revealed that gene block B was the most widely conserved, whereas the remaining gene blocks exhibited family or lineage-specific conservation patterns, reflecting extensive mitochondrial gene rearrangements during the evolution of the Thysanoptera. Moreover, 228 synapomorphic and 68 autapomorphic gene boundaries were identified across thysanopteran mitogenomes. Phylogenies indicated that the family Uzelothripidae is in a sister relationship with Stenurothripidae, and the Uzelothripidae + Stenurothripidae clade is sister to Thripidae. This study provides the first mitogenomic insights into Uzelothripidae and highlights the need for broader taxon sampling and nuclear genomic data to resolve deep evolutionary relationships within Thysanoptera.

Comparative analysis↗

PubMed2026

Evidence for the prognostic value of TP53 mutations in circulating tumor DNA across solid malignancies: a systematic review and meta-analysis.

BACKGROUND: The purpose of this meta-analysis study is to provide evidence for the clinical utility of TP53 mutations in circulating tumor DNA (ctDNA) as a prognostic biomarker. METHODS: We searched the PubMed, Embase, Cochrane, and Web of Science databases (last update May 2025) for studies on TP53 mutations in ctDNA or cfDNA as prognosis overall survival and in solid tumors. A total of 21 studies that met the criteria were utilized and data was collected regarding the authors, year of publication, study design, site of the study, number of patients, detection, mutation sample size and outcome measures were collected. The Newcastle-Ottawa Scale (NOS) was used to evaluate the quality of the study, and meta-analysis was done by using STATA 16.0. Effect sizes were in the form of hazard ratios (HR) that had 95% confidence intervals (CI). The models used were fixed-effects and random-effects based on heterogeneity. Funnel plots, and Egger's test was used to measure publication bias, and sensitivity analysis conducted through a leave-one-out method. RESULTS: A total of 21 studies (2,685 TP53-mutated patients, one unreported) showed: Mutated patients had worse progression-free survival (PFS) (HR=2.10, p=0.000; 12 studies, heterogeneity resolved after excluding Yoshida 2023), shorter OS (HR=1.74, p=0.014; 9 studies), and reduced DFS (HR=1.73, p=0.007; 3 studies), but RFS (2 items) showed no statistically significant differences. Subgroup analyses revealed: Prospective studies showed stronger PFS (HR=2.14 vs retrospective 1.90) with Japanese subgroup HR=4.90; Lung/liver cancers had higher HRs than breast. Prospective OS HR=2.25 (lung 3.14, endometrial 0.75). Retrospective DFS HR=1.89 vs Japanese breast RFS HR=4.00. Heterogeneity originated from study design, region, and cancer type variations, with no significant publication bias (Egger's test p>0.05). CONCLUSION: Current evidence suggests that TP53 mutations detected in ctDNA are significantly associated with poor prognosis in various solid tumors, particularly lung cancer. The association is robust for PFS and OS, though high heterogeneity and biological complexity warrant cautious interpretation. These findings support the potential incorporation of ctDNA-based TP53 mutation status into clinical prognostic assessment systems as an adjunctive parameter; however, further standardization of detection protocols, functional annotation of mutation types (e.g., LOF vs. GOF), incorporation of VAF and clonality analysis, and validation in large prospective multicenter cohorts are needed before routine clinical implementation. PROSPERO REGISTRATION NUMBER: CRD420251021095.

Humans↗

PubMed2026

Integrated signatures define mutational processes in prostate cancer.

Prostate cancer follows a long and heterogeneous disease course with incompletely understood aetiology1. Here we dissect the mutational processes shaping the genomes of 959 donors from the Pan Prostate Cancer Group and assess their clinical relevance. By integrating de novo extracted single-base substitution, insertion-deletion and copy-number signatures with six novel complex structural variant signatures, we identify eight integrated mutational footprints (IMFs) that collectively explain the mutational processes in 85% of primary prostate cancer genomes. IMFs were strongly influenced by regional biases in the genome, most prevalently androgen receptor-mediated mutagenesis and replication stress. Four IMFs, present in 37% of primary tumours, were significantly associated with shorter time to metastasis. These included reactive oxygen-species-driven mutagenesis and both canonical and non-canonical homologous recombination deficiency, the latter being enriched in patients of African ancestry. Extending to the metastatic setting, we found that IMFs predicted sensitivity to androgen receptor pathway inhibitors. Taken together, our study delineates the aetiologies and mutational processes that drive the genomic and clinical heterogeneity of prostate cancer, introduces IMFs as a unifying framework, and highlights their potential to improve both risk stratification and biomarker-guided treatment selection.

Journal Article↗