- 1Academician Workstation, Changsha Medical University, Changsha, China
- 2Hunan Key Laboratory of the Research and Development of Novel Pharmaceutical Preparations, Changsha Medical University, Changsha, China
- 3Geneis Beijing Co., Ltd., Beijing, China
- 4School of Chemistry and Chemical Engineering, Hainan Normal University, Haikou, China
Background: Linc00312 is dysregulated in nasopharyngeal carcinoma (NPC) and participates in the initiation and progression of NPC. Our previous studies suggested that linc00312 was able to enhance the sensitivity of NPC cells to irradiation and NPC patients with higher expression of linc00312 was associated with better short-term curative effect and overall survival. The single nucleotide polymorphisms (SNPs) of lncRNAs may influence the disease course and outcome by affecting the expression, secondary structure or function of lncRNAs. However, the role of SNPs in linc00312 on the occurrence and survival of NPC remains unknown.
Methods: We recruited 684 NPC patients and 823 healthy controls to evaluate the association between linc00312 SNPs and NPC susceptibility by using multivariate logistic regression analysis. Kaplan-Meier analysis and Cox proportional hazards regression were applied to assess the effect of linc00312 SNPs on the survival of NPC patients. The relative expression of linc00312 in NPC tissues was determined by real-time PCR. The interaction between linc00312 and mir-411-3p was explored by luciferase reporter assay. In silico prediction of the changes on linc00312 folding structure was conducted by RNAfold WebServer.
Result: We demonstrated that rs12497104 (G > A) GA genotype carriers had a higher risk than others for suffering from NPC (GA vs GG, OR = 1.437, P = 0.003). Besides, patients with rs12497104 AA genotype showed a poorer overall survival in contrast to GG genotype (AA vs GG, HR = 2.117, P = 0.011). In addition, the heterozygous carriers of rs15734 (G > A) and rs164966 (A > G) were correlated with decreased risk of NPC (GA vs GG, OR = 0.778, P = 0.031; GA vs AA, OR = 0.781, P = 0.033, respectively). We found that the three SNPs might influence the expression of linc00312 in a genotype specific feature. The local centroid secondary structure as well as the minimum free energy of linc00312 were changed following the candidate SNPs alterations. Besides, we revealed that the G to A alteration at rs12497104 disrupted the binding between mir-411-3p and linc00312.
Conclusion: Our results indicated genetic polymorphisms of linc00312 might serve as potential biomarkers for NPC carcinogenesis and prognosis.
Introduction
Nasopharyngeal carcinoma has a remarkable geographic distribution and is relatively predominant in Southern China. According to the statistics report of nasopharyngeal carcinoma (NPC) in 2018, more than 70% of new cases are in east and southeast Asia, with an age-standardized rate (world) of 3.0/100,000 in China to 0.4/100,000 in white populations (Bray et al., 2018). The etiology of NPC is multifactorial and is widely suspected to be complex interactions of genetic predisposition, Epstein–Barr virus infection, and environmental factors (Paul et al., 2018).
Distinguished from other types of head and neck cancer, NPC is more aggressive due to extensive local infiltration, early lymphatic spread, and high tendency of hematogenous dissemination. More than 70% of patients are classified as locoregionally advanced stages at the time of diagnosis (Zhang et al., 2019b). Currently, radiotherapy and platinum-based chemotherapy are the mainstay treatment modality for NPC (Chen et al., 2019). The 3-years failure-free survival is around 80% for NPC patients receiving induction chemotherapy plus concurrent chemoradiotherapy (Sun et al., 2016; Zhang et al., 2019a).
It is now well recognized that genetic aberrations play a vital role in the pathogenesis, progression and prognosis of NPC (Lin et al., 2014; Tang et al., 2018). Previous studies have reported several susceptibility and recurrent associated loci of NPC, for example, CYLD, TRAF3, MST1R, TNFRSF19, MECOM, and CDKN2A–CDKN2B (Bei et al., 2010; Dai et al., 2016; Li et al., 2017). As we can see, most of the identified loci are located within coding regions. Noteworthy, results from GWAS showed only 7% of diseases associated loci was located in protein-coding regions, indicating genetic variations in non-coding regions are likely to exert important functions in disease (Freedman et al., 2011).
Long non-coding RNAs (lncRNAs) lack the capability of protein-coding and are once considered to be simply transcriptional “noise”. Recently, emerging evidence have shown that lncRNAs participates in the regulation of many cellular biological processes (Lin and Yang, 2018). Deregulation of lncRNAs are closely related to tumorigenesis and prognosis of cancers (Wu and Hann, 2018). Single-nucleotide polymorphism (SNP) is one of the most commonly occurred type of genetic variants within lncRNA genes. The SNPs may exert its effect by affecting the expression, secondary structure or function of lncRNAs, thus influencing the disease course and outcome (Minotti et al., 2018).
In the present study, we focus on linc00312, which is dysregulated in NPC. The expression of linc00312 in NPC tissues is negatively correlated with tumor size (Zhang et al., 2013). Linc00312 exerts dual role on NPC cells as it inhibits cells proliferation, but increases cell adhesion and invasion (Huang et al., 2009). Our previous study revealed that overexpression of linc00312 enhanced the sensitivity of NPC cells to irradiation through regulating DNA damage response. NPC patients with higher expression of linc00312 were associated with longer overall survival (Guo et al., 2021). All of these indicates linc00312 might serve as a potential biomarker for NPC.
To date, the functional significance of SNPs in linc00312 locus is still unclear. Given the fact mentioned above, we speculate that the functional genetic variations in linc00312 may affect its expression and/or function, which finally influence the carcinogenesis and prognosis of NPC. In order to verify the conjecture, we selected three SNPs with potential regulatory feature of linc00312 and performed a case-control study to explore the role of linc00312 variations on NPC risk and survival.
Materials and Methods
Study Population
684 histopathology-confirmed NPC patients were enrolled at Hunan Provincial Cancer Hospital between 2014 and 2015. Meanwhile, 823 cancer-free volunteers taking physical examination at Xiangya Hospital were enrolled as controls. Moreover, we collected 82 NPC tissue biopsies from patients taking nasopharyngoscopy inspection. The patients who did not receive IMRT radiation technique or abandoned treatment or were out of contact were excluded in the survival analysis. All of the patients underwent platinum-based chemoradiotherapy. A total of 68–74 Gy irradiation was administered to the primary tumor for 7 weeks. The chemotherapy regimens in the present study including: platinum + paclitaxel (TP); platinum + 5-fluorouracil (FP); platinum + docetaxel (DP); platinum + gemcitabine (GP); cisplatin alone (DDP); nedaplatin alone (NDP). This study was approved by the Independent Ethical Committee of Institute of Clinical Pharmacology, Central South University (CTXY-140007-2) and all the participants signed the informed consent at the time of enrollment.
Candidate SNPs Selection
Three online bioinformatic databases, including ENCODE, ENSEMBL and lncRNASNP were applied to select potentially functional SNPs that located in miRNA binding site/3′-untranslated region (3′-UTR)/enhancer region/5′-untranslated region (5′-UTR)/alternative splice region/open chromatin region. Finally, we selected three SNPs of linc00312 [rs12497104 (G > A), rs15734 (G > A), rs164966 (A > G)] that have never been reported for further study.
Genotyping
Genomic DNA was isolated from peripheral blood lymphocytes and NPC tissues using the QIAamp DNA Mini and Blood Mini Kit (Qiagen Inc., Valencia, CA, United States). The concentrations and purity of DNA were measured using NanoDropTM 1000 spectrophotometer. Genotypes of the candidate SNPs were determined by using Sequenom MassARRAY iPLEX (Sequenom, Inc., San Diego, CA, United States). The call rate threshold was at least 95%.
Real-Time PCR
Total RNA was extracted from NPC tissues and reverse-transcribed to cDNA by using a PrimeScriptTM RT kit (Takara, Japan). LightCycler 480 system was utilized for amplification of cDNA using a SYBR Premix Ex TaqTM kit (Takara, Japan). The relative expression of linc00312 was calculated by 2–ΔΔCt method. The primers are as the following: GAPDH forward: 5′-ACAACTTTGGTATCGTGGAAGG-3′ and reverse: 5′-GCCATC ACGCCACAGTTTC-3′; linc00312 forward: 5′-GATCTATG GCCCATCATTCTTT-3′ and reverse: 5′-GTCCATCATGTAGC AAGCAGT-3′.
Dual Luciferase Reporter Assay
The human NPC cell lines HONE1 and HNE1 were cultured in RPMI 1640 medium (Invitrogen, Carlsbad, CA, United States) containing 10% FBS (Invitrogen, Carlsbad, CA, United States) and 1% penicillin-streptomycin (100 units/mL) at 37°C in a 5% CO2 incubator. The 3′-UTR of linc00312 covering wild-type of rs12497104 was directly synthesized by Genechem (Genechem, Shanghai, China) and cloned into XbaI site of pGL3 promoter vector to get the wild-type vector (WT). Site-directed mutagenesis at rs12497104 was utilized to get the mutant-type vector (MUT). DNA sequencing was performed in order to make sure of the sequence of these vectors. HONE1 and HNE1 cells were co-transfected with Renilla vector, WT\MUT reporter vector, mir-411-3p mimics\mimics control or mir-411-3p inhibitor\inhibitor control using Lipofectamine 3000 (Invitrogen, Carlsbad, CA, United States). Cells were harvested 48 h post transfection and assayed for luciferase activities using the Dual Luciferase Reporter Assay Kit (Promega, Madison, WI, United States).
Statistical Analysis
Chi-square test and Student’s t-test were used to compare the distributions of categorical variables and continuous variables between patients and controls. The odd ratios (ORs) and 95% confidence intervals (95% CI) for each genotype were calculated by multivariate logistic regression to determine the association between the selected SNPs and NPC susceptibility. Kaplan-Meier analysis and Cox proportional hazards regression were applied to evaluate the overall survival of NPC patients. The SPSS 19.0 was used for statistical analyses (Chicago, IL, United States) and P < 0.05 was set as the threshold of statistically significant.
Results
Population Characteristics
The detailed information of the NPC cases and controls were summarized in Supplementary Table 1, as we have previously reported (Guo et al., 2019). There is no distribution difference in the cases and controls for age and BMI (P > 0.05). However, the male-to-female ratio was higher in the cases than controls, which was in line with NPC epidemiology. The cases were more likely to have ever smoked (P < 0.05) and drunk patients (P < 0.05). To eliminate the influence of these factors, the covariant were further adjusted in multivariate logistics regression analysis.
Linc00312 SNPs Were Associated With NPC Susceptibility
The genotyping call rate of the selected SNPs were greater than 95%, and the genotype frequency of the candidate SNPs were in accordance with Hardy-Weinberg equilibrium, as shown in Table 1. We used multivariate logistic regression analysis to estimate the association between selected SNPs and NPC susceptibility. We found that individuals with rs12497104 GA genotype were associated with significantly increased risk of NPC (GA vs GG, OR = 1.437, P = 0.003) (Table 1). Individuals with rs15734 GA genotype showed a 22.2% lower risk of NPC than those with rs15734 GG genotype (GA vs GG, OR = 0.778, P = 0.031). As for rs164966, a 21.9% decreased risk was observed in rs164966 GA genotype carriers in comparison with rs164966 AA genotype carriers (GA vs AA, OR = 0.781, P = 0.033). Since gender might be a risk factor for NPC, we conducted the stratified analysis to eliminate the influence of gender. As shown in Supplementary Table 2, the ORs of each SNP were even more remarkable in the stratified group than in the whole cohort.
Table 1. Association between the candidate single nucleotide polymorphisms (SNPs) of linc00312 and nasopharyngeal carcinoma (NPC) susceptibility was analyzed by multivariate logistic regression.
Linc00312 rs12497104 AA Genotype Was Correlated With Poorer Overall Survival of NPC
We included 650 NPC patients in the survival analysis and the detailed therapeutic information was listed in Table 2. The median follow-up duration was 49.27 months (range, 3.5–60 months). The 3-years OS rate of patients with rs12497104 GG genotype were 85.6% and patients with AA genotype were 74.6%. The univariate analysis indicated that the covariates correlated with OS were rs12497104, age, gender, clinical stage, irradiation dose and chemotherapeutic regimen (P < 0.05). Multivariate analysis adjusting for these covariates by Cox proportional hazards regression model was used to assess the role of rs12497104. As shown in Figure 1, rs12497104 AA genotype carriers showed a poorer overall survival than GG genotype carriers (AA vs GG, HR = 2.117, P = 0.011). The mean survival time was 53.56 months for rs12497104 GG genotype carriers and 48.71 months for AA genotype patients. No significant difference in survival was noted in rs15734 and rs164966.
Figure 1. Linc00312 single nucleotide polymorphisms (SNPs) are associated with OS of nasopharyngeal carcinoma (NPC) patients. The survival curve of NPC patients with different genotypes of rs12497104 (A), rs13734 (B), and rs164966 (C).
Genotype-Specific Expression Effect of the Candidate SNPs
The three candidate SNPs are located in the regulatory region of linc00312 that could cause miRNA–lncRNA binding site gain or loss based on the bioinformatic prediction (Table 3). We speculated the SNPs might regulate the expression of linc00312 by affecting the interaction with miRNAs and we detected the expression level of linc00312 in NPC tissues with different genotypes. The result demonstrated that rs12497104 AA genotype was associated with a significantly lower expression of linc00312 in comparison with rs12497104 GG genotype (P < 0.05) (Figure 2A). The expression of linc00312 in rs15734 GA genotype carriers were higher than rs15734 GG genotype carriers (P < 0.05) (Figure 2B). Moreover, patients with rs164966 GA and GG genotypes all displayed a higher expression of linc00312 compared with rs164966 AA genotype carriers (P < 0.05) (Figure 2C).
Figure 2. Genotype-specific expression effect of the candidate SNPs. (A) The relative expression of linc00312 in different rs12497104 genotype carriers. (B) The relative expression of linc00312 in different rs15734 genotype carriers. (C) The relative expression of linc00312 in different rs164966 genotype carriers. *P < 0.05; **P < 0.01.
Actually, we verified that rs12497104 could destroy the interaction between linc00312 and mir-411-3p by luciferase reporter assay. The relative luciferase activity was significantly decreased when mir-411-3p mimics was co-transfected with the wild-type linc00312 (WT) reporter (P < 0.01). On the contrary, the luciferase activity was almost unchanged when mir-411-3p mimics was co-transfected with the rs12497104 mutant linc00312 (MUT) reporter (P > 0.05) (Figures 3A,B). What’s more, cotransfection of mir-411-3p inhibitor caused a notable gain on the luciferase activity of the WT reporter (P < 0.05) while the luciferase activity of MUT reporter displayed no visible change (P > 0.05) (Figures 3C,D).
Figure 3. rs12497104 destroys the binding between mir-411-3p and linc00312. (A) Relative luciferase activity of the wild type linc00312 (WT) and rs12497104 mutant type linc00312 (MUT) reporter vector that co-transfected with mir-411-3p mimics or mir-411-3p control in HONE1 cells. (B) Relative luciferase activity of the WT and MUT reporter vector that co-transfected with mir-411-3p mimics or mir-411-3p control in HNE1 cells. (C) Relative luciferase activity of the WT and MUT reporter vector that co-transfected with mir-411-3p inhibitor or mir-411-3p inhibitor control in HONE1 cells. (D) Relative luciferase activity of the WT and MUT reporter vector that co-transfected with mir-411-3p inhibitor or mir-411-3p inhibitor control in HNE1 cells. *P < 0.05; **P < 0.01.
In silico Prediction of the Candidate SNPs on Linc00312 Folding Structure
By in silico analysis, we predicted weather the candidate SNPs could possibly affecting the folding structure of linc00312. As shown in Figure 4, the local centroid secondary structure and minimum free energy of linc00312 were changed following the candidate SNPs alterations. The minimum free energy of the local centroid secondary structure of linc00312 changed from −121.60 to −98.42, −127.20, and −113.50 kcal/mol for rs12497104, rs15734, and rs164966 variations, respectively.
Figure 4. In silico prediction of the candidate SNPs on linc00312 folding structure. (A) The local centroid secondary structure of the wild type linc00312. (B) The local centroid secondary structure of rs12497104 mutant linc00312. (C) The local centroid secondary structure of rs15734 mutant linc00312. (D) The local centroid secondary structure of rs164966 mutant linc00312. The color of each base represents base-pair probabilities. (E) The mountain plot of the wild type linc00312. (F) The mountain plot of rs12497104 mutant linc00312. (G) The mountain plot of rs15734 mutant linc00312. (H) The mountain plot of rs164966 mutant linc00312. The mountain plot presents the MFE structure, the thermodynamic ensemble of RNA structures, the centroid structure and the positional entropy for each position. All of the results were computed by RNAfold WebServer.
Discussion
Considering the critical role of linc00312 in NPC, we wondered if the SNPs in linc00312 might serve as potential biomarkers for NPC. In the present study, we identified three SNPs of linc00312 (rs12497104, rs15734, and rs164966) were associated with susceptibility of NPC. In the following survival analysis, we revealed rs12497104 was an independent risk factor for the prognosis of NPC.
Not only for NPC, linc00312 also displayed a vital role in other types of cancers via regulating diverse cellular processes. As for non-small cell lung cancer (NSCLC), patients with lower expression of linc00312 were correlated with larger tumor size, advanced clinical stages as well as shorter overall survival. Overexpression of linc00312 was able to repress cell growth and induce apoptosis of NSCLC cells (Zhu et al., 2017). The other study demonstrated that linc00312 could promote metastasis and angiogenesis in lung cancer via interacting with the transcription factor YBX1 (Peng et al., 2018). Linc00312 also regulated oral fibrogenesis by binding to YBX1 in oral cavity (Yu et al., 2020). Evidence have shown that linc00312 inhibited cell growth and migration of thyroid cancer cells through suppressing the PI3K/Akt and MMP-9 (Min et al., 2018). By targeting miR-197-3p and miR-21, linc00312 functioned as a competing endogenous RNA (ceRNA) and repressed cell migration and invasion of bladder cancer, thyroid cancer and colorectal cancer (Wang et al., 2016; Liu et al., 2017; Li et al., 2018). Linc00312 could induce cell cycle arrest of hepatocellular carcinoma cell and leading to the suppression of proliferation (Wu et al., 2018). In addition, linc00312 participated in drug resistance as it sensitized ovarian cancer cells to cisplatin by activating the Bcl-2/Caspase-3 pathway (Zhang et al., 2018).
A great number of studies have indicated SNPs of lncRNA have great potential as predictive markers for carcinogenesis, prognosis, and drug resistance (Gao and Wei, 2017; Minotti et al., 2018). Our previous findings have suggested NPC patients with lncRNA GAS5 rs2067079 CT genotypes were linked with an obviously increased risk of severe myelosuppression during concurrent radiochemotherapy period compared with CC genotype carriers (Guo et al., 2017). In addition, NPC patients with lncRNA MEG3 rs10132552 CC genotype were vulnerable to suffering chemoradiotherapy induced anemia (Wang et al., 2017). Nevertheless, nothing has been known about linc00312 SNPs and NPC so far.
Unveiling the molecular mechanism of lncRNAs SNPs that contributing to disease risk is quite important for better understanding the pathogenesis of disease. The SNPs may regulate the function of lncRNA through interfering the interaction between lncRNAs and transcription factors, miRNAs or other binding protein partners. Evidence have proven that linc00673 rs11655237 G > A alteration created a binding site for miR-1231, which diminished the function of linc00673 and lead to increased risk of pancreatic cancer (Zheng et al., 2016). The risk variant of rs11672691 in lncRNA PCAT19 suppresses binding of transcription factor NKX3.1 to the promoter of PCAT19-short, resulting in PCAT19-long activation and prostate cancer growth and metastasis (Hua et al., 2018). Another prostate cancer-associated variant at lncRNA PCAT1 rs7463708 increases binding of transcription factor ONECUT2 to the PCAT1 promoter, resulting in upregulation of PCAT1 and prostate transformation (Guo et al., 2016). In colon cancer, by interacting with the CFIm complex with allele specific affinities, lncRNA CCAT2 rs6983267 regulates the alternative splicing of glutaminase, resulting in reprogramming of cancer metabolism (Redis et al., 2016). To data, literature that clearly clarify the molecular mechanism of cancer related SNPs in lncRNA is very scarce.
Bioinformatics approaches have usually been applied to predict the potential role of SNPs in lncRNAs. In the present study, we applied LncRNASNP, ENSEMBL, and RNAfold for SNP selection and function prediction (Gruber et al., 2008; Miao et al., 2018). Actually, bioinformatics approach is a cost-effective way to screening causal SNPs. The prediction result indicated all of the three SNPs of linc00312 (rs12497104, rs15734, and rs164966) had an eQTL trait and were likely to create or destroy the binding sites of miRNAs with linc00312. In light of this, we wondered if the SNPs could affect the function/expression of linc00312. It has been reported that 75% of the lncRNA SNPs affected the expression level of lncRNA (Kumar et al., 2013). So, we examined the expression of linc00312 with different genotypes and found the genotype specific expression feature. The rs12497104 AA genotype carriers had a lower expression of linc00312 and were subjected to enhanced risk of NPC and poorer survival. Likewise, the patients with GA genotype of rs15734 and rs164966 showed a lower expression of linc00312 and were correlated with decreased risk of NPC. In addition, the local secondary structure and minimum free energy of linc00312 was affected by the SNPs.
Increasing evidences have indicated that lncRNA could function as ceRNA to regulate mRNA by sponging miRNA. For the first time, we proved that rs12497104 could interfere the binding between linc00312 and mir-411-3p. This finding provided one probable mechanism by which rs12497104 affected the function of linc00312. As far as we know, there was only two miRNAs (mir-197-3p, mir-21) have been identified to interact with linc00312 (Wang et al., 2016; Liu et al., 2017; Li et al., 2018). To date, the biological function of mir-411-3p is still largely unknown. A recent study reported mir-411-3p interacted with lncRNA ANRIL and inhibited the malignant proliferation and tumor stem cell like property of multiple myeloma (Wang et al., 2020). Another study demonstrated lncRNA TTN-AS1 could function as the mir-411-3p sponge in oral squamous cell carcinoma (OSCC) and mir-411-3p exerted the inhibitory functions on OSCC growth (Fu et al., 2020). In addition, lncRNA CDKN2B-AS1 could interact with mir-411-3p and contribute to carcinogenesis in ovarian cancer (Wang et al., 2019). The role of mir-411-3p in NPC still needs further investigation.
In conclusion, we demonstrated the SNPs of linc00312 were associated with NPC susceptibility and survival possibly by influencing the expression of linc00312. Our findings may shed some light on the biomarkers for predicting NPC risk and prognosis.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by the Independent Ethical Committee of Institute of Clinical Pharmacology, Central South University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
ZG performed the experiments, analyzed the data, and drafted the manuscript. M-HB helped for sample collection. Y-XF, YZ, and H-YL performed the quality control of data. X-LZ and BW provided the technical assistance. B-SH and Q-QL participated in manuscript editing. J-YL participated in manuscript review. All authors read and approved the final manuscript.
Funding
This work was supported by the Outstanding Youth Project of Hunan Provincial Education Department (No. 20B075), Hunan Provincial Innovation Platform and Open Fund Project (No. 20K015), and Project of Improving the Basic Scientific Research ability of Young College Teachers in Guangxi (Nos. 2021KY1942 and 2021KY1953).
Conflict of Interest
Q-QL was employed by Genesis Beijing Co., Ltd.
The remaining 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 W. Zhang for helping us collect specimens.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.698558/full#supplementary-material
References
Bei, J. X., Li, Y., Jia, W. H., Feng, B. J., Zhou, G., Chen, L. Z., et al. (2010). A genome-wide association study of nasopharyngeal carcinoma identifies three new susceptibility loci. Nat. Genet. 42, 599–603.
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. doi: 10.3322/caac.21492
Chen, Y. P., Chan, A. T. C., Le, Q. T., Blanchard, P., Sun, Y., and Ma, J. (2019). Nasopharyngeal carcinoma. Lancet 394, 64–80.
Dai, W., Zheng, H., Cheung, A. K., Tang, C. S., Ko, J. M., Wong, B. W., et al. (2016). Whole-exome sequencing identifies MST1R as a genetic susceptibility gene in nasopharyngeal carcinoma. Proc. Natl. Acad. Sci. U. S. A. 113, 3317–3322.
Freedman, M. L., Monteiro, A. N., Gayther, S. A., Coetzee, G. A., Risch, A., Plass, C., et al. (2011). Principles for the post-GWAS functional characterization of cancer risk loci. Nat. Genet. 43, 513–518. doi: 10.1038/ng.840
Fu, S. W., Zhang, Y., Li, S., Shi, Z. Y., Zhao, J., and He, Q. L. (2020). LncRNA TTN-AS1 promotes the progression of oral squamous cell carcinoma via miR-411-3p/NFAT5 axis. Cancer Cell Int. 20:415.
Gao, P., and Wei, G. H. (2017). Genomic Insight into the Role of lncRNA in Cancer Susceptibility. Int. J. Mol. Sci. 18:1239. doi: 10.3390/ijms18061239
Gruber, A. R., Lorenz, R., Bernhart, S. H., Neuböck, R., and Hofacker, I. L. (2008). The Vienna RNA websuite. Nucleic Acids Res. 36, W70–W74.
Guo, H., Ahmed, M., Zhang, F., Yao, C. Q., Li, S., Liang, Y., et al. (2016). Modulation of long noncoding RNAs by risk SNPs underlying genetic predispositions to prostate cancer. Nat. Genet. 48, 1142–1150. doi: 10.1038/ng.3637
Guo, Z., Wang, Y. H., Xu, H., Yuan, C. S., Zhou, H. H., Huang, W. H., et al. (2021). LncRNA linc00312 suppresses radiotherapy resistance by targeting DNA-PKcs and impairing DNA damage repair in nasopharyngeal carcinoma. Cell Death Dis. 12:69.
Guo, Z., Wang, Y., Zhao, Y., Jin, Y., An, L., Wu, B., et al. (2017). Genetic polymorphisms of long non-coding RNA GAS5 predict platinum-based concurrent chemoradiotherapy response in nasopharyngeal carcinoma patients. Oncotarget 8, 62286–62297. doi: 10.18632/oncotarget.19725
Guo, Z., Wang, Y., Zhao, Y., Jin, Y., An, L., Xu, H., et al. (2019). A Functional 5′-UTR Polymorphism of MYC Contributes to Nasopharyngeal Carcinoma Susceptibility and Chemoradiotherapy Induced Toxicities. J. Cancer 10, 147–155. doi: 10.7150/jca.28534
Hua, J. T., Ahmed, M., Guo, H., Zhang, Y., Chen, S., Soares, F., et al. (2018). Risk SNP-Mediated Promoter-Enhancer Switching Drives Prostate Cancer through lncRNA PCAT19. Cell 174, 564–575.e18.
Huang, C., Wu, M., Tang, Y., Li, X., Ouyang, J., Xiao, L., et al. (2009). NAG7 promotes human nasopharyngeal carcinoma invasion through inhibition of estrogen receptor alpha and up-regulation of JNK2/AP-1/MMP1 pathways. J. Cell Physiol. 221, 394–401. doi: 10.1002/jcp.21867
Kumar, V., Westra, H. J., Karjalainen, J., Zhernakova, D. V., Esko, T., Hrdlickova, B., et al. (2013). Human disease-associated genetic variation impacts large intergenic non-coding RNA expression. PLoS Genet. 9:e1003201. doi: 10.1371/journal.pgen.1003201
Li, G., Wang, C., Wang, Y., Xu, B., and Zhang, W. (2018). LINC00312 represses proliferation and metastasis of colorectal cancer cells by regulation of miR-21. J. Cell. Mol. Med. 22, 5565–5572. doi: 10.1111/jcmm.13830
Li, Y. Y., Chung, G. T., Lui, V. W., To, K. F., Ma, B. B., Chow, C., et al. (2017). Exome and genome sequencing of nasopharynx cancer identifies NF-κB pathway activating mutations. Nat. Commun. 8:14121.
Lin, C., and Yang, L. (2018). Long Noncoding RNA in Cancer: wiring Signaling Circuitry. Trends Cell Biol. 28, 287–301. doi: 10.1016/j.tcb.2017.11.008
Lin, D. C., Meng, X., Hazawa, M., Nagata, Y., Varela, A. M., Xu, L., et al. (2014). The genomic landscape of nasopharyngeal carcinoma. Nat. Genet. 46, 866–871.
Liu, K., Huang, W., Yan, D. Q., Luo, Q., and Min, X. (2017). Overexpression of long intergenic noncoding RNA LINC00312 inhibits the invasion and migration of thyroid cancer cells by down-regulating microRNA-197-3p. Biosci. Rep. 37:BSR20170109.
Miao, Y. R., Liu, W., Zhang, Q., and Guo, A. Y. (2018). lncRNASNP2: an updated database of functional SNPs and mutations in human and mouse lncRNAs. Nucleic Acids Res. 46:D276–D280.
Min, X., Liu, K., Zhu, H., and Zhang, J. (2018). Long Noncoding RNA LINC003121 Inhibits Proliferation and Invasion of Thyroid Cancer Cells by Suppression of the Phosphatidylinositol-3-Kinase (PI3K)/Akt Signaling Pathway. Med. Sci. Monit. 24, 4592–4601. doi: 10.12659/msm.908652
Minotti, L., Agnoletto, C., Baldassari, F., Corrà, F., and Volinia, S. (2018). SNPs and Somatic Mutation on Long Non-Coding RNA: new Frontier in the Cancer Studies? High Throughput 7:34. doi: 10.3390/ht7040034
Paul, P., Deka, H., Malakar, A. K., Halder, B., and Chakraborty, S. (2018). Nasopharyngeal carcinoma: understanding its molecular biology at a fine scale. Eur. J. Cancer Prev. 27, 33–41. doi: 10.1097/cej.0000000000000314
Peng, Z., Wang, J., Shan, B., Li, B., Peng, W., Dong, Y., et al. (2018). The long noncoding RNA LINC00312 induces lung adenocarcinoma migration and vasculogenic mimicry through directly binding YBX1. Mol. Cancer 17:167.
Redis, R. S., Vela, L. E., Lu, W., Ferreira de Oliveira, J., Ivan, C., Rodriguez-Aguayo, C., et al. (2016). Allele-Specific Reprogramming of Cancer Metabolism by the Long Non-coding RNA CCAT2. Mol. Cell 61:640.
Sun, Y., Li, W. F., Chen, N. Y., Zhang, N., Hu, G. Q., Xie, F. Y., et al. (2016). Induction chemotherapy plus concurrent chemoradiotherapy versus concurrent chemoradiotherapy alone in locoregionally advanced nasopharyngeal carcinoma: a phase 3, multicentre, randomised controlled trial. Lancet Oncol. 17, 1509–1520.
Tang, X. R., Li, Y. Q., Liang, S. B., Jiang, W., Liu, F., Ge, W. X., et al. (2018). Development and validation of a gene expression-based signature to predict distant metastasis in locoregionally advanced nasopharyngeal carcinoma: a retrospective, multicentre, cohort study. Lancet Oncol. 19, 382–393. doi: 10.1016/s1470-2045(18)30080-9
Wang, M., Zhao, H. Y., Zhang, J. L., Wan, D. M., Li, Y. M., and Jiang, Z. X. (2020). Dysregulation of LncRNA ANRIL mediated by miR-411-3p inhibits the malignant proliferation and tumor stem cell like property of multiple myeloma via hypoxia-inducible factor 1α. Exp. Cell Res. 396:112280. doi: 10.1016/j.yexcr.2020.112280
Wang, Y. Y., Wu, Z. Y., Wang, G. C., Liu, K., Niu, X. B., Gu, S., et al. (2016). LINC00312 inhibits the migration and invasion of bladder cancer cells by targeting miR-197-3p. Tumour Biol. 37, 14553–14563. doi: 10.1007/s13277-016-5303-8
Wang, Y., Guo, Z., Zhao, Y., Jin, Y., An, L., Wu, B., et al. (2017). Genetic polymorphisms of lncRNA-p53 regulatory network genes are associated with concurrent chemoradiotherapy toxicities and efficacy in nasopharyngeal carcinoma patients. Sci. Rep. 7:8320.
Wang, Y., Huang, Y., Liu, H., Su, D., Luo, F., and Zhou, F. (2019). Long noncoding RNA CDKN2B-AS1 interacts with miR-411-3p to regulate ovarian cancer in vitro and in vivo through HIF-1a/VEGF/P38 pathway. Biochem. Biophys. Res. Commun. 514, 44–50. doi: 10.1016/j.bbrc.2019.03.141
Wu, J., and Hann, S. S. (2018). Functions and Roles of Long-Non-Coding RNAs in Human Nasopharyngeal Carcinoma. Cell Physiol. Biochem. 45, 1191–1204. doi: 10.1159/000487451
Wu, J., Zhou, X., Fan, Y., Cheng, X., Lu, B., and Chen, Z. (2018). Long non-coding RNA 00312 downregulates cyclin B1 and inhibits hepatocellular carcinoma cell proliferation in vitro and in vivo. Biochem. Biophys. Res. Commun. 497, 173–180. doi: 10.1016/j.bbrc.2018.02.049
Yu, C. H., Fang, C. Y., Yu, C. C., Hsieh, P. L., Liao, Y. W., Tsai, L. L., et al. (2020). LINC00312/YBX1 Axis Regulates Myofibroblast Activities in Oral Submucous Fibrosis. Int. J. Mol. Sci. 21:2979. doi: 10.3390/ijms21082979
Zhang, C., Wang, M., Shi, C., Shi, F., and Pei, C. (2018). Long non-coding RNA Linc00312 modulates the sensitivity of ovarian cancer to cisplatin via the Bcl-2/Caspase-3 signaling pathway. Biosci. Trends 12, 309–316. doi: 10.5582/bst.2018.01052
Zhang, W., Huang, C., Gong, Z., Zhao, Y., Tang, K., Li, X., et al. (2013). Expression of LINC00312, a long intergenic non-coding RNA, is negatively correlated with tumor size but positively correlated with lymph node metastasis in nasopharyngeal carcinoma. J. Mol. Histol. 44, 545–554. doi: 10.1007/s10735-013-9503-x
Zhang, Y., Chen, L., Hu, G. Q., Zhang, N., Zhu, X. D., Yang, K. Y., et al. (2019a). Gemcitabine and Cisplatin Induction Chemotherapy in Nasopharyngeal Carcinoma. N. Engl. J. Med. 381, 1124–1135.
Zhang, Y., Sun, Y., and Ma, J. (2019b). Induction gemcitabine and cisplatin in locoregionally advanced nasopharyngeal carcinoma. Cancer Commun. 39:39. doi: 10.1186/s40880-019-0385-5
Zheng, J., Huang, X., Tan, W., Yu, D., Du, Z., Chang, J., et al. (2016). Pancreatic cancer risk variant in LINC00673 creates a miR-1231 binding site and interferes with PTPN11 degradation. Nat. Genet. 48, 747–757. doi: 10.1038/ng.3568
Keywords: nasopharyngeal carcinoma, linc00312, polymorphism, susceptibility, survival
Citation: Guo Z, Bao M-H, Fan Y-X, Zhang Y, Liu H-Y, Zhou X-L, Wu B, Lu Q-Q, He B-S, Nan X-Y and Lu J-Y (2021) Genetic Polymorphisms of Long Non-coding RNA Linc00312 Are Associated With Susceptibility and Predict Poor Survival of Nasopharyngeal Carcinoma. Front. Cell Dev. Biol. 9:698558. doi: 10.3389/fcell.2021.698558
Received: 21 April 2021; Accepted: 07 June 2021;
Published: 16 July 2021.
Edited by:
Tao Huang, Shanghai Institute of Nutrition and Health (CAS), ChinaReviewed by:
Guohua Huang, Shaoyang University, ChinaBing Wang, Anhui University of Technology, China
Copyright © 2021 Guo, Bao, Fan, Zhang, Liu, Zhou, Wu, Lu, He, Nan and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jiao-Yang Lu, lulujiaoyang@163.com; Xu-Ying Nan, nanxuying@163.com