A Four−Gene-Based Risk Score With High Prognostic Value in Gastric Cancer

Abstract

Background:

Gastric adenocarcinoma is an important contributor to cancer mortality and morbidity. This study aimed to explore the prognostic value of mutation patterns in gastric adenocarcinoma.

Materials and Methods:

We extracted somatic mutation data for 437 gastric adenocarcinoma samples from The Cancer Genome Atlas (TCGA) Stomach Adenocarcinoma (STAD) cohort. Kaplan–Meier survival in the R package maftools was used to analyze associations between mutations and survival. Multivariate Cox proportional model was used to establish risk formula. A four-gene-based risk score was developed to predict the overall survival of patients with gastric adenocarcinoma. We used the Tianjin cohort dataset with survival information to further evaluate the clinical value of this mutation signature.

Results:

Forty-five survival-related mutated genes were identified and verified, most of which were co-occurring in their mutation pattern and co-occurring with MLH3 and polymerase ϵ (POLE) mutations. Gastric adenocarcinoma samples with the 45 mutated genes had a significantly higher mutation count. Four-gene [UTRN, MUC16, coiled-coil domain-containing protein 178 (CCDC178), and HYDIN] mutation status was used to build a prognostic risk score that could be translated into the clinical setting. The association between the four-gene-based signature and overall survival remained statistically significant after controlling for age, sex, TNM stage, and POLE mutation status in the multivariate model [hazard ratio (HR), 1.88; 95% CI, 1.33–2.7; p < 0.001]. The prognostic significance of the four-gene-based risk score identified in TCGA cohort was validated in the Tianjin cohort.

Conclusion:

A four-mutated gene risk formula was developed that correlated with the overall survival of patients with gastric adenocarcinoma using a multivariable Cox regression model. In two independent genomic datasets from TCGA and Tianjin cohorts, low risk scores were associated with higher tumor mutation loads and improved outcome in patients with gastric adenocarcinoma. This finding may have implications for prognostic prediction and therapeutic guidance for gastric adenocarcinoma.

Introduction

Gastric adenocarcinoma is an important contributor to cancer mortality and morbidity, and its molecular mechanism remains largely incomprehensible. Next-generation sequencing (NGS) technology could provide genomic-level information about the mechanism of cancer. Numbers of large-scale genomic analyses on gastric adenocarcinoma have been completed, including The Cancer Genome Atlas (TCGA) project.

Recent research has shown that gastric adenocarcinoma is a heterogeneous disease. Surgical resection is still the main means of curative treatment for gastric adenocarcinoma. However, a portion of patients with advanced gastric adenocarcinoma developed local recurrences and distant metastases and had a poor prognosis after resection (1). Patients would have received the best treatment if their prognosis was depicted in advance. However, different prognoses of patients with similar clinical stages or pathologic grades remain unpredictable (2–5). Profiling the genetic mutation of gastric adenocarcinoma that influences the prognosis and accurate risk assessment based on genetic screening will lead to more effective clinical strategies in precision medicine.

In this study, we identified and verified 45 survival-related mutated genes with bioinformatics analysis from TCGA Stomach Adenocarcinoma (STAD) cohort. We investigated the function of these genes via Gene Ontology (GO) analysis. Through random survival forest algorithm, we ranked these mutated genes by importance and constructed a four−gene-based risk score with multivariable Cox regression model. Using the Tianjin cohort dataset with survival information, we evaluated the clinical value of the risk score.

Materials and Methods

Stomach Adenocarcinoma Datasets

Genomic data of gastric adenocarcinoma somatic mutation and gene expression data for 437 gastric adenocarcinoma samples in TCGA data portal (level 3) were downloaded from Genomic Data Commons (https://portal.gdc.cancer.gov) (Release Date: August 23, 2018). The Tianjin cohort contained data from 78 patients from northern China (6). Frozen tissue samples derived from surgical resection specimens of primary gastric adenocarcinoma from 294 northern Chinese patients without preoperative chemotherapy or radiotherapy were obtained from the Tianjin Medical University Cancer Institute and Hospital-National Foundation for Cancer Research Joint Tissue Banking Facility. Whole-exome sequencing was performed on 78 samples. Germline DNA was obtained from matching blood samples and used as a reference sequence to detect somatic mutations. Histopathologic diagnoses were independently reviewed by at least two experienced pathologists. Clinical follow-up data were complete for 78 participants with 25.08 months of median follow-up (32 deceased, 41.03%) (Supplementary Table 1). TCGA cohort also has the follow-up and vital status of patients (Supplementary Table 2). This study was approved by the Chinese PLA General Hospital (Beijing), which waived additional informed consent because all data used in this study were obtained from public databases. This study met the publication guidelines provided by TCGA (http://cancergenome.nih.gov/publications/publicationguidelines). All data were processed and analyzed by Excel 2010 and R (version 3.5.0).

Prognosis

Kaplan–Meier survival analyses implemented in the R package maftools and survival were used to analyze the correlation between mutations and survival (7). The log-rank test was used to determine significant differences of survival curves stratified by mutations. A two-sided p < 0.05 was considered statistically significant. Correlation between mutations and survival was also explored by multivariate Cox regression analyses by the R package survival.

Gene Ontology Pathway Analysis

The GO pathway analysis mutated genes were annotated by the R package of clusterProfiler (8). The cutoff p.adjust value was 0.01.

Random Survival Forest Algorithm

Random survival forest algorithm implemented in the R randomForestSRC was used to rank the survival-related genes by their importance.

Results

Survival in The Cancer Genome Atlas Stomach Adenocarcinoma Cohort

Somatic mutation data for 437 gastric adenocarcinoma samples from TCGA STAD cohort were extracted by maftools version 1.6.15 (http://www.bioconductor.org/packages/release/bioc/vignettes/maftools/inst/doc/maftools.html). These somatic mutations included point mutations and insertions/deletions (indels). We annotated 17,431 protein-coding genes with somatic mutations that mostly consist of missense mutations (Figure 1A), and single-nucleotide polymorphism (SNP) was more common than insertion–deletion (InDel) (Figure 1B). Furthermore, the mutational contexts are derived from combinations of six mutational types (i.e., T>G, T>C, T>A, C>T, C>G, and C>A) (Figure 1C). The somatic mutation rates varied considerably among the samples (Figures 1D, E), though an average of 108 mutations occurred in each sample (Figure 1D). The top 10 mutated genes were TTN, TP53, MUC16, LRP1B, SYNE1, ARID1A, CSMD3, FAT4, FLG, and PCLO (Figure 1F).

Figure 1

In order to select the most weighted genes associated with survival outcome, Kaplan–Meier survival in the R package maftools was used to analyze associations between mutations and survival for each of the 17,431 protein-coding genes (Figure 1B). A two-sided p < 0.05 was considered statistically significant. The most studied tumor repressor gene TP53 was analyzed for example. All samples were categorized into two groups representing the wild-type and mutated TP53 respectively, and Kaplan–Meier survival analysis showed that there were no significant associations between TP53 mutations and survival in TCGA STAD cohort (Figure 1G). Forty-four mutated genes, such as MUC16 and GLI3, occurring in more than 5% of the patients, were significantly associated with a better survival outcome (Figures 1H, I). Only coiled-coil domain-containing protein 178 (CCDC178) mutations were significantly associated with a poor survival outcome (Figure 1J). The 45 mutated genes were listed in Table 1.

To elucidate the function of these survival-related mutated genes, we conducted GO analysis and revealed that many genes play an important role in “cell–cell adhesion via plasma–membrane adhesion”, “extracellular matrix component”, and “alpha-catenin binding”, which were highly correlated to cancer metastasis and invasiveness (Table 2).

Clinical Features of Patients in The Cancer Genome Atlas Stomach Adenocarcinoma Cohort With the 45 Mutated Genes

To explore the relationship between these survival-related mutations, we performed pair-wise Fisher’s exact test to detect significant pairs of mutated genes in the 45 genes and DNA mismatch repair (MMR)-related genes [PMS2, MSH2, MLH1, MSH3, MLH3, MSH6, polymerase ϵ (POLE)]. Interestingly, most of the 45 mutated genes were co-occurring in their mutation pattern and co-occurring with MLH3 and POLE mutations (Figure 2A). Tumor mutation burden (TMB) is an important determinant for molecular subtyping of gastric adenocarcinoma in TCGA (9). Recent studies have shown that gastric adenocarcinoma with POLE mutations or microsatellite instability–high (MSI-H) had DNA MMR signatures and higher TMBs (10). Gastric adenocarcinoma samples with the 45 mutated genes had a significantly higher tumor count (Figure 2B; Mann–Whitney test p-value <0.0001). Since these mutations tended to occur simultaneously, we explored the relationship between all these mutations and prognosis. Kaplan–Meier survival analysis showed that the 45 mutated genes were significantly associated with a better survival outcome in TCGA STAD cohort (Figure 2C; log-rank test, p < 0.0001). Multivariate Cox regression analyses showed that the correlation remained statistically significant after controlling for confounding factors such as sex, age, and TNM stage (Figure 2D). So, the 45 mutated genes were important prognostic indicators associated with tumor mutation load.

Figure 2

Construction of Risk Score Formula

The most common mutations of the 45 genes are missense substitutions. Other alterations include silent mutations, frameshift insertions and deletions, nonsense mutations, and other infrequent mutations. Both base substitution mutation and frameshift mutation can change the composition or sequence of amino acids in the polypeptide chain. According to the impact of different type mutations on DNA composition or sequence of amino acids, we scored these mutations, as follows: No mutation or Silent mutation: 0; Missense mutation, In-frame insertion or deletion: +1; Splice site mutation: +2; Nonsense mutation, frameshift insertion or deletion: +3; Multiple mutation: the maximum score of all (Figure 3A). There are four levels of mutations in total. Finally, we get the mutation score matrix (Supplementary Table 3). In order to select the most weighted genes, we used random survival forest algorithm (Ntree = 1,000, default parameters of Hemant Ishwaran algorithm) and set the 45 mutated genes score as variables in this model. We ranked these 45 genes by their importance from the random survival forest model (Figure 3B).

Figure 3

Since the mutation status is divided into four levels, we selected four genes to construct the model in order to preserve the effect of all mutation states on the prognosis and take into account the minimalist principle. Combining Kaplan–Meier survival analysis log-rank p-value, mutation coexistence pattern, and the result of random survival forest algorithm, we selected UTRN, MUC16, CCDC178, and HYDIN as candidates for an accurate prediction of survival in gastric adenocarcinoma patients. UTRN ranked first among the candidates in random survival forest model (Figure 3B). MUC16 had the highest mutation frequency among the 45 genes (Table 1). CCDC178 mutations seem to be an independent factor (Figure 2A). MUC16 and CCDC178 also ranked in the top 10 important candidates (Figure 3B). UTRN and CCDC178 mutations were not significantly co-occurring or exclusive with each other in their mutation pattern (Figure 3C). Based on the selection of the three genes, using a computer to build the multivariate Cox model repeatedly, HYDIN was selected among the genes that had no obvious co-occurring or exclusive mutations with UTRN and CCDC178.

Of the 437 patients in TCGA cohort, the four genes were altered in 177 patients (40.5%). Consistently, gastric adenocarcinoma samples with the four mutated genes had a significantly higher mutation load (Figures 3D, E; Mann–Whitney test, p < 0.0001).

The mutated genes chosen from the previous step were constructed into the multivariate Cox proportional model to calculate the coefficients in TCGA cohort, thereby establishing the risk formula by which a risk score for each patient was calculated. Risk score = −0.1445* (mutation score of MUC16) − 0.459* (mutation score of UTRN) - 0.332* (mutation score of HYDIN) + 0.3102* (mutation score of CCDC178). Cutting off by 0, we defined risk score <0 as low-risk group, risk score = 0 as medium-risk group, and risk score >0 as high-risk group. Patients in the low-risk group had a markedly longer overall survival than those in the medium-risk group, and high-risk group had the shortest overall survival (Figure 3F; p < 0.0001, by log rank). The correlation between the four-gene-based signature and overall survival remained statistically significant after controlling for age, sex, TNM stage, and POLE mutation status in the multivariate model [hazard ratio (HR), 1.88; 95% CI, 1.33–2.7; p < 0.001] (Figure 3G). A significantly higher mutation count was also observed in gastric adenocarcinoma samples within the low-risk group (Figure 3H; Mann–Whitney test p-value <0.0001). So, the signature of the four-gene mutation would*nbsp;be a good prediction for survival of gastric adenocarcinoma patients.

Independent Validation of Four−Gene-Based Risk Score in the Tianjin Cohort

To further evaluate the clinical value of this four-gene mutation signature, we used Tianjin cohort dataset with survival information. Kaplan–Meier survival analyses showed low-risk scores were significantly associated with better survival outcomes (Figure 4A; log-rank test, p = 0.036). Significantly higher mutation count was also observed in Tianjin cohort gastric adenocarcinoma samples in the low-risk group (Figure 4B; Mann–Whitney test p-value <0.0001). Multivariable Cox regression analysis also showed that the association of the risk score with overall survival was statistically significant after controlling for age, sex, TNM stage, and POLE mutation status (Figure 4C; HR, 4.33; 95% CI, 1.29–14.5; p = 0.018).

Figure 4

Discussion

We analyzed 437 gastric adenocarcinoma samples from TCGA cohort and 78 gastric adenocarcinoma samples from a Tianjin cohort for survival prediction genes. We have identified and verified 45 mutated genes related to survival from 17,431 protein-coding genes in STAD. The GO enrichment showed that these genes play an important role in “cell–cell adhesion via plasma–membrane adhesion”, “extracellular matrix component”, and “alpha-catenin binding”. Obviously, these genes played a pivotal role in cancer metastasis. After that, we ranked the 45 survival−related mutated genes by random survival forest algorithm. Whole-genome sequencing is expensive, and data analysis takes a long time. The aim of this study is to develop a cheap and practical prognostic tool that can be accomplished by PCR. Combining Kaplan–Meier survival analysis log-rank p-value, mutation coexistence pattern, and the result of random survival forest algorithm, we selected four-gene mutation status to build a prognostic risk score that could be transformed into the clinical setting. Gastric adenocarcinoma samples classified into low-risk group had a significantly higher tumor mutation load and better survival outcome. The association between the four−gene-based risk score and overall survival was independent of mutations in POLE mutation status, age, sex, and TNM stage. Many studies have reported the association between these genes and tumors.

The deletion of chromosome 6q has been extensively mapped in a variety of tumors (11, 12). UTRN is located in this region, which encodes dystrophin. UTRN is a tumor suppressor inducing cell transformation when expressed in an antisense orientation. Studies showed decreased expression and inactivation mutations of UTRN in tumors. Expression of a wild-type UTRN in breast cancer cells inhibited tumor cell growth in vitro and reduced their tumor potential in nude mice (13). HYDIN is a gene whose impaired function has been linked to abnormal ciliary function, dyskinesia, and brain abnormalities (14, 15). HYDIN-derived sequences are targeted by the adaptive immunity in patients with cancer (16). Somatic mutations in HYDIN were found in breast cancer samples (17–19). MUC16, encoding a type I transmembrane mucin protein (20, 21), is frequently mutated in multiple types of human cancer (22). MUC16 was reported to modulate immune response to cancer (23–25). CA125 is a repeating peptide epitope of the mucin MUC16 (26, 27). MUC16 mutations were found to be associated with higher tumor mutation load, better survival outcomes, and immune response and cell cycle pathways in gastric adenocarcinoma (28). The CCDC178 is an 867-amino acid polypeptide and belongs to the superfamily of coiled-coil domain-containing protein. CCDC178 was reported to be mutated in hepatocellular carcinoma (29) and gastric carcinoma (30). CCDC178 associated with BRCA1-associated protein 2 (BRAP2), a negative regulator of extracellular signal-regulated kinase (ERK) pathway, and promoted its degradation (31). CCDC178 deficiency impaired the ERK activation in hepatocellular carcinoma (31). In our study, CCDC178 mutations were significantly associated with a poor survival outcome in gastric adenocarcinoma. The relative transcriptional level of CCDC178 was significantly downregulated in several types of carcinoma compared with adjacent non-cancerous tissues in TCGA cohorts (Supplementary Figure 1).

Compared with mutation detection, measurement of gene expression in cross-platform is unstable. Due to the lack of reproducibility and standardization, its clinical application may be limited. Recently, high-throughput sequencing technologies have been widely utilized in clinical cancer research. Compared to normal tissues, many high/low expressed proteins and mutated genes in tumor cells were identified. Combining several altered genes together may be feasible in predicting gastric adenocarcinoma risk and prognosis. In our study, the four genes (UTRN, MUC16, CCDC178, and HYDIN) were mutated in 177 patients (40.5%) in TCGA STAD cohort. The risk score is powerful and accurate in prognostic stratification. Our work provided an advanced method toward clinical applications of gene mutation profiling in STAD, especially in future personalized prediction and precision medicine.

However, our study has several limitations. The number of samples with follow-up data in the Tianjin cohort was limited. No CCDC178 mutation was detected in the Tianjin cohort. Gastric adenocarcinoma is a heterogeneous disease. Molecular subtyping can encompass this heterogeneity and provide useful clinical information. Prognostic tool constructed on the basis of anatomic site, histopathology, and molecular subtype may be more powerful and accurate. Considering the number of samples, anatomic site, histopathology, and molecular subtype were not included in this study. So, more prospective studies are necessary to further validate the reliability and stability of this risk score.

Conclusions

A four-mutated gene risk formula was developed that correlated with the overall survival of patients with gastric adenocarcinoma using a multivariable Cox regression model. In two independent genomic datasets from TCGA and Tianjin cohorts, low risk scores were associated with higher tumor mutation loads and improved outcome in patients with gastric adenocarcinoma. This finding may have implications for prognostic prediction and therapeutic guidance for gastric adenocarcinoma.

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.

References

• 1

  SongZWuYYangJYangDFangX. Progress in the Treatment of Advanced Gastric Cancer. Tumour Biol: J Int Soc Oncodevelopmental Biol Med (2017) 39(7):1010428317714626. doi: 10.1177/1010428317714626

  • CrossRef
  • Google Scholar

• 2

  GayedDTWodeyarJWangZXWeiXYaoYYChenXXet al. Prognostic Values of Inhibitory kappaB Kinases mRNA Expression in Human Gastric Cancer. Biosci Rep (2019) 39(1):BSR20180617. doi: 10.1042/BSR20180617

  • CrossRef
  • Google Scholar

• 3

  JinCLinJRMaLSongYShiYXJiangPet al. Elevated Spondin-2 Expression Correlates With Progression and Prognosis in Gastric Cancer. Oncotarget (2017) 8(6):10416–24. doi: 10.18632/oncotarget.14423

  • CrossRef
  • Google Scholar

• 4

  ErichsenRKelshMAOlinerKSNielsenKBFroslevTLaenkholmAVet al. Prognostic Impact of Tumor MET Expression Among Patients With Stage IV Gastric Cancer: A Danish Cohort Study. Ann Epidemiol (2016) 26(7):500–3. doi: 10.1016/j.annepidem.2016.05.002

  • CrossRef
  • Google Scholar

• 5

  LiuXXuDHuangCGuoYWangSZhuCet al. Regulatory T Cells and M2 Macrophages Present Diverse Prognostic Value in Gastric Cancer Patients With Different Clinicopathologic Characteristics and Chemotherapy Strategies. J Trans Med (2019) 17(1):192. doi: 10.1186/s12967-019-1929-9

  • CrossRef
  • Google Scholar

• 6

  ChenKYangDLiXSunBSongFCaoWet al. Mutational Landscape of Gastric Adenocarcinoma in Chinese: Implications for Prognosis and Therapy. Proc Natl Acad Sci USA (2015) 112(4):1107–12. doi: 10.1073/pnas.1422640112

  • CrossRef
  • Google Scholar

• 7

  MayakondaALinDCAssenovYPlassCKoefflerHP. Maftools: Efficient and Comprehensive Analysis of Somatic Variants in Cancer. Genome Res (2018) 28(11):1747–56. doi: 10.1101/gr.239244.118

  • CrossRef
  • Google Scholar

• 8

  YuGWangLGHanYHeQY. Clusterprofiler: An R Package for Comparing Biological Themes Among Gene Clusters. Omics: J Integr Biol May (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118

  • CrossRef
  • Google Scholar

• 9

  Cancer Genome Atlas ResearchN.Comprehensive Molecular Characterization of Gastric Adenocarcinoma. Nature (2014) 513(7517):202–9. doi: 10.1038/nature13480

  • CrossRef
  • Google Scholar

• 10

  LiXWuWKXingRWongSHLiuYFangXet al. Distinct Subtypes of Gastric Cancer Defined by Molecular Characterization Include Novel Mutational Signatures With Prognostic Capability. Cancer Res (2016) 76(7):1724–32. doi: 10.1158/0008-5472.CAN-15-2443

  • CrossRef
  • Google Scholar

• 11

  MertensFJohanssonBHoglundMMitelmanF. Chromosomal Imbalance Maps of Malignant Solid Tumors: A Cytogenetic Survey of 3185 Neoplasms. Cancer Res (1997) 57(13):2765–80.

  • Google Scholar

• 12

  OffitKParsaNZGaidanoGFilippaDALouieDPanDet al. 6q Deletions Define Distinct Clinico-Pathologic Subsets of Non-Hodgkin’s Lymphoma. Blood (1993) 82(7):2157–62. doi: 10.1182/blood.V82.7.2157.2157

  • CrossRef
  • Google Scholar

• 13

  LiYHuangJZhaoYLHeJWangWDaviesKEet al. UTRN on Chromosome 6q24 Is Mutated in Multiple Tumors. Oncogene (2007) 26(42):6220–8. doi: 10.1038/sj.onc.1210432

  • CrossRef
  • Google Scholar

• 14

  Brunetti-PierriNBergJSScagliaFBelmontJBacinoCASahooTet al. Recurrent Reciprocal 1q21.1 Deletions and Duplications Associated With Microcephaly or Macrocephaly and Developmental and Behavioral Abnormalities. Nat Genet (2008) 40(12):1466–71. doi: 10.1038/ng.279

  • CrossRef
  • Google Scholar

• 15

  OlbrichHSchmidtsMWernerCOnoufriadisALogesNTRaidtJet al. Recessive HYDIN Mutations Cause Primary Ciliary Dyskinesia Without Randomization of Left-Right Body Asymmetry. Am J Hum Genet (2012) 91(4):672–84. doi: 10.1016/j.ajhg.2012.08.016

  • CrossRef
  • Google Scholar

• 16

  LaskeKShebzukhovYVGrosse-HovestLKuprashDVKhlgatianSVKorolevaEPet al. Alternative Variants of Human HYDIN Are Novel Cancer-Associated Antigens Recognized by Adaptive Immunity. Cancer Immunol Res (2013) 1(3):190–200. doi: 10.1158/2326-6066.CIR-13-0079

  • CrossRef
  • Google Scholar

• 17

  StephensPJTarpeyPSDaviesHVan LooPGreenmanCWedgeDCet al. The Landscape of Cancer Genes and Mutational Processes in Breast Cancer. Nature (2012) 486(7403):400–4. doi: 10.1038/nature11017

  • CrossRef
  • Google Scholar

• 18

  BanerjiSCibulskisKRangel-EscarenoCBrownKKCarterSLFrederickAMet al. Sequence Analysis of Mutations and Translocations Across Breast Cancer Subtypes. Nature (2012) 486(7403):405–9. doi: 10.1038/nature11154

  • CrossRef
  • Google Scholar

• 19

  ZhangYCaiQShuXOGaoYTLiCZhengWet al. Whole-Exome Sequencing Identifies Novel Somatic Mutations in Chinese Breast Cancer Patients. J Mol Genet Med Int J Biomed Res (2015) 9(4):183. doi: 10.4172/1747-0862.1000183

  • CrossRef
  • Google Scholar

• 20

  XiangXFengMFelderMConnorJPManYGPatankarMSet al. HN125: A Novel Immunoadhesin Targeting MUC16 With Potential for Cancer Therapy. J Cancer (2011) 2:280–91. doi: 10.7150/jca.2.280

  • CrossRef
  • Google Scholar

• 21

  FelderMKapurAGonzalez-BosquetJHoribataSHeintzJAlbrechtRet al. MUC16 (CA125): Tumor Biomarker to Cancer Therapy, a Work in Progress. Mol Cancer (2014) 13:129. doi: 10.1186/1476-4598-13-129

  • CrossRef
  • Google Scholar

• 22

  LawrenceMSStojanovPPolakPKryukovGVCibulskisKSivachenkoAet al. Mutational Heterogeneity in Cancer and the Search for New Cancer-Associated Genes. Nature (2013) 499(7457):214–8. doi: 10.1038/nature12213

  • CrossRef
  • Google Scholar

• 23

  RumpAMorikawaYTanakaMMinamiSUmesakiNTakeuchiMet al. Binding of Ovarian Cancer Antigen CA125/MUC16 to Mesothelin Mediates Cell Adhesion. J Biol Chem (2004) 279(10):9190–8. doi: 10.1074/jbc.M312372200

  • CrossRef
  • Google Scholar

• 24

  HattrupCLGendlerSJ. Structure and Function of the Cell Surface (Tethered) Mucins. Annu Rev Physiol (2008) 70:431–57. doi: 10.1146/annurev.physiol.70.113006.100659

  • CrossRef
  • Google Scholar

• 25

  ComamalaMPinardMTheriaultCMatteIAlbertABoivinMet al. Downregulation of Cell Surface CA125/MUC16 Induces Epithelial-to-Mesenchymal Transition and Restores EGFR Signalling in NIH : OVCAR3 Ovarian Carcinoma Cells. Br J Cancer (2011) 104(6):989–99. doi: 10.1038/bjc.2011.34

  • CrossRef
  • Google Scholar

• 26

  O’BrienTJBeardJBUnderwoodLJDennisRASantinADYorkL. The CA 125 Gene: An Extracellular Superstructure Dominated by Repeat Sequences. Tumour Biol: J Int Soc Oncodevelopmental Biol Med (2001) 22(6):348–66. doi: 10.1159/000050638

  • CrossRef
  • Google Scholar

• 27

  YinBWLloydKO. Molecular Cloning of the CA125 Ovarian Cancer Antigen: Identification as a New Mucin, MUC16. J Biol Chem (2001) 276(29):27371–5. doi: 10.1074/jbc.M103554200

  • CrossRef
  • Google Scholar

• 28

  LiXPascheBZhangWChenK. Association of MUC16 Mutation With Tumor Mutation Load and Outcomes in Patients With Gastric Cancer. JAMA Oncol (2018) 4(12):1691–8. doi: 10.1001/jamaoncol.2018.2805

  • CrossRef
  • Google Scholar

• 29

  HuangJDengQWangQLiKYDaiJHLiNet al. Exome Sequencing of Hepatitis B Virus-Associated Hepatocellular Carcinoma. Nat Genet (2012) 44(10):1117–21. doi: 10.1038/ng.2391

  • CrossRef
  • Google Scholar

• 30

  ZhangJHuangJYChenYNYuanFZhangHYanFHet al. Whole Genome and Transcriptome Sequencing of Matched Primary and Peritoneal Metastatic Gastric Carcinoma. Sci Rep (2015) 5:13750. doi: 10.1038/srep13750

  • CrossRef
  • Google Scholar

• 31

  HuXZhaoYWeiLZhuBSongDWangJet al. CCDC178 Promotes Hepatocellular Carcinoma Metastasis Through Modulation of Anoikis. Oncogene (2017) 36(28):4047–59. doi: 10.1038/onc.2017.10

  • CrossRef
  • Google Scholar