Next Article in Journal
Dysfunction of Endocytic Kinase AAK1 in ALS
Next Article in Special Issue
Associations between Variants in IL-33/ST2 Signaling Pathway Genes and Coronary Heart Disease Risk
Previous Article in Journal
LRP-1: Functions, Signaling and Implications in Kidney and Other Diseases
Previous Article in Special Issue
Risk-Association of DNMT1 Gene Polymorphisms with Coronary Artery Disease in Chinese Han Population
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Associations of NR5A2 Gene Polymorphisms with the Clinicopathological Characteristics and Survival of Gastric Cancer

1
Department of Oncology, Nanjing First Hospital, Nanjing Medical University, 68 Changle Road, Nanjing 210006, China
2
Department of Oncology, Nantong Tumor Hospital, Nantong 226300, China
3
Department of Environmental Genomics, Jiangsu Key Laboratory of Cancer Biomarkers, Prevention and Treatment, Cancer Center, Nanjing Medical University, Nanjing 211166, China
4
Department of Genetic Toxicology, the Key Laboratory of Modern Toxicology of Ministry of Education, School of Public Health, Nanjing Medical University, Nanjing 211166, China
5
Department of Surgery, Yixing People's Hospital, Yixing 214200, China
6
Department of General Surgery, Yixing Tumor Hospital, Yixing 214200, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2014, 15(12), 22902-22917; https://doi.org/10.3390/ijms151222902
Submission received: 22 September 2014 / Revised: 15 October 2014 / Accepted: 16 November 2014 / Published: 10 December 2014
(This article belongs to the Special Issue Human Single Nucleotide Polymorphisms and Disease Diagnostics)

Abstract

:
The orphan nuclear receptor (NR5A2), which belongs to the NR5A subfamily of nuclear receptors, is expressed in developing and adult tissues of endodermal origin, and can contribute to the development of several cancers through regulating cell proliferation. NR5A2 (rs3790843 and rs3790844) single nucleotide polymorphisms (SNPs) genotyping were examined in DNA samples, extracted from paraffin-embedded cancer tissue. Clinicopathologic and follow-up data were collected from 944 patients with gastric cancer (GC). Associations of the 2 SNPs with the progression and prognosis in gastric cancer patients were analyzed using the SPSS version 18.0. We found that NR5A2 rs3790843 polymorphism was significantly associated with the risk of GC which had regional lymph node metastasis (p = 0.044) or distant metastasis (p = 0.020). Our results also indicated that rs3790844 polymorphism was associated with the increased overall survival (OS) of GC patients in the dominant model (GG vs. GA/AA, HR (hazard ratio) = 0.823, 95% CI (confidence interval) = 0.679–0.997), suggesting a potential protective role of the variant A allele. Additionally, in the stratified analysis, both NR5A2 rs3790843 and rs3790844 polymorphism were associated with significantly lower risk of death in the groups of female, tumor size >5 cm in a dominant model. Our results represent the first demonstration that the NR5A2 rs3790844 polymorphism is associated with increased OS of GC patients in the dominant model, and similar results were found among the female group and tumor size >5 cm group for NR5A2 rs3790843 polymorphism. Further validation in other larger studies with different ethnic populations and functional evaluations are needed.

Graphical Abstract

1. Introduction

Gastric cancer (GC), as the third leading cause of cancer death in men and the fifth in women, remains one of the major public health problems worldwide [1]. About two-thirds of the cases occur in developing countries and 42% in China alone, while the estimated high-risk areas in developed countries include Eastern Europe, and parts of Central and South America [2]. In recent decades, its incidence rate has declined, and remarkable progress has been achieved in comprehensive treatment strategies of combined therapy. However, the prognosis for patients still remains poor, with 5-year overall survival rates of 30% [3]. The development of gastric cancer is a multi-step, sequential process, which initiates from chronic gastritis, atrophy, intestinal metaplasia, dysplasia, and finally malignant transformation to invasive gastric cancer [4]. Moreover, multiple genetic and epigenetic alterations are implicated in the multi-step process of human stomach carcinogenesis and development [5]. Therefore, discovery of these biomarkers could be helpful in the improvement of early diagnosis, screening of high-risk individuals, as well as patient care [6,7]. In recent years, increased studies have focused on the detection of genetic variants that could play roles in the development, progression and prognoses of gastric cancer [8].
The orphan nuclear receptor (NR5A2), also known as liver receptor homologue-1 (LRH-1) and fetoprotein transcription factor (FTF), belongs to the NR5A or FTZ-F1 subfamily of nuclear receptors. NR5A2 is expressed in developing and adult tissues of endodermal origin, including liver, pancreas, intestine and the ovary [9,10]. Functionally, NR5A2 has been implicated in the regulation of bile acid and cholesterol homeostasis [9,10] and the regulation of inflammatory responses in the liver and gut [11]. Aromatase is an NR5A2 target gene, which catalyses the conversion of androgens (testosterone and primarily androstenedione) to oestrogens. NR5A2 may aid breast cancer progression in postmenopausal women by promoting local oestrogen biosynthesis [12,13,14]. The expression of NR5A2 is also elevated in pancreatic cancer and promotes pancreatic cancer cell growth through stimulation of cyclin D1, cyclin E1 and c-Myc [15]. A recent genome-wide association study (GWAS) identified rs3790844 and rs3790843, located in the first intron of NR5A2 in 1q32.11, as associated with pancreatic cancer susceptibility [16]. In the gastrointestinal tract, NR5A2 has been shown to participate in intestinal cell renewal [17] and is expressed in the stomach epithelium [18]. Therefore, we conducted this study to examine whether NR5A2 rs3790844 and rs3790843 polymorphisms are associated with clinical outcomes of gastric cancer. The two SNPs may serve as potential molecular prognostic markers for gastric cancer, which will promote further defined sub-populations at higher risk of the disease. Consequently, these sub-populations may require more rigorous treatment and postoperative follow up.

2. Results

2.1. Associations between Clinicopathological Variables and Overall Survival

The patients’ characteristics and clinical information are summarized in Table 1. In the follow-up period of 119 months, 442 patients died. There were 727 males (77.0%) and 217 females (23.0%), with a median age of 62 years ranging from 28 to 83 years. Clinicopathological characteristics including tumor size, depth of invasion, lymph node metastasis, distant metastasis, Tumor, Node and Metastasis (TNM) stage and Lauren classification were significantly associated with survival time (log-rank p < 0.05). Specifically, patients with tumor size >5 cm (median survival time (MST), 48 months) had a 42% significantly higher risk of death (HR (hazard ratio) = 1.42, 95% CI (confidence interval) = 1.178–1.716), compared with those with tumor size <5 cm (MST, 98 months), and the patients with lymph node metastasis or distant metastasis had significantly higher risk of death than those patients without lymph node metastasis or distant metastasis (log-rank p < 0.05). In addition, as the depth of invasion and TNM stage increased, the risk of death for gastric cancer showed a significant increase in a dose-response manner (log-rank p < 0.05).
Table 1. Associations between clinicopathological variables and overall survival.
Table 1. Associations between clinicopathological variables and overall survival.
VariablePatients, n = 944MST (Months)Log-Rank pHR (95% CI) d
Age (years)
≤60442880.2391.000
>6050259 1.118 (0.927–1.349)
Sex
Male727700.5361.000
Female21763 1.071 (0.860–1.335)
Tumor Size
≤5 cm58198<0.0011.000
>5 cm36348 1.422 (1.178–1.716)
Location
Non-Cardia626780.2531.000
Cardia31863 0.891 (0.730–1.088)
Histological Types a
Intestinal370740.4361.000
Diffuse49760 1.069 (0.918–1.244)
Differentiation a
Well to Moderate305800.5181.000
Poorly49359 1.165 (0.942–1.443)
Mucinous/Signet-Ring Cell6962 1.199 (0.825–1.742)
Lauren a
139976<0.0011.000
254150 1.467 (1.211–1.776)
Depth of Invasion b
T118284<0.0011.000
T213878 0.540 (0.408–0.714)
T3870 0.800 (0.608–1.054)
T459451 1.008 (0.416–2.442)
Lymph Node Metastasis c
N037881<0.0011.000
N1/N2/N356643 1.814 (1.480–2.222)
Distant Metastasis
M0886740.0041.000
M15826 1.646 (1.165–2.326)
TNM Stage
I25083<0.0011.000
II20388 1.231 (0.909–1.666)
III45841 1.949 (1.524–2.492)
IV2544 2.046 (1.163–3.600)
Chemotherapy
No638610.6841.000
Yes30669 1.043 (0.852–1.276)
Mean survival time was presented when the median survival time could not be measured. TNM, Tumor, Node and Metastasis; MST, median survival time; HR, hazard ratio; CI, confidence interval; AJCC, American Joint Commission on Cancer. a: Partial data were not available, and statistics were based on available data; b: Invaded depth of tumor was classified according to the criteria of AJCC 7th; c: Lymph nodes were staged according to tumor node metastasis classification of the 7th edition of AJCC in which the number of lymph nodes with a metastasis of 1–2, 3–6, and ≥7 were classified as N1, N2, and N3, respectively; d: Adjusted by age and sex.

2.2. Association Analyses of NR5A2 rs3790843 and rs3790844 Genotypes with Clinicopathological Features

Among 944 specimens from patients with gastric cancer, NR5A2 rs3790843 was successfully genotyped in 907 specimens, and NR5A2 rs3790844 was successfully genotyped in 912 specimens. The frequency of each NR5A2 rs3790843 genotypes was 45.6% (414 specimens) for the TT variant, 45.2% (410 specimens) for the CT variant, and 9.2% (83 specimens) for the CC variant. Table 2 indicates that the risk of GC for the CT and CC variants compared with that for the TT variant in NR5A2 rs3790843 was associated significantly with regional lymph node metastasis (p = 0.044) and distant metastasis (p = 0.020) but not with age, sex, tumor size or location, histological type, differentiation, depth of invasion, chemotherapy history, American Joint Committee on Cancer staging, or Lauren classification. In reference to NR5A2 rs3790844, the GG, GA, and AA genotypes were identified in 399 specimens (43.8%), 416 specimens (45.6%), and 97 specimens (10.6%), respectively. Compared with the GG genotype, the GA/AA genotypes were not associated with any clinicopathologic data.
Table 2. Association analyses of NR5A2 rs3790843 and rs3790844 genotypes with clinicopathological features.
Table 2. Association analyses of NR5A2 rs3790843 and rs3790844 genotypes with clinicopathological features.
Variable rs3790843 (n = 907)rs3790844 (n = 912)
TTTC/CCpGGGA/AAp
Age (years)≤601922370.6101822490.380
>60222256217264
SexMale3213750.6013113910.539
Female9311888122
LocationNon-Cardia Cancer2793200.4322643380.930
Cardia Cancer135173135175
Tumor Size≤5 cm2613020.5812433240.486
>5 cm153191156189
Lymph Node Metastasis aN01801820.0441681970.257
N1/N2/N3234311231316
Distant MetastasisM03984560.020 3824780.098
M11637935
Histological Types bIntestinal1681870.2731562050.541
Diffuse207270206270
Differentiation bWell to Moderate1371550.5561271700.687
Poorly208267205270
Mucinous/Signet-Ring Cell30353035
Lauren b11821990.2531692180.450
2229293227294
ChemotherapyNo2753330.7212663480.709
Yes139160133165
Depth of Invasion cT186930.854761010.929
T260695971
T33324
T4257320254326
TNM StageI1201180.3691111280.533
II8311178119
III198248198250
IV1013913
HR, hazard ratio; CI, confidence interval; AJCC, American Joint Commission on Cancer; TNM, Tumor, Node and Metastasis. a: Lymph nodes were staged according to tumor node metastasis classification of the 7th edition of AJCC in which the number of lymph nodes with a metastasis of 1–2, 3–6, and ≥7 were classified as N1, N2, and N3, respectively; b: Partial data were not available, and statistics were based on available data; c: Invaded depth of tumor was classified according to the criteria of AJCC 7th.
The distribution of genotypes in the population was consistent with Hardy–Weinberg equilibrium (p = 0.20 for rs3790843 SNP, p = 0.46 for rs3790844 SNP), while strong linkage disequilibrium (LD) exists between alleles of the two loci with the following R2 values: rs3790843/rs3790844 = 0.798. Haplotype analyses were performed for all patients using the SHEsis software (http://analysis. bio-x.cn), and the six possible haplotype frequencies of the two SNPs are shown in Table 3. According to the results, the TTGG and CTGA haplotypes are the most common and represent 41.5% and 40.8%, respectively, in all groups. After haplotype analyses, however, no significant differences in these haplotype frequencies were found in any group.
Table 3. Haplotype analyses of NR5A2 rs3790843 and rs3790844 genotypes.
Table 3. Haplotype analyses of NR5A2 rs3790843 and rs3790844 genotypes.
HaplotypeFrequenciesPatients/DeathsHR (95% CI) a
TTGG0.415379/1841.00
CTGA0.408373/1700.90 (0.73–1.10)
CCAA0.07972/300.79 (0.54–1.17)
TTGA0.03633//110.54 (0.30–1.00)
CTAA0.02422//60.47 (0.21–1.07)
CTGG0.01514//101.55 (0.81–2.94)
a: Adjusted by age and sex.

2.3. Associations of NR5A2 rs3790843 and rs3790844 with Clinical Outcomes of OS

The Kaplan-Meier survival analysis and Cox proportional hazard models were used to assess the prognostic effect of NR5A2 rs3790843 and rs3790844 on GC patients in different genetic models (Table 4). We found rs3790844 polymorphism was associated with the OS of GC patients in the dominant model (GG vs. GA/AA, log-rank p = 0.045, Figure 1). Patients with the GA or AA genotype (MST 88 months) had a 17.7% lower risk of death than that of patients with the GG genotype (MST 52 months) (HR = 0.823, 95% CI = 0.679–0.997), suggesting a potential protective role of the variant A allele. No significant associations were observed between the rs3790843 genotypes and OS of GC patients in any genetic models. We further evaluated the associations by stratified analysis of variants of clinicopathological features (Table 5). In the dominant model, NR5A2 rs3790844 polymorphism was associated with significantly lower risk of death in the groups of patients with lymph node metastasis, no distant metastasis, diffuse type and no chemotherapy history (log-rank p < 0.05). Additionally, both NR5A2 rs3790843 and rs3790844 polymorphism were associated with significantly better prognosis among the female patients group and tumor size >5 cm group in a dominant model. Cox stepwise regression analysis was conducted to evaluate the independent effect of clinicopathological variables, rs3790843 and rs3790844 SNP on the OS of the patients with GC. As shown in Table 6, two variables (lymph node metastasis and NR5A2 rs3790844) were included in the regression model by stepwise selection of the covariant variables and rs3790844 SNP was shown to be an independent protective factor for GC with a 21.6% decreased risk (HR = 0.784, 95% CI = 0.646–0.951, p = 0.014).
Table 4. Associations of NR5A2 rs3790843 and rs3790844 with clinical outcomes of overall survival.
Table 4. Associations of NR5A2 rs3790843 and rs3790844 with clinical outcomes of overall survival.
Genetic ModelGenotypesPatientsDeathsMST (Months)Log-Rank pHR (95%CI) a
rs3790843
Codominant ModelTT414196590.6681.000
CT410187700.934 (0.765–1.142)
CC8336880.872 (0.611–1.244)
Dominant ModelTT414196590.4151.000
CT or CC493223700.924 (0.762–1.119)
Recessive ModelTT or CT824383700.5511.000
CC8336880.902 (0.641–1.269)
rs3790844
Codominant ModelGG399197520.0751.000
GA416186780.853 (0.698–1.043)
AA9737980.699 (0.492–0.993)
Dominant ModelGG399197520.0451.000
GA or AA513223880.823 (0.679–0.997)
Recessive ModelGG or GA815383650.1031.000
AA9737980.757 (0.540–1.061)
MST, median survival time; HR, hazard ratio; CI, confidence interval. Mean survival time was presented when the median survival time could not be measured. a: Adjusted by age and sex.

3. Discussion

In the present study, we investigated the effects of two SNPs (rs3790843 and rs3790844) of the NR5A2 gene on the progression and survival of GC in Chinese populations. We found that NR5A2 rs3790843 polymorphism was significantly associated with the risk of GC when compared with regional lymph node metastasis and distant metastasis. Our results also indicated that rs3790844 polymorphism was associated with the increased OS of GC patients in the dominant model, suggesting a potential protective role of the variant A allele. No significant associations were observed between the rs3790843 genotypes and OS of GC patients in any genetic models. Additionally, in the stratified analysis, NR5A2 rs3790844 polymorphism was associated with significantly lower risk of death in the groups of female, tumor size >5 cm, lymph node metastasis, no distant metastasis, diffuse type and no chemotherapy history in the dominant model. Similar results were found among the female patients group and tumor size >5 cm group for the NR5A2 rs3790843 polymorphism in a dominant model.
Figure 1. Overall survival curve in relation to NR5A2 rs3790844 polymorphism in patients with gastric cancer in dominant model.
Figure 1. Overall survival curve in relation to NR5A2 rs3790844 polymorphism in patients with gastric cancer in dominant model.
Ijms 15 22902 g001
Nuclear receptor 5 subtype A2 (NR5A2), also known as Liver receptor homolog-1 (LRH-1), is a member of the orphan family of nuclear receptors that participates in a wide range of developmental processes [9,10]. Functionally, NR5A2 has been implicated in the regulation of bile acid and cholesterol homeostasis, and is expressed in developing and adult tissues of endodermal origin, including liver, pancreas, intestine and the ovary [9,10]. NR5A2 was first isolated as a transcriptional activator of the alpha-fetoprotein (AFP) gene, which encodes a factor that plays a crucial rolein hepatic specification [19]. Further, NR5A2 expression is elevated in pancreatic cancer andpromotes pancreatic cancer cell growth through stimulation of cyclin D1, cyclin E1 and c-Myc [15], while genome-wide association studies implicate mutations in the NR5A2 gene in pancreatic ductal adenocarcinoma [16]. In the gastrointestinal tract, NR5A2 has been shown to participate in intestinal cell renewal [17] and is expressed in the stomach epithelium [18]. However, it seemed no study focused on the association between NR5A2 and the progression and prognoses of gastric cancer. Thus, this study was conducted, and it indicates that the NR5A2 polymorphism is significantly associated with the development and overall survival of gastric cancer.
In the gastrointestinal tract, NR5A2 regulates proliferation through at least two different mechanisms: direct activation of the cyclin E1 promoter and, in cooperation with β-catenin, induction of cyclin D1 and c-Myc transcription [20]. A previous study showed that NR5A2, which interacts with β-catenin, induces cell proliferation through the concomitant induction of cyclin D1 and E1, which are downstream NR5A2 effectors whose induction may account for the diminished threshold for G/S transition [20]. Whereas β-catenin coactivates NR5A2 on the cyclin E1 promoter, NR5A2 acts as a potent tissue-restricted coactivator of β-catenin on the cyclin D1 promoter.
Table 5. Stratified analysis of rs3790843 and rs3790844 polymorphism among GC patients.
Table 5. Stratified analysis of rs3790843 and rs3790844 polymorphism among GC patients.
rs3790843rs3790844
TTTC/CCHR (95% CI) epGGGA/AAHR (95% CI) ep
Age (years)
≤6088/192108/2370.977 (0.737–1.294)0.86990/182107/2490.806 (0.609–1.067)0.132
>60108/222115/2560.873 (0.672–1.136)0.313107/217116/2640.837 (0.644–1.089)0.185
Sex
Male144/321175/3751.044 (0.838–1.302)0.700149/311172/3910.877 (0.705–1.093)0.243
Female52/9348/1180.619 (0.417–0.916)0.01748/8851/1220.660 (0.445–0.981)0.040
Location
Non–CardiaCancer133/279148/3200.948 (0.750–1.198)0.652133/264148/3380.810 (0.641–1.024)0.077
CardiaCancer63/13575/1730.861 (0.615–1.204)0.38164/13575/1750.840 (0.601–1.174)0.308
Tumor Size
≤5 cm108/261127/3021.044 (0.807–1.349)0.745106/243132/3240.927 (0.718–1.197)0.559
>5 cm88/15396/1910.747 (0.558–0.999)0.04991/15691/1890.693 (0.517–0.928)0.014
Lymph Node Metastasis a
N064/18060/1820.947 (0.666–1.347)0.76363/16863/1970.865 (0.610–1.226)0.414
N1/N2/N3132/234163/3110.829 (0.659–1.043)0.109134/231160/3160.744 (0.591–0.937)0.012
Distant Metastasis
M0188/398200/4560.886 (0.726–1.081)0.232188/382202/4780.793 (0.650–0.968)0.022
M18/1623/371.340 (0.599–2.998)0.4767/921/351.189 (0.544–2.599)0.664
Histological Types b
Intestinal71/16880/1871.041 (0.756–1.434)0.80569/15684/2050.952 (0.692–1.309)0.762
Diffuse108/207127/2700.826 (0.639–1.068)0.145111/206123/2700.735 (0.568–0.950)0.019
Differentiation b
Well to Moderate57/13767/1551.099 (0.772–1.565)0.60156/12769/1700.962 (0.676–1.368)0.828
Poorly104/208126/2670.882 (0.680–1.144)0.344105/205125/2700.822 (0.634–1.067)0.140
Others c18/3014/350.522 (0.259–1.056)0.0719/3013/350.397 (0.195–0.808)0.011
Lauren b
169/18278/1991.049 (0.759–1.451)0.77169/16980/2180.899 (0.652–1.241)0.655
2125/229144/2930.824 (0.648–1.047)0.113126/227142/2940.763 (0.600–0.970)0.027
Chemotherapy
No132/275151/3330.908 (0.719–1.147)0.416136/266149/3480.767 (0.608–0.968)0.025
Yes64/13972/1600.957 (0.683–1.340)0.79961/13374/1650.953 (0.679–1.338)0.783
Depth of Invasion d
T125/8634/931.297 (0.774–2.173)0.32424/7633/1011.051 (0.621–1.779)0.852
T231/6024/690.669 (0.393–1.141)0.140 31/5925/710.648 (0.382–1.099)0.107
T31/32/31.405 (0.125–15.838)0.7821/22/40.809 (0.071–9.157)0.864
T4134/257159/3200.880 (0.699–1.108)0.276135/254157/3260.809 (0.642–1.019)0.072
TNM Stage
I40/12039/1180.992 (0.638–1.542)0.97240/11140/1280.862 (0.556–1.336)0.506
II38/8338/1110.743 (0.474–1.166)0.19637/7841/1190.711 (0.456–1.109)0.132
III111/198135/2480.861 (0.670–1.108)0.245113/198132/2500.800 (0.622–1.029)0.082
IV4/108/131.855 (0.557–6.171)0.3144/97/131.257 (0.368–4.302)0.715
HR, hazard ratio; CI, confidence interval; AJCC, American Joint Commission on Cancer; TNM, Tumor, Node and Metastasis. a: Lymph nodes were staged according to tumor node metastasis classification of the 7th edition of AJCC in which the number of lymph nodes with a metastasis of 1–2, 3–6, and ≥7 were classified as N1, N2, and N3, respectively; b: Partial data were not available, and statistics were based on available data; c: Others: mucinous carcinoma and Signet-ring cell carcinoma; d: Invaded depth of tumor was classified according to the criteria of AJCC 7th; e: Adjusted by age and sex.
Table 6. Stepwise Cox regression analysis on the survival of GC.
Table 6. Stepwise Cox regression analysis on the survival of GC.
VariablesβSEHR b95% CIp Value
Age a0.0670.1001.070(0.878–1.302)0.503
Sex0.040.1181.041(0.825–1.313)0.743
Lymph Node Metastasis0.5670.1121.763(1.414–2.198)<0.001
rs3790844 (GG vs. GA/AA)−0.2430.0990.784(0.646–0.951)0.014
rs3790843 (TT vs. CT/CC)−0.1520.1010.859(0.705–1.046)0.131
β, regression coefficient; SE, standard error; HR, hazard ratio; CI, confidence interval. a: Age was included as a continuous variable in the Cox stepwise regression analysis; b: Adjusted by age and sex.
Tumor growth is the result of uncontrolled cell proliferation or a defective cell death program. For the proliferation of cells, the phosphorylation of pRb by cyclin–CDK complexes to release the transcription factor E2F is an essential step. The cyclin D1, cyclin E and c-Myc are important parts of expression peaks at the G1–S transition, and then decreases as cells proceed through the S phase [21,22,23]. The G1–S transition is the major regulation point of the cell cycle, and during the G1–S transition, alterations in cell cycle regulators lead to the deregulation of the cell cycle, which can cause unbridled cell division, contributing to cancer development. Previously, overexpression of cyclin E was demonstrated in many tumors and correlated with prognosis [24,25,26,27] including gastric cancer [28]. Wang et al. [29] found that the mRNA expression of NR5A2 was significantly upregulated in gastric cancer, as compared with self-paired normal control. In addition, overexpression of NR5A2 was shown to promote the proliferation of gastric adenocarcinoma SGC-7901 cells via induction of cyclin E1, which may lead to the tumorigenesis of gastric cancer [29]. Botrugno et al. [20] found that cyclin E was a direct NR5A2 target gene, which could partly explain the mechanism of NR5A2 polymorphism in the progression and prognoses of gastric cancer. On the other hand, β-catenin has been proposed to act as a docking protein that assembles both general and specific factors required for the activation of target genes. Botrugno et al. also found that NR5A2 could act as a coactivator for β-Catenin/Tcf4 to drive the expression of cyclin D1 and other β-catenin/Tcf target genes, such as c-Myc [30,31,32]. As reported, cyclin D1 is a proto-oncogene that belongs to the family of G1 cyclins, and plays an important role in cell cycle G1 to S transition by binding its partners cyclin dependent kinase 4 and 6 to phosphorylate and inactivate the Rb protein [33]. In addition, cyclin D1 over-expression is usually an early event in carcinogenesis and a prognostic indicator associated with poor survival in cancers [34]. Given above-mentioned evidences, NR5A2 polymorphism could also affect the carcinogenesis and prognoses of cancer through acting as a coactivator for β-catenin on the cyclin D1 promoter.
Interestingly, in the stratified analyses by sex, our results indicated that both the NR5A2 rs3790843 and rs3790844 polymorphisms were associated with significantly better prognosis in the female patient group. NR5A2 is a direct estrogen receptor α (ERα) target gene [35,36], its expression correlates with ERα in breast tumours [37] and it promotes breast cancer proliferation and invasion [38]. Previous study showed that NR5A2 is an important regulator of ERα target genes and it shares many binding sites with ERα. Importantly, at shared sites, NR5A2 promotes ERα recruitment and vice versa, ERα stimulates NR5A2 recruitment [39]. As a result, in the female patients with more ERα, more NR5A2 polymorphisms, which act as a protective factor, contribute to a better prognosis.
Taken together, the presented findings of the potential involvement of the NR5A2 gene in anti-tumorigenesis prompts us to further characterize its structure, biological function and interaction with other partners by in vitro or in vivo studies. However, some limitations of the present study should be addressed. First, Helicobacter pylorus, a known crucial factor in gastric carcinogenesis, was not considered due to the lack of related follow-up information. Second, only two SNPs in NR5A2 are evaluated, and it is possible that some other important SNPs are neglected or the observed associations may be due to other polymorphisms in linkage disequilibrium with the rs3790843 and rs3790844 polymorphisms. Finally, for validation of the genotype–phenotype relationship, further investigation is underway to clarify the association between rs3790843 and rs3790844 polymorphisms and expression levels of NR5A2 protein in gastric cancer tissues, and will be reported separately.

4. Experimental Section

4.1. Study Subjects

A retrospective cohort of 944 patients with gastric cancer who underwent a surgical resection at the Yixing People’s Hospital (Yixing, Jiangsu Province, China) between 1999 and 2006, who had a median follow-up of 35 months (range, 0–119 months) were recruited for this study [40]. The demographic features and clinicopathologic data were summarized in Table 1. All patients were diagnosed with gastric carcinoma histopathologically. None had received neoadjuvant radiochemotherapy or postoperative radiotherapy, and some had received adjuvant chemotherapy. Formalin-fixed and paraffin-embedded blocks of these patients were collected from the department of pathology of this hospital. The samples used for genotyping were reviewed and classified by two independent pathologists. Death dates were confirmed via review of death certificates of inpatient and outpatient records or obtained from patients’ families through follow-up telephone calls. The study protocol was approved by the Institutional Review Board of Nanjing Medical University (Nanjing, China), and all patients provided written informed consent on the use of paraffin specimens for gene polymorphism analyses.

4.2. Genotyping

Genomic DNA was extracted from tumor specimens using proteinase K digestion, followed by isopropanol extraction and ethanol precipitation. The NR5A2 (rs3790844 and rs3790843) SNPs were examined by multiplex snapshot technology using an ABI fluorescence-based assay allelic discrimination method (Applied Biosystems, Foster City, CA, USA), which has been described indetail in a previous study [41]. The primers were designed to anneal immediately adjacent to the nucleotideat the mutation site: rs3790844: forward, 5'-TTCCGTGTGGAAACACAGGTCA-3'; reverse, 5'-TCGACTGGAGCCCAAGATCA-3'; rs3790843: forward, 5'-TCTTTGCCCCGATGAGTTCG-3'; reverse, 5'-TCCCAGATGCTCTGGTGCAG-3'. The primers for extension were as follows: rs3790844, 5'-TTTTTTTTTTTTT…TTTTTTTTTTTTTTTTTGGAAACACAGGTCACTAAAACTGG-3'; and rs3790843, 5'-TCTGGTGCAGCCGAAGTAG-3'. The snapshot products were analyzed by using an ABI 3130 genetic analyzer (Applied Biosystems) and the genotypes were determined by GeneMapper Analysis Software version 4 (Applied Biosystems). Genotype analysis was performed by two investigators blinded to the survival end points. Genotyping was validated by sequencing a randomly selected 10% of samples, and the results were 100% concordant. However, 37 cases of rs3790843 and 32 cases of rs3790844 failed in genotyping because of poor DNA quality, which were excluded in further analysis. Finally, 907 cases of rs3790843 and 912 cases of rs3790844 were included in the analysis.

4.3. Statistical Methods

The associations of each genotype or combinations of genotypes with clinicopathologic features were compared using the Pearson chi-square test for categorical variables and the Student’s t-test for continuous data. Associations of genetic variants and clinicopathological features with the overall survival (OS) were estimated using Kaplan-Meier method and comparisons between different groups of patients were performed with the log-rank test. Univariate or multivariate Cox proportional hazard models was used to estimate the adjusted HR and their 95% CI. Moreover, Cox stepwise regression analysis was conducted to assess the independent impacts of SNP or clinicopathological features on the overall survival after adjusting for other covariates, with a significance level of p < 0.05 for entering and p > 0.10 for removing the respective explanatory variables. All tests were two-sided and p < 0.05 was considered statistically significant. All the statistical analyses were carried out using SPSS version 16.0 for Windows (SPSS Inc., Chicago, IL, USA).

5. Conclusions

In conclusion, our results represent the first demonstration that NR5A2 rs3790844 polymorphism was associated with the increased OS of GC patients in the dominant model, and similar results were found among the female patients group and tumor size >5 cm group for the NR5A2 rs3790843 polymorphism, suggesting that the mutant alleles may serve as suitable markers for predicting the survival of gastric cancer patients, especially in a Chinese population. Large population-based prospective studies with ethnically diverse populations are warranted to verify these findings.

Acknowledgments

The work was partly supported by grants from the National Natural Science Foundation of China (Grant No. 81272469), the National 973 Basic Research Program of China (Grant No. 2013CB911300) and the clinical special project for Natural Science Foundation of Jiangsu Province (Grant No. BL2012016), and a grant from The Project of Plans for the Development of Science and Technology of Nanjing, China (Grant No. 201208020), and the grant from Nanjing 12th Five-Year Key Scientific Project of Medicine.

Author Contributions

Xunlei Zhang, Mulong Du and Meilin Wang designed the study and applied for Research Ethics Board Approval; Chunxiang Cao, Lili Shen, Meng Kuang, Yongfei Tan and Xinying Huo recruited the patients and collected the data; Yongfei Tan, Weida Gong and Zhi Xu analyzed the data and prepared draft figures and tables. Zhi Xu and Dongying Gu prepared the manuscript draft. All authors approved the final article. Zhengdong Zhang, Cuiju Tang and Jinfei Chen had complete access to the study data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jemal, A.; Bray, F.; Center, M.M.; Ferlay, J.; Ward, E.; Forman, D. Global cancer statistics. CA Cancer J. Clin. 2011, 61, 69–90. [Google Scholar] [CrossRef] [PubMed]
  2. Parkin, D.M.; Bray, F.; Ferlay, J.; Pisani, P. Global cancer statistics, 2002. CA Cancer J. Clin. 2005, 55, 74–108. [Google Scholar] [CrossRef]
  3. Santoro, R.; Carboni, F.; Lepiane, P.; Ettorre, G.M.; Santoro, E. Clinicopathological features and prognosis of gastric cancer in young european adults. Br. J. Surg. 2007, 94, 737–742. [Google Scholar] [CrossRef]
  4. Correa, P. Human gastric carcinogenesis: A multistep and multifactorial process—First american cancer society award lecture on cancer epidemiology and prevention. Cancer Res. 1992, 52, 6735–6740. [Google Scholar] [PubMed]
  5. Resende, C.; Ristimaki, A.; Machado, J.C. Genetic and epigenetic alteration in gastric carcinogenesis. Helicobacter 2010, 15 (S1), 34–39. [Google Scholar] [CrossRef]
  6. Ludwig, J.A.; Weinstein, J.N. Biomarkers in cancer staging, prognosis and treatment selection. Nat. Rev. Cancer 2005, 5, 845–856. [Google Scholar] [CrossRef] [PubMed]
  7. Panani, A.D. Cytogenetic and molecular aspects of gastric cancer: Clinical implications. Cancer Lett. 2008, 266, 99–115. [Google Scholar] [CrossRef] [PubMed]
  8. Becker, K.F.; Keller, G.; Hoefler, H. The use of molecular biology in diagnosis and prognosis of gastric cancer. Surg. Oncol. 2000, 9, 5–11. [Google Scholar] [CrossRef] [PubMed]
  9. Fayard, E.; Auwerx, J.; Schoonjans, K. LRH-1: An orphan nuclear receptor involved in development, metabolism and steroidogenesis. Trends Cell Biol. 2004, 14, 250–260. [Google Scholar] [CrossRef] [PubMed]
  10. Fernandez-Marcos, P.J.; Auwerx, J.; Schoonjans, K. Emerging actions of the nuclear receptor LRH-1 in the gut. Biochim. Biophys. Acta 2011, 1812, 947–955. [Google Scholar] [CrossRef] [PubMed]
  11. Venteclef, N.; Jakobsson, T.; Ehrlund, A.; Damdimopoulos, A.; Mikkonen, L.; Ellis, E.; Nilsson, L.M.; Parini, P.; Janne, O.A.; Gustafsson, J.A.; et al. GPS2-dependent corepressor/SUMO pathways govern anti-inflammatory actions of LRH-1 and LXRβ in the hepatic acute phase response. Genes Dev. 2010, 24, 381–395. [Google Scholar] [CrossRef] [PubMed]
  12. Clyne, C.D.; Kovacic, A.; Speed, C.J.; Zhou, J.; Pezzi, V.; Simpson, E.R. Regulation of aromatase expression by the nuclear receptor LRH-1 in adipose tissue. Mol. Cell. Endocrinol. 2004, 215, 39–44. [Google Scholar] [CrossRef] [PubMed]
  13. Clyne, C.D.; Speed, C.J.; Zhou, J.; Simpson, E.R. Liver receptor homologue-1 (LRH-1) regulates expression of aromatase in preadipocytes. J. Biol. Chem. 2002, 277, 20591–20597. [Google Scholar] [CrossRef] [PubMed]
  14. Lazarus, K.A.; Wijayakumara, D.; Chand, A.L.; Simpson, E.R.; Clyne, C.D. Therapeutic potential of liver receptor homolog-1 modulators. J. Steroid Biochem. Mol. Biol. 2012, 130, 138–146. [Google Scholar] [CrossRef] [PubMed]
  15. Benod, C.; Vinogradova, M.V.; Jouravel, N.; Kim, G.E.; Fletterick, R.J.; Sablin, E.P. Nuclear receptor liver receptor homologue 1 (LRH-1) regulates pancreatic cancer cell growth and proliferation. Proc. Natl. Acad. Sci. USA 2011, 108, 16927–16931. [Google Scholar] [CrossRef] [PubMed]
  16. Petersen, G.M.; Amundadottir, L.; Fuchs, C.S.; Kraft, P.; Stolzenberg-Solomon, R.Z.; Jacobs, K.B.; Arslan, A.A.; Bueno-de-Mesquita, H.B.; Gallinger, S.; Gross, M.; et al. A genome-wide association study identifies pancreatic cancer susceptibility loci on chromosomes 13q22.1, 1q32.1 and 5p15.33. Nat. Genet. 2010, 42, 224–228. [Google Scholar] [CrossRef] [PubMed]
  17. Schoonjans, K.; Dubuquoy, L.; Mebis, J.; Fayard, E.; Wendling, O.; Haby, C.; Geboes, K.; Auwerx, J. Liver receptor homolog 1 contributes to intestinal tumor formation through effects on cell cycle and inflammation. Proc. Natl. Acad. Sci. USA 2005, 102, 2058–2062. [Google Scholar] [CrossRef] [PubMed]
  18. Rausa, F.M.; Galarneau, L.; Belanger, L.; Costa, R.H. The nuclear receptor fetoprotein transcription factor is coexpressed with its target gene HNF-3β in the developing murine liver, intestine and pancreas. Mech. Dev. 1999, 89, 185–188. [Google Scholar] [CrossRef]
  19. Galarneau, L.; Pare, J.F.; Allard, D.; Hamel, D.; Levesque, L.; Tugwood, J.D.; Green, S.; Belanger, L. The α1-fetoprotein locus is activated by a nuclear receptor of the drosophila FTZ-F1 family. Mol. Cell. Biol. 1996, 16, 3853–3865. [Google Scholar] [PubMed]
  20. Botrugno, O.A.; Fayard, E.; Annicotte, J.S.; Haby, C.; Brennan, T.; Wendling, O.; Tanaka, T.; Kodama, T.; Thomas, W.; Auwerx, J.; et al. Synergy between LRH-1 and β-catenin induces G1 cyclin-mediated cell proliferation. Mol. Cell 2004, 15, 499–509. [Google Scholar] [CrossRef] [PubMed]
  21. Lew, D.J.; Dulic, V.; Reed, S.I. Isolation of three novel human cyclins by rescue of G1 cyclin (cln) function in yeast. Cell 1991, 66, 1197–1206. [Google Scholar] [CrossRef] [PubMed]
  22. Sherr, C.J. Mammalian G1 cyclins. Cell 1993, 73, 1059–1065. [Google Scholar] [CrossRef] [PubMed]
  23. Wimmel, A.; Lucibello, F.C.; Sewing, A.; Adolph, S.; Muller, R. Inducible acceleration of G1 progression through tetracycline-regulated expression of human cyclin E. Oncogene 1994, 9, 995–997. [Google Scholar] [PubMed]
  24. Florenes, V.A.; Faye, R.S.; Maelandsmo, G.M.; Nesland, J.M.; Holm, R. Levels of cyclin D1 and D3 in malignant melanoma: Deregulated cyclin D3 expression is associated with poor clinical outcome in superficial melanoma. Clin. Cancer Res. 2000, 6, 3614–3620. [Google Scholar] [PubMed]
  25. Richter, J.; Wagner, U.; Kononen, J.; Fijan, A.; Bruderer, J.; Schmid, U.; Ackermann, D.; Maurer, R.; Alund, G.; Knonagel, H.; et al. High-throughput tissue microarray analysis of cyclin E gene amplification and overexpression in urinary bladder cancer. Am. J. Pathol. 2000, 157, 787–794. [Google Scholar] [CrossRef] [PubMed]
  26. Kamai, T.; Takagi, K.; Asami, H.; Ito, Y.; Oshima, H.; Yoshida, K.I. Decreasing of p27kip1 and cyclin E protein levels is associated with progression from superficial into invasive bladder cancer. Br. J. Cancer 2001, 84, 1242–1251. [Google Scholar] [CrossRef] [PubMed]
  27. Heah, K.G.; Hassan, M.I.; Huat, S.C. p53 Expression as a marker of microinvasion in oral squamous cell carcinoma. Asian Pac. J. Cancer Prev. 2011, 12, 1017–1022. [Google Scholar] [PubMed]
  28. Ahn, M.J.; Kim, B.H.; Jang, S.J.; Hong, E.K.; Lee, W.M.; Baik, H.K.; Park, H.K.; Lee, C.B.; Ki, M. Expression of cyclin D1 and cyclin E in human gastric carcinoma and its clinicopathologic significance. J. Korean Med. Sci. 1998, 13, 513–518. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, S.L.; Zheng, D.Z.; Lan, F.H.; Deng, X.J.; Zeng, J.; Li, C.J.; Wang, R.; Zhu, Z.Y. Increased expression of hLRH-1 in human gastric cancer and its implication in tumorigenesis. Mol. Cell. Biochem. 2008, 308, 93–100. [Google Scholar] [CrossRef] [PubMed]
  30. He, T.C.; Sparks, A.B.; Rago, C.; Hermeking, H.; Zawel, L.; da Costa, L.T.; Morin, P.J.; Vogelstein, B.; Kinzler, K.W. Identification of c-myc as a target of the APC pathway. Science 1998, 281, 1509–1512. [Google Scholar] [CrossRef] [PubMed]
  31. Shtutman, M.; Zhurinsky, J.; Simcha, I.; Albanese, C.; D’Amico, M.; Pestell, R.; Ben-Ze’ev, A. The cyclin D1 gene is a target of the β-catenin/LEF-1 pathway. Proc. Natl. Acad. Sci. USA 1999, 96, 5522–5527. [Google Scholar] [CrossRef] [PubMed]
  32. Tetsu, O.; McCormick, F. β-catenin regulates expression of cyclin D1 in colon carcinoma cells. Nature 1999, 398, 422–426. [Google Scholar] [CrossRef] [PubMed]
  33. Kato, J.; Matsushime, H.; Hiebert, S.W.; Ewen, M.E.; Sherr, C.J. Direct binding of cyclin D to the retinoblastoma gene product (pRb) and pRb phosphorylation by the cyclin D-dependent kinase CKD4. Genes Dev. 1993, 7, 331–342. [Google Scholar] [CrossRef] [PubMed]
  34. Kim, J.K.; Diehl, J.A. Nuclear cyclin D1: An oncogenic driver in human cancer. J. Cell. Physiol. 2009, 220, 292–296. [Google Scholar] [CrossRef] [PubMed]
  35. Annicotte, J.S.; Chavey, C.; Servant, N.; Teyssier, J.; Bardin, A.; Licznar, A.; Badia, E.; Pujol, P.; Vignon, F.; Maudelonde, T.; et al. The nuclear receptor liver receptor homolog-1 is an estrogen receptor target gene. Oncogene 2005, 24, 8167–8175. [Google Scholar] [PubMed]
  36. Thiruchelvam, P.T.; Lai, C.F.; Hua, H.; Thomas, R.S.; Hurtado, A.; Hudson, W.; Bayly, A.R.; Kyle, F.J.; Periyasamy, M.; Photiou, A.; et al. The liver receptor homolog-1 regulates estrogen receptor expression in breast cancer cells. Breast Cancer Res. Treat. 2011, 127, 385–396. [Google Scholar] [CrossRef] [PubMed]
  37. Miki, Y.; Clyne, C.D.; Suzuki, T.; Moriya, T.; Shibuya, R.; Nakamura, Y.; Ishida, T.; Yabuki, N.; Kitada, K.; Hayashi, S.; et al. Immunolocalization of liver receptor homologue-1 (LRH-1) in human breast carcinoma: Possible regulator of insitu steroidogenesis. Cancer Lett. 2006, 244, 24–33. [Google Scholar] [CrossRef] [PubMed]
  38. Chand, A.L.; Herridge, K.A.; Thompson, E.W.; Clyne, C.D. The orphan nuclear receptor LRH-1 promotes breast cancer motility and invasion. Endocr. Relat. Cancer 2010, 17, 965–975. [Google Scholar] [CrossRef] [PubMed]
  39. Lai, C.F.; Flach, K.D.; Alexi, X.; Fox, S.P.; Ottaviani, S.; Thiruchelvam, P.T.; Kyle, F.J.; Thomas, R.S.; Launchbury, R.; Hua, H.; et al. Co-regulated gene expression by oestrogen receptor alpha and liver receptor homolog-1 is a feature of the oestrogen response in breast cancer cells. Nucleic Acids Res. 2013, 41, 10228–10240. [Google Scholar] [CrossRef] [PubMed]
  40. Wang, M.; Bai, J.; Tan, Y.; Wang, S.; Tian, Y.; Gong, W.; Zhou, Y.; Gao, Y.; Zhou, J.; Zhang, Z. Genetic variant in PSCA predicts survival of diffuse-type gastric cancer in a Chinese population. Int. J. Cancer 2011, 129, 1207–1213. [Google Scholar] [CrossRef] [PubMed]
  41. Zhang, Y.; Wang, M.; Gu, D.; Wu, D.; Zhang, X.; Gong, W.; Tan, Y.; Zhou, J.; Wu, X.; Tang, C.; et al. Association of XRCC1 gene polymorphisms with the survival and clinicopathological characteristics of gastric cancer. DNA Cell Biol. 2013, 32, 111–118. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Zhang, X.; Gu, D.; Du, M.; Wang, M.; Cao, C.; Shen, L.; Kuang, M.; Tan, Y.; Huo, X.; Gong, W.; et al. Associations of NR5A2 Gene Polymorphisms with the Clinicopathological Characteristics and Survival of Gastric Cancer. Int. J. Mol. Sci. 2014, 15, 22902-22917. https://doi.org/10.3390/ijms151222902

AMA Style

Zhang X, Gu D, Du M, Wang M, Cao C, Shen L, Kuang M, Tan Y, Huo X, Gong W, et al. Associations of NR5A2 Gene Polymorphisms with the Clinicopathological Characteristics and Survival of Gastric Cancer. International Journal of Molecular Sciences. 2014; 15(12):22902-22917. https://doi.org/10.3390/ijms151222902

Chicago/Turabian Style

Zhang, Xunlei, Dongying Gu, Mulong Du, Meilin Wang, Chunxiang Cao, Lili Shen, Meng Kuang, Yongfei Tan, Xinying Huo, Weida Gong, and et al. 2014. "Associations of NR5A2 Gene Polymorphisms with the Clinicopathological Characteristics and Survival of Gastric Cancer" International Journal of Molecular Sciences 15, no. 12: 22902-22917. https://doi.org/10.3390/ijms151222902

APA Style

Zhang, X., Gu, D., Du, M., Wang, M., Cao, C., Shen, L., Kuang, M., Tan, Y., Huo, X., Gong, W., Xu, Z., Chen, J., Zhang, Z., & Tang, C. (2014). Associations of NR5A2 Gene Polymorphisms with the Clinicopathological Characteristics and Survival of Gastric Cancer. International Journal of Molecular Sciences, 15(12), 22902-22917. https://doi.org/10.3390/ijms151222902

Article Metrics

Back to TopTop