Skip to main content

Polymorphism rs2327430 in TCF21 predicts the risk and prognosis of gastric cancer by affecting the binding between TFAP2A and TCF21

Abstract

Background

Single nuclear polymorphisms (SNPs) have been published to be correlated with multiple diseases. Transcription Factor 21 (TCF21) is a critical transcription factor involved in various types of cancers. However, the association of TCF21 genetic polymorphisms with gastric cancer (GC) susceptibility and prognosis remains unclear.

Methods

A case-control study comprising 890 patients diagnosed with GC and an equal number of cancer-free controls was conducted. After rigorous statistical analysis, molecular experiments were carried out to elucidate the functional significance of the SNPs in the context of GC.

Results

TCF21 rs2327430 (OR = 0.78, P = 0.026) provides protection against GC, while rs4896011 (OR = 1.39, P = 0.005) exhibit significant associations with GC risk. Furthermore, patients with the (TC + CC) genotype of rs2327430 demonstrate a relatively favorable prognosis (OR = 0.47, P = 0.012). Mechanistically, chromatin immunoprecipitation assay and luciferase reporter assay revealed that the C allele of rs2327430 disrupts the binding of Transcription Factor AP-2 Alpha (TFAP2A) to the promoter region of TCF21, resulting in increased expression of TCF21 and inhibition of malignant behaviors in GC cells.

Conclusion

Our findings highlight the significant role of TCF21 SNPs in both the risk and prognosis of GC and provide valuable insights into the underlying molecular mechanisms. Specifically, the disruptive effect of rs2327430 on TCF21 expression and its ability to modulate malignant cell behaviors suggest that rs2327430 may serve as a potential predictive marker for GC risk and prognosis.

Introduction

Despite the availability of numerous potential therapies, GC remains the fifth most common cancer in terms of incidence and the fourth leading cause of cancer-related death worldwide, posing a significant threat to human health [1]. According to recent data, GC was responsible for 397,000 new cases and approximately 290,000 deaths in China in 2016 [2], most of whom were diagnosed at advanced stages. Subsequent studies have consistently demonstrated that early treatment could have a positive impact on the prognosis of GC patients [3]. Therefore, early screening for GC among populations with high risk could contribute to improving the patient’s prognosis and overall survival [4, 5].

Single nucleotide polymorphism (SNP) is a common form of heritable DNA variation linked to tumor development and progression [6, 7]. Numerous studies have suggested that SNPs could assist in detecting malignant tumors at an early stage, predicting the risk in key populations, and assessing the prognosis of patients [8]. For instance, the rs3807213 variant genotypes of IFRD1 have been associated with increased susceptibility to GC [9]. Moreover, the rs1141023 variant located in TRIM59 has been linked to a rising predisposition to GC, potentially correlated with early-stage GC [10]. The effects of transcription factors on the malignant cell behaviors of GC have also been widely reported [11, 12]. However, comprehensive empirical support regarding the effects of SNPs located in transcription factors in GC is still lacking. Thus, our study aimed to explore the role of TCF21 SNPs in GC.

The TCF21 gene encodes a basic helix-loop-helix transcription factor at 6q23.2 [13, 14]. TCF21 is downregulated and plays protective roles in various malignant tumors, such as hepatocellular carcinoma, ovarian cancer, and renal tumors [13, 14]. Yang et al. demonstrated that TCF21 was downregulated and suppressed malignant behaviors through the AKT pathway in GC [15]. Furthermore, researchers have suggested that rs12190287 located in TCF21 is associated with a poor prognosis in breast cancer [16]. To investigate the function of SNPs located in TCF21, we selected four tagSNPs and examined their relationship with GC.

Materials and methods

Study subjects

All participants included in the study comprised 890 GC cases histopathologically confirmed by at least two pathologists and 890 cancer-free controls without any history of malignancy. All patients received medical treatment at the First Affiliated Hospital of Nanjing Medical University between February 2012 and March 2020. Exclusion criteria and detailed definitions are provided in a previous study [1].

SNP selection

Several studies suggest that selecting tagSNPs instead of all SNPs in a given interval does not result in a substantial loss of power, and it is reasonable to choose potentially functional SNPs and a selection of tagSNPs as objects of a study [17, 18]. All tagSNPs of the TCF21 gene full-length and upstream 2000 bp in chromosome 6, position 133,887,113–133,895,553 in the Chinese population were selected on the grounds of the 1000 Genomes Project and HaploView 4.2 software under the following conditions: r2 > 0.8 and a minor allele frequency (MAF) > 0.05. Previous studies have suggested that SNPs in nonsynonymous mutated regions, splice site regions, and untranslated regions (UTRs) have a higher potential to have functional effects than SNPs in intronic regions, as introns are typically removed during the formation of mature mRNA [19, 20]. Therefore, we excluded SNPs in the intron region and identified five SNPs [rs12190287 (3’-UTR), rs2327429 (promoter), rs2327430 (promoter), rs3734281 (3’-UTR), rs4896011 (3’-UTR)] that potentially influence the expression of TCF21. Based on respective genomic locations, the JASPAR database and miRNASNP website were used to analyze the function of sites situated in promoter regions or 3’UTRs, respectively (Fig. S1A-C). Due to the lack of predicted transcription factors associated with rs2327429, which is located in the promoter region of TCF21, the four remaining sites were considered for inclusion in the study.

Extraction and genotyping of DNA

Extraction of DNA was performed following our previous introductions [1]. The genotypes of the chosen tagSNPs were determined by the SNaPshot genotyping method, which has been validated in previous studies [21, 22]. Here are the primer sequences used (Table S1).

Tissue samples and cells

A total of 120 pairs of tumorous and adjacent tissues were collected from radical gastric resection specimens of patients with GC in the First Affiliated Hospital of Nanjing Medical University. The study employed two human GC cell lines, MKN45 and AGS, obtained from the cell bank of the Chinese Academy of Sciences (Shanghai, China). The preservation of tissues and cell culture conditions were performed concerning a previous study [23]. All cells used in this study were authenticated using STR profiling within the last three years, and all experiments were performed with mycoplasma-free cells.

RNA extraction and quantitative real-time PCR (qRT-PCR)

RNA was extracted from cells as well as 120 out of 890 pairs of tumorous tissues and adjacent normal tissues using TRIzol Reagent (Invitrogen, USA). Following reverse transcription, the expression of mRNA was measured by quantitative real-time PCR (ABI 7300). The primers used in this research are shown in the supplementary material. All experiments were performed in triplicate.

Plasmid construction, siRNA interference, and transfection

We constructed luciferase reporter plasmids and TFAP2A-overexpressing plasmids in which the sequences contained the 2 kb upstream of the transcription start sites (TSS) of TCF21 with different alleles of rs2327430. The siRNA sequences designed against TFAP2A and TCF21 whose efficiency has been validated in other research are listed in the supplementary material [24, 25]. The rs2327430 wild (T) or mutant (C) vectors were obtained from GenePharma (GenePharma, Shanghai, China). When the cells were grown to 40–60% confluencein six-well plates, Lipofectamine 3000 (Invitrogen) and p3000 were used to transfect small interfering RNA (siRNA) and plasmids into the aforementioned cells.

Western blotting

The protein extraction method and detailed steps of western blotting were adapted from a previous study [23]. The following antibodies were used: TCF21 and cleaved caspase 3 from Abcam (Britain) and p-AKT, AKT, caspase 3, BCL-xL, TFAP2A, and GAPDH from Proteintech (Wuhan, China).

Immunohistochemistry (IHC)

After fixing and then embedding all specimens, we de-waxed the paraffin-embedded sections in xylene and incubated them with corresponding antibodies.

Luciferase assay

In total, 0.4 µg luciferase reporter plasmids containing different alleles of rs2327430 were transfected into MKN45 and AGS cells, and equivalent amounts of TFAP2A-overexpressing plasmids were cotransfected for the luciferase assay. The empty vector served as a negative control. The luciferase assay was performed on the Dual-Luciferase Reporter Assay System (Promega, USA), and the outcomes were evaluated based on firefly luciferase activities normalized to Renilla luciferase.

Chromatin immunoprecipitation (ChIP)

For the ChIP assay, we employed a ChIP assay kit (Beyotime, China). Broadly, we sonicated the cross-linked chromatin DNA into 200–1000 bp fragments that were immunoprecipitated with an anti-TFAP2A antibody (Proteintech, China), while normal rabbit IgG was used as the reference group. Quantification of the immunoprecipitated DNA was assessed by qRT-PCR with SYBR Green Mix (Vazyme Biotech Co., Ltd, China), and the results are presented relative to input (% input). The TCF21 promoter primers used are as follows: forward 5′-AGATGGACAGAACATGCTGC-3′ and reverse 5′-AGGGAAACTCAATGCACAGA CC-3′.

Cell proliferation, wound healing, and transwell assays

As described in our previous study [23], cell counting kit 8 (CCK-8) and colony formation assays were performed to validate the influence of TFAP2A on the proliferative capacity of pre-treated GC cells altered by TCF21, while transwell assays and wound healing assays were used to investigate the effect on cell migration and invasion. All the above experiments were conducted in triplicate.

Flow cytometric analysis of apoptosis

We conducted the apoptosis assay using an Annexin V-APC/PI Apoptosis Detection Kit (Multisciences, China). The obtained data were analyzed with flow cytometry (BD Biosciences, USA). The total apoptosis rate was calculated based on the ratio of early and late apoptotic cells.

Animal experiment

To construct the tumor xenograft model, a total of 30 female nude mice were randomly assigned. Pre-treated cells (1 × 106 cells/100 µl of PBS) were injected into the flank of each nude mouse (BALB/c) in the different groups, respectively. We measured tumor volume every 2 days and weighed xenografts after sacrificing the mice.

False-positive report probability

According to previously reported studies [26], false-positive report probability (FPRP) was provided, and a false-positive rate between genetic polymorphisms and GC hazard was available (Table S2).

Statistical analysis

SPSS22.0 software (SPSS Inc., Chicago, IL, USA) was chosen for performing the statistical analysis in our research. Statistical methods were used as described in our previous studies [1].

Results

Demographic information of the study subjects

The baseline characteristics of the 890 cases and 890 cancer-free controls are summarized in Table 1. These characteristics were well-matched between the two groups (P > 0.05). Among the GC patients, there were 490 cases of non-cardia cancer (55.1%). Cases with lymph node metastasis accounted for a greater proportion (62.0%), while a lower proportion had lymphovascular invasion (39.6%) and perineural invasion (39.6%). The distributions of the depth of invasion are presented below: 24.7% (T1), 11.8% (T2), 25.5% (T3), and 38.0% (T4). By the 8th edition of the AJCC cancer staging manual, the TNM stages of 890 patients were categorized as follows: stage I (27.9%), stage II (23.9%), stage III (42.7%), and stage IV (5.5%).

Table 1 Demographic information

Associations between TCF21 polymorphisms and the risk of GC

Due to the absence of allele C (n = 0) in the genotyping results, rs3734281 was not investigated. Table 2 shows the frequencies of genotypes and alleles at each selected SNP locus. However, the examination of HWE showed disequilibrium in rs12190287 (Table S3); therefore, further analysis was not conducted for this SNP. The results of rs2327430 and rs4896011 were the focus of our study. As presented in Table 2, the presence of the C allele of rs2327430 could decrease the risk of GC occurrence (OR = 0.78, 95% CI = 0.63–0.97), which was also observed in the dominant model (OR = 0.74, 95% CI = 0.58–0.93). In addition, the TC genotype in the codominant model played an active role in reducing the risk of GC (OR = 0.72, 95% CI = 0.56–0.92). Conversely, the A allele in rs4896011 was associated with an increased risk of GC in the codominant mode (OR = 1.36, 95% CI = 1.07–1.75) as well as the dominant model (OR = 1.39, 95% CI = 1.09–1.77).

Table 2 The association between TCF21 gene polymorphisms (rs2327430, rs4896011) and the risk of gastric cancer

Stratification analysis of TCF21 SNPs and GC risk

To clarify the potential effects of risk factors and comorbidities related to GC, a subgroup analysis was conducted using dominant models of rs2327430 and rs4896011. Table 3 presents the results for rs2327430, indicating that the protective effects of mutant genotypes TC + CC were particularly significant among females (OR = 0.52, 95% CI = 0.33–0.82) and patients residing in urban areas (OR = 0.67, 95% CI = 0.47–0.95), without hypertension (OR = 0.67, 95% CI = 0.50–0.90), without diabetes (OR = 0.76, 95% CI = 0.59–0.98), without smoking history (OR = 0.68, 95% CI = 0.56–0.94), and without drinking history (OR = 0.73, 95% CI = 0.56–0.95). For rs4896011, the AT + AA variant genotypes were found to significantly increase the risk of GC in males (OR = 1.56, 95% CI = 1.17–2.08), individuals without hypertension (OR = 1.67, 95% CI = 1.22–2.29), individuals without diabetes (OR = 1.40, 95% CI = 1.07–1.83), nonsmokers (OR = 1.48, 95% CI = 1.12–1.97), nondrinkers (OR = 1.40, 95% CI = 1.07–1.84), and individuals living in rural areas (OR = 1.40, 95% CI = 1.00-1.97). No obvious correlations were observed in the opposite subgroups.

Table 3 Stratified analysis for rs2327430 and rs4896011 genotype in cases and controls

Further analyses were conducted to assess the correlations between the aforementioned tagSNPs and various clinicopathological characteristics. However, obvious proof of relevance between them was rarely obtained (Table S4).

Associations of TCF21 polymorphisms with clinicopathological characteristics and the prognosis of GC patients

Analysis of the correlation between the clinicopathological features and the prognosis of GC patients indicated that larger tumor size, positive lymphovascular invasion, more advanced TNM stage (III and IV), and the diffuse type based on Lauren classification were independent predictors for poor outcome of GC patients (Table 4) (Fig. 1A-F), which aligns with their recognized roles. Based on the preliminary findings of this study, we explored the potential effects of the two aforementioned SNPs on the overall survival of GC patients. Unfortunately, analyses between wild-type homozygotes and mutant homozygotes (with SNPs affecting both alleles) could not be performed, and it was not possible to make predictions based on the MST data due to the limitation of the small sample size and single survival outcome of mutant homozygote samples. According to the results obtained thus far, the TC genotype (HR = 0.50, 95% CI = 0.27–0.91) and TC + CC genotype (HR = 0.47, 95% CI = 0.26–0.85) in rs2327430 could have a favorable impact on prognosis (Fig. 1G-H), while the polymorphisms in rs4896011 did not correlate with overall survival in GC patients when tumor size, lymphovascular invasion, TNM stage, and Lauren classification were considered (Table 5).

Table 4 The association of patients’ clinical features and overall survival time
Fig. 1
figure 1

The associations of various factors with the prognosis of GC patients. Notes: (A–F) Kaplan-Meier survival curves for the overall survival by characteristics and clinical features of gastric cancer patients. (tumor size, tumor differentiation, TMN stage, LVI, PNI, Lauren classification.) (G and H) Kaplan-Meier survival curves of rs2327430 polymorphism for the overall survival in patients with gastric cancer. (TT vs. TC, TT vs. CC, TT vs. TC + CC.)

Table 5 Associations between rs2327430 and rs4896011 genotypes and GC overall survival in 627 GC patients

The C allele of rs2327430 modulates malignant behaviors in GC cells

Given the protective roles of rs2327430 in both the risk and prognosis of GC, we further explored the underlying mechanism. Bioinformatic analyses showed a high abundance of promoter and enhancer histone marks in the 2000 bp upstream genomic region of TCF21, especially around rs2327430 (Fig. S1A-B). The impact of rs2327430 on the open chromatin state was also investigated by single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) (Fig. S1D). The 3DSNP database and UCSC database were utilized to visualize the structure of rs2327430 through circular and linear plots, respectively (Fig. S1E-F). In addition, we predicted the role of rs2327430 in modulating RNA secondary structure (Fig. S1G). Preliminary exploratory analysis indicated that rs2327430 may regulate promoter activity by altering the binding of a specific transcription factor. First, we quantified the expression of TCF21 in 120 GC tissues by qRT-PCR. The results indicated that cases carrying the rs2327430 C allele had a higher expression level of TCF21 than those with the TT genotype (Fig. 2A-B). Based on the above observations, we examined the expression level of TCF21 in cells transfected with plasmids carrying different alleles of rs2327430, and the rs2327430 C allele caused an increase in TCF21 expression at both the RNA level and protein level (Fig. 2C-D). The CCK-8 assay and colony formation assay validated that the C allele could inhibit the proliferation of GC cells (Fig. 2E-H). The migration and invasion capacity of GC cells may be weakened by the C allele of rs2327430 (Fig. 2I-M). The C allele may also be associated with a higher apoptosis rate in GC cells than the T allele (Fig. 2N-O). A previous study validated that TCF21 plays a role in GC through the AKT-Bcl-xL signaling pathway [15]. As a result, we explored whether there was a correlation between rs2327430 and the AKT-Bcl-xL signaling pathway, and we obtained a positive result (Fig. 2P).

Fig. 2
figure 2

The C allele in rs232730 inhibited the malignant behavior of GC cell lines relative to the T allele by affecting the expression level of TCF21 Notes: (A and B) Correlation between different rs2327430 genotypes and TCF21 expression levels in tumor tissues. (C and D) The expression level of TCF21 was analyzed in GC cell lines transfected with different alleles of rs2327430 by qRT-PCR and western blots. (E–H) CCK-8 and colony formation assay showed the C allele in rs2327430 could inhibit the proliferation ability of GC cell lines. (I–M) The wound healing assay and transwell assay were performed to show the effect of rs2327430 on the migration and invasion of GC cell lines. (N and O) The apoptotic rate of GC cell lines transfected with different alleles in rs2327430 was examined by flow cytometry. (P) The effect of different alleles in rs2327430 on the AKT/Bcl-xL signaling pathway was examined by western blots. (Q) Images of the subcutaneous xenografts were presented. (R) tumor growth curves were drawn and the tumor weight was calculated. (S) IHC was used to examine the expression of TCF21 and Ki67 in xenografts

In vivo, subcutaneous xenograft tumor models were constructed, and the expression of TCF21 as well as Ki67 was examined by IHC, whose results validated the aforementioned conclusion as well (Fig. 2Q-S).

The C allele in rs2327430 blocks the binding of TFAP2A to the TCF21 promoter

To understand the underlying mechanisms, we utilized the JASPAR database and found that there existed a TFAP2A binding motif in rs2327430 (Fig. 3A-B). Analysis of molecular correlates also showed that the expression level of TCF21 was negatively correlated with the expression level of TFAP2A (Fig. 3C). To further validate this conjecture, various types of luciferase reporter vectors with promoters containing either the rs2327430 allele (T or C) were constructed and transfected into MKN45 (rs2327430-TT) and AGS (rs2327430-TC) cells for the luciferase assay. Additionally, a plasmid overexpressing TFAP2A was cotransfected into these cells. The cells transfected with a vector carrying the C allele exhibited significantly increased luciferase activity compared to those with the T allele, and the overexpression of TFAP2A negatively affected luciferase activity (Fig. 3D). The study revealed a negative correlation between the expression levels of TFAP2A and the transcriptional activity of TCF21. Furthermore, it was observed that the T allele inhibited the transcription of TCF21. The direct binding of TFAP2A to rs2327430 in TCF21 was confirmed by ChIP allele-specific qPCR assay in the two cell lines (Fig. 3E). To perform the following rescue assays, the efficiency of siRNAs in knockdown of TCF21 and TFAP2A was determined by qRT‒PCR and western blots, respectively (Fig. 3F-G). TFAP2A promotes malignant cell behaviors in GC, while TCF21 plays a protective role in GC [15, 25]. To further substantiate these findings, rescue experiments were conducted in vitro. Quantitative qRT-PCR and western blotting were employed to validate the impact of TFAP2A on the expression of TCF21 (Fig. 3H). In accordance with the above results, the knockdown of TFAP2A partially rescued several functional phenotypes, such as proliferation, migration, invasion, and apoptosis when TCF21 was also knocked down (Fig. 3I-O). A similar phenomenon was also observed in vivo experiments (Fig. 3P-Q). In summary, the polymorphisms in rs2327430 influenced malignant behaviors and the overall survival of GC patients by altering the binding with TFAP2A (Fig. 4).

Fig. 3
figure 3

The C allele of rs2327430 played a role in GC by blocking the binding of TFAP2A to the TCF21 promoter. Notes: (A and B) The potential of the binding between tagSNPs and transcription factor was analyzed and displayed according to the JASPAR database. (C) Correlation analysis between TFAP2A and TCF21 mRNA was performed in 120 GC patients and analyzed with Spearman’s correlation test. (D) Schematic illustration of the reporter gene containing rs2327430 C or T allele constructs and reporter plasmids with the different alleles of rs2327430 were transfected into AGS and MKN45 cells. Then, the results were expressed as relative luciferase activity (Firefly luciferase/Renilla luciferase). Both cells were also co-transfected with TFAP2A plasmid, respectively. (E) CHIP assay evaluated by qRT-PCR was performed with control IgG or antibody against TFAP2A in the AGS and MKN45 cells, results of which are normalized to the input group and shown as means ± SEM in 3 independent experiments. (F and G) The knock-down efficiency of siRNAs on TCF21 and TFAP2A was determined by qRT-PCR and western blots, respectively. (H) The rescue experiment was carried out to verify whether the expression of TCF21 is affected by TFAP2A through qRT-PCR and western blot. (I–N) It is shown that TCF21 knockdown could partially rescue functional phenotypes caused by TFAP2A knockdown including proliferation, invasion, metastasis, and apoptosis. (O) The effect of TFAP2A knockdown on AKT/Bcl-xL signaling pathway could be rescued partially by TCF21 knockdown. (P) Images of the subcutaneous xenografts were presented. tumor growth curves were drawn and the tumor weight was calculated. (Q) IHC was used to examine the expression of TCF21 and Ki67 in xenografts

Fig. 4
figure 4

The role of rs2327430 in modulating gastric cancer. Notes: The allele C played a protective role in gastric cancer development by inhibiting the binding with TFAP2A

Discussion

SNPs are considered to be closely related to environmental, ethnic, and hereditary factors [27]. Specific SNP sites could result in changes in gene expression profiles and gene functions, thereby influencing susceptibilities to various diseases, including tumors [28, 29]. Genome-wide association studies (GWASs) and high-throughput sequencing have identified numerous unstudied SNP sites, providing new perspectives on the mechanisms of tumorigenesis and methods for cancer screening [30, 31]. Additionally, there have been publications linking genetic polymorphisms of TCF21 to numerous diseases [16, 32]. Considering the above background, we focused on the association between polymorphisms of TCF21 and GC. After conducting statistical analyses, allele C in rs2327430 was found to be associated with a lower incidence of GC, and the TC + CC genotype was associated with a better prognosis.

As reported, SNPs occurring at any location in the genome have the potential to influence the trait. For instance, in the coding DNA sequence (CDS) region, nonsynonymous mutations could alter the amino acid encoded by influencing codons, while synonymous mutations may affect RNA secondary structures, protein folding, and cellular localization. Mutations located in UTRs play a role in regulating the expression of related genes, highlighting the significant implications of studying SNP sites in different parts of the genome [33,34,35]. SNP positioning in promoter regions can directly modulate gene transcriptional activity, attracting considerable attention from researchers [36, 37]. We took notice of the effects on binding with a transcription factor, mainly because rs2327430 is positioned 214 bp upstream of the first exon of TCF21. The conclusion that TFAP2A may act as a transcriptional repressor in this process was tentatively proposed by us. Although little is known about the mechanism by which transcriptional repressors exert their effects, while much is known about transcriptional activators, numerous studies have revealed the function of allele-specific transcriptional repressor binding sites [32, 38,39,40], which provided a theoretical basis for this study.

During the statistical analysis, the association of the C allele in rs2327430 with GC seemed to be driven by the heterozygous genotype, whereas the homozygous CC genotype had an OR of 1.26. However, the sample size of patients with the homozygous CC genotype may have resulted in a discrepancy between the direction of the correlations for the heterozygous TC and homozygous CC genotypes. The observation of a rare association between polymorphisms and clinicopathological features in the stratified analysis was noteworthy. Interestingly, one study recently reported a similar phenomenon [41]. This phenomenon may be attributed to many factors in vivo that differ from those in vitro. First, TCF21 is part of a complex regulatory network, while tumors are regulated by multiple molecules and axes. Additionally, the heterogeneity of the tumor may contribute to the limited impact of SNPs on clinicopathological characteristics at the population level. A larger sample size and a wider population range may be beneficial to reducing the effect of individual differences. Furthermore, strong immune-mediated surveillance and clearance could play an unpredictable role in vivo [42, 43]. Moreover, the interindividual variation in tumor metabolism and the complex tumor microenvironment in vivo may weaken the efficacy of SNPs [44, 45]. The protective effects on GC patients’ overall survival may indicate that these SNPs not only inhibit malignant properties but also alter the subsequent chemoradiotherapy sensitivity of GC [46, 47].

The mechanism by which TCF21 affects malignant behaviors in GC is also noteworthy. TCF21 has been reported to influence multiple pathways, including the PI3K-AKT pathway and MAPK pathway [48, 49]. However, further validation is needed to identify the direct target genes of TCF21 that regulate these downstream signaling pathways. This validation is of great value for future research.

Despite the significant discoveries revealed by our study, we must acknowledge its limitations. Above all, the modes by which transcriptional repressors work are diverse; some compete against transcriptional activators, and others may mediate chemical modifications. However, the mechanism of the transcriptional inhibition of TCF21 exerted by TFAP2A awaits further exploration [50]. Moreover, the reason that rs4896011 is associated with the occurrence but not the development of GC is still unclear. In addition, H. pylori infection, one of the major risk factors for GC in China, was not included in our study due to a lack of relevant medical records for a significant number of patients. Finally, the estimates of MST and the analysis of SNP correlation with clinicopathological data were imperfect due to limitations imposed by the small sample size and missing pathological information and survival results.

Data availability

No datasets were generated or analysed during the current study.

References

  1. Shen K, Liu KH, Wang YH, Ni PD, Xiao J, Hao F, Zhou XY, Xu ZK, Yang L. Polymorphisms of an oncogenic gene, mesothelin, predict the risk and prognosis of gastric cancer in a Chinese Han population. Arch Toxicol. 2022;96(7):2097–111.

    Article  CAS  PubMed  Google Scholar 

  2. Chen WQ, Zheng RS, Baade PD, Zhang SW, Zeng HM, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA-Cancer J Clin. 2016;66(2):115–32.

    Article  PubMed  Google Scholar 

  3. Smyth EC, Nilsson M, Grabsch HI, van Grieken NCT, Lordick F. Gastric cancer. Lancet. 2020;396(10251):635–48.

    Article  CAS  PubMed  Google Scholar 

  4. Ding P, Zheng C, Cao GH, Gao ZM, Lei YY, Deng P, Hou B, Li K. Combination of preoperative plasma fibrinogen and AJCC staging improves the accuracy of survival prediction for patients with stage I-II gastric cancer after curative gastrectomy. Cancer Med. 2019;8(6):2919–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Hundahl SA. Staging, stage migration, and patterns of spread in gastric cancer. Semin Radiat Oncol. 2002;12(2):141–9.

    Article  PubMed  Google Scholar 

  6. Briones-Orta MA, Avendano-Vazquez S, Aparicio-Bautista DI, Coombes JD, Weber GF, Syn WK. Osteopontin splice variants and polymorphisms in cancer progression and prognosis. Biochim Biophys Acta-Rev Cancer. 2017;1868(1):93–U347.

    Article  CAS  PubMed  Google Scholar 

  7. Jiao JH, Wu JY, Huang DS, Liu L. Lack of association of the iNOS gene polymorphism with risk of cancer: a systematic review and Meta-analysis. Sci Rep. 2015;5:7.

    Article  Google Scholar 

  8. Gao CD, Zhuang J, Zhou C, Li HY, Liu C, Liu LJ, Feng FB, Liu RJ, Sun CG. SNP mutation-related genes in breast cancer for monitoring and prognosis of patients: a study based on the TCGA database. Cancer Med. 2019;8(5):2303–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Xu R, Peng CB, Xiao SM, Zhuang W. IFRD1 polymorphisms and gastric cancer risk in a Chinese population. Med Oncol. 2014;31(9):4.

    Article  Google Scholar 

  10. Luo DK, Wang YN, Huan XK, Huang C, Yang C, Fan H, Xu ZK, Yang L. Identification of a synonymous variant in TRIM59 gene for gastric cancer risk in a Chinese population. Oncotarget. 2017;8(7):11507–16.

    Article  PubMed  Google Scholar 

  11. Guo ZJ, Zhou K, Wang Q, Huang YS, Ji J, Peng Y, Zhang XY, Zheng TH, Zhang Z, Chong DC, et al. The transcription factor RUNX2 fuels YAP1 signaling and gastric cancer tumorigenesis. Cancer Sci. 2021;112(9):3533–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zhao X, Chen LL, Wu JX, You J, Hong QQ, Ye F. Transcription factor KLF15 inhibits the proliferation and migration of gastric cancer cells via regulating the TFAP2A-AS1/NISCH axis. Biol Direct. 2021;16(1):11.

    Article  Google Scholar 

  13. Gooskens SL, Klasson TD, Gremmels H, Logister I, Pieters R, Perlman EJ, Giles RH, van den Heuvel-Eibrink MM. TCF21 hypermethylation regulates renal tumor cell clonogenic proliferation and migration. Mol Oncol. 2018;12(2):166–79.

    Article  CAS  PubMed  Google Scholar 

  14. Hussain A, Voisin V, Poon S, Karamboulas C, Bui NHB, Meens J, Dmytryshyn J, Ho VW, Tang KH, Paterson J, et al. Distinct fibroblast functional states drive clinical outcomes in ovarian cancer and are regulated by TCF21. J Exp Med. 2020;217(8):26.

    Article  Google Scholar 

  15. Yang Z, Jiang XD, Li DM, Dong QZ, Zhao HY, Jiang XF. TCF21 inhibits proliferation and chemoresistance through the AKT pathway in human gastric cancer. Gene. 2019;682:42–9.

    Article  CAS  PubMed  Google Scholar 

  16. Bakshi D, Nagpal A, Sharma V, Sharma I, Shah R, Sharma B, Bhat A, Verma S, Bhat GR, Abrol D, et al. MassARRAY-based single nucleotide polymorphism analysis in breast cancer of north Indian population. BMC Cancer. 2020;20(1):8.

    Article  Google Scholar 

  17. Pettersson FH, Anderson CA, Clarke GM, Barrett JC, Cardon LR, Morris AP, Zondervan KT. Marker selection for genetic case-control association studies. Nat Protoc. 2009;4(5):743–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Gopalakrishnan S, Qin ZS. TagSNP selection based on pairwise LD criteria and power analysis in association studies. Pacific Symposium on Biocomputing Pacific Symposium on Biocomputing 2006:511–522.

  19. Wang YN, Zhang F, Mukiibi R, Chen LH, Vinsky M, Plastow G, Basarab J, Stothard P, Li CX. Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: II: carcass merit traits. BMC Genomics. 2020;21(1):22.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Zhang F, Wang YN, Mukiibi R, Chen LH, Vinsky M, Plastow G, Basarab J, Stothard P, Li CX. Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: I: feed efficiency and component traits. BMC Genomics. 2020;21(1):22.

    CAS  Google Scholar 

  21. Ben ST, Cooper-DeHoff RM, Flaten HK, Evero O, Ferrara TM, Spritz RA, Monte AA. Multiplex SNaPshot-a new simple and efficient CYP2D6 and ADRB1 genotyping method. Hum Genomics. 2016;10:9.

    Article  Google Scholar 

  22. Lai GQ, Zhang WL, Tang H, Zhao TT, Wei LW, Tao Y, Wang ZC, Huang AL. A SNaPshot assay for the rapid and simple detection of hepatitis B virus genotypes. Mol Med Rep. 2014;10(3):1245–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Lin LL, Xiao J, Shi L, Chen WW, Ge YG, Jiang MK, Li ZL, Fan H, Yang L, Xu ZK. STRA6 exerts oncogenic role in gastric tumorigenesis by acting as a crucial target of miR-873. J Exp Clin Cancer Res. 2019;38(1):16.

    Article  Google Scholar 

  24. Akama T, Chun TH. Transcription factor 21 (TCF21) promotes proinflammatory interleukin 6 expression and extracellular matrix remodeling in visceral adipose stem cells. J Biol Chem. 2018;293(17):6603–10.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Cao L, Wang SY, Zhang YQ, Wong KC, Nakatsu G, Wang XH, Wong S, Ji JF, Yu J. Zinc-finger protein 471 suppresses gastric cancer through transcriptionally repressing downstream oncogenic PLS3 and TFAP2Avol 37, pg 3601, (2018). Oncogene 2022, 41(23):3298–3299.

  26. Xiao J, Zhu CM, Ni PD, Chen WW, Li ZL, Fan H, Ma X, Xu ZK, Yang L. Correlations of SNHG5 genetic polymorphisms with susceptibility and prognosis to gastric cancer in a Chinese population. Genomics. 2021;113(4):1754–60.

    Article  CAS  PubMed  Google Scholar 

  27. Edwards DRV, Naj AC, Monda K, North KE, Neuhouser M, Magvanjav O, Kusimo I, Vitolins MZ, Manson JE, O’Sullivan MJ, et al. Gene-environment interactions and obesity traits among postmenopausal African-American and hispanic women in the women’s Health Initiative SHARe Study. Hum Genet. 2013;132(3):323–36.

    Article  PubMed  Google Scholar 

  28. Cheng TC, Gamage SMK, Hewage D, Lu CT, Aktar S, Gopalan V, Lam AK. AHR gene expression and the polymorphism rs2066853 are associated with clinicopathological parameters in colorectal carcinoma. Hum Pathol. 2022;122:50–9.

    Article  CAS  PubMed  Google Scholar 

  29. Cui J, Li GQ, Yin J, Li LW, Tan Y, Wei HR, Liu B, Deng LH, Tang JL, Chen YL, et al. GSTP1 and cancer: expression, methylation, polymorphisms and signaling (review). Int J Oncol. 2020;56(4):867–78.

    CAS  PubMed  Google Scholar 

  30. Klijn C, Settleman J, Seshagiri S, Zhang ZM. A comprehensive transcriptional portrait of human cancer cell lines. Cancer Res. 2015;75(22):2.

    Google Scholar 

  31. Wang MH, Cordell HJ, Van Steen K. Statistical methods for genome-wide association studies. Semin Cancer Biol. 2019;55:53–60.

    Article  PubMed  Google Scholar 

  32. Suriano AR, Sanford AN, Kim N, Oh M, Kennedy S, Henderson MJ, Dietzmann K, Sullivan KE. GCF2/LRRFIP1 represses tumor necrosis factor alpha expression. Mol Cell Biol. 2005;25(20):9073–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Chen R, Davydov EV, Sirota M, Butte AJ. Non-synonymous and synonymous coding SNPs show similar likelihood and effect size of Human Disease Association. PLoS ONE. 2010;5(10):6.

    Article  Google Scholar 

  34. Kim S, Misra A. SNP genotyping: technologies and biomedical applications. Annu Rev Biomed Eng. 2007;9:289–320.

    Article  CAS  PubMed  Google Scholar 

  35. Lukas J, Hermann A, Giese AK. Silent but significant - A synonymous SNV alters prognosis in pompe disease. EBioMedicine. 2019;43:20–1.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Madureira MWS, Queiroz MAF, Lima SS, Pereira LMS, da Costa CA, de Sousa MS, Feitosa RNM, Monteiro JC, Ishak R, Vallinoto ACR et al. The FOXP3-924A/G Single Nucleotide Polymorphism May Be Associated with Predictive Factors for Human T Lymphotropic Virus 1 Associated Myelopathy. Viral immunology 2023.

  37. Wang Y, Ma RN, Liu B, Kong JY, Lin HY, Yu X, Wang RY, Li L, Gao M, Zhou BS, et al. SNP rs17079281 decreases lung cancer risk through creating an YY1-binding site to suppress DCBLD1 expression. Oncogene. 2020;39(20):4092–102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Cavalli M, Pan G, Nord H, Arzt EW, Wallerman O, Wadelius C. Genetic prevention of hepatitis C virus-induced liver fibrosis by allele-specific downregulation of MERTK. Hepatol Res. 2017;47(8):826–30.

    Article  CAS  PubMed  Google Scholar 

  39. Ma SL, Tang NLS, Tam CWC, Lui VWC, Lam LCW, Chiu HFK, Driver JA, Pastorino L, Lu KP. A PIN1 polymorphism that prevents its suppression by AP4 associates with delayed onset of Alzheimer’s disease. Neurobiol Aging. 2012;33(4):804–13.

    Article  CAS  PubMed  Google Scholar 

  40. Syddall CM, Reynard LN, Young DA, Loughlin J. The identification of Trans-acting factors that regulate the expression of GDF5 via the Osteoarthritis susceptibility SNP rs143383. PLoS Genet. 2013;9(6):14.

    Article  Google Scholar 

  41. Lv JC, Song Q, Bai KX, Han J, Yu H, Li K, Zhuang JT, Yang X, Yang HW, Lu Q. N6-methyladenosine-related single-nucleotide polymorphism analyses identify oncogene RNFT2 in bladder cancer. Cancer Cell Int. 2022;22(1):14.

    Article  Google Scholar 

  42. Sattler S. The Role of the Immune System Beyond the Fight Against Infection. In: Immunology of Cardiovascular Homeostasis and Pathology Edited by Sattler S, KennedyLydon T, vol. 1003. Cham: Springer International Publishing Ag; 2017: 3–14.

  43. Certo M, Tsai CH, Pucino V, Ho PC, Mauro C. Lactate modulation of immune responses in inflammatory versus tumour microenvironments. Nat Rev Immunol. 2021;21(3):151–61.

    Article  CAS  PubMed  Google Scholar 

  44. Li JQ, Eu JQ, Kong LR, Wang LZ, Lim YC, Goh BC, Wong ALA. Targeting metabolism in Cancer cells and the Tumour Microenvironment for Cancer Therapy. Molecules. 2020;25(20):40.

    Article  Google Scholar 

  45. Shi J, Wu Z, Wu X, Huangfu L, Guo T, Cheng X, Han J, Li Z, Xing X, Ji J. Characterization of glycometabolism and tumor immune microenvironment for predicting clinical outcomes in gastric cancer. iScience. 2023;26(3):106214.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Bukowska-Strakova K, Wlodek J, Pitera E, Kozakowska M, Konturek-Ciesla A, Ciesla M, Gonka M, Nowak W, Wieczorek A, Pawinska-Wasikowska K, et al. Role of HMOX1 promoter genetic variants in Chemoresistance and Chemotherapy Induced Neutropenia in Children with Acute Lymphoblastic Leukemia. Int J Mol Sci. 2021;22(3):22.

    Article  Google Scholar 

  47. Zhang ZH, Liu XP, Li L, Yang Y, Yang JG, Wang Y, Wu JJ, Wu XD, Shan L, Pei F, et al. SNP rs4971059 predisposes to breast carcinogenesis and chemoresistance via TRIM46-mediated HDAC1 degradation. Embo J. 2021;40(19):22.

    Article  CAS  Google Scholar 

  48. Duan HX, Li BW, Zhuang X, Wang LT, Cao Q, Tan LH, Qu GF, Xiao S. TCF21 inhibits tumor-associated angiogenesis and suppresses the growth of cholangiocarcinoma by targeting PI3K/Akt and ERK signaling. Am J Physiol-Gastroint Liver Physiol. 2019;316(6):G763–73.

    Article  CAS  Google Scholar 

  49. Zhang XY, Cheng BH, Jiang HX, Liu C, Cao ZP, Luan P, Wang N, Li H. Transcription factor 21 promotes Chicken Adipocyte differentiation at least in part via activating MAPK/JNK signaling. Genes. 2021;12(12):11.

    Article  CAS  Google Scholar 

  50. Elowitz MB, Leibler S. A synthetic oscillatory network of transcriptional regulators. Nature. 2000;403(6767):335–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was financially funded by the National Natural Science Foundation of China (Grant No. 81874219), Jiangsu Province Capability Improvement Project through Science, Technology and Education (Jiangsu Provincial Medical Key Discipline, ZDXK202222), and the Project of Cultivating Innovation in Science and Technology Plan of Liyang City (Grant No. LC2021001).

Author information

Authors and Affiliations

Authors

Contributions

Xinyi Zhou: Conceptualization, Methodology, Software, Writing-Reviewing and Editing. Kuan Shen: Data curation, Writing-Original. Shuqing Cao: draft preparation. Pengyu Li: Visualization. Jian Xiao: Investigation. Jiacheng Dong: Data curation. Quan Cheng: Software. Li Hu: Validation. Zekuan Xu: Project administration. Li Yang: Supervision, Resources, Funding acquisition.

Corresponding author

Correspondence to Li Yang.

Ethics declarations

Ethics approval and consent to participate

All subjects have signed informed consent forms, following the ethical approval of the local ethical committee.

Consent for publication

All authors have agreed with the publication of this article.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, X., Shen, K., Cao, S. et al. Polymorphism rs2327430 in TCF21 predicts the risk and prognosis of gastric cancer by affecting the binding between TFAP2A and TCF21. Cancer Cell Int 24, 159 (2024). https://doi.org/10.1186/s12935-024-03343-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-024-03343-z

Keywords