Mathematical Modeling of Impacts of Patient Differences on COVID-19 Lung Fibrosis Outcomes

Genomic landscape of adult testicular germ cell tumours in the 100,000 Genomes Project

Testicular germ cell tumours (TGCT), which comprise seminoma and non-seminoma subtypes, are the most common cancers in young men. In this study, we present a comprehensive whole genome sequencing analysis of adult TGCTs. Leveraging samples from participants recruited via the UK National Health Service and data from the Genomics England 100,000 Genomes Project, our results provide an extended description of genomic elements underlying TGCT pathogenesis. This catalogue offers a comprehensive, high-resolution map of copy number alterations, structural variation, and key global genome features, including mutational signatures and analysis of extrachromosomal DNA amplification. This study establishes correlations between genomic alterations and histological diversification, revealing divergent evolutionary trajectories among TGCT subtypes. By reconstructing the chronological order of driver events, we identify a subgroup of adult TGCTs undergoing relatively late whole genome duplication. Additionally, we present evidence that human leukocyte antigen loss is a more prevalent mechanism of immune disruption in seminomas. Collectively, our findings provide valuable insights into the developmental and immune modulatory processes implicated in TGCT pathogenesis and progression.

Mutational spectra are associated with bacterial niche

As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

Centriolar subdistal appendages promote double-strand break repair through homologous recombination

The centrosome is a cytoplasmic organelle with roles in microtubule organization that has also been proposed to act as a hub for cellular signaling. Some centrosomal components are required for full activation of the DNA damage response. However, whether the centrosome regulates specific DNA repair pathways is not known. Here, we show that centrosome presence is required to fully activate recombination, specifically to completely license its initial step, the so-called DNA end resection. Furthermore, we identify a centriolar structure, the subdistal appendages, and a specific factor, CEP170, as the critical centrosomal component involved in the regulation of recombination and resection. Cells lacking centrosomes or depleted for CEP170 are, consequently, hypersensitive to DNA damaging agents. Moreover, low levels of CEP170 in multiple cancer types correlate with an increase of the mutation burden associated with specific mutational signatures and a better prognosis, suggesting that changes in CEP170 can act as a mutation driver but could also be targeted to improve current oncological treatments.

Genomic and evolutionary characteristics of metastatic gastric cancer by routes

BACKGROUND

In gastric cancer (GC) patients, metastatic progression through the lymphatic, hematogenous, peritoneal, and ovarian routes, is the ultimate cause of death. However, the genomic and evolutionary characteristics of metastatic GC have not been widely evaluated.

METHODS

Whole-exome sequencing data were analyzed for 99 primary and paired metastatic gastric cancers from 15 patients who underwent gastrectomy and metastasectomy.

RESULTS

Hematogenous metastatic tumors were associated with increased chromosomal instability and de novo gain/amplification in cancer driver genes, whereas peritoneal/ovarian metastasis was linked to sustained chromosomal stability and de novo somatic mutations in driver genes. The genomic distance of the hematogenous and peritoneal metastatic tumors was found to be closer to the primary tumors than lymph node (LN) metastasis, while ovarian metastasis was closer to LN and peritoneal metastasis than the primary tumor. Two migration patterns for metastatic GCs were identified; branched and diaspora. Both molecular subtypes of the metastatic tumors, rather than the primary tumor, and their migration patterns were related to patient survival.

CONCLUSIONS

Genomic characteristics of metastatic gastric cancer is distinctive by routes and associated with patients' prognosis along with genomic evolution pattenrs, indicating that both primary and metastatic gastric cancers require genomic evaluation.

"Evolution of the mutation spectrum across a mammalian phylogeny"

Little is known about how the spectrum and etiology of germline mutagenesis might vary among mammalian species. To shed light on this mystery, we quantify variation in mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k -mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clocklike mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these clocklike signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the mutation spectrum's phylogenetic signal when fit to non-context-dependent mutation spectrum data in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.

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.

A Comprehensive Analysis of Cutaneous Melanoma Patients in Greece Based on Multi-Omic Data

Cutaneous melanoma (CM) is the most aggressive type of skin cancer, and it is characterised by high mutational load and heterogeneity. In this study, we aimed to analyse the genomic and transcriptomic profile of primary melanomas from forty-six Formalin-Fixed, Paraffin-Embedded (FFPE) tissues from Greek patients. Molecular analysis for both germline and somatic variations was performed in genomic DNA from peripheral blood and melanoma samples, respectively, exploiting whole exome and targeted sequencing, and transcriptomic analysis. Detailed clinicopathological data were also included in our analyses and previously reported associations with specific mutations were recognised. Most analysed samples (43/46) were found to harbour at least one clinically actionable somatic variant. A subset of samples was profiled at the transcriptomic level, and it was shown that specific melanoma phenotypic states could be inferred from bulk RNA isolated from FFPE primary melanoma tissue. Integrative bioinformatics analyses, including variant prioritisation, differential gene expression analysis, and functional and gene set enrichment analysis by group and per sample, were conducted and molecular circuits that are implicated in melanoma cell programmes were highlighted. Integration of mutational and transcriptomic data in CM characterisation could shed light on genes and pathways that support the maintenance of phenotypic states encrypted into heterogeneous primary tumours.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Postmortem Human Dura Mater Cells Exhibit Phenotypic, Transcriptomic and Genetic Abnormalities that Impact their Use for Disease Modeling

Patient-derived cells hold great promise for precision medicine approaches in human health. Human dermal fibroblasts have been a major source of cells for reprogramming and differentiating into specific cell types for disease modeling. Postmortem human dura mater has been suggested as a primary source of fibroblasts for in vitro modeling of neurodegenerative diseases. Although fibroblast-like cells from human and mouse dura mater have been previously described, their utility for reprogramming and direct differentiation protocols has not been fully established. In this study, cells derived from postmortem dura mater are directly compared to those from dermal biopsies of living subjects. In two instances, we have isolated and compared dermal and dural cell lines from the same subject. Notably, striking differences were observed between cells of dermal and dural origin. Compared to dermal fibroblasts, postmortem dura mater-derived cells demonstrated different morphology, slower growth rates, and a higher rate of karyotype abnormality. Dura mater-derived cells also failed to express fibroblast protein markers. When dermal fibroblasts and dura mater-derived cells from the same subject were compared, they exhibited highly divergent gene expression profiles that suggest dura mater cells originated from a mixed mural lineage. Given their postmortem origin, somatic mutation signatures of dura mater-derived cells were assessed and suggest defective DNA damage repair. This study argues for rigorous karyotyping of postmortem derived cell lines and highlights limitations of postmortem human dura mater-derived cells for modeling normal biology or disease-associated pathobiology.

Bioinformatics roadmap for therapy selection in cancer genomics

Tumour heterogeneity is one of the main characteristics of cancer and can be categorised into inter- or intratumour heterogeneity. This heterogeneity has been revealed as one of the key causes of treatment failure and relapse. Precision oncology is an emerging field that seeks to design tailored treatments for each cancer patient according to epidemiological, clinical and omics data. This discipline relies on bioinformatics tools designed to compute scores to prioritise available drugs, with the aim of helping clinicians in treatment selection. In this review, we describe the current approaches for therapy selection depending on which type of tumour heterogeneity is being targeted and the available next-generation sequencing data. We cover intertumour heterogeneity studies and individual treatment selection using genomics variants, expression data or multi-omics strategies. We also describe intratumour dissection through clonal inference and single-cell transcriptomics, in each case providing bioinformatics tools for tailored treatment selection. Finally, we discuss how these therapy selection workflows could be integrated into the clinical practice.

Integrated genomic analyses of acral and mucosal melanomas nominate novel driver genes

BACKGROUND

Acral and mucosal melanomas are aggressive subtypes of melanoma, which have a significantly lower burden of somatic mutations than cutaneous melanomas, but more frequent copy number variations, focused gene amplifications, and structural alterations. The landscapes of their genomic alterations remain to be fully characterized.

METHODS

We compiled sequencing data of 240 human acral and mucosal melanoma samples from 11 previously published studies and applied a uniform pipeline to call tumor cell content, ploidy, somatic and germline mutations, as well as CNVs, LOH, and SVs. We identified genes that are significantly mutated or recurrently affected by CNVs and implicated in oncogenesis. We further examined the difference in the frequency of recurrent pathogenic alterations between the two melanoma subtypes, correlation between pathogenic alterations, and their association with clinical features.

RESULTS

We nominated PTPRJ, mutated and homozygously deleted in 3.8% (9/240) and 0.8% (2/240) of samples, respectively, as a probable tumor suppressor gene, and FER and SKP2, amplified in 3.8% and 11.7% of samples, respectively, as probable oncogenes. We further identified a long tail of infrequent pathogenic alterations, involving genes such as CIC and LZTR1. Pathogenic germline mutations were observed on MITF, PTEN, ATM, and PRKN. We found BRAF V600E mutations in acral melanomas with fewer structural variations, suggesting that they are distinct and related to cutaneous melanomas. Amplifications of PAK1 and GAB2 were more commonly observed in acral melanomas, whereas SF3B1 R625 codon mutations were unique to mucosal melanomas (12.9%). Amplifications at 11q13-14 were frequently accompanied by fusion to a region on chromosome 6q12, revealing a recurrent novel structural rearrangement whose role remains to be elucidated.

CONCLUSIONS

Our meta-analysis expands the catalog of driver mutations in acral and mucosal melanomas, sheds new light on their pathogenesis and broadens the catalog of therapeutic targets for these difficult-to-treat cancers.

Computational methods and translational applications for targeted next-generation sequencing platforms

During the past decade, next-generation sequencing (NGS) technologies have become widely adopted in cancer research and clinical care. Common applications within the clinical setting include patient stratification into relevant molecular subtypes, identification of biomarkers of response and resistance to targeted and systemic therapies, assessment of heritable cancer risk based on known pathogenic variants, and longitudinal monitoring of treatment response. The need for efficient downstream processing and reliable interpretation of sequencing data has led to the development of novel algorithms and computational pipelines, as well as structured knowledge bases that link genomic alterations to currently available drugs and ongoing clinical trials. Cancer centers around the world use different types of targeted solid-tissue and blood based NGS assays to analyze the genomic and transcriptomic profile of patients as part of their routine clinical care. Recently, cross-institutional collaborations have led to the creation of large pooled datasets that can offer valuable insights into the genomics of rare cancers.

Signatures of copy number alterations in human cancer

Gains and losses of DNA are prevalent in cancer and emerge as a consequence of inter-related processes of replication stress, mitotic errors, spindle multipolarity and breakage-fusion-bridge cycles, among others, which may lead to chromosomal instability and aneuploidy1,2. These copy number alterations contribute to cancer initiation, progression and therapeutic resistance3-5. Here we present a conceptual framework to examine the patterns of copy number alterations in human cancer that is widely applicable to diverse data types, including whole-genome sequencing, whole-exome sequencing, reduced representation bisulfite sequencing, single-cell DNA sequencing and SNP6 microarray data. Deploying this framework to 9,873 cancers representing 33 human cancer types from The Cancer Genome Atlas6 revealed a set of 21 copy number signatures that explain the copy number patterns of 97% of samples. Seventeen copy number signatures were attributed to biological phenomena of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, chromothripsis and haploidization. The aetiologies of four copy number signatures remain unexplained. Some cancer types harbour amplicon signatures associated with extrachromosomal DNA, disease-specific survival and proto-oncogene gains such as MDM2. In contrast to base-scale mutational signatures, no copy number signature was associated with many known exogenous cancer risk factors. Our results synthesize the global landscape of copy number alterations in human cancer by revealing a diversity of mutational processes that give rise to these alterations.

Mutational signatures are markers of drug sensitivity of cancer cells

Genomic analyses have revealed mutational footprints associated with DNA maintenance gone awry, or with mutagen exposures. Because cancer therapeutics often target DNA synthesis or repair, we asked if mutational signatures make useful markers of drug sensitivity. We detect mutational signatures in cancer cell line exomes (where matched healthy tissues are not available) by adjusting for the confounding germline mutation spectra across ancestries. We identify robust associations between various mutational signatures and drug activity across cancer cell lines; these are as numerous as associations with established genetic markers such as driver gene alterations. Signatures of prior exposures to DNA damaging agents - including chemotherapy - tend to associate with drug resistance, while signatures of deficiencies in DNA repair tend to predict sensitivity towards particular therapeutics. Replication analyses across independent drug and CRISPR genetic screening data sets reveal hundreds of robust associations, which are provided as a resource for drug repurposing guided by mutational signature markers.

A natural mutator allele shapes mutation spectrum variation in mice

Although germline mutation rates and spectra can vary within and between species, common genetic modifiers of the mutation rate have not been identified in nature1. Here we searched for loci that influence germline mutagenesis using a uniquely powerful resource: a panel of recombinant inbred mouse lines known as the BXD, descended from the laboratory strains C57BL/6J (B haplotype) and DBA/2J (D haplotype). Each BXD lineage has been maintained by brother-sister mating in the near absence of natural selection, accumulating de novo mutations for up to 50 years on a known genetic background that is a unique linear mosaic of B and D haplotypes2. We show that mice inheriting D haplotypes at a quantitative trait locus on chromosome 4 accumulate C>A germline mutations at a 50% higher rate than those inheriting B haplotypes, primarily owing to the activity of a C>A-dominated mutational signature known as SBS18. The B and D quantitative trait locus haplotypes encode different alleles of Mutyh, a DNA repair gene that underlies the heritable cancer predisposition syndrome that causes colorectal tumors with a high SBS18 mutation load3,4. Both B and D Mutyh alleles are present in wild populations of Mus musculus domesticus, providing evidence that common genetic variation modulates germline mutagenesis in a model mammalian species.

Substitution mutational signatures in whole-genome-sequenced cancers in the UK population

Whole-genome sequencing (WGS) permits comprehensive cancer genome analyses, revealing mutational signatures, imprints of DNA damage and repair processes that have arisen in each patient's cancer. We performed mutational signature analyses on 12,222 WGS tumor-normal matched pairs, from patients recruited via the UK National Health Service. We contrasted our results to two independent cancer WGS datasets, the International Cancer Genome Consortium (ICGC) and Hartwig Foundation, involving 18,640 WGS cancers in total. Our analyses add 40 single and 18 double substitution signatures to the current mutational signature tally. Critically, we show for each organ, that cancers have a limited number of 'common' signatures and a long tail of 'rare' signatures. We provide a practical solution for utilizing this concept of common versus rare signatures in future analyses.

Spectrum of DNA mismatch repair failures viewed through the lens of cancer genomics and implications for therapy

Genome sequencing can be used to detect DNA repair failures in tumors and learn about underlying mechanisms. Here, we synthesize findings from genomic studies that examined deficiencies of the DNA mismatch repair (MMR) pathway. The impairment of MMR results in genome-wide hypermutation and in the 'microsatellite instability' (MSI) phenotype-occurrence of indel mutations at short tandem repeat (microsatellite) loci. The MSI status of tumors was traditionally assessed by molecular testing of a selected set of MS loci or by measuring MMR protein expression levels. Today, genomic data can provide a more complete picture of the consequences on genomic instability. Multiple computational studies examined somatic mutation distributions that result from failed DNA repair pathways in tumors. These include analyzing the commonly studied trinucleotide mutational spectra of single-nucleotide variants (SNVs), as well as of other features such as indels, structural variants, mutation clusters and regional mutation rate redistribution. The identified mutation patterns can be used to rigorously measure prevalence of MMR failures across cancer types, and potentially to subcategorize the MMR deficiencies. Diverse data sources, genomic and pre-genomic, from human and from experimental models, suggest there are different ways in which MMR can fail, and/or that the cell-type or genetic background may result in different types of MMR mutational patterns. The spectrum of MMR failures may direct cancer evolution, generating particular sets of driver mutations. Moreover, MMR affects outcomes of therapy by DNA damaging drugs, antimetabolites, nonsense-mediated mRNA decay (NMD) inhibitors, and immunotherapy by promoting either resistance or sensitivity, depending on the type of therapy.

Discovery of an unusually high number of de novo mutations in sperm of older men using duplex sequencing

De novo mutations (DNMs) are important players in heritable diseases and evolution. Of particular interest are highly recurrent DNMs associated with congenital disorders that have been described as selfish mutations expanding in the male germline, thus becoming more frequent with age. Here, we have adapted duplex sequencing (DS), an ultradeep sequencing method that renders sequence information on both DNA strands; thus, one mutation can be reliably called in millions of sequenced bases. With DS, we examined ∼4.5 kb of the FGFR3 coding region in sperm DNA from older and younger donors. We identified sites with variant allele frequencies (VAFs) of 10-4 to 10-5, with an overall mutation frequency of the region of ∼6 × 10-7 Some of the substitutions are recurrent and are found at a higher VAF in older donors than in younger ones or are found exclusively in older donors. Also, older donors harbor more mutations associated with congenital disorders. Other mutations are present in both age groups, suggesting that these might result from a different mechanism (e.g., postzygotic mosaicism). We also observe that independent of age, the frequency and deleteriousness of the mutational spectra are more similar to COSMIC than to gnomAD variants. Our approach is an important strategy to identify mutations that could be associated with a gain of function of the receptor tyrosine kinase activity, with unexplored consequences in a society with delayed fatherhood.

Mapping clustered mutations in cancer reveals APOBEC3 mutagenesis of ecDNA

Clustered somatic mutations are common in cancer genomes and previous analyses reveal several types of clustered single-base substitutions, which include doublet- and multi-base substitutions1-5, diffuse hypermutation termed omikli6, and longer strand-coordinated events termed kataegis3,7-9. Here we provide a comprehensive characterization of clustered substitutions and clustered small insertions and deletions (indels) across 2,583 whole-genome-sequenced cancers from 30 types of cancer10. Clustered mutations were highly enriched in driver genes and associated with differential gene expression and changes in overall survival. Several distinct mutational processes gave rise to clustered indels, including signatures that were enriched in tobacco smokers and homologous-recombination-deficient cancers. Doublet-base substitutions were caused by at least 12 mutational processes, whereas most multi-base substitutions were generated by either tobacco smoking or exposure to ultraviolet light. Omikli events, which have previously been attributed to APOBEC3 activity6, accounted for a large proportion of clustered substitutions; however, only 16.2% of omikli matched APOBEC3 patterns. Kataegis was generated by multiple mutational processes, and 76.1% of all kataegic events exhibited mutational patterns that are associated with the activation-induced deaminase (AID) and APOBEC3 family of deaminases. Co-occurrence of APOBEC3 kataegis and extrachromosomal DNA (ecDNA), termed kyklonas (Greek for cyclone), was found in 31% of samples with ecDNA. Multiple distinct kyklonic events were observed on most mutated ecDNA. ecDNA containing known cancer genes exhibited both positive selection and kyklonic hypermutation. Our results reveal the diversity of clustered mutational processes in human cancer and the role of APOBEC3 in recurrently mutating and fuelling the evolution of ecDNA.

The breast pre-cancer atlas illustrates the molecular and micro-environmental diversity of ductal carcinoma in situ

Microenvironmental and molecular factors mediating the progression of Breast Ductal Carcinoma In Situ (DCIS) are not well understood, impeding the development of prevention strategies and the safe testing of treatment de-escalation. We addressed methodological barriers and characterized the mutational, transcriptional, histological, and microenvironmental landscape across 85 multiple microdissected regions from 39 cases. Most somatic alterations, including whole-genome duplications, were clonal, but genetic divergence increased with physical distance. Phenotypic and subtype heterogeneity was frequently associated with underlying genetic heterogeneity and regions with low-risk features preceded those with high-risk features according to the inferred phylogeny. B- and T-lymphocytes spatial analysis identified three immune states, including an epithelial excluded state located preferentially at DCIS regions, and characterized by histological and molecular features of immune escape, independently from molecular subtypes. Such breast pre-cancer atlas with uniquely integrated observations will help scope future expansion studies and build finer models of outcomes and progression risk.

Accuracy of mutational signature software on correlated signatures

Mutational signatures are characteristic patterns of mutations generated by exogenous mutagens or by endogenous mutational processes. Mutational signatures are important for research into DNA damage and repair, aging, cancer biology, genetic toxicology, and epidemiology. Unsupervised learning can infer mutational signatures from the somatic mutations in large numbers of tumors, and separating correlated signatures is a notable challenge for this task. To investigate which methods can best meet this challenge, we assessed 18 computational methods for inferring mutational signatures on 20 synthetic data sets that incorporated varying degrees of correlated activity of two common mutational signatures. Performance varied widely, and four methods noticeably outperformed the others: hdp (based on hierarchical Dirichlet processes), SigProExtractor (based on multiple non-negative matrix factorizations over resampled data), TCSM (based on an approach used in document topic analysis), and mutSpec.NMF (also based on non-negative matrix factorization). The results underscored the complexities of mutational signature extraction, including the importance and difficulty of determining the correct number of signatures and the importance of hyperparameters. Our findings indicate directions for improvement of the software and show a need for care when interpreting results from any of these methods, including the need for assessing sensitivity of the results to input parameters.

MSA: reproducible mutational signature attribution with confidence based on simulations

BACKGROUND

Mutational signatures proved to be a useful tool for identifying patterns of mutations in genomes, often providing valuable insights about mutagenic processes or normal DNA damage. De novo extraction of signatures is commonly performed using Non-Negative Matrix Factorisation methods, however, accurate attribution of these signatures to individual samples is a distinct problem requiring uncertainty estimation, particularly in noisy scenarios or when the acting signatures have similar shapes. Whilst many packages for signature attribution exist, a few provide accuracy measures, and most are not easily reproducible nor scalable in high-performance computing environments.

RESULTS

We present Mutational Signature Attribution (MSA), a reproducible pipeline designed to assign signatures of different mutation types on a single-sample basis, using Non-Negative Least Squares method with optimisation based on configurable simulations. Parametric bootstrap is proposed as a way to measure statistical uncertainties of signature attribution. Supported mutation types include single and doublet base substitutions, indels and structural variants. Results are validated using simulations with reference COSMIC signatures, as well as randomly generated signatures.

CONCLUSIONS

MSA is a tool for optimised mutational signature attribution based on simulations, providing confidence intervals using parametric bootstrap. It comprises a set of Python scripts unified in a single Nextflow pipeline with containerisation for cross-platform reproducibility and scalability in high-performance computing environments. The tool is publicly available from https://gitlab.com/s.senkin/MSA .

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

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

Germline and Somatic Whole-Exome Sequencing Identifies New Candidate Genes Involved in Familial Predisposition to Serrated Polyposis Syndrome

The serrated polyposis syndrome (SPS) is the most common and yet underdiagnosed colorectal polyposis syndrome. It is characterized by multiple and/or large colonic serrated polyps and a higher associated risk for colorectal cancer (CRC). The main objective of this study was to identify new candidate genes involved in the germline predisposition to SPS/CRC. Thirty-nine SPS patients from 16 families (≥2 patients per family) were recruited without alterations in well-known hereditary CRC genes, and germline and somatic whole-exome sequencing were performed. Germline rare variants with plausible pathogenicity, located in genes involved in cancer development, senescence and epigenetic regulation were selected. Somatic mutational profiling and signature analysis was pursued in one sample per family, when possible. After data filtering, ANXA10, ASXL1, CFTR, DOT1L, HIC1, INO80, KLF3, MCM3AP, MCM8, PDLIM2, POLD1, TP53BP1, WNK2 and WRN were highlighted as the more promising candidate genes for SPS germline predisposition with potentially pathogenic variants shared within families. Somatic analysis characterized mutational profiles in advanced serrated polyps/tumors, revealing a high proportion of hypermutated samples, with a prevalence of clock-like mutational signatures in most samples and the presence of DNA mismatch repair-defective signatures in some cases. In conclusion, we identified new candidate genes to be involved in familial SPS. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.