A Case-Control Study of the Association Between the SPP1 Gene SNPs and the Susceptibility to Breast Cancer in Guangxi, China

Secreted phosphoprotein-1 (SPP1) has been reported to be involved in the pathogenesis of breast cancer (BRC), but the influence of SPP1 single nucleotide polymorphisms on the BRC susceptibility has been rarely reported. In this study, we explored the association between rs11730582, rs2853750, and rs35893069 in the SPP1 gene and the BRC susceptibility. We used Snapshot assay to detect SPP1 single nucleotide polymorphisms in 471 BRC patients and 471 controls. The plasma SPP1 level was measured by ELISA. We found that the CC genotype and C allele of rs11730582 were associated with a significantly decreased BRC risk compared with the TT genotype and T allele, respectively [CC vs. TT: odds ratio (OR) = 0.59, 95% CI = 0.37–0.94, P = 0.026; C vs. T: OR = 0.79, 95% CI = 0.65–0.96, P = 0.022]. In addition, BRC patients and controls with the rs11730582 CC genotype had a lower plasma SPP1 level than did BRC patients and controls with TT genotype (P = 0.007 and P = 0.011, respectively). Moreover, the proportions of rs11730582 CC genotype and C allele were decreased in BRC patients with clinical stages I–III compared with those with clinical stage IV (P = 0.012 and P = 0.003, respectively). Besides, the C-G-T haplotype was associated with a significantly decreased BRC risk compared with the T-A-T haplotype (OR = 0.69, 95% CI = 0.52–0.93, P = 0.015). However, there was no significant association between rs2853750 or rs35893069 and the BRC risk. In summary, our study found the association between rs11730582 and the risk of BRC and suggested that rs11730582 may promote the occurrence and development of BRC by regulating SPP1 expression.

Introduction

Breast cancer (BRC) is the most common cancer type and a leading cause of mortality among females in the world (1). According to the statistics from the American Cancer Society, the BRC incidence in the United States was 46.3 per 100 thousand females and the mortality was 13.0 per 100 thousand females (2). In India, BRC ranked first among different types of cancer, with the incidence of 25.8 per 100 thousand females and the mortality of 12.7 per 100 thousand females (3). China has the largest number of BRC patients in the world. About 278 thousand females were diagnosed as having BRC in 2014; the incidence was 30.6 per 100 thousand females and the mortality was 6.5 per 100 thousand females (4, 5). Because of the high incidence and mortality, BRC has brought numerous families in the world a great economic and health burden (6, 7). The etiology of BRC is involved in multiple factors including smoking, diet, estrogen exposure, irregularity of menstruation periods, family history of BRC, and so on (8, 9). However, these risk factors just clarify some parts of the BRC pathogenesis, and the underlying etiology of BRC still needs to be further elucidated.

Secreted phosphoprotein-1 (SPP1), also known as osteopontin, is an extracellular matrix protein that is produced in multiple human tissues and involved in the control of biomineralization and calcification (10, 11). Some studies suggested that SPP1 played an important role in the progression and metastasis of several types of cancers. In a study about ovarian cancer, SPP1 was found to be overexpressed in ovarian cancer tissues and promote ovarian cancer progression (12). The overexpression of SPP1 was also found in colorectal cancer tissues and related to poor prognosis (13). Recently, Pio et al. reported that SPP1 was up-regulated in the BRC cells and promoted their migration and stem-like behavior, potentially through activation of the WNK-1 and PRAS40-related pathways (14). Similarly, SPP1 was also found to be increased in normal human breast tissue at a high risk of developing BRC (15) and promote inflammation and tumor growth by reprograming normal mammary fibroblasts (16). The studies above suggested that SPP1 was a potential biomarker or a novel therapeutic target for BRC.

Previous works had shown that single nucleotide polymorphisms (SNPs) in the SPP1 gene were associated with the risk or metastasis of some human cancers, such as oral carcinogenesis (17), hepatocellular carcinoma (18), and BRC (19). Particularly, the association between SNPs in the SPP1 gene and the susceptibility to BRC had been evaluated among the population in Egypt, but the sample size was small (60 cases and 60 controls) (20). A further study with a larger sample size should be conducted immediately to reveal the association between SNPs in the SPP1 gene and the susceptibility to BRC. Therefore, we conducted this study among the population in Guangxi, China. In our study, we genotyped SNPs in the SPP1 gene, measured its expression level, and assessed the impact of SNPs in the SPP1 gene on the expression level of SPP1 and other BRC-related characteristics.

Materials and Methods

Study Population

Four hundred seventy-one histologically confirmed BRC patients were recruited in the Department of Breast and Thyroid Surgery, Affiliated Hospital of Youjiang Medical University for Nationalities from February 2017 to November 2018. Clinical information was collected from medical records including age at diagnosis, age at menarche, estrogen receptor status, progesterone receptor status, and clinical stages. Patients who had a family history of cancer, recurring BRC, or BRC combined with other types of cancers were excluded. Four hundred seventy-one controls were healthy volunteers visiting the same hospital from December 2017 to December 2018. We excluded those controls who had a family history of cancer, breast lesions, or breast mastitis. All the subjects came from Guangxi, China. About 2–3 ml anticoagulant whole blood sample was taken from each participant. The study was approved by the Ethics Committee of Youjiang Medical University for Nationalities, conformed to the Declaration of Helsinki, and received written informed consent from all subjects before the study.

SNP Selection and Genotyping

We used Promoter 2.0 Prediction Server (http://www.cbs.dtu.dk/services/Promoter/) and The Eukaryotic Promoter Database (http://www.epd.isb-sib.ch/) to predict promoter. Subsequently, SNPs with minor allele frequency (MAF) >30% among the Southern Han Chinese population were selected by searching in the promoter and open reading frame of the SPP1 gene in the National Center for Biotechnology Information SNP database (http://www.ncbi.nlm.nih.gov/SNP) and 1000 Genomes Project (http://browser.1000genomes.org). Finally, we evaluated the association between three SNPs (rs11730582, rs2853750, and rs35893069) in the SPP1 gene and the susceptibility to BRC. The basic information of the three SNPs is listed in Supplementary Table 1.

DNA extraction kit (Tiangen, China) was used to extract DNA from blood samples and then stored the DNA at −80°C. Primers were designed by Primer 5 (PREMIER Biosoft) and synthesized by Sangon Company, and are listed in Supplementary Table 2. After amplifying the SPP1 gene by polymerase chain reaction, the three SNPs were genotyped by Snapshot assay.

Plasma SPP1 Level Determination

Plasma was separated from whole blood samples of patients and controls at room temperature and stored at −70°C until analysis. Plasma SPP1 level was determined according to the instruction of ELISA kit (BioVision, Catalog: K7335). The developed color reaction was measured at 450-nm wavelength on the ELISA reader (Bio-Rad, USA), and detection range was from 0.2 to 70 ng/ml.

Statistics

Statistics were mainly proceeded by Statistical Package for Social Science 21. The percentages were calculated for clinical features and different genotypes and alleles, and the means and SDs were calculated for age at diagnosis and age at menarche. The data normality for continuous variables was assessed by Shapiro-Wilk tests, then the continuous variables were analyzed by the independent-sample t-test. We performed the Hardy-Weinberg equation for the three SNPs in our study. Two-sided chi-square test was used to assess the conformity to the Hardy-Weinberg rule of the three SNPs genotypes. Categorical variables comparisons were also proceeded by two-sided chi-square test. The associations of the three SNPs with BRC susceptibility were assessed in genotypic, allelic, recessive (rs11730582: CC + TC vs. TT; rs2853750: GG + AG vs. AA; rs35893069: TT + AT vs. AA), and dominant model (rs11730582: CC vs. TC + TT; rs2853750: GG vs. AG + AA; rs35893069: TT vs. AT + AA). Logistic regression was used for exploring the association between SPP1 SNPs and the susceptibility to BRC, and P values, odds ratio (OR), and 95% CI were adjusted by age at diagnosis and age at menarche. Genetical distribution of rs11730582 among 11 different populations were collected on the website below (http://www.internationalgenome.org). Haplotype analysis was proceeded by SHEsis Main (http://analysis.bio-x.cn/). P < 0.05 was regarded as statistically significant. The GAS Power Calculator was used for power calculation, which was proceeded by rs11730582 with the minimum MAF among the three SNPs. Under the conditions of sample size (cases = 471, controls = 471), disease prevalence = 0.00041 (4), significance level = 0.05, MAF = 0.31, and OR = 1.5), we had 98% statistical power in this study.

Results

Clinical Features of Cases and Controls

Clinical features of cases and controls are shown in Table 1. As for the distribution of age at diagnosis and age at menarche, there was no significant difference between BRC patients and the controls (P = 0.630 and P = 0.596, respectively). The proportions of cases with positive estrogen receptor, positive progesterone receptor, and III–IV clinical stage were 61.9, 55.0, and 37.1%, respectively.

TABLE 1

Table 1. Clinical features of cases and controls.

Association Between the Three SPP1 SNPs and the Risk of BRC

As shown in Supplementary Table 1, the genotypes of rs11730582, rs2853750, and rs35893069 conformed to the Hardy-Weinberg rule (P = 0.685, P = 0.973, and P = 0.966, respectively). The association between the three SPP1 SNPs (rs11730582, rs2853750, and rs35893069) and the risk of BRC is shown in Table 2. We found that the rs11730582 CC genotype was associated with a significantly decreased BRC risk compared with the TT genotype (OR = 0.59, 95% CI = 0.37–0.94, P = 0.026). Similarly, we also observed a significantly decreased risk in allele analysis (OR = 0.79, 95% CI = 0.65–0.96, P = 0.022). However, no significant association of the rs2853750 and rs35893069 with the risk of BRC was found (P > 0.05).

TABLE 2

Table 2. Association between SPP1 SNPs and the risk of BRC.

Association Between rs11730582 and Plasma SPP1 Level

Association between rs11730582 and plasma SPP1 level is shown in Figure 1. Plasma SPP1 level was significantly up-regulated in BRC patients compared with controls (P < 0.001). Moreover, we found that BRC patients with the rs11730582 CC genotype had a lower plasma SPP1 level than did those with the TT genotype (P = 0.007). Similarly, controls with the rs11730582 CC genotype also had a lower plasma SPP1 level than did those with the TT genotype (P = 0.011). Nevertheless, there was no significant difference between individuals with the rs11730582 CC genotype and those with the TC genotype on plasma SPP1 level.

FIGURE 1

Figure 1. Plasma secreted phosphoprotein-1 (SPP1) level in controls and cases. (A) Decreased plasma SPP1 level was found in breast cancer (BRC) patients compared to controls (P < 0.001). (B) Decreased plasma SPP1 level was found in controls carrying the rs11730582 CC genotype compared with those carrying the rs11730582 TT (P = 0.011). (C) Decreased plasma SPP1 level was found in BRC patients carrying the rs11730582 CC genotype compared with those carrying the rs11730582 TT (P = 0.007). The upper and lower boundaries of the boxes represented the 25th and 75th percentiles of plasma SPP1 level value, and the horizontal line within the box represents the median value. All P-values were adjusted by age at diagnosis and age at menarche.

Association Between rs11730582 and BRC Clinical Features

Association between rs11730582 and BRC clinical features is shown in Table 3, and the proportion of rs11730582 CC genotype was decreased in BRC patients with clinical stages I–III compared with those with clinical stage IV (P = 0.012). Besides, the rs11730582 A allele was also decreased in BRC patients with clinical stages I–III compared with those with clinical stage IV (P = 0.003). However, there was no significant difference between rs11730582 and estrogen receptor status or progesterone receptor status of BRC patients.

TABLE 3

Table 3. The associations between rs11730582 and clinical features of BRC.

Haplotype Analysis

Haplotype analysis of rs11730582-rs2853750-rs35893069 with the BRC risk is shown in Table 4; six haplotypes were listed, rs11730582 was in a strong linkage disequilibrium with rs2853750 (D′ = 0.981) and rs35893069 (D′ = 0.810), and rs2853750 was also in a strong linkage disequilibrium with the rs35893069 (D′ = 0.849). Moreover, T-A-T and T-A-A were two major haplotypes in both controls (30.5 and 24.6%, respectively) and BRC patients (31.7 and 24.3%, respectively). In addition, we observed that the C-G-T haplotype was associated with a significantly decreased BRC risk compared with the T-A-T haplotype (OR = 0.69, 95% CI = 0.52–0.93, P = 0.015).

TABLE 4

Table 4. Haplotype analysis of rs11730582-rs2853750-rs35893069 with the BRC risk.

Genetical Distribution of rs11730582 Among 11 Different Populations

Genetical distribution of rs11730582 among 11 different populations is shown in Table 5. The genetical distribution of rs11730582 in this study was significantly different from those among 1KGP-ACB, 1KGP-CLM, 1KGP-ASW, 1KGP-PEL, 1KGP-TSI, 1KGP-PJL, and 1KGP-CEU populations (P < 0.05). In contrast, there was no significant difference after comparing with 1KGP-CHB, 1KGP-CHS, and 1KGP-JPT populations (P > 0.05).

TABLE 5

Table 5. Genotypic and allelic distributions of rs11730582 among different populations.

Discussion

There are millions of females suffering from BRC in the world; it has brought numerous families an unaffordable economic burden (6, 7). SPP1, a tumor-related extracellular matrix protein, has been reported to be involved in the pathogenesis of BRC (10, 11). Pio et al. reported that SPP1 was up-regulated in BRC cells and promoted their migration and stem-like behavior (14). Similar result was found in a study conducted by Lindahl et al. (15). In our study, the upregulation of plasma SPP1 level in BRC patients was also observed and consistent with the results of previous studies.

SNPs have been broadly reported to cause the abnormal functions of various genes. Until now, a few studies had focused on the association of SPP1 SNPs with the susceptibility to some human cancers. Chiu et al. reported that a less prevalence of the rs11730582 CC genotype was found in oral squamous cell carcinoma patients from Taiwan (17). Likewise, Dong et al. reported that hepatocellular carcinoma patients with the rs11730582 CC genotype had a longer overall survival and time to recurrence compared to those with the rs11730582 TT/TC genotype (18). Notably, Ramchandani and Weber corroborated that the rs11730582 C allele was associated with lower aggressiveness of BRC (19). In addition, in a study performed by Zakhary et al., the rs11730582 CC genotype was found to be significantly lower in BRC patients compared to controls among the Egyptian population, but the sample size was too small (60 cases and 60 controls) (20). In our study, we found that the rs11730582 CC genotype and C allele were associated with a significantly decreased BRC risk compared with the TT genotype and T allele, respectively. The results in our study were consistent with the studies above and indicated that rs11730582 may be involved in the BRC etiology. The association of rs2853750 and rs11730582 with the BRC risk was firstly reported. However, we found a negative association of rs2853750 or rs35893069 with the BRC risk in this study. There are some explanations for the negative result. Both rs2853750 and rs35893069 are the SNPs locating in the introns of SPP1 gene. Previous studies showed that an SNP in the intron could influence RNA splicing process during the maturation of mRNA (21, 22), but sometimes this effect was lost (23, 24). In addition, few important gene regulation-related elements locate in the introns of a gene.

In consideration of that rs11730582 located in the promoter of SPP1 gene, we further assessed SPP1 expression with the rs11730582 genotypes. In the present study, BRC patients with the rs11730582 CC genotype had a lower plasma SPP1 level than did those with the TT genotype. Similarly, controls with the rs11730582 CC genotype also had a lower plasma SPP1 level than did those with the TT genotype. Some previous studies supported our results. In a study performed by Ramchandani and Weber, the proportion of rs11730582 CC genotype was found to be decreased in BRC patients and controls with a high SPP1 level (19). Analogously, Zakhary et al. also found that the rs11730582 CC genotype was associated with a lower SPP1 level compared to the TT genotype in BRC patients (20). As we have known, SNPs in the promoter of a gene could alter the promoter transcription activity; we deduced that the rs11730582 genotypes might alter plasma SPP1 level by changing the transcription activity of SPP1 gene. The deduction was verified in a study conducted by Dong et al. (18). In their study, SPP1 level in hepatocellular carcinoma cell transfected with SPP1 promoter containing the rs11730582 CC genotype was significantly lower than that in hepatocellular carcinoma cell transected with the TT genotype.

We further investigated the association between rs11730582 and BRC clinical features. In the current study, the proportions of rs11730582 CC genotype and C allele were decreased in BRC patients with clinical stages I–III compared with those with clinical stage IV. Based on the published studies mentioned above, a significantly increased SPP1 level had been regarded as an important role to promote the development of BRC, and the individuals with rs11730582 CC genotype had a lower plasma SPP1 level than did those with TT genotype. Herein, we concluded that the rs11730582 CC genotype might downregulate the transcription activity of SPP1 gene and plasma SPP1 level, finally promoting the occurrence and development of BRC.

In the present study, we performed a haplotype analysis for rs11730582-rs2853750-rs35893069. We observed that the C-G-T haplotype was associated with a significantly decreased BRC risk compared with the major haplotype T-A-T. Similarly, in a study about sarcoidosis, Maver et al. found that a haplotype with three SPP1 SNPs including rs11730582 significantly decreased the sarcoidosis risk (25). In addition, Glas et al. reported that a haplotype with eight SPP1 SNPs including rs11730582 was strongly associated with the susceptibility to Crohn's disease (26). These studies above indicated that rs11730582-related SPP1 haplotypes might be closely associated with the BRC risk, and also prove the importance of rs11730582 in the BRC pathogenesis.

There are two limitations in our study. First, all cases and controls in our study came from Guangxi; the subjects from other regions of the world should be recruited. As shown in Table 5, the genetical distribution of rs11730582 among in the population of Guangxi was similar to those among 1KGP-CHB, 1KGP-CHS, and 1KGP-JPT populations, but significantly different from those among 1KGP-ACB, 1KGP-CLM, 1KGP-ASW, 1KGP-PEL, 1KGP-TSI, 1KGP-PJL, and 1KGP-CEU populations. The difference in genetical distribution among different populations caused the different susceptibility to BRC. Although our study firstly reported the susceptibility to BRC in the population of Guangxi and offered important evidence for BRC pathogenesis, it was still very important to investigate the susceptibility to BRC in other populations. In the future, further studies in other populations will be performed. Second, we investigated SNPs in the SPP1 gene in our study; the results provided a partial insight into the BRC pathogenesis. More BRC risk factors like epigenetics and so on should be investigated to understand the BRC pathogenesis thoroughly in future studies.

In conclusion, the present study found the association between SPP1 rs11730582 and the BRC risk among the population of Guangxi China, and suggested that rs11730582 may promote the occurrence and development of BRC by regulating SPP1 expression. The results may provide important evidence for the etiology of BRC and a potential BRC biomarker.

Data Availability Statement

The datasets generated for this study can be found in the datasets in the current study are available from the corresponding author on reasonable request.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Youjiang Medical University for Nationalities. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

LinL and GW designed the study and wrote the manuscript. LinL, GL, and GP performed experiments. YD, JLia, and LimL performed statistical analysis. JLiu and YT helped to revise the manuscript.

Funding

This work was supported by the Guangxi Education Department Basic Ability Improvement project (No. 2016YB348), the Baise Scientific Research and Technology Development project (Nos. 20140617 and 20192101), Project of Guangxi Education Department (Nos. KY2015YB232 and KY2016YB342), Project of Guangxi Health committee (No. Z2010092), the First Batch of High-level Talent Scientific Research Projects of the Affiliated Hospital of Youjiang Medical University for Nationalities in 2019 (No. R20196307), and Natural Science Foundation of Guangxi (No. 2019JJB140033).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Zhizhai Luo (Department of Breast and Thyroid Surgery, Affiliated Hospital of Youjiang Medical University for Nationalities) for samples collection.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.01415/full#supplementary-material

References

1. Harbeck N, Gnant M. Breast cancer. Lancet. (2017) 389:1134–50. doi: 10.1016/s0140-6736(16)31891-8

PubMed Abstract | CrossRef Full Text | Google Scholar

2. 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) 61:69–90. doi: 10.3322/caac.21492

CrossRef Full Text | Google Scholar

3. Malvia S, Bagadi SA, Dubey US, Saxena S. Epidemiology of breast cancer in Indian women. Asia Pac J Clin Oncol. (2017) 13:289–95. doi: 10.1111/ajco.12661

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Li H, Zheng R, Zhang S, Zeng H, Sun K, Xia C, et al. Incidence and mortality of female breast cancer in China, 2014. Zhonghua Yu Fang Yi Xue Za Zhi. (2018) 40:166–71. doi: 10.3760/cma.j.issn.0253-3766.2018.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Sun K, Zheng R, Gu X, Zhang S, Zeng H, Zou X, et al. Incidence trend and change in the age distribution of female breast cancer in cancer registration areas of China from 2000 to 2014. Zhonghua Yu Fang Yi Xue Za Zhi. (2018) 52:567–75. doi: 10.3760/cma.j.issn.0253-9624.2018.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ye M, Lu J, Yang F, Wu B. Economic evaluation of letrozole for early breast cancer in a health resource-limited setting. Biomed Res Int. (2018) 2018:e9282646. doi: 10.1155/2018/9282646

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Tartari F, Santoni M, Pistelli M, Berardi R. Healthcare cost of HER2-positive and negative breast tumors in the United States (2012-2035). Cancer Treat Rev. (2017) 60:12–7. doi: 10.1016/j.ctrv.2017.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Sun YS, Zhao Z, Yang ZN, Xu F, Lu HJ, Zhu ZY, et al. Risk factors and preventions of breast cancer. Int J Biol Sci. (2017) 13:1387–97. doi: 10.7150/ijbs.21635

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Winters S, Martin C, Murphy D, Shokar NK. Breast cancer epidemiology, prevention, and screening. Prog Mol Biol Transl Sci. (2017) 151:1–32. doi: 10.1016/bs.pmbts.2017.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Icer MA, Gezmen M. The multiple functions and mechanisms of osteopontin. Clin Biochem. (2018) 59:17–24. doi: 10.1016/j.clinbiochem.2018.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Hao C, Cui Y, Owen S, Li W, Cheng S, Jiang WG. Human osteopontin: potential clinical applications in cancer. Int J Mol Med. (2017) 39:1327–37. doi: 10.3892/ijmm.2017.2964

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zeng B, Zhou M, Wu H, Xiong Z. SPP1 promotes ovarian cancer progression via Integrin beta1/FAK/AKT signaling pathway. Oncol Targets Ther. (2018) 11:1333–43. doi: 10.2147/ott.s154215

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Choe E, Yi J, Chai Y, Park K. Upregulation of the adipokine genes ADIPOR1 and SPP1 is related to poor survival outcomes in colorectal cancer. J Surg Oncol. (2018) 117:1833–40. doi: 10.1002/jso.25078

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Pio G, Xia Y, Piaseczny M, Chu J, Allan A. Soluble bone-derived osteopontin promotes migration and stem-like behavior of breast cancer cells. PLoS ONE. (2017) 12:e0177640. doi: 10.1371/journal.pone.0177640

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lindahl G, Rzepecka A, Dabrosin C. Increased extracellular osteopontin levels in normal human breast tissue at high risk of developing cancer and its association with inflammatory biomarkers in situ. Front Oncol. (2019) 9:e746. doi: 10.3389/fonc.2019.00746

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Sharon Y, Raz Y, Cohen N, Ben-Shmuel A, Schwartz H, Geiger T, et al. Tumor-derived osteopontin reprograms normal mammary fibroblasts to promote inflammation and tumor growth in breast cancer. Cancer Res. (2015) 75:963–73. doi: 10.1158/0008-5472.can-14-1990

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Chiu Y, Tu H, Wang I, Wu C, Chang K, Liu T, et al. The implication of osteopontin (OPN) expression and genetic polymorphisms of OPN promoter in oral carcinogenesis. Oral Oncol. (2010) 46:302–6. doi: 10.1016/j.oraloncology.2010.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Dong Q, Zhang X, Zhao Y, Jia H, Zhou H, Dai C, et al. Osteopontin promoter polymorphisms at locus−443 significantly affect the metastasis and prognosis of human hepatocellular carcinoma. Hepatology. (2013) 57:1024–34. doi: 10.1002/hep.26103

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ramchandani D, Weber G. An osteopontin promoter polymorphism is associated with aggressiveness in breast cancer. Oncol Rep. (2013) 30:1860–8. doi: 10.3892/or.2013.2632

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Zakhary M, Mahmoud A, Hashim M. Role of osteopontin and its rs11730582 gene polymorphism in breast cancer. Meta Gene. (2018) 17:124–7. doi: 10.1016/j.mgene.2018.05.006

CrossRef Full Text | Google Scholar

21. Liu X, Liu H, Wang M, Li R, Zeng J, Mo D, et al. Disruption of the ZBED6 binding site in intron 3 of IGF2 by CRISPR/Cas9 leads to enhanced muscle development in Liang Guang Small Spotted pigs. Transgenic Res. (2019) 28:141-50. doi: 10.1007/s11248-018-0107-9

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Su K, Chen S, Ye J, Kuang L, Zhang T, Wang H, et al. A functional indel polymorphism rs34396413 in TFAP2A intron-5 significantly increases female encephalocele risk in Han Chinese population. Childs Nerv Syst. (2019) 35:965–72. doi: 10.1007/s00381-019-04131-1

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Selinski S, Blaszkewicz M, Ickstadt K, Gerullis H, Otto T, Roth E, et al. Identification and replication of the interplay of four genetic high-risk variants for urinary bladder cancer. Carcinogenesis. (2017) 38:1167–79. doi: 10.1093/carcin/bgx102

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Mahdavipour M, Zarei S, Fatemi R, Edalatkhah H, Heidarivala H, Jedditehrani M, et al. Polymorphisms in the estrogen receptor beta gene and the risk of unexplained recurrent spontaneous abortion. Avicenna J Med Biotechnol. (2017) 9:150–4.

PubMed Abstract | Google Scholar

25. Maver A, Medica I, Salobir B, Tercelj M, Peterlin B. Genetic variation in osteopontin gene is associated with susceptibility to sarcoidosis in Slovenian population. Dis Mark. (2009) 27:295–302. doi: 10.3233/dma-2009-0675

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Glas J, Seiderer J, Bayrle C, Wetzke M, Fries C, Tillack C, et al. The role of osteopontin (OPN/SPP1) haplotypes in the susceptibility to Crohn's disease. PLoS ONE. (2011) 6:e29309. doi: 10.1371/journal.pone.0029309

PubMed Abstract | CrossRef Full Text | Google Scholar