SNP rs12982687 affects binding capacity of lncRNA UCA1 with miR-873-5p: involvement in smoking-triggered colorectal cancer progression

Background This investigation was arranged to elucidate whether single nucleotide polymorphisms (SNPs) of lncRNA UCA1 was implicated in elevating colorectal cancer (CRC) risk by interacting with environmental exposures. Methods LncRNASNP database was firstly adopted to predict SNPs that possibly affected binding of UCA1 with miRNAs and then the interactive effect of SNPs and environmental exposure on CRC risk was evaluated by recurring to type 2 gene-environment interactions (GEI) model. Besides, MTT assay, colony formation assay, transwell assay and wound healing assay were performed to assess the activity of CRC cell lines which carried distinct genotypes of specific SNPs. The impact of nicotine on activity of CRC cells was also appraised. Results SNP rs12982687 of UCA1 intervened in the binding capacity of UCA1 with several miRNAs, especially miR-873-5p. MiRNAs regulated by UCA1, as predicted by mirPath software, shared genes that were enriched in HIF1 signaling pathway. Moreover, homozygote TT of rs12982687 reduced CRC risk among smokers, and CRC cells that carried rs12982687 (CC) displayed strong migration and invasion. By contrast, miR-873-5p mimic, which reduced UCA1 expression, delayed metastasis of CRC cells (all P < 0.05). Additionally, nicotine not merely elevated UCA1 and HIF-1α expressions in CRC cells, but also facilitated proliferation and metastasis of CRC cells (P < 0.05). Conclusions SNP rs12982687 was involved in smoking-triggered CRC progression, given its influence on UCA1's binding with miR-873-5p and HIF-1 signaling.

pathogenesis and figure out effective strategies for CRC diagnosis and treatment.
In addition to possessing protein-binding sites, lncRNAs could also dynamically and specifically acted upon DNAs and RNAs based on their sophisticated secondary structure [5], which altogether caused expressional change of pathogenic genes and ultimately disease progression [6]. Of note, there were huge potential for lncRNAs, despite their inability to encode proteins, to boost or decelerate neoplastic development (e.g. CRC) [7][8][9][10]. For instance, highly-expressed lncRNA CRNDE was an indicator of advanced CRC, and knockout of lncRNA CRNDE could induce in-vivo and in-vitro apoptosis of CRC cells [11,12]. Moreover, incremental expression of CCAT1 and CCAT2 might enable aggressive growth and intrusion of CRC cells [13,14], and lncRNA MIR100HG was able to reinforce cetuximab-tolerance of CRC cells via activation of Wnt/βcatenin signaling [15]. With regard to lncRNA urothelial carcinoma associated 1 (UCA1), its high expression could reflect poor survival of CRC patients, and artificially raising its intracellular expression might powerfully enhance epithelial-mesenchymal transition (EMT) and radioresistance of CRC cells [16,17]. Apart from CRC, the oncogenic function of UCA1 was also embodied in nonsmall cell lung cancer [18], bladder cancer [19], breast cancer [20], tongue squamous cell carcinoma [21] and esophageal cancer [22], a portion of which were achieved through its interplay with specific miRNAs.
There were several regulatory manners between lncRNAs and miRNAs [23,24]. An emerging theory proposed that single nucleotide polymorphisms (SNPs) were competent in altering the secondary structure of lncRNAs [25], which produced sizeable effects on the interaction of lncRNA with corresponding miRNAs and ultimately heightened cancer risk. For instance, SNP rs11752942 (A > G) was estimated to interfere with lincRNA-uc003opf.1 expression through binding with miR-149, which altogether facilitated progression of esophageal squamous cell carcinoma [26]. This rule might also be applicable to lncRNA UCA1, which suggested that certain SNPs of UCA1 counted in regulating cancer (e.g. CRC) risk by interfering with binding of UCA1 to miRNAs [25]. Besides the genetic factors, environmental exposures were also estimated to enlarge cancer risk [27,28]. And there was a speculation that SNPs could interact with environmental exposures to multiply disease risk [29][30][31][32][33]. It was thus hypothesized that risky SNPSs in UCA1 might intensify risk via its interaction with environmental hazards, although relevant molecular explanations remained far from complete.
Hence, this investigation was aimed at locating SNPs that might disturb the binding capacity of UCA1 to CRC-regulating miRNAs. Besides, environmental parameters that aggravated or relieved genetic impacts on CRC risk were also explored, which might conduce to developing effective strategies for diagnosing and treating CRC.

Collection of clinical samples
From April 2012 to July 2013, 762 patients pathologically examined as CRC were recruited from the First Affiliated Hospital of Zhengzhou University. Their severity was classified according to the TNM staging system revised by Union for International Cancer Control (UICC). Both cancer tissues and corresponding normal mucosa (≥10 cm distant from tumor tissues) were excised from the patients during surgery. The included CRC patients all met following requirements: 1) their age ranged from 18 years old to 75 years old; 2) they received none of chemotherapy, targeted therapy or radiotherapy pre-operatively; 3) their follow-up information was complete; and 4) their RNA samples were qualified. Correspondingly, those diagnosed with following diseases would be excluded from this project, including familial adenomatous polyps, other colonic disorders (e.g. ulcerative colitis, colon tuberculosis and Crohn's disease), other primary cancers and immune-deficient disorders (e.g. HIV and AIDS). For another, 1408 healthy Chinese who took medical examinations in the First Affiliated Hospital of Zhengzhou University were included into the control group, and they were without any history of tumors, liver/kidney dysfunction or any other smoking-induced disorders. Last but not the least, this study was approved by the First Affiliated Hospital of Zhengzhou University and its ethics committee, and all the participants have signed informed consents.

Follow-up of CRC patients
The CRC patients were followed up for 4 years through household survey, questionnaire-assisted investigation and telephone. Forty-one CRC patients quitted midway from this project, owing to loss of contact. The followup period started from the time point when diagnosis was confirmed until death of the patient. The survival time of CRC patients ranged from 0.34 month to 50 months, with a median value of 34.29 months.

Selection of SNPs
We screened SNPs of UCAl mainly following four procedures. In the first place, UCAl was positioned through retrieving UCSC human GRCh37hg19 database (http:// genome.ucsc.edu). Secondly, SNPs within UCAl were filtrated from database of 1000 Genomes Project, and ones that accorded with following three conditions were selected: 1) their genetic typing rate was > 95%; b) their minor allele frequency (MAF) was > 0.05; and c) their Hardy-Weinberg equilibrium (HWE) was > l × 10 − 5 . Thirdly, SNPs situated in the intron region of UCAl, which were marked based on NCBI dbSNP database (https//www.ncbi.nlm.nih.gov/projects/SNP), were removed. Finally, SNPs that might affect binding of UCA1 with miRNAs were screened, assisted by lncRNASNP online database (http://bioinfo.1ife.hust.edu.cn/lncRNASNP).

Cell culture
A couple of human CRC cell lines (i.e. SW480 and SW620) were purchased from the Chinese Academy of Sciences (Shanghai, China). They were cultured in 5% CO 2 and 95% humidity at 37°C, and the culture solution was managed as RPMI 1640 medium (Gibo, USA), which was inclusive of 10% fetal bovine serum (FBS), 100 U/ml penicillin and 100 μg/ml streptomycin. Until cells grew to 90% confluence, they were prepared for following experiments.

Luciferase reporter gene assay
The UCA1 fragments that separately carried allele T and allele C, synthesized by Generay Biotech (Shanghai, China), were firstly digested by XhoI/NotI, and were then ligated to the psiCHECK-2 vector. Then the two vectors were, respectively, co-transfected with miR-873-5p mimic, miR-1207-5p mimic and miR-584 mimic (Ribobio, China) into CRC cells. Twenty-four hours later, luciferase activity of the cells was examined, as per the guidance of dual reporter gene detection kit (Promega, USA).

Western blotting
The CRC cells, having been washed with pre-cooled PBS for twice, were dissociated by protease inhibitorcontaining RIPA lysate. Then the concentration of extracted total protein was examined by referring to bicinchoninic acid (BCA) method. Exactly 25 μg protein of each sample was taken to carry out 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), after which the samples were transferred onto the polyvinylidene fluoride (PVDF) membrane. After 2-h blockage of the samples with 5% skimmed milk powder, mouse anti-human primary antibodies against E-cadherin, Ncadherin and Vimentin (1:1000, Abcam, USA) were supplemented to incubate samples at 4°C for overnight. With TBST applied to rinse proteins for 3 times, goat anti-mouse secondary antibodies labeled by horseradish peroxidase (HRP) (1:1000, Abcam, USA) were arranged to incubate proteins for another 1 h. Finally, development was accomplished via adoption of electrochemical luminescence, and gray values of the protein bands were analyzed by employing Image J software.

MTT assay
Precisely 3 × 10 3 cells of each treatment were, respectively, incubated for 0 h, 24 h, 48 h and 72 h, and then 50 μl MTT (5 g/L, Sigma, German) was prepared to treat the cells. After 4-h incubation at 37°C, the reactions were terminated by addition of 200 μl DMSO. At last, the absorbance (A) value of each sample was measured on the ELISA-Reader (Bio-Rad, USA) at the wavelength of 490 nm.

Colony formation assay
Exactly 3 × 10 3 cells per well were cultured in 5% CO 2 at 37°C for 9 days After being fixated with 95% methanol for 10 min, the cells were stained with 0.1% methylene blue. Finally, colonies that included > 50 cells were recorded under the microscope.

Cell apoptosis assay
Exactly 1 × 10 5 cells of each group were re-suspended by addition of 500 μl binding buffer, in strict accordance with the guidance of Annexin V-FITC/PI doublestaining kit (KeyGen Biotech, NanJing, China). Then 5 μl Annexin V-fluorescein isothiocyanate (FITC) and 5 μl propidium iodide (PI) were successively added to incubate cells. After incubation in the dark for 15 min, flow cytometry (model: FACS Calibur, BD, USA) was operated to determine apoptotic rate of the cells.

Transwell invasion assay
The polycarbonate membrane of transwell chambers was coated with a layer of Matrigel at 37°C for 45 min. Then CRC cells adjusted to a concentration of 4 × 10 5 /ml were added to the upper layer of Transwell chamber (Corning, USA), while 700 μl 10% FBS-containing medium was supplemented to the lower chamber. After incubating the cells in 5% CO 2 for 24 h, the cells in the lower chamber were fixated by pre-cooled 4% paraformaldehyde for 30 min and then stained by Giemsa for 20 min. A total of 8 views were randomly selected to calculate the average cell number.

Wound healing assay
When CRC cells inoculated into 6-well plates grew to nearly 80% confluence, a line was drawn on the back of cell plates utilizing the tip of a microsyringe (model: 10 μl). After on-going cultivation in 5% CO 2 for 24 h, the migratory distance of the cells was analyzed utilizing Image Pro-Plus 6.0 software (Media Cybernetics, USA).

Statistical analyses
SPSS software ver. 18 (Chicago, USA) was employed to perform all statistical analyses. The quantitative data were analyzed by virtue of independent two-tailed t test, and chi-square test was applied to compare categorical data. Logistic regression was implemented to evaluate the association of SNPs with CRC risk and prognosis. And Kaplan-Meier curves were devised to mirror survival of CRC patients, with log-rank test adopted to estimate between-group distinctions. Additionally, type 2 gene-environment interaction (GEI) model [34,35] and multifactor dimensionality reduction (MDR) model [36] were established, so as to assess the interactive contribution of SNPs in UCA1 and environmental factors to CRC risk. Statistical significance was identified in case of P < 0.05.

Results
Screening of SNPs that affected the secondary structure of lncRNA UCA1 According to the retrieval results from UCSC human GRCh37/hg19 database, UCA1, with a length of 7375 bp, was located within the segment of chr19:15,939,757-15, 947,131. Through usage of lncRNASNP database, a couple of SNPs (i.e. rs11085996 and rs12982687) were expected to affect binding of UCA1 to certain miRNAs. Subsequently, RNAfold software was employed to single out SNPs whose mutation was capable of altering the secondary structure of UCA1. It was indicated that the secondary structure of UCA1 differed when rs12982687 (allele C vs. allele T, Fig. 1a) or rs11085996 (allele G vs. allele A, Fig. 1b) carried distinct alleles.

Association of significant SNPs in UCA1 with CRC risk and prognosis
The CRC patients and healthy controls were well matched in terms of mean age (P = 0.140) and sex ratio (P = 0.069) ( Table 1). Nevertheless, respondents characterized by smoking/alcohol consumption, low intake of fruit/vegetable and uptake of hot-temperature food/beverage were associated with incremental susceptibility to CRC, as compared with those without such habits (all P < 0.05). Furthermore, homozygote GG of rs11085996 appeared as a protector against CRC in comparison to genotype AA/AG (OR = 0.84, 95% CI: 0.72-0.98, P = 0.026), yet barely any significant association was discoverable under the allelic or dominant model (both P > 0.05). Regarding rs12982687, its allele T seemed capable of reducing CRC risk as relative to allele C (OR = 0.78, 95% CI: 0.67-0.90, P = 0.001), and carriers of genotypes TC/TT were less susceptible to CRC than those carrying homozygote CC (OR = 0.78, 95% CI: 0.62-0.97, P = 0.025) ( Table 2). For another, CRC patients in possession of genotypes TT/CT of rs12982687 enjoyed more favorable prognosis than those with homozygote CC (P < 0.05), and homozygote GG of rs11085996 was associated with longer lifespan than homozygote AA (P = 0.008) (Fig. 1c). Considering an outstanding association of rs12982687 with CRC risk, this SNP was managed for later cellular experiments.
Stratified analyses regarding association of candidate SNPs in UCA1 with susceptibility to CRC (See Additional file 6) Evaluating the interactive contribution of candidate SNPs in UCA1 and smoking to CRC risk by establishing GEI and MDR models (See Additional file 6) The results of GEI model demonstrated that homozygote CC of rs12982687 could amplify the promoting effect of advanced age (≥ 50 years old) (γ = 2.80), smoking (γ = 1.98) and up-take of hot-temperature food/beverage (γ = 2.75) on CRC susceptibility (Table 3). With respect to rs11085996, smoking (γ > 1) and hottemperature food intake (γ > 1) both made AG/AA carriers more readily plagued by CRC, yet intake of fruit largely reduced the risk of developing CRC (γ = − 81.58) ( Table 4).
Mutation of rs12982687 disordered binding of UCA1 to certain CRC-specific miRNAs As predicted by LncRNASNP2 database, rs12982687C > T might disturb the binding capacity of UCA1 to some miRNAs (http://bioinfo.life.hust.edu.cn/lncRNASNP#!/ snp_info?snp=rs12982687). Among the miRNAs, miR-873-5p, miR-1207-5p, miR-584, miR-1910-3p and miR-940 The secondary structure of UCA1 differed between allele A (left) and allele G (right) of rs11085996. The red block was indicative of changes in the secondary structure of UCA1. c Four-year survival was compared among CRC patients who carried distinct genotypes of rs11085996 (up) and rs12982687 (down). d MiRNAs that possibly bound to UCA1 were quantified within CRC tissues and adjacent non-tumor tissues. The miRNAs were deemed as down-regulated when their -log10 (false discovery rate) was > 1.2 and log2 (fold change) was <− 0.5. e Correlations of UCA1 expression with miRNA levels were displayed by establishing correlation matrix. f Survival curves were depicted to assess the association of UCA1, miR-873-5p, miR-1207-5p and miR-584 expressions with CRC prognosis were evidently down-regulated within CRC tissues as opposed to adjacent normal tissues, with -log10(false discovery rate) > 1.2 and log2(fold change) < − 0.5 as the screening criterion (Fig. 1d). More than that, expressions of miR-873-5p (r s = − 0.36), miR-1207-5p (r s = − 0.21) and miR-584 (r s = − 0.20) were separately negatively correlated with UCA1 expression among the recruited CRC patients, as illustrated by the correlation matrix (Fig. 1e). Nevertheless, merely UCA1 and miR-873-5p were significant biomarkers for predicting the prognosis of CRC patients (both P < 0.05) (Fig. 1f).

Mutation of rs12982687 in UCA1 influenced EMT, proliferation and apoptosis of CRC cells
Transfection of miR-873-5p mimic was found to impair the trans-membrane potency of SW620 cell line (P < 0.05), rather than SW480 cell line (P > 0.05) (Fig. 3a). Correspondingly, SW620 cell line transfected by miR-873-5p inhibitor displayed stronger migratory potential than those transfected with plasmid vector (P < 0.05). Results of wound healing also presented that migration of SW620 cell line (P < 0.05), instead of SW480 cell line (P > 0.05), was weakened by miR-873-5p mimic (Fig. 3b). In stark contrast, SW620 cell line was prone to invade aggressively after being transfected by miR-873-5p inhibitor in comparison to the vector group (P < 0.05).  What's more, SW620 cell line (P < 0.05), but not SW480 cell line (P > 0.05), in the miR-873-5p mimic group expressed higher E-cadherin level and lower vimentin/ N-cadherin level than NC group (Fig. 3c). Conversely, managing the SW620 cell line with miR-873-5p inhibitor rendered them to produce low E-cadherin level and yet high vimentin/N-cadherin level (P < 0.05).
In addition, both proliferation and activity of SW620 cell line were hindered by miR-873-5p mimic (P < 0.05) ( Fig. 4a-b), and the apoptotic rate of SW620 cell line was also boosted under the influence of miR-873-5p mimic (P < 0.05) (Fig. 4c). Treatment of miR-873-5p inhibitor, on the contrary, boosted proliferation/activity and prohibited apoptosis of SW620 cell lines (P < 0.05) (Fig. 4). Nonetheless, the activity of SW480 cell line was nearly unaltered, regardless of proliferation, activity or apoptosis (all P > 0.05). Hence, it was implied that rs12982687 might indirectly regulate UCA1 expression by affecting the binding of UCA1 to miR-873-5p, which ultimately regulated the malignant phenotype of CRC cells.

Nicotine promoted malignancy of CRC cells and altered UCA1, miR-873 and HIF-1α expressions within CRC cells
Under treatment of 10 μmol/L nicotine, viability, migration and invasion of both SW620 and SW480 cell lines were elevated remarkably (P < 0.05), although the intensifying degree of SW620 cell line was more pronounced than SW480 cell line (P < 0.05) (Fig. 5). At the meantime, SW620 cell line was associated with lower E-cadherin expression and higher vimentin/N-cadherin expression than SW480 cell line (P < 0.05) (Fig. 6a). Concurrently, the apoptotic tendency of SW620 cell line was not as obvious as that of SW480 cell line (P < 0.05) (Fig. 6b). However, treatment of 50 μmol/L nicotine failed to significantly facilitate proliferation (Fig. 5a-b), invasion (Fig. 5c) and migration (Fig. 5d) of CRC cells, or prohibited their apoptosis (Fig. 6b) with statistical significance (P > 0.05). In addition, driven by 10 μmol/L nicotine, UCA1 expression was promoted, yet miR-873 expression was suppressed in both SW620 and SW480 cell lines (P < 0.05) (Fig. 6c). Further comparisons manifested that expressions of UCA1 and HIF-1α were altered less significantly in SW480 cell line than in SW620 cell line (P < 0.05).

Discussion
CRC onset was implied as the comprehensive outcome of genetic disorders and exposure to hazardous surroundings, and here we attempted to clarify the interactive role of SNPs in UCA1 and environmental risks in enhancing CRC risk. In the first place, the GEI model was established, which disclosed that genotypes TT/CT of rs12982687 could make smokers more vulnerable to CRC risk than homozygote CC (Table 3). Then as speculated by the RNAfold software, rs12982687 (C > T) seemed capable of changing the secondary structure of UCA1 (Fig. 1a), which could result in disordered binding of UCA1 to several miRNAs (http://bioinfo.life.hust.edu.cn/lncRN ASNP#!/snp_info?snp=rs12982687). Among the miR-NAs, it was intriguing to discover that miR-873-5p, miR-1207-5p and miR-873-5p were pronounced biomarkers in symbolizing CRC onset (Fig. 1d), and the miR-873-5p even displayed a tight linkage with poor survival of CRC patients (Fig. 1f). There appeared a possibility that mutation of rs12982687 was responsible for altering CRC susceptibility, due to its modifying binding of UCA1 with some miRNAs, especially miR-873-5p.
To date, a broad array of evidences have corroborated the inherent association of UCA1 with etiologies of diverse neoplasms, including bladder carcinoma, breast cancer, CRC, esophageal squamous cell carcinoma, gastric cancer, hepatocellular carcinoma, melanoma, ovarian cancer and squamous-cell carcinoma [16,21,[43][44][45][46][47][48]. Actually, UCA1 expression dwindled in adult tissues in comparison to embryonic tissues [49,50], which, from another angle, implied the adverse effect of highlyexpressed UCA1 on physiology of adults. Within this investigation, we strived to explain the CRC-promoting role of UCA1, from the perspective of its interaction with miRNAs [51]. In effect, the reciprocity of lncRNAs and miRNAs has been introduced in various disorders [25]. For instance, lncRNA ANRIL was reported to encourage proliferation of gastric carcinoma cells by silencing miR-99a and miR-449a [52], while miR-372 down-regulated expression of lncRNA HULC in hepatoma carcinoma cells by disaggregating transcriptional factors from the promoter of HULC [53]. Here it was suggested that homozygote CC of *The interaction coefficient (i.e. γ) was the result of βeg divided by βe. In case of γ > 1, genetic factors amplified the effect of environmental exposures; in case of γ > 1, genetic factors weakened the effect of environmental exposures; in case of γ = 1, no interaction was present between genetic factors and environmental exposures. When the environmental exposure was detrimental, genetic factors were deemed to be strongly protective in case of βe > 0 and γ < 0. The bold results indicated there is a statistically significance when its p value <0.05  The interaction coefficient (i.e. γ) was the result of βeg divided by βe. In case of γ > 1, genetic factors amplified the effect of environmental exposures; in case of γ > 1, genetic factors weakened the effect of environmental exposures; in case of γ = 1, no interaction was present between genetic factors and environmental exposures. When the environmental exposure was detrimental, genetic factors were deemed to be strongly protective in case of βe > 0 and γ < 0. The bold results indicated there is a statistically significance when its p value <0.05 Fig. 2 Single nucleotide polymorphism (SNP) rs12982687 (C > T) affected binding of UCA1 with miRNAs in colorectal cancer (CRC) cells. a The luciferase activity of CRC cells was compared between allele T + miR-873-5p/miR-1207-5p/miR-584 mimic group and allele C + miR-873-5p/miR-1207-5p/miR-584 mimic group. *: P < 0.05 when compared with allele T + miR-873-5p/miR-1207-5p/miR-584 mimic group. b UCA1 expression was determined when miR-873-5p mimic, miR-873-5p inhibitor or miR-NC was transfected into CRC cells, and miR-873-5p expression was measured among CRC cells of pcDNA3.1, pcDNA3.1-UCA1, si-UCA1 and NC groups. *: P < 0.05 when compared with NC group. c The genotype of rs12982687 (C > T) was determined within SW620 cell line and within SW480 cell line. d Expressional change of HIF-1α in SW620 and SW480 cell lines was determined in case of altered expression of UCA1 and miR-873-5p, and expressions of UCA1 and miR-873-5p were monitored when HIF-1α level was intentionally changed in CRC cells *: P < 0.05 when compared with NC group rs12982687 might prevent UCA1 expression from being suppressed by miR-873-5p (Fig. 2b), owing to which the anti-tumor role of miR-873-5p was hindered [54,55]. Intriguingly, our study demonstrated that homozygote CC of rs12982687 could facilitate proliferation, migration and invasion of CRC cells (Figs. 3 and 4), which might be supportive of this suggestion. Additionally, HIF-1 signaling seemed as a noteworthy pathway that mediated the role of UCA1 and miRNAs (i.e. miR-584, miR-1207-5p and miR-873-5p) in modifying CRC development, according to the results of KEGG enrichment analysis (Fig. 1d-e). There were convincing proofs that over-activation of HIF-1 could lead to aggressive proliferation and invasion of tumor cells, which was achieved via alteration of MMP9, VEGF and NDRG1 expressions [56][57][58]. In this investigation, we also observed that HIF-1α expression was boosted within CRC cells when UCA1 expression was up-regulated, and interestingly, a rise of UCA1 expression was detectable in case of over-expressed HIF-1α (Fig. 2d). There might be a feedback regulatory loop between UCA1 and HIF-1α, which was similar to the conclusion drawn by Li et al. [41,42], yet mechanisms in this regard demanded in-depth exploration. Single nucleotide polymorphism (SNP) rs12982687 (C > T) affected invasion, migration and epithelial-mesenchymal transition (EMT) of colorectal cancer (CRC) cells. a Invasive capability of SW620 and SW480 cell lines was assessed when miR-873-5p expression was up-regulated and down-regulated. *: P < 0.05 when compared with miR-NC group. b The migratory potency of SW620 and SW480 cell lines was compared among miR-NC, miR-873-5p mimic and miR-873-5p inhibitor groups. *: P < 0.05 when compared with miR-NC group. c Expressions of EMT-related proteins were detected within SW620 and SW480 cell lines of miR-NC, miR-873-5p mimic and miR-873-5p inhibitor groups. *: P < 0.05 when compared with miR-NC group Besides the above, we proposed that smoking, which created a hypoxic micro-environment that was quite common in solid tumors [59][60][61][62], might cooperate with homozygote CC of rs12982687 to engender a superimposed effect on CRC progression. On one hand, smoking was found to enlarge the likelihood of CRC risk among the population who carried rs12982687 (CC) (Figs. 5 and 6, Table 3). On the other hand, in-vitro experiments indicated that rs12982687 (CC) was associated with intensified proliferation, migration and invasion of Fig. 4 Single nucleotide polymorphism rs12982687 (C > T) affected viability, proliferation and apoptosis of colorectal cancer (CRC) cells. a Viability of SW620 and SW480 cell lines was measured among miR-NC, miR-873-5p mimic and miR-873-5p inhibitor groups. *: P < 0.05 when compared with miR-NC group. b Multiplication of SW620 and SW480 cell lines was contrasted among miR-NC, miR-873-5p mimic and miR-873-5p inhibitor groups. *: P < 0.05 when compared with miR-NC group. c The apoptotic rates of SW620 and SW480 cell lines were appraised among miR-NC, miR-873-5p mimic and miR-873-5p inhibitor groups. *: P < 0.05 when compared with miR-NC group CRC cells (Figs. 3 and 4). The UCA1 expression was significantly raised in SW620 and SW480 cell lines after supplementation of nicotine (Fig. 6c), although the expressional change of UCA1 was more profound in SW620 cell lines than in SW480 cell line. There was an implication that the genotype of rs12982687 was indeed correlative with the proliferative capability of CRC cells. Nonetheless, nicotine could promote proliferation of cancer cells through regulation of multifarious genetic pathways, and the expressional change of UCA1 might not reflect all carcinogenic actions of nicotine.
As a matter of fact, oxygen-poor was distinctly seen within CRC tissues [63]. For one thing, structure confusion of local tissues deformed tumor neovascularization and damaged normal blood supply of colon tissues. For another, extensive proliferation of tumor cells and expansion of tumor volume both aggravated hypoxia of CRC tissues. To survive the oxygen-deficient microenvironment, CRC cells generated a series of adaptive changes, such as variation of genomes [64]. Notably, HIF-1 signaling was among the genetic changes that happened during CRC exacerbation [65], and activation of HIF-1 signaling could powerfully enhance metastasis and drug resistance of CRC cells [66][67][68]. Moreover, HIF-1 signaling was also involved in the etiology of smokingtriggered cancer. In particular, knockdown of HIF-1α could interdict oncogenesis triggered by 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), a specific Nicotine regulated viability, proliferation, migration and invasion of colorectal cancer (CRC) cells. a Viability of SW620 and SW480 cell lines was assessed under treatments of 0, 5, 10 and 50 μmol/L nicotine. *: P < 0.05 when compared with treatment of 0 μmol/L nicotine. b Proliferation of SW620 and SW480 cell lines was measured under treatments of 0, 5, 10 and 50 μmol/L nicotine. *: P < 0.05 when compared with treatment of 0 μmol/L nicotine. c Invasive ability of SW620 and SW480 cell lines was evaluated under treatments of 0, 5, 10 and 50 μmol/L nicotine. *: P < 0.05 when compared with treatment of 0 μmol/L nicotine. d Nicotine at different concentrations (i.e. 0, 5, 10 and 50 μmol/L) facilitated migratory potential of SW620 and SW480 cell lines. *: P < 0.05 when compared with treatment of 0 μmol/L nicotine component of tobacco [37]. These evidences might, from another aspect, explain the interactive effect of UCA1 and smoking on CRC risk, owing to their sharing HIF-1 signaling.

Conclusion
SNP rs12982687 of UCA1 might influence the binding capacity of UCA1 with miR-873 and thereby affected function of HIF-1 signaling, which comprehensively altered risk and development of CRC. Moreover, smoking was also involved with CRC aggravation by modulating UCA1, miR-873 and HIF-1 signaling. Nonetheless, there were several flaws in the experimental design. Firstly, besides genotypes of rs12982687, the SW480 and SW620 cell lines utilized also differed in other genetic traits, so it might not be convincing enough to estimate the influence Fig. 6 Nicotine affected epithelial-mesenchymal transition (EMT) and apoptosis of colorectal cancer (CRC) cells. a Expressions of EMT-specific proteins were determined under treatments of 0, 5, 10 and 50 μmol/L nicotine. *: P < 0.05 when compared with treatment of 0 μmol/L nicotine. b Nicotine at different concentrations (i.e. 0, 5, 10 and 50 μmol/L) prohibited apoptosis of SW620 and SW480 cell lines. *: P < 0.05 when compared with treatment of 0 μmol/L nicotine. c UCA1 and miR-873-5p expressions were determined under the treatment of 10 μmol/L nicotine. *: P < 0.05 when compared with NC group of rs12982687 genotype on CRC risk with the couple of cell lines. It was highly recommended that the identical CRC cell line with two genotypes of rs12982687 could be employed, which was unfortunately unfinished here due to technical deficiencies. Secondly, how rs12982687 affected miR-873's regulation of UCA1 expression remained unclear, and if other mechanisms were involved in this modulation was still an uncertainty. Last but not the least, mice models were not established here to highlight the role of UCA1, miR-873 and HIF-1 signaling in CRC development. These deficiencies, after all, demanded further improvements.