Establishment and characterization of six human gastric carcinoma cell lines, including one naturally infected with Epstein-Barr virus

Epigenetic promoter alterations in GI tumour immune-editing and resistance to immune checkpoint inhibition

OBJECTIVES

Epigenomic alterations in cancer interact with the immune microenvironment to dictate tumour evolution and therapeutic response. We aimed to study the regulation of the tumour immune microenvironment through epigenetic alternate promoter use in gastric cancer and to expand our findings to other gastrointestinal tumours.

DESIGN

Alternate promoter burden (APB) was quantified using a novel bioinformatic algorithm (proActiv) to infer promoter activity from short-read RNA sequencing and samples categorised into APB_high, APB_int and APB_low. Single-cell RNA sequencing was performed to analyse the intratumour immune microenvironment. A humanised mouse cancer in vivo model was used to explore dynamic temporal interactions between tumour kinetics, alternate promoter usage and the human immune system. Multiple cohorts of gastrointestinal tumours treated with immunotherapy were assessed for correlation between APB and treatment outcomes.

RESULTS

APB_high gastric cancer tumours expressed decreased levels of T-cell cytolytic activity and exhibited signatures of immune depletion. Single-cell RNAsequencing analysis confirmed distinct immunological populations and lower T-cell proportions in APB_high tumours. Functional in vivo studies using 'humanised mice' harbouring an active human immune system revealed distinct temporal relationships between APB and tumour growth, with APB_high tumours having almost no human T-cell infiltration. Analysis of immunotherapy-treated patients with GI cancer confirmed resistance of APB_high tumours to immune checkpoint inhibition. APB_high gastric cancer exhibited significantly poorer progression-free survival compared with APB_low (median 55 days vs 121 days, HR 0.40, 95% CI 0.18 to 0.93, p=0.032).

CONCLUSION

These findings demonstrate an association between alternate promoter use and the tumour microenvironment, leading to immune evasion and immunotherapy resistance.

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.

Highly recurrent CBS epimutations in gastric cancer CpG island methylator phenotypes and inflammation

Background

CIMP (CpG island methylator phenotype) is an epigenetic molecular subtype, observed in multiple malignancies and associated with the epigenetic silencing of tumor suppressors. Currently, for most cancers including gastric cancer (GC), mechanisms underlying CIMP remain poorly understood. We sought to discover molecular contributors to CIMP in GC, by performing global DNA methylation, gene expression, and proteomics profiling across 14 gastric cell lines, followed by similar integrative analysis in 50 GC cell lines and 467 primary GCs.

Results

We identify the cystathionine beta-synthase enzyme (CBS) as a highly recurrent target of epigenetic silencing in CIMP GC. Likewise, we show that CBS epimutations are significantly associated with CIMP in various other cancers, occurring even in premalignant gastroesophageal conditions and longitudinally linked to clinical persistence. Of note, CRISPR deletion of CBS in normal gastric epithelial cells induces widespread DNA methylation changes that overlap with primary GC CIMP patterns. Reflecting its metabolic role as a gatekeeper interlinking the methionine and homocysteine cycles, CBS loss in vitro also causes reductions in the anti-inflammatory gasotransmitter hydrogen sulfide (H₂S), with concomitant increase in NF-κB activity. In a murine genetic model of CBS deficiency, preliminary data indicate upregulated immune-mediated transcriptional signatures in the stomach.

Conclusions

Our results implicate CBS as a bi-faceted modifier of aberrant DNA methylation and inflammation in GC and highlights H₂S donors as a potential new therapy for CBS-silenced lesions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02375-2.

Establishment and characterization of a new gastric cancer cell line, XGC-1

Background

To establish a primary human gastric cancer cell line.

Methods

Fresh gastric cancer tissue samples were separated into a cell suspension, and DMEM/F12 medium containing 10% foetal bovine serum was used for primary culture and subculture. The morphology of the cells was observed under a light microscope, and the cell growth curve was plotted. A soft agar colony formation assay was used to detect the colony formation ability of the cell line. Immunohistochemical methods were used to detect cytokeratin, vimentin and Ki-67, the chromosome G banding method was used to analyse the karyotype of the cells, and the tumourigenic ability of the cells was detected by subcutaneous inoculation of BALB/C nude mice.

Results

We established a gastric cancer cell line from a 68-year-old male patient. This gastric cancer cell line was named XGC-1 and had a doubling time of approximately 48 h. The cell line displayed strong colony formation ability and tumourigenicity in BALB/C nude mice and had complicated chromosomal abnormalities. When nutrients were insufficient, the cells shed and floated in the medium, but adherent growth was observed in nutrient-rich conditions.

Conclusions

The XGC-1 cell line will be useful for future studies of gastric cancer development, progression, metastasis and therapy.

Epstein-Barr Virus miR-BART17-5p Promotes Migration and Anchorage-Independent Growth by Targeting Kruppel-Like Factor 2 in Gastric Cancer

Epstein-Barr virus (EBV) infects more than 90% of the global population and is associated with a variety of tumors including nasopharyngeal carcinoma, Hodgkin lymphoma, natural killer/T lymphoma, and gastric carcinoma. In EBV-associated gastric cancer (EBVaGC), highly expressed EBV BamHI A rightward transcripts (BART) miRNAs may contribute to tumorigenesis with limited viral antigens. Despite previous studies on the targets of BART miRNAs, the functions of all 44 BART miRNAs have not been fully clarified. Here, we used RNA sequencing data from the Cancer Genome Atlas to find genes with decreased expression in EBVaGC. Furthermore, we used AGS cells infected with EBV to determine whether expression was reduced by BART miRNA. We showed that the expression of Kruppel-like factor 2 (KLF2) is lower in AGS-EBV cells than in the AGS control. Using bioinformatics analysis, four BART miRNAs were selected to check whether they suppress KLF2 expression. We found that only miR-BART17-5p directly down-regulated KLF2 and promoted gastric carcinoma cell migration and anchorage-independent growth. Our data suggest that KLF2 functions as a tumor suppressor in EBVaGC and that miR-BART17-5p may be a valuable target for effective EBVaGC treatment.

Heterogeneous Responses of Gastric Cancer Cell Lines to Tenovin-6 and Synergistic Effect with Chloroquine

Gastric cancer (GC) is the fifth most frequently diagnosed cancer and the third leading cause of cancer death. Approximately 15% of GC is associated with Epstein-Barr virus (EBV). GC is largely incurable with a dismal five-year survival rate. There is an urgent need to identify new therapeutic agents for the treatment of GC. Tenovin-6 was initially identified as a p53 activator, but it was later found to inhibit autophagy flux, and the protein deacetylase activity of sirtuins. Tenovin-6 shows promising therapeutic effect in various malignancies. However, it remains unknown whether Tenovin-6 is effective for GC. In this study, we found that EBV-positive and -negative GC cell lines were sensitive to Tenovin-6 but with different response times and doses. Tenovin-6 suppressed anchorage-independent growth of GC cells. Tenovin-6 induced different levels of apoptosis and phases of cell-cycle arrest depending on the cell lines with some manifesting gap 1 (G1) and others showing synthesis (S) phase cell-cycle arrest. Mechanistically, Tenovin-6 induced autophagy or p53 activation in GC cells depending on the status of TP53 gene. However, initiation of autophagy following treatment with Tenovin-6 conferred some protective effect on numerous cells. Combined treatment with Tenovin-6 and autophagy inhibitor chloroquine increased the cytotoxic effect by inducing microtubule-associated protein 1 light chain 3B (LC3B)-II accumulation, and by enhancing apoptosis and cell-cycle arrest. These results indicated that Tenovin-6 can be used as a potential therapeutic agent for GC, but the genetic background of the cancer cells might determine the response and mechanism of action. Treatment with Tenovin-6 alone or in combination with chloroquine could be a promising therapeutic approach for GC.

Screen technical noise in single cell RNA sequencing data

We proposed a data cleaning pipeline for single cell (SC) RNA-seq data, where we first screen genes (gene-wise screening) followed by screening cell libraries (library-wise screening). Gene-wise screening is based on the expectation that for a gene with a low technical noise, a gene's count in a library will tend to increase with the increase of library size, which was tested using negative binomial regression of gene count (as dependent variable) against library size (as independent variable). Library-wise screening is based on the expectation that across-library correlations for housekeeping (HK) genes is expected to be higher than the correlations for non-housekeeping (NHK) genes in those libraries with low technical noise. We removed those libraries, whose mean pairwise correlation for HK genes is NOT significantly higher than that for NHK genes. We successfully applied the pipeline to two large SC RNA-seq datasets. The pipeline was also developed into an R package.

EBV-associated gastric cancer evades T-cell immunity by PD-1/PD-L1 interactions

BACKGROUND

Epstein-Barr virus (EBV) is an oncogenic human herpesvirus involved in the development of around 10% of gastric cancers. The overexpression of PD-L1 is one of the features of EBV-associated gastric cancer (EBVaGC); however, the function of PD-L1 has not been studied in EBVaGC.

METHODS

We used three EBVaGC cell lines, SNU719 cells, NCC24 cells, and YCCEL1 cells, to evaluate the PD-L1 expression and function in EBVaGC. Jurkat T-lymphocytes expressing PD-1 were co-cultured with NCC24 and YCCEL1 cells and the cell cycles were analyzed. To study the regulatory mechanism for PD-L1 expression, the 3'UTR of PD-L1 was sequenced, and the effect of inhibitors of the IFN-γ signaling pathway was evaluated.

RESULTS

All of the EBVaGC cell lines expressed PD-L1, and its expression was further enhanced by stimulation with IFN-γ. In Jurkat T-cells co-cultured with IFN-γ-stimulated NCC24 and YCCEL1 cells, the number of cells in the G0/G1 phase was significantly increased. This G0/G1 arrest was partially released by administration of anti-PD-L1 antibody. We found SNPs in PD-L1 3'UTR nucleotide sequences that were located at seed regions for microRNAs. Treatment of EBVaGC cell lines with JAK2-inhibitor, PI3K-inhibitor, and mTOR inhibitor reduced the level of PD-L1 expression to the same level as cells without IFN-γ stimulation.

CONCLUSIONS

EBVaGC cells expressing high levels of PD-L1 suppress T-cell proliferation, and the IFN-γ signaling pathway is involved in the expression of PD-L1.

A Comprehensive Human Gastric Cancer Organoid Biobank Captures Tumor Subtype Heterogeneity and Enables Therapeutic Screening

Gastric cancer displays marked molecular heterogeneity with aggressive behavior and treatment resistance. Therefore, good in vitro models that encompass unique subtypes are urgently needed for precision medicine development. Here, we have established a primary gastric cancer organoid (GCO) biobank that comprises normal, dysplastic, cancer, and lymph node metastases (n = 63) from 34 patients, including detailed whole-exome and transcriptome analysis. The cohort encompasses most known molecular subtypes (including EBV, MSI, intestinal/CIN, and diffuse/GS, with CLDN18-ARHGAP6 or CTNND1-ARHGAP26 fusions or RHOA mutations), capturing regional heterogeneity and subclonal architecture, while their morphology, transcriptome, and genomic profiles remain closely similar to in vivo tumors, even after long-term culture. Large-scale drug screening revealed sensitivity to unexpected drugs that were recently approved or in clinical trials, including Napabucasin, Abemaciclib, and the ATR inhibitor VE-822. Overall, this new GCO biobank, with linked genomic data, provides a useful resource for studying both cancer cell biology and precision cancer therapy.

Induction of MAPK- and ROS-dependent autophagy and apoptosis in gastric carcinoma by combination of romidepsin and bortezomib

Proteasome inhibitors and histone deacetylase (HDAC) inhibitors can synergistically induce apoptotic cell death in certain cancer cell types but their combinatorial effect on the induction of autophagy remains unknown. Here, we investigated the combinatorial effects of a proteasome inhibitor, bortezomib, and an HDAC inhibitor, romidepsin, on the induction of apoptotic and autophagic cell death in gastric carcinoma (GC) cells. Isobologram analysis showed that low nanomolar concentrations of bortezomib/romidepsin could synergistically induce killing of GC cells. The synergistic killing was due to the summative effect of caspase-dependent intrinsic apoptosis and caspase-independent autophagy. The autophagic cell death was dependent on the activation of MAPK family members (ERK1/2 and JNK), and generation of reactive oxygen species (ROS), but was independent of Epstein-Barr virus infection. In vivo, bortezomib/romidepsin also significantly induced apoptosis and autophagy in GC xenografts in nude mice. This is the first report demonstrating the potent effect of combination of HDAC and proteasome inhibitors on the induction of MAPK- and ROS-dependent autophagy in addition to caspase-dependent apoptosis in a cancer type.

Establishment and characterization of an expanding-type gastric cancer cell line by Ming's classification

According to Ming's classification, gastric cancer (GC) can be divided into two types: expanding and infiltrative. The two types are readily recognizable by histology: expanding carcinomas grow en masse and by expansion, resulting in the formation of discrete tumour nodules, whereas in infiltrative carcinoma, tumour cells invade individually. Both types show varying degrees of cell maturation. The two types of carcinomas have vastly different pathological and clinical features. However, little is known concerning the mechanisms underlying these differences since no GC cell line models are available. For comprehensive and insightful analyses of mechanisms and treatment methods, new cell lines derived from expanding- and infiltrative-type gastric tumours are urgently needed. In the present study, we established an expanding-type GC cell line from a 72-year-old male patient. Different in vitro and in vivo methods were used to characterize the phenotypes of this cell line. This GC cell line was named XGC-2 and had an ~60 h doubling time. The cell line displayed strong colony formation and tumourigenicity in nude mice and had complicated chromosomal abnormalities. XGC-2 cells showed some markers of epithelial-to-mesenchymal transition (EMT), with decreased E-cadherin expression levels and increased vimentin expression levels. The XGC-2 cell line may be useful for future studies of GC development, progression, metastasis and therapy.

Improving gastric cancer preclinical studies using diverse in vitro and in vivo model systems

BACKGROUND

"Biomarker-driven targeted therapy," the practice of tailoring patients' treatment to the expression/activity levels of disease-specific genes/proteins, remains challenging. For example, while the anti-ERBB2 monoclonal antibody, trastuzumab, was first developed using well-characterized, diverse in vitro breast cancer models (and is now a standard adjuvant therapy for ERBB2-positive breast cancer patients), trastuzumab approval for ERBB2-positive gastric cancer was largely based on preclinical studies of a single cell line, NCI-N87. Ensuing clinical trials revealed only modest patient efficacy, and many ERBB2-positive gastric cancer (GC) patients failed to respond at all (i.e., were inherently recalcitrant), or succumbed to acquired resistance.

METHOD

To assess mechanisms underlying GC insensitivity to ERBB2 therapies, we established a diverse panel of GC cells, differing in ERBB2 expression levels, for comprehensive in vitro and in vivo characterization. For higher throughput assays of ERBB2 DNA and protein levels, we compared the concordance of various laboratory quantification methods, including those of in vitro and in vivo genetic anomalies (FISH and SISH) and xenograft protein expression (Western blot vs. IHC), of both cell and xenograft (tissue-sectioned) microarrays.

RESULTS

The biomarker assessment methods strongly agreed, as did correlation between RNA and protein expression. However, although ERBB2 genomic anomalies showed good in vitro vs. in vivo correlation, we observed striking differences in protein expression between cultured cells and mouse xenografts (even within the same GC cell type). Via our unique pathway analysis, we delineated a signaling network, in addition to specific pathways/biological processes, emanating from the ERBB2 signaling cascade, as a potential useful target of clinical treatment. Integrated analysis of public data from gastric tumors revealed frequent (10 - 20 %) amplification of the genes NFKBIE, PTK2, and PIK3CA, each of which resides in an ERBB2-derived subpathway network.

CONCLUSION

Our comprehensive bioinformatics analyses of highly heterogeneous cancer cells, combined with tumor "omics" profiles, can optimally characterize the expression patterns and activity of specific tumor biomarkers. Subsequent in vitro and in vivo validation, of specific disease biomarkers (using multiple methodologies), can improve prediction of patient stratification according to drug response or nonresponse.

Inhibition of class I histone deacetylases by romidepsin potently induces Epstein-Barr virus lytic cycle and mediates enhanced cell death with ganciclovir

Pan-histone deacetylase (HDAC) inhibitors, which inhibit 11 HDAC isoforms, are widely used to induce Epstein-Barr virus (EBV) lytic cycle in EBV-associated cancers in vitro and in clinical trials. Here, we hypothesized that inhibition of one or several specific HDAC isoforms by selective HDAC inhibitors could potently induce EBV lytic cycle in EBV-associated malignancies such as nasopharyngeal carcinoma (NPC) and gastric carcinoma (GC). We found that inhibition of class I HDACs, particularly HDAC-1, -2 and -3, was sufficient to induce EBV lytic cycle in NPC and GC cells in vitro and in vivo. Among a panel of selective HDAC inhibitors, the FDA-approved HDAC inhibitor romidepsin was found to be the most potent lytic inducer, which could activate EBV lytic cycle at ∼0.5 to 5 nM (versus ∼800 nM achievable concentration in patients' plasma) in more than 75% of cells. Upregulation of p21(WAF1) , which is negatively regulated by class I HDACs, was observed before the induction of EBV lytic cycle. The upregulation of p21(WAF1) and induction of lytic cycle were abrogated by a specific inhibitor of PKC-δ but not the inhibitors of PI3K, MEK, p38 MAPK, JNK or ATM pathways. Interestingly, inhibition of HDAC-1, -2 and -3 by romidepsin or shRNA knockdown could confer susceptibility of EBV-positive epithelial cells to the treatment with ganciclovir (GCV). In conclusion, we demonstrated that inhibition of class I HDACs by romidepsin could potently induce EBV lytic cycle and mediate enhanced cell death with GCV, suggesting potential application of romidepsin for the treatment of EBV-associated cancers.

Latent Expression of the Epstein-Barr Virus (EBV)-Encoded Major Histocompatibility Complex Class I TAP Inhibitor, BNLF2a, in EBV-Positive Gastric Carcinomas

The Epstein-Barr virus (EBV) BNLF2a gene product provides immune evasion properties to infected cells through inhibition of transporter associated with antigen processing (TAP)-mediated transport of antigen peptides. Although BNLF2a is considered to be a lytic gene, we demonstrate that it is expressed in nearly half of the EBV-associated gastric carcinomas analyzed. Further, we show that BNLF2a expression is dissociated from lytic gene expression. BNLF2a is therefore expressed in this latency setting, potentially helping protect the infected tumor cells from immunosurveillance.

Deregulation of immune response genes in patients with Epstein-Barr virus-associated gastric cancer and outcomes

BACKGROUND & AIMS

Patients with Epstein-Barr virus-associated gastric carcinoma (EBVaGC) have a better prognosis than those with gastric cancer not associated with EBV infection (EBVnGC). This is partly because EBV infection recruits lymphocytes, which infiltrate the tumor. A high degree of tumor heterogeneity is likely to be associated with poor response. We investigated differences in gene expression patterns between EBVaGC and EBVnGC.

METHODS

We used gene expression profile analysis to compare tumor and nontumor gastric tissues from 12 patients with EBVaGC and 14 patients with EBVnGC. Findings were validated by whole transcriptome RNAseq and real-time quantitative polymerase chain reaction analyses. CD3(+) primary T cells were isolated from human blood samples; migration of these cells and of Jurkat cells were measured in culture with EBV-infected and uninfected gastric cancer cells.

RESULTS

Based on Pearson correlation matrix analysis, EBVaGCs had a higher degree of homogeneity than EBVnGCs. Although 4550 genes were differentially expressed between tumor and nontumor gastric tissues of patients with EBVnGC, only 186 genes were differentially expressed between tumor and nontumor gastric tissues of patients with EBVaGC (P < .001). This finding supports the concept that EBVaGCs have fewer genetic and epigenetic alterations than EBVnGCs. Expression of major histocompatibility complex class II genes and genes that regulate chemokine activity were more often deregulated in EBVaGCs compared with nontumor tissues. In culture, more T cells migrated to EBV-infected gastric cancer cells than to uninfected cells; migration was blocked with a neutralizing antibody against CXCR3 (a receptor for many chemokines).

CONCLUSIONS

Fewer genes are deregulated in EBVaGC than in EBVnGC. Most changes in EBVaGCs occur in immune response genes. These changes might allow EBVaGC to recruit reactive immune cells; this might contribute to the better outcomes of these patients compared with those with EBVnGC.

Epstein-barr virus in gastric carcinoma

The Epstein-Barr virus (EBV) is detected in about 10% of gastric carcinoma cases throughout the world. In EBV-associated gastric carcinoma, all tumor cells harbor the clonal EBV genome. Gastric carcinoma associated with EBV has distinct clinicopathological features, occurs predominately in men and in younger-aged individuals, and presents a generally diffuse histological type. Most cases of EBV-associated gastric carcinoma exhibit a histology rich in lymphocyte infiltration. The immunological reactiveness in the host may represent a relatively preferable prognosis in EBV-positive cases. This fact highlights the important role of EBV in the development of EBV-associated gastric carcinoma. We have clearly proved direct infection of human gastric epithelialcells by EBV. The infection was achieved by using a recombinant EBV. Promotion of growth by EBV infection was observed in the cells. Considerable data suggest that EBV may directly contribute to the development of EBV-associated GC. This tumor-promoting effect seems to involve multiple mechanisms, because EBV affects several host proteins and pathways that normally promote apoptosis and regulate cell proliferation.

Establishing a peritoneal dissemination xenograft mouse model for survival outcome assessment of experimental gastric cancer

BACKGROUND

Peritoneal dissemination of gastric cancer is a common reason for unresectability, a frequent recurrence mechanism, and a common cause for death. The present study was performed to test peritoneal dissemination gastric cancer xenografts mouse models that would support survival outcome analyses.

MATERIALS AND METHODS

Human gastric cancer cell lines AGS, NCI-N87, and SNU-16 were intraperitoneally injected into nude mice and severe combined immunodeficiency (SCID) mice. The peritoneal tumor formation and mouse survival were compared among different groups. Mice were treated with oxaliplatin (5 mg/kg) and NVP-BEZ235 (10 mg/kg).

RESULTS

The formation rate of peritoneal cancer after intraperitoneal injection of 5 × 10(6) SNU16, NCI-N87, and AGS cells was 2/8, 6/8, and 0/8 in nude mice, and 6/6, 6/6, and 0/6 in SCID mice, respectively. Median animal survival with peritoneal dissemination was 74 d for NCI-N87 cells (10 × 10(6)), 95 d for SNU16 cells (10 × 10(6)), 78 d for SNU16 cells (20 × 10(6)), and 44 d for SNU16 cells (40 × 10(6)). In a therapeutic experiment with 40 × 10(6) SNU16 cells, animal survival was significantly improved by oxaliplatin treatment compared with the control group (58.5 d versus 45 d, P < 0.001), but not by NVP-BEZ235 (48 d versus 45 d, P = 0.249) treatment. In the accompanying subcutaneous SNU16 mouse model, relative tumor volume compared with controls was not significantly decreased by oxaliplatin treatment (P = 0.151) but by NVP-BEZ235 therapy (P = 0.008).

CONCLUSIONS

Peritoneal gastric cancer xenografts were successfully established after intraperitoneal injection NCI-N87 and SNU16 cells. These findings provide a useful survival outcome assessment model for experimental gastric cancer research.