- Primary research
- Open Access
Polymorphism in lncRNA AC016683.6 and its interaction with smoking exposure on the susceptibility of lung cancer
Cancer Cell International volume 18, Article number: 91 (2018)
Long non-coding RNAs play pivotal roles in the carcinogenesis of multiple types of cancers. This study is firstly to evaluate influence of rs4848320 and rs1110839 polymorphisms in long non-coding RNA AC016683.6 on the susceptibility of lung cancer.
The present study was a hospital-based case–control study with 434 lung cancer patients and 593 cancer-free controls. Genotyping of the two SNPs detected by Taqman® allelic discrimination method.
There were no statistically significant associations between rs4848320 and rs1110839 polymorphisms in AC016683.6 and risk of lung cancer in overall population. However, in the smoking population, rs4848320 and rs1110839 polymorphisms significantly increased the risk of lung cancer in dominant and homozygous models (Rs4848320: P = 0.029; Rs1110839: P = 0.034), respectively. In male population, rs1110839 genetic variant was related to the risk of lung cancer in all genetic models (GG vs. TT: P = 0.008; Dominant model: P = 0.029; Recessive model: P = 0.027) rather than heterozygous model. The crossover analyses provided rs4848320 and rs1110839 risk genotypes carriers combined with smoking exposure 2.218-fold, 1.755-fold increased risk of lung cancer (Rs4848320: P = 0.005; Rs1110839: P = 0.017). Additionally, there were significantly positive multiplicative interaction of rs4848320 polymorphism with smoking status, with adjusted OR of 2.244 (1.162–4.334), but rs1110839 polymorphism did not exist.
Rs4848320 and rs1110839 polymorphisms may be associated with lung cancer susceptibility. Interaction of rs4848320 risk genotypes with smoking exposure may strengthen the risk effect on lung cancer.
Lung cancer is one of the most frequent malignant neoplasms in the world according to previous data gathered from different countries including China [1, 2]. The crude incidence and mortality of lung cancer were increasing rapidly over the past 30 years in China . The occurrence and development of lung cancer may result from complicated interactions of environmental and genetic factors through previous epidemiological studies, which has considered tobacco exposure as a major environmental risk factor. With the in-depth exploration of genome-wide association studies (GWAS), multiple studies for underlying associations between single nucleotide polymorphisms (SNPs) in non-coding RNA genes and various cancers has gradually emerged with respect to genetic risk factors.
Long non-coding RNAs (lncRNAs), a kind of non-coding RNAs of more than 200 nucleotides in length, possesses the properties of lack transcriptional protein-coding function and open-reading frames [4,5,6]. On account of cellular function, lncRNAs could also be characteristically classified into tumor-suppressive and oncogenic lncRNAs in a way similar to the protein-coding genes . LncRNAs may act as gene regulators through intricate mechanisms, including post-transcriptional, trans- and cis- gene regulation in carcinogenic pathways [8,9,10,11]. Most of previous studies have focused on revealing the critical role of lncRNAs for cancer risk, prognosis, diagnosis and targeted therapy, implying lncRNAs have shifted to the forefront of human cancer research.
Paired-box gene 8 (PAX8), one of members of the pair-box (PAX) gene family, is mapped on chromosome 2q13, which has multiple oncogenic properties such as inhibition of apoptosis to promote the survival of tumor cells, activation of BCl2 transcription and suppression of p53 [12,13,14]. A large number of studies have shown that PAX8 strikingly expressed in several types of tumors, exerts its biological functions and ultimately affects the development of cancer . PAX2/5/8 could promote cell proliferation, but its transcription factors play conservative roles in affecting cell apoptosis [16,17,18]. Besides, specific protein–protein interactions could regulate the activity of PAX transcription factors. For example, activation of PAX8 transcriptional activity may be occur on condition that PAX2/5/8 interact with the tumor suppressor pRb . Additionally, several studies had identified that the most common genetic event is the PAX8-PPARγ rearrangement except for the RAS gene mutations in follicular thyroid gland tumour [20,21,22,23,24]. PAX8-PPARγ rearrangement may affect the function of PPARγ protein and leads to impairment of normal PPARγ function eventually, which loss of the functions would result in uncontrolled cell growth [25, 26]. PAX8 may also bring about the elimination of the PTEN gene’s inhibitory effect by acting on the AKT to immunologically activate the PI3K signaling pathway . Simultaneous activation of the signal pathways is highly carcinogenic, which may cause distant metastasis of tumour cells and the occurrence of locally infiltrating follicular carcinoma. Particularly, Kanteti et al.  reported that PAX8 is overexpressed in non-small cell lung cancer (NSCLC), whereas PAX5 protein is predominantly expressed in small cell lung cancer (SCLC) by analyzing the expression profile of the PAX gene family, and speculated that PAX8 is likely to interfere with p53 transcription and contribute to the development of tumor. Silencing of PAX8 expression may impair viability and motility of NSCLC cells in A549 cells . Previously these studies had laid a solid foundation for gathering evidence to examine the hypothesis.
Recently, two studies has evaluated the association of expression quantitative trait loci (eQTL) SNPs in lncRNA AC016683.6, which located in the upstream region of PAX8 gene and could alter the expression of PAX8, with the susceptibility of cervical cancer  and the prognosis of hepatocellular carcinoma . In this study, we hypothesized that rs4848320 and rs1110839 polymorphisms in AC016683.6 may shape the risk of lung cancer by supporting of the above evidence. To verify this hypothesis, we implemented a hospital-based case–control study consisting of 434 lung cancer cases and 593 cancer-free controls to evaluate the associations between the two SNPs and susceptibility of lung cancer. We also investigated the potential interaction of the two SNPs with smoking exposure on lung cancer risk, which was valuable for exploring cancer etiology and proposing improvement of environmental risk factors for cancer prevention.
Materials and methods
Study subjects and sample data collection
The hospital-based case–control study was conducted in Shenyang city of Liaoning Province, which located in the northeast China. All enrolled subjects were unrelated ethnic Han Chinese. Totally, 434 newly diagnosed lung cancer patients (111 males and 323 females) and 593 cancer-free subjects (180 males and 413 females) enrolled in the case and control group, respectively. The included criteria of the selected cases were as follows: (a) patients who examined through histopathological confirmation and newly diagnosed lung cancer; (b) The patients had no previous cancers; (c) These patients had not experienced radiotherapy or chemotherapy before. The included criteria of controls was consistent with the (b) (c) of the above criteria for cases. A small number of controls suffered mainly from coronary heart disease, hypertension and diabetes. Meanwhile, the recruited controls were from medical examination centers of the hospitals in the same period and matched to the selected patients on age (± 5 years) to ensure comparability among subjects. The study has acquired the approval of the Institutional Review Board of China Medical University. All enrolled participants or every participant’s representatives signed the informed consent according to relevant regulations. Each participants donated 10 ml peripheral blood and possess capability for undergoing a 1.5 h interview to complete the demographic data by filling out questionnaire. Additionally, a subject who smoked more than 100 cigarettes in the whole lifetime is considered as a smoker, otherwise, he or she is regarded as a non-smoker.
SNP selection and genotyping
Based upon the data from Regulome Database, the study selected the rs4848320 and rs1110839 SNPs in PAX8-AS1 (http://www.regulomedb.org) , which are meet the criteria of minor allele frequency (MAF) > 0.05 in Han Chinese population [32, 33]. Following the previous method , we isolated genomic DNA samples from the vein blood of each participants by Phenol–chloroform extraction. Genotyping of the two SNPs conducted by using an Applied Biosystems 7500 FAST Real-Time PCR System (Foster City, CA, USA) utilizing Taqman® allelic discrimination (Applied Biosystems, Foster City, CA) with commercial primer–probe set. The reaction condition of real-time PCR was set at 10 min for 95 °C, 30 s at 92 °C and 1 min at 60 °C for all 47 cycles.
T test, χ2 test and logistic regression analysis tested the differences between the case and control group according to type of demographic variables (including categorical and continuous variables). A goodness-of-fit χ2 test was performed to estimate the Hardy–Weinberg equilibrium (HWE) of the control group. The linear-by-linear association of χ2 test was conducted to assess trend analysis in a similar dose-dependent manner. The odds ratios (ORs) and their 95% confidence intervals (CIs) were used to evaluate the relationships between rs4848320, rs1110839 genotypes and lung cancer risk by means of logistic regression analysis. Crossover analysis was to examine interaction between the SNPs and smoking exposure. The multiplicative interaction was evaluated by OR and its 95% CI with the unconditional logistic regression model. Assessment of the addictive interaction applied relative excess risk due to interaction (RERI), synergy index (S) and attributable proportion due to interaction (AP). If the 95% CI of RERI/AP did not contain 0 and the 95% CI of S did not contain 1, the addictive interaction was existed . All statistical analyses was obtained by performing SPSS 22.0 software (IBM SPSS, Inc. Chicago, IL, USA), and P < 0.05 was regarded as a statistically significant result under conditions of the all tests with two-sided.
Table 1 clearly presented the demographic characteristics of 593 controls and 434 lung cancer patients, which covers 331 non-small cell lung cancer (NSCLC), 233 adenocarcinoma (AD), 89 squamous cell carcinoma (SQ) and 61 small cell lung cancer (SCLC), respectively. The age mean ± SD of case group and control group were 57.97 ± 11.85 and 56.67 ± 15.88, respectively. There were no statistical differences on age and gender of the case and control groups. However, the distribution of smoking status in both groups presented the significant difference (P = 0.044). Therefore, all the further statistical analyses adjusted by age, gender and smoking status. Genotype frequency distributions of the two SNPs in control group was in agreement with HWE (P = 0.556 for rs4848320, P = 0.643 for rs1110839), which indicated selected subjects were a randomly representative sample from the overall population.
Genotype distribution and lung cancer susceptibility
Table 2 summarized the associations between the SNPs genotypes and the susceptibility to lung cancer as well as NSCLC. No significantly statistical significance was in all genetic (CT vs. CC: OR = 1.036, 95% CI = 0.788–1.362, P = 0.800; TT vs. CC: OR = 1.041, 95% CI = 0.482–2.247, P = 0.918; Dominant model: OR = 1.036, 95% CI = 0.795–1.352, P = 0. 0.792; Recessive model: OR = 1.030, 95% CI = 0.479–2.212, P = 0.940, adjusted by age, gender and smoking) models for rs4848320 polymorphism. Rs1110839 genetic variation was not significantly statistical significance in all (TG vs. TT: OR = 1.083, 95% CI = 0.830–1.413, P = 0.559; GG vs. TT: OR = 1.079, 95% CI = 0.714 –1.632, P = 0.717; Dominant model: OR = 1.082, 95% CI = 0.839–1.394, P = 0.543; Recessive model: OR = 1.034, 95% CI = 0.702–1.525, P = 0.864, adjusted by age, gender and smoking) models. Moreover, similar results occurred in NSCLC, lung adenocarcinoma and squamous cell carcinoma subgroups (Tables 2 and 3). Subsequently, we estimated the combined effect of variant alleles (Rs4848320-T and Rs1110839-G) on the risk of lung cancer and NSCLC (shown in Table 4). We obtained that lung cancer and NSCLC risks was not remarkable significance with the increasing number of risk genotypes of the rs4848320 and rs1110839 in a similar dose-dependent manner (LC: P = 0.893; NSCLC: P = 0.964). The results of stratified analyses by smoking exposure, gender and age presented in Tables 5, 6 and 7, respectively. In the smoking population, patients with CT and TT genotypes were significantly increased with the risk of lung cancer in dominant (Rs4848320: OR = 1.929, 95% CI = 1.070–3.479, P = 0.029, adjusted by age and gender) model rather than other three models. Rs110839 GG genotype was significantly associated with increased risk of lung cancer (OR = 3.224, 95% CI = 1.089–9.544, P = 0.034). In male population, rs1110839 genotypes were related to significantly increased risk of lung cancer in all genetic models (GG vs. TT: OR = 3.096, 95% CI = 1.343–7.134, P = 0.008; Dominant model: OR = 1.804, 95% CI = 1.061–3.067, P = 0.029; Recessive model: OR = 2.376, 95% CI = 1.102–5.121, P = 0.027) except for heterozygous model. However, the results of non-smoking, age and female subgroups did not reveal significant effect on lung cancer susceptibility.
Interaction between SNPs and smoking exposure
Table 8 provided the results of the crossover analysis. TT and CT genotypes carriers with smoking exposure were significantly associated with 2.218-fold increased risk to lung cancer (P = 0.005) for rs4848320. Additionally, rs1110839 TG and GG genotypes carriers with smoking exposure 1.755-fold increased risk to lung cancer (P = 0.017). However, the parameters of additive interaction presented that rs4848320 and rs1110839 polymorphisms had no additive interaction with smoking exposure (shown in Table 9). The OR value of multiplicative interaction of rs4848320 risk genotypes with smoking exposure more than 1 indicated that the interaction was positive multiplicative interaction, with a P value of 0.016 (P < 0.05) (Table 10). In other words, when the rs4848320 TT and CT genotypes coexist with smoking exposure, the interaction between the two risk factors becomes stronger, and its biological significance is a synergistic effect.
The current case–control study initially investigated the associations between polymorphisms in AC016683.6 and lung cancer susceptibility. We obtained that rs4848320 and rs1110839 polymorphisms significantly increased the risk of lung cancer in the smoking and male population. Subsequently, the crossover analyses provided the two SNPs risk genotypes carriers with smoking exposure increased the risk of lung cancer. Moreover, the interaction between rs4848320 risk genotypes and smoking was positive multiplicative interaction.
Genome-wide cancer mutation analyses are uncovering how the particular properties of lncRNAs to contribute crucial cancer phenotypes through the regulation of lncRNAs transcription , suggesting that lncRNAs are powerful modulators, involving chromatin interaction, transcriptional regulation, and DNA, Protein, RNA processing . McGinnis et al. initially discovered the paired box gene family in Drosophila . Each member among the family is a very conservative transcriptional regulator and actively participates in the sophisticated regulation of intracellular signal transduction. The cooperative regulation among the PAX genes plays a crucial role in different biological processes such as the transcription of development-related genes, the promotion of cell proliferation and differentiation, the inhibition of apoptosis, the induction of directional metastasis of precursor cells as well as selective processing of different splice sites of the precursor mRNA [38,39,40]. PAX8, a product of the pair-box gene family, plays a complex role in the developmental of various disease including cancers . Overwhelming evidence elucidated that PAX8 is a efficacious marker to distinguish serous ovarian tumor [41, 42], metastatic müllerian carcinomas , renal cell carcinoma . The expression of PAX5 and PAX8 are different in lung cancer tissues in order to further identify potential therapeutic targets, while PAX5 inclined to expressed in SCLC cells, overexpression of PAX8 in most NSCLC cell lines . Moreover, the results of Muratovska et al. presented that PAX genes influence cell survival in specific types of cancer .
Previous studies demonstrated that AC016683.6 eQTLs SNPs might contribute to susceptibility or prognosis to different cancers [29, 30, 45]. Han et al. firstly conducted a case–control study, which revealed that rs4848320 TT and rs1110839 GG decreased risk of cervical cancer, with the corresponding adjusted ORs (95% CIs) of 0.58 (0.36–0.93) and 0.75 (0.58–0.97), respectively . However, another study found that rs4848320 polymorphism may also be a risk factor for the risk of childhood acute lymphoblastic leukemia . Moreover, Ma et al. demonstrated the SNPs in AC016683.6 were significantly associated with the prognosis of hepatocellular carcinoma . However, the study did not provide both rs4848320 and rs1110839 polymorphisms associated with lung cancer risk in overall population. We identified that rs4848320 and rs1110839 polymorphisms increased the risk of lung cancer in the smoking and male population. Besides, we founded that there was positive multiplicative interaction between rs4848320 risk genotypes and smoking exposure, indicating that the interaction increased the susceptibility to lung cancer.
Previous evidence identified that SNP mutations in lncRNAs may alter their structural stability, generate alternative splicing, impair the translation of target mRNAs, and ultimately affect the risk of various cancers . For example, rs35592567 polymorphism affect the expression of TP63 by interfering miR-140, which may serve as a reasonable explanation for the increased susceptibility of gastric cancer . Additionally, Xue et al. reported that rs7958904 G > C strikingly altered the secondary structure of HOTAIR, suggesting the SNPs affect the susceptibility to colorectal cancer. LncRNA AC016683.6 is located in the intron region of PAX8, which mapped on chromosome 2q13. Rs1110839 and rs4848320 polymorphisms in AC016683.6 may be the eQTLs SNPs for PAX8 gene. Consequently, it is likely that the two SNPs could influence the specific interaction between AC016683.6 and PAX8, directly regulating the expression of PAX8. Therefore, we need to explore the functional properties of SNPs in oncogenic lncRNAs, which may be beneficial for developing the potential of lncRNAs in the diagnosis and treatment of cancer.
The study should emphatically highlight several implicated limitations. Firstly, the subjects from Shenyang city may lack of representative to some extent. Secondly, sample size of the study is limited, especially stratified subgroup. Thirdly, functional verification of two SNPs in AC016683.6 did not in our current study. Therefore, more large-scale studies confirm these results of the study across different ethnicities in the future.
Rs4848320 and rs1110839 genetic variants may be associated with the increased risk of lung cancer. The rs4848320 risk genotypes combined with smoking status strengthened the risk effect on lung cancer risk. However, the rs4848320 and rs1110839 genetic variants in AC016683.6 may serve as potential genetic risk factors for lung cancer susceptibility, which needs to be verified in a larger population. Additionally, biological function of the two SNPs in AC016683.6 need further to clarify in lung cancer. Previous studies and the current study provide novel clues for functional analysis of cancer-related susceptibility loci.
long non-coding RNA
paired-box gene 8
single nucleotide polymorphism
non-small cell lung cancer
squamous cell carcinoma
small cell lung cancer
- 95% CI:
95% confidence interval
real-time polymerase chain reaction
Hong QY, Wu GM, Qian GS, Hu CP, Zhou JY, Chen LA, Li WM, Li SY, Wang K, Wang Q, et al. Prevention and management of lung cancer in China. Cancer. 2015;121(Suppl 17):3080–8.
Fang JY, Dong HL, Wu KS, Du PL, Xu ZX, Lin K. Characteristics and prediction of lung cancer mortality in China from 1991 to 2013. Asian Pac J Cancer Prev (APJCP). 2015;16(14):5829–34.
Zhou C. Lung cancer molecular epidemiology in China: recent trends. Transl Lung Cancer Res. 2014;3(5):270–9.
Zhu H, Lv Z, An C, Shi M, Pan W, Zhou L, Yang W, Yang M. Onco-lncRNA HOTAIR and its functional genetic variants in papillary thyroid carcinoma. Sci Rep. 2016;6:31969.
Orom UA, Derrien T, Beringer M, Gumireddy K, Gardini A, Bussotti G, Lai F, Zytnicki M, Notredame C, Huang Q, et al. Long noncoding RNAs with enhancer-like function in human cells. Cell. 2010;143(1):46–58.
Schmitt AM, Chang HY. Long Noncoding RNAs in cancer pathways. Cancer Cell. 2016;29(4):452–63.
Murugan AK, Munirajan AK, Alzahrani AS. Long noncoding RNAs: emerging players in thyroid cancer pathogenesis. Endocr Relat Cancer. 2018;25(2):R59–82.
Li J, Xuan Z, Liu C. Long non-coding RNAs and complex human diseases. Int J Mol Sci. 2013;14(9):18790–808.
Huarte M. The emerging role of lncRNAs in cancer. Nat Med. 2015;21(11):1253–61.
Zhang A, Xu M, Mo YY. Role of the lncRNA-p53 regulatory network in cancer. J Mol Cell Biol. 2014;6(3):181–91.
Tao H, Yang JJ, Zhou X, Deng ZY, Shi KH, Li J. Emerging role of long noncoding RNAs in lung cancer: current status and future prospects. Respir Med. 2016;110:12–9.
Lang D, Powell SK, Plummer RS, Young KP, Ruggeri BA. PAX genes: roles in development, pathophysiology, and cancer. Biochem Pharmacol. 2007;73(1):1–14.
Hewitt SM, Hamada S, Monarres A, Kottical LV, Saunders GF, McDonnell TJ. Transcriptional activation of the bcl-2 apoptosis suppressor gene by the paired box transcription factor PAX8. Anticancer Res. 1997;17(5a):3211–5.
Stuart ET, Haffner R, Oren M, Gruss P. Loss of p53 function through PAX-mediated transcriptional repression. EMBO J. 1995;14(22):5638–45.
Heidarpour M, Tavanafar Z. Diagnostic utility of PAX8 in differentiation of mullerian from non-mullerian tumors. Adv Biomed Res. 2014;3:96.
Park D, Jia H, Rajakumar V, Chamberlin HM. Pax2/5/8 proteins promote cell survival in C. elegans. Development. 2006;133(21):4193–202.
Muratovska A, Zhou C, He S, Goodyer P, Eccles MR. Paired-box genes are frequently expressed in cancer and often required for cancer cell survival. Oncogene. 2003;22(39):7989–97.
Buttiglieri S, Deregibus MC, Bravo S, Cassoni P, Chiarle R, Bussolati B, Camussi G. Role of Pax2 in apoptosis resistance and proinvasive phenotype of Kaposi’s sarcoma cells. J Biol Chem. 2004;279(6):4136–43.
Miccadei S, Provenzano C, Mojzisek M, Natali PG, Civitareale D. Retinoblastoma protein acts as Pax 8 transcriptional coactivator. Oncogene. 2005;24(47):6993–7001.
Romitti M, Ceolin L, Siqueira DR, Ferreira CV, Wajner SM, Maia AL. Signaling pathways in follicular cell-derived thyroid carcinomas (review). Int J Oncol. 2013;42(1):19–28.
Cheung L, Messina M, Gill A, Clarkson A, Learoyd D, Delbridge L, Wentworth J, Philips J, Clifton-Bligh R, Robinson BG. Detection of the PAX8-PPAR gamma fusion oncogene in both follicular thyroid carcinomas and adenomas. J Clin Endocrinol Metab. 2003;88(1):354–7.
Marques AR, Espadinha C, Catarino AL, Moniz S, Pereira T, Sobrinho LG, Leite V. Expression of PAX8-PPAR gamma 1 rearrangements in both follicular thyroid carcinomas and adenomas. J Clin Endocrinol Metab. 2002;87(8):3947–52.
Klemke M, Drieschner N, Belge G, Burchardt K, Junker K, Bullerdiek J. Detection of PAX8-PPARG fusion transcripts in archival thyroid carcinoma samples by conventional RT-PCR. Genes Chromosomes Cancer. 2012;51(4):402–8.
Nikiforova MN, Lynch RA, Biddinger PW, Alexander EK, Dorn GW 2nd, Tallini G, Kroll TG, Nikiforov YE. RAS point mutations and PAX8-PPAR gamma rearrangement in thyroid tumors: evidence for distinct molecular pathways in thyroid follicular carcinoma. J Clin Endocrinol Metab. 2003;88(5):2318–26.
Lui WO, Foukakis T, Liden J, Thoppe SR, Dwight T, Hoog A, Zedenius J, Wallin G, Reimers M, Larsson C. Expression profiling reveals a distinct transcription signature in follicular thyroid carcinomas with a PAX8-PPARγ fusion oncogene. Oncogene. 2005;24(8):1467–76.
Reddi HV, McIver B, Grebe SK, Eberhardt NL. The paired box-8/peroxisome proliferator-activated receptor-gamma oncogene in thyroid tumorigenesis. Endocrinology. 2007;148(3):932–5.
Kanteti R, Nallasura V, Loganathan S, Tretiakova M, Kroll T, Krishnaswamy S, Faoro L, Cagle P, Husain AN, Vokes EE, et al. PAX5 is expressed in small-cell lung cancer and positively regulates c-Met transcription. Lab Invest. 2009;89(3):301–14.
Kanteti R, El-Hashani E, Dhanasingh I, Tretiakova M, Husain AN, Sharma S, Sharma J, Vokes EE, Salgia R. Role of PAX8 in the regulation of MET and RON receptor tyrosine kinases in non-small cell lung cancer. BMC Cancer. 2014;14:185.
Han J, Zhou W, Jia M, Wen J, Jiang J, Shi J, Zhang K, Ma H, Liu J, Ren J, et al. Expression quantitative trait loci in long non-coding RNA PAX8-AS1 are associated with decreased risk of cervical cancer. Mol Genet Genom (MGG). 2016;291(4):1743–8.
Ma S, Yang J, Song C, Ge Z, Zhou J, Zhang G, Hu Z. Expression quantitative trait loci for PAX8 contributes to the prognosis of hepatocellular carcinoma. PLoS ONE. 2017;12(3):e0173700.
Boyle AP, Hong EL, Hariharan M, Cheng Y, Schaub MA, Kasowski M, Karczewski KJ, Park J, Hitz BC, Weng S, et al. Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 2012;22(9):1790–7.
Veyrieras JB, Kudaravalli S, Kim SY, Dermitzakis ET, Gilad Y, Stephens M, Pritchard JK. High-resolution mapping of expression-QTLs yields insight into human gene regulation. PLoS Genet. 2008;4(10):e1000214.
Stranger BE, Nica AC, Forrest MS, Dimas A, Bird CP, Beazley C, Ingle CE, Dunning M, Flicek P, Koller D, et al. Population genomics of human gene expression. Nat Genet. 2007;39(10):1217–24.
Yin Z, Cui Z, Ren Y, Xia L, Wang Q, Zhang Y, He Q, Zhou B. Association between polymorphisms in pre-miRNA genes and risk of lung cancer in a Chinese non-smoking female population. Lung Cancer. 2016;94:15–21.
Wang PH, Shen HB, Chen F, Zhao JK. Study on the significance and application of crossover analysis in assessing gene-environmental interaction. Zhonghua liu xing bing xue za zhi (Zhonghua liuxingbingxue zazhi). 2005;26(1):54–7.
Guttman M, Donaghey J, Carey BW, Garber M, Grenier JK, Munson G, Young G, Lucas AB, Ach R, Bruhn L, et al. lincRNAs act in the circuitry controlling pluripotency and differentiation. Nature. 2011;477(7364):295–300.
McGinnis W, Levine MS, Hafen E, Kuroiwa A, Gehring WJ. A conserved DNA sequence in homoeotic genes of the Drosophila Antennapedia and bithorax complexes. Nature. 1984;308(5958):428–33.
Fabian P, Kozmikova I, Kozmik Z, Pantzartzi CN. Pax2/5/8 and Pax6 alternative splicing events in basal chordates and vertebrates: a focus on paired box domain. Front Genet. 2015;6:228.
Li L, Yang Y, Xue L. Regulatory functions of Pax gene family in Drosophila development. Yi chuan (Hereditas). 2010;32(2):115–21.
Paixao-Cortes VR, Salzano FM, Bortolini MC. Origins and evolvability of the PAX family. Semin Cell Dev Biol. 2015;44:64–74.
Nonaka D, Chiriboga L, Soslow RA. Expression of pax8 as a useful marker in distinguishing ovarian carcinomas from mammary carcinomas. Am J Surg Pathol. 2008;32(10):1566–71.
Laury AR, Hornick JL, Perets R, Krane JF, Corson J, Drapkin R, Hirsch MS. PAX8 reliably distinguishes ovarian serous tumors from malignant mesothelioma. Am J Surg Pathol. 2010;34(5):627–35.
Waters L, Crumley S, Truong L, Mody D, Coffey D. PAX2 and PAX8: useful markers for metastatic effusions. Acta Cytol. 2014;58(1):60–6.
Tacha D, Zhou D, Cheng L. Expression of PAX8 in normal and neoplastic tissues: a comprehensive immunohistochemical study. Appl Immunohistochem Mol Morphol (AIMM). 2011;19(4):293–9.
Bahari G, Hashemi M, Naderi M, Sadeghi-Bojd S, Taheri M. Long non-coding RNA PAX8-AS1 polymorphisms increase the risk of childhood acute lymphoblastic leukemia. Biomed Rep. 2018;8(2):184–90.
Lu C, Yang Y, Ma S. A functional variant (Rs35592567) in TP63 at 3q28 is associated with gastric cancer risk via modifying its regulation by MicroRNA-140. Cell Physiol Biochem. 2018;47(1):235–44.
JL and ZHY conceived and wrote the paper; JL and XTL gathered and analyzed the data; HL, MG, YHB, ZTY revised the whole paper. All authors have reviewed the final version of the manuscript and approved to submit to your journal. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The data of the study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
The study has acquired approval of the Institutional Review Board of China Medical University, and all enrolled participants or their representatives signed the informed consent according to relevant regulations. All participants signed informed consent in the study.
This study was funded by National Natural Science Foundation of China (No. 81673261).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.