Genomic alterations in gastric cancers discovered via whole-exome sequencing

Mutational signatures in 175 Chinese gastric cancer patients

BACKGROUND

Gastric cancer (GC), a molecularly heterogeneous disease, is the third leading cause of cancer death worldwide. The majority of GC cases worldwide occur in East Asia, predominantly China. Mutational Signature Framework offers an elegant approach to identify mutational processes present in tumors.

METHODS

To identify mutational signature patterns, we conducted whole exome sequencing (WES) analysis in Chinese patients with GC. Mutect2 and MutsigCV were used to identify significantly mutated genes in 175 Chinese GC cases using paired tumor-normal tissues. We investigated mutational signatures using Catalogue of Somatic Mutations in Cancer (COSMIC) Version 2 (V2) and Version 3 (V3).

RESULTS

We identified 104 mutated genes with P < 0.01. Seven genes (OR6B1, B2M, ELF3, RHOA, RPL22, TP53, ARIDIA) had q < 0.0001, including six previously associated with GC. Mutational signatures (COSMIC-V3) observed include 14 single base substitutions (SBS), one doublet base substitution (DBS) Signature A, and one InDel (ID2). The most frequent SBS signatures (SBS05, SBS01, SBS15, SBS20, SBS40) were also observed in 254 White GC cases from The Cancer Genome Atlas (TCGA) Project. However, SBS01 and SBS20 showed significant differences between Whites vs. All Asians (19.3% vs. 11.3% for SBS 1 (P = 0.012) and 11.4% vs. 5.9% for SBS20 (P = 0.025), respectively). Using COSMIC V2, signatures 6, 15, and 1 were the most frequent in Chinese GC cases. Further, most Chinese GC cases carried multiple signatures.

CONCLUSIONS

This effort represents the most detailed mutational signatures analysis of GC cases from China to date. Results hold promise for new insights in understanding risk and prognosis factors in GC.

Bioinformatics-based analysis of the relationship between disulfidptosis and prognosis and treatment response in pancreatic cancer

Tumor formation is closely associated with disulfidptosis, a new form of cell death induced by disulfide stress-induced. The exact mechanism of action of disulfidptosis in pancreatic cancer (PCa) is not clear. This study analyzed the impact of disulfidptosis-related genes (DRGs) on the prognosis of PCa and identified clusters of DRGs, and based on this, a risk score (RS) signature was developed to assess the impact of RS on the prognosis, immune and chemotherapeutic response of PCa patients. Based on transcriptomic data and clinical information from PCa tissue and normal pancreatic tissue samples obtained from the TCGA and GTEx databases, differentially expressed and differentially surviving DRGs in PCa were identified from among 15 DRGs. Two DRGs clusters were identified by consensus clustering by merging the PCa samples in the GSE183795 dataset. Analysis of DRGs clusters about the PCa tumor microenvironment and differential analysis to obtain differential genes between the two DRG clusters. Patients were then randomized into the training and testing sets, and a prognostic prediction signature associated with disulfidptosis was constructed in the training set. Then all samples were divided into high-disulfidptosis-risk (HDR) and low-disulfidptosis-risk (LDR) subgroups based on the RS calculated from the signature. The predictive efficacy of the signature was assessed by survival analysis, nomograms, correlation analysis of clinicopathological characteristics, and the receiver operating characteristic (ROC) curves. To assess differences between different risk subgroups in immune cell infiltration, expression of immune checkpoint molecules, somatic gene mutations, and effectiveness of immunotherapy and chemotherapy. The GSE57495 dataset was used as external validation, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was used to detect the expression levels of DRGs. A total of 12 DRGs with differential expression and prognosis in PCa were identified, based on which a risk-prognosis signature containing five differentially expressed genes (DEGs) was developed. The signature was a good predictor and an independent risk factor. The nomogram and calibration curve shows the signature's excellent clinical applicability. Functional enrichment analysis showed that RS was associated with tumor and immune-related pathways. RS was strongly associated with the tumor microenvironment, and analysis of response to immunotherapy and chemotherapy suggests that the signature can be used to assess the sensitivity of treatments. External validation further demonstrated the model's efficacy in predicting the prognosis of PCa patients, with RT-qPCR and immunohistochemical maps visualizing the expression of each gene in PCa cell lines and the tissue. Our study is the first to apply the subtyping model of disulfidptosis to PCa and construct a signature based on the disulfidptosis subtype, which can provide an accurate assessment of prognosis, immunotherapy, and chemotherapy response in PCa patients, providing new targets and directions for the prognosis and treatment of PCa.

Analysis of genomic and immune intratumor heterogeneity in linitis plastica via multiregional exome and T-cell receptor sequencing

The molecular landscape and the intratumor heterogeneity (ITH) architecture of gastric linitis plastica (LP) are poorly understood. We performed whole-exome sequencing (WES) and T-cell receptor (TCR) sequencing on 40 tumor regions from four LP patients. The landscape and ITH at the genomic and immunological levels in LP tumors were compared with multiple cancers that have previously been reported. The lymphocyte infiltration was further assessed by immunohistochemistry (IHC) in LP tumors. In total, we identified 6339 non-silent mutations from multi-samples, with a median tumor mutation burden (TMB) of 3.30 mutations per Mb, comparable to gastric adenocarcinoma from the Cancer Genome Atlas (TCGA) cohort (P = 0.53). An extremely high level of genomic ITH was observed, with only 12.42%, 5.37%, 5.35%, and 30.67% of mutations detectable across 10 regions within the same tumors of each patient, respectively. TCR sequencing revealed that TCR clonality was substantially lower in LP than in multi-cancers. IHC using antibodies against CD4, CD8, and PD-L1 demonstrated scant T-cell infiltration in the four LP tumors. Furthermore, profound TCR ITH was observed in all LP tumors, with no T-cell clones shared across tumor regions in any of the patients, while over 94% of T-cell clones were restricted to individual tumor regions. The Morisita overlap index (MOI) ranged from 0.21 to 0.66 among multi-regions within the same tumors, significantly lower than that of lung cancer (P = 0.002). Our results show that LP harbored extremely high genomic and TCR ITH and suppressed T-cell infiltration, suggesting a potential contribution to the frequent recurrence and poor therapeutic response of this adenocarcinoma.

KLF4 Mutation Shapes Pathologic Characteristics of Foveolar-Type Gastric Adenoma in Helicobacter pylori-Naive Patients

Along with a recent remarkable decrease in Helicobacter pylori-infected individuals, reports of gastric neoplasms such as sporadic foveolar-type gastric adenoma (FGA) in H. pylori-naive patients have been increasing. This tumor, with its raspberry-like appearance, is common in H. pylori-naive gastric mucosa. The current study investigated the genomic features of sporadic FGA. Fresh-frozen sporadic FGA tissue samples from H. pylori-naive patients were subjected to whole genome analysis using a next-generation sequencer. Proliferation ability and apoptotic profiles of human gastric epithelial cells, along with plasmid transfection of candidate variants, were examined. A mean of 6.65 × 10⁸ total reads were obtained for each sample. Common genetic abnormalities in well-known proliferation driver genes of conventional gastric dysplasia/cancer were not found. However, a common single-nucleotide variation (SNV) was noted within the DNA-binding domain of the tumor suppressor gene KLF4. This novel SNV was located in the zinc finger 2 region. Additional experiments showed that it significantly suppressed proliferation of gastric epithelial cells compared with wild-type KLF4 plasmid-transfected cells, although suppression was reduced in early apoptotic phase-related genes. A novel SNV in the KLF4 zinc finger 2 region was commonly found in sporadic FGA tissue samples, which may explain the slow-growing properties of this neoplasm.

Application of next-generation sequencing in the diagnosis of gastric cancer

Objectives: Gastric cancer (GC) is a disease with high mortality, poor prognosis and numerous risk factors. GC has an asymptomatic nature in early stages of the diseases, making timely diagnosis complicated using common conventional approaches, namely pathological examinations and imaging tests. Recently, molecular profiling of GC using next generation sequencing (NGS) has opened new doors to efficient prognostic, diagnostic, and therapeutic strategies. The current review aims to thoroughly discuss and compare the current NGS techniques and commercial platforms utilized for GC diagnosis and treatment, highlighting the most recent NGS-based GC studies. Furthermore, this review addresses the challenges of clinical implementation of NGS in GC.Materials and methods: This review was conducted according to the eligible studies identified via search of Web of Science, PubMed, Scopus, Embase and the Cochrane Library. In the present study, data on gastric cancer patients and NGS methods used to diagnose the disease were reviewed.Conclusion: Given the ever-rising advancements in NGS technologies, bioinformatics, healthcare guidelines and refined classifications, it is hoped that these technologies can actualize their advantages and optimize GC patients' experience.

Cooperative participation of epigenomic and genomic alterations in the clinicopathological diversity of gastric adenocarcinomas: significance of cell adhesion and epithelial-mesenchymal transition-related signaling pathways

The present study was conducted to clarify the cooperative significance of epigenomic and genomic abnormalities during gastric carcinogenesis. Using 21 samples of normal control gastric mucosa (C), 109 samples of non-cancerous gastric mucosa (N) and 105 samples of cancerous tissue (T) from 109 patients with primary gastric adenocarcinomas, genome-wide DNA methylation analysis was performed using Infinium assay. Among these samples, 66 paired N and corresponding T samples were subjected to whole-exome and single nucleotide polymorphism array analyses. As had been shown in our previous study, 109 patients were clustered clinicopathologically into least aggressive Cluster A (n = 20), most aggressive Cluster B1 (n = 20) and Cluster B2 (n = 69). Most DNA methylation alterations in each cluster had already occurred even in N samples compared with C samples, and DNA methylation alterations at the precancerous N stage were inherited by the established cancers themselves. Recurrent single nucleotide variants and insertions/deletions resulting in functional disruption of the proteins encoded by the ABCA10, BNC2, CDH1, CTNNB1, SMAD4 and VAV2 genes were specific to Cluster B1, whereas those of the APC, EGFR, ERBB2, ERBB3, MLH1 and MUC6 genes were specific to Cluster A. MetaCore pathway analysis revealed that the epigenomically affected TWIST1 gene and genomically affected CDH1, CTNNB1, MMP9, TLN2, ROCK1 and SMAD4 genes were accumulated in signaling pathways related to cell adhesion, cytoskeleton remodeling and epithelial–mesenchymal transition in Cluster B1. These data indicate that epigenomic alterations at the precancerous stage are important in gastric carcinogenesis and that epigenomic and genomic alterations cooperatively underlie the aggressiveness of gastric adenocarcinomas.

Somatic alterations and mutational burden are potential predictive factors for metachronous development of early gastric cancer

The risk of developing metachronous gastric cancer (MGC) following curative endoscopic submucosal dissection (ESD) of early gastric cancer (EGC) remains even after eradicating Helicobacter pylori (HP) successfully. We screened initial EGC and adjacent non-cancerous mucosa ESD-resected specimens for somatic variants of 409 cancer-related genes, assessing their mutational burden (MB) to predict molecular markers for metachronous post-ESD development. We compared variants between ten patients diagnosed with MGC more than 3 years after ESD and ten age-matched patients who did not have MGC developments after successful HP eradication. We found no significant background differences between the two groups. In adjacent non-cancerous mucosa, the MB tended to be higher in the patients with metachronous developments than in the others. Somatic genomic alterations of RECQL4, JAK3, ARID1A, and MAGI1 genes were significantly associated with MGC development. The criteria including both the MB and their variants, which had potential significant values for predicting MGC. In conclusion, combined of assessing specific somatic variants and MB may be useful for predicting MGC development. This study included a limited number of subjects; however, our novel findings may encourage further exploration of the significance of the molecular features of EGC that predict MGC development, thereby promoting focused follow-up strategies and helping elucidate the mechanisms.

Deciphering the impact of missense mutations on structure and dynamics of SMAD4 protein involved in pathogenesis of gall bladder cancer

Gall bladder cancer (GBC) is the most common malignancy of biliary tract cancer associated with high mortality rate and poor prognosis due to lack of suitable biomarkers. In this study, we explored the structural and functional effects of different missense mutations occurs in SMAD4 that was associated with the development of GBC. We utilized in silico methods to predict the harmful effects of nonsynonymous missense mutations and monitored the stability of protein. We found that all mutations (D351N, G352E, R361C, R361H, E526Q) associated with SMAD4 were deleterious in nature resulting in the formation of deformed or unstable protein structure. Molecular dynamics simulation studies revealed how these mutations affect protein stability, structure, conformation and function. We observed, different mutants increase the compactness and rigidity of SMAD4 protein, alter secondary structure composition, decrease the surface area and protein-ligand interaction and affect its conformation. Findings of current work indicated that the analyzed mutations might affect the structure of protein and its caliber to interact with other molecules, which probably related to functional impairment of SMAD4 upon D351N, G352E, R361C, R361H, E526Q mutations and their involvement in cancer. Hence, the present study has significance of rational drug design and further increase our understanding of GBC development.Communicated by Ramaswamy H. Sarma.

Genome-Wide Study Updates in the International Genetics and Translational Research in Transplantation Network (iGeneTRAiN)

The prevalence of end-stage renal disease (ESRD) and the number of kidney transplants performed continues to rise every year, straining the procurement of deceased and living kidney allografts and health systems. Genome-wide genotyping and sequencing of diseased populations have uncovered genetic contributors in substantial proportions of ESRD patients. A number of these discoveries are beginning to be utilized in risk stratification and clinical management of patients. Specifically, genetics can provide insight into the primary cause of chronic kidney disease (CKD), the risk of progression to ESRD, and post-transplant outcomes, including various forms of allograft rejection. The International Genetics & Translational Research in Transplantation Network (iGeneTRAiN), is a multi-site consortium that encompasses >45 genetic studies with genome-wide genotyping from over 51,000 transplant samples, including genome-wide data from >30 kidney transplant cohorts (n = 28,015). iGeneTRAiN is statistically powered to capture both rare and common genetic contributions to ESRD and post-transplant outcomes. The primary cause of ESRD is often difficult to ascertain, especially where formal biopsy diagnosis is not performed, and is unavailable in ∼2% to >20% of kidney transplant recipients in iGeneTRAiN studies. We overview our current copy number variant (CNV) screening approaches from genome-wide genotyping datasets in iGeneTRAiN, in attempts to discover and validate genetic contributors to CKD and ESRD. Greater aggregation and analyses of well phenotyped patients with genome-wide datasets will undoubtedly yield insights into the underlying pathophysiological mechanisms of CKD, leading the way to improved diagnostic precision in nephrology.