- Research article
- Open Access
- Open Peer Review
Characterization and risk association of polymorphisms in Aurora kinases A, B and C with genetic susceptibility to gastric cancer development
BMC Cancer volume 19, Article number: 919 (2019)
Single nucleotide polymorphisms (SNPs) in genes encoding mitotic kinases could influence development and progression of gastric cancer (GC).
Case-control study of nine SNPs in mitotic genes was conducted using qPCR. The study included 116 GC patients and 203 controls. In silico analysis was performed to evaluate the effects of polymorphisms on transcription factors binding sites.
The AURKA rs1047972 genotypes (CT vs. CC: OR, 1.96; 95% CI, 1.05–3.65; p = 0.033; CC + TT vs. CT: OR, 1.94; 95% CI, 1.04–3.60; p = 0.036) and rs911160 (CC vs. GG: OR, 5.56; 95% CI, 1.24–24.81; p = 0.025; GG + CG vs. CC: OR, 5.26; 95% CI, 1.19–23.22; p = 0.028), were associated with increased GC risk, whereas certain rs8173 genotypes (CG vs. CC: OR, 0.60; 95% CI, 0.36–0.99; p = 0.049; GG vs. CC: OR, 0.38; 95% CI, 0.18–0.79; p = 0.010; CC + CG vs. GG: OR, 0.49; 95% CI, 0.25–0.98; p = 0.043) were protective. Association with increased GC risk was demonstrated for AURKB rs2241909 (GG + AG vs. AA: OR, 1.61; 95% CI, 1.01–2.56; p = 0.041) and rs2289590 (AC vs. AA: OR, 2.41; 95% CI, 1.47–3.98; p = 0.001; CC vs. AA: OR, 6.77; 95% CI, 2.24–20.47; p = 0.001; AA+AC vs. CC: OR, 4.23; 95% CI, 1.44–12.40; p = 0.009). Furthermore, AURKC rs11084490 (GG + CG vs. CC: OR, 1.71; 95% CI, 1.04–2.81; p = 0.033) was associated with increased GC risk. A combined analysis of five SNPs, associated with an increased GC risk, detected polymorphism profiles where all the combinations contribute to the higher GC risk, with an OR increased 1.51-fold for the rs1047972(CT)/rs11084490(CG + GG) to 2.29-fold for the rs1047972(CT)/rs911160(CC) combinations. In silico analysis for rs911160 and rs2289590 demonstrated that different transcription factors preferentially bind to polymorphic sites, indicating that AURKA and AURKB could be regulated differently depending on the presence of particular allele.
Our results revealed that AURKA (rs1047972 and rs911160), AURKB (rs2241909 and rs2289590) and AURKC (rs11084490) are associated with a higher risk of GC susceptibility. Our findings also showed that the combined effect of these SNPs may influence GC risk, thus indicating the significance of assessing multiple polymorphisms, jointly. The study was conducted on a less numerous but ethnically homogeneous Bosnian population, therefore further investigations in larger and multiethnic groups and the assessment of functional impact of the results are needed to strengthen the findings.
Gastric cancer (GC) represents one of the major causes of tumor-linked death, with geographical and ethnical variations in incidence . Accurate chromosomal segregation in rapidly dividing tumor cells and defects during the spindle assembly checkpoint may contribute to tumorigenesis . Genetic alterations in mitotic genes could enhance susceptibility to malignant transformation through modifications of gene expression profiles [3, 4]. Aurora kinases are members of serine-threonine kinases family essential for cell cycle control . Aurora kinase A (AURKA) is involved in regulation of a several oncogenic signaling processes, including mitotic entry, cytokinesis, functions of centrosome, chromosome segregation, and chromosome alignment [6, 7]. Aurora kinase B (AURKB) assists in chromatin modification, spindle checkpoint regulation, cytokinesis and plays a significant role in establishment of the correct kinetochore/microtubule binding . Aurora kinase C (AURKC) acts as a chromosomal passenger protein, participating in the proper centrosome functioning . Polo-like kinase 1 (PLK1) is essential for cell division and regulates various cellular events including centrosome maturation, mitotic checkpoint activation, spindle assembly, kinetochore/microtubule attachment, exit from the mitosis, and cytokinesis .
In this study, using a case-control approach, we estimated the impact of rs2273535, rs1047972, rs911160 and rs8173 in AURKA, rs2241909 and rs2289590 in AURKB, rs758099 and rs11084490 in AURKC and rs42873 in PLK1 mitotic checkpoint genes on GC susceptibility in Bosnia and Herzegovina population. In addition, the associations between single nucleotide polymorphisms and the histological types of gastric cancer (intestinal and diffuse types) have been investigated. By conducting in silico analysis of SNPs, we evaluated the impact of the studied polymorphisms in introns and untranslated regions (UTRs) within candidate genes (AURKA, AURKB, AURKC and PLK1) on transcription factors binding sites.
Study design and populations
Our examined population consisted of 116 GC patients with diagnosed gastric adenocarcinoma from the Clinical Pathology and Cytology at the University Clinical Center Sarajevo, Bosnia and Herzegovina. General status of gastric cancer patients is given in Table 1. Gastric cancer patients in the case group were not subjected to any type of treatment (radiotherapy or chemotherapy).The formalin fixed paraffin embedded (FFPE) cancer tissue sections were collected during surgical procedures. Simultaneously, 203 healthy blood donors (controls) of Bosnian origin (matched to cases for ethnicity) were randomly selected and signed up for the present study. Individuals in the control group had no history of any neoplastic formation, were not related to each other and to the patients group. Three ml of blood was sampled from each control individual and stored at − 80 °C. The study was approved by the Ethical Committee at the University Clinical Centre Sarajevo (No. 0302–36,765). Personal information was encrypted to provide maximum anonymity in compliance with the Helsinki Declaration.
Genomic DNA from FFPE GC tissues was isolated using the Chemagic FFPE DNA Kit special (PerkinElmer Inc., Waltham, MA, USA), according to manufacturer’s recommendations. Automated DNA washing and elution was conducted on Chemagic Magnetic Separation Module I robot (PerkinElmer Inc., Waltham, MA, USA), following manufacturer’s standard programme. All sample transfers were performed with 4-eye principle to avoid sample mix-ups. DNA from lymphocytes (control DNA) was extracted using the Promega™ Wizard™ Genomic DNA Purification Kit Protocol (Promega Corp., Fitchburg, WI, USA), in concordance with the manufacturer’s recommendations. The qualitative and quantitative analysis of extracted DNA was conducted by use of the DropSense96 photometer (Trinean, Gentbrugge, Belgium) and Synergy™ 2 Multi Mode Reader (BioTek, Inc., Winooski, VT, USA).
Selection of polymorphisms
We selected nine polymorphisms in mitotic genes, namely rs2273535, rs1047972, rs911160 and rs8173 (AURKA), rs2241909 and rs2289590 (AURKB), rs758099 and rs11084490 (AURKC) and rs42873 (PLK1). The positions of selected genetic variants in mitotic genes are presented in Fig. 1. For this purpose, gene structures were extracted from the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank (PDB) . Selection of the polymorphisms for this study was conducted in accordance with the parameters described below: (a) previously demonstrated association with respect to certain cancer types; (b) minor allele frequency (MAF) of less than or equal to 10% in the population of Utah residents with Northern and Western European ancestry (CEU), as stated by the Phase 31,000 Genomes; and (c) tagging polymorphisms (tagSNPs) status, which was anticipated in silico by use of LD Tag Selection of SNP (tagSNP) (https://snpinfo.niehs.nih.gov) , with the following parameters: 1 kb of the sequences upstream/downstream from gene was selected, linkage disequilibrium (LD) lower limit of 0.8, and MAF range 0.05–0.5 for CEU subpopulation (Table 2 and Fig. 2).
Genotyping was conducted using TaqMan SNP genotyping assays (Applied Biosystems, Foster City, CA). The assay ID numbers are presented in Table 2. The reaction mixtures, GC samples (5 μl) and controls (10 μl), were composed of 20X TaqMan® assay with 2X Master Mix (Applied Biosystems, Foster City, CA), and 20 nanograms of DNA. The polymerase chain reaction (PCR) profile was carried out in concordance with the manufacturer’s recommendations (Initial denaturation at 95 °C for 10 min, 45 cycles at 92 °C for 15 s and 60 °C for 90 s, using the ViiA 7 Real Time PCR System (Applied Biosystems, Foster City, CA). In each plate, at least two negative controls were included. PCR results were analyzed using TaqMan® Genotyper Software (Applied Biosystems, Foster City, CA, USA).
The genotype frequencies of the investigated variants were tested for Hardy-Weinberg equilibrium (HWE) in the case/control groups separately, using Michael H. Court’s online HWE calculator (http://www.tufts.edu) . The differences in genotype frequencies amongst GC cases and controls were calculated by use of the Chi-square test or Fisher’s exact test. Association between examined polymorphisms and the GC risk was estimated by multinomial logistic regression. Odds ratio (OR) with 95% confidence interval (CI) were computed in order to evaluate the relative risk. For the assessment of each genotype, risk estimates were computed for dominant, overdominant and recessive models using the most frequent homozygote as the reference. Akaike information criterion (AIC) was calculated to define which of the models best fits the data. A combined analysis was performed to evaluate synergistic effect of the studied polymorphisms. All statistical calculations were conducted using SPSS 20.0 software package (SPSS, Chicago, IL, USA). P ≤ 0.05 was chosen as threshold value in significance testing. MAF plot was created by use of the PAST software package, version 3.18 (http://folk.uio.no/ohammer/past/) .
Determination of the haplotype block structure and haplotype analysis, which encompassed subsequent corrections for multiple comparisons by 10,000 permutations, were evaluated using the Haploview software, version 4.2 ., and SNP tools V1.80 (MS Windows, Microsoft Excel). To construct the haplotype block, the solid spine of the linkage disequilibrium algorithm with a minimum Lewontin’s D′ value of 0.8 was selected.
In silico analysis of SNPs
Impact of the polymorphic DNA sequences (SNPs in introns and untranslated regions (UTRs)) on transcription factors binding sites (TFBSs) was estimated in silico. Bioinformatic functional evaluation was carried out using PROMO software (ALGGEN web-server), which is utilizing data from TRANSFAC database V8.3 [15, 16]. FASTA sequences for the investigated genetic variants were downloaded from Ensembl 90 (www.ensembl.org/index.html) . Identification of transcription factor binding sites was performed with the following criteria: human species, all sites and factors.
Genotype distributions for examined SNPs
For all of the 9 studied variants, rs2273535 (AURKA), rs1047972 (AURKA), rs911160 (AURKA), rs8173 (AURKA), rs2241909 (AURKB), rs2289590 (AURKB), rs758099 (AURKC), rs11084490 (AURKC), rs42873 (PLK1) was determined to be in HWE in both, case and control populations (P > 0.05). When chi-square test and Fisher exact test were conducted for the frequency distributions at the genotypic level, a significant differences for rs911160 in AURKA (P = 0.044), rs8173 in AURKA (P = 0.018), rs2289590 in AURKB (P < 0.001) and rs11084490 in AURKC (P = 0.009) between the cases and controls for all types of GC were observed (summarized in Table 3).
Effect of studied polymorphisms on gastric cancer risk
Patients with rs1047972 (AURKA) CT genotype had a higher risk of GC development in comparison with the reference CC genotype (OR = 1.96, 95% CI = 1.05–3.65, P = 0.033) (Table 4). Genotypes (TT + CT) vs. reference CC genotype in dominant model (OR = 1.92, 95% CI = 1.06–3.48, P = 0.030) and CT vs. reference (CC + TT) genotypes in overdominant model (OR = 1.94, 95% CI = 1.04–3.60, P = 0.036) were associated with higher disease risk (Table 4). Based on Akaike information criterion (AIC), the overdominant model was selected as the model that best fits the data. The rs911160 (AURKA) CC genotype was positively associated with an increased gastric cancer risk in comparison with the reference GG genotype (OR = 5.56, 95% CI = 1.24–24.81, P = 0.025). Also, CC genotype was associated with disease risk in the recessive genetic model (GG + CG) vs. CC genotypes, (OR = 5.26, 95% CI = 1.19–23.22, P = 0.028). However, the confidence intervals in those two cases were wide; therefore, other factors might play a significant role in GC risk in interaction with this polymorphism. Comparison of genotype distributions for rs8173 (AURKA) showed that patients with GG genotype (OR = 0.38, 95% CI = 0.18–0.79, P = 0.010), and CG genotype (OR = 0.60, 95% CI = 0.36–0.99, P = 0.049) had decreased risk of gastric cancer. Analysis of genetic models showed that GG + CG genotypes in comparison with the reference CC genotype in dominant model (OR = 0.54, 95% CI = 0.33–0.87, P = 0.012) and GG vs. reference (CC + CG) (OR = 0.49, 95% CI = 0.25–0.98, P = 0.043) genotypes (recessive genetic model) were associated with decreased GC risk. According to the calculated AIC values, (CC + CG):GG recessive model had more statistical power than dominant model CC:(GG + CG). Analysis of rs2241909 (AURKB) demonstrated that G allele (dominant model: (GG + AG) vs. common AA genotype) was associated with higher risk of GC development (OR = 1.61, 95%CI = 1.01–2.56, P = 0.041). Comparison of the reference AA genotype with AC (OR = 2.41, 95% CI = 1.47–3.98, P = 0.001) and CC (OR = 6.77, 95% CI = 2.24–20.47, P = 0.001) genotypes of rs2289590 (AURKB) also revealed a significant effect of these two genotypes on increased GC risk. CC and AC genotypes in dominant model (OR = 2.78, 95% CI = 1.71–4.51, P < 0.001) as well as CC genotype in recessive model (OR = 4.23, 95% CI = 1.44–12.40, P = 0.009) and AC genotype in overdominant genetic model (OR = 1.77, 95% CI = 1.10–2.85, P = 0.017) were associated with an elevated disease risk. Since recessive genetic model had the lowest AIC value, when compared to the dominant and overdominant models, it was considered to be preferred model. However, in this model the confidence interval was wide, therefore, other factors could influence its effect. For rs11084490 (AURKC) polymorphism, (GG + CG) vs. CC genotypes in dominant model demonstrated statistically significant effect on higher GC risk (OR = 1.71, 95% CI = 1.04–2.81, P = 0.033). Additionally, the five polymorphisms rs1047972, rs911160, rs2241909, rs2289590 and rs11084490, associated with an increased GC risk individually in this study, were subjected to the combined analysis in order to determine polymorphism profiles related to the higher risk of this disease. The results of the synergistic effects of these SNPs are summarized in Table 5. By analyzing various combinations of risk genotypes (two to five combined SNPs), we demonstrated that the additive effect of all combinations significantly affected the risk of GC development, with an odds ratio ranging from (OR = 1.51, 95% CI = 1.03–2.22, P = 0.034) for the combined rs1047972(CT)/rs11084490(CG + GG) risk genotypes to (OR = 2.29, 95% CI = 1.32–3.97, P = 0.003) for the rs1047972(CT)/rs911160(CC) combination. Another interesting combined effect was demonstrated for five-polymorphisms combination rs1047972(CT)/rs911160 (CC)/rs2241909 (AG + GG)/rs2289590(AC + CC)/rs11084490 (CG + GG). In this case, this combination was significantly associated with an increased GC risk (OR = 1.83 95% CI = 1.46–2.29, P < 0.001). No significant effects on gastric cancer susceptibility were revealed for rs2273535 (AURKA), rs758099 (AURKC) and rs42873 (PLK1) polymorphisms (P > 0.05), when patients with both types of GC, intestinal and diffuse, were taken into account.
Next, we estimated the effects of genotypes on GC subtypes (presented in Table 4). CT genotype of rs1047972 (AURKA) was more frequent in patients with intestinal type (OR = 2.53, 95% CI = 1.02–6.30, P = 0.045) in comparison with the reference CC genotype. Likewise, (TT + CT) genotypes vs. reference CC (OR = 2.39, 95% CI = 1.02–5.63, P = 0.045) and CT vs. common (CC + TT) genotypes (OR = 2.50, 95%CI = 1.01–6.22, P = 0.047) were associated with higher risk for the development of intestinal subtype. According to the AIC values, (CC + TT):CT overdominant genetic model displayed stronger statistical confidence than dominant model CC:(TT + CT). The rs8173 (AURKA), GG genotype, in comparison with the reference CC genotype, was underrepresented in patients with diffuse GC type (OR = 0.32, 95% CI = 0.13–0.77, P = 0.012). Furthermore, both (GG + CG) genotypes as compared to its common CC genotype in dominant model (OR = 0.49, 95% CI = 0.27–0.89, P = 0.021) and GG vs. reference (CC + CG) genotypes in recessive model (OR = 0.44, 95% CI = 0.20–0.98, P = 0.044) were associated with the decreased diffuse type GC risk. In order to discriminate between these two competing models, in accordance with AIC, recessive model represents the preferred model in comparison with the dominant model. In stratified analysis for rs2241909 (AURKB), we found that carriers of AG genotype had elevated risk of developing intestinal type GC as compared to its reference AA genotype (OR = 2.23, 95% CI = 1.16–4.27, P = 0.016). Carriers of (GG + AG) genotypes had more frequently intestinal type of GC when compared to the carriers of the more common AA genotype in dominant model (OR = 2.38, 95% CI = 1.27–4.46, P = 0.007). In overdominant model (OR = 1.93, 95%CI = 1.02–3.67, P = 0.042) individuals with AG genotype had more frequently intestinal type GC in comparison with reference genotypes (AA+GG). According to the calculated AIC values, overdominant model had more statistical power than dominant, therefore it represents the model that better fitted the data. The higher risk for intestinal type GC development was also detected for the patients with CC genotype of rs2289590 (AURKB) (OR = 5.19, 95% CI = 1.14–23.56, P = 0.033). Dominant genetic model revealed that patients with (CC + AC) genotypes when compared to the AA genotype (OR = 2.04, 95% CI = 1.07–3.88, P = 0.028) had significantly more frequently intestinal GC subtype. AC genotype (OR = 3.12, 95% CI = 1.68–5.80, P < 0.001) was more frequently observed in patients with diffuse subtype. Regarding genetic models, (CC + AC) genotypes in dominant model (OR = 3.58, 95% CI = 1.96–6.52, P < 0.001), CC genotype in recessive model (OR = 4.72, 95%CI = 1.09–20.43, P = 0.038) and AC genotype in overdominant model (OR = 2.27, 95% CI = 1.24–4.13, P = 0.007) were associated with the increased risk of diffuse subtype, with a recessive model being the one that best suited the data (according to the AIC value), however, the confidence interval in this model was also the largest. For rs11084490 (AURKC), dominant model (GG + CG) vs. CC (ref.) genotypes reveled a significant effect of GG and CG genotypes on the higher risk of intestinal subtype (OR = 2.03, 95% CI = 1.02–4.04, P = 0.043).
For genotypes of rs2273535 (AURKA), rs911160 (AURKA), rs758099 (AURKC) and rs42873 (PLK1) no significant effect on any of the GC histological subtypes was noted (P > 0.05).
Raw genotyping data for the studied polymorphisms rs2273535, rs1047972, rs911160 and rs8173 in AURKA gene were used to perform haplotype analysis. Using the Haploview software, our results showed that no haplotype block was created with an average Lewontin’s D < 0.8 (Fig. 3) thus, no haplotypes were available for the analysis of their potential association with GC risk.
Bioinformatic SNP analysis
Our in silico analysis suggested that polymorphic sequences in transcription factors binding sites (TFBSs), within AURKA, AURKB, AURKC and PLK1 genes, bind various transcription factors (TFs). In this regard, the region comprising G allele of rs911160 in AURKA was linked with C/EBPalpha, C/EBPbeta and GR-beta proteins, whereas for C allele, additional binding sites for NF-Y, NFI-CTF and NF-1 were identified (Table 6). For rs2289590 in AURKB, an additional motif for YY1 binding was recognized when C allele was present. The region near C allele of rs758099 was associated with binding sites for NF-1, NF-Y, XBP-1, ENKTF-1, CTF, PEA3 and POU2F1, whereas in the presence of T allele NF-1, NF-Y, GATA-1 and TFII-I sequence-specific DNA-binding factors were recorded. Only in the case of rs11084490 in AURKC there were no changes in transcription factor binding site motif (XBP-1), if different alleles, either C or G, were present. The G allele of rs42873 in PLK1 was linked with an additional recognition motif for c-Jun transcription factor.
In this study, SNPs rs2273535, rs1047972, rs911160 and rs8173 (AURKA), rs2241909 and rs2289590 (AURKB), rs758099 and rs11084490 (AURKC), and rs42873 (PLK1) mitotic kinases were screened for associations with the genetic susceptibility to gastric cancer (GC) in Bosnian population. We also examined genotype effects of the investigated polymorphisms for each GC subtype.
In our study, a significant association between AURKA rs1047972 CT genotype with the overall GC susceptibility was found. Similarly, in stratified analysis established on Lauren’s classification , this genotype has affected intestinal GC subtype, whereas association was lost in patients with diffuse type GC. Furthermore, for rs911160 in AURKA, analysis showed that its CC genotype showed effect on increased disease risk. Our results also revealed that AURKA rs8173 GG genotype could be associated with a decreased GC risk. In stratified analysis of GC types, the association was significant in patients with the diffuse type GC. These findings could underlie different epidemiological and clinical patterns observed in intestinal and diffuse subtypes .
Bioinformatic analysis of transcription binding sites reveled that in the case of rs911160 C allele, an extra NF-Y, NFI-CTF and NF-1 transcription factors were detected in comparison with G allele. NF-Y regulates some of the genes enrolled in regulation of cell cycle, which are also deregulated in certain human diseases . NF-1 family of sequence-specific TFs affect the rate of transcription, either through repression or activation . NFI-CTF corresponds to the protein family involved in transcription activation, which is guided by the RNA polymerase II . Single nucleotide polymorphisms in TFBSs, can alter gene expression through linkage of different TFs, by removing existing or creating new binding motifs . Also, it has been demonstrated that introns, particularly long ones, harboring more functional cis-acting elements, could accommodate sites for binding several TFs, and consequently regulate transcription . Thus, our results suggest that rs911160 alleles in TFBS regions could bind various transcription factors which might affect the rate of AURKA expression, resulting in distinctions in exposure to the risk of GC development. In our previous study conducted in Slovenian population, we reported AURKA rs911160 association with an increased GC risk , and our findings from this study are supportive to these findings. Polymorphisms in 3′ untranslated regions (3’UTRs) of genes might affect mRNA stability, translation and overall level of post-transcriptional expression through effects on polyadenylation and/or changing binding sites for regulatory proteins as well as for microRNAs (miRNAs) . Recent study has demonstrated that 3’UTR variant in high mobility group box-1 (HMGB1) gene have a protective effect on overall survival in GC patients through decreased HMGB1 mRNA expression levels . Thus, it is reasonable to believe that protective effect of GG genotype of SNP rs8173 in AURKA 3’UTR, evaluated in our study, could be associated with an aberrant AURKA expression.
AURKA confers major contribution to the processes, such as centrosome duplication, entry into mitosis and in spindle assembly checkpoint . Several studies have suggested that AURKA overexpression leads to malignant transformation . A number of polymorphisms in the AURKA have also been reported to exhibit an effect on the risk of cancer onset. Genetic variant rs2273535 was associated with colorectal and lung cancer [29, 30]. In our study no significant association was observed between rs2273535 (AURKA) and GC risk. Polymorphism rs1047972, one of the most investigated variants in AURKA gene, showed significant association with the increased esophagus cancer risk as well as with gastric cancer risk and progression [31,32,33]. Our results from the present study confirm these previous findings. SNP rs1047972 might increase relative kinase activity of AURKA . AURKA is involved in phosphorylation of p53, which is followed by MDM2 induced degradation of p53, or resulting in silencing of the p53 transcriptional function . The absence of p53 can result in mitotic checkpoint dysfunction and subsequent chromosomal instability . Moreover, by suppressing p53 and p73 pro-apoptotic functions, AURKA enables a mechanism for cancer cells to evade apoptosis . Thus, it could be expected that slightly higher kinase activity could be involved in cancer development as well as cancer cell survival. In AURKA gene, rs1047972 and rs2273535 variants are located in exon 3 with high LD amongst them, suggesting that phenotypic effects of both polymorphisms could be consequence of a synergistic act. In addition, it was suggested that rs1047972 could possess a noticeable role in carcinogenesis by alteration of rs2273535 secondary structure and/or function . Our findings, regarding evaluated genetic variants in AURKA gene, suggest that rs1047972 and rs911160 polymorphisms could act as factors which contribute to GC susceptibility, whereas rs8173 variant might be protective factor for GC development.
Aurora kinase B (AURKB) is a subunit of chromosomal passenger complex (CPC), involved in the segregation of chromatids, cytokinesis and modification of histones  and has been overexpressed in different types of cancers encompassing prostate, thyroid and brain . It has been proposed that AURKB overexpression causes defects in chromosome segregation, aneuploidy and tumor development . We examined rs2241909 SNP in AURKB and found a significant association between (AG/GG) genotypes and increased susceptibility to GC. In addition to this, in analyses of genetic models, AG genotype demonstrated an effect on a higher risk of intestinal type GC growth. In an earlier study, rs2241909 showed association with familial breast cancer risk . The rs2241909 variant is a silent variant positioned on C terminal end of aurora kinase B. This amino acid change does not abolish or create splice site, nor affects exonic splicing enhancers/silencers motifs, and it has also been demonstrated that it does not change AURKB mRNA secondary structure . Therefore, the observed risk between GC risk and rs2241909 could be due to its linkage with another unidentified functional genetic variant. The analysis of the second polymorphism in AURKB, rs2289590, demonstrated that CC genotype was associated with higher risk of GC onset. In stratified analysis of GC types, both CC and AC genotypes had an effect on diffuse type GC risk, whereas CC genotype was related to the increased risk of developing intestinal GC subtype. In silico analysis of rs2289590 region revealed binding of additional YY1 transcription factor, if C allele was present.
The YY1 TF is associated with a cell cycle progression and it has been demonstrated that YY1 expression is with uncontrolled cell proliferation, apoptosis resistance and metastasis, thus acting as an initiator of carcinogenesis . Transcription factors (TFs) are important gene regulators with specific roles in cell cycle, thus when improperly regulated, they contribute to the failure in proper cellular functioning, instability and malignant transformation [41, 42]. SNPs in regulatory regions can moderate expression of genes through potential disruption of sequence specific DNA-binding motifs, which consequently alters the binding of the appropriate TFs . Our data for intronic rs2289590 in AURKB suggest that additional binding of the YY1 sequence-specific DNA-binding factor, when C allele is present within TF binding site, could modify AURKB expression level, which might result in higher susceptibility to gastric cancer occurrence. Important roles of introns in regulation of transcription have been reported in cell cycle and apoptosis genes, highlighting the significance of intronic genetic variants in tumorigenesis . More importantly, our findings from this study for rs2289590 (AURKB) association with an increased GC risk, are in accordance with the findings from our previous study conducted in Slovenian population .
Aurora kinase C (AURKC) represents a catalytic chromosomal passenger protein, similarly as Aurora kinase B, which plays essential role mitotic events, segregation and centrosome function throughout meiosis [8, 44]. AURKC overexpression has been described in malignant thyroid cell lines and tissues . It has been shown that overexpression of AURKC induces centrosome amplification, multinucleation and that its abnormal expression in somatic cells has an oncogenic potential . We examined rs11084490 in AURKC and its potential relationship with gastric cancer risk. A link between CG and GG genotypes and increased gastric cancer risk was observed. Stratified analyses revealed that these genotypes were more common in patients with intestinal type of GC. Polymorphism rs11084490 is situated within the 5’UTR region of AURKC. Eukaryotic 5’UTR various elements and structures e.g. hairpins, RNA G-quadruplexes (RG4s), Kozak sequences around the initiation codons, upstream open reading frames (uORFs) and start codons AUGs, internal ribosome entry sites (IRESs) and iron responsive elements (IREs) greatly influence mRNA translation . It has been demonstrated that 5′ uORF-altering polymorphisms and mutations significantly silence expression of the downstream protein . Additionally, genetic variations such as mutations and SNPs, by disrupting motifs within 5’UTR, are capable of causing damaging effects on human health, and could be associated with diseases such as multiple myeloma, esophageal cancer and many others . Therefore, observed association of the rs11084490 (AURKC) polymorphism with the increased GC risk in our study could be due to altered AURKC translation mediated by risk genotypes affecting the above mentioned functional motifs in AURKC 5’UTR. Our results demonstrated that rs758099 (AURKC) polymorphism exhibited no effect on GC susceptibility.
As reported above, the results of our study demonstrated involvement of the rs1047972 (AURKA), rs911160 (AURKA), rs2241909 (AURKB), rs2289590 (AURKB) and rs11084490 (AURKC) polymorphisms in gastric tumorigenesis. However, considering different genes included in chromosome segregation process, it is difficult to explain the association of gastric cancer development with an individual polymorphism. Therefore, a combined analysis spanning various gene polymorphisms enables the assessment of gene-gene interactions, and consequently determination of genetic profiles associated with a risk of GC.
In this study, combined analysis of the five polymorphisms and their risk genotypes associated with an increased susceptibility to gastric cancer, rs1047972(CT)/rs911160(CC)/rs2241909(AG + GG)/rs2289590(AC + CC)/rs11084490(CG + GG, revealed polymorphism profiles where all the combinations (two to five combined risk genotypes) influence the higher risk of GC, with an OR increased 1.51-fold for the rs1047972(CT)/rs11084490(CG + GG) to 2.29-fold for the rs1047972(CT)/rs911160(CC) combinations. It is also important to highlight that five-polymorphisms combination rs1047972(CT)/rs911160 (CC)/rs2241909(AG + GG)/rs2289590(AC + CC)/rs11084490 (CG + GG) showed significant effect on an increased GC risk (OR = 1.83 95%CI = 1.46–2.29, P < 0.001).
Several studies have conducted combined analysis of polymorphisms in gastric cancer. In one of them, it has been demonstrated that the risk of noncardia gastric cancer increased 27.3-fold with increasing number of proinflammatory genotypes for three or four polymorphisms . Similarly, another study revealed that combination of polymorphisms in genes involved in the inflammatory process could affect the increased risk of gastric cancer . These findings may be explained by an additive effect of the polymorphisms in inflammatory genes. Therefore, based on these results, we could assume that particular combinations of genetic variants in aurora kinases A, B and C, could act synergistically, in mediating aberrations in the process of chromosome segregation, leading to aneuploidy and consequently to gastric cancer development.
Polo-like kinase 1 (PLK1) is essential for cell division and it has been demonstrated that PLK1 with other signal proteins is responsible for mitotic progression and has also been linked to cellular proliferation . Moreover, it has been demonstrated that polymorphisms in PLK1 influence its expression, therefore they could potentially affect cancer risk and progression . We selected rs42873 (PLK1) polymorphism for the assessment of its possible effect on an increased gastric cancer risk, however, our results showed no significant association between rs42873 genetic variant and GC risk.
The results of this study revealed that AURKA (rs1047972 and rs911160), AURKB (rs2241909 and rs2289590) and AURKC (rs11084490) polymorphisms could affect the risk of gastric cancer, both individually and synergistically. Contrary, we found that AURKA (rs8173) polymorphism appeared to be associated with decreased GC risk. Collectively, these findings indicated the existence of the plausible roles of genetic variations in AURKA, AURKB and AURKC in stomach carcinogenesis. Our results could be beneficial in the further investigations of the functional impact of these polymorphisms. The present study is based on a reduced number of cases which represents its limitation, therefore it is important that larger prospective studies confirm our findings.
Availability of data and materials
The data used in this study contain personal information and are not publicly available, but can be requested from the Clinical Pathology and Cytology at the University Clinical Center Sarajevo, subject to ethical approvals.
Akaike information criterion
Algorithmics and genetics group
All phase 3 individuals
- AP-2 alpha A:
Activating enhancer binding Protein 2 alpha
- AURKA :
Aurora kinase A
- AURKB :
Aurora kinase B
- AURKC :
Aurora kinase C
- C/EBP alpha:
CCAAT/enhancer-binding protein alpha
- C/EBP beta:
CCAAT/enhancer-binding protein beta
Utah residents with Northern and Western European ancestry population
Transcription factor c-Jun
Chromosomal passenger complex
CCAAT box-binding transcription factor
Enkephalin transcription factor 1
Formalin fixed paraffin-embedded
GATA binding factor 1
Glucocorticoid receptor alpha
Glucocorticoid receptor beta
High mobility group box-1
Internal ribosome entry site
Logarithm of odds
Minor allele frequency
Mouse double minute 2 homolog
Nuclear factor 1
Nuclear factor I-CCAAT-binding transcription factor
Nuclear transcription factor Y
Tumor protein p53
Tumor protein p73
Polymerase chain reaction
Protein Data Bank
Polyoma enhancer activator 3
- PLK1 :
Polo-like kinase 1
POU domain, class 2, transcription factor 1
Research Collaboratory for Structural Bioinformatics
Single nucleotide polymorphism
Statistical package for social sciences
Thyroid hormone receptor beta 1
Tagging single nucleotide polymorphism
Transcription factor binding site
General transcription factor II-I
Upstream open reading frame
X-box binding protein 1
Yin yang 1 transcription factor
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
Vaclavicek A, Bermejo JL, Wappenschmidt B, Meindl A, Sutter C, Schmutzler RK, et al. Genetic variation in the major mitotic checkpoint genes does not affect familial breast cancer risk. Breast Cancer Res Treat. 2007;106:205–13.
Tomonaga T, Nomura F. Chromosome instability and kinetochore dysfunction. Histol Histopathol. 2007;22:191–7.
McLean MH, El-Omar EM. Genetics of gastric cancer. Nat Rev Gastroenterol Hepatol. 2014;11:664–74.
Glover DM, Leibowitz MH, McLean DA, Parry H. Mutations in aurora prevent centrosome separation leading to the formation of monopolar spindles. Cell. 1995;81:95–105.
Gavriilidis P, Giakoustidis A, Giakoustidis D. Aurora kinases and potential medical applications of Aurora kinase inhibitors: a review. J Clin Med Res. 2015;10:742–51.
Katsha A, Belkhiri A, Goff L, El-Rifai W. Aurora kinase a in gastrointestinal cancers: time to targer. Mol Cancer. 2015;14:106.
Sasai K, Katayama H, Stenoien DL, Fujii S, Honda R, Kimura M, et al. Aurora-C kinase is a novel chromosomal passenger protein that can complement Aurora-B kinase function in mitotic cells. Cell Motil Cytoskeleton. 2004;59:249–63.
Strebhardt K. Multifaceted polo-like kinases. Drug targets and antitargets for cancer therapy. Nat Rev Drug Discov. 2010;9:643–60.
Burley SK, Berman HM, Christie C, Duarte JM, Feng Z, Westbrook J, et al. RCSB protein data Bank: sustaining a living digital data resource that enables breakthroughs in scientific research and biomedical education. Protein Sci. 2018;27:316–30.
Xu Z, Taylor JA. SNPinfo: integrating GWAS and candidate gene information into functional SNP selection for genetic association studies. Nucleic Acids Res. 2009;37:W600–5.
Chahil JK, Munretnam K, Samsudin N, Lye SH, Hashim NAN, Ramzi NH, et al. Genetic polymorphisms associated with breast Cancer in Malaysian cohort. Indian J Clin Biochem. 2015;30:134–9.
Hammer Ø, Harper DAT, Ryan PD. PAST: paleontological statistics software package for education and data analysis. Palaeontol Electron. 2001;4:9pp.
Barrett JC. Haploview: Visualization and analysis of SNP genotype data. Cold Spring Harb Protoc. 2009;4:1-6.
Messeguer X, Escudero R, Farre D, Nunez O, Martinez J, Alba MM. PROMO: detection of known transcription regulatory elements using species-tailored searches. Bioinformatics. 2002;18:333–4.
Farre D, Roset R, Huerta M, Adsuara JE, Rosello L, Alba MM, et al. Identification of patterns in biological sequences at the ALGGEN server: PROMO and MALGEN. Nucleic Acids Res. 2003;31:3651–3.
Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, et al. Ensembl 2018. Nucleic Acids Res. 2018;46:D754–61.
Lauren P. The two histological main types of gastric carcinoma: diffuse and so called intestinal-type carcinoma. An attempt at a histo-clinical classification. Acta Pathol Microbiol Immunol Scand. 1965;64:31–9.
Vauhkonen M, Vauhkonen H, Sipponen P. Pathology and molecular biology of gastric cancer. Best Pract Res Clin Gastroenterol. 2006;20:651–74.
Ly LL, Yoshida H, Yamaguchi M. Nuclear transcription factor Y and its roles in cellular processes related to human disease. Am J Cancer Res. 2003;3:339–46.
Gronostajski RM. Roles of the NFI/CTF gene family in transcription and development. Gene. 2000;249:31–45.
Wenzelides S, Altmann H, Wendler W, Winnacker EL. CTFP-a new transcriptional activator of the NFI/CTF family. Nucleic Acids Res. 1996;24:2416–21.
Wang X, Tomso DJ, Liu X, Bell DA. Single nucleotide polymorphisms in transcriptional regulatory regions and expression of environmentally responsive genes. Toxicol Appl Pharmacol. 2005;207:84–90.
Li H, Chen D, Zhang J. Analysis of intron sequence features associated with transcriptional regulation in human genes. PLoS One. 2012;7:e46784.
Mesic A, Markocic E, Rogar M, Juvan R, Hudler P, Komel R. Single nucleotide polymorphisms rs911160 in AURKA and rs2289590 in AURKB mitotic checkpoint genes contribute to gastric cancer susceptibility. Environ Mol Mutagen. 2017;58:701–11.
Skeeles LE, Fleming JL, Mahler KL, Toland AE. The impact of 3’UTR variants on differential expression of candidate cancer susceptibility genes. PLoS One. 2013;8:e58609.
Bao G, Qu F, He L, Zhao H, Wang N, Ji G, et al. Prognostic significance of tag SNP rs1045411 in HMGB1 of the aggressive gastric cancer in a Chinese population. PLoS One. 2016;11:e0154378.
Wang X, Zhou YX, Qiao W, Tominaga Y, Ouchi M, Ouchi T, et al. Overexpression of aurora kinase a in mouse mammary epithelium induces genetic instability preceding mammary tumor formation. Oncogene. 2006;25:7148–58.
Hienonen T, Salovaara R, Mecklin JP, Jarvinen H, Karhu A, Aaltonen LA. Preferential amplification of AURKA 91A (Ile31) in familial colorectal cancers. Int J Cancer. 2006;118:505–8.
Gu J, Gong Y, Huang M, Lu C, Spitz MR, Wu X. Polymorphisms of STK15 (Aurora-a) gene and lung cancer risk in Caucasians. Carcinogenesis. 2007;28:350–5.
Ju H, Cho H, Kim YS, Kim WH, Ihm C, Noh SM, et al. Functional polymorphism 57Val>Ile of Aurora kinase a associated with increased risk of gastric cancer progression. Cancer Lett. 2006;242:273–9.
Jaboin JJ, Ausborn NL, Hwang M, Chen H, Niermann KJ, Giacalone NJ, et al. The Aurora kinase a polymorphisms are not associated with recurrence-free survival in prostate cancer patients. J Cancer Sci Ther. 2012;4:016–22.
Mesic A, Rogar M, Hudler P, Juvan R, Komel R. Association of the AURKA and AURKC gene polymorphisms with an increased risk of gastric cancer. IUBMB Life. 2016;68:634–44.
Marumoto T, Zhang D, Sava H. Aurora-A-a guardian of poles. Nat Rev Cancer. 2005;5:42–50.
Sehdev V, Katsha A, Arras J, Peng D, Soutto M, Ecsedy J, et al. HDM2 regulation by AURKA promotes cell survival in gastric cancer. Clin Cancer Res. 2014;20:76–86.
Tang W, Qiu H, Jiang H, Wang L, Sun B, Gu H. Aurora-a V57I (rs1047972) polymorphism and cancer susceptibility: a meta-analysis involving 27,269 subjects. PLoS One. 2014;9:e90328.
Vader G, Medema RH, Lens SM. The chromosomal passenger complex: guiding aurora-B through mitosis. J Cell Biol. 2006;173:833–7.
Gautschi O, Heighway J, Mack PC, Purnell PR, Lara PN Jr, Gandara DR. Aurora kinases as anticancer drug targets. Clin Cancer Res. 2008;14:1639–48.
González-Loyola A, Fernández-Miranda G, Trakala M, Partida D, Samejima K, Ogawa H, et al. Aurora B overexpression causes aneuploidy and p21Cip1 repression during tumor development. Mol Cell Biol. 2015;35:3566–78.
Tchatchou S, Wirtenberger M, Hemminki K, Sutter C, Meindl A, Wappenschmidt B, et al. Aurora kinases a and B and familial breast cancer risk. Cancer Lett. 2007;247:266–72.
Gordon S, Akopyan G, Garban H, Bonavida B. Transcription factor YY1: structure, function, and therapeutic implications in cancer biology. Oncogene. 2006;25:1125–42.
Broos S, Hulpiau P, Galle J, Hooghe B, Van Roy F, De Bleser P. ConTra v2: a tool to identify transcription factors binding sites across species, update 2011. Nucleic Acids Res. 2011;39:W74–8.
Kumar S, Ambrosini G, Bucher P. SNP2TFBS – a database of regulatory SNPs affecting predicted trabscription factor binding site affinity. Nucleic Acids Res. 2017;45:D139–44.
Fellmeth JE, Gordon D, Robins CE, Scott RT Jr, Treff NR, Schindler K. Expression and characterization of three Aurora kinase C splice variants found in human oocytes. Mol Hum Reprod. 2015;21:633–44.
Ulisse S, Delcros JG, Baldini E, Toller M, Curcio F, Giacomelli L, et al. Expression of Aurora kinases in human thyroid carcinoma cell lines and tissues. Int J Cancer. 2006;119:275–82.
Khan J, Ezan F, Cremet JY, Fautrel A, Gilot D, Lambert M, et al. Overexpression of active aurora-C kinase results in cell transformation and tumour formation. PLoS One. 2011;6:e26512.
Leppek K, Das R, Barna M. Functional 5’UTR mRNA structures in eukaryotic translation regulation and how to find them. Nat Rev Mol Cell Biol. 2018;19:158–74.
Calvo SE, Pagliarini DJ, Mootha VK. Upstream open reading frames cause widespread reduction of protein expression and are polymorphic among humans. Proc Natl Acad Sci U S A. 2009;106:7507–12.
Chatterjee S, Berwal SK, Pal JK. Pathological mutations in 5’ untranslated regions of human genes. In: Encyclopedia of Life Sciences (ELS). Chichester: Wiley; 2010. p. 1-8.
El-Omar EM, Rabkin CS, Gammon MD, Vaughan TL, Risch HA, Schoenberg JB, et al. Increased risk of noncardia gastric cancer associated with proinflammatory cytokine gene polymorphisms. Gastroenterology. 2003;124:1193–201.
De Oliveira JG, Rosi AFT, Nizato DM, Miyasaki K, Silva AE. Profiles of gene polymorphisms in cytokines and tool-like receptors with higher risk for gastric cancer. Dig Dis Sci. 2013;58:978–88.
Lens SM, Voest EE, Medema RH. Shared and separate functions of polo-like kinases and aurora kinases in cancer. Nat Rev Cancer. 2010;12:825–41.
Akdeli N, Riemann K, Westphal J, Hess J, Siffert W, Bachmann HS. A 3’UTR polymorphism modulates mRNA stability of the oncogene and drug target Polo-like kinase 1. Mol Cancer. 2014;13:87.
This work was supported by Slovenian Research Agency (ARRS) (P1–0390; J3–5504; BI-BA/14–15-010) and Federal Ministry of Education and Science of Bosnia and Herzegovina (FMON) (05–39–116-23/14). The funding body had no role in the design of the study and collection, analysis and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Ethical approval for this study and for the use of patient’s personal and medical data from the Clinical Pathology and Cytology at the University Clinical Center Sarajevo was issued by the Ethical Committee at the University Clinical Centre Sarajevo (No. 0302–36765). The data were analyzed anonymously. Written informed consents from all the patients for using their samples for the purposes of this study were obtained prior to surgery. The requirement to obtain the consents from the control individuals had been waived due to the following reasons: the control participants were healthy blood donors, randomly selected upon regular medical examinations whereby only the remains of their samples were used for the purposes of the present study, no additional risks to the subjects existed, no individually identifiable as well as medical or any other sensitive information were used for this study, the study would not adversely affect the rights and welfare of the subjects, and the results of the study would have no effect on the subjects.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.