Abstract

Wilms tumor is the most common type of renal malignancy in children. Previous studies have demonstrated that single nucleotide polymorphisms (SNPs) in the AURKA gene could predispose to several human malignancies. We recruited 145 cases and 531 cancer-free controls to investigate whether AURKA gene variants modify Wilms tumor susceptibility. Three AURKA SNPs (rs1047972 C>T, rs2273535 T>A, and rs8173 G>C) were genotyped by the Taqman methodology. Odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the strength of association between AURKA SNPs and Wilms tumor risk. We found that only the rs8173 G>C polymorphism was significantly associated with Wilms tumor risk (GC vs. GG: adjusted OR (AOR) = 0.50, 95% CI = 0.35–0.73, ; GC/CC vs. GG: AOR = 0.60, 95% CI = 0.42–0.88, ). Stratification analysis revealed that rs8173 GC/CC genotypes were associated with Wilms tumor risk among children aged >18 months (AOR = 0.56, 95% CI = 0.34–0.93, ), male children (AOR = 0.54, 95% CI = 0.33–0.90, ), and children with clinical stage III + IV diseases (AOR = 0.56, 95% CI = 0.35–0.90, ). Haplotype analysis indicated that the CAG haplotype was significantly associated with increased Wilms tumor risk. In conclusion, our findings indicated that the AURKA rs8173 G>C polymorphism was associated with decreased Wilms tumor risk in Chinese children.

1. Introduction

Wilms tumor (WT), also known as nephroblastoma, is the most common renal malignancy in children [1, 2]. It accounts for 6% to 7% of malignant tumors in children under the age of 15 years, with an incidence rate of about 7–10 cases per million in Western countries [3]. The incidence rate of WT in Chinese is around 3.3 cases per million [4]. Over the past 20 years, because of the application of multimodality treatment, including surgery, radiotherapy, chemotherapy, and autologous stem cell transplantation, the overall survival rate of 5 years has increased from 30% to 90% [2, 57]. Despite the great achievements in the treatment of WT, the prognosis of nearly 25% of patients with high-risk diseases remains unsatisfying [8]. In addition, about 25% of survivors suffer from the high cost of treatment and the physical torment of some chronic conditions [9, 10].

Previous studies have found that genetic factors contribute to the risk of Wilms tumor and identified a number of genes associated with WT, including Wilms’ tumor protein 1 (WT1), β-catenin, tumor protein 53 (TP53), catenin beta 1 (CTNNB1), and AMER1 [1113]. Although recent genome-wide studies on WT have revealed some previously unknown gene mutations implicated with WT, current known genetic variants are not adequate to fully elucidate the pathogenesis of WT [14]. Therefore, it is particularly important to identify more causal genetic variants for WT.

Aurora kinase A (AURKA) is located at chromosome 20q13.2, consisting of 12 exons [15]. It encodes a centrosome-related serine/threonine kinase found to be overexpressed in many human cancers, including primary colorectal carcinoma, esophageal squamous cell carcinoma, neuroblastoma, and breast and ovarian cancers [1619]. Overexpression of AURKA leads to abnormal centrosome amplification, thereby affecting the stability of the genome and inducing tumorigenesis [16, 20]. Furthermore, inhibition of AURKA expression could result in abnormal cell mitosis, eventually leading to cell death [2123]. Therefore, AURKA might be an attractive target for anticancer therapies [24].

In addition, previous studies had found that the AURKA gene single nucleotide polymorphism (SNP) was closely related to the risk of several human malignancies [2527]. Early research demonstrated that the AA genotype of AURKA rs2273535 T>A could increase the risk of oral cancer [28]. Studies in recent years have also found that AURKA rs2273535 T>A polymorphism could enhance the susceptibility of gastric cancer [15]. Furthermore, there were related studies reported that Caucasians carrying AURKA rs1047972 T>C had a relatively low risk of developing breast cancer [26]. However, studies regarding their association with WT risk are lacking. Therefore, we performed this case-control study with 145 cases of WT and 531 control subjects to assess the association between AURKA gene polymorphisms and WT susceptibility in Chinese children.

2. Materials and Methods

2.1. Study Subjects

In this study, a total of 145 cases of WT and 531 controls were collected from the Guangzhou Women and Children Medical Center from 2001 to 2016, as we described previously [29]. Controls were selected from those receiving routine physical examination during the same period. The cases included in this study were all sporadic and the control groups were confirmed free of history of malignancy. Both WT patients and controls were frequency-matched based on the age, gender, and ethnicity. The frequency distribution information of selected variables for WT cases and cancer-free controls was summarized in Supplementary Table 1. We obtained written informed consents from the parents or legal guardians of all subjects. The study protocol was approved by the Institutional Review Board of the Guangzhou Women and Children’s Center (Ethical Approval Number: 2018022102).

2.2. SNP Selection and Genotyping

The selection criteria for potentially functional polymorphisms have been described in detail previously [30]. Briefly, we choose potential functional polymorphisms in the 5′-flanking region, 5′UTR, 3′UTR, and exons of the AURKA gene. Three SNPs in the AURKA gene were finally selected (rs1047972 C>T (Assay ID: AN322ZF), rs2273535 T>A (Assay ID: C__25623289_10), and rs8173 G>C (Assay ID: C__8947675_10)). The rs1047972 C>T polymorphism may lead to amino acid alteration from Val to Ile at codon 57; the rs2273535 T>A polymorphism may lead to splicing alteration, whereas the rs8173 G>C polymorphism may lead to miRNA binding alteration. These three SNPs could capture additional fourteen SNPs (Supplementary Table 2). Moreover, there was no significant linkage disequilibrium between paired polymorphisms in Supplementary Figure 1 (rs8173/rs1047972, R2 = 0.119; rs8173/rs2273535, R2 = 0.527; rs1047972/rs2273535, R2 = 0.291) [31]. Human genomic DNA was mainly extracted from the peripheral blood samples using the TIANamp Blood DNA Kit (TianGen Biotech, Beijing, China). The mean yield of DNA was 12.21 μg (median = 6.25 μg; range = 1.08 to 60 μg), with a mean A260/A280 ratio of 1.89 (median = 1.90; range = 1.78 to 2.05). AURKA SNPs were genotyped using the ABI-7900 Sequence Detection System (Applied Biosystem, Foster City, CA) and Taqman real-time PCR, as described previously [30]. The information on all the genotyped samples for all three SNPs was summarized in Supplementary Figure 2.

2.3. Genotype and Gene Expression Correlation Analysis

The FAM210B mRNA expression data were retrieved for transformed fibroblasts cells with different genotypes of AURKA SNPs from the GTEx Portal database (https://www.gtexportal.org/home/). These data were used to assess the correlation between rs8173 G>C genotypes and mRNA expression-level alteration [32].

2.4. Statistical Analysis

All statistical analysis of data was processed by SAS (version 9.4; SAS Institute, Cary, NC, USA). The distribution of subject characteristics between cases and controls was examined by a bilateral χ2 test. The goodness-of-chi-squared test was used to validate whether the genotype frequencies in the controls were consistent with the Hardy–Weinberg equilibrium (HWE). Multivariate logistic regression analysis was adopted to calculate the odds ratios (ORs) and 95% confidence intervals (CIs). The adjusted ORs and 95% CIs were calculated using an unconditional logistic regression model with adjustment for age and gender to assess the association between AURKA polymorphisms and WT risk. Inferred haplotypes of AURKA gene were based on observed genotypes, and adjusted ORs and 95% CIs were obtained by logistic regression models with adjustment for age and gender. was considered to be statistically significant.

3. Results

3.1. Association between AURKA Gene Polymorphisms and Wilms Tumor Susceptibility

In this study, 143 cases and 531 controls were successfully genotyped. All selected SNP genotypes had a frequency distribution consistent with HWE (rs1047972 C>T, ; rs2273535 T>A, ; rs8173 G>C, ). The genotype frequencies of the case group and control group are shown in Table 1. We found rs8173 G>C polymorphism was significantly associated with decreased WT risk (GC vs. GG: adjusted OR (AOR) = 0.50, 95% CI = 0.35–0.73, ; GC/CC vs. GG: AOR = 0.60, 95% CI = 0.42–0.88, ).

3.2. Stratification Analysis

Subsequently, we explored the association between rs8173 G>C polymorphism and WT risk by stratified analysis. The results of the stratified analysis were based on the age, gender, and clinical stage (Table 2). We found that carriers of rs8173 GC/CC genotypes had a decreased WT risk when compared with GG genotype carriers, among children older than 18 months (AOR = 0.56, 95% CI = 0.34–0.93, ), male children (AOR = 0.54, 95% CI = 0.33–0.90, ), and those with the clinical stage III + IV diseases (AOR = 0.56, 95% CI = 0.35–0.90, ).

3.3. Haplotype Analysis

Based on a combined analysis of the three SNP polymorphisms in the AURKA gene, eight haplotypes were inferred (Table 3). Compared with the reference CAC haplotype, only CAG haplotype was associated with significant increased WT risk (AOR = 1.99, 95% CI = 1.05–3.77, ).

3.4. Genotype-Based mRNA Expression Analysis

We found that the rs8173 G>C polymorphism was significantly associated with altered AURKA gene expression in transformed fibroblast cells () using data from GTEx portal (Figure 1).

4. Discussion

In this case-control study, we analyzed 143 WT patients and 531 controls to investigate the association between three AURKA SNPs (rs1047972, rs2273535, and rs8173) and WT risk of in Chinese children. To the best of our knowledge, this is the first study to investigate the association between AURKA gene polymorphisms and WT susceptibility.

Previous studies have found that AURKA was overexpressed in several common human malignancies, which in turn promoted cell proliferation and tumor progression and metastasis [3335]. The abnormal expression of AURKA could lead to abnormal chromosome segregation, decreased chromosome stability, and finally increase the susceptibility to malignant transformation [36]. Studies in recent years have shown that AURKA SNPs were closely related to cancer risk. An early study by Lee et al. indicated that the AURKA gene rs2273535 was associated with oral cancer risk [28]. A research by Guo et al. [37] illustrated that the rs2273535 polymorphism was closely related to the increased risk of breast cancer, especially in Asian populations. Interestingly, studies have also found that the Phe allelic variation in AURKA rs2273535 appears to prevent breast cancer in Malaysian Chinese [38]. In addition, Dai et al. [26] reported that the AURKA rs1047972 polymorphism was associated with reduced incidence of breast cancer in Caucasians. These findings suggest that the effects of AURKA SNPs might be tissue dependent and ethnicity dependent.

We previously did not find a significant association between AURKA SNPs and neuroblastoma susceptibility [31]. Interestingly, in this study, our results showed that AURKA rs8173 G>C significantly reduced the WT risk, although no significant association was found for AURKA rs1047972 C>T and rs2273535 T>A. In addition, stratified analysis revealed that individuals carrying AURKA rs8173 GC/CC genotypes had significantly decreased susceptibility to WT in several subgroups, including older than 18 months, male, and clinical stages III + IV.

A previous study showed that the AURKA rs2273535 polymorphism had a strong LD with the rs1047972 genotype. Patients with AURKA haplotype variants exhibited high kinase activity and tended to develop advanced gastric cancer more readily [39]. In this study, we found that the WT risk in individuals with rs1047972/rs2273535/rs8173 CAG haplotype almost doubled, when compared with individuals with CAC haplotype. However, there were no significant association between the other six haplotypes and the risk of WT.

AURKA has been extensively investigated in neuroblastoma. ShRNA-mediated AURKA gene silencing assays have demonstrated that reducing AURKA expression would inhibit cell proliferation in neuroblastoma [4042]. It is worth mentioning that studies also have found that overexpression or expansion of AURKA is associated with poor prognosis in a variety of cancer patients and inhibition of AURKA expression can trigger tumor cell death [2123]. Previous studies have shown that the LIN28B-RAN-AURKA axis was involved in the development of neuroblastoma, and AURKA, as a confluence of LIN28B-RAN signaling, could further promote cell cycle progression by phosphorylating many cell cycle regulators and stabilizing N-myc protein (encoded by MYCN gene) [41, 42]. The LIN28B-RAN-AURKA-MYCN signaling cascade in the development of neuroblastoma provides a new insight into the molecular mechanism by which AURKA rs8173 reduced the risk of WT. The role of AURKA in the development of Wilms tumor is currently lacking. However, based on literature search, we found that AURKA is involved in a variety of tumors. Because of the universal contributions of AURKA in tumors, we hypothesized that AURKA might be also implicated in WT and therefore performed the current study. Based on previous reports, we speculated that the AURKA rs8173 might reduce the risk of WT by inhibiting the proliferation and migration of tumor cells through the LIN28B-RAN-AURKA-MYCN signaling pathway.

Some shortcomings should be mentioned in our research. First, the study included 145 patients and 531 controls. Because of the limited sample size, some important results may be accidental. Studies with large samples are indispensible to verify our results. Second, only three AURKA SNPs were investigated in the study. Other potentially functional AURKA gene polymorphisms should be studied in the future. Third, the genotype distribution in this hospital-based study might not represent genotype distribution across the general population, which may bias the case-control study to some extent. Finally, functional experiments should be conducted to strengthen the findings in the current study.

5. Conclusions

Overall, our study confirmed that the AURKA rs8173 G>C polymorphism is associated with reduced WT risk in Chinese children. In the future, multicenter collaboration is needed to further expand the sample size from different regions and different races to clarify the impact of AURKA SNPs on the risk of WT more accurately.

Abbreviations

WT:Wilms tumor
WT1:Wilms’ tumor protein 1
TP53:Tumor protein 53
CTNNB1:Catenin beta 1
AURKA:Aurora kinase A
SNP:Single nucleotide polymorphism
HWE:Hardy–Weinberg equilibrium
OR:Odds ratio
CI:Confidence interval.

Data Availability

All the data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Tongyi Lu, Li Li, and Jinhong Zhu contributed equally to this work. All authors contributed significantly to this work. TL, JL, AL, WF, and GL performed the research study and collected the samples and data. JH analyzed the data. HX, JH, and TZ designed the research study. TL, LL, and JZ wrote the paper. JH prepared all the tables and figures. In addition, all authors read, reviewed, and approved the final manuscript to be published.

Acknowledgments

This study was funded by grants from the Pearl River S&T Nova Program of Guangzhou (no. 201710010086), the National Natural Science Foundation of China (nos. 81803320 and 81560262), the Science and Technology Project of Guangzhou (no. 201804010037), the Guangdong Provincial Key Laboratory of Research in Structural Birth Defect Disease (no. 2019B030301004), and the Basic Applied Study Planning Projects of Yunnan Province (no. 2018FB130).

Supplementary Materials

Supplemental Table 1: frequency distribution of selected variables for Wilms tumor cases and cancer-free controls. Supplemental Table 2: polymorphisms captured by the three included AURKA polymorphisms as predicted by SNPinfo online software. Supplemental Figure 1: linkage disequilibrium (LD) analysis for the three included polymorphisms in Han Chinese population consisting of CHB (Han Chinese in Beijing, China) and CHS (Southern Han Chinese) subjects. Supplemental Figure 2: diagram showing all the genotyped samples for the three included polymorphisms. (Supplementary Materials)