Genetic Variants of HOTAIR Associated With Colorectal Cancer Susceptibility and Mortality

In colorectal carcinogenesis, the unique molecular and genetic changes that occur within cells result in specific CRC phenotypes. The involvement of the long non-coding RNA, HOTAIR, in cancer development, progression, and metastasis is well-established. Various studies have reported on the contribution of HOTAIR to cancer pathogenesis. Therefore, we selected four HOTAIR polymorphisms (rs7958904G>C, rs1899663G>T, rs4759314A>G, and rs920778T>C) to evaluate the association of each variant with CRC prevalence and prognosis. We conducted a case–control study of 850 individuals to identify the genotype frequencies of each polymorphism. The study population included 450 CRC patients and 400 control individuals that were randomly selected following a health screening. Notably, rs7958904 and rs1899663, their hetero genotype, and the dominant model were significantly different when compared to the healthy control group (rs7958904; AOR = 1.392, 95% CI = 1.052–1.843, P = 0.021). To evaluate the effect of HOTAIR polymorphisms on the survival rate, we analyzed patient mortality and relapse occurrence within 3 and 5 years with Cox-regression analysis. The rs7958904 CC polymorphism mortality rate was significantly higher than the GG polymorphism mortality rate (adjusted HR = 2.995, 95% CI = 1.189–7.542, P = 0.021). In addition, the rs920778 CC genotype was significantly different than the TT genotype (adjusted HR = 3.639, 95% CI = 1.435–9.230, P = 0.007). In addition, this study confirmed that genetic variants of HOTAIR alter the mRNA expression level (P < 0.01). We suggest that HOTAIR rs7958904G>C which is associated with CRC prevalence and mortality is a potential biomarker for CRC. The association between HOTAIR gene polymorphisms and CRC prevalence were reported for the first time.

Introduction

Colorectal cancer (CRC) is the third most common type of cancer and the second leading cause of cancer-related mortalities in Western countries (1, 2). The prognosis of patients with CRC depends on the tumor stage at the time of diagnosis. Alarmingly, over 57% of patients already exhibit regional or distant tumor cell spreading when the cancer is first diagnosed (1). The pathogenesis of CRC usually follows a stepwise progression that begins with a benign polyp and advances to invasive adenocarcinoma. In colorectal carcinogenesis, the unique molecular and genetic changes that occur within cells characterize the specific CRC phenotype. CRC phenotypes are associated with variable tumor behaviors that are relevant to disease prognosis and therapeutic efficacy. It was also reported that many various factors had an effect on CRC susceptibility and prognosis. Accordingly, an active area of research within the field is the assessment of biomarkers that can be used to predict disease prognosis or the response to therapy, which will contribute to individualized disease management.

The involvement of HOTAIR, a long non-coding RNA (lncRNA) in cancer development, progression, and metastasis is well-known. The oncogenic lncRNA is involved in the progression of multiple human cancers including breast, gastric, pancreatic, liver, hepatocellular, colon, lung, colorectal, and ovarian cancer. HOTAIR levels are also a predictive cancer biomarker, and HOTAIR expression is highly correlated with patient prognosis. HOTAIR promotes different processes including tumor growth, metastasis, invasion and migration, the epithelial-to mesenchymal-transition (EMT), and stemness via cancer-type specific pathways. These cancer phenotypes predominantly occur through HOTAIR-mediated epigenetic changes, which illustrates that lncRNA-guided mechanisms can be hijacked in the context of cancer (3–6).

Various studies have reported on the contribution of HOTAIR to cancer pathogenesis. A gene silencing mechanism that is activated by histone protein modification is caused by HOTAIR overexpression. HOTAIR is also directly involved in the translational regulation of target mRNA. It has also been suggested that HOTAIR functions as a molecular decoy in tumors where it sequesters several microRNAs (miRNAs) and RNA binding proteins (RBPs). For example, in esophageal and epithelial cancer, HOTAIR acts as a competing endogenous RNA to negatively regulate miR-148a. Suppression of miR-148a promotes the expression of Snail2, which is a key driver of the EMT in cancer, to enhance cell invasion and metastasis (7, 8). Interestingly, a subset of miRNAs, including miR-141, can regulate HOTAIR expression by targeting HOTAIR to the RNA-induced silencing complex for subsequent degradation by Ago2-induced cleavage (7, 9). In addition, HOTAIR can accelerate colon cancer development by down-regulating miRNA-34a (10). In renal carcinoma, HOTAIR can bind with the RBP human antigen R, which then directs HOTAIR to the Let7 miRNA-Ago2 complex and leads to microRNA-mediated suppression of HOTAIR through degradation (11, 12). Notably, proteins or miRNAs that bind to HOTAIR can repress its function, which adds another layer of complexity to the regulation of HOTAIR activity.

Several studies (13–15) report that HOTAIR is upregulated in esophageal squamous-cell carcinoma (ESCC), and that elevated HOTAIR expression is related to advanced TNM stage and poor histological differentiation. Molecular studies (16) revealed that HOTAIR combined with polycomb repressive complex 2 (PRC2) promoted histone H3K27 methylation of the define (WIF1) promoter, which reduced the WIF1 protein synthesis and expression. Consequently, the reduced ß-catenin degradation and the increased T-cell factor/lymphoid enhancer-binding factor levels activated the Wnt/β-catenin signaling pathway (17). This process eventually increased the expression of target genes and enhanced tumor cell proliferation, invasion, and metastasis. HOTAIR may also promote EMT. HOTAIR silencing increases E-cadherin expression and is accompanied by a decreased vimentin and differential matrix metalloproteinase (MMP) nine expression in colon cancer cells. Collectively, these downstream events suggest that HOTAIR may function as a novel pleiotropic regulator of the EMT (18). Furthermore, Xu et al. reported that HOTAIR silencing in gastric cancer cells downregulated MMP1 and MMP3 expression and inhibited invasion; whereas upregulation of E-cadherin and differential zonula occludens-1 expression reversed the EMT (7). HOTAIR can be influenced by miRNA regulation and EMT promotion in carcinogenesis, cancer development, and metastasis suggested by many previous studies. Factors that can regulate the expression of HOTAIR can provide new clues to cancer treatment and coping.

Recent studies demonstrate an association between HOTAIR polymorphisms and CRC (19, 20). Recent studies suggested that SNPs of HOTAIR (such as rs920778, rs4759314, rs1899663, rs12826786, rs874945, rs7958904, and rs10783618) acted as potential cancer susceptibility loci and were significantly associated with the increased risk of various cancers. However, there are no studies on the association of these five SNPs in Koreans. Therefore, we selected four HOTAIR polymorphisms (rs7958904G>C, rs1899663G>T, rs4759314A>G, and rs920778T>C) to evaluate the association between HOTAIR variants and CRC prevalence and prognosis. The minor allele frequencies of the selected polymorphisms are >5% in the Korean population.

Materials and Methods

Ethics Statement

All study protocols of participants were reviewed and approved by The Institutional Review Board of CHA Bundang Medical Center and followed the recommendations of the Declaration of Helsinki. Study subjects were recruited from the South Korean provinces of Seoul and Gyeonggi-do between 1996 and 2009. The Institutional Review Board of CHA Bundang Medical Center approved this genetic study in June 2009 (IRB No. 2009-08-077) and informed consent was obtained from study participants.

Study Population

We conducted a case–control study of 890 individuals. Four hundred and fifty patients that were diagnosed with CRC at the CHA Bundang Medical Center (Seongnam, South Korea) were enrolled from June 1996 to January 2009. The study only included CRC patients who had undergone surgical resection with a curative intent and who had histologically-proven adenocarcinoma. Collectively, there were 272 consecutive patients with colon cancer, 189 consecutive patients with rectal cancer, and 13 consecutive patients with unclassified CRC that underwent primary surgery. We retrospectively obtained patient information including the diagnosis date, the pathological stage, relapse events, and survival. Tumor staging of CRCs was performed according to the Sixth Edition of the American Joint Committee on Cancer (AJCC) staging manual. The control group was comprised of 416 individuals that were randomly selected following a health screening. This screening excluded patients with a history of thrombotic diseases or cancer. Patients with a high baseline blood pressure (systolic ≥140 mmHg or diastolic ≥90 mmHg) on more than one occasion or a history of antihypertensive medication were classified as having hypertension (HTN). Patients with a high fasting plasma glucose level (≥126 mg/dL), individuals who took oral hypoglycemic agents, or those with a history of insulin treatment were classified as having diabetes mellitus (DM). All study subjects were of Korean ethnicity and provided written informed consent. The study protocol was approved by the Institutional Review Board of the CHA Bundang Medical Center, Seongnam, South Korea.

Phenotype Measurements

Anthropometric measurements, including body mass index (BMI), were collected for each individual. Systolic and diastolic blood pressures were measured in the seated position after 10 min of rest. To measure physiological parameters, 3 mL of blood were collected after overnight fasting. Plasma glucose levels were measured in duplicate by adapting the hexokinase method for use with an automated analyzer (TBA 200FR NEO, Toshiba Medical Systems, Tokyo, Japan). High density lipoprotein-cholesterol (HDL-C) levels were measured with enzymatic colorimetric methods using commercial reagent sets (TBA 200FR NEO, Toshiba Medical Systems). The homocysteine (Hcy) plasma concentration was measured with a fluorescent polarizing immunoassay (FPIA) with IMx (Abbott Laboratories, Chicago, IL, USA). The folate (FA) plasma concentration was measured with a radioimmunoassay kit (ACS:180; Bayer, Tarrytown, NY, USA).

Genotyping

DNA was extracted from leukocytes using a G-DEX II Genomic DNA Extraction kit (Intron Biotechnology, Seongnam, Korea) according to the manufacturer's directions. To analyze HOTAIR genotypes, polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and TaqMan allele discrimination analysis were chosen because these procedures were more economical and cost-effective when compared to entire gene sequencing. Detailed conditions for PCR-RFLP, real-time PCR, and qRT-PCR methods are presented in Supplementary Materials. To validate the RFLP findings for each polymorphism, 30% of the PCR assays were randomly selected and duplicated, and followed by DNA sequencing. Sequencing was performed with an ABI 3730xl DNA Analyzer (Applied Biosystems, Foster City, CA, USA). The concordance of the quality control samples was 100%.

Statistical Analysis

To analyze baseline characteristics of CRC and control patients, we used chi-square tests for categorical data and Student's t-tests for continuous data (21, 22). We estimated the association of VEGF and KDR polymorphisms with colorectal cancer incidence using adjusted odds ratios (AORs) and 95% confidence intervals (95% CIs) with multivariate logistic regression that was adjusted for age, gender, HTN, DM, BMI, and HDL-C. We chose HTN, DM, BMI, and HDL-C as adjustment variables because the risk factors for metabolic syndrome are closely associated with colorectal cancer. Cox-regression models were used to analyze the independent prognostic importance of various markers, and the results, which excluded 100 CRC patients that had insufficient medical history, were adjusted for age, gender, tumor differentiation, tumor site, chemotherapy, and cancer stage. Overall survival (OS) was defined as the time period from surgery until death or the last follow-up. Relapse-free survival (RFS) was defined as the time period from surgery until cancer reoccurred or the last follow-up. Hazard ratios (HRs) are presented with a 95% CI. Participants were followed for a median of 34 months (range, 4–173 months). The estimated 3-years OS and RFS rates for all patients were 82.6 and 81.7%, respectively. Analyses were performed using GraphPad Prism 4.0 (GraphPad Software Inc., San Diego, CA, USA) and Medcalc Version 12.7.1.0 (Medcalc Software, Mariakerke, Belgium). Haplotypes for multiple loci were estimated using the expectation-maximization algorithm with SNPAlyze (Version 5.1; DYNACOM Co, Ltd, Yokohama, Japan).

Results

We evaluated HOTAIR gene polymorphisms (rs7958904, rs1899663, rs4759314, and rs920778) in CRC patients. We first classified the characteristics of patients in each polymorphism group, which included dividing the CRC patients into two groups based on the cancer location. Healthy control and colorectal patients were also age and gender matched. The detail information for sex ratios between the control and CRC patients were 41.6 and 45.1%, respectively, and the mean age was 61.2 and 62.1 years. There was no statistical significance between the two groups. The hypertension proportion in participants was higher, 40.4%, in the control group, compared with 33.1% in CRC patients (P = 0.030). Also, the total subject ratio of patients with diabetes mellitus was 13.2% in the control group and 16.0% of CRC patients showed a similar frequency between the two groups. In comparing the mean levels of plasma folate and plasma homocysteine, folate was found to be 9.0 in the control group and 7.8 in the patient group (P < 0.0001). However, homocysteine did not show any difference between the two groups (Table 1).

TABLE 1

Table 1. Baseline characteristics between CRC patients and healthy control subjects.

In Table 2, we identified the genotype frequencies of each polymorphism. The rs7958904 and rs1899663 polymorphisms, their hetero genotypes, and dominant models were significantly different when compared to the healthy control group (HOTAIR rs7958904 GG vs. GC: AOR = 1.352, 95% CI = 1.014–1.801, P = 0.040; GG vs. GC+CC: AOR = 1.351, 95% CI = 1.027–1.777, P = 0.032; HOTAIR rs1899663 GG vs. GT: AOR = 1.3385, 95% CI = 1.004–1.784, P = 0.047; GG vs. GT+TT: AOR = 1.378, 95% CI = 1.043–1.822, P = 0.024). However, the rs4759314 and rs920778 polymorphisms were not statistically different when compared with the control group. Additionally, rectal cancer patients (189 samples) exhibited a significant difference with respect to both the rs7958904 GC type polymorphism (AOR = 1.559, 95% CI = 1.076–2.258, P = 0.019) and the dominant model (AOR = 1.547, 95% CI = 1.085–2.205, P = 0.016) when compared with the control patients. Although the rs1899663 GT type polymorphism was not statistically significant, the dominant model was statistically different when compared to control patients (AOR = 1.481, 95% CI = 1.035–2.120, P = 0.032). Interestingly, patients with localized colon cancer did not exhibit any differences when compared with healthy controls (P > 0.05).

TABLE 2

Table 2. Comparison of HOTAIR polymorphism genotype frequencies between CRC patients and healthy controls.

To evaluate the effect HOTAIR polymorphisms on survival rate, we analyzed mortality and 3-years relapse occurrence with Cox-regression analysis. The HOTAIR rs7958904 CC genotype mortality was significantly higher than the GG genotype (adjusted HR = 2.995, 95 %CI = 1.189–7.542, P = 0.021). Additionally, the HOTAIR rs920778 CC genotype was significantly different when compared to the TT genotype (adjusted HR = 3.639, 95% CI = 1.435–9.230, P = 0.007; Table 3 and Figure 1). To complete the subgroup analysis, we analyzed the mortality of colon cancer patients, and found that patients with HOTAIR rs7958904 and rs920778 polymorphisms had greater survival rates. Interestingly, the HOTAIR rs1899663 TT genotype was associated with increased colon cancer mortality (adjusted HR = 4.507, 95% CI = 1.167–17.413, P = 0.030), and the association with increased mortality only held for colon cancer (Table 3).

TABLE 3

Table 3. HOTAIR polymorphism genotype frequencies and CRC patient 3-years mortality in overall, colon, and rectum cancer.

FIGURE 1

Figure 1. Survival plot from a Cox proportiona hazards model with HOTAIR rs7958904G>C and rs920778T>C polymorphisms in colorectal cancer (CRC). Survival curve of patients grouped by 3 years mortality in colorectal cancer based on (A) HORAIR rs7958904GG vs. HOTAIR rs7958904CC genotypes (B) HORAIR rs920778TT vs. HOTAIR rs920778CC genotypes. HR, hazard ratio.

We evaluated synergic effect between clinical parameters and HOTAIR polymorphisms (Supplementary Table 1). In order to confirm the synergistic effects of clinical factors and genotypes, we performed stratified analysis based on the classification of sex, age, hypertension, diabetes mellitus, plasma folate, and homocysteine levels analyzed the interactions with HOTAIR polymorphisms. Especially, when analyzing the synergistic effect of increasing the risk of CRC with genotypes, according to the concentration of folate in plasma, it was found that the association of CRC risk was increased at low folate concentrations (≤ 3.08 nmol/L) group in HOTAIR rs7958904 GC+CC and rs1899663 GT+TT types (Figure 2).

FIGURE 2

Figure 2. The effects of HOTAIR rs7958904 G>C and rs1899663 variant on colorectal cancer (CRC) risk modulated by clinical factors. (A) The synergic effect for CRC susceptibility combined between HOTAIR rs7958904 GC+CC and low plasma folate group (AOR, 3.191) than rs7958904 GC+CC with high plasma folate group (AOR, 1.285) or rs7958904 GG with low plasma folate group (AOR, 2.545). (B) The HOTAIR rs1899663 G>T were show to elevated CRC prevalence, in the case of HOTAIR rs1988663 GT+TT with low plasma folate group (AOR, 3.196) than rs1988663 GT+TT with high plasma folate group (AOR, 1.353) or rs1899663 GG with low plasma folate group (AOR, 2.601).

We measured HOTAIR expression levels in colon cancer and normal tissues to determine the expression pattern of specific HOTAIR polymorphisms. The overall expression of HOTAIR was 2-fold greater in CRC when compared to control tissues (Supplementary Figure 1, P < 0.01). Next, to identify the altered expression patterns of the specific HOTAIR genotype, rs7958904G>C, we analyzed the expression of the HOTAIR rs7958904G>C polymorphism compared with normal tissue and tumor tissue. The HOTAIR rs7958904 CC genotype was highly expressed in tumor tissue than normal tissues (P < 0.01). However, GG and C genotypes do not have significantly different expression levels compared with normal tissue and tumor tissue (Supplementary Figure 2).

LD block and haplotype analyses were performed to identify linkage disequilibrium between SNPs. When analyzed using Haploview software (https://www.broadinstitute.org/haploview/haploview), it was confirmed that the R² value was the highest at 0.810 between rs7958904 and rs1899663 in the analysis between four SNPs, but no LD block was generated (Supplementary Figure 3). The haplotype analysis (Supplementary Tables 2–4) revealed that the various haplotypes were at a higher risk. The data suggests that the C-T (HOTAIR rs7958904G>C/rs1899663G>T) haplotype is more prevalent in CRC when compared to the baseline (OR = 1.584, 95% CI = 1.117–2.247, P = 0.010). The C-G-A (HOTAIR rs7958904G>C/rs1899663G>T/rs4759314A>G) haplotype analysis revealed that the C-G-A haplotype was specifically associated with reduced CRC occurrence (OR = 0.212, 95%CI = 0.060–0.750, P = 0.008). Furthermore, we identified genotype combinations (Supplementary Tables 5, 6) where the HOTAIR rs1899663 GT genotype, when combined with other polymorphisms, had a synergistic effect that increased CRC risk. However, several genotype combinations demonstrated that the occurrence of cancer was location-dependent. The GC+GT (HOTAIR rs7598904+rs1899663) and AG+TT (HOTAIR rs4759314+rs920778) genotypes were significantly different in rectal cancer patients, but there was not any statistical significance in colon cancer patients. Supplementary Table 5 summarizes the statistical powers of positive genetic associations in the present study.

Discussion

In this study, we analyzed the association between the risk development of four SNPs of HOTAIR in Korean CRC patients and the correlation between CRC mortality by survival analysis. In addition, we confirmed that the HOTAIR expression for rs7958904G> C showed the strongest statistical correlation. Interestingly, our results confirm the association between rs7958904G> C and rs1899663G> T rather than the previously known rs4759314A> G that is related to CRC susceptibility. It was confirmed that similar correlations were found in the colon and rectum regions.

The lncRNAs participate in gene regulation (23, 24). Numerous lncRNAs have been shown to promote cell invasion and metastasis (25, 26). HOTAIR is an lncRNA located in the HOXC locus, and it can interact with PRC2. This interaction mediates the methylation of lysine 27 and the demethylation of lysine 4 on histone H3 in the HOXD locus, where enhancer of zeste homolog 2 also plays a regulatory role (27, 28). HOTAIR can alter the state of chromosomes, which affects gene expression. Importantly, HOTAIR expression was upregulated in cancer tissue samples from patients with breast, pancreatic, liver, gastric, and non-small cell lung cancers, and elevated HOTAIR expression was even more pronounced in metastatic tissue. Both in vivo and in vitro studies confirm that upregulated expression of HOTAIR enhances tumors invasion and metastasis (28–30). For example, prior studies report that HOTAIR is significantly upregulated and closely related to invasion and metastasis in endothelial cells (15, 31). Studies with CRC patients revealed that HOTAIR expression levels are higher in CRC tissue when compared with corresponding non-cancerous tissue. Similar findings were reported for ovarian (32) and laryngeal squamous-cell cancers. High expression levels of HOTAIR are correlated with the presence of liver metastases, and CRC patients with elevated HOTAIR expression have a worse prognosis than patients with tumors exhibiting low HOTAIR expression (33). Based on an in-situ hybridization assay, 91 of 160 (56.87%) paraffin embedded nasopharyngeal carcinoma (NPC) biopsy specimens had elevated HOTAIR levels. NPC patients with high tumor HOTAIR expression levels had a worse OS prognosis than patients with low HOTAIR expression (34). Furthermore, elevated HOTAIR expression was detected in patients that had high histological grade tumors or advanced clinical stage cancer (35). Depletion of HOTAIR by short interfering RNA, decreased invasion and increased the apoptosis of human epithelial type 2 cells. In addition, tumor growth was significantly inhibited in mice injected with HOTAIR-deficient cells (35). Overexpression of HOTAIR was also strongly associated with high-grade tumors and metastasis in gastrointestinal stromal tumor specimens (36). ESCC patients with elevated levels of HOTAIR had significantly lower 5-years survival rates than HOTAIR-negative patients (15). Finally, high HOTAIR levels in primary sarcoma correlated with a high metastasis probability (37).

Recent studies identified that lncRNAs bind to miRNAs to “communicate” with other RNA targets (38). Interestingly, lncRNAs and miRNA exhibit reciprocal regulation. As many lncRNAs control important physiological functions, the abundance and binding of each miRNA and lncRNA directly alters cell function. By sharing common miRNA binding sites with mRNA targets, lncRNAs sequester and compete with miRNAs to inhibit miRNA functionality and alleviate mRNA repression (39). Interestingly, 40% of miRNAs are found in the introns of protein coding genes (40). Furthermore, analyses indicate that 10% of lncRNA genes also host an miRNA, either in an intron or exon (41). In addition, miRNAs bind to lncRNAs, and when combined with other RBPs, regulate lncRNA stability and miRNA-mediated decay. Recently, it was reported that RNA interference-mediated knockdown of HOTAIR altered the expression of HOTAIR target genes and suppressed the invasion of gastrointestinal stromal tumor cells. Another study revealed that HOTAIR functions as a “miRNA sponge” that silenced miRNAs (tumor suppressor), and thereby induced the overexpression of oncogenic genes (8).

At this point, HOTAIR expression has been well-studied. The mechanisms of cancer onset and metastasis, and HOTAIR SNP have also actively been studied. However, there are insufficient studies to establish the link between HOTAIR polymorphism and CRC, and some target SNPs of the previously reported HOTAIR polymorphism are limited. In particular, according to genotypes that vary by SNP, there are no existing reported studies on the changes in HOTAIR expression extending from the linkage study between SNP and CRC. Functional studies through studies on statistically related polymorphisms may be a way to confirm the contribution of HOTAIR regulation in cancer onset/metastasis and maybe a major guideline for cancer onset and prognosis. In our results, the C allele of rs7958904G>C was found to be the association for CRC susceptibility, mortality, and altered expression level. In previous studies, rs7958904 is a SNP present in 3′UTR of HOTAIR and it has been reported that SNPs located in 3′UTR. Regulation mechanisms of gene expression including mRNA stability, miRNA binding, cis-regulation. However, in some cases, the results of the analysis regarding the association of SNPs is not clear. In our results, there was not an LD block between the four SNPs we selected in the haplotype analysis, and the allele combination analysis showed that the CRC risk increased in some cases, even with the combination of minor alleles had shown to reduce the risk. This result is considered to be a variation that could have been caused by classifying the haplotype with a small sample size- Performing statistical progress based on these results, it is necessary to collect additional samples and analyze the population with more people.

CRC is associated with a variety of factors. Representative factors include hypertension (42, 43), folate (44, 45), homocysteine (46), and diabetes milieus (47). Thus, treatment strategies and prognostic management methods change. Therefore, it is well-known to check the association with clinical parameters, and it is necessary to perform disease management based on this. In addition, due to the development of many studies and techniques, the pathogenesis and diagnosis methods of CRC are developing, and the rate of CRC treatment by early diagnosis is increasing (25, 26, 48, 49). However, despite advances in technology, people in blind spots, those with the burden of medical expenses and the difficulties of living conditions do not receive these benefits (50–52). Therefore, there is an urgent need to find new therapeutic targets. In this regard, it is necessary to study markers for early diagnosis and treatment, and the HOTAIR polymorphism proposed in our study is expected to satisfy this need.

There are some limitations to our research. First, case-control studies were conducted only in the Korean population. However, the frequency of genotypes was consistent with studies in other populations previously reported and confirmed the frequency of meeting the HWE p-value. Although the size of our research sample was small, statistical verification was confirmed through statistical power analysis. Secondly, the duration of survival analysis in cancer research is generally 5 years, but our analysis is 3 years. This was done because when the data was compiled for a 5-years survival analysis, the explanatory power of the analysis decreased because there were too many samples that lost the survival information. Finally, although our study suggested a link between the HOTAIR SNP risk of CRC and survival analysis, it confirmed differences in HOTAIR expression for the genotype of rs7958904, RNA second structure changes, RNA stability regulation, and RNA of these SNPs. The mechanism of expression regulation such as the role of binding protein sites, was not elucidated.

In conclusion, this study analyzed the expression levels of HOTAIR and miRNAs in genetic variants. In addition, the association between HOTAIR gene polymorphisms and colorectal cancer prevalence, including the altered expression patterns of HOTAIR polymorphisms, were reported for the first time.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Ethics Statement

All study protocols of participants were reviewed and approved by The Institutional Review Board of HA Bundang Medical Center and followed the recommendations of the Declaration of Helsinki. Study subjects were recruited from the South Korean provinces of Seoul and Gyeonggi-do between 1996 and 2009. The Institutional Review Board of CHA Bundang Medical Center approved this genetic study in June 2009 (IRB No. 2009-08-077) and informed consent was obtained from study participants.

Author Contributions

NK and JWK: onceptualization. JOK, EK, and JL: methodology. HP, CR, SK, and DO: validation. JOK, HJ, JL, HP, and CR: formal analysis. EK and SK: investigation. JWK and DO: resources. NK: data curation. JOK, HJ, and EK: writing—original draft preparation. JOK: writing—review and editing. NK and JWK: project administration.

Funding

This research was funded by a National Research Foundation of Korea (NRF) Grant (2018R1D1A1B07047604), funded by the Korean Government and was supported by a grant of the Korea Healthcare technology R&D project, Ministry for Health, Welfare & Family Affairs (HI18C19990200). The funding bodies were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

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.

Supplementary Material

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

References

1. Figueredo A, Coombes ME, Mukherjee S. Adjuvant therapy for completely resected stage II colon cancer. Cochrane Database Syst Rev. (2008) 16:CD005390. doi: 10.1002/14651858.CD005390.pub2

CrossRef Full Text | Google Scholar

2. Jemal A, Siegel R, Ward E, Hao Y, Xu J, Murray T, et al. Cancer statistics, 2008. CA Cancer J Clin. (2008) 58:71–96. doi: 10.3322/CA.2007.0010

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Loewen G, Zhuo Y, Zhuang Y, Jayawickramarajah J, Shan B. lincRNA HOTAIR as a novel promoter of cancer progression. J Can Res Updates. (2014) 3:134–40. doi: 10.6000/1929-2279.2014.03.03.3

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Zhang J, Zhang P, Wang L, Piao HL, Ma L. Long non-coding RNA HOTAIR in carcinogenesis and metastasis. Acta Biochim Biophys Sin. (2014) 46:1–5. doi: 10.1093/abbs/gmt117

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Hajjari M, Salavaty A. HOTAIR: an oncogenic long non-coding RNA in different cancers. Cancer Biol Med. (2015) 12:1–9. doi: 10.7497/j.issn.2095-3941.2015.0006

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Schmitt AM, Chang HY. Long noncoding RNAs in cancer pathways. Cancer Cell. (2016) 29:452–63. doi: 10.1016/j.ccell.2016.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Xu ZY, Yu QM, Du YA, Yang LT, Dong RZ, Huang L, et al. Knockdown of long non-coding RNA HOTAIR suppresses tumor invasion and reverses epithelial-mesenchymal transition in gastric cancer. Int J Biol Sci. (2013) 9:587–97. doi: 10.7150/ijbs.6339

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Xu F, Zhang J. Long non-coding RNA HOTAIR functions as miRNA sponge to promote the epithelial to mesenchymal transition in esophageal cancer. Biomed Pharmacother. (2017) 90:888–96. doi: 10.1016/j.biopha.2017.03.103

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Padua Alves C, Fonseca AS, Muys BR, de Barros ELBR, Burger MC, de Souza JE, et al. Brief report: The lincRNA Hotair is required for epithelial-to-mesenchymal transition and stemness maintenance of cancer cell lines. Stem Cells. (2013) 31:2827–32. doi: 10.1002/stem.1547

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Peng CL, Zhao XJ, Wei CC, Wu JW. LncRNA HOTAIR promotes colon cancer development by down-regulating miRNA-34a. Eur Rev Med Pharmacol Sci. (2019) 23:5752–61. doi: 10.26355/eurrev_201907_18312

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Yoon JH, Abdelmohsen K, Kim J, Yang X, Martindale JL, Tominaga-Yamanaka K, et al. Scaffold function of long non-coding RNA HOTAIR in protein ubiquitination. Nat Commun. (2013) 4:2939. doi: 10.1038/ncomms3939

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Loewen G, Jayawickramarajah J, Zhuo Y, Shan B. Functions of lncRNA HOTAIR in lung cancer. J Hematol Oncol. (2014) 7:90. doi: 10.1186/s13045-014-0090-4

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chen FJ, Sun M, Li SQ, Wu QQ, Ji L, Liu ZL, et al. Upregulation of the long non-coding RNA HOTAIR promotes esophageal squamous cell carcinoma metastasis and poor prognosis. Mol Carcinog. (2013) 52:908–15. doi: 10.1002/mc.21944

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Li X, Wu Z, Mei Q, Li X, Guo M, Fu X, et al. Long non-coding RNA HOTAIR, a driver of malignancy, predicts negative prognosis and exhibits oncogenic activity in oesophageal squamous cell carcinoma. Br J Cancer. (2013) 109:2266–78. doi: 10.1038/bjc.2013.548

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lv XB, Lian GY, Wang HR, Song E, Yao H, Wang MH. Long noncoding RNA HOTAIR is a prognostic marker for esophageal squamous cell carcinoma progression and survival. PLoS ONE. (2013) 8:e63516. doi: 10.1371/journal.pone.0063516

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Ge XS, Ma HJ, Zheng XH, Ruan HL, Liao XY, Xue WQ, et al. HOTAIR, a prognostic factor in esophageal squamous cell carcinoma, inhibits WIF-1 expression and activates Wnt pathway. Cancer Sci. (2013) 104:1675–82. doi: 10.1111/cas.12296

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Salmerón-Bárcenas EG, Illades-Aguiar B, Del Moral-Hernández O, Ortega-Soto A, Hernández-Sotelo D. HOTAIR knockdown decreased the activity Wnt/β-catenin signaling pathway and increased the mRNA levels of its negative regulators in Hela Cells. Cell Physiol Biochem. (2019) 53:948–60. doi: 10.33594/000000188

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wu ZH, Wang XL, Tang HM, Jiang T, Chen J, Lu S, et al. Long non-coding RNA HOTAIR is a powerful predictor of metastasis and poor prognosis and is associated with epithelial-mesenchymal transition in colon cancer. Oncol Rep. (2014) 32:395–402. doi: 10.3892/or.2014.3186

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Xue Y, Gu D, Ma G, Zhu L, Hua Q, Chu H, et al. Genetic variants in lncRNA HOTAIR are associated with risk of colorectal cancer. Mutagenesis. (2015) 30:303–10. doi: 10.1093/mutage/geu076

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Zhang J, Liu X, You LH, Zhou RZ. Significant association between long non-coding RNA HOTAIR polymorphisms and cancer susceptibility: a meta-analysis. Onco Targets Ther. (2016) 9:3335–43. doi: 10.2147/OTT.S107190

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Kim TN, Kim WT, Jeong MS, Mun MH, Kim MH, Lee JZ, et al. Short rare minisatellite variant of BORIS-MS2 is related to bladder cancer susceptibility. Genes Genomics. (2019) 41:249–56. doi: 10.1007/s13258-018-0771-4

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ryu CS, Sakong JH, Ahn EH, Kim JO, Ko D, Kim JH, et al. Association study of the three functional polymorphisms (TAS2R46G>A, OR4C16G>A, and OR4X1A>T) with recurrent pregnancy loss. Genes Genomics. (2019) 41:61–70. doi: 10.1007/s13258-018-0738-5

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Khalil AM, Guttman M, Huarte M, Garber M, Raj A, Rivea Morales D, et al. Many human large intergenic noncoding RNAs associate with chromatin-modifying complexes and affect gene expression. Proc Natl Acad Sci USA. (2009) 106:11667–72. doi: 10.1073/pnas.0904715106

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ponting CP, Oliver PL, Reik W. Evolution and functions of long noncoding RNAs. Cell. (2009) 136:629–41. doi: 10.1016/j.cell.2009.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Lloyd JM, McIver CM, Stephenson SA, Hewett PJ, Rieger N, Hardingham JE. Identification of early-stage colorectal cancer patients at risk of relapse post-resection by immunobead reverse transcription-PCR analysis of peritoneal lavage fluid for malignant cells. Clin Cancer Res. (2006) 12:417–23. doi: 10.1158/1078-0432.CCR-05-1473

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Lien GS, Lin CH, Yang YL, Wu MS, Chen BC. Ghrelin induces colon cancer cell proliferation through the GHS-R, Ras, PI3K, Akt, and mTOR signaling pathways. Eur J Pharmacol. (2016) 776:124–31. doi: 10.1016/j.ejphar.2016.02.044

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. (2010) 464:1071–6. doi: 10.1038/nature08975

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Tsai MC, Manor O, Wan Y, Mosammaparast N, Wang JK, Lan F, et al. Long noncoding RNA as modular scaffold of histone modification complexes. Science. (2010) 329:689–93. doi: 10.1126/science.1192002

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Ishibashi M, Kogo R, Shibata K, Sawada G, Takahashi Y, Kurashige J, et al. Clinical significance of the expression of long non-coding RNA HOTAIR in primary hepatocellular carcinoma. Oncol Rep. (2013) 29:946–50. doi: 10.3892/or.2012.2219

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kim K, Jutooru I, Chadalapaka G, Johnson G, Frank J, Burghardt R, et al. HOTAIR is a negative prognostic factor and exhibits pro-oncogenic activity in pancreatic cancer. Oncogene. (2013) 32:1616–25. doi: 10.1038/onc.2012.193

CrossRef Full Text | Google Scholar

31. Cao W, Shi F, Wu L, Yang K, Tian F, Chen G, et al. [Genome wide screening and characterization of long non-coding RNAs in esophageal cancer]. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. (2014) 31:587–90. doi: 10.3760/cma.j.issn.1003-9406.2014.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Cui L, Xie XY, Wang H, Chen XL, Liu SL, Hu LN. [Expression of long non-coding RNA HOTAIR mRNA in ovarian cancer]. Sichuan Da Xue Xue Bao Yi Xue Ban. (2013) 44:57–9.

PubMed Abstract | Google Scholar

33. Kogo R, Shimamura T, Mimori K, Kawahara K, Imoto S, Sudo T, et al. Long noncoding RNA HOTAIR regulates polycomb-dependent chromatin modification and is associated with poor prognosis in colorectal cancers. Cancer Res. (2011) 71:6320–6. doi: 10.1158/0008-5472.CAN-11-1021

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Nie Y, Liu X, Qu S, Song E, Zou H, Gong C. Long non-coding RNA HOTAIR is an independent prognostic marker for nasopharyngeal carcinoma progression and survival. Cancer Sci. (2013) 104:458–64. doi: 10.1111/cas.12092

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Li D, Feng J, Wu T, Wang Y, Sun Y, Ren J, et al. Long intergenic noncoding RNA HOTAIR is overexpressed and regulates PTEN methylation in laryngeal squamous cell carcinoma. Am J Pathol. (2013) 182:64–70. doi: 10.1016/j.ajpath.2012.08.042

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Niinuma T, Suzuki H, Nojima M, Nosho K, Yamamoto H, Takamaru H, et al. Upregulation of miR-196a and HOTAIR drive malignant character in gastrointestinal stromal tumors. Cancer Res. (2012) 72:1126–36. doi: 10.1158/0008-5472.CAN-11-1803

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Milhem MM, Knutson T, Yang S, Zhu D, Wang X, Leslie KK, et al. Correlation of MTDH/AEG-1 and HOTAIR expression with metastasis and response to treatment in sarcoma patients. J Cancer Sci Ther S. (2011) 5:4. doi: 10.4172/1948-5956.S5-004

CrossRef Full Text | Google Scholar

38. Ebert MS, Sharp PA. Emerging roles for natural microRNA sponges. Curr Biol. (2010) 20:R858–61. doi: 10.1016/j.cub.2010.08.052

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Xia T, Liao Q, Jiang X, Shao Y, Xiao B, Xi Y, et al. Long noncoding RNA associated-competing endogenous RNAs in gastric cancer. Sci Rep. (2014) 4:6088. doi: 10.1038/srep06088

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Carthew RW, Sontheimer EJ. Origins and mechanisms of miRNAs and siRNAs. Cell. (2009) 136:642–55. doi: 10.1016/j.cell.2009.01.035

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Xiong W, Qu Y, Chen H, Qian J. Insight into long noncoding RNA-miRNA-mRNA axes in myocardial ischemia-reperfusion injury: the implications for mechanism and therapy. Epigenomics. (2019) 11:1733–48. doi: 10.2217/epi-2019-0119

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Sud S, O'Callaghan C, Jonker C, Karapetis C, Price T, Tebbutt N, et al. Hypertension as a predictor of advanced colorectal cancer outcome and cetuximab treatment response. Curr Oncol. (2018) 25:e516–26. doi: 10.3747/co.25.4069

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Wang SC, Schulman-Marcus J, Fantauzzi J, Bevington T, Sayegh A, Lee E, et al. Colon cancer laterality is associated with atherosclerosis and coronary artery disease. J Gastrointest Oncol. (2019) 10:30–6. doi: 10.21037/jgo.2018.09.18

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Catala GN, Bestwick CS, Russell WR, Tortora K, Giovannelli L, Moyer MP, et al. Folate, genomic stability and colon cancer: the use of single cell gel electrophoresis in assessing the impact of folate in vitro, in vivo and in human biomonitoring. Mutat Res. (2019) 843:73–80. doi: 10.1016/j.mrgentox.2018.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Agarwal S, Behring M, Hale K, Al Diffalha S, Wang K, Manne U, et al. MTHFD1L, a folate cycle enzyme, is involved in progression of colorectal cancer. Transl Oncol. (2019) 12:1461–7. doi: 10.1016/j.tranon.2019.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Weinstein SJ, Albanes D, Selhub J, Graubard B, Lim U, Taylor PR, et al. One-carbon metabolism biomarkers and risk of colon and rectal cancers. Cancer Epidemiol Biomarkers Prev. (2008) 17:3233–40. doi: 10.1158/1055-9965.EPI-08-0459

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Soltani G, Poursheikhani A, Yassi M, Hayatbakhsh A, Kerachian M, Kerachian MA, et al. Obesity, diabetes and the risk of colorectal adenoma and cancer. BMC Endocr Disord. (2019) 19:113. doi: 10.1186/s12902-019-0444-6

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Shahbazian H, Nasuri Y, Hosseini SM, Arvandi S, Razzaghi S. A report of the frequency of colorectal carcinoma and involved lymph nodes in South-West Iran. Indian J Med Paediatr Oncol. (2016) 37:38–41. doi: 10.4103/0971-5851.177014

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Zaenker P, Ziman MR. Serologic autoantibodies as diagnostic cancer biomarkers–a review. Cancer Epidemiol Biomarkers Prev. (2013) 22:2161–81. doi: 10.1158/1055-9965.EPI-13-0621

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Ahmed D, Eide PW, Eilertsen IA, Danielsen SA, Eknaes M, Hektoen M, et al. Epigenetic and genetic features of 24 colon cancer cell lines. Oncogenesis. (2013) 2:e71. doi: 10.1038/oncsis.2013.35

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Corso G, Pascale V, Flauti G, Ferrara F, Marrelli D, Roviello F. Oncogenic mutations and microsatellite instability phenotype predict specific anatomical subsite in colorectal cancer patients. Eur J Hum Genet. (2013) 21:1383–8. doi: 10.1038/ejhg.2013.66

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Mouradov D, Domingo E, Gibbs P, Jorissen RN, Li S, Soo PY, et al. Survival in stage II/III colorectal cancer is independently predicted by chromosomal and microsatellite instability, but not by specific driver mutations. Am J Gastroenterol. (2013) 108:1785–93. doi: 10.1038/ajg.2013.292

CrossRef Full Text | Google Scholar