Mutation hotspots at CTCF binding sites coupled to chromosomal instability in gastrointestinal cancers

ERCC2 mutations alter the genomic distribution pattern of somatic mutations and are independently prognostic in bladder cancer

Excision repair cross-complementation group 2 (ERCC2) encodes the DNA helicase xeroderma pigmentosum group D, which functions in transcription and nucleotide excision repair. Point mutations in ERCC2 are putative drivers in around 10% of bladder cancers (BLCAs) and a potential positive biomarker for cisplatin therapy response. Nevertheless, the prognostic significance directly attributed to ERCC2 mutations and its pathogenic role in genome instability remain poorly understood. We first demonstrated that mutant ERCC2 is an independent predictor of prognosis in BLCA. We then examined its impact on the somatic mutational landscape using a cohort of ERCC2 wild-type (n = 343) and mutant (n = 39) BLCA whole genomes. The genome-wide distribution of somatic mutations is significantly altered in ERCC2 mutants, including T[C>T]N enrichment, altered replication time correlations, and CTCF-cohesin binding site mutation hotspots. We leverage these alterations to develop a machine learning model for predicting pathogenic ERCC2 mutations, which may be useful to inform treatment of patients with BLCA.

Machine learning enables pan-cancer identification of mutational hotspots at persistent CTCF binding sites

CCCTC-binding factor (CTCF) is an insulator protein that binds to a highly conserved DNA motif and facilitates regulation of three-dimensional (3D) nuclear architecture and transcription. CTCF binding sites (CTCF-BSs) reside in non-coding DNA and are frequently mutated in cancer. Our previous study identified a small subclass of CTCF-BSs that are resistant to CTCF knock down, termed persistent CTCF binding sites (P-CTCF-BSs). P-CTCF-BSs show high binding conservation and potentially regulate cell-type constitutive 3D chromatin architecture. Here, using ICGC sequencing data we made the striking observation that P-CTCF-BSs display a highly elevated mutation rate in breast and prostate cancer when compared to all CTCF-BSs. To address whether P-CTCF-BS mutations are also enriched in other cell-types, we developed CTCF-INSITE-a tool utilising machine learning to predict persistence based on genetic and epigenetic features of experimentally-determined P-CTCF-BSs. Notably, predicted P-CTCF-BSs also show a significantly elevated mutational burden in all 12 cancer-types tested. Enrichment was even stronger for P-CTCF-BS mutations with predicted functional impact to CTCF binding and chromatin looping. Using in vitro binding assays we validated that P-CTCF-BS cancer mutations, predicted to be disruptive, indeed reduced CTCF binding. Together this study reveals a new subclass of cancer specific CTCF-BS DNA mutations and provides insights into their importance in genome organization in a pan-cancer setting.

Correlation of syndecan gene amplification with metastatic potential and clinical outcomes in carcinomas

Cell surface receptors play crucial roles in cellular responses to extracellular ligands, helping to modulate the functions of a cell based on information coming from outside the cell. Syndecan refers to a family of cell adhesion receptors that regulate both extracellular and cytosolic events. Alteration of syndecan expression disrupts regulatory mechanisms in a cell type-specific fashion, often leading to serious diseases, notably cancer. Given the multifaceted functions and distinct tissue distributions of syndecan, it will be important to unravel the gene-level intricacies of syndecan expression and thereby further understand its involvement in various carcinogenic processes. Although accumulating evidence indicates that the protein expression patterns of syndecan family members are significantly altered in cancer cells, the underlying gene-level mechanisms remain largely unknown. This review endeavors to explore syndecan gene expression levels across different cancer types by scrutinizing extensive cancer genome datasets using tools such as cBioPortal. Our analysis unveils that somatic mutations in SDC genes are rare occurrences, whereas copy number alterations are frequently observed across diverse cancers, particularly in SDC2 and SDC4. Notably, amplifications of SDC2 and SDC4 correlate with heightened metastatic potential and dismal prognosis. This underscores the recurrent nature of SDC2 and SDC4 amplifications during carcinogenesis and sheds light on their role in promoting cancer activity through augmented protein expression. The identification of these amplifications not only enriches our understanding of carcinogenic mechanisms but also hints at the potential therapeutic avenue of targeting SDC2 and SDC4 to curb cancer cell proliferation and metastasis.

Hidden secrets of the cancer genome: unlocking the impact of non-coding mutations in gene regulatory elements

Discoveries in the field of genomics have revealed that non-coding genomic regions are not merely "junk DNA", but rather comprise critical elements involved in gene expression. These gene regulatory elements (GREs) include enhancers, insulators, silencers, and gene promoters. Notably, new evidence shows how mutations within these regions substantially influence gene expression programs, especially in the context of cancer. Advances in high-throughput sequencing technologies have accelerated the identification of somatic and germline single nucleotide mutations in non-coding genomic regions. This review provides an overview of somatic and germline non-coding single nucleotide alterations affecting transcription factor binding sites in GREs, specifically involved in cancer biology. It also summarizes the technologies available for exploring GREs and the challenges associated with studying and characterizing non-coding single nucleotide mutations. Understanding the role of GRE alterations in cancer is essential for improving diagnostic and prognostic capabilities in the precision medicine era, leading to enhanced patient-centered clinical outcomes.

Efficacy of immunotherapy in ARID1A-mutant solid tumors: a single-center retrospective study

BACKGROUND

Immune checkpoint inhibitors (ICIs), especially those targeting programmed cell death-1 (PD-1) and programmed cell death ligand-1 (PD-L1), have introduced a new treatment landscape for many types of tumors. However, they only achieve a limited therapeutic response. Hence, identifying patients who may benefit from ICIs is currently a challenge.

METHODS

47 tumor patients harboring ARID1A mutations were retrospectively studied. The genomic profiling data through next-generation sequencing (NGS) and relevant clinical information were collected and analyzed. Additionally, bioinformatics analysis of the expression of immune checkpoints and immune cell infiltration levels was conducted in ARID1A-mutant gastric cancer (GC).

RESULTS

ARID1A mutations frequently co-occur with mutations in DNA damage repair (DDR)-associated genes. Among the 35 ARID1A-mutant patients who received immunotherapy, 27 were evaluable., with the objective response rate (ORR) was 48.15% (13/27), and the disease control rate (DCR) was 92.59% (25/27). Moreover, survival assays revealed that ARID1A-mutant patients had longer median overall survival (mOS) after immunotherapy. In ARID1A-mutated GC patients, receiving ICIs treatment indicated longer progressive-free survival (PFS). Additionally, the incidence of microsatellite instability-high (MSI-H), high tumor mutation burden (TMB-H) and Epstein‒Barr virus (EBV) infection was elevated. Bioinformatic analysis showed significant enrichment of immune response and T cell activation pathway within differentially expressed genes in ARID1A-mutant GC group. Finally, ARID1A mutations status was considered to be highly correlated with the level of tumor infiltrating lymphocytes (TILs) and high expression of immune checkpoints.

CONCLUSIONS

Patients with tumors harboring ARID1A mutations may achieve better clinical outcomes from immunotherapy, especially in GC. ARID1A mutations can lead to genomic instability and reshape the tumor immune microenvironment (TIME), which can be used as a biomarker for immunotherapy.

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation - a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

Prediction of novel biomarkers for gastric intestinal metaplasia and gastric adenocarcinoma using bioinformatics analysis

Background & aim

The histologic and molecular changes from intestinal metaplasia (IM) to gastric cancer (GC) have not been fully characterized. The present study sought to identify potential alterations in signaling pathways in IM and GC to predict disease progression; these alterations can be considered therapeutic targets.

Materials & methods

Seven gene expression profiles were selected from the GEO database. Discriminate differentially expressed genes (DEGs) were analyzed by EnrichR. The STRING database, Cytoscape, Gene Expression Profiling Interactive Analysis (GEPIA), cBioPortal, NetworkAnalyst, MirWalk database, OncomiR, and bipartite miRNA‒mRNA correlation network was used for downstream analyses of selected module genes.

Results

Analyses revealed that extracellular matrix-receptor interactions (ITGB1, COL1A1, COL1A2, COL4A1, FN1, COL6A3, and THBS2) in GC and PPAR signaling pathway interactions (FABP1, APOC3, APOA1, HMGCS2, and PPARA and PCK1) in IM may play key roles in both the carcinogenesis and progression of underlying GC from intestinal metaplasia. IM enrichment indicated that this is closely related to digestion and absorption. The TF-hub gene regulatory network revealed that AR, TCF4, SALL4, and ESR1 were more important for hub gene expression. It was revealed that the development and prediction of GC may be affected by hsa-miR-29. It was found that PTGR1, C1orf115, CRYL1, ALDOB, and SULT1B1 were downregulated in GC and upregulated in IM. Therefore, they might have tumor suppressor activity in GC progression.

Conclusion

New potential biomarkers and pathways involved in GC and IM were identified that are important for the transformation of GC from IM to adenocarcinoma and can be therapeutic targets for GC.

Extended Enrichment for Ultrasensitive Detection of Low-Frequency Mutations by Long Blocker Displacement Amplification

Detecting low-frequency DNA mutations hotspots cluster is critical for cancer diagnosis but remains challenging. Quantitative PCR (qPCR) is constrained by sensitivity, and allele-specific PCR is restricted by throughput. Here we develop a long blocker displacement amplification (LBDA) coupled with qPCR for ultrasensitive and multiplexed variants detection. By designing long blocker oligos to perfectly match wildtype sequences while mispairing with mutants, long blockers enable 14-44 nt enrichment regions which is 2-fold longer than normal BDA in the experiments. For wild template with a specific nucleotide, LBDA can detect different mutation types down to 0.5 % variant allele frequency (VAF) in one reaction, with median enrichment fold of 1,000 on 21 mutant DNA templates compared to the wild type. We applied LBDA-qPCR to detect KRAS and NRAS mutation hotspots, utilizing a single plex assay capable of covering 81 mutations and tested in synthetic templates and colorectal cancer tissue samples. Moreover, the mutation types were verified through Sanger sequencing, demonstrating concordance with results obtained from next generation sequencing. Overall, LBDA-qPCR provides a simple yet ultrasensitive approach for multiplexed detection of low VAF mutations hotspots, presenting a powerful tool for cancer diagnosis and monitoring.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

Dissection of a CTCF topological boundary uncovers principles of enhancer-oncogene regulation

Enhancer-gene communication is dependent on topologically associating domains (TADs) and boundaries enforced by the CCCTC-binding factor (CTCF) insulator, but the underlying structures and mechanisms remain controversial. Here, we investigate a boundary that typically insulates fibroblast growth factor (FGF) oncogenes but is disrupted by DNA hypermethylation in gastrointestinal stromal tumors (GISTs). The boundary contains an array of CTCF sites that enforce adjacent TADs, one containing FGF genes and the other containing ANO1 and its putative enhancers, which are specifically active in GIST and its likely cell of origin. We show that coordinate disruption of four CTCF motifs in the boundary fuses the adjacent TADs, allows the ANO1 enhancer to contact FGF3, and causes its robust induction. High-resolution micro-C maps reveal specific contact between transcription initiation sites in the ANO1 enhancer and FGF3 promoter that quantitatively scales with FGF3 induction such that modest changes in contact frequency result in strong changes in expression, consistent with a causal relationship.

Enhancer-promoter specificity in gene transcription: molecular mechanisms and disease associations

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

Hotspot propensity across mutational processes

The sparsity of mutations observed across tumours hinders our ability to study mutation rate variability at nucleotide resolution. To circumvent this, here we investigated the propensity of mutational processes to form mutational hotspots as a readout of their mutation rate variability at single base resolution. Mutational signatures 1 and 17 have the highest hotspot propensity (5-78 times higher than other processes). After accounting for trinucleotide mutational probabilities, sequence composition and mutational heterogeneity at 10 Kbp, most (94-95%) signature 17 hotspots remain unexplained, suggesting a significant role of local genomic features. For signature 1, the inclusion of genome-wide distribution of methylated CpG sites into models can explain most (80-100%) of the hotspot propensity. There is an increased hotspot propensity of signature 1 in normal tissues and de novo germline mutations. We demonstrate that hotspot propensity is a useful readout to assess the accuracy of mutation rate models at nucleotide resolution. This new approach and the findings derived from it open up new avenues for a range of somatic and germline studies investigating and modelling mutagenesis.

Involvement of CCCTC-binding factor in epigenetic regulation of cancer

A major global health burden continues to be borne by the complex and multifaceted disease of cancer. Epigenetic changes, which are essential for the emergence and spread of cancer, have drawn a huge amount of attention recently. The CCCTC-binding factor (CTCF), which takes part in a wide range of cellular processes including genomic imprinting, X chromosome inactivation, 3D chromatin architecture, local modifications of histone, and RNA polymerase II-mediated gene transcription, stands out among the diverse array of epigenetic regulators. CTCF not only functions as an architectural protein but also modulates DNA methylation and histone modifications. Epigenetic regulation of cancer has already been the focus of plenty of studies. Understanding the role of CTCF in the cancer epigenetic landscape may lead to the development of novel targeted therapeutic strategies for cancer. CTCF has already earned its status as a tumor suppressor gene by acting like a homeostatic regulator of genome integrity and function. Moreover, CTCF has a direct effect on many important transcriptional regulators that control the cell cycle, apoptosis, senescence, and differentiation. As we learn more about CTCF-mediated epigenetic modifications and transcriptional regulations, the possibility of utilizing CTCF as a diagnostic marker and therapeutic target for cancer will also increase. Thus, the current review intends to promote personalized and precision-based therapeutics for cancer patients by shedding light on the complex interplay between CTCF and epigenetic processes.

High CTCF expression mediated by FGD5-AS1/miR-19a-3p axis is associated with immunosuppression and pancreatic cancer progression

The most common reason for cancer-related death globally is predicted to be pancreatic cancer (PC), one of the deadliest cancers. The CCCTC-binding factor (CTCF) regulates the three-dimensional structure of chromatin, was reported to be highly regulated in various malignancies. However, the underlying biological functions and possible pathways via which CTCF promotes PC progression remain unclear. Herein, we examined the CTCF function in PC and discovered that CTCF expression in PC tissues was significantly raised compared to neighboring healthy tissues. Additionally, Kaplan-Meier survival analysis demonstrated a strong connection between elevated CTCF expression and poor patient prognosis. A study of the ROC curve (receiver operating characteristic) revealed an AUC value for CTCF of 0.968. Subsequent correlation analysis exhibited a strong relationship between immunosuppression and CTCF expression in PC. CTCF knockdown significantly inhibited the malignant biological process of PC in vitro and in vivo, suggesting that CTCF may be a potential PC treatment target. We also identified the FGD5 antisense RNA 1 (FGD5-AS1)/miR-19a-3p axis as a possible upstream mechanism for CTCF overexpression. In conclusion, our data suggest that ceRNA-mediated CTCF overexpression contributes to the suppression of anti-tumor immune responses in PC and could be a predictive biomarker and potential PC treatment target.

Enhancers dysfunction in the 3D genome of cancer cells

Eukaryotic genomes are spatially organized inside the cell nucleus, forming a threedimensional (3D) architecture that allows for spatial separation of nuclear processes and for controlled expression of genes required for cell identity specification and tissue homeostasis. Hence, it is of no surprise that mis-regulation of genome architecture through rearrangements of the linear genome sequence or epigenetic perturbations are often linked to aberrant gene expression programs in tumor cells. Increasing research efforts have shed light into the causes and consequences of alterations of 3D genome organization. In this review, we summarize the current knowledge on how 3D genome architecture is dysregulated in cancer, with a focus on enhancer highjacking events and their contribution to tumorigenesis. Studying the functional effects of genome architecture perturbations on gene expression in cancer offers a unique opportunity for a deeper understanding of tumor biology and sets the basis for the discovery of novel therapeutic targets.

Predicting regional somatic mutation rates using DNA motifs

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

ERBB2-Mutant Gastrointestinal Tumors Represent Heterogeneous Molecular Biology, Particularly in Microsatellite Instability, Tumor Mutation Burden, and Co-Mutated Genes: An In Silico Study

During recent years, activating mutations in ERBB2 have been reported in solid tumors of various organs, and clinical trials targeting ERBB2-mutant tumors have been conducted. However, no effective treatment has been established for gastrointestinal tumors targeting ERBB2 mutations. ERBB2-mutant tumors have a higher tumor mutation burden (TMB) and microsatellite instability (MSI) than ERBB2 non-mutant tumors, but not all ERBB2-mutant tumors are TMB- and MSI-high. Thus, a more detailed classification of ERBB2-mutant tumors based on the underlying molecular mechanisms is required. Herein, we classified ERBB2 mutations into three groups-group 1: both ERBB2 mutations and amplifications; group 2: ERBB2 mutations annotated as putative driver mutations but without amplifications; group 3: ERBB2 mutations annotated as non-driver mutations (passenger mutations or unknown significance) and those that were not amplified in gastrointestinal tumors. Esophageal adenocarcinoma, gastric cancer, and colorectal cancer presented significantly higher MSI and TMB in the ERBB2-mutant group than in the ERBB2-wild-type group. The proportions of TMB- and MSI-high tumors and frequency of co-mutated downstream genes differed among the groups. We identified TMB- and MSI-high groups; this classification is considered important for guiding the selection of drugs for ERBB2-mutant tumors with downstream genetic mutations.

The epigenetic modifier lysine methyltransferase 2C is frequently mutated in gastric remnant carcinoma

Gastric remnant carcinoma (GRC), which occurs in the stomach after partial gastrectomy, is a rare and aggressive form of gastric adenocarcinoma (GAC). Comprehensive profiling of genomic mutations in GRC could provide the basis for elucidating the origin and characteristics of this cancer. Herein, whole-exome sequencing (WES) was performed on 36 matched tumor-normal samples from patients with GRC and identified recurrent mutations in epigenetic modifiers, notably KMT2C, ARID1A, NSD1, and KMT2D, in 61.11% of cases. Mutational signature analysis revealed a low frequency of microsatellite instability (MSI) in GRC, which was further identified by MSIsensor, MSI-polymerase chain reaction, and immunohistochemistry analysis. Comparative analysis demonstrated that GRC had a distinct mutation spectrum compared to that of GAC in The Cancer Genome Atlas samples, with a significantly higher mutation rate of KMT2C. Targeted deep sequencing (Target-seq) of an additional 25 paired tumor-normal samples verified the high mutation frequency (48%) of KMT2C in GRC. KMT2C mutations correlated with poor overall survival in both WES and Target-seq cohorts and were independent prognosticators in GRC. In addition, KMT2C mutations were positively correlated with favorable outcomes in immune checkpoint inhibitor-treated pan-cancer patients and associated with higher intratumoral CD3+ , CD8+ tumor-infiltrating lymphocyte counts, and PD-L1 expression in GRC samples (p = 0.018, 0.092, 0.047, 0.010, and 0.034, respectively). Our dataset provides a platform for information and knowledge mining of the genomic characteristics of GRC and helps to frame new therapeutic approaches for this disease.

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

Human papillomavirus integration transforms chromatin to drive oncogenesis

BACKGROUND

Human papillomavirus (HPV) drives almost all cervical cancers and up to 70% of head and neck cancers. Frequent integration into the host genome occurs predominantly in tumorigenic types of HPV. We hypothesize that changes in chromatin state at the location of integration can result in changes in gene expression that contribute to the tumorigenicity of HPV.

RESULTS

We find that viral integration events often occur along with changes in chromatin state and expression of genes near the integration site. We investigate whether introduction of new transcription factor binding sites due to HPV integration could invoke these changes. Some regions within the HPV genome, particularly the position of a conserved CTCF binding site, show enriched chromatin accessibility signal. ChIP-seq reveals that the conserved CTCF binding site within the HPV genome binds CTCF in 4 HPV^+ cancer cell lines. Significant changes in CTCF binding pattern and increases in chromatin accessibility occur exclusively within 100 kbp of HPV integration sites. The chromatin changes co-occur with out-sized changes in transcription and alternative splicing of local genes. Analysis of The Cancer Genome Atlas (TCGA) HPV⁺ tumors indicates that HPV integration upregulates genes which have significantly higher essentiality scores compared to randomly selected upregulated genes from the same tumors.

CONCLUSIONS

Our results suggest that introduction of a new CTCF binding site due to HPV integration reorganizes chromatin state and upregulates genes essential for tumor viability in some HPV⁺ tumors. These findings emphasize a newly recognized role of HPV integration in oncogenesis.

A correlation graph attention network for classifying chromosomal instabilities from histopathology whole-slide images

The chromosome instability (CIN) is one of the hallmarks of cancer and is closely related to tumor metastasis. However, the sheer size and resolution of histopathology whole-slide images (WSIs) already challenges the capabilities of computational pathology. In this study, we propose a correlation graph attention network (MLP-GAT) that can construct graphs for classifying multi-type CINs from the WSIs of breast cancer. We construct a WSIs dataset of breast cancer from the Cancer Genome Atlas Breast Invasive Carcinoma (TCGA-BRCA). Extensive experiments show that MLP-GAT far outperforms accepted state-of-the-art methods and demonstrate the advantages of the constructed graph networks for analyzing WSI data. The visualization shows the difference among the tiles in a WSI. Furthermore, the generalization performance of the proposed method was verified on the stomach cancer. This study provides guidance for studying the relationship between CIN and cancer from the perspective of image phenotype.

Integrative identification of non-coding regulatory regions driving metastatic prostate cancer

Large-scale sequencing efforts of thousands of tumor samples have been undertaken to understand the mutational landscape of the coding genome. However, the vast majority of germline and somatic variants occur within non-coding portions of the genome. These genomic regions do not directly encode for specific proteins, but can play key roles in cancer progression, for example by driving aberrant gene expression control. Here, we designed an integrative computational and experimental framework to identify recurrently mutated non-coding regulatory regions that drive tumor progression. Application of this approach to whole-genome sequencing (WGS) data from a large cohort of metastatic castration-resistant prostate cancer (mCRPC) revealed a large set of recurrently mutated regions. We used (i) in silico prioritization of functional non-coding mutations, (ii) massively parallel reporter assays, and (iii) in vivo CRISPR-interference (CRISPRi) screens in xenografted mice to systematically identify and validate driver regulatory regions that drive mCRPC. We discovered that one of these enhancer regions, GH22I030351, acts on a bidirectional promoter to simultaneously modulate expression of U2-associated splicing factor SF3A1 and chromosomal protein CCDC157. We found that both SF3A1 and CCDC157 are promoters of tumor growth in xenograft models of prostate cancer. We nominated a number of transcription factors, including SOX6, to be responsible for higher expression of SF3A1 and CCDC157. Collectively, we have established and confirmed an integrative computational and experimental approach that enables the systematic detection of non-coding regulatory regions that drive the progression of human cancers.

Whole Exome-Trio Analysis Reveals Rare Variants Associated with Congenital Pouch Colon

Anorectal malformations (ARM) are individually common, but Congenital Pouch Colon (CPC) is a rare anorectal anomaly that causes a dilated pouch and communication with the genitourinary tract. In this work, we attempted to identify de novo heterozygous missense variants, and further discovered variants of unknown significance (VUS) which could provide insights into CPC manifestation. From whole exome sequencing (WES) performed earlier, the trio exomes were analyzed from those who were admitted to J.K. Lon Hospital, SMS Medical College, Jaipur, India, between 2011 and 2017. The proband exomes were compared with the unaffected sibling/family members, and we sought to ask whether any variants of significant interest were associated with the CPC manifestation. The WES data from a total of 64 samples including 16 affected neonates (11 male and 5 female) with their parents and unaffected siblings were used for the study. We examined the role of rare allelic variation associated with CPC in a 16 proband/parent trio family, comparing the mutations to those of their unaffected parents/siblings. We also performed RNA-Seq as a pilot to find whether or not the genes harboring these mutations were differentially expressed. Our study revealed extremely rare variants, viz., TAF1B, MUC5B and FRG1, which were further validated for disease-causing mutations associated with CPC, further closing the gaps of surgery by bringing intervention in therapies.

Benefiting from the intrinsic role of epigenetics to predict patterns of CTCF binding

Motivation

One of the most relevant mechanisms involved in the determination of chromatin structure is the formation of structural loops that are also related with the conservation of chromatin states. Many of these loops are stabilized by CCCTC-binding factor (CTCF) proteins at their base. Despite the relevance of chromatin structure and the key role of CTCF, the role of the epigenetic factors that are involved in the regulation of CTCF binding, and thus, in the formation of structural loops in the chromatin, is not thoroughly understood.

Results

Here we describe a CTCF binding predictor based on Random Forest that employs different epigenetic data and genomic features. Importantly, given the ability of Random Forests to determine the relevance of features for the prediction, our approach also shows how the different types of descriptors impact the binding of CTCF, confirming previous knowledge on the relevance of chromatin accessibility and DNA methylation, but demonstrating the effect of epigenetic modifications on the activity of CTCF. We compared our approach against other predictors and found improved performance in terms of areas under PR and ROC curves (PRAUC-ROCAUC), outperforming current state-of-the-art methods.

When 3D genome changes cause disease: the impact of structural variations in congenital disease and cancer

Large structural variations (SV) are a class of mutations that have long been known to cause a wide range of genetic diseases, from rare congenital disease to cancer. Many of these SVs do not directly disrupt disease-related genes and determining causal genotype-phenotype relationships has been challenging to disentangle in the past. This has started to change with our increased understanding of the 3D genome folding. The pathophysiologies of the different types of genetic diseases influence the type of SVs observed and their genetic consequences, and how these are connected to 3D genome folding. We propose guiding principles for interpreting disease-associated SVs based on our current understanding of 3D chromatin architecture and the gene-regulatory and physiological mechanisms disrupted in disease.

3D genome alterations and editing in pathology

The human genome is folded into a multi-level 3D structure that controls many nuclear functions including gene expression. Recently, alterations in 3D genome organization were associated with several genetic diseases and cancer. As a consequence, experimental approaches are now being developed to modify the global 3D genome organization and that of specific loci. Here, we discuss emerging experimental approaches of 3D genome editing that may prove useful in biomedicine.

Myoepithelial carcinoma of the parotid gland with a novel CTCF::NCOA2 fusion

We describe a case of a myoepithelial carcinoma of the superficial parotid gland in a 46-year-old male harboring a novel CTCF::NCOA2 gene fusion. To our knowledge, this novel gene fusion has not been described previously in myoepithelial carcinoma. A 46-year-old male patient presented with a mass involving the superficial left parotid gland with extension into the external auditory canal (EAC) and erosion of the conchal cartilage. Histologically, the neoplasm was composed of uniform spindled, epithelioid/ovoid cells arranged in cords and nests within hyalinized to myxoid stroma. On immunohistochemistry (IHC), the tumor cells demonstrated patchy and variable staining for low molecular weight cytokeratin (CAM5.2), pan-cytokeratin (OSCAR), and S-100. Overall, the morphological and immunohistochemical attributes supported a locally aggressive tumor of myoepithelial differentiation consistent with myoepithelial carcinoma. Molecular analysis using a custom 115-gene gene panel by targeted RNA sequencing, showed an in-frame CTCF::NCOA2 fusion. In addition to reporting this novel fusion in myoepithelial carcinoma, we also discuss relevant differential diagnosis, and provide a brief review of NCOA2 gene function in both normal and neoplastic contexts.

Gastric adenocarcinoma with high‑level microsatellite instability: A case report

Gastric cancer (GC) ranks fifth on the list of the most common malignancies worldwide. In Peru, gastric neoplasms are considered the second leading cause of mortality among males. Among the molecular subgroups of GC, microsatellite instability presents a favorable prognosis due to its hypermutated phenotype, which activates immunosurveillance. The present study describes the case of a 75-year-old patient, who was admitted in the hospital with a history of upper gastrointestinal bleeding and recurrent hospital admission, due to severe anemia. The patient presented with pale skin, normal vital functions, slight swelling of the lower extremities, and abdominal distention and bloating upon a physical examination. An endoscopic examination revealed an infiltrating circular ulcerated lesion. The histopathological analysis identified a moderately differentiated intestinal-type adenocarcinoma with pathological stage T3N0M0. Tumor genomic profiling demonstrated alterations in 15 different genes with a tumor mutational burden of 28 mutations/Mb. Finally, the patient underwent a partial gastrectomy without pre-operative chemotherapy. After 4 days, the patient presented with post-operative complications for which he was re-operated on. The patient did not survive. To the best of our knowledge, in the present case, pernicious anemia was an early sign of GC and a gastroscopy had to be performed. Furthermore, MutS homolog 3 alterations probably conditioned the presence of multiple frame-shift mutations.

Alterations of cohesin complex genes in acute myeloid leukemia: differential co-mutations, clinical presentation and impact on outcome

Functional perturbations of the cohesin complex with subsequent changes in chromatin structure and replication are reported in a multitude of cancers including acute myeloid leukemia (AML). Mutations of its STAG2 subunit may predict unfavorable risk as recognized by the 2022 European Leukemia Net recommendations, but the underlying evidence is limited by small sample sizes and conflicting observations regarding clinical outcomes, as well as scarce information on other cohesion complex subunits. We retrospectively analyzed data from a multi-center cohort of 1615 intensively treated AML patients and identified distinct co-mutational patters for mutations of STAG2, which were associated with normal karyotypes (NK) and concomitant mutations in IDH2, RUNX1, BCOR, ASXL1, and SRSF2. Mutated RAD21 was associated with NK, mutated EZH2, KRAS, CBL, and NPM1. Patients harboring mutated STAG2 were older and presented with decreased white blood cell, bone marrow and peripheral blood blast counts. Overall, neither mutated STAG2, RAD21, SMC1A nor SMC3 displayed any significant, independent effect on clinical outcomes defined as complete remission, event-free, relapse-free or overall survival. However, we found almost complete mutual exclusivity of genetic alterations of individual cohesin subunits. This mutual exclusivity may be the basis for therapeutic strategies via synthetic lethality in cohesin mutated AML.

Nuclear envelope, chromatin organizers, histones, and DNA: The many achilles heels exploited across cancers

In eukaryotic cells, the genome is organized in the form of chromatin composed of DNA and histones that organize and regulate gene expression. The dysregulation of chromatin remodeling, including the aberrant incorporation of histone variants and their consequent post-translational modifications, is prevalent across cancers. Additionally, nuclear envelope proteins are often deregulated in cancers, which impacts the 3D organization of the genome. Altered nuclear morphology, genome organization, and gene expression are defining features of cancers. With advances in single-cell sequencing, imaging technologies, and high-end data mining approaches, we are now at the forefront of designing appropriate small molecules to selectively inhibit the growth and proliferation of cancer cells in a genome- and epigenome-specific manner. Here, we review recent advances and the emerging significance of aberrations in nuclear envelope proteins, histone variants, and oncohistones in deregulating chromatin organization and gene expression in oncogenesis.

Accurate somatic variant detection using weakly supervised deep learning

Identification of somatic mutations in tumor samples is commonly based on statistical methods in combination with heuristic filters. Here we develop VarNet, an end-to-end deep learning approach for identification of somatic variants from aligned tumor and matched normal DNA reads. VarNet is trained using image representations of 4.6 million high-confidence somatic variants annotated in 356 tumor whole genomes. We benchmark VarNet across a range of publicly available datasets, demonstrating performance often exceeding current state-of-the-art methods. Overall, our results demonstrate how a scalable deep learning approach could augment and potentially supplant human engineered features and heuristic filters in somatic variant calling.

Multilevel view on chromatin architecture alterations in cancer

Chromosomes inside the nucleus are not located in the form of linear molecules. Instead, there is a complex multilevel genome folding that includes nucleosomes packaging, formation of chromatin loops, domains, compartments, and finally, chromosomal territories. Proper spatial organization play an essential role for the correct functioning of the genome, and is therefore dynamically changed during development or disease. Here we discuss how the organization of the cancer cell genome differs from the healthy genome at various levels. A better understanding of how malignization affects genome organization and long-range gene regulation will help to reveal the molecular mechanisms underlying cancer development and evolution.

Deep learning approaches for noncoding variant prioritization in neurodegenerative diseases

Determining how noncoding genetic variants contribute to neurodegenerative dementias is fundamental to understanding disease pathogenesis, improving patient prognostication, and developing new clinical treatments. Next generation sequencing technologies have produced vast amounts of genomic data on cell type-specific transcription factor binding, gene expression, and three-dimensional chromatin interactions, with the promise of providing key insights into the biological mechanisms underlying disease. However, this data is highly complex, making it challenging for researchers to interpret, assimilate, and dissect. To this end, deep learning has emerged as a powerful tool for genome analysis that can capture the intricate patterns and dependencies within these large datasets. In this review, we organize and discuss the many unique model architectures, development philosophies, and interpretation methods that have emerged in the last few years with a focus on using deep learning to predict the impact of genetic variants on disease pathogenesis. We highlight both broadly-applicable genomic deep learning methods that can be fine-tuned to disease-specific contexts as well as existing neurodegenerative disease research, with an emphasis on Alzheimer's-specific literature. We conclude with an overview of the future of the field at the intersection of neurodegeneration, genomics, and deep learning.

Untangling the roles of TOP2A and TOP2B in transcription and cancer

Type II topoisomerases (TOP2) are conserved regulators of chromatin topology that catalyze reversible DNA double-strand breaks (DSBs) and are essential for maintaining genomic integrity in diverse dynamic processes such as transcription, replication, and cell division. While controlled TOP2-mediated DSBs are an elegant solution to topological constraints of DNA, DSBs also contribute to the emergence of chromosomal translocations and mutations that drive cancer. The central importance of TOP2 enzymes as frontline chemotherapeutic targets is well known; however, their precise biological functions and impact in cancer development are still poorly understood. In this review, we provide an updated overview of TOP2A and TOP2B in the regulation of chromatin topology and transcription, and discuss the recent discoveries linking TOP2 activities with cancer pathogenesis.

ERR-activated GPR35 promotes immune infiltration level of macrophages in gastric cancer tissues

Enhancer release and retargeting (ERR) events could activate disease-causing gene promoters for increasing the expression level of oncogenes. Meanwhile, class A orphan GPCRs (oGPCRs) are known as potential biomarkers or drug targets for various cancers, such as gastric cancer (GC). Hence, systemic investigation of ERR events for class A oGPCRs in GC could help to explore biomarkers for GC. In this study, ENCODE and GTEx eQTL data were utilized to define ERR events in GC. Only GPR35 was then detected that could be activated by ERR in GC based on these data and ChIP-seq. Then, activated GPR35 functional in GC cells were explored by flow cytometry, cell-based wound healing assay, Transwell migration assay, and M2 polarization of macrophages assay. Meanwhile, according to TCGA and GEO database, overall survival, immune-related gene expression, and immune cell infiltration level in different GPR35 expressions were calculated. Here, we found ERR event activate GPR35 results in GC cells proliferation and migration, and partly immune cells significance exhaustion (CD8 + T-cells and CD4 + memory T-cells) and/or infiltration (T-cells and macrophage). Meanwhile, high GRP35 level leads to a poor prognosis in GC patients, probably partly due to it promoting the immune infiltration level of macrophages and then inducing polarization of M2 macrophages. Notably, GPR35’s high expression in CTSB+ and CD68 + macrophage could be a genetic indicator for early warning of primary GC. Hence, our findings provide a novel activation approach for oGPCRs, and GPR35 could be determined as a new drugable receptor and early genetic indicator for GC.

Phosphorylated histone variant γH2Av is associated with chromatin insulators in Drosophila

Chromatin insulators are responsible for orchestrating long-range interactions between enhancers and promoters throughout the genome and align with the boundaries of Topologically Associating Domains (TADs). Here, we demonstrate an association between gypsy insulator proteins and the phosphorylated histone variant H2Av (γH2Av), normally a marker of DNA double strand breaks. Gypsy insulator components colocalize with γH2Av throughout the genome, in polytene chromosomes and in diploid cells in which Chromatin IP data shows it is enriched at TAD boundaries. Mutation of insulator components su(Hw) and Cp190 results in a significant reduction in γH2Av levels in chromatin and phosphatase inhibition strengthens the association between insulator components and γH2Av and rescues γH2Av localization in insulator mutants. We also show that γH2Av, but not H2Av, is a component of insulator bodies, which are protein condensates that form during osmotic stress. Phosphatase activity is required for insulator body dissolution after stress recovery. Together, our results implicate the H2A variant with a novel mechanism of insulator function and boundary formation.

Somatic Mutational Landscape in Mexican Patients: CDH1 Mutations and chr20q13.33 Amplifications Are Associated with Diffuse-Type Gastric Adenocarcinoma

The Hispanic population, compared with other ethnic groups, presents a more aggressive gastric cancer phenotype with higher frequency of diffuse-type gastric adenocarcinoma (GA); this could be related to the mutational landscape of GA in these patients. Using whole-exome sequencing, we sought to present the mutational landscape of GA from 50 Mexican patients who were treated at The Instituto Nacional de Cancerología from 2019 to 2020. We performed a comprehensive statistical analysis to explore the relationship of the genomic variants and clinical data such as tumor histology and presence of signet-ring cell, H. pylori, and EBV. We describe a potentially different mutational landscape between diffuse and intestinal GA in Mexican patients. Patients with intestinal-type GA tended to present a higher frequency of NOTCH1 mutations, copy number gains in cytobands 13.14, 10q23.33, and 12q25.1, and copy number losses in cytobands 7p12, 14q24.2, and 11q13.1; whereas patients with diffuse-type GA tended to present a high frequency of CDH1 mutations and CNV gains in cytobands 20q13.33 and 22q11.21. This is the first description of a mutational landscape of GA in Mexican patients to better understand tumorigenesis in Hispanic patients and lay the groundwork for discovering potential biomarkers and therapeutic targets.

Zinc finger protein 280C contributes to colorectal tumorigenesis by maintaining epigenetic repression at H3K27me3-marked loci

This study uncovered the role of ZNF280C, a known DNA damage response protein, as a tumorigenic transcription regulator that contributes to colorectal tumorigenesis and metastasis through maintaining an epigenetic repression program at key cancer gene loci. These findings identified a contributor with potential prognostic value to colorectal pathogenesis and provide mechanistic insight to the essential function of transcription factor in fine-tuning the activity of chromatin regulators for proper transcription control.

Dysregulated epigenetic and transcriptional programming due to abnormalities of transcription factors (TFs) contributes to and sustains the oncogenicity of cancer cells. Here, we unveiled the role of zinc finger protein 280C (ZNF280C), a known DNA damage response protein, as a tumorigenic TF in colorectal cancer (CRC), required for colitis-associated carcinogenesis and Apc deficiency–driven intestinal tumorigenesis in mice. Consistently, ZNF280C silencing in human CRC cells inhibited proliferation, clonogenicity, migration, xenograft growth, and liver metastasis. As a C2H2 (Cys2-His2) zinc finger-containing TF, ZNF280C occupied genomic intervals with both transcriptionally active and repressive states and coincided with CCCTC-binding factor (CTCF) and cohesin binding. Notably, ZNF280C was crucial for the repression program of trimethylation of histone H3 at lysine 27 (H3K27me3)-marked genes and the maintenance of both focal and broad H3K27me3 levels. Mechanistically, ZNF280C counteracted CTCF/cohesin activities and condensed the chromatin environment at the cis elements of certain tumor suppressor genes marked by H3K27me3, at least partially through recruiting the epigenetic repressor structural maintenance of chromosomes flexible hinge domain-containing 1 (SMCHD1). In clinical relevance, ZNF280C was highly expressed in primary CRCs and distant metastases, and a higher ZNF280C level independently predicted worse prognosis of CRC patients. Thus, our study uncovered a contributor with good prognostic value to CRC pathogenesis and also elucidated the essence of DNA-binding TFs in orchestrating the epigenetic programming of gene regulation.

CTCF: a misguided jack-of-all-trades in cancer cells

The emergence and progression of cancers is accompanied by a dysregulation of transcriptional programs. The three-dimensional (3D) organization of the human genome has emerged as an important multi-level mediator of gene transcription and regulation. In cancer cells, this organization can be restructured, providing a framework for the deregulation of gene activity. The CTCF protein, initially identified as the product from a tumor suppressor gene, is a jack-of-all-trades for the formation of 3D genome organization in normal cells. Here, we summarize how CTCF is involved in the multi-level organization of the human genome and we discuss emerging insights into how perturbed CTCF function and DNA binding causes the activation of oncogenes in cancer cells, mostly through a process of enhancer hijacking. Moreover, we highlight non-canonical functions of CTCF that can be relevant for the emergence of cancers as well. Finally, we provide guidelines for the computational identification of perturbed CTCF binding and reorganized 3D genome structure in cancer cells.

Can aggressive cancers be identified by the "aggressiveness" of their chromatin?

Phenotypic plasticity is a crucial feature of aggressive cancer, providing the means for cancer progression. Stochastic changes in tumor cell transcriptional programs increase the chances of survival under any condition. I hypothesize that unstable chromatin permits stochastic transitions between transcriptional programs in aggressive cancers and supports non-genetic heterogeneity of tumor cells as a basis for their adaptability. I present a mechanistic model for unstable chromatin which includes destabilized nucleosomes, mobile chromatin fibers and random enhancer-promoter contacts, resulting in stochastic transcription. I suggest potential markers for "unsettled" chromatin in tumors associated with poor prognosis. Although many of the characteristics of unstable chromatin have been described, they were mostly used to explain changes in the transcription of individual genes. I discuss approaches to evaluate the role of unstable chromatin in non-genetic tumor cell heterogeneity and suggest using the degree of chromatin instability and transcriptional noise in tumor cells to predict cancer aggressiveness.

Pan-cancer analysis of mutations in open chromatin regions and their possible association with cancer pathogenesis

BACKGROUND

Open chromatin is associated with gene transcription. Previous studies have shown that the density of mutations in open chromatin regions is lower than that in flanking regions because of the higher accessibility of DNA repair machinery. However, in several cancer types, open chromatin regions show an increased local density of mutations in activated regulatory regions. Although the mutation distribution within open chromatin regions in cancer cells has been investigated, only few studies have focused on their functional implications in cancer. To reveal the impact of highly mutated open chromatin regions on cancer, we investigated the association between mutations in open chromatin regions and their possible functions.

METHODS

Whole-genome sequencing data of 18 cancer types were downloaded from the PanCancer Analysis of Whole Genomes and Catalog of Somatic Mutations in Cancer. We quantified the mutations located in open chromatin regions defined by The Cancer Genome Atlas and classified open chromatin regions into three categories based on the number of mutations. Then, we investigated the chromatin state, amplification, and possible target genes of the open chromatin regions with a high number of mutations. We also analyzed the association between the number of mutations in open chromatin regions and patient prognosis.

RESULTS

In some cancer types, the proportion of promoter or enhancer chromatin state in open chromatin regions with a high number of mutations was significantly higher than that in the regions with a low number of mutations. The possible target genes of open chromatin regions with a high number of mutations were more strongly associated with cancer than those of other open chromatin regions. Moreover, a high number of mutations in open chromatin regions was significantly associated with a poor prognosis in some cancer types.

CONCLUSIONS

These results suggest that highly mutated open chromatin regions play an important role in cancer pathogenesis and can be effectively used to predict patient prognosis.

Genome-wide analysis of somatic noncoding mutation patterns in cancer

We established a genome-wide compendium of somatic mutation events in 3949 whole cancer genomes representing 19 tumor types. Protein-coding events captured well-established drivers. Noncoding events near tissue-specific genes, such as ALB in the liver or KLK3 in the prostate, characterized localized passenger mutation patterns and may reflect tumor-cell-of-origin imprinting. Noncoding events in regulatory promoter and enhancer regions frequently involved cancer-relevant genes such as BCL6, FGFR2, RAD51B, SMC6, TERT, and XBP1 and represent possible drivers. Unlike most noncoding regulatory events, XBP1 mutations primarily accumulated outside the gene's promoter, and we validated their effect on gene expression using CRISPR-interference screening and luciferase reporter assays. Broadly, our study provides a blueprint for capturing mutation events across the entire genome to guide advances in biological discovery, therapies, and diagnostics.

High-resolution mapping demonstrates inhibition of DNA excision repair by transcription factors

DNA base damage arises frequently in living cells and needs to be removed by base excision repair (BER) to prevent mutagenesis and genome instability. Both the formation and repair of base damage occur in chromatin and are conceivably affected by DNA-binding proteins such as transcription factors (TFs). However, to what extent TF binding affects base damage distribution and BER in cells is unclear. Here, we used a genome-wide damage mapping method, N-methylpurine-sequencing (NMP-seq), and characterized alkylation damage distribution and BER at TF binding sites in yeast cells treated with the alkylating agent methyl methanesulfonate (MMS). Our data show that alkylation damage formation was mainly suppressed at the binding sites of yeast TFs ARS binding factor 1 (Abf1) and rDNA enhancer binding protein 1 (Reb1), but individual hotspots with elevated damage levels were also found. Additionally, Abf1 and Reb1 binding strongly inhibits BER in vivo and in vitro, causing slow repair both within the core motif and its adjacent DNA. Repair of ultraviolet (UV) damage by nucleotide excision repair (NER) was also inhibited by TF binding. Interestingly, TF binding inhibits a larger DNA region for NER relative to BER. The observed effects are caused by the TF–DNA interaction, because damage formation and BER can be restored by depletion of Abf1 or Reb1 protein from the nucleus. Thus, our data reveal that TF binding significantly modulates alkylation base damage formation and inhibits repair by the BER pathway. The interplay between base damage formation and BER may play an important role in affecting mutation frequency in gene regulatory regions.

Single-Cell Atlas of Lineage States, Tumor Microenvironment, and Subtype-Specific Expression Programs in Gastric Cancer

Gastric cancer heterogeneity represents a barrier to disease management. We generated a comprehensive single-cell atlas of gastric cancer (>200,000 cells) comprising 48 samples from 31 patients across clinical stages and histologic subtypes. We identified 34 distinct cell-lineage states including novel rare cell populations. Many lineage states exhibited distinct cancer-associated expression profiles, individually contributing to a combined tumor-wide molecular collage. We observed increased plasma cell proportions in diffuse-type tumors associated with epithelial-resident KLF2 and stage-wise accrual of cancer-associated fibroblast subpopulations marked by high INHBA and FAP coexpression. Single-cell comparisons between patient-derived organoids (PDO) and primary tumors highlighted inter- and intralineage similarities and differences, demarcating molecular boundaries of PDOs as experimental models. We complemented these findings by spatial transcriptomics, orthogonal validation in independent bulk RNA-sequencing cohorts, and functional demonstration using in vitro and in vivo models. Our results provide a high-resolution molecular resource of intra- and interpatient lineage states across distinct gastric cancer subtypes.

SIGNIFICANCE

We profiled gastric malignancies at single-cell resolution and identified increased plasma cell proportions as a novel feature of diffuse-type tumors. We also uncovered distinct cancer-associated fibroblast subtypes with INHBA-FAP-high cell populations as predictors of poor clinical prognosis. Our findings highlight potential origins of deregulated cell states in the gastric tumor ecosystem. This article is highlighted in the In This Issue feature, p. 587.

G-Quadruplex Matters in Tissue-Specific Tumorigenesis by BRCA1 Deficiency

How and why distinct genetic alterations, such as BRCA1 mutation, promote tumorigenesis in certain tissues, but not others, remain an important issue in cancer research. The underlying mechanisms may reveal tissue-specific therapeutic vulnerabilities. Although the roles of BRCA1, such as DNA damage repair and stalled fork stabilization, obviously contribute to tumor suppression, these ubiquitously important functions cannot explain tissue-specific tumorigenesis by BRCA1 mutations. Recent advances in our understanding of the cancer genome and fundamental cellular processes on DNA, such as transcription and DNA replication, have provided new insights regarding BRCA1-associated tumorigenesis, suggesting that G-quadruplex (G4) plays a critical role. In this review, we summarize the importance of G4 structures in mutagenesis of the cancer genome and cell type-specific gene regulation, and discuss a recently revealed molecular mechanism of G4/base excision repair (BER)-mediated transcriptional activation. The latter adequately explains the correlation between the accumulation of unresolved transcriptional regulatory G4s and multi-level genomic alterations observed in BRCA1-associated tumors. In summary, tissue-specific tumorigenesis by BRCA1 deficiency can be explained by cell type-specific levels of transcriptional regulatory G4s and the role of BRCA1 in resolving it. This mechanism would provide an integrated understanding of the initiation and development of BRCA1-associated tumors.

Multi-omic profiling of peritoneal metastases in gastric cancer identifies molecular subtypes and therapeutic vulnerabilities

Peritoneal metastasis, a hallmark of incurable advanced gastric cancer (GC), presently has no curative therapy and its molecular features have not been examined extensively. Here we present a comprehensive multi-omic analysis of malignant ascitic fluid samples and their corresponding tumor cell lines from 98 patients, including whole-genome sequencing, RNA sequencing, DNA methylation and enhancer landscape. We identify a higher frequency of receptor tyrosine kinase and mitogen-activated protein kinase pathway alterations compared to primary GC; moreover, approximately half of the gene alterations are potentially treatable with targeted therapy. Our analyses also stratify ascites-disseminated GC into two distinct molecular subtypes: one displaying active super enhancers (SEs) at the ELF3, KLF5 and EHF loci, and a second subtype bearing transforming growth factor-β (TGF-β) pathway activation through SMAD3 SE activation and high expression of transcriptional enhancer factor TEF-1 (TEAD1). In the TGF-β subtype, inhibition of the TEAD pathway circumvents therapy resistance, suggesting a potential molecular-guided therapeutic strategy for this subtype of intractable GC.

Mapping oncogenic protein interactions for precision medicine

Normal protein‐protein interactions (normPPIs) occur with high fidelity to regulate almost every physiological process. In cancer, this highly organised and precisely regulated network is disrupted, hijacked or reprogrammed resulting in oncogenic protein‐protein interactions (oncoPPIs). OncoPPIs, which can result from genomic alterations, are a hallmark of many types of cancers. Recent technological advances in the field of mass spectrometry (MS)‐based interactomics, structural biology and drug discovery have prompted scientists to identify and characterise oncoPPIs. Disruption of oncoPPI interfaces has become a major focus of drug discovery programs and has resulted in the use of PPI‐specific drugs clinically. However, due to several technical hurdles, studies to build a reference oncoPPI map for various cancer types have not been undertaken. Therefore, there is an urgent need for experimental workflows to overcome the existing challenges in studying oncoPPIs in various cancers and to build comprehensive reference maps. Here, we discuss the important hurdles for characterising oncoPPIs and propose a three‐phase multidisciplinary workflow to identify and characterise oncoPPIs. Systematic identification of cancer‐type‐specific oncogenic interactions will spur new opportunities for PPI‐focused drug discovery projects and precision medicine.

Mapping the genomic diaspora of gastric cancer

Gastric cancer (GC) is a leading contributor to global cancer incidence and mortality. Pioneering genomic studies, focusing largely on primary GCs, revealed driver alterations in genes such as ERBB2, FGFR2, TP53 and ARID1A as well as multiple molecular subtypes. However, clinical efforts targeting these alterations have produced variable results, hampered by complex co-alteration patterns in molecular profiles and intra-patient genomic heterogeneity. In this Review, we highlight foundational and translational advances in dissecting the genomic cartography of GC, including non-coding variants, epigenomic aberrations and transcriptomic alterations, and describe how these alterations interplay with environmental influences, germline factors and the tumour microenvironment. Mapping of these alterations over the GC life cycle in normal gastric tissues, metaplasia, primary carcinoma and distant metastasis will improve our understanding of biological mechanisms driving GC development and promoting cancer hallmarks. On the translational front, integrative genomic approaches are identifying diverse mechanisms of GC therapy resistance and emerging preclinical targets, enabled by technologies such as single-cell sequencing and liquid biopsies. Validating these insights will require specifically designed GC cohorts, converging multi-modal genomic data with longitudinal data on therapeutic challenges and patient outcomes. Genomic findings from these studies will facilitate 'next-generation' clinical initiatives in GC precision oncology and prevention.