The repertoire of mutational signatures in human cancer

From Genes to Clinical Practice: Exploring the Genomic Underpinnings of Endometrial Cancer

Endometrial cancer (EC), a prevalent gynecological malignancy, presents significant challenges due to its genetic complexity and heterogeneity. The genomic landscape of EC is underpinned by genetic alterations, such as mutations in PTEN, PIK3CA, and ARID1A, and chromosomal abnormalities. The identification of molecular subtypes-POLE ultramutated, microsatellite instability (MSI), copy number low, and copy number high-illustrates the diverse genetic profiles within EC and underscores the need for subtype-specific therapeutic strategies. The integration of multi-omics technologies such as single-cell genomics and spatial transcriptomics has revolutionized our understanding and approach to studying EC and offers a holistic perspective that enhances the ability to identify novel biomarkers and therapeutic targets. The translation of these multi-omics findings into personalized medicine and precision oncology is increasingly feasible in clinical practice. Targeted therapies such as PI3K/AKT/mTOR inhibitors have demonstrated the potential for improved treatment efficacy tailored to specific genetic alterations. Despite these advancements, challenges persist in terms of variability in patient responses, the integration of genomic data into clinical workflows, and ethical considerations. This review explores the genomic underpinnings of EC, from genes to clinical practice. It highlights the ongoing need for multidisciplinary research and collaboration to address the complexities of EC and improve diagnosis, treatment, and patient outcomes.

Genomic and transcriptomic signatures of sequential carcinogenesis from papillary neoplasm to biliary tract cancer

BACKGROUND & AIMS

Papillary neoplasms of the biliary tree, including intraductal papillary neoplasms (IPN) and intracholecystic papillary neoplasms (ICPN), are recognized as precancerous lesions. However, the genetic characteristics underlying sequential carcinogenesis remain unclear.

METHODS

Whole-exome sequencing was performed on 166 neoplasms (33 intrahepatic IPNs, 44 extrahepatic IPNs, and 89 ICPNs), and 41 associated carcinomas. Nine available cases were also subjected to spatial transcriptomic analysis.

RESULTS

Mutations in the MAPK (48%), genomic integrity maintenance (42%), and Wnt/β-catenin (33%) pathways were prevalent in intrahepatic IPNs, extrahepatic IPNs, and ICPNs, respectively. KRAS mutations were enriched in intrahepatic IPN (42%, P<0.001), whereas SMAD4 mutations were enriched in extrahepatic IPN (21%, P=0.005). ICPNs frequently exhibit CTNNB1 mutations, particularly in low-grade lesions. Mutational signature analysis revealed that SBS1 and SBS5 signatures were homogeneously enriched in intrahepatic IPN, in contrast to the heterogeneous distribution of SBS1, SBS2, SBS5, SBS13, SBS7b, and SBS23 in extrahepatic IPN and ICPN. Copy number aberrations gradually increased from low- to high-grade intraepithelial neoplasia and eventually to carcinoma. Phylogenetic analysis revealed that 89% of carcinomas were derived from IPN/ICPN through sequential carcinogenesis, with the majority sharing driver mutations between IPN/ICPN and carcinoma. Furthermore, multifocal, independent carcinogenesis events were observed in IPNs/ICPNs, resulting in mutationally distinct carcinoma lesions. Carcinogenesis of IPN/ICPN occurs in multiple subclones through mutational accumulation and transcriptomic alterations that affect vascular development, cell morphogenesis, extracellular matrix organization, and growth factor response.

CONCLUSIONS

With the largest IPN/ICPN cohort reported to date, our study provides a genome- and spatial transcriptome-level portrait of sequential carcinogenesis and differences in the anatomical location of biliary papillary neoplasms.

IMPACT AND IMPLICATIONS

Biliary tract cancer is a fatal malignancy. However, its genome-level sequential carcinogenesis from intraepithelial neoplasia to carcinoma has not yet been evaluated in a sufficiently large cohort. Papillary lesions of the bile duct and gallbladder are collectively termed intraductal papillary neoplasms (IPN) of the bile duct and intracholecystic papillary neoplasms (ICPN), respectively. They are primarily diagnosed based on histopathological studies. This study provides a comprehensive mutational and spatial transcriptomic landscape of papillary neoplasms of the bile duct and gallbladder. The results of this study offer insights into the mechanism of sequential carcinogenesis in papillary biliary tract tumors, pathology-genomics correlation, and potential therapeutic targets.

Molecular characterization of EBV-associated primary pulmonary lymphoepithelial carcinoma by multiomics analysis

BACKGROUND

Primary pulmonary lymphoepithelial carcinoma (pLEC) is a subtype of non-small cell lung cancer (NSCLC) characterized by Epstein-Barr virus (EBV) infection. However, the molecular pathogenesis of pLEC remains poorly understood.

METHODS

In this study, we explored pLEC using whole-exome sequencing (WES) and RNA-whole-transcriptome sequencing (RNA-seq) technologies. Datasets of normal lung tissue, other types of NSCLC, and EBV-positive nasopharyngeal carcinoma (EBV+-NPC) were obtained from public databases. Furthermore, we described the gene signatures, viral integration, cell quantification, cell death and immune infiltration of pLEC.

RESULTS

Compared with other types of NSCLC and EBV+-NPC, pLEC patients exhibited a lower somatic mutation burden and extensive copy number deletions, including 1p36.23, 3p21.1, 7q11.23, and 11q23.3. Integration of EBV associated dysregulation of gene expression, with CNV-altered regions coinciding with EBV integration sites. Specifically, ZBTB16 and ERRFI1 were downregulated by CNV loss, and the FOXD family genes were overexpressed with CNV gain. Decreased expression of the FOXD family might be associated with a favorable prognosis in pLEC patients, and these patients exhibited enhanced cytotoxicity.

CONCLUSION

Compared with other types of NSCLC and NPC, pLEC has distinct molecular characteristics. EBV integration, the aberrant expression of genes, as well as the loss of CNVs, may play a crucial role in the pathogenesis of pLEC. However, further research is needed to assess the potential role of the FOXD gene family as a biomarker.

Somatic mutation as an explanation for epigenetic aging

DNA methylation marks have recently been used to build models known as epigenetic clocks, which predict calendar age. As methylation of cytosine promotes C-to-T mutations, we hypothesized that the methylation changes observed with age should reflect the accrual of somatic mutations, and the two should yield analogous aging estimates. In an analysis of multimodal data from 9,331 human individuals, we found that CpG mutations indeed coincide with changes in methylation, not only at the mutated site but with pervasive remodeling of the methylome out to ±10 kilobases. This one-to-many mapping allows mutation-based predictions of age that agree with epigenetic clocks, including which individuals are aging more rapidly or slowly than expected. Moreover, genomic loci where mutations accumulate with age also tend to have methylation patterns that are especially predictive of age. These results suggest a close coupling between the accumulation of sporadic somatic mutations and the widespread changes in methylation observed over the course of life.

Implementing Mutational Epidemiology on a Global Scale: Lessons from Mutographs

The Mutographs Cancer Grand Challenge team aimed to discover unknown causes of cancer through mutational epidemiology, an alliance of cancer epidemiology and somatic genomics. By generating whole-genome sequences from thousands of cancers and normal tissues from more than 30 countries on five continents, it discovered unsuspected mutagenic exposures affecting millions of people, raised the possibility that some carcinogens act by altering forces of selection in tissue microenvironments rather than by mutagenesis, and demonstrated changes to the direction of somatic evolution in normal cells of the human body in response to exogenous exposures and noncancer diseases. See related article by Bressan et al., p. 16 See related article by Bhattacharjee et al., p. 28 See related article by Goodwin et al., p. 34.

Single-nucleotide-resolution genomic maps of O6-methylguanine from the glioblastoma drug temozolomide

Temozolomide kills cancer cells by forming O6-methylguanine (O6-MeG), which leads to cell cycle arrest and apoptosis. However, O6-MeG repair by O6-methylguanine-DNA methyltransferase (MGMT) contributes to drug resistance. Characterizing genomic profiles of O6-MeG could elucidate how O6-MeG accumulation is influenced by repair, but there are no methods to map genomic locations of O6-MeG. Here, we developed an immunoprecipitation- and polymerase-stalling-based method, termed O6-MeG-seq, to locate O6-MeG across the whole genome at single-nucleotide resolution. We analyzed O6-MeG formation and repair across sequence contexts and functional genomic regions in relation to MGMT expression in a glioblastoma-derived cell line. O6-MeG signatures were highly similar to mutational signatures from patients previously treated with temozolomide. Furthermore, MGMT did not preferentially repair O6-MeG with respect to sequence context, chromatin state or gene expression level, however, may protect oncogenes from mutations. Finally, we found an MGMT-independent strand bias in O6-MeG accumulation in highly expressed genes. These data provide high resolution insight on how O6-MeG formation and repair are impacted by genome structure and nucleotide sequence. Further, O6-MeG-seq is expected to enable future studies of DNA modification signatures as diagnostic markers for addressing drug resistance and preventing secondary cancers.

Translational and clinical comparison of whole genome and transcriptome to panel sequencing in precision oncology

Precision oncology offers new cancer treatment options, yet sequencing methods vary in type and scope. In this study, we compared whole-exome/whole-genome (WES/WGS) and transcriptome sequencing (TS) with broad panel sequencing by resequencing the same tumor DNA and RNA as well as normal tissue DNA for germline assessment, from 20 patients with rare or advanced tumors, who were originally sequenced by WES/WGS ± TS within the DKFZ/NCT/DKTK MASTER program from 2015 to 2020. Molecular analyses resulted in a median number of 2.5 (gene panel) to 3.5 (WES/WGS ± TS) treatment recommendations per patient. Our results showed that approximately half of the therapy recommendations (TRs) of both sequencing programs were identical, while approximately one-third of the TRs in WES/WGS ± TS relied on biomarkers not covered by the panel. Eight of 10 molecularly informed therapy implementations were supported by the panel, the remaining two were based on biomarkers absent from the panel, highlighting the potential additional clinical benefit of WGS and TS.

A cost-effective and scalable approach for DNA extraction from FFPE tissues

Genomic profiling of cancer plays an increasingly vital role for diagnosis and therapy planning. In addition, research of novel diagnostic applications such as DNA methylation profiling requires large training and validation cohorts. Currently, most diagnostic cases processed in pathology departments are stored as formalin-fixed and paraffin embedded tissue blocks (FFPE). Consequently, there is a growing demand for high-throughput extraction of nucleic acids from FFPE tissue samples. While proprietary kits are available, they are expensive and offer little flexibility. Here, we present ht-HiTE, a high-throughput implementation of a recently published and highly efficient DNA extraction protocol. This approach enables manual and automated processing of 96-well plates with a liquid handler, offers two options for purification and utilizes off-the-shelf reagents. Finally, we show that NGS and DNA methylation microarray data obtained from DNA processed with ht-HiTE are of equivalent quality as compared to a manual, kit-based approach.

TP53 Mutations and PD-L1 Amplification in Vulvar Adenocarcinoma of the Intestinal Type: Insights From Whole Exome Sequencing of 2 Cases

Vulvar adenocarcinoma of the intestinal type (VAIt) is a rare subtype of primary vulvar carcinoma, with ∼30 cases documented in the English literature. This study presents 2 new cases of HPV-independent VAIt with lymph node metastasis and discusses their clinical presentation, histopathologic features, and whole exome sequencing (WES) analysis. Both cases exhibited histologic features consistent with VAIt, including tubular, papillary, and mucinous carcinoma components. Immunohistochemical analysis showed p16 patchy staining, CDX2, CK20, and SATB2 positivity, while being negative for ER, PAX8, and CK7. WES revealed pathogenic TP53 mutations in both cases, accompanied by distinct additional mutations (GRIN2A and KDM6A in Case #1; CHD4 in Case #2). Common copy number alterations (CNAs) included TP53 loss of heterozygosity and CD274/PD-L1 amplification. However, other CNAs varied between the cases. Immunohistochemistry for p53 suggests the presence of both wild-type and mutant subclones, indicating that TP53 abnormalities may be acquired during tumor progression. Both tumors showed mutational signatures SBS1 and SBS5, associated with aging and DNA damage. Our findings deepen the understanding of the genetic events involved in the tumorigenesis of HPV-independent VAIt. Given the TP53 abnormalities and CD274/PD-L1 amplification, emerging p53-based therapies and immune checkpoint inhibitors may represent potential treatment targets. While these findings contribute to the understanding of VAIt tumorigenesis, further research is required to validate these observations in a larger cohort.

Genomic characterization reveals distinct mutational landscape of acral melanoma in East Asian

Acral melanoma, the most common melanoma subtype in East Asia, is associated with a poor prognosis. This study aims to comprehensively analyze the genomic characteristics of acral melanoma in East Asians. We conduct whole-genome sequencing of 55 acral melanoma tumors and perform data mining with relevant clinical data. Our findings reveal a unique mutational profile in East Asian acral melanoma, characterized by fewer point mutations and structural variations, a higher prevalence of NRAS mutations, and a lower frequency of BRAF mutations compared to patients of European descent. Notably, we identify previously underestimated ultraviolet radiation signatures and their significant association with BRAF and NRAS mutations. Structural rearrangement signatures indicate distinct mutational processes in BRAF-driven versus NRAS-driven tumors. We also find that homologous recombination deficiency with MAPK pathway mutations correlated with poor prognosis. The structural variations and amplifications in EP300, TERT, RAC1, and LZTR1 point to potential novel therapeutic targets tailored to East Asian populations. The high prevalence of whole-genome duplication events in BRAF/NRAS-mutated tumors suggests a synergistic carcinogenic effect that warrants further investigation. In summary, our study provides important insights into the genetic underpinnings of acral melanoma in East Asians, creating opportunities for targeted therapies.

Deep Learning for Biomarker Discovery in Cancer Genomes

Background

Genomic data is essential for clinical decision-making in precision oncology. Bioinformatic algorithms are widely used to analyze next-generation sequencing (NGS) data, but they face two major challenges. First, these pipelines are highly complex, involving multiple steps and the integration of various tools. Second, they generate features that are human-interpretable but often result in information loss by focusing only on predefined genetic properties. This limitation restricts the full potential of NGS data in biomarker extraction and slows the discovery of new biomarkers in precision oncology.

Methods

We propose an end-to-end deep learning (DL) approach for analyzing NGS data. Specifically, we developed a multiple instance learning DL framework that integrates somatic mutation sequences to predict two compound biomarkers: microsatellite instability (MSI) and homologous recombination deficiency (HRD). To achieve this, we utilized data from 3,184 cancer patients obtained from two public databases: The Cancer Genome Atlas (TCGA) and the Clinical Proteome Tumor Analysis Consortium (CPTAC).

Results

Our proposed deep learning method demonstrated high accuracy in identifying clinically relevant biomarkers. For predicting MSI status, the model achieved an accuracy of 0.98, a sensitivity of 0.95, and a specificity of 1.00 on an external validation cohort. For predicting HRD status, the model achieved an accuracy of 0.80, a sensitivity of 0.75, and a specificity of 0.86. Furthermore, the deep learning approach significantly outperformed traditional machine learning methods in both tasks (MSI accuracy, p-value = 5.11×10-18; HRD accuracy, p-value = 1.07×10-10). Using explainability techniques, we demonstrated that the model's predictions are based on biologically meaningful features, aligning with key DNA damage repair mutation signatures.

Conclusion

We demonstrate that deep learning can identify patterns in unfiltered somatic mutations without the need for manual feature extraction. This approach enhances the detection of actionable targets and paves the way for developing NGS-based biomarkers using minimally processed data.

Investigating the origins of the mutational signatures in cancer

Most of the risk factors associated with chronic and complex diseases, such as cancer, stem from exogenous and endogenous exposures experienced throughout an individual's life, collectively known as the exposome. These exposures can modify DNA, which can subsequently lead to the somatic mutations found in all normal and tumor tissues. Understanding the precise origins of specific somatic mutations has been challenging due to multitude of DNA adducts (i.e. the DNA adductome) and their diverse positions within the genome. Thus far, this limitation has prevented researchers from precisely linking exposures to DNA adducts and DNA adducts to subsequent mutational outcomes. Indeed, many common mutations observed in human cancers appear to originate from error-prone endogenous processes. Consequently, it remains unclear whether these mutations result from exposure-induced DNA adducts, or arise indirectly from endogenous processes or are a combination of both. In this review, we summarize approaches that aim to bridge our understanding of the mechanism by which exposure leads to DNA damage and then to mutation and highlight some of the remaining challenges and shortcomings to fully supporting this paradigm. We emphasize the need to integrate cellular DNA adductomics, long read-based mapping, single-molecule duplex sequencing of native DNA molecules and advanced computational analysis. This proposed holistic approach aims to unveil the causal connections between key DNA modifications and the mutational landscape, whether they originate from external exposures, internal processes or a combination of both, thereby addressing key questions in cancer biology.

Comprehensive whole-genome sequencing reveals origins of mutational signatures associated with aging, mismatch repair deficiency and temozolomide chemotherapy

In a comprehensive study to decipher the multi-layered response to the chemotherapeutic agent temozolomide (TMZ), we analyzed 427 genomes and determined mutational patterns in a collection of ∼40 isogenic DNA repair-deficient human TK6 lymphoblast cell lines. We first demonstrate that the spontaneous mutational background is very similar to the aging-associated mutational signature SBS40 and mainly caused by polymerase zeta-mediated translesion synthesis (TLS). MSH2-/- mismatch repair (MMR) knockout in conjunction with additional repair deficiencies uncovers cryptic mutational patterns. We next report how distinct mutational signatures are induced by TMZ upon sequential inactivation of DNA repair pathways, mirroring the acquisition of chemotherapy resistance by glioblastomas. The most toxic adduct induced by TMZ, O6-meG, is directly repaired by the O6-methylguanine-DNA methyltransferase (MGMT). In MGMT-/- cells, MMR leads to cell death and limits mutagenesis. MMR deficiency results in TMZ resistance, allowing the accumulation of ∼105 C > T substitutions corresponding to signature SBS11. Under these conditions, N3-methyladenine (3-meA), processed by base excision repair (BER), limits cell survival. Without BER, 3-meA is read through via error-prone TLS, causing T > A substitutions but not affecting survival. Blocking BER after abasic site formation results in large deletions and TMZ hypersensitization. Our findings reveal potential vulnerabilities of TMZ-resistant tumors.

The Distinct Roles of NEIL1 and XPA in Limiting Aflatoxin B1-Induced Mutagenesis in Mice

Dietary exposure to aflatoxin B1 (AFB1) is a risk factor for the development of hepatocellular carcinomas. Following metabolic activation, AFB1 reacts with guanines to form covalent DNA adducts, which induce high-frequency G > T transversions. The molecular signature associated with these mutational events aligns with the single-base substitution signature 24 (SBS24) in the Catalog of Somatic Mutations in Cancer database. Deficiencies in either base excision repair due to the absence of Nei-like DNA glycosylase 1 (NEIL1) or nucleotide excision repair due to the absence of xeroderma complementation group A protein (XPA) contribute to hepatocellular carcinomas in murine models. In the current study, ultra-low error duplex sequencing was used to characterize mutational profiles in liver DNAs of NEIL1-deficient, XPA-deficient, and DNA repair-proficient mice following neonatal injection of 1 mg/kg AFB1. Analyses of AFB1-induced mutations showed high cosine similarity to SBS24 regardless of repair proficiency status. The absence of NEIL1 resulted in an approximately 30% increase in the frequency of mutations, with the distribution suggesting preferential NEIL1-dependent repair of AFB1 lesions in open chromatin regions. A trend of increased mutagenesis was also observed in the absence of XPA. Consistent with the role of XPA in transcription-coupled repair, mutational profiles in XPA-deficient mice showed disruption of the transcriptional bias in mutations associated with SBS24. Implications: Our findings define the roles of DNA repair pathways in AFB1-induced mutagenesis and carcinogenesis in murine models, with these findings having implications in human health for those with base excision repair and nucleotide excision repair deficiencies.

Innovation in cancer pharmacotherapy through integrative consideration of germline and tumor genomes

Precision cancer medicine is widely established, and numerous molecularly targeted drugs for various tumor entities are approved or are in development. Personalized pharmacotherapy in oncology has so far been based primarily on tumor characteristics, for example, somatic mutations. However, the response to drug treatment also depends on pharmacological processes summarized under the term ADME (absorption, distribution, metabolism, and excretion). Variations in ADME genes have been the subject of intensive research for >5 decades, considering individual patients' genetic makeup, referred to as pharmacogenomics (PGx). The combined impact of a patient's tumor and germline genome is only partially understood and often not adequately considered in cancer therapy. This may be attributed, in part, to the lack of methods for combined analysis of both data layers. Optimized personalized cancer therapies should, therefore, aim to integrate molecular information, which derives from both the tumor and the germline genome, and taking into account existing PGx guidelines for drug therapy. Moreover, such strategies should provide the opportunity to consider genetic variants of previously unknown functional significance. Bioinformatic analysis methods and corresponding algorithms for data interpretation need to be developed to integrate PGx data in cancer therapy with a special meaning for interdisciplinary molecular tumor boards, in which cancer patients are discussed to provide evidence-based recommendations for clinical management based on individual tumor profiles. SIGNIFICANCE STATEMENT: The era of personalized oncology has seen the emergence of drugs tailored to genetic variants associated with cancer biology. However, the full potential of targeted therapy remains untapped owing to the predominant focus on acquired tumor-specific alterations. Optimized cancer care must integrate tumor and patient genomes, guided by pharmacogenomic principles. An essential prerequisite for realizing truly personalized drug treatment of cancer patients is the development of bioinformatic tools for comprehensive analysis of all data layers generated in modern precision oncology programs.

Characterizing the evolutionary dynamics of cancer proliferation in single-cell clones with SPRINTER

Proliferation is a key hallmark of cancer, but whether it differs between evolutionarily distinct clones co-existing within a tumor is unknown. We introduce the Single-cell Proliferation Rate Inference in Non-homogeneous Tumors through Evolutionary Routes (SPRINTER) algorithm that uses single-cell whole-genome DNA sequencing data to enable accurate identification and clone assignment of S- and G2-phase cells, as assessed by generating accurate ground truth data. Applied to a newly generated longitudinal, primary-metastasis-matched dataset of 14,994 non-small cell lung cancer cells, SPRINTER revealed widespread clone proliferation heterogeneity, orthogonally supported by Ki-67 staining, nuclei imaging and clinical imaging. We further demonstrated that high-proliferation clones have increased metastatic seeding potential, increased circulating tumor DNA shedding and clone-specific altered replication timing in proliferation- or metastasis-related genes associated with expression changes. Applied to previously generated datasets of 61,914 breast and ovarian cancer cells, SPRINTER revealed increased single-cell rates of different genomic variants and enrichment of proliferation-related gene amplifications in high-proliferation clones.

FBP1 controls liver cancer evolution from senescent MASH hepatocytes

Hepatocellular carcinoma (HCC) originates from differentiated hepatocytes undergoing compensatory proliferation in livers damaged by viruses or metabolic-dysfunction-associated steatohepatitis (MASH)1. While increasing HCC risk2, MASH triggers p53-dependent hepatocyte senescence3, which we found to parallel hypernutrition-induced DNA breaks. How this tumour-suppressive response is bypassed to license oncogenic mutagenesis and enable HCC evolution was previously unclear. Here we identified the gluconeogenic enzyme fructose-1,6-bisphosphatase 1 (FBP1) as a p53 target that is elevated in senescent-like MASH hepatocytes but suppressed through promoter hypermethylation and proteasomal degradation in most human HCCs. FBP1 first declines in metabolically stressed premalignant disease-associated hepatocytes and HCC progenitor cells4,5, paralleling the protumorigenic activation of AKT and NRF2. By accelerating FBP1 and p53 degradation, AKT and NRF2 enhance the proliferation and metabolic activity of previously senescent HCC progenitors. The senescence-reversing and proliferation-supportive NRF2-FBP1-AKT-p53 metabolic switch, operative in mice and humans, also enhances the accumulation of DNA-damage-induced somatic mutations needed for MASH-to-HCC progression.

Deep CRISPR mutagenesis characterizes the functional diversity of TP53 mutations

The mutational landscape of TP53, a tumor suppressor mutated in about half of all cancers, includes over 2,000 known missense mutations. To fully leverage TP53 mutation status for personalized medicine, a thorough understanding of the functional diversity of these mutations is essential. We conducted a deep mutational scan using saturation genome editing with CRISPR-mediated homology-directed repair to engineer 9,225 TP53 variants in cancer cells. This high-resolution approach, covering 94.5% of all cancer-associated TP53 missense mutations, precisely mapped the impact of individual mutations on tumor cell fitness, surpassing previous deep mutational scan studies in distinguishing benign from pathogenic variants. Our results revealed even subtle loss-of-function phenotypes and identified promising mutants for pharmacological reactivation. Moreover, we uncovered the roles of splicing alterations and nonsense-mediated messenger RNA decay in mutation-driven TP53 dysfunction. These findings underscore the power of saturation genome editing in advancing clinical TP53 variant interpretation for genetic counseling and personalized cancer therapy.

[(Over-)living with cancer: secondary malignancies (incl. genetics)]

Secondary malignancies (secondary cancers) are malignant diseases that occur at a certain time after cancer treatment. The malignant neoplasms can occur anywhere from 2 months to decades after cancer treatment. In addition, multiple tumor diseases can also develop due to a hereditary tendency to tumors. This article provides an overview of the causes, early detection and individual treatment.

Tumor Mutation Signature Reveals the Risk Factors of Lung Adenocarcinoma with EGFR or KRAS Mutation

INTRODUCTION

EGFR and KRAS mutations are frequently detected in lung adenocarcinoma (LUAD). Tumor mutational signature (TMS) determination is an approach to identify somatic mutational patterns associated with pathogenic factors. In this study, through the analysis of TMS, the underlying pathogenic factors of LUAD with EGFR and KRAS mutations were traced.

METHODS

This was a retrospective study. TMS of LUAD with KRAS and EGFR mutations from the TCGA, OncoSG, and MSK datasets was determined by two bioinformatics tools, namely the "MutationalPatterns" and "FitMS" packages. Elevated microsatellite alterations at selected tetranucleotide repeats (EMAST) of LUAD clinical specimens was analyzed using capillary electrophoresis.

RESULTS

In LUAD with KRAS mutations, TMS analysis indicated that the smoking-related SBS4 signature was enriched. For LUAD with EGFR L858R mutation, the smoking-related SBS4 signature was enriched in the Western population from the TCGA database; however, the smoking-related SBS4 signature was not obvious in Asian LUAD patients. LUAD with EGFR exon19 deletion (19Del) exhibited stronger SBS15 signature, which was related to defective DNA mismatch repair. Capillary electrophoresis analysis showed that an EMAST locus was frequently instable in LUAD with EGFR 19Del. Different from the Western population, Asian LUAD patients with EGFR mutations exhibited the enrichment of SBS1, SBS2, and SBS13 signatures, which were associated with the endogenous mutation process of cytidine deamination.

CONCLUSIONS

TMS analysis reveals that smoking is associated with LUAD with KRAS mutations. Defective DNA mismatch repair and endogenous cytidine deamination are associated with LUAD with EGFR mutations, especially for the EGFR 19Del. The endogenous mutational process is stronger in Asian LUAD patients than Western LUAD patients.

Mutational signatures and their association with cancer survival and gene expression in multiple cancer types

Different endogenous and exogenous mutational processes cause specific patterns of somatic mutations and mutational signatures. Although their biological research has been intensive, there are only rare studies assessing the possible prognostic role of mutational signatures. We used data from The Cancer Genome Atlas to study the associations between the activity of the mutational signatures and four survival endpoints in 18 types of malignancies. We further explored the prognostic differences according to, for example, the HPV status in head and neck squamous cell carcinomas and smoking status in lung cancers. The predictive power of the signatures over time was evaluated with a dynamic area under the curve model, and the links between mutational signature activities and differences in gene expression patterns were analyzed. In 12 of 18 studied cancer types, we identified at least one mutational signature whose activity predicted survival outcomes after adjusting for the established prognostic factors. For example, overall survival was associated with the activity of mutational signatures in nine cancer types and disease-specific survival in seven cancer types. The clock-like signatures SBS5 and SBS40 were most commonly associated with survival endpoints. The genes of the myosin binding protein and melanoma antigen families were among the most substantially dysregulated genes between the signatures of low and high activity. The differences in gene expression also revealed various enriched pathways. Based on these data, specific mutational signatures associate with the gene expression and have the potential to serve as strong prognostic factors in several cancer types.

TRACERx analysis identifies a role for FAT1 in regulating chromosomal instability and whole-genome doubling via Hippo signalling

Chromosomal instability (CIN) is common in solid tumours and fuels evolutionary adaptation and poor prognosis by increasing intratumour heterogeneity. Systematic characterization of driver events in the TRACERx non-small-cell lung cancer (NSCLC) cohort identified that genetic alterations in six genes, including FAT1, result in homologous recombination (HR) repair deficiencies and CIN. Using orthogonal genetic and experimental approaches, we demonstrate that FAT1 alterations are positively selected before genome doubling and associated with HR deficiency. FAT1 ablation causes persistent replication stress, an elevated mitotic failure rate, nuclear deformation and elevated structural CIN, including chromosome translocations and radial chromosomes. FAT1 loss contributes to whole-genome doubling (a form of numerical CIN) through the dysregulation of YAP1. Co-depletion of YAP1 partially rescues numerical CIN caused by FAT1 loss but does not relieve HR deficiencies, nor structural CIN. Importantly, overexpression of constitutively active YAP15SA is sufficient to induce numerical CIN. Taken together, we show that FAT1 loss in NSCLC attenuates HR and exacerbates CIN through two distinct downstream mechanisms, leading to increased tumour heterogeneity.

Molecular dependencies and genomic consequences of a global DNA damage tolerance defect

BACKGROUND

DNA damage tolerance (DDT) enables replication to continue in the presence of fork stalling lesions. In mammalian cells, DDT is regulated by two independent pathways, controlled by the polymerase REV1 and ubiquitinated PCNA, respectively.

RESULTS

To determine the molecular and genomic impact of a global DDT defect, we studied PcnaK164R/-;Rev1-/- compound mutants in mouse cells. Double-mutant cells display increased replication stress, hypersensitivity to genotoxic agents, replication speed, and repriming. A whole-genome CRISPR-Cas9 screen revealed a strict reliance of double-mutant cells on the CST complex, where CST promotes fork stability. Whole-genome sequencing indicated that this double-mutant DDT defect favors the generation of large, replication-stress inducible deletions of 0.4-4.0 kbp, defined as type 3 deletions. Junction break sites of these deletions reveal microhomology preferences of 1-2 base pairs, differing from the smaller type 1 and type 2 deletions. These differential characteristics suggest the existence of molecularly distinct deletion pathways. Type 3 deletions are abundant in human tumors, can dominate the deletion landscape, and are associated with DNA damage response status and treatment modality.

CONCLUSIONS

Our data highlight the essential contribution of the DDT system to genome maintenance and type 3 deletions as mutational signature of replication stress. The unique characteristics of type 3 deletions implicate the existence of a novel deletion pathway in mice and humans that is counteracted by DDT.

DNA replication initiation drives focal mutagenesis and rearrangements in human cancers

The rate and pattern of mutagenesis in cancer genomes is significantly influenced by DNA accessibility and active biological processes. Here we show that efficient sites of replication initiation drive and modulate specific mutational processes in cancer. Sites of replication initiation impede nucleotide excision repair in melanoma and are off-targets for activation-induced deaminase (AICDA) activity in lymphomas. Using ductal pancreatic adenocarcinoma as a cancer model, we demonstrate that the initiation of DNA synthesis is error-prone at G-quadruplex-forming sequences in tumours displaying markers of replication stress, resulting in a previously recognised but uncharacterised mutational signature. Finally, we demonstrate that replication origins serve as hotspots for genomic rearrangements, including structural and copy number variations. These findings reveal replication origins as functional determinants of tumour biology and demonstrate that replication initiation both passively and actively drives focal mutagenesis in cancer genomes.

Mutational signatures define immune and Wnt-associated subtypes of ampullary carcinoma

BACKGROUND AND OBJECTIVE

Ampullary carcinoma (AMPAC) taxonomy is based on morphology and immunohistochemistry. This classification lacks prognostic reliability and unique genetic associations. We applied an approach of integrative genomics characterising patients with AMPAC exploring molecular subtypes that may guide personalised treatments.

DESIGN

We analysed the mutational landscapes of 170 patients with AMPAC. The discovery included 110 tumour/normal pairs and the validation comprised 60 patients. In a tumour subset, we interrogated the transcriptomes and DNA methylomes. Patients were stratified based on mutational signatures and associated with molecular and clinical features. To evaluate tumour and immune cellularity, 22 tumours were independently assessed histomorphologically and by digital pathology.

RESULTS

We defined three patient clusters by mutational signatures independent of histomorphology. Cluster 1 (C1) was defined by spontaneous deamination of DNA 5-methylcytosine and defective mismatch repair. C2 and C3 were related to the activity of transcription-coupled nucleotide excision repair but C3 was further defined by the polymerase eta mutational process. C1-2 showed enrichment of Wnt pathway alterations, aberrant DNA methylation profiles, immune cell exclusion and patients with poor prognosis. These features were associated with a hypermutator phenotype caused by C>T alterations at CpGs. C3 patients with improved overall survival were associated with activation of immune-related pathways, immune infiltration and elevated expression of immunoinhibitory checkpoint genes.

CONCLUSION

Immunogenicity and Wnt pathway associations, emphasised by the mutational signatures, defined patients with prospective sensitivity to either immunotherapy or Wnt pathway inhibitors. This emphasises a novel mutational signature-based AMPAC classification with prognostic potential, suggesting prospective implications for subgroup-specific management of patients with AMPAC.

Advances in DNA/RNA Sequencing and Their Applications in Acute Myeloid Leukemia (AML)

Acute myeloid leukemia (AML) is an aggressive malignancy that poses significant challenges due to high rates of relapse and resistance to treatment, particularly in older populations. While therapeutic advances have been made, survival outcomes remain suboptimal. The evolution of DNA and RNA sequencing technologies, including whole-genome sequencing (WGS), whole-exome sequencing (WES), and RNA sequencing (RNA-Seq), has significantly enhanced our understanding of AML at the molecular level. These technologies have led to the discovery of driver mutations and transcriptomic alterations critical for improving diagnosis, prognosis, and personalized therapy development. Furthermore, single-cell RNA sequencing (scRNA-Seq) has uncovered rare subpopulations of leukemia stem cells (LSCs) contributing to disease progression and relapse. However, widespread clinical integration of these tools remains limited by costs, data complexity, and ethical challenges. This review explores recent advancements in DNA/RNA sequencing in AML and highlights both the potential and limitations of these techniques in clinical practice.

Contributing factors to the oxidation-induced mutational landscape in human cells

8-oxoguanine (8-oxoG) is a common oxidative DNA lesion that causes G > T substitutions. Determinants of local and regional differences in 8-oxoG-induced mutability across genomes are currently unknown. Here, we show DNA oxidation induces G > T substitutions and insertion/deletion (INDEL) mutations in human cells and cancers. Potassium bromate (KBrO3)-induced 8-oxoGs occur with similar sequence preferences as their derived substitutions, indicating that the reactivity of specific oxidants dictates mutation sequence specificity. While 8-oxoG occurs uniformly across chromatin, 8-oxoG-induced mutations are elevated in compact genomic regions, within nucleosomes, and at inward facing guanines within strongly positioned nucleosomes. Cryo-electron microscopy structures of OGG1-nucleosome complexes indicate that these effects originate from OGG1's ability to flip outward positioned 8-oxoG lesions into the catalytic pocket while inward facing lesions are occluded by the histone octamer. Mutation spectra from human cells with DNA repair deficiencies reveals contributions of a DNA repair network limiting 8-oxoG mutagenesis, where OGG1- and MUTYH-mediated base excision repair is supplemented by the replication-associated factors Pol η and HMCES. Transcriptional asymmetry of KBrO3-induced mutations in OGG1- and Pol η-deficient cells also demonstrates transcription-coupled repair can prevent 8-oxoG-induced mutation. Thus, oxidant chemistry, chromatin structures, and DNA repair processes combine to dictate the oxidative mutational landscape in human genomes.

Somatic mutation in autosomal dominant polycystic kidney disease revealed by deep sequencing human kidney cysts

Autosomal Dominant Polycystic Kidney Disease (ADPKD) results in progressive cysts that lead to kidney failure, and is caused by heterozygous germline variants in PKD1 or PKD2. Cyst pathogenesis is not definitively understood. Somatic second-hit mutations have been implicated in cyst pathogenesis, though technical sequencing challenges have limited investigation. We used unique molecular identifiers, high-depth massively parallel sequencing and custom analysis techniques to identify somatic second-hit mutations in 24 whole cysts from disparate regions of six human ADPKD kidneys, utilising replicate samples and orthogonal confirmation. Average depth of coverage of 1166 error-corrected reads for PKD1 and 539 reads for PKD2 was obtained. 58% (14/24) of cysts had a detectable PKD1 somatic variant, with 5/6 participants having at least one cyst with a somatic variant. We demonstrate that low-frequency somatic mutations are detectable in a proportion of cysts from end-stage ADPKD human kidneys. Further studies are required to understand the drivers of this somatic mutation.

ERBB2/HOXB13 co-amplification with interstitial loss of BRCA1 defines a unique subset of breast cancers

BACKGROUND

The HOXB13/IL17RB gene expression biomarker has been shown to predict response to adjuvant and extended endocrine therapy in patients with early-stage ER+ HER2- breast tumors. HOXB13 gene expression is the primary determinant driving the prognostic and endocrine treatment-predictive performance of the biomarker. Currently, there is limited data on HOXB13 expression in HER2+ and ER- breast cancers. Herein, we studied the expression of HOXB13 in large cohorts of HER2+ and ER- breast cancers.

METHODS

We investigated gene expression, genomic copy number, mutational signatures, and clinical outcome data in the TGGA and METABRIC breast cancer cohorts. Genomic-based gene amplification data was validated with tri-colored fluorescence in situ hybridization.

RESULTS

In the TCGA breast cancer cohort, HOXB13 gene expression was significantly higher in HER2+ versus HER2- breast cancers, and its expression was also significantly higher in the ER- versus ER+ breast cancers. HOXB13 is frequently co-gained or co-amplified with ERBB2. Joint copy gains of HOXB13 and ERBB2 occurred with low-level co-gains or high-level co-amplifications (co-amp), the latter of which is associated with an interstitial loss that includes the tumor suppressor BRCA1. ERBB2/HOXB13 co-amp tumors with interstitial BRCA1 loss exhibit a mutational signature associated with APOBEC deaminase activity and copy number signatures associated with chromothripsis and genomic instability. Among ERBB2-amplified tumors of different tissue origins, ERBB2/HOXB13 co-amp with a BRCA1 loss appeared to be enriched in breast cancer compared to other tumor types. Lastly, patients with ERBB2/HOXB13 co-amplified and BRCA1 lost tumors displayed a significantly shorter progression-free survival (PFS) than those with ERBB2-only amplifications. The difference in PFS was restricted to the ER- subset patients and this difference in PFS was not solely driven by HOXB13 gene expression.

CONCLUSIONS

HOXB13 is frequently co-gained with ERBB2 at both low-copy number level or as complex high-level amplification with relative BRCA1 loss. ERBB2/HOXB13 amplified, BRCA1-lost tumors are strongly enriched in breast cancer, and patients with such breast tumors experience a shortened PFS.

Novel structural variants that impact cell cycle genes are elucidated in metastatic gastrointestinal stromal tumors

Gastrointestinal stromal tumors (GISTs) are the most common mesenchymal neoplasm of the digestive tract. Despite multiple therapeutic advances, patients with advanced disease frequently develop resistance to tyrosine kinase inhibitors (TKIs), and therefore represent a therapeutic challenge. We employed whole genome sequencing (WGS) on three metastatic GISTs refractory to various TKIs and explored a publicly available cohort of 499 GISTs. This study sheds light on the clinical importance of alterations in cell cycle genes such as cyclin-dependent kinase 2 A (CDKN2A), and cyclin-dependent kinase 2B (CDKN2B), their frequent alteration in metastatic GISTs and their potential role in tumor progression of this neoplasm. Likewise, new structural variations were identified in cyclin-dependent kinase 12 (CDK12). Whole genome profiling of metastatic GIST provides new insights to advance precision care of the disease, focusing on new therapeutic possibilities, especially for emerging targets such as CDK12.

The role of the JunD-RhoH axis in the pathogenesis of hairy cell leukemia and its ability to identify existing therapeutics that could be repurposed to treat relapsed or refractory disease

Hairy cell leukemia (HCL) is an indolent malignancy of mature B-lymphocytes. While existing front-line therapies achieve excellent initial results, a significant number of patients relapse and become increasingly treatment resistant. A major molecular driver of HCL is aberrant interlocking expression of the transcription factor JunD and the intracellular signaling molecule RhoH. Here we discuss the molecular basis of how the JunD-RhoH axis contributes to HCL pathogenesis. We also discuss how leveraging the JunD-RhoH axis identifies CD23, CD38, CD66a, CD115, CD269, integrin β7, and MET as new potential therapeutic targets. Critically, preclinical studies have already demonstrated that targeting CD38 with isatuximab effectively treats preexisiting HCL. Isatuximab and therapeutics directed against each of the other six new HCL targets are currently in clinical use to treat other disorders. Consequently, leveraging the JunD-RhoH axis has identified a battery of therapies that could be repurposed as new means of treating relapsed or refractory HCL.

The theory of massively repeated evolution and full identifications of cancer-driving nucleotides (CDNs)

Tumorigenesis, like most complex genetic traits, is driven by the joint actions of many mutations. At the nucleotide level, such mutations are cancer-driving nucleotides (CDNs). The full sets of CDNs are necessary, and perhaps even sufficient, for the understanding and treatment of each cancer patient. Currently, only a small fraction of CDNs is known as most mutations accrued in tumors are not drivers. We now develop the theory of CDNs on the basis that cancer evolution is massively repeated in millions of individuals. Hence, any advantageous mutation should recur frequently and, conversely, any mutation that does not is either a passenger or deleterious mutation. In the TCGA cancer database (sample size n=300-1000), point mutations may recur in i out of n patients. This study explores a wide range of mutation characteristics to determine the limit of recurrences (i*) driven solely by neutral evolution. Since no neutral mutation can reach i*=3, all mutations recurring at i≥3 are CDNs. The theory shows the feasibility of identifying almost all CDNs if n increases to 100,000 for each cancer type. At present, only <10% of CDNs have been identified. When the full sets of CDNs are identified, the evolutionary mechanism of tumorigenesis in each case can be known and, importantly, gene targeted therapy will be far more effective in treatment and robust against drug resistance.

Validation of the APOBEC3A-Mediated RNA Single Base Substitution Signature and Proposal of Novel APOBEC1, APOBEC3B, and APOBEC3G RNA Signatures

Mutational signature analysis gained significant attention for providing critical insights into the underlying mutational processes for various DNA single base substitution (SBS) signatures and their associations with different cancer types. Recently, RNA single base substitution (RNA-SBS) signatures were defined and described by decomposing RNA variants found in non-small cell lung cancer. Through statistical association, they attributed Apolipoprotein B mRNA Editing Enzyme, Catalytic Polypeptide 3A (APOBEC3A) mutagenesis to the RNA-SBS2 signature. Here, we provide the first validation of an RNA-SBS mutational signature by decomposing novel exogenous and endogenous APOBEC3A RNA editing signatures into COSMICv3.4 RNA-SBS reference signatures. Additionally, we have identified novel RNA-SBS signatures for APOBEC1, APOBEC3B, and APOBEC3G.

Whole-Exome Sequencing, Mutational Signature Analysis, and Outcome in Multiple Myeloma-A Pilot Study

The complex and heterogeneous genomic landscape of multiple myeloma (MM) and many of its clinical and prognostic implications remains to be understood. In other cancers, such as breast cancer, using whole-exome sequencing (WES) and molecular signatures in clinical practice has revolutionized classification, prognostic prediction, and patient management. However, such integration is still in its early stages in MM. In this study, we analyzed WES data from 35 MM patients to identify potential mutational signatures and driver mutations correlated with clinical and cytogenetic characteristics. Our findings confirm the complex mutational spectrum and its impact on previously described ontogenetic and epigenetic pathways. They show TYW1 as a possible new potential driver gene and find no significant associations of mutational signatures with clinical findings. Further studies are needed to strengthen the role of mutational signatures in the clinical context of patients with MM to improve patient management.

Pancreatic Neuroendocrine Tumor: The Case Report of a Patient with Germline FANCD2 Mutation and Tumor Analysis Using Single-Cell RNA Sequencing

Background: Neuroendocrine neoplasms are a rare and heterogeneous group of neoplasms. Small-sized (≤2 cm) pancreatic neuroendocrine tumors (PanNETs) are of particular interest as they are often associated with aggressive behavior, with no specific prognostic or progression markers.

METHODS

This article describes a clinical case characterized by a progressive growth of nonfunctional PanNET requiring surgical treatment in a patient with a germline FANCD2 mutation, previously not reported in PanNETs. The patient underwent whole exome sequencing and single-cell RNA sequencing.

RESULTS

The patient underwent surgical treatment. We confirmed the presence of the germline mutation FANCD2 and also detected the germline mutation WNT10A. The cellular composition of the PanNET was analyzed using single-cell sequencing, and the main cell clusters were identified. We analyzed the tumor genomics, and used the data to define the effect the germline FANCD2 mutation had.

CONCLUSIONS

Analysis of the mutational status of patients with PanNET may provide additional data that may influence treatment tactics, refine the plan for monitoring such patients, and provide more information about the pathogenesis of PanNET. PanNET research using scRNA-seq data may help in predicting the effect of therapy on neuroendocrine cells with FANCD2 mutations.

Exploring the Complexity of Pan-Cancer: Gene Convergences and in silico Analyses

Cancer is a complex and multifaceted group of diseases characterized by highly intricate mechanisms of tumorigenesis and tumor progression, which complicates diagnosis, prognosis, and treatment. In recent years, targeted therapies have gained prominence by focusing on specific mutations and molecular features unique to each tumor type, offering more effective and personalized treatment options. However, it is equally critical to explore the genetic commonalities across different types of cancer, which has led to the rise of pan-cancer studies. These approaches help identify shared therapeutic targets across various tumor types, enabling the development of broader and potentially more widely applicable treatment strategies. This review aims to provide a comprehensive overview of key concepts related to tumors, including tumorigenesis processes, the tumor microenvironment, and the role of extracellular vesicles in tumor biology. Additionally, we explore the molecular interactions and mechanisms driving tumor progression, with a particular focus on the pan-cancer perspective. To achieve this, we conducted an in silico analysis using publicly available datasets, which facilitated the identification of both common and divergent genetic and molecular patterns across different tumor types. By integrating these diverse areas, this review offers a clearer and deeper understanding of the factors influencing tumorigenesis and highlights potential therapeutic targets.

Regulatory interactions between APOBEC3B N- and C-terminal domains

APOBEC3B (A3B) is implicated in DNA mutations that facilitate tumor evolution. Although structures of its individual N- and C-terminal domains (NTD and CTD) have been resolved through X-ray crystallography, the full-length A3B (fl-A3B) structure remains elusive, limiting understanding of its dynamics and mechanisms. In particular, the APOBEC3B C-terminal domain (A3Bctd) active site is frequently closed in models and structures. In this study, we built several new models of fl-A3B using integrative structural biology methods and selected a top model for further dynamical investigation. We compared dynamics of the truncated (A3Bctd) to the fl-A3B via conventional and Gaussian accelerated molecular dynamics (MD) simulations. Subsequently, we employed weighted ensemble methods to explore the fl-A3B active site opening mechanism, finding that interactions at the NTD-CTD interface enhance the opening frequency of the fl-A3B active site. Our findings shed light on the structural dynamics of fl-A3B, which may offer new avenues for therapeutic intervention in cancer.

Utilizing insights of DNA repair machinery to discover MMEJ deletions and novel mechanisms

We developed Del-read, an algorithm targeting medium-sized deletions (6-100 bp) in short-reads, which are challenging for current variant callers relying on alignment. Our focus was on Micro-Homolog mediated End Joining deletions (MMEJ-dels), prevalent in myeloid malignancies. MMEJ-dels follow a distinct pattern, occurring between two homologies, allowing us to generate a comprehensive list of MMEJ-dels in the exome. Using Del-read, we identified numerous novel germline and somatic MMEJ-dels in BEAT-AML and TCGA-breast datasets. Validation in 672 healthy individuals confirmed their presence. These novel MMEJ-dels were linked to genomic features associated with replication stress, like G-quadruplexes and minisatellite. Additionally, we observed a new category of MMEJ-dels with an imperfect-match at the flanking sequences of the homologies, suggesting a mechanism involving mispairing in homology alignment. We demonstrated robustness of the repair system despite CRISPR/Cas9-induced mismatches in the homologies. Further analysis of the canonical ASXL1 deletion revealed a diverse array of these imperfect-matches. This suggests a potentially more flexible and error-prone MMEJ repair system than previously understood. Our findings highlight Del-read's potential in uncovering previously undetected deletions and deepen our understanding of repair mechanisms.

Dose-Related Mutagenic and Clastogenic Effects of Benzo[b]fluoranthene in Mouse Somatic Tissues Detected by Duplex Sequencing and the Micronucleus Assay

Polycyclic aromatic hydrocarbons (PAHs) are common environmental pollutants that originate from the incomplete combustion of organic materials. We investigated the clastogenicity and mutagenicity of benzo[b]fluoranthene (BbF), one of 16 priority PAHs, in MutaMouse males after a 28 day oral exposure. BbF causes robust dose-dependent increases in micronucleus frequency in peripheral blood, indicative of chromosome damage. Duplex sequencing (DS), an error-corrected sequencing technology, reveals that BbF induces dose-dependent increases in mutation frequencies in bone marrow (BM) and liver. Mutagenicity is increased in intergenic relative to genic regions, suggesting a role for transcription-coupled repair of BbF-induced DNA damage. At higher doses, the maximum mutagenic response to BbF is higher in liver, which has a lower mitotic index but higher metabolic capacity than BM; however, mutagenic potency is comparable between the two tissues. BbF induces primarily C:G > A:T mutations, followed by C:G > T:A and C:G > G:C, indicating that BbF metabolites mainly target guanines and cytosines. The mutation spectrum of BbF correlates with cancer mutational signatures associated with tobacco exposure, supporting its contribution to the carcinogenicity of combustion-derived PAHs in humans. Overall, BbF's mutagenic effects are similar to benzo[a]pyrene, a well-studied mutagenic PAH. Our work showcases the utility of DS for effective mutagenicity assessment of environmental pollutants.

Genomic landscape of circulating tumor DNA in HER2-low metastatic breast cancer

The large population of HER2-low breast cancer patients necessitates further research to provide enhanced clinical guidance. In this study, we retrospectively analyzed 1071 metastatic breast cancer (MBC) patients and the circulating tumor DNA (ctDNA) to investigate clinicopathological and genetic alterations of HER2-low MBC patients. The effect of HER2-low status on different treatment modalities was explored in two prospective clinical trials (NCT03412383, NCT01917279) and a retrospective study. Our findings suggest TP53, PIK3CA, and ESR1 are frequently mutated genes in HER2-low MBC. Compared to the HER2-0 group, mutations observed in the HER2-low group are more closely associated with metabolic pathway alterations. Additionally, among patients with ERBB2 mutations and treated with pyrotinib, the HER2-low group may experience superior prognosis when compared to the HER2-0 group. Notably, we did not find any statistically significant disparity in the response to chemotherapy, endocrine therapy, or CDK4/6 inhibitor therapy between HER2-0 and HER2-low breast cancer patients. Interestingly, within the subgroup of individuals with metabolic pathway-related gene mutations, we found that HER2-low group exhibited a more favorable response to these treatments compared to HER2-0 group. Additionally, dynamic analysis showed the HER2-low MBC patients whose molecular tumor burden index decreased or achieved early clearance of ctDNA after the initial two treatment cycles, exhibited prolonged survival. Moreover, we classified HER2-low MBC into three clusters, providing a reference for subsequent treatment with enhanced precision. Our study offers valuable insights into the biology of HER2-low MBC and may provide reference for personalized treatment strategies.

Elucidating acquired PARP inhibitor resistance in advanced prostate cancer

PARP inhibition (PARPi) has anti-tumor activity against castration-resistant prostate cancer (CRPC) with homologous recombination repair (HRR) defects. However, mechanisms underlying PARPi resistance are not fully understood. While acquired mutations restoring BRCA genes are well documented, their clinical relevance, frequency, and mechanism of generation remain unclear. Moreover, how resistance emerges in BRCA2 homozygously deleted (HomDel) CRPC is unknown. Evaluating samples from patients with metastatic CRPC treated in the TOPARP-B trial, we identify reversion mutations in most BRCA2/PALB2-mutated tumors (79%) by end of treatment. Among reversions mediated by frameshift deletions, 60% are flanked by DNA microhomologies, implicating POLQ-mediated repair. The number of reversions and time of their detection associate with radiological progression-free survival and overall survival (p < 0.01). For BRCA2 HomDels, selection for rare subclones without BRCA2-HomDel is observed following PARPi, confirmed by single circulating-tumor-cell genomics, biopsy fluorescence in situ hybridization (FISH), and RNAish. These data support the need for restored HRR function in PARPi resistance.

Genetic evolution of keratinocytes to cutaneous squamous cell carcinoma

We performed multi-omic profiling of epidermal keratinocytes, precancerous actinic keratoses, and squamous cell carcinomas to understand the molecular transitions during skin carcinogenesis. Single-cell mutational analyses of normal skin cells showed that most keratinocytes have remarkably low mutation burdens, despite decades of sun exposure, however keratinocytes with TP53 or NOTCH1 mutations had substantially higher mutation burdens. These observations suggest that wild-type keratinocytes (i.e. without pathogenic mutations) are able to withstand high dosages of cumulative UV radiation, but certain pathogenic mutations break these adaptive mechanisms, priming keratinocytes for transformation by increasing their mutation rate. Mutational profiling of squamous cell carcinomas adjacent to actinic keratoses revealed TERT promoter and CDKN2A mutations emerging in actinic keratoses, whereas additional mutations inactivating ARID2 and activating the MAPK-pathway delineated the transition to squamous cell carcinomas. Surprisingly, actinic keratoses were often not related to their neighboring squamous cell carcinoma, indicating that collisions of unrelated neoplasms are common in the skin. Spatial variation in gene expression patterns was common in both tumor and immune cells, with high expression of checkpoint molecules at the invasive front of tumors. In conclusion, this study catalogues the key events during the evolution of cutaneous squamous cell carcinoma.

Optimizing design of genomics studies for clonal evolution analysis

Motivation

Genomic biotechnology has rapidly advanced, allowing for the inference and modification of genetic and epigenetic information at the single-cell level. While these tools hold enormous potential for basic and clinical research, they also raise difficult issues of how to design studies to deploy them most effectively. In designing a genomic study, a modern researcher might combine many sequencing modalities and sampling protocols, each with different utility, costs, and other tradeoffs. This is especially relevant for studies of somatic variation, which may involve highly heterogeneous cell populations whose differences can be probed via an extensive set of biotechnological tools. Efficiently deploying genomic technologies in this space will require principled ways to create study designs that recover desired genomic information while minimizing various measures of cost.

Results

The central problem this paper attempts to address is how one might create an optimal study design for a genomic analysis, with particular focus on studies involving somatic variation that occur most often with application to cancer genomics. We pose the study design problem as a stochastic constrained nonlinear optimization problem. We introduce a Bayesian optimization framework that iteratively optimizes for an objective function using surrogate modeling combined with pattern and gradient search. We demonstrate our procedure on several test cases to derive resource and study design allocations optimized for various goals and criteria, demonstrating its ability to optimize study designs efficiently across diverse scenarios.

Availability and implementation

https://github.com/CMUSchwartzLab/StudyDesignOptimization.

Comprehensive analysis of somatic mutations and structural variations in domestic pig

Understanding somatic mutations and structural variations in domestic pigs (Sus scrofa domestica) is critical due to their increasing importance as model organisms in biomedical research. In this study, we conducted a comprehensive analysis through whole-genome sequencing of skin, organs, and blood samples. By examining two pig pedigrees, we investigated the inheritance and sharedness of structural variants among fathers, mothers, and offsprings. Utilizing single-cell clonal expansion techniques, we observed significant variations in the number of somatic mutations across different tissues. An in-house developed pipeline enabled precise filtering and analysis of these mutations, resulting in the construction of individual phylogenetic trees for two pigs. These trees explored the developmental relationships between different tissues, revealing insights into clonal expansions from various anatomical locations. This study enhances the understanding of pig genomes, affirming their increasing value in clinical and genomic research, and provides a foundation for future studies in other animals, paralleling previous studies in mice and humans. This approach not only deepens our understanding of mammalian genomic variations but also strengthens the role of pigs as a crucial model in human health and disease research.

Detection of human papillomavirus (HPV) in malignant melanoma

The most common type of melanoma is cutaneous melanoma (CM). The predominant mutational signature is that of ultraviolet radiation (UVR) exposure. The Cancer Genome Atlas (TCGA) molecular classification includes four major subtypes of CM based on common genetic alterations involving the following genes: BRAF, NRAS, and NF1, with a small fraction being "triple" wild-type. The two main signaling pathway abnormalities in CM are the mitogen-activated protein kinase (MAPK) pathway and the phosphoinositol-3-kinase (PI3K) pathway. Other less common types include mucosal melanomas (MM) and uveal melanoma (UM), which have a significantly different genomic landscape. Although few studies reported rare cases with HPV-positive (HPV+) melanoma, the clinicopathological and molecular characteristic of this entity has not been well-described. Among the 2084 melanoma cases queried at our institution, we identified seven patients diagnosed with HPV+ melanoma (prevalence 0.03 %), including five instances of CM and two of MM. The majority of cases were positive for HPV16 (n = 6). Most of the patients were elderly and with advanced disease (n = 6), although this finding may be attributed to the relative frequency of our institution testing advanced-stage tumors. Histologically, most cases showed high degree of pleomorphism and high mitotic count (5 or more mitoses/mm2) (n = 6). UVR signature was present in the CM, but not in the MM cases. Alterations in either MAPK and/or PI3K pathways were detected in the majority of cases (n = 6). The most common genetic abnormalities detected in this study occurred in the TERT promoter (TERTp) (n = 5), a finding that has been reported to be associated with aggressive disease. Our data shows that while HPV+ melanoma is rare, identifying this disease entity could help guide therapy given the demonstrated genomic alterations.

Assessing pathogenicity of mismatch repair variants of uncertain significance by molecular tumor analysis

Functional analyses are the main method to classify mismatch repair (MMR) gene variants of uncertain significance (VUSs). However, the pathogenicity remains unclear for many variants because of conflicting results between clinical, molecular, and functional data. In this study, we evaluated whether whole exome sequencing (WES) could add another layer of evidence to elucidate the pathogenicity of MMR variants with conflicting interpretations. WES was performed on formalin-fixed paraffin-embedded tumor tissue of eight patients with a constitutional MMR VUS (seven families), including eight colorectal and two endometrial carcinomas and one ovarian carcinoma. Cell-free CIMRA assays were performed to assign Odds of Pathogenicity to these VUSs. In four families, seven tumors showed MMR deficiency-associated mutational signatures, supporting the pathogenicity of the VUS. Moreover, somatic (second) MMR hits identified in the WES data were found to explain MMR staining patterns when the MMR staining was discordant with the reported germline MMR gene variant. In conclusion, WES did not significantly reclassify VUS in these cases but clarified some phenotypic aspects such as age of onset and explanations in case of discordant MMR stainings.

Evolution of genome and immunogenome in esophageal squamous cell carcinomas driven by neoadjuvant chemoradiotherapy

Neoadjuvant chemoradiotherapy (NCRT) followed by surgery is a standard treatment for locally advanced esophageal squamous cell carcinomas (ESCCs). However, the evolution of genome and immunogenome in ESCCs driven by NCRT remains incompletely elucidated. We performed whole-exome sequencing of 51 ESCC tumors collected before and after NCRT, 36 of which were subjected to transcriptome sequencing. Clonal analysis identified clonal extinction in 13 ESCC patients wherein all pre-NCRT clones disappeared after NCRT, and clonal persistence in 9 patients wherein clones endured following NCRT. The clone-persistent patients showed higher pre-NCRT genomic intratumoral heterogeneity and worse prognosis than the clone-extinct ones. In contrast to the clone-extinct patients, the clone-persistent patients demonstrated a high proportion of subclonal neoantigens within pre-treatment specimens. Transcriptome analysis revealed increased immune infiltrations and up-regulated immune-related pathways after NCRT, especially in the clone-extinct patients. The number of T cell receptor-neoantigen interactions was higher in the clone-extinct patients than in the clone-persistent ones. The decrease in T cell repertoire evenness positively correlated to the decreased number of clonal neoantigens after NCRT, especially in the clone-extinct patients. In conclusion, we identified two prognosis-related clonal dynamic modes driven by NCRT in ESCCs. This study extended our knowledge of the ESCC genome and immunogenome evolutions driven by NCRT.

Whole genome profiling of rare pediatric thoracic tumors elucidates a YAP1::LEUTX fusion in an unclassified biphasic embryonal neoplasm

Malignant biphasic tumors of the lungs are rare, more so in the pediatric population. Here, we present the whole-genome characterization of a pleuropulmonary blastoma Type III and an unclassified biphasic thoracic embryonal neoplasm. The pleuropulmonary blastoma harbored pathogenic DICER1 germline and somatic mutations, and additional somatic variants in TP53 and BCOR. The other malignant tumor demonstrated a t(11;19) balanced translocation with a YAP1::LEUTX fusion that was confirmed by fluorescence in situ hybridization. No DICER1 germline or somatic mutation was present. YAP1 and LEUTX have been implicated in tumorigenesis of various neoplasms, and YAP1 fusion genes are an emerging oncogenic entity in a variety of malignancies. In this study we highlight the importance of whole-genome characterization of rare and unclassified tumors to identify biologic mechanisms and potential therapeutic targets.

Good response of stage IV melanoma to high‑dose radiation therapy combined with immunotherapy: A case report

Patients with advanced malignant melanoma (MM) often do not receive satisfactory treatment. The present study reports the case of a 51-year-old female patient with stage IV MM of unknown primary. After undergoing immune checkpoint inhibitor therapy, the patient received multiple doses of hypofractionated radiotherapy (HFRT) for the left inguinal lymph node and single-fraction high-dose-rate brachytherapy for the left and right lung metastases. After combination treatment, the patient experienced almost complete remission of the inguinal target area, significant relief of pain and discomfort and an improved quality of life. The time of lung radiotherapy lesion control was 8 months. Meanwhile, the observed lesions (observation lesions 1, 2, 3 and 5) adjacent to the target lesion received lower doses of scattering (0.9-1.8 Gy) and the time of control for these lung observation lesions was 9 months. In addition, restarting targeted therapy after cessation of other treatments due to myelosuppression resulted in a progression-free survival time of 6 months. Nevertheless, the patient developed new metastases in the brain and abdomen. The present case report demonstrates that high-dose radiotherapy combined with immunotherapy may be effective for local lesions and that multiple doses of HFRT may be superior to single-fraction high-dose-rate brachytherapy for certain patients. Low-dose scattering also shows improvement for local lesions. Furthermore, restarting targeted therapy may be effective in the presence of target sites. Thus, the present case report provides a possible therapeutic option for the treatment of advanced melanoma.

Methods and applications of genome-wide profiling of DNA damage and rare mutations

DNA damage is a threat to genome integrity and can be a cause of many human diseases, owing to either changes in the chemical structure of DNA or conversion of the damage into a mutation, that is, a permanent change in DNA sequence. Determining the exact positions of DNA damage and ensuing mutations in the genome are important for identifying mechanisms of disease aetiology when characteristic mutations are prevalent and probably causative in a particular disease. However, this approach is challenging particularly when levels of DNA damage are low, for example, as a result of chronic exposure to environmental agents or certain endogenous processes, such as the generation of reactive oxygen species. Over the past few years, a comprehensive toolbox of genome-wide methods has been developed for the detection of DNA damage and rare mutations at single-nucleotide resolution in mammalian cells. Here, we review and compare these methods, describe their current applications and discuss future research questions that can now be addressed.