Skip to main content

Genome-wide association study of serum tumor markers in Southern Chinese Han population



Serum indicators AFP, CA50, CA125, CA153, CA19-9, CEA, f-PSA, SCC-Ag have been confirmed as tumor markers (TMs). We conducted a genome-wide association study on 8 tumor markers of our 427 Han population in southern China, in order to identify genetic loci that are significantly associated with the level of 8 tumor markers.


We use Gene Titan multi-channel instrument and Axiom Analysis Suite 6.0 software for genotyping. We used IMPUTE2 software for imputation, and 1000 Genomes Project (Phase 3) was used as haplotype reference. After necessary quality control and statistical analysis, genetic loci genome-wide associated with TMs (p < 5E-8) will be identified. Finally, we selected Top SNPs (p < 5E-7) from the GWAS results for replication test. We used SPSS software to draw the distribution box plots of serum TMs under different genotypes of significant loci.


The results showed that there were only MUC1 (mucin 1)-rs4072037 significantly genome-wide associated with CA153 (p = 1.28E-18). However, we found that a total of 30 genetic loci have a suggestively significant genome-wide association with the level of 8 serum tumor markers (p < 5E-6). Then 3 Top SNPs (p < 5E-7) were selected for replication verification. The results showed that MUC1-rs4072037 was still significantly associated with CA153 in another population (p = 3.73E-08). Comparing with the TT genotype of rs4072037, the CA153 level was higher under CC or CT genotype of rs4072037.


MUC1-rs4072037 is significantly genome-wide associated with CA153 level. There are 30 genetic loci suggestively genome-wide associated with level of tumor markers among the Han population from Southern China.

Peer Review reports


Tumor marker (TM) is a series of substances that reflect the existence or development of tumors. It not only plays a role in the disease detection and target treatment of tumor patients, but also has important practical value in the anti-cancer screening of healthy people. Numerous studies have confirmed that the levels of certain related indicators in serum are closely associated with the occurrence and development of one or more tumors, and can be regarded as tumor markers. Such as alpha fetoprotein (AFP): liver cancer [1, 2]; CA19-9: pancreatic cancer; gastric cancer; colorectal cancer, etc. [3,4,5,6]; CA125: ovarian cancer [7]; CA153: breast cancer [8]; CA50: pancreatic cancer [9]; carcino embryonic antigen (CEA): colorectal cancer [10]; free Prostate Specific Antigen (f-PSA): prostate cancer [11]; squamous cell carcinoma antigen (SCC-Ag): squamous cell carcinomas [12], etc. In addition, previous studies suggested that the level of tumor markers such as CA19-9, CEA and AFP may differ between healthy individuals and cancer patients due to individual genetic differences and environmental factors [13,14,15,16,17]. Therefore, the detection of these tumor marker levels is helpful to monitor the occurrence and development of specific tumors in specific clinical populations.

Genome-wide association study (GWAS) has become the main common method for studying complex diseases and their susceptibility genes because of its ability to cover single nucleotide polymorphisms in the whole genome [18]. It can more effectively find genetic loci associated with the occurrence and development of diseases. GWAS can perform population-level statistical analysis of genotype and phenotype, so as to helping us identify genetic loci related to phenotypic changes. Up to now, numerous SNPs associated with levels of tumor makers among different populations have been reported [13, 17, 19,20,21]. SNPs significantly associated with levels of tumor markers in specific populations can be identified through GWAS. However, there are few studies few studies on tumor markers GWAS. He, M. et al. identified several loci associated with CA19-9, CEA and AFP levels through GWAS studies among 3451 healthy participants [13]. The above reminds us that it is of great significance to identify genetic loci associated with tumor marker levels through GWAS, which provides new ideas and data supplement for further understanding the genetic mechanism of serum TMs in specific populations.

Therefore, we conducted a genome-wide association study on 8 serum tumor markers (AFP, CA125, CA19-9, CA153, CA50, CEA, f-PSA, SCC-Ag) in the Han population from southern China. With a view to discovering genetic loci that are significantly associated with serum tumor markers in the Han population from southern China. Our study will provide valuable reference for clinical monitoring and diagnosis of the occurrence and development of tumors among Han population from southern China.

Materials and method

Study subjects and DNA extraction

The subjects consisted of 427 healthy individuals (221 males and 206 females) from health examination center. The inclusion criteria for all participants included: age mainly 40-50 years old, no family history of tumor, no disease (at least within 2 weeks), and no medication (at least within 2 weeks). Then the whole bloods of the participants were collected, and we extracted the whole genome DNA. The specific operation procedure was carried out according to the kit instructions (GoldMag, Xi’an). Subsequently, a GWAS study was conducted on 8 serum tumor markers in the participating population, which are AFP, CA50, CA125, CA153, CA19-9, CEA, f-PSA, SCC-Ag.

Genotyping and quality control

In our study, Thermo Fisher Genotyping Chip was used (Applied Biosystems™ Axiom™ Precision Medicine Diversity Array, PMDA). We use Gene Titan multi-channel instrument and Axiom Analysis Suite 6.0 software for genotyping. We performed a full gene scan through Axiom, and the results showed that within our sample range (427 participants), there were a total of 874,190 loci. After excluding Indel, copy number variation, sex chromosomes and duplicate sites, we then perform the necessary quality control on the remaining sites (sample call rate> 0.95, maker call rate> 0.90, HWE > 5 × 10-6). In the end, there are 796,288 loci left before we perform the imputation.

Imputation and quality control

We used IMPUTE2 software and used the haplotype of the 1000 Genomes Project (Phase 3) as a reference for imputation. After imputation, there are many poor quality loci, which need to be filtered out. Loci need to be filtered out include loci with MAF=0, loci with MAF<0.01 and info>0.3, and loci with excessive deletion rates (98% or more). Finally, the loci that meet the following conditions will be retained. These conditions include: sample call rate > 95%, marker call rate > 90%, HWE > 5 × 10-6, Allele = 2. In the end, a total of 6423,076 SNPs were used for subsequent analysis.

Statistical analysis

We use Gold Helix SNP & Variation Suite 8.7 version for association analysis. We mainly use the Mixed Linear Model-additive genetic model, and add the IBD matrix to the model to detect the SNPs associated with each indicator. All data in our study were adjusted by age and gender. At the same time, we also constructed Manhattan plots and Quantile-quantile plots related to serum TMs. The p-value is less than 5E-8, which means the genetic polymorphism is genome-wide significantly associated with the levels of tumor markers. Genetic polymorphism with p-value less than 5E-6 suggests that it may have a suggestively significant genome-wide association with tumor markers level [22].

Repeat detection

SNPs with p < 5E-7 were selected as Top SNPs from all SNPs identified by GWAS, and then they were repeatedly verified among 398 participants from the health examination center of Hainan General Hospital (different from the participants in the GWAS analysis). We removed the extreme measurement values (mean > ± 3 sd) of each serum TMs, and then normalized them with rankbaseINTs. Genotyping was performed by Agena MassARRAY. Then, the association analysis between Top SNPs and serum TMs was analyzed by Gold Helix SNP & Variation Suite 8.7. Through the mixed linear model-additive genetic model, while adding the IBD matrix to the model, the correlation between each tumor marker and Top SNPs will be analyzed. The final results were all adjusted by age and gender. And we also use False Discovery Rate analysis to perform correction on positive results. Finally, SPSS software was used to draw the distribution box plots of serum TMs under different genotypes of significant loci. In this part, p < 0.05 indicates that the candidate SNPs are significant associated with the 8 serum TMs.


GWAS identified 31 single nucleotide polymorphisms

A total of 427 Southern Chinese Han people participated in the GWAS assessment of this study. The number of participants, and ‘mean ± standard’ of each serum TMs were summarized in Table 1. There were differences in sample sizes for each indicator, which was due to data missing during sample collection.

Table 1 Sample information divided by different tumor markers

The results showed that there were only MUC1 (mucin 1)-rs4072037 significantly genome-wide associated with CA153 (p = 1.28E-18). However, we found that a total of 30 genetic loci have a suggestively significant genome-wide association with the level of 8 serum tumor markers. (p < 5E-6). Specifically, the genetic loci with suggestively genome-wide significances with 8 serum tumor markers are as follows: AFP: rs1524622 (LOC100287704/LOC100287834), rs143687479 (SGCZ), rs10425983 (ZNF444/GALP); CA50: rs112401204 (RBFOX2/APOL3), rs9295778 (LINC00533/LINC01623); CA125: rs75332508 (GCFC2/LRRTM4), rs12476806 (LINC01790/SLC39A10), rs7563672 (D2HGDH/GAL3ST2), rs6554915 (LOC728613/MIR4277), rs2028385 (CDK17/CFAP54); CA153: rs112804057 (PCDH7), rs12113394 (LOC100128317), rs7838015 (PTDSS1/LOC102724804), rs1053878 (ABO), rs8015097 (KCNH5), rs62406978 (RASGEF1C); CA19-9: rs67300776 (DPP10), rs1482587 (ADAMTS9-AS2/MAGI1), rs147309725 (ADRA1B/TTC1), rs4716999 (LOC389602/LOC285889); CEA: rs11060471 (TMEM132D), rs17568726 (FARP1); f-PSA: rs72641143 (PDE4B), rs466235 (PAM), rs78255505 (PCAT5/ANKRD30A), rs113984953 (PABPN1L/CBFA2T3), rs4816598 (ERG); SCC-Ag: rs111241781 (MIR6072/LINC00701), rs117089993 (HNF4G/LINC01111), rs72773580 (SHISA9/ERCC4).levelThe specific information was summarized in Table 2. We constructed Manhattan (Fig. 1) and Quantile–quantile (Supplementary Figure 1) plots. The horizontal line in Fig. 1 represents a cut-off value with suggestively genome-wide significance (p <5E-6). Supplementary Figure 2 showed a map of the regions associated with each serum tumor marker on different chromosomes.

Table 2 Genetic loci significantly associated with tumor markers (8 indicators) identified by the GWAS after imputation analysis in the Han population from Southern China
Fig. 1
figure 1

Manhattan graph of the results of the genome-wide association study. A AFP; B CA50; C CA125; D CA153; E CA19-9; F CEA; G f-PSA; H SCC-Ag. The chromosomes are displayed on the x-axis, while the y-axis represents the −log10 of the p-value. The red line represents the cut-off value of the suggestively genome-wide significance (5.0×10 -6)

Repeat verification

We selected TOP SNPs (p < 5E-7) from the 31 SNPs identified by GWAS to narrow the scope of verification. Finally, 3 SNPs were selected for the verification test. The results showed (Table 3) that MUC1-rs4072037 was still significantly associated with CA153 in different participants (p = 3.73E-08, FDR = 2.31E-06). The average level of CA153 under different genotypes of rs4072037 is summarized in Table 4. Compared with the wild genotype TT of MUC1-rs4072037, the level of CA153 under the CC/CT genotype was significantly increased (p < 0.001). Fig. 2 showed a box plot of the CA153 level changes under different genotypes of rs4072037.

Table 3 According to the results of GWAS, SNPs with p-value < 5E-07 were selected for replication testing in another population
Table 4 The average concentration of CA153 under different genotypes of rs4072037
Fig. 2
figure 2

The box plot of CA153 level under different genotypes of rs4072037


GWAS results of our study showed that only MUC1-rs4072037 significantly genome-wide associated with CA153 (p = 1.28E-18). However, a total of 30 genetic loci were identified in the Han population from Southern China, which were suggestively genome-wide associated with level of TMs. In addition, the results of repeated verification showed that only MUC1 rs4072037 still had a significant association with CA153 level.

The level of CA15-3 is increased in the serum of patients with malignant breast cancer and ovarian cancer [23, 24]. In this study, we found that the association between rs4072037 (MUC1) and CA15-3 levels reached the whole genome level. In addition, there were seven genetic loci that have a suggestive genome-wide significant association with CA153 level. Among the genetic loci identified in this study, rs4072037 (MUC1) [25, 26] and rs1053878 (ABO) [27] have been reported to be associated with the occurrence and development of cancers. However, their specific mechanism in cancer remains unclear. Gu, X. et al. found the important evidence that MUC1 rs4072037 may be used as a tumor marker in a large sample of ‘case Vs. control’ study [28]. Combining the results of this study, we speculate that MUC1-rs4072037 and ABO-rs1053878 may affect the expression level of CA153, thereby affecting the occurrence and development of cancer.

We only conducted GWAS of CA153 among 165 participants due to data missing in sample collection. Nevertheless, MUC1 rs4072037 was still detected to be significantly associated with CA153 levels during the repeated verification, which was passed FDR correction. It shows that this study result is worth believing. And compared with the TT genotype, the CA153 level under the CC or CT genotype was significantly increased. Recent studies have shown that the median expression of CA153 is significantly higher in a variety of cancers (lung cancer, breast cancer, ovarian cancer, cervical cancer) [29]. Combined with the results of our study, the appearance of the minor allele ‘C’ will significantly increase the level of CA153. Based on the above, MUC1 rs4072037 can be used as a new genetic signal for tumor prevention in southern Chinese Han population. To our knowledge, our study is the first to report a significant association between MUC1 rs4072037 and CA153 level. Our study will provide a valuable reference for tumor prevention.

The median survival rate of hepatocellular carcinoma patients with high serum AFP level is low [30]. A Japanese GWAS study found that the AT/TT genotype of rs17047200 on the TLL1 intron was associated with the e level of AFP [31]. Perhaps due to the different study methods, sample size or the genetic backgrounds of the participating populations, our results are different from the previous studies. In our study, the GWAS results showed that the genetic loci that have genome-wide significance with the AFP level in the Han population of southern China are rs1524622 (LOC100287704/LOC100287834), rs143687479 (SGCZ), rs10425983 (ZNF444/GALP). As far as we know, these three SNPs were discovered for the first time that may be associated with changes in AFP level.

High CEA level are associated with tumorigenesis and are often used to assist in the diagnosis of cancers caused by the gastrointestinal tract [32]. Carcinoembryonic antigen (CEA) is also highly expressed in malignant breast cancer [23]. Up to now, there were studies have identified 5 genetic variants of the ABO gene associated with CEA level through GWAS [13]. However, the association between the two SNPs (TMEM132D-rs11060471; FARP1-rs17568726) identified by GWAS and CEA level changes in our study was reported for the first time. They may be new genetic signals to monitor changes of CEA level in cancer patients, but this requires further experiments to verify.

Serum f-PSA is an important indicator for early detection of prostate cancer. Jeannette et al. found that KLK3 SNPs may be associated with f-PSA level in the GWAS assessment conducted among African Americans and Europeans [11]. Our results are different from those of previous studies. We found that the SNPs of PDE4B, PAM, PCAT5/ANKRD30A, PABPN1L/CBFA2T3 were significantly associated with changes in f-PSA level among the Chinese Southern Han population. We speculate that the reason for this difference may be influenced by different genetic background or environment. Nevertheless, our results still provide data supplements for GWAS studies on f-PSA in different populations.

Detection of serum SCC-Ag level can be used in the early diagnosis of cervical cancer recurrence [33]. Up to now, there is no study on the SCC-Ag GWAS has been reported. We conducted the SCC-Ag GWAS study for the first time in the Han population of southern China. And three genetic loci (rs111241781, rs117089993, rs72773580) have been identified that may be associated with the SCC-Ag level. They may be used as new genetic signals for early diagnosis of cervical cancer recurrence.

Carbohydrate antigen (CA) plays a key role in tumor progression, CA15-3, CA125, and CA19-9 are widely used in various cancer screening [34]. Studies have reported that CA19-9 is the most sensitive tumor marker for pancreatic cancer [3]. And the expression of CA19-9 is related to the recurrence [35]. CA125 is a reliable marker for clinical diagnosis of ovarian cancer [7], and its increasing level is a signal of recurrence of female germline tumors [36]. The expression of CA50 in the serum of pancreatic cancer patients was significantly higher than that of healthy controls [9]. Although some SNPs that are significantly associated to the level of CA19-9/CA125 have been reported in previous GWAS studies [13, 37], the genetic loci identified in our study have never been reported in Southern Chinese Han population. Our results provide reliable evidence that these SNPs can be used as new genetic signals for one or some tumor markers.

The characteristics and changes of tumor markers has important guiding significance for health assessment management, cancer prevention and tumor detection. Our study provides important evidence for in-depth understanding of the mechanism of serum tumor markers. The small sample size of this study is a limitation that cannot be ignored. In subsequent studies, we will expand the sample size or replace populations with different genetic backgrounds as study subjects to conduct validation tests. In any case, our study have provided a valuable reference for clinical monitoring and diagnosis of the occurrence and development of tumors among Han population from southern China.


In summary, MUC1-rs4072037 is significantly genome-wide associated with CA153 level (p = 1.28E-18) and 30 genetic loci were suggestively genome-wide associated with level of tumor markers (AFP, CA50, CA125, CA153, CA19-9, CEA, f-PSA, SCC-Ag) among the Han population from Southern China. Our study provides a valuable reference for clinical tumor prevention, and also provides a new theoretical basis for the phenotype research of tumor markers in healthy populations and individual health management.

Availability of data and materials

The datasets used and analyzed in the current study are available from the corresponding author on reasonable request.



Tumor markers


Alpha fetoprotein


Carcino embryonic antigen


Free Prostate Specific Antigen


Squamous cell carcinoma antigen


Genome-wide association study


  1. Cai J, Chen L. Genome-wide mapping of 5-hydroxymethylcytosines in circulating cell-free DNA as a non-invasive approach for early detection of hepatocellular carcinoma. Gut. 2019;68(12):2195–205.

    Article  CAS  Google Scholar 

  2. Chen W, Peng J, Ye J, Dai W, Li G, He Y. Aberrant AFP expression characterizes a subset of hepatocellular carcinoma with distinct gene expression patterns and inferior prognosis. J Cancer. 2020;11(2):403–13.

    Article  CAS  Google Scholar 

  3. Scarà S, Bottoni P, Scatena R. CA 19-9: Biochemical and Clinical Aspects. Adv Exp Med Biol. 2015;867:247–60.

    Article  Google Scholar 

  4. Yajima H, Omura N, Matai K, Mitsumori N, Yoshida K, Yanaga K. Clinicopathological features of CA19-9-producing gastric cancer. Hepatogastroenterology. 2014;61(129):221–5.

    CAS  PubMed  Google Scholar 

  5. Vanderbeck K, Voutsadakis IA. Elevation of Ca19-9 tumor antigen in colorectal cancer: an in silico investigation of pathogenesis. Int J Color Dis. 2016;31(5):1105–7.

    Article  Google Scholar 

  6. Chen Z, Liu Z. Combination of CA19-9 and the neutrophil-to-lymphocyte ratio for the differential diagnosis of gallbladder carcinoma. Cancer Manag Res. 2020;12:4475–82.

    Article  CAS  Google Scholar 

  7. Felder M, Kapur A, Gonzalez-Bosquet J, Horibata S, Heintz J, Albrecht R, et al. MUC16 (CA125): tumor biomarker to cancer therapy, a work in progress. Mol Cancer. 2014;13:129.

    Article  Google Scholar 

  8. Bon GG, Kenemans P, Verstraeten AA, Go S, Philipi PA, van Kamp GJ, et al. Maternal serum Ca125 and Ca15-3 antigen levels in normal and pathological pregnancy. Fetal Diagn Ther. 2001;16(3):166–72.

    Article  CAS  Google Scholar 

  9. Lei XF, Jia SZ, Ye J, Qiao YL, Zhao GM, Li XH, et al. Application values of detection of serum CA199, CA242 and CA50 in the diagnosis of pancreatic cancer. J Biol Regul Homeost Agents. 2017;31(2):383–8.

    CAS  PubMed  Google Scholar 

  10. Bajenova O, Gorbunova A, Evsyukov I, Rayko M, Gapon S, Bozhokina E, et al. The genome-wide analysis of carcinoembryonic antigen signaling by colorectal cancer cells using RNA sequencing. PLoS One. 2016;11(9):e0161256.

    Article  Google Scholar 

  11. Bensen JT, Xu Z, Smith GJ, Mohler JL, Fontham ET, Taylor JA. Genetic polymorphism and prostate cancer aggressiveness: a case-only study of 1,536 GWAS and candidate SNPs in African-Americans and European-Americans. Prostate. 2013;73(1):11–22.

    Article  CAS  Google Scholar 

  12. Hamada K, Hanakawa Y, Hashimoto K, Iwamoto M, Kihana T, Hirose S, et al. Gene expression of human squamous cell carcinoma antigens 1 and 2 in human cell lines. Oncol Rep. 2001;8(2):347–54.

    CAS  PubMed  Google Scholar 

  13. He M, Wu C, Xu J, Guo H, Yang H, Zhang X, et al. A genome wide association study of genetic loci that influence tumour biomarkers cancer antigen 19-9, carcinoembryonic antigen and α fetoprotein and their associations with cancer risk. Gut. 2014;63(1):143–51.

    Article  CAS  Google Scholar 

  14. Zhang G, Xu Y, Lu X, Huang H, Zhou Y, Lu B, et al. Diagnosis value of serum B7-H3 expression in non-small cell lung cancer. Lung Cancer. 2009;66(2):245–9.

    Article  Google Scholar 

  15. Evdokimova VN, Liu Y, Potter DM, Butterfield LH. AFP-specific CD4+ helper T-cell responses in healthy donors and HCC patients. J Immunother. 2007;30(4):425–37.

    Article  CAS  Google Scholar 

  16. Hatate K, Yamashita K, Hirai K, Kumamoto H, Sato T, Ozawa H, et al. Liver metastasis of colorectal cancer by protein-tyrosine phosphatase type 4A, 3 (PRL-3) is mediated through lymph node metastasis and elevated serum tumor markers such as CEA and CA19-9. Oncol Rep. 2008;20(4):737–43.

    PubMed  Google Scholar 

  17. Motzer RJ, Hutson TE, Hudes GR, Figlin RA, Martini JF, English PA, et al. Investigation of novel circulating proteins, germ line single-nucleotide polymorphisms, and molecular tumor markers as potential efficacy biomarkers of first-line sunitinib therapy for advanced renal cell carcinoma. Cancer Chemother Pharmacol. 2014;74(4):739–50.

    Article  CAS  Google Scholar 

  18. Hirschhorn JN, Daly MJ. Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005;6(2):95–108.

    Article  CAS  Google Scholar 

  19. Berković MC, Jokić M, Marout J, Radosević S, Zjacić-Rotkvić V, Kapitanović S. IL-2 -330 T/G SNP and serum values-potential new tumor markers in neuroendocrine tumors of the gastrointestinal tract and pancreas (GEP-NETs). J Mol Med (Berl). 2010;88(4):423–9.

    Article  Google Scholar 

  20. Zhang J, Chi H, Xiao H, Tian X, Wang Y, Yun X, et al. Interleukin 6 (IL-6) and Tumor Necrosis Factor α (TNF-α) Single Nucleotide Polymorphisms (SNPs), Inflammation and Metabolism in Gestational Diabetes Mellitus in Inner Mongolia. Med Sci Monit. 2017;23:4149–57.

    Article  Google Scholar 

  21. Grigorova AA, Trenova AG, Stanilova SA. Association of polymorphism -308G/A in tumor necrosis factor-alpha gene (TNF-α) and TNF-α serum levels in patients with relapsing-remitting multiple sclerosis. Neurol Res. 2021;43(4):291–8.

    Article  CAS  Google Scholar 

  22. Son HY, Hwangbo Y, Yoo SK, Im SW, Yang SD, Kwak SJ, et al. Genome-wide association and expression quantitative trait loci studies identify multiple susceptibility loci for thyroid cancer. Nat Commun. 2017;8:15966.

    Article  CAS  Google Scholar 

  23. Dolscheid-Pommerich RC, Keyver-Paik M, Hecking T, Kuhn W, Hartmann G, Stoffel-Wagner B, et al. Clinical performance of LOCI™-based tumor marker assays for tumor markers CA 15-3, CA 125, CEA, CA 19-9 and AFP in gynecological cancers. Tumour Biol. 2017;39(10):1010428317730246.

    Article  Google Scholar 

  24. Williams KA, Terry KL, Tworoger SS, Vitonis AF, Titus LJ, Cramer DW. Polymorphisms of MUC16 (CA125) and MUC1 (CA15.3) in relation to ovarian cancer risk and survival. PLoS One. 2014;9(2):e88334.

    Article  Google Scholar 

  25. Duan F, Song C, Dai L, Cui S, Zhang X, Zhao X. The effect of MUC1 rs4072037 functional polymorphism on cancer susceptibility: evidence from published studies. PLoS One. 2014;9(4):e95651.

    Article  Google Scholar 

  26. Kim BS, Lee I. Association between the MUC1 rs4072037 Polymorphism and Risk of Gastric Cancer and. Clin Outcomes. 2020;20(2):127–38.

    Google Scholar 

  27. Cozzi GD, Levinson RT, Toole H, Snyder MR, Deng A, Crispens MA, et al. Blood type, ABO genetic variants, and ovarian cancer survival. PLoS One. 2017;12(4):e0175119.

    Article  Google Scholar 

  28. Gu X, Feng J, Liu L, Lu M, Ma X, Cao Y, et al. Association of MUC1 rs4072037 Functional Polymorphism and Cancer Risk: Evidence from 12551 Cases and 13436 Controls. J Cancer. 2018;9(18):3343–51.

    Article  Google Scholar 

  29. Wang Y, Lei L, Chi YG, Liu LB, Yang BP. A comprehensive understanding of ovarian carcinoma survival prognosis by novel biomarkers. Eur Rev Med Pharmacol Sci. 2019;23(19):8257–64.

    CAS  PubMed  Google Scholar 

  30. Tangkijvanich P, Anukulkarnkusol N, Suwangool P, Lertmaharit S, Hanvivatvong O, Kullavanijaya P, et al. Clinical characteristics and prognosis of hepatocellular carcinoma: analysis based on serum alpha-fetoprotein levels. J Clin Gastroenterol. 2000;31(4):302–8.

    Article  CAS  Google Scholar 

  31. Matsuura K, Sawai H, Ikeo K, Ogawa S, Iio E, Isogawa M, et al. Genome-Wide Association Study Identifies TLL1 Variant Associated With Development of Hepatocellular Carcinoma After Eradication of Hepatitis C Virus Infection. Gastroenterology. 2017;152(6):1383–94.

    Article  CAS  Google Scholar 

  32. Goldenberg DM, Neville AM, Carter AC, Go VL, Holyoke ED, Isselbacher KJ, et al. CEA (carcinoembryonic antigen): its role as a marker in the management of cancer. J Cancer Res Clin Oncol. 1981;101(3):239–42.

    Article  CAS  Google Scholar 

  33. Salvatici M, Achilarre MT, Sandri MT, Boveri S, Vanna Z, Landoni F. Squamous cell carcinoma antigen (SCC-Ag) during follow-up of cervical cancer patients: Role in the early diagnosis of recurrence. Gynecol Oncol. 2016;142(1):115–9.

    Article  CAS  Google Scholar 

  34. Sturgeon C. Practice guidelines for tumor marker use in the clinic. Clin Chem. 2002;48(8):1151–9.

    Article  CAS  Google Scholar 

  35. Sakamoto T, Saito H, Uchinaka EI, Morimoto M, Amisaki M, Tokuyasu N, et al. The Combination of Neutrophil-to-lymphocyte Ratio and Serum Carbohydrate Antigen 19-9 Level as a Prognostic Indicator in Patients with Recurrent Pancreatic Cancer. Anticancer Res. 2018;38(9):5497–503.

    Article  CAS  Google Scholar 

  36. Xiao C, Zhao J, Guo P, Wang D, Zhao D, Ren T, et al. Clinical analysis of primary primitive neuroectodermal tumors in the female genital tract. Int J Gynecol Cancer. 2014;24(3):404–9.

    Article  Google Scholar 

  37. Rupp C, Friedrich K, Folseraas T, Wannhoff A, Bode KA, Weiss KH, et al. Fut2 genotype is a risk factor for dominant stenosis and biliary candida infections in primary sclerosing cholangitis. Aliment Pharmacol Ther. 2014;39(8):873–82.

    Article  CAS  Google Scholar 

Download references


We thank all the participants from Hainan General Hospital for providing blood samples and all people involved in this study.


This work was supported by the National Natural Science Foundation of China [grant numbers No. 81660013, No.81860015] and the Key Research and Development Plan of Hainan Province [grant number No. ZDYF2018116].

Author information

Authors and Affiliations



YPD and PH conceived and designed the experiments; XKL and FHB performed the experiments; XWW, TBJ, CL and YTZ collected samples; YTZ, XWW, ML and XLZ analyzed the data; YTZ, YFX, CYH and QL contributed reagents/materials/analysis tools; XKL and FHB drafted the paper. SYC revised and reviewed the paper. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Ping He, Shuyuan Chu or Yipeng Ding.

Ethics declarations

Competing of interests

The authors declared that they have no conflicts of interest.

Ethics approval and consent to participate

This study was conducted under the standard approved by the Ethics Committee of Hainan General Hospital, and conformed to the ethical principles for medical research involving humans of the World Medical Association Declaration of Helsinki. All participants signed informed consent forms before participating in this study.

Consent for publication

All authors agreed to publish the manuscript.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplemental Figure 1.

Quantile–quantile plots about the results of the GWAS. (A) AFP; (B) CA50; (C) CA125; (D) CA153; (E) CA19-9; (F) CEA; (G) f-PSA; (H) SCC-Ag.

Additional file 2: Supplemental Figure 2.

A map of the SNPs associated with each serum tumor marker on different chromosomes. (A) AFP; (B) CA50; (C) CA125; (D) CA153; (E) CA19-9; (F) CEA; (G) f-PSA; (H) SCC-Ag.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, X., Bai, F., Wei, X. et al. Genome-wide association study of serum tumor markers in Southern Chinese Han population. BMC Cancer 22, 160 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Tumor markers
  • GWAS
  • Imputation
  • Southern Chinese Han population