- 1Basic Medical College, Shanxi Medical University, Taiyuan, Shanxi, China
- 2Department of Medical Laboratory, Jincheng People’s Hospital, Jincheng, Shanxi, China
- 3Laboratory Animal Center, Shanxi Province Cancer Hospital/Shanxi Hospital Affiliated to Cancer Hospital, Chinese Academy of Medical Sciences/Cancer Hospital Affiliated to Shanxi Medical University, Taiyuan, Shanxi, China
- 4Department of Tumor Biobank, Shanxi Province Cancer Hospital/Shanxi Hospital Affiliated to Cancer Hospital, Chinese Academy of Medical Sciences/Cancer Hospital Affiliated to Shanxi Medical University, Taiyuan, Shanxi, China
- 5Central Laboratory, Shanxi Province Cancer Hospital/Shanxi Hospital Affiliated to Cancer Hospital, Chinese Academy of Medical Sciences/Cancer Hospital Affiliated to Shanxi Medical University, Taiyuan, Shanxi, China
- 6Department of Laboratory Medicine, Shanxi Province Cancer Hospital/Shanxi Hospital Affiliated to Cancer Hospital, Chinese Academy of Medical Sciences/Cancer Hospital Affiliated to Shanxi Medical University, Taiyuan, Shanxi, China
Background: This study aimed to identify novel mutations associated with the progression of gastric cancer by establishing patient-derived xenograft (PDX) models and performing comprehensive genomic characterization of these PDX models and their corresponding primary tumors.
Methods: Fresh gastric cancer tissue samples were collected from 20 patients who underwent surgical resection at Shanxi Cancer Hospital and were subsequently implanted into NOD-SCID mice to establish PDX models. Histopathological features were evaluated using hematoxylin and eosin (H&E) staining. Whole-exome sequencing (WES) was performed on both primary tumors and their corresponding F1-PDX and F3-PDX tumors, focusing on mutations within 559 cancer-related genes. Predictive tools, including SIFT, Polyphen2_HVAR, Polyphen2_HDIV, and Mutation Taster, were utilized to identify potentially deleterious mutations, while I-Mutant and MUpro were employed to assess protein stability.
Results: Nine gastric cancer PDX models were successfully established, with seven models propagated to the third generation (F3-PDX), achieving an initial engraftment success rate of 45%. The latency of tumor establishment significantly decreased with each successive generation. The histological characteristics of the primary tumors were well preserved in the PDX models. WES of the three selected models revealed key mutated genes in primary tumors (F0), including IRS2, BLM, PDE4DIP, NUMA1, MYH9, TP53, PIK3CD, ERCC5, and ASXL1. A total of 28 somatic mutations were conserved across all three generations (F0, F1-PDX, and F3-PDX) in these models, representing a conservation rate of 43.75% (28/64). Among these conserved mutations, 10 were identified as potentially deleterious by multiple bioinformatics algorithms. Mutations in PTPRK (p.L988S), PIK3CB (p.F934L), LRP1B (p.A1912T), and IGF2R (p.G2052R) were predicted to significantly decrease protein stability.
Conclusion: This study demonstrated that PDX models effectively preserve the biological and genetic characteristics of primary gastric tumors, underscoring their utility in studying tumor heterogeneity. The integrated analysis of longitudinal WES data from primary tumors and matched PDXs enabled the identification of a core set of conserved, potentially deleterious mutations. The four prioritized mutations (PTPRK, PIK3CB, LRP1B, and IGF2R) provide new insights into the genetic landscape of gastric cancer and represent promising candidates for the development of targeted therapeutic strategies.
Highlights
• Successful establishment and longitudinal propagation of gastric cancer patient-derived xenograft (PDX) models in NOD-SCID mice, with histological features preserved across passages (F0 → F3).
• Integrated whole-exome sequencing (WES) analysis of matched primary tumors (F0) and PDXs (F1/F3) from three patients, revealing a core set of 28 somatic mutations conserved across all generations.
• Identification and prioritization of four novel, potentially deleterious missense mutations—PTPRK (p.L988S), PIK3CB (p.F934L), LRP1B (p.A1912T), and IGF2R (p.G2052R)—predicted to disrupt protein stability and key oncogenic pathways in gastric cancer.
Introduction
Gastric cancer remains one of the most common malignancies worldwide, with over one million new cases diagnosed annually. Despite advances in early detection and multimodal therapy, it continues to be a leading cause of cancer-related mortality, responsible for 768,793 deaths globally in 2020 (Sung et al., 2021; Chandarana et al., 2024). The clinical management of advanced disease is particularly challenging due to pronounced inter- and intra-tumoral heterogeneity and frequent resistance to conventional chemotherapy (Cho, 2020; Marin et al., 2016). Targeted therapies represent a promising strategy to overcome these limitations by exploiting tumor-specific molecular targets. To date, however, trastuzumab—approved for first-line treatment of HER2-positive advanced gastric cancer—is the only targeted agent with established clinical utility (Zhu et al., 2021). Its efficacy is often compromised by primary or rapidly acquired resistance, and broader molecularly guided approaches have yet to translate into substantial survival gains for most patients (Wang et al., 2022; Reddavid et al., 2021). These shortcomings highlight an urgent need to deepen our understanding of gastric cancer genomics and to identify novel, actionable driver alterations.
Patient-derived xenograft (PDX) models, generated by direct engraftment of tumor tissue into immunodeficient mice, have emerged as powerful preclinical tools that faithfully recapitulate the histopathological architecture, genomic landscape, and therapeutic responses of original patient tumors across serial passages (Meehan, 2019). In gastric cancer, PDX models have been successfully employed for drug screening and molecular profiling (Peille et al., 2020; Zhu et al., 2015). Nevertheless, their potential to uncover evolutionarily conserved, potentially deleterious somatic mutations—particularly those driving clonal adaptation during engraftment and propagation—remains underexploited. Whole-exome sequencing (WES) enables comprehensive characterization of coding-region mutations, facilitating the discovery of driver genes, elucidation of mutational processes, and identification of therapeutic targets (Nemtsova et al., 2020).
Here, we established a cohort of gastric cancer PDX models and performed WES on primary tumors (F0) and their corresponding early- (F1) and late-passage (F3) PDX derivatives. Our goal was to identify a core set of somatically mutated genes conserved across passages—thereby enriched for functionally relevant alterations—and to prioritize novel, potentially deleterious variants with implications for gastric cancer pathogenesis and targeted therapy development.
Materials and methods
Human gastric cancer tissue specimens
Fresh gastric cancer tissue samples were collected from 20 patients who underwent surgical resection at Shanxi Cancer Hospital in 2018. All procedures were conducted in accordance with approved protocols, and written informed consent was obtained from all participants. The study protocol and tissue collection procedures were reviewed and approved by the Medical Ethics Committee of Shanxi Cancer Hospital (Approval No.202009).
Immune-deficient mice
Female NOD-SCID (non-obese diabetic severe combined immunodeficiency) mice, aged 6–7 weeks, were purchased from Beijing Vital River Laboratory Animal Technology Co., Ltd. (Beijing, China). All animals were housed in a specific pathogen-free facility at the Laboratory Animal Center of Shanxi Cancer Hospital. All animal experiments were performed in compliance with institutional ethical guidelines and were approved by the Experimental Animal Welfare Ethics Review Committee of Shanxi Cancer Hospital (Approval No.202009).
Establishment of gastric cancer PDX models
Fresh gastric cancer tissue specimens (designated as F0) were collected immediately after surgical resection and immersed in ice-cold DMEM. Tissues were minced into fragments of approximately 2–3 mm3 and surgically implanted into the subcutaneous space of the right or left axilla of recipient NOD-SCID mice. Tumor fragments from each patient were engrafted into 3–4 recipient mice, which constituted the first passage (F1) cohort. Tumor volumes were measured weekly using digital caliper and calculated according to the formula: volume = 0.5× (width2 × length) (Faustino-Rocha et al., 2013). Successful engraftment was defined as the formation of a measurable tumor reaching a volume of ≥1,500 mm3 within 180 days post-implantation.
Upon reaching predefined endpoint (tumor volume ≥1,500 mm3), F1 tumors were harvested and sectioned into uniform fragments (about 2–3 mm3). A portion of each fragment was serially passaged into new recipient mice to generate subsequent generations (F2, F3, etc.), while the remainder was cryopreserved in liquid nitrogen for downstream molecular analyses (Figure 1). Differences in tumor latency among passages were assessed using the Kruskal–Wallis test with Dunn’s post hoc multiple comparisons test.
Figure 1. Workflow for whole-exome sequencing (WES) and bioinformatics analysis. The left panel illustrates the experimental pipeline: genomic DNA was extracted from tumor samples, quantified (≥0.5 μg), fragmented to 180–280 bp, and used for library preparation. Exome enrichment was performed using capture-based kits, followed by paired-end sequencing (PE150) on a next-generation sequencing platform. The right panel depicts the computational workflow: raw sequencing reads underwent quality control (QC) assessment and trimming to generate clean data, which were then aligned to the human reference genome (GRCh37). PCR duplicates were marked and removed, and variants—including single-nucleotide polymorphisms (SNPs) and insertions/deletions (INDELs)—were called and annotated. Functional impact of nonsynonymous variants was predicted using tools such as SIFT and PolyPhen-2, with potentially deleterious mutations identified based on combined prediction Scores.
Histological analysis
Primary tumors (F0), F1-PDX, and F3-PDX tissues were fixed in 4% paraformaldehyde for 24 h. Following fixation, tissues were processed, embedded in paraffin, and sectioned at 4 μm. Sections were stained with hematoxylin and eosin (H&E) and examined under a optical microscope. Histopathological concordance between primary tumors and matched PDXs was independently assessed by two board-certified pathologists blinded to sample identity, based on three criteria: 1) tissue architecture, 2) cellular morphology, and 3) nuclear atypia. Each feature was scored on a 4-point ordinal scale: 0, no similarity; 1, mild similarity; 2, moderate similarity; 3, near-identical appearance. Initial scoring was performed independently, and any discrepancies were resolved by consensus.
Whole-exome sequencing (WES) of 559 cancer-related genes
Among the successfully established PDX models, we selected a subset for longitudinal WES based on the following criteria: (i) robust engraftment and propagation to at least the F3 generation; (ii) high-quality DNA yield (>0.5μg, concentration ≥20 ng/μL) from both F0 and PDX tissues; and (iii) high histopathological fidelity to the original tumor.
Genomic DNA was extracted from frozen F0, F1-PDX, and F3-PDX tissues using a standard phenol-chloroform or commercial kit method, and quantified with a Qubit 2.0 fluorometer (Thermo Fisher Scientific). Samples with a DNA concentration of at least 20 ng/μL and a total yield of 0.5 μg or more were selected for library preparation. Genomic DNA was sheared to a target size of 180–280 bp using a Covaris M220 focused-ultrasonicator (Covaris, Inc.). Indexed sequencing libraries were prepared and enriched using the Agilent SureSelect XT Custom Target Enrichment Kit (Agilent Technologies) following the manufacturer’s instructions. Libraries were hybridized with biotinylated RNA baits designed to capture all exons of 559 cancer-associated genes (Supplementary Table S1). Magnetic beads coated with streptavidin were used to capture biotinylated probe–DNA hybrids. Captured libraries were PCR-amplified, and library quality was assessed using a Bioanalyzer (Agilent Technologies). Libraries were sequenced on an Illumina HiSeq PE150 platform using 150-bp paired-end reads (Illumina, Inc., San Diego, CA, USA). An average sequencing depth of ≥100× was achieved across all samples, with >99% of the target region covered at ≥20× (Figure 1).
Sequencing data analysis
Raw sequencing reads were assessed for quality using FastQC (v0.11.9) and trimmed to remove adapters and low-quality bases using Trim Galore (v0.6.6). For PDX samples, murine stromal reads were removed using Xenome (v1.0.1) with GRCh37 (human) and GRCm38 (mouse) reference genomes. High-quality reads were aligned to the human reference genome (hg19/GRCh37) using BWA-MEM (v0.7.17) (Li and Durbin, 2009). PCR duplicates were marked and removed using Sambamba (v0.7.1) (Tarasov et al., 2015). Somatic single-nucleotide variants (SNVs) and small insertions/deletions (indels) were called using GATK Mutect2 (v4.1.20) following the GATK Best Practices workflow.
A Panel of Normals (PON) was generated from 12 in-house non-tumor gastric mucosal samples and publicly available WES data from East Asian individuals, and used in GATK Mutect2 to filter out germline variants and technical artifacts. Only variants with a read depth ≥20× and VAF ≥5% in at least one sample were retained for downstream analysis. Variant calls were annotated using ANNOVAR (2019Oct24) and converted to Mutation Annotation Format (MAF). Cross-sample comparison of mutational profiles across F0, F1-PDX, and F3-PDX was performed using the maftools R package (v2.10.0) (Mayakonda et al., 2018).
Bioinformatics analysis
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed using the DAVID bioinformatics database (v6.8; https://david.ncifcrf.gov/). Enrichment significance was assessed using the hypergeometric test, and p-values were adjusted for multiple testing using the Benjamini–Hochberg procedure to control the false discovery rate (FDR). Pathways with an FDR <0.05 were considered significantly enriched.
Variant allele frequencies (VAF) across passages was visualized as a heatmap using the ggplot2 package (v3.4.0) in R (v4.3.1). For genes with multiple somatic variants, only the variant with the highest VAF in each sample was retained to represent the gene-level alteration, under the assumption that the dominant clone most likely drives functional consequences. Variants with VAF <0.1 were excluded from visualization to minimize the influence of potential sequencing artifacts. A red–white color scale was used (white: VAF = 0; dark red: VAF = 1). Samples were ordered chronologically (F0 → F1 → F3), and genes were sorted alphabetically. Plot parameters were optimized for clarity in publication.
The functional and clinical relevance of mutations was assessed by cross-referencing the dbSNP (build 155) and COSMIC (v99) databases. The functional impact of missense variants was predicted using SIFT (Kumar et al., 2009), Polyphen-2 HDIV and HVAR (Adzhubei et al., 2010), Mutation Taster (Schwarz et al., 2010), I-Mutant (Capriotti et al., 2005), and MUpro (Cheng et al., 2006).
Results
Establishment of gastric cancer patient-derived xenograft (PDX) models
Patient-derived xenograft (PDX) models of gastric cancer were established by subcutaneously implanting tumor fragments from 20 patients into 3–4 NOD-SCID mice per donor specimen. A total of 9 F1, 8 F2, and 7 F3 PDX lines were successfully established. Clinical characteristics of the 9 patients whose tumors yielded F1 PDX models are summarized in Supplementary Table S2.
Engraftment success rates were 45% (9/20) for F0-to-F1, 88.9% (8/9) for F1-to-F2, and 87.5% (7/8) for F2-to-F3 passages. Tumor growth kinetics accelerated significantly across passages, with median time to reach 1,500 mm3 decreasing from 28 to 65 days in F1 to 10–25 days in F3. These intergenerational differences were statistically significant overall (Kruskal–Wallis test, p < 0.001), with pairwise comparisons revealing significant reductions in latency between F1 and F2 (p = 0.003) and between F2 and F3 (p = 0.012) (Figure 2).
Figure 2. Flowchart of patient-derived xenograft (PDX) model establishment in gastric cancer. Primary tumors from 20 gastric cancer patients (F0) were implanted into female NOD-SCID mice (3–4 fragments per mouse). Engraftment success was defined as tumor volume ≥1,500 mm3 within 180 days. Of the initial 20 cases, 9 (45%) successfully engrafted at F1; 8 of these 9 (88.9%) propagated to F2, and 7 of 8 (87.5%) progressed to F3. The latency period—the time from implantation to tumor volume reaching ≥1,500 mm3—significantly decreased across passages: F1, 28–65 days; F2, 13–28 days (p = 0.003 vs. F1); F3, 10–25 days (p = 0.012 vs. F2), as determined by Dunn’s test with Benjamini–Hochberg correction for multiple comparisons. Three PDX lines were selected for whole-exome sequencing based on high DNA quality, sufficient yield, and histopathological fidelity (>80% tumor cellularity confirmed by H&E staining).
Histopathological assessment demonstrated high morphological concordance between primary tumors and their matched PDX models at both early (F1) and late (F3) passages. In both GC002 and GC006 PDX lines, key architectural and cytological features of the original tumors were preserved, with similarity scores consistently ranging from 2 to 3 across all three histological domains: tissue architecture, cellular morphology, and nuclear atypia (Figure 3; Supplementary Table S3). These results indicate that PDX models retain the histological features of their parental tumors across multiple passages, supporting their utility as preclinical models for translational research.
Figure 3. Hematoxylin and eosin (H&E) staining of primary F0 tumors and their corresponding F1-PDX and F3-PDX models from two representative cases (GC002 and GC006). Histological features—including irregular glandular architecture, branching luminal spaces (GC002), and marked nuclear atypia (enlarged, hyperchromatic nuclei with irregular contours)—are consistently preserved across generations. In GC002, the F3-PDX shows increased cellular density and reduced glandular differentiation, possibly reflecting clonal selection during serial passaging; however, nuclear morphology remains faithful to the original tumor. In GC006, tubular gland formation and nuclear pleomorphism are maintained throughout multiple passages. These findings demonstrate high morphological fidelity between primary tumors and PDX models, confirming that key histopathological characteristics are retained over successive engraftments. This preservation supports the use of these models in preclinical studies of gastric cancer biology and therapy. Scale bars: 100 μm (×20 magnification).
Somatic mutations in primary gastric tumors (F0)
Whole-exome sequencing (WES) of primary tumors (F0) identified 64 somatic mutations across 53 of the 559 targeted cancer-related genes. Mutation annotation revealed that missense variants constituted the predominant mutation type (Figures 4A,E), Single-nucleotide polymorphisms (SNPs) were the most frequent class (Figure 4B), and C>T transitions represented the most common substitution (Figure 4C). Among the cohort, GC005-F0 exhibited the highest mutational burden (Figure 4D).
Figure 4. Molecular characterization of somatic mutations in primary gastric tumors based on whole-exome sequencing of 559 cancer-related genes. (A) Variant classification distribution: missense mutations (green) are the most frequent, followed by nonsense (light red), in-frame insertions (dark red), and frameshift deletions (blue). Missense variants account for the majority of coding alterations, reflecting their common role in functional modulation. (B) Distribution of variant types: single-nucleotide polymorphisms (SNPs; purple) constitute ∼60% of all variants, with insertions (INS; yellow) and deletions (DEL; cyan) occurring less frequently. (C) Spectrum of single-nucleotide variants (SNVs): C>T transitions (red) are the most prevalent substitution type, consistent with known mutational signatures in gastric cancer, followed by T>C (yellow), C>G (blue), and T>A (green) changes. (D) Number of somatic variants per tumor sample: each bar represents an individual tumor (GC002-F0, GC006-F0, GC005-F0), with the dashed line indicating the median variant burden (18 variants/sample). (E) Summary of variant classifications across samples: box plot showing the distribution of somatic variant counts by type. Missense mutations exhibit the highest median count, while frameshift deletions and nonsense mutations are rare. (F) Top 10 most frequently mutated genes in primary tumors: horizontal bar chart displaying total somatic variant counts per gene (X-axis), with color-coded mutation types (red = frameshift, green = missense). Notable recurrently altered genes include IRS2, BLM, PDE4DP, NUMA1, and TP53. In panels (A), (B), (C), and (F), the X-axis represents somatic variant counts; in panel (D), the Y-axis shows variant counts per sample; in panel (E), the Y-axis indicates variant counts by classification.
Nine genes harbored two or more somatic mutations across the cohort: IRS2, BLM, PDE4DIP, NUMA1, MYH9, TP53, PIK3CD, ERCC5, and ASXL1. Notably, recurrent mutations in TP53, ERCC5, PIK3CD, BLM, and ASXL1 were observed in multiple tumors (Figure 4F; Supplementary Table S4).
In GC002-F0, NUMA1 carried two missense mutations (p.R1667C and p. A794G), and MYH9 harbored two missense variants (p.V615A and p. D7N). In GC005-F0, IRS2 harbored three distinct in-frame indels (p.H29_N30delinsNH, p. N23_T24delinsTN) and a missense variant (p.G11V), while PDE4DIP carried two missense mutations (p.P30S and p. P2355H) (Supplementary Table S4).
To investigate the biological relevance of the 53 mutated genes, KEGG pathway enrichment analysis was performed (Supplementary Table S5). Of these, 40 genes were mapped to 89 KEGG pathways. Notably, 13 of these genes were annotated to the PI3K–Akt signaling pathway: PDGFRA, CSF3R, SYK, FLT3, PIK3CD, PIK3CB, MTOR, PIK3CG, ERBB3, KIT, ERBB2, FGFR4, and TP53.
Somatic mutations in patient-derived xenografts (PDX)
The somatic mutation burden increased progressively across passages, with 64, 81, and 121 mutations identified in F0, F1, and F3, respectively. A mutation profile heatmap revealed that a core set of driver alterations was stably maintained across generations, whereas additional mutations emerged during serial passaging (Figure 5). Specifically, 28 genes—including ARID1B, ASXL1, ERBB3, IGF2R, IRS2, KIT, PIK3CB, and PTPRK—retained mutations across all three generations, reflecting the genomic stability of core driver events in PDX models. Notably, the variant allele frequencies (VAFs) of mutations in ARID1B, IRS2, PIK3CB, PTPRK, and TAL1 increased consistently from F0 to F3, suggestive of clonal expansion under in vivo selective pressure.
Figure 5. Heatmap of somatic mutation profiles across primary gastric tumors (F0) and their corresponding first- (F1) and third-generation (F3) PDX models. Each row represents a gene, and each column represents a sample (F0: primary tumor; F1/F3: first and third-generation PDX models). The color intensity reflects the VAF of somatic mutations, ranging from 0.00 (white) to 1.00 (dark red). Mutations are detected in 53 cancer-related genes across three representative gastric cancer cases (GC002, GC005, GC006).
Conversely, multiple F0-specific mutations—including those in ABCC4, CTNNA1, ERBB2, EPHB6, FAH, NOTCH4, NUP214, PIK3CG, RNF213, RNF43, ROS1, SDHA, SLC34A2, SYK, TET2, and TRIP11—were undetectable in later passages. This loss may reflect either true clonal extinction driven by selective pressures in the murine microenvironment or reduced detection sensitivity due to changes in tumor purity or sequencing coverage. De novo somatic mutations not present in F0 emerged in F1 and/or F3, including alterations in ARID1A, BCL9, BCOR, DICER1, PLAG1, and SNCAIP.
As shown in Table 1, ARID1A was wild-type in all F0 samples but acquired three mutations in F1 and six in F3. Similarly, ARID1B harbored a single mutation in F0 but expanded to 16 and 20 mutations in F1 and F3, respectively. Mapping of mutation sites revealed that most variants in both genes clustered within their N-terminal regions, which contain intrinsically disordered domains essential for protein–protein interactions (Supplementary Figure S1).
Prediction of potential deleterious mutations in gastric cancer
We identified 28 somatic mutations conserved between primary tumors and matched PDXs across three cases: GC002 (n = 9), GC005 (n = 14), and GC006 (n = 6). Among these, the ASXL1 p. G652S variant was shared by GC005 and GC006, whereas the remaining mutations were case-specific. The overall conservation rate of somatic mutations between F0 and PDX models was 43.8% (28/64).
The functional impact of these 28 conserved mutations was assessed using four in silico prediction algorithms: SIFT, PolyPhen-2 HDIV, PolyPhen-2 HVAR, and MutationTaster. SIFT predicted 13 mutations as deleterious, PolyPhen-2 (both HDIV and HVAR) classified 12 as probably damaging, and MutationTaster predicted 16 as disease-causing (Table 2). Ten mutations were concordantly predicted as deleterious by all four algorithms. The genomic positions of these ten mutations are depicted in a lollipop plot (Supplementary Figure S2).
Notably, four of these mutations—LRP1B p. A1912T, PDGFRA p. H425Y, NUMA1 p. R1667C, and UGT1A3 p. R45W—are catalogued in both dbSNP and COSMIC. Additionally, six novel deleterious mutations were identified in MTOR (p.K1606Q), NOTCH2 (p.C772Y), PTPRK (p.L988S), IGF2R (p.G2052R), ERBB3 (p.R1213W), and PIK3CB (p.F934L). To evaluate recurrence in gastric cancer, we queried the TCGA-STAD (Firehose Legacy) dataset (n = 441) via cBioPortal (https://www.cbioportal.org). Of these 10 genes, 8 (80%) were altered in ≥5% of TCGA-STAD cases, supporting their status as recurrently mutated genes in gastric adenocarcinoma (Supplementary Figure S3).
We further evaluated the impact of these ten mutations on protein stability using I-Mutant 2.0 and MUpro (Table 3). Both algorithms predicted decreased protein stability for nine of the ten mutations. Notably, PTPRK p. L988S, PIK3CB p. F934L, LRP1B p. A1912T, and IGF2R p. G2052R exhibited substantial destabilization, with ΔΔG values of −3.44 (I-Mutant)/−1.95 (MUpro), −2.70/−1.21, −1.41/−1.38, and −1.34/−0.76, respectively.
Discussion
Patient-derived xenograft (PDX) models have emerged as powerful platforms for recapitulating the biological and genomic fidelity of human tumors, offering significant advantages over conventional cell line–based systems in preclinical oncology research. In this study, we successfully established seven F3-PDX models from 20 surgically resected gastric cancer specimens, yielding an engraftment success rate of 35%—substantially higher than previously reported rates of 15.1% and 24.4% (Kuwata et al., 2019; Choi et al., 2016). This improved efficiency likely stems from optimized tumor acquisition protocols, selection of treatment-naïve specimens, and preferential engraftment of tumors from anatomically favorable primary sites, all of which are known to influence PDX take rates (Katsiampoura et al., 2017; Byrne et al., 2017). Furthermore, we observed a progressive shortening of tumor latency across successive passages (F0 → F1 → F3), consistent with prior observations in gastric cancer PDX models (Zhu et al., 2015; Yao et al., 2022), suggesting adaptation and clonal selection during serial propagation.
Despite their growing utility, most existing gastric cancer PDX resources lack comprehensive molecular annotation, limiting their application in precision oncology and biomarker discovery. To address this gap, we performed whole-exome sequencing (WES) on three primary tumors (F0) and their matched PDX passages (F1–F3), enabling longitudinal tracking of somatic alterations throughout engraftment and in vivo expansion. This approach allowed us to distinguish tumor-intrinsic, stably conserved mutations from acquired or context-dependent variants, thereby refining the identification of putative driver events.
Our WES analysis of the three F0 tumors revealed somatic alterations in IRS2, BLM, PDE4DIP, NUMA1, MYH9, TP53, PIK3CD, ERCC5, and ASXL1. Notably, this mutational profile only partially overlaps with large-scale genomic studies of gastric adenocarcinoma. For instance, a large-scale genomic study of 529 Chinese gastric adenocarcinomas highlighted TP53, ARID1A, LRP1B, PIK3CA, and CDH1 as top mutated genes (Yu et al., 2021). whereas another NGS-based analysis of 45 Chinese patients emphasized mutations in MLL4, ERBB3, FBXW7, MTOR, and NOTCH1 (Pan et al., 2018). The discrepancies observed here are almost certainly attributable to our extremely limited sample size (n = 3). With such a small cohort, even moderately recurrent drivers (e.g., ARID1A or CDH1, altered in >20% of cases in larger datasets) may be missed by chance, while rare or passenger variants unique to individual tumors may appear spuriously enriched. Moreover, given the well-established molecular heterogeneity of gastric cancer—including chromosomal instability (CIN), microsatellite instability (MSI), Epstein–Barr virus (EBV)-positive, and genomically stable (GS) subtypes—our three cases may represent a non-representative subset, further skewing the observed mutational landscape. Consequently, the specific gene alterations identified in our primary tumors should be interpreted as hypothesis-generating rather than reflective of population-level recurrence patterns.
Nonetheless, the longitudinal design of our study provides a unique advantage: by tracking variant allele frequencies (VAFs) across F0 → F1 → F3 passages, we operationally classified somatic mutations into conserved (present in all passages) and acquired (emerging or expanding during passaging) categories (Figure 5). Among the 28 conserved mutations, several reside in canonical cancer driver genes—including KIT, LRP1B, PDGFRA, ERBB3, MTOR, NOTCH2, PTPRD, PIK3CB, and PIK3CD—suggesting their functional importance in maintaining core oncogenic programs. In contrast, early subclonal mutations such as TP53, ERBB2, EPHA3, and SDHA were progressively lost during PDX propagation, implying that the corresponding subclones either failed to engraft or became dispensable under murine selective pressures (Ben-David et al., 2017).
Strikingly, we observed the emergence of late-acquired variants—including ARID1A/B, BCL9, BCOR, DICER1, KMT2A, and PLAG1—which were undetectable in F0 but reached VAFs ≥15% by F3. This rapid clonal expansion strongly suggests positive selection within the murine microenvironment. Notably, ARID1A and ARID1B mutations showed progressive enrichment across passages. As core subunits of the SWI/SNF (cBAF) chromatin remodeling complex, loss-of-function alterations in these genes disrupt chromatin architecture, impair DNA repair, and promote epigenetic plasticity (Krishnamurthy et al., 2022; Liang et al., 2024). Such changes may facilitate metabolic reprogramming and immune evasion, conferring a fitness advantage during xenograft adaptation (Dong et al., 2022; Chen et al., 2023; Caumanns et al., 2025). Importantly, concurrent ARID1A/ARID1B mutations have been reported in gastric, endometrial, and liver cancers (Kelso et al., 2017) and may act synergistically to drive tumorigenesis, positioning ARID1B as a potential therapeutic vulnerability in ARID1A-deficient contexts (Helming et al., 2014).
Intriguingly, the ARID1A and ARID1B mutations in our models clustered within their N-terminal intrinsically disordered regions (IDRs)—a domain critical for cBAF condensate formation, chromatin targeting, and protein–protein interactions (Patil et al., 2023). Although these variants were classified as “benign” by conventional in silico predictors (e.g., SIFT, PolyPhen-2)—which are primarily trained on structured domains—recent evidence indicates that IDR mutations can subtly alter biophysical properties (e.g., charge, hydrophobicity, linear motifs) to disrupt higher-order biomolecular organization without affecting global protein stability (Garcia et al., 2022). The recurrent localization of these variants within the IDR, coupled with their progressive VAF increase across passages, strongly argues against them being neutral passengers. Instead, they likely fine-tune cBAF complex dynamics to enhance epigenetic adaptability—a mechanism distinct from complete loss-of-function and potentially less susceptible to negative selection.
To prioritize functionally relevant conserved mutations, we integrated predictions from SIFT, PolyPhen-2, and MutationTaster, identifying ten high-confidence deleterious variants. Acknowledging the limitations of in silico tools—including training bias toward common variants and inter-algorithm discordance (Donner et al., 2021; Wang et al., 2024) —we adopted a conservative consensus approach to minimize false positives. Subsequent protein stability analysis using I-Mutant 2.0 and MUpro (Capriotti et al., 2005; Cheng et al., 2006) highlighted four mutations with pronounced destabilizing effects: PTPRK (p.L988S), PIK3CB (p.F934L), LRP1B (p.A1912T), and IGF2R (p.G2052R).
PTPRK, a candidate tumor suppressor encoding a receptor-type tyrosine phosphatase, harbors its L988S mutation within the catalytically active D1 domain. This alteration may impair dephosphorylation of oncogenic substrates, potentially acting as a gain-of-function driver. Similarly, the PIK3CB p. F934L mutation resides in the same kinase domain as the experimentally validated oncogenic E1051K variant (Whale et al., 2017), which enhances p110β catalytic activity and drives PI3K signaling independent of upstream activation. Given that p. F934L similarly localizes to this critical functional region, it represents a compelling candidate for isoform-selective PI3Kβ inhibition. LRP1B, one of the most frequently mutated tumor suppressors in gastric cancer (Wang et al., 2023), carries the A1912T substitution—a variant documented in COSMIC and dbSNP but lacking functional validation. Finally, the IGF2R G2052R mutation lies within the 14th extracellular repeat domain responsible for IGF-II binding (Brown et al., 2008); its presence may compromise ligand sequestration, thereby attenuating IGF2R’s tumor-suppressive function (Martin-Kleiner and Koraljka, 2010).
Collectively, these destabilizing mutations likely perturb key cellular processes—including signal transduction, proliferation, and apoptosis—and may serve as actionable targets. Of particular therapeutic interest is PIK3CB p. F934L, which localizes to the catalytic core and may confer resistance to pan-PI3K inhibitors while remaining sensitive to β-isoform–selective agents (e.g., GSK2636771 or AZD8186). To validate this hypothesis, we are currently generating isogenic cell lines and transgenic mouse models expressing this variant to assess drug sensitivity. Concurrently, we plan to use CRISPR interference (CRISPRi) to knock down mutant ARID1B in PDX-derived organoids, evaluating its impact on SWI/SNF complex integrity, chromatin accessibility, and transcriptional networks driving gastric cancer progression.
Finally, while our KEGG pathway enrichment analysis of the 53 somatic mutations in F0 tumors indicated significant enrichment in “Pathways in cancer” and the “PI3K-Akt signaling pathway” (Supplementary Table S4), these findings must be interpreted with caution. The extremely small input gene set (n = 53) and limited cohort size (n = 3) severely constrain statistical power and increase susceptibility to overrepresentation artifacts—particularly for broad, highly connected pathways like PI3K-Akt. Thus, these results should be regarded as exploratory and require validation in larger, molecularly stratified cohorts.
In summary, our study not only demonstrates the feasibility of establishing high-fidelity gastric cancer PDX models with enhanced engraftment efficiency but also provides a longitudinally annotated molecular resource that distinguishes stable driver events from microenvironment-driven adaptations. By integrating genomic tracking, in silico prioritization, and functional hypotheses, we lay the groundwork for future mechanistic studies and biomarker-driven therapeutic development in gastric cancer.
Conclusion
In conclusion, our study establishes a molecularly annotated panel of gastric cancer PDX models with an enhanced engraftment rate and demonstrates the power of longitudinal genomic tracking to distinguish tumor-intrinsic driver alterations from microenvironment-driven adaptive changes. By integrating exome sequencing across serial passages, functional prediction, and pathway analysis, we identify conserved deleterious mutations—including PIK3CB p. F934L and PTPRK p. L988S—as potential therapeutic targets, and reveal progressive selection of ARID1A/B mutations in intrinsically disordered regions that may fine-tune SWI/SNF complex function during xenograft adaptation. Although limited by a small cohort size, our findings provide a robust preclinical framework for validating genotype–phenotype relationships and advancing precision oncology strategies in gastric cancer.
Data availability statement
The original contributions presented in the study are publicly available. Raw sequence data reported in this study have been deposited in the Genome Sequence Archive (GSA) at the National Genomics Data Center, China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences, under accession numbers HRA016162 (GSA-Human) and CRA037299 (GSA). The datasets are publicly accessible at https://ngdc.cncb.ac.cn/gsa-human and https://ngdc.cncb.ac.cn/gsa, respectively.
Ethics statement
The studies involving humans were approved by the Medical Ethics Committee of Shanxi Cancer Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from a previous study for which ethical approval was obtained. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements. The animal study was approved by the Medical Ethics Committee of Shanxi Cancer Hospital. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
LK: Conceptualization, Data curation, Writing – original draft, Writing – review and editing. JW: Writing – original draft, Software. JZ: Software, Visualization, Writing – original draft. XY: Methodology, Writing – original draft. RS: Methodology, Writing – original draft. JK: Software, Visualization, Writing – original draft. YY: Software, Visualization, Writing – original draft. FL: Methodology, Writing – original draft. FW: Conceptualization, Formal Analysis, Writing – original draft. SG: Conceptualization, Funding acquisition, Writing – review and editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research was funded by the National Natural Science Foundation of China (81972840), the Science and Education Cultivation Fund of the National Cancer Regional Medical Center of Shanxi Provincial Cancer Hospital (BD2023005 and SD2023004), and Jincheng People’s Hospital Institutional Research Project (JSY-2025D004).
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2025.1571535/full#supplementary-material
References
Adzhubei, I. A., Schmidt, S., Peshkin, L., Ramensky, V. E., Gerasimova, A., Bork, P., et al. (2010). A method and server for predicting damaging missense mutations. Nat. Methods 7 (4), 248–249. doi:10.1038/nmeth0410-248
Ben-David, U., Ha, G., Tseng, Y.-Yi, Greenwald, N. F., Oh, C., Shih, J., et al. (2017). Patient-derived xenografts undergo Mouse-Specific tumor evolution. Nat. Genet. 49 (11), 1567–1575. doi:10.1038/ng.3967
Brown, J., Delaine, C., Zaccheo, O. J., Siebold, C., Gilbert, R. J., van Boxel, G., et al. (2008). Structure and functional analysis of the IGF-II/IGF2R interaction. EMBO J. 27 (1), 265–276. doi:10.1038/sj.emboj.7601938
Byrne, A. T., Alférez, D. G., Amant, F., Annibali, D., Arribas, J., Biankin, A. V., et al. (2017). Interrogating open issues in cancer precision medicine with patient-derived xenografts. Nat. Rev. Cancer 17 (4), 254–268. doi:10.1038/nrc.2016.140
Capriotti, E., Fariselli, P., and Casadio, R. (2005). I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 33 (Web Server), W306–W310. doi:10.1093/nar/gki375
Caumanns, J. J., Shang, Li, Meersma, G. J., Duiker, E. W., van der Zee, A. G. J., Wisman, G. B. A., et al. (2025). Establishment and characterization of ovarian clear cell carcinoma patient-derived xenografts. Sci. Rep. 15 (1), 4434. doi:10.1038/s41598-025-86384-2
Chandarana, C. V., Mithani, N. T., Singh, D. V., and Kikani, U. B. (2024). Vibrational spectrophotometry: a comprehensive review on the diagnosis of gastric and liver cancer. Curr. Pharm. Anal. 20 (7), 453–465. doi:10.2174/0115734129322567240821052326
Chen, X.-Yu, Li, B., Wang, Ye, Jin, J., Yang, Y., Huang, L. H., et al. (2023). Low level of ARID1A contributes to adaptive immune resistance and sensitizes triple-negative breast cancer to immune checkpoint inhibitors. Cancer Commun. 43 (9), 1003–1026. doi:10.1002/cac2.12465
Cheng, J., Randall, A., and Baldi, P. (2006). Prediction of protein stability changes for single-site mutations using support vector machines. Proteins Struct. Funct. Bioinforma. 62 (4), 1125–1132. doi:10.1002/prot.20810
Cho, S.-Y. (2020). Patient-Derived xenografts as compatible models for precision oncology. Laboratory Animal Res. 36 (1), 14. doi:10.1186/s42826-020-00045-1
Choi, Y. Y., Lee, J. E., Kim, H., Sim, M. H., Kim, K.-K., Lee, G., et al. (2016). Establishment and characterisation of patient-derived xenografts as paraclinical models for gastric cancer. Sci. Rep. 6 (1), 22172. doi:10.1038/srep22172
Dong, X., Song, S., Yuan, Li, Fan, Y., Wang, L., Wang, R., et al. (2022). Loss of ARID1A activates mTOR signaling and SOX9 in gastric Adenocarcinoma—Rationale for targeting ARID1A deficiency. Gut 71 (3), 467–478. doi:10.1136/gutjnl-2020-322660
Donner, I., Kiviluoto, T., Ristimaki, A., Aaltonen, L. A., and Vahteristo, P. (2021). Exome sequencing reveals three novel candidate predisposition genes for Diffuse gastric cancer.
Faustino-Rocha, A., Oliveira, P. A., Pinho-Oliveira, J., Teixeira-Guedes, C., Soares-Maia, R., da Costa, R. G., et al. (2013). Estimation of Rat mammary tumor volume using caliper and ultrasonography measurements. Lab. Anim. 42 (6), 217–224. doi:10.1038/laban.254
Garcia, F. A.De O., Andrade, E. S.De, and Palmero, E. I. (2022). Insights on variant analysis in silico tools for pathogenicity prediction. Front. Genet. 13 (November), 1010327. doi:10.3389/fgene.2022.1010327
Helming, K. C., Wang, X., Wilson, B. G., Vazquez, F., Haswell, J. R., Manchester, H. E., et al. (2014). ARID1B is a specific vulnerability in ARID1A-Mutant cancers. Nat. Med. 20 (3), 251–254. doi:10.1038/nm.3480
Katsiampoura, A., Raghav, K., Jiang, Z.-Q., Menter, D. G., Varkaris, A., Morelli, M. P., et al. (2017). Modeling of patient-derived xenografts in colorectal cancer. Mol. Cancer Ther. 16 (7), 1435–1442. doi:10.1158/1535-7163.MCT-16-0721
Kelso, T. W. R., Porter, D. K., Amaral, M. L., Shokhirev, M. N., Benner, C., and Hargreaves, D. C. (2017). Chromatin accessibility underlies synthetic lethality of SWI/SNF subunits in ARID1A-Mutant cancers. eLife 6 (October), e30506. doi:10.7554/eLife.30506
Krishnamurthy, N., Kato, S., Lippman, S., and Kurzrock, R. (2022). Chromatin remodeling (SWI/SNF) complexes, cancer, and response to immunotherapy. J. Immunother. Cancer 10 (9), e004669. doi:10.1136/jitc-2022-004669
Kumar, P., Henikoff, S., and Ng, P. C. (2009). Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 4 (7), 1073–1081. doi:10.1038/nprot.2009.86
Kuwata, T., Yanagihara, K., Iino, Y., Komatsu, T., Ochiai, A., Sekine, S., et al. (2019). Establishment of novel gastric cancer patient-derived xenografts and cell lines: pathological comparison between primary tumor, patient-Derived, and cell-line derived xenografts. Cells 8 (6), 585. doi:10.3390/cells8060585
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with burrows-wheeler transform. Bioinforma. Oxf. Engl. 25 (14), 1754–1760. doi:10.1093/bioinformatics/btp324
Liang, H., Zheng, X., Zhang, X., Zhang, Y., and Zheng, J. (2024). The role of SWI/SNF complexes in digestive System neoplasms. Med. Oncol. 41 (5), 119. doi:10.1007/s12032-024-02343-3
Marin, J. J. G., Al-Abdulla, R., Lozano, E., Briz, O., Bujanda, L., Banales, J. M., et al. (2016). Mechanisms of resistance to chemotherapy in gastric cancer. Anti-Cancer Agents Med. Chem. 16 (3), 318–334. doi:10.2174/1871520615666150803125121
Martin-Kleiner, I., and Koraljka, G. T. (2010). Mannose-6-Phosphate/Insulin-like growth factor 2 receptor (M6P/IGF2R) in carcinogenesis. Cancer Lett. 289 (1), 11–22. doi:10.1016/j.canlet.2009.06.036
Mayakonda, A., Lin, D.-C., Assenov, Y., Plass, C., and Phillip Koeffler, H. (2018). Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 28 (11), 1747–1756. doi:10.1101/gr.239244.118
Meehan, T. F. (2019). Know thy PDX model. Cancer Res. 79 (17), 4324–4325. doi:10.1158/0008-5472.CAN-19-2023
Nemtsova, M. V., Kalinkin, A. I., Kuznetsova, E. B., Bure, I. V., Alekseeva, E. A., Bykov, I. I., et al. (2020). Clinical relevance of somatic mutations in main driver genes detected in gastric cancer patients by next-generation DNA sequencing. Sci. Rep. 10 (1), 504. doi:10.1038/s41598-020-57544-3
Pan, X., Ji, X., Zhang, R., Zhou, Z., Zhong, Y., Peng, W., et al. (2018). Landscape of somatic mutations in gastric cancer assessed using next-generation sequencing analysis. Oncol. Lett. 16, 4863–4870. doi:10.3892/ol.2018.9314
Patil, A., Amy, R. S., Paulo, J. A., Collings, C. K., Ruff, K. M., Shinn, M. K., et al. (2023). A disordered Region controls cBAF activity via condensation and partner recruitment. Cell 186 (22), 4936–4955.e26. doi:10.1016/j.cell.2023.08.032
Peille, A.-L., Vincent, V., Wong, S.-S., Ting, J., Klingner, K., Zeitouni, B., et al. (2020). Evaluation of molecular subtypes and clonal selection during establishment of patient-derived tumor xenografts from Gastric adenocarcinoma. Commun. Biol. 3 (1), 367. doi:10.1038/s42003-020-1077-z
Reddavid, R., Dagatti, S., Franco, C., Puca, L., Tomatis, M., Corso, S., et al. (2021). Molecularly targeted therapies for gastric cancer. State of the art. Cancers 13 (16), 4094. doi:10.3390/cancers13164094
Schwarz, J. M., Rödelsperger, C., Schuelke, M., and Seelow, D. (2010). MutationTaster evaluates disease-causing potential of sequence alterations. Nat. Methods 7 (8), 575–576. doi:10.1038/nmeth0810-575
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA A Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Tarasov, A., Vilella, A. J., Cuppen, E., Nijman, I. J., and Prins, P. (2015). Sambamba: fast processing of NGS alignment formats. Bioinforma. Oxf. Engl. 31 (12), 2032–2034. doi:10.1093/bioinformatics/btv098
Wang, J., Huang, Q., Hu, X., Zhang, S., Jiang, Y., Yao, G., et al. (2022). Disrupting circadian rhythm via the PER1–HK2 axis reverses trastuzumab resistance in gastric cancer. Cancer Res. 82 (8), 1503–1517. doi:10.1158/0008-5472.CAN-21-1820
Wang, R., Zhang, G., Zhu, X., Xu, Y., Cao, N., Li, Z., et al. (2023). Prognostic implications of LRP1B and its relationship with the tumor-infiltrating immune cells in gastric cancer. Cancers 15 (24), 5759. doi:10.3390/cancers15245759
Wang, X., Ouyang, X., Zhang, D., Zhu, Y., Wu, L., Xiao, Z., et al. (2024). Two novel and three recurrent mutations in the mevalonate pathway genes in Chinese patients with porokeratosis. Clin. Cosmet. Investigational Dermatology 17 (January), 191–197. doi:10.2147/ccid.s444985
Whale, A. D., Colman, L., Lensun, L., Rogers, H. L., and Shuttleworth, S. J. (2017). Functional characterization of a novel somatic oncogenic mutation of PIK3CB. Signal Transduct. Target. Ther. 2 (1), 17063. doi:10.1038/sigtrans.2017.63
Yao, Y., Wang, Y., Chen, L., Tian, Z., Yang, G., Wang, R., et al. (2022). Clinical utility of PDX cohorts to reveal biomarkers of intrinsic resistance and clonal Architecture changes underlying acquired resistance to cetuximab in HNSCC. Signal Transduct. Target. Ther. 7 (1), 73. doi:10.1038/s41392-022-00908-0
Yu, P., Wang, Y., Yu, Y., Wang, A., Huang, L., Zhang, Y., et al. (2021). Deep targeted sequencing and its potential implication for cancer therapy in Chinese patients with gastric adenocarcinoma. Oncologist 26, e756–e768. doi:10.1002/onco.13695
Zhu, Y., Tian, T., Li, Z., Tang, Z., Wang, L., Wu, J., et al. (2015). Establishment and characterization of patient-derived tumor xenograft using gastroscopic biopsies in gastric cancer. Sci. Rep. 5 (1), 8542. doi:10.1038/srep08542
Keywords: bioinformatics tools, deleterious mutations, gastric cancer, patient-derived xenografts (PDX), whole-exome sequencing (WES)
Citation: Kong L, Wang J, Zheng J, Yang X, Sun R, Kou J, Yao Y, Li F, Wang F and Guo S (2026) Identification of potentially deleterious mutations in gastric cancer using patient-derived xenograft models. Front. Genet. 16:1571535. doi: 10.3389/fgene.2025.1571535
Received: 05 February 2025; Accepted: 17 December 2025;
Published: 29 January 2026.
Edited by:
Brian A. Boone, West Virginia University, United StatesReviewed by:
Md Salimullah, National Institute of Biotechnology (NIB), BangladeshXiao-Peng Duan, Shanghai Changzheng Hospital, China
Copyright © 2026 Kong, Wang, Zheng, Yang, Sun, Kou, Yao, Li, Wang and Guo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fuhua Wang, d2FuZ2Z1aHVhLjE2OEAxNjMuY29t; Sutang Guo, c3V0YW5nZ3VvQGFsaXl1bi5jb20=
†These authors have contributed equally to this work and share first authorship
Jie Wang1†