Association of Long Non-Coding RNAs (lncRNAs) ANRIL and MALAT1 Polymorphism with Cervical Cancer.

Background: Long non-coding RNAs (lncRNAs) and their polymorphisms play crucial roles in the development of different cancers.
Methods: Eight single-nucleotide polymorphisms (SNPs) in ANRIL and MALAT1 (rs1333045, rs4977574, rs1333048, and rs10757278 in ANRIL and rs11227209, rs619586, rs664589, and rs3200401 in MALAT1) were enrolled and genotyped in a total of 1248 samples, including 587 patients with cervical cancer (CC) and 661 healthy individuals using in TaqMan assay. The association of these SNPs with CC was then evaluated.
Results: Our results showed that the allele and genotype frequencies of rs3200401 in MALAT1 were significantly different between the control and CC groups after Bonferroni correction (P = 0.001 and P = 0.004, respectively), indicating that the C allele is a protective factor against CC (OR = 0.70; 95% CI = 0.57-0.87). In addition, the allele and genotype frequencies of rs4977574 in ANRIL were significantly different between the control and CC groups after Bonferroni correction (P = 0.004 and P = 0.014, respectively), and the A allele might be a protective factor for CC (OR = 0.80; 95% CI = 0.68-0.93). For subgroup analysis, the alleles of rs3200401 in MALAT1 showed significant differences between the control and adenocarcinoma (AC) and control and squamous cell carcinoma (SCC) groups (P = 0.005 and P = 0.004, respectively). The rs3200401C allele could be a protective factor for AC and SCC development (OR = 0.57; 95% CI = 0.38-0.85; OR = 0.72; 95% CI = 0.58-0.90). Moreover, the rs3200401C allele could be a protective factor for cervical cancer stage I development (OR = 0.67; 95% CI = 0.53-0.86).
Conclusion: Our results indicate that rs3200401 in MALAT1 and rs4977574 in ANRIL could play key roles in the CC development.

Introduction

Cervical cancer (CC) is the fourth common cancer among women. China has contributed to an increasing global cervical cancer burden, with 106,000 cases and 48,000 deaths.1 Cervical cancer includes two main pathological types: squamous cell carcinoma (SCC) and adenocarcinoma (AC). SCC is the most common pathological type which covers 75% to 80% of cases, followed by AC, which accounts for 10.0% to 25.0%.2

Recently, a genome-wide association study found that host genetic factors could play an important role in the CC development.3–7 Many studies have focused on the dysregulation of long non-coding RNAs (lncRNAs) in host genes, because they play critical regulatory roles in cancers as either tumour suppressors or oncogenic lncRNAs.8,9 lncRNA antisense non-coding RNA in the INK4 locus (ANRIL) is located in the chromosome 9p21, which is associated with different cancers.10–15 Besides the ANRIL, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), located on chromosome 11q13.1, participates in epigenetic changes, changes in gene expression, and alternative splicing, which is also associated with different cancers.16,17

Many studies have reported that the polymorphisms of lncRNAs could influence the interaction between lncRNAs and other molecules and are associated with cancer development.18–20 In 2017, Taheri et al investigated the association between four ANRIL SNPs (rs1333045, rs4977574, rs1333048, and rs10757278) and benign prostate hyperplasia and prostate cancer.21 Their results showed that rs4977574, rs1333048, and rs10757278 were associated with benign prostate hyperplasia and prostate cancer risk.21 In 2019, Qu et al investigated the association between MALAT1 tagSNPs (rs11227209, rs619586, rs664589, and rs3200401) and oesophageal squamous cell carcinoma (ESCC).22 Their results showed that rs3200401C is associated with an increased risk of ESCC.22

In the current study, we genotyped four SNPs (rs1333045, rs4977574, rs1333048, and rs10757278) in ANRIL and four tagSNPs (rs11227209, rs619586, rs664589, and rs3200401) in MALAT1 in patients with CC and healthy individuals. We then analyzed the association of these eight SNPs with CC development.

Materials and Methods

Ethics Declarations

This study was approved by the Institutional Review Board of the Third Affiliated Hospital of Kunming Medical University. The study protocol was in accordance with the principles of the Helsinki Declaration of 1975, as revised in 2008.

Subjects

A total of 1248 subjects, including 661 healthy control and 587 patients with CC were recruited for this study. All patients were diagnosed at the Third Affiliated Hospital of Kunming Medical University between 2013 and 2018. The diagnosis of CC was based on the World Health Organization Comprehensive Cervical Cancer Control: A Guide to Essential Practice and the International Federation of Gynaecology and Obstetrics. The inclusion criteria for healthy control was as follows: (1) lack of cancer history and lesion; (2) HPV negative. All subjects were Han Chinese from Yunnan Province, and all signed informed consent forms.

DNA Extraction

Venous blood (5 mL) was collected from each subject, all of whom had fasted. Genomic DNA from peripheral blood was extracted using genomic DNA mini kit (QIAamp DNA Blood Mini Kit). An ultramicro UV-visible spectrophotometer (ND-2000, Thermo Scientific, USA) was used to detect the concentration and purity of the DNA.

SNP Genotyping

The probes and primers used for genotyping these eight SNPs were all purchased from ABI (). Primers and probes for the genotyping were commercially available. The array ID were C__32062431_10 (rs11227209), C__1060479_10 (rs619586), C__1060482_20 (rs664589), C__3246069_10 (rs3200401), C__1754667_10 (rs1333048), C__1754681_10 (rs4977574), C__8766826_10 (rs1333045), C__11841860_10 (rs10757278). The TaqMan Genotyping Master Mix used in the genotyping was purchased from ABI, and the genotyping process was the same as that used in our previous study.23 The PCR experiment data were analyzed using TaqMan Genotyper Software (Version 1.3.1). Several samples of each genotype were sequenced to identify the accuracy of SNP genotyping. Each genotype sample was then used as a positive sample in TaqMan assay.

Statistical Analysis

Differences in age between the CC and control groups were analyzed using the Student’s t-test. Hardy-Weinberg equilibrium (HWE) was calculated using SHEsis program.24,25 The associations of eight SNPs with CC were calculated by SHEsis program.24,25 The inheritance model analysis is an attempt to identify the mode of inheritance of association between eight SNPs genotypes and CC which was used by SNPstats software.26 The inheritance models included codominant, dominant, recessive, overdominant, and log-additive models. The best fit model for each SNP was evaluated using the Akaike information criterion (AIC) and the Bayesian information criterion (BIC), and the smallest AIC and BIC were the best fit models.26 Bonferroni correction was applied for multiple comparisons. The differences were statistically significant at P<0.006 (0.05/8).

Results

Subject Characteristics

Table 1 shows the clinical data of the participants in this study. The CC group contained 587 patients, including 488 with squamous cell carcinoma (SCC), 76 with adenocarcinoma (AC), and 23 with other types of cancer. According to CC staging, 349 patients were in stage I, 206 in stage II, and 32 in stage III+IV. The average ages for CC and control groups were 46.12±9.34 and 46.80±6.91 respectively. No significant difference in age was found between the two groups (P=0.153).

Table 1

Characteristics of the Subjects Enrolled in the Current Study

		CC	Control	t	P value
N		587	661	1.430	0.153
Ages (year)		46.12±9.34	46.80±6.91
Pathologic types	SCC (n)	488
	AC(n)	76
	Others (n)	23
Stages of CC	I (n)	349
	II (n)	206
	III and IV (n)	32

Abbreviations: CC, cervical cancer; SCC, squamous cell carcinoma; AC, adenocarcinoma.

Association Analysis of Eight SNPs in the MALAT1 and ANRIL in Control and CC Groups

The alleles and genotype frequencies of the eight SNPs in MALAT1 and ANRIL between the control and CC groups are shown in Table 2. The results of HWE show that the four SNPs in MALAT1 and four SNPs in ANRIL were all in HWE in both the control and CC groups (P>0.05) (Table 2). For MALAT1, the allele and genotype frequencies of rs3200401 were significantly different between the control and CC groups after Bonferroni correction (P=0.001 and P=0.004, respectively), and the C allele might be a protective factor against CC (OR=0.70; 95% CI=0.57–0.87). For ANRIL, the allele and genotype frequencies of rs4977574 were significantly different between the control and CC groups after Bonferroni correction (P=0.004 and P=0.014, respectively), and the A allele might be a protective factor for CC (OR=0.80; 95% CI=0.68–0.93). For the other six SNPs, there were no difference in alleles and genotypes between these two groups after Bonferroni correction (P>0.006).

Table 2

The Allele and Genotype Distribution of the Eight SNPs in Control and Cervical Cancer Groups

SNPs	Alleles/Genotypes	Control n (%)	CC n (%)	OR [95% CI]	P value	HWE (P value)
						Control	CC
rs11227209	C	1257(95.1)	1115(95.0)	0.98 [0.68–1.40]	0.901	0.738	0.199
	G	65(4.9)	59(5.0)
	C/C	598(90.5)	528(89.9)		0.367
	C/G	61(9.2)	59(10.1)
	G/G	2(0.3)	0(0.0)
rs619586	A	1210(91.5)	1069(91.1)	0.94[0.71–1.25]	0.676	0.529	0.094
	G	112(8.5)	105(8.9)
	A/A	555(84.0)	490(83.5)		0.747
	A/G	100(15.1)	89(15.2)
	G/G	6(0.9)	8(1.3)
rs664589	C	1218(92.1)	1086(92.5)	1.05[0.78–1.42]	0.728	0.118	0.439
	G	104(7.9)	88(7.5)
	C/C	564(85.3)	501(85.3)		0.312
	C/G	90(13.6)	84(14.3)
	G/G	7(1.1)	2(0.3)
rs3200401	C	1131(85.6)	946(80.6)	0.70[0.57–0.87]	0.001	0.949	0.451
	T	191(14.4)	228(19.4)
	C/C	484(73.2)	384(65.4)		0.004
	C/T	163(24.7)	178(30.3)
	T/T	14(2.1)	25(4.3)
rs1333048	A	673(50.9)	556(47.4)	0.87[0.74–1.02]	0.077	0.109	0.660
	C	649(49.1)	618(52.6)
	A/A	161(24.4)	129(22.0)		0.144
	A/C	351(53.1)	298(50.8)
	C/C	149(22.5)	160(27.3)
rs4977574	A	713(53.9)	566(48.2)	0.80[0.68–0.93]	0.004	0.368	0.195
	G	609(46.1)	608(51.8)
	A/A	184(27.8)	131(22.3)		0.014
	A/G	345(52.2)	304(51.8)
	G/G	132(20.0)	152(25.9)
rs1333045	T	665(50.3)	546(46.5)	1.16[1.00–1.36]	0.058	0.080	0.410
	C	657(49.7)	628(53.5)
	T/T	156(23.6)	122(20.8)		0.126
	T/C	353(53.4)	302(51.4)
	C/C	152(23.0)	163(27.8)
rs10757278	A	681(51.5)	567(48.3)	0.88[0.75–1.03]	0.109	0.076	0.989
	G	641(48.5)	607(51.7)
	A/A	164(24.8)	137(23.3)		0.124
	A/G	353(53.4)	293(49.9)
	G/G	144(21.8)	157(26.7)

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviations: CC, cervical cancer; HWE, Hardy-Weinberg equilibrium.

Inheritance Model Analysis of Eight SNPs in the MALAT1 and ANRIL in Control and CC Groups

The inheritance models were used to analyze the association of these SNP genotypes with CC. The results of the inheritance model analysis of eight SNPs in MALAT1 and ANRIL between the control and CC groups are shown in Tables 3 and 4. Comparison between CC and control groups showed that the log-additive model was the best fit model for rs3200401 in MALAT1. In this model, the 2T/T+C/T genotype was associated with an increased CC risk (P=0.001, OR=1.42; 95% CI=1.15–1.75). The best fit for rs4977574 in ANRIL was the log-additive inheritance model. In this model, the 2G/G+A/G genotype was associated with an increased CC risk (P=0.003, OR=1.28; 95% CI=1.09–1.50). There were no significant differences of other SNPs between the control and CC groups in the inheritance models (P>0.006).

Table 3

The Inheritance Model Analysis for the Four SNPs in MALAT1 Between CC and Control Groups

SNPs	Models	Genotypes	Control n (%)	CC n (%)	OR (95% CI)	P value	AIC	BIC
rs11227209	Codominant	C/C	598 (90.5)	528 (90.0)	1	0.260	1728.7	1749.3
		C/G	61 (9.2)	59 (10.1)	1.10 (0.76–1.61)
		G/G	2 (0.3)	0 (0.0)	0.00 (0.00-NA)
	Dominant	C/C	598 (90.5)	528 (90.0)	1	0.730	1729.3	1744.7
		C/G-G/G	63 (9.5)	59 (10.1)	1.07 (0.73–1.55)
	Recessive	C/C-C/G	659 (99.7)	587 (100.0)	1	0.120	1727.0	1742.4
		G/G	2 (0.3)	0 (0.0)	0.00 (0.00-NA)
	Overdominant	C/C-G/G	600 (90.8)	528 (90.0)	1	0.600	1729.2	1744.6
		C/G	61 (9.2)	59 (10.1)	1.11 (0.76–1.61)
	Log-additive				1.03 (0.72–1.49)	0.870	1729.4	1744.8
rs619586	Codominant	A/A	555 (84)	490 (83.5)	1	0.740	1730.8	1751.3
		G/A	100 (15.1)	89 (15.2)	1.01 (0.74–1.37)
		G/G	6 (0.9)	8 (1.4)	1.52 (0.52–4.43)
	Dominant	A/A	555 (84.0)	490 (83.5)	1	0.810	1729.4	1744.8
		G/A-G/G	106 (16.0)	97 (16.5)	1.04 (0.77–1.40)
	Recessive	A/A-G/A	655 (99.1)	579 (98.6)	1	0.440	1728.8	1744.2
		G/G	6 (0.9)	8 (1.4)	1.52 (0.52–4.42)
	Overdominant	A/A-G/G	561 (84.9)	498 (84.8)	1	0.990	1729.4	1744.8
		G/A	100 (15.1)	89 (15.2)	1.00 (0.73–1.37)
	Log-additive				1.06 (0.81–1.39)	0.680	1729.3	1744.7
rs664589	Codominant	C/C	564 (85.3)	501 (85.3)	1	0.290	1729.0	1749.5
		C/G	90 (13.6)	84 (14.3)	1.05 (0.76–1.45)
		G/G	7 (1.1)	2 (0.3)	0.32 (0.07–1.55)
	Dominant	C/C	564 (85.3)	501 (85.3)	1	0.980	1729.4	1744.8
		C/G-G/G	97 (14.7)	86 (14.7)	1.00 (0.73–1.36)
	Recessive	C/C-C/G	654 (98.9)	585 (99.7)	1	0.120	1727.0	1742.4
		G/G	7 (1.1)	2 (0.3)	0.32 (0.07–1.54)
	Overdominant	C/C-G/G	571 (86.4)	503 (85.7)	1	0.730	1729.3	1744.7
		C/G	90 (13.6)	84 (14.3)	1.06 (0.77–1.46)
	Log-additive				0.95 (0.71–1.27)	0.720	1729.3	1744.7
rs3200401	Codominant	C/C	484 (73.2)	384 (65.4)	1	0.004	1720.5	1741.0
		C/T	163 (24.7)	178 (30.3)	1.37 (1.07–1.77)
		T/T	14 (2.1)	25 (4.3)	2.26 (1.16–4.42)
	Dominant	C/C	484 (73.2)	384 (65.4)	1	0.003	1720.6	1736.0
		C/T-T/T	177 (26.8)	203 (34.6)	1.44 (1.13–1.84)
	Recessive	C/C-C/T	647 (97.9)	562 (95.7)	1	0.028	1724.6	1740.0
		T/T	14 (2.1)	25 (4.3)	2.07 (1.06–4.02)
	Overdominant	C/C-T/T	498 (75.3)	409 (69.7)	1	0.026	1724.5	1739.9
		C/T	163 (24.7)	178 (30.3)	1.33 (1.03–1.70)
	Log-additive				1.42 (1.15–1.75)	0.001	1718.7	1734.1

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviation: CC, cervical cancer.

Table 4

The Inheritance Model Analysis for the Four SNPs in ANRIL Between CC and Control Groups

SNPs	Models	Genotypes	Control n (%)	CC n (%)	OR (95% CI)	P value	AIC	BIC
rs1333048	Codominant	C/C	149 (22.5)	160 (27.3)	1	0.130	1727.3	1747.8
		A/C	351 (53.1)	298 (50.8)	0.78 (0.60–1.03)
		A/A	161 (24.4)	129 (22)	0.74 (0.54–1.02)
	Dominant	C/C	149 (22.5)	160 (27.3)	1	0.046	1725.5	1740.9
		A/C-A/A	512 (77.5)	427 (72.7)	0.77 (0.59–1.00)
	Recessive	C/C-A/C	500 (75.6)	458 (78.0)	1	0.320	1728.5	1743.8
		A/A	161 (24.4)	129 (22.0)	0.88 (0.67–1.14)
	Overdominant	C/C-A/A	310 (46.9)	289 (49.2)	1	0.380	1728.7	1744.0
		A/C	351 (53.1)	298 (50.8)	0.90 (0.72–1.13)
	Log-additive				0.86 (0.73–1.01)	0.065	1726.0	1741.4
rs4977574	Codominant	A/A	184 (27.8)	131 (22.3)	1	0.012	1722.6	1743.1
		A/G	345 (52.2)	304 (51.8)	1.24 (0.94–1.62)
		G/G	132 (20.0)	152 (25.9)	1.63 (1.18–2.26)
	Dominant	A/A	184 (27.8)	131 (22.3)	1	0.024	1724.3	1739.7
		A/G-G/G	477 (72.2)	456 (77.7)	1.35 (1.04–1.74)
	Recessive	A/A-A/G	529 (80.0)	435 (74.1)	1	0.011	1722.9	1738.3
		G/G	132 (20.0)	152 (25.9)	1.41 (1.08–1.84)
	Overdominant	A/A-G/G	316 (47.8)	283 (48.2)	1	0.860	1729.4	1744.8
		A/G	345 (52.2)	304 (51.8)	0.98 (0.78–1.22)
	Log-additive				1.28 (1.09–1.50)	0.003	1720.6	1736.0
rs1333045	Codominant	C/C	152 (23.0)	163 (27.8)	1	0.110	1727.1	1747.6
		C/T	353 (53.4)	302 (51.5)	0.79 (0.60–1.03)
		T/T	156 (23.6)	122 (20.8)	0.72 (0.52–1.00)
	Dominant	C/C	152 (23.0)	163 (27.8)	1	0.045	1725.4	1740.8
		C/T-T/T	509 (77.0)	424 (72.2)	0.77 (0.60–0.99)
	Recessive	C/C-C/T	505 (76.4)	465 (79.2)	1	0.240	1728.0	1743.4
		T/T	156 (23.6)	122 (20.8)	0.85 (0.65–1.11)
	Overdominant	C/C-T/T	308 (46.6)	285 (48.5)	1	0.450	1728.9	1744.2
		C/T	353 (53.4)	302 (51.5)	0.92 (0.73–1.15)
	Log-additive				0.85 (0.72–1.00)	0.047	1725.5	1740.9
rs10757278	Codominant	A/A	164 (24.8)	137 (23.3)	1	0.100	1726.9	1747.4
		A/G	353 (53.4)	293 (49.9)	0.99 (0.75–1.30)
		G/G	144 (21.8)	157 (26.8)	1.32 (0.95–1.81)
	Dominant	A/A	164 (24.8)	137 (23.3)	1	0.550	1729.1	1744.5
		A/G-G/G	497 (75.2)	450 (76.7)	1.08 (0.83–1.41)
	Recessive	A/A-A/G	517 (78.2)	430 (73.2)	1	0.034	1724.9	1740.3
		G/G	144 (21.8)	157 (26.8)	1.33 (1.02–1.72)
	Overdominant	A/A-G/G	308 (46.6)	294 (50.1)	1	0.190	1727.7	1743.1
		A/G	353 (53.4)	293 (49.9)	0.86 (0.69–1.08)
	Log-additive				1.15 (0.98–1.35)	0.093	1726.6	1742.0

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviation: CC, cervical cancer.

Association Analysis of Eight SNPs in the MALAT1 and ANRIL with Different CC Pathological Types

A comparison between patients with different CC pathological types is shown in Tables 5 and 6. However, the alleles and genotypes of rs3200401 in MALAT1 were significantly different between the control and AC groups (P=0.005 and P=0.005, respectively). The rs3200401C allele could be a protective factor for AC development (OR=0.57; 95% CI=0.38–0.85). In addition, the alleles of rs3200401 in MALAT1 were significantly different between the control and SCC groups (P=0.004). The rs3200401C allele could be a protective factor for SCC development (OR=0.72; 95% CI=0.58–0.90). The allelic frequency and genotype frequency of eight SNPs in MALAT1 and ANRIL were not significantly different between the SCC and AC groups (P>0.006).

Table 5

The Allele and Genotype Distribution of the Eight SNPs in Control and SCC Groups

SNPs	Alleles/Genotypes	Control n (%)	SCC n (%)	OR [95% CI]	P value
rs11227209	C	1257(95.1)	929(95.2)	1.02[0.70–1.50]	0.911
	G	65(4.9)	47(4.8)
	C/C	598(90.5)	441(90.4)		0.466
	C/G	61(9.2)	47(9.6)
	G/G	2(0.3)	0(0.0)
rs619586	A	1210(91.5)	889(91.1)	0.95[0.71–1.27]	0.710
	G	112(8.5)	87(8.9)
	A/A	555(84.0)	407(83.4)		0.861
	A/G	100(15.1)	75(15.4)
	G/G	6(0.9)	6(1.2)
rs664589	C	1218(92.1)	905(92.7)	1.09[0.80–1.49]	0.597
	G	104(7.9)	71(7.3)
	C/C	564(85.3)	417(85.5)		0.069
	C/G	90(13.6)	71(14.5)
	G/G	7(1.1)	0(0.0)
rs3200401	C	1131(85.6)	791(81.0)	0.72[0.58–0.90]	0.004
	T	191(14.4)	185(19.0)
	C/C	484(73.2)	322(66.0)		0.069
	C/T	163(24.7)	147(30.1)
	T/T	14(2.1)	19(3.9)
rs1333048	A	673(50.9)	467(47.8)	0.89[0.750–1.04]	0.147
	C	649(49.1)	509(52.2)
	A/A	161(24.4)	109(22.3)		0.264
	A/C	351(53.1)	249(51.0)
	C/C	149(22.5)	130(26.6)
rs4977574	A	713(53.9)	475(48.7)	0.81[0.69–0.96]	0.013
	G	609(46.1)	501(51.3)
	A/A	184(27.8)	110(22.5)		0.038
	A/G	345(52.2)	255(52.3)
	G/G	132(20.0)	123(25.2)
rs1333045	T	665(50.3)	460(47.1)	0.88[0.75–1.04]	0.133
	C	657(49.7)	516(52.9)
	T/T	156(23.6)	105(21.5)		0.241
	T/C	353(53.4)	250(51.2)
	C/C	152(23.0)	133(27.3)
rs10757278	A	681(51.5)	477(48.9)	0.90[0.76–1.06]	0.211
	G	641(48.5)	499(51.1)
	A/A	164(24.8)	116(23.8)		0.245
	A/G	353(53.4)	245(50.2)
	G/G	144(21.8)	127(26.0)

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviation: SCC, squamous cell carcinoma.

Table 6

The Allele and Genotype Distribution of the Eight SNPs in Control and AC Groups

SNPs	Alleles/Genotypes	Control n (%)	AC n (%)	OR [95% CI]	P value
rs11227209	C	1257(95.1)	145(95.4)	1.07[0.48–2.38]	0.866
	G	65(4.9)	7(4.6)
	C/C	598(90.5)	69(90.8)		0.891
	C/G	61(9.2)	7(9.2)
	G/G	2(0.3)	0(0.0)
rs619586	A	1210(91.5)	142(93.4)	0.95[0.705–1.27]	0.710
	G	112(8.5)	10(6.6)
	A/A	555(84.0)	66(86.8)		0.627
	A/G	100(15.1)	10(13.2)
	G/G	6(0.9)	0(0.0)
rs664589	C	1218(92.1)	141(92.8)	1.09[0.57–2.09]	0.784
	G	104(7.9)	11(7.2)
	C/C	564(85.3)	66(86.8)		0.896
	C/G	90(13.6)	9(11.8)
	G/G	7(1.1)	1(1.3)
rs3200401	C	1131(85.6)	117(77.0)	0.57[0.38–0.85]	0.005
	T	191(14.4)	35(23.0)
	C/C	484(73.2)	47(61.8)		0.005
	C/T	163(24.7)	23(30.3)
	T/T	14(2.1)	6(7.9)
rs1333048	A	673(50.9)	67(44.1)	0.76[0.54–1.07]	0.111
	C	649(49.1)	85(55.9)
	A/A	161(24.4)	14(18.4)		0.247
	A/C	351(53.1)	39(51.3)
	C/C	149(22.5)	23(30.3)
rs4977574	A	713(53.9)	68(44.7)	0.69[0.49–0.97]	0.031
	G	609(46.1)	84(55.3)
	A/A	184(27.8)	15(19.7)		0.076
	A/G	345(52.2)	38(50.0)
	G/G	132(20.0)	23(30.3)
rs1333045	T	665(50.3)	66(43.4)	0.76[0.54–1.06]	0.108
	C	657(49.7)	86(56.6)
	T/T	156(23.6)	13(17.1)		0.249
	T/C	353(53.4)	40(52.6)
	C/C	152(23.0)	23(30.3)
rs10757278	A	681(51.5)	68(44.7)	0.76[0.54–1.07]	0.114
	G	641(48.5)	84(55.3)
	A/A	164(24.8)	15(19.7)		0.219
	A/G	353(53.4)	38(50.0)
	G/G	144(21.8)	23(30.3)

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviation: AC, adenocarcinoma.

Inheritance Model Analysis of Eight SNPs in the MALAT1 and ANRIL with Different CC Pathological Types

The results of the analysis of different pathological types of cervical cancer under different inheritance models are shown in Tables 7 and 8. Comparison between the control and SCC groups showed that the log-additive model was the best fit model for rs3200401 in MALAT1. In this model, the 2T/T+C/T genotype was associated with an increased risk of cervical cancer (P=0.004, OR=1.38; 95% CI=1.11–1.73). There were no significant differences between the AC and control groups, and between AC and SCC groups of the eight SNPs in the inheritance model (P>0.006).

Table 7

The Inheritance Model Analysis Between Control and SCC Groups

SNPs	Model	Genotype	Control n (%)	SCC n (%)	OR (95% CI)	P value	AIC	BIC
rs11227209	Codominant	C/C	598 (90.5)	441 (90.4)	1	0.330	1571.5	1591.7
		C/G	61 (9.2)	47 (9.6)	1.05 (0.70–1.57)
		G/G	2 (0.3)	0 (0.0)	0.00 (0.00-NA)
	Dominant	C/C	598 (90.5)	441 (90.4)	1	0.940	1571.7	1586.8
		C/G-G/G	63 (9.5)	47 (9.6)	1.02 (0.68–1.51)
	Recessive	C/C-C/G	659 (99.7)	488 (100.0)	1	0.140	1569.6	1584.7
		G/G	2 (0.3)	0 (0.0)	0.00 (0.00-NA)
	Overdominant	C/C-G/G	600 (90.8)	441 (90.4)	1	0.800	1571.6	1586.8
		C/G	61 (9.2)	47 (9.6)	1.05 (0.71–1.57)
	Log-additive				0.98 (0.67–1.45)	0.930	1571.7	1586.8
rs619586	Codominant	A/A	555 (84.0)	407 (83.4)	1	0.860	1573.4	1593.6
		G/A	100 (15.1)	75 (15.4)	1.02 (0.74–1.42)
		G/G	6 (0.9)	6 (1.2)	1.37 (0.44–4.28)
	Dominant	A/A	555 (84.0)	407 (83.4)	1	0.800	1571.6	1586.8
		G/A-G/G	106 (16.0)	81 (16.6)	1.04 (0.76–1.43)
	Recessive	A/A-G/A	655 (99.1)	482 (98.8)	1	0.590	1571.4	1586.6
		G/G	6 (0.9)	6 (1.2)	1.37 (0.44–4.26)
	Overdominant	A/A-G/G	561 (84.9)	413 (84.6)	1	0.920	1571.7	1586.8
		G/A	100 (15.1)	75 (15.4)	1.02 (0.74–1.41)
	Log-additive				1.05 (0.79–1.41)	0.720	1571.6	1586.7
rs664589	Codominant	C/C	564 (85.3)	417 (85.5)	1	0.019	1565.8	1586.0
		C/G	90 (13.6)	71 (14.6)	1.06 (0.76–1.49)
		G/G	7 (1.1)	0 (0.0)	0.00 (0.00-NA)
	Dominant	C/C	564 (85.3)	417 (85.5)	1	0.940	1571.7	1586.8
		C/G-G/G	97 (14.7)	71 (14.6)	0.99 (0.71–1.37)
	Recessive	C/C-C/G	654 (98.9)	488 (100.0)	1	0.005	1563.9	1579.1
		G/G	7 (1.1)	0 (0.0)	0.00 (0.00-NA)
	Overdominant	C/C-G/G	571 (86.4)	417 (85.5)	1	0.670	1571.5	1586.7
		C/G	90 (13.6)	71 (14.6)	1.08 (0.77–1.51)
	Log-additive				0.92 (0.67–1.25)	0.580	1571.4	1586.5
rs3200401	Codominant	C/C	484 (73.2)	322 (66.0)	1	0.015	1565.3	1585.5
		T/C	163 (24.7)	147 (30.1)	1.36 (1.04–1.77)
		T/T	14 (2.1)	19 (3.9)	2.07 (1.02–4.18)
	Dominant	C/C	484 (73.2)	322 (66.0)	1	0.008	1564.6	1579.8
		T/C-T/T	177 (26.8)	166 (34.0)	1.41 (1.10–1.82)
	Recessive	C/C-T/C	647 (97.9)	469 (96.1)	1	0.072	1568.5	1583.6
		T/T	14 (2.1)	19 (3.9)	1.90 (0.94–3.82)
	Overdominant	C/C-T/T	498 (75.3)	341 (69.9)	1	0.039	1567.5	1582.6
		T/C	163 (24.7)	147 (30.1)	1.32 (1.01–1.71)
	Log-additive				1.38 (1.11–1.73)	0.004	1563.4	1578.5
rs1333048	Codominant	C/C	149 (22.5)	130 (26.6)	1	0.250	1570.9	1591.1
		A/C	351 (53.1)	249 (51.0)	0.81 (0.61–1.08)
		A/A	161 (24.4)	109 (22.3)	0.77 (0.55–1.08)
	Dominant	C/C	149 (22.5)	130 (26.6)	1	0.100	1569.0	1584.2
		A/C-A/A	512 (77.5)	358 (73.4)	0.80 (0.61–1.05)
	Recessive	C/C-A/C	500 (75.6)	379 (77.7)	1	0.430	1571.1	1586.2
		A/A	161 (24.4)	109 (22.3)	0.89 (0.68–1.18)
	Overdominant	C/C-A/A	310 (46.9)	239 (49.0)	1	0.460	1571.2	1586.3
		A/C	351 (53.1)	249 (51.0)	0.92 (0.72–1.16)
	Log-additive				0.88 (0.74–1.04)	0.130	1569.4	1584.6
rs4977574	Codominant	A/A	184 (27.8)	110 (22.5)	1	0.035	1567.0	1587.2
		A/G	345 (52.2)	255 (52.2)	1.24 (0.93–1.65)
		G/G	132 (20.0)	123 (25.2)	1.57 (1.11–2.20)
	Dominant	A/A	184 (27.8)	110 (22.5)	1	0.041	1567.5	1582.7
		A/G-G/G	477 (72.2)	378 (77.5)	1.33 (1.01–1.74)
	Recessive	A/A-A/G	529 (80.0)	365 (74.8)	1	0.032	1567.1	1582.3
		G/G	132 (20.0)	123 (25.2)	1.36 (1.03–1.80)
	Overdominant	A/A-G/G	316 (47.8)	233 (47.8)	1	1.000	1571.7	1586.8
		A/G	345 (52.2)	255 (52.2)	1.00 (0.79–1.26)
	Log-additive				1.25 (1.05–1.48)	0.010	1565.0	1580.2
rs1333045	Codominant	C/C	152 (23.0)	133 (27.2)	1	0.230	1570.7	1590.9
		C/T	353 (53.4)	250 (51.2)	0.80 (0.60–1.07)
		T/T	156 (23.6)	105 (21.5)	0.77 (0.55–1.08)
	Dominant	C/C	152 (23.0)	133 (27.2)	1	0.090	1568.8	1584.0
		C/T-T/T	509 (77.0)	355 (72.8)	0.79 (0.60–1.04)
	Recessive	C/C-C/T	505 (76.4)	383 (78.5)	1	0.410	1571.0	1586.2
		T/T	156 (23.6)	105 (21.5)	0.89 (0.67–1.18)
	Overdominant	C/C-T/T	308 (46.6)	238 (48.8)	1	0.440	1571.1	1586.2
		C/T	353 (53.4)	250 (51.2)	0.91 (0.72–1.15)
	Log-additive				0.87 (0.74–1.03)	0.120	1569.2	1584.4
rs10757278	Codominant	A/A	164 (24.8)	116 (23.8)	1	0.220	1570.7	1590.9
		A/G	353 (53.4)	245 (50.2)	0.98 (0.73–1.31)
		G/G	144 (21.8)	127 (26.0)	1.25 (0.90–1.76)
	Dominant	A/A	164 (24.8)	116 (23.8)	1	0.690	1571.5	1586.7
		A/G-G/G	497 (75.2)	372 (76.2)	1.06 (0.80–1.39)
	Recessive	A/A-A/G	517 (78.2)	361 (74.0)	1	0.085	1568.7	1583.9
		G/G	144 (21.8)	127 (26.0)	1.27 (0.97–1.67)
	Overdominant	A/A-G/G	308 (46.6)	243 (49.8)	1	0.260	1570.5	1585.6
		A/G	353 (53.4)	245 (50.2)	0.87 (0.69–1.11)
	Log-additive				1.12 (0.95–1.33)	0.190	1570.0	1585.1

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviation: SCC, squamous cell carcinoma.

Table 8

The Inheritance Model Analysis Between Control and AC Groups

SNPs	Model	Genotype	Control n (%)	AC n (%)	OR (95% CI)	P value	AIC	BIC
rs11227209	Codominant	C/C	598 (90.5)	69 (90.8)	1	0.840	489.3	507.7
		C/G	61 (9.2)	7 (9.2)	0.97 (0.42–2.21)
		G/G	2 (0.3)	0 (0.0)	0.00 (0.00-NA)
	Dominant	C/C	598 (90.5)	69 (90.8)	1	0.890	487.6	501.4
		C/G-G/G	63 (9.5)	7 (9.2)	0.94 (0.41–2.15)
	Recessive	C/C-C/G	659 (99.7)	76 (100.0)	1	0.560	487.3	501.1
		G/G	2 (0.3)	0 (0.0)	0.00 (0.00-NA)
	Overdominant	C/C-G/G	600 (90.8)	69 (90.8)	1	0.940	487.6	501.4
		C/G	61 (9.2)	7 (9.2)	0.97 (0.42–2.21)
	Log-additive				0.92 (0.41–2.06)	0.840	487.6	501.4
rs619586	Codominant	A/A	555 (84.0)	66 (86.8)	1	0.460	488.0	506.5
		G/A	100 (15.1)	10 (13.2)	0.83 (0.41–1.68)
		G/G	6 (0.9)	0 (0.0)	0.00 (0.00-NA)
	Dominant	A/A	555 (84.0)	66 (86.8)	1	0.490	487.1	500.9
		G/A-G/G	106 (16.0)	10 (13.2)	0.79 (0.39–1.58)
	Recessive	A/A-G/A	655 (99.1)	76 (100.0)	1	0.260	486.3	500.1
		G/G	6 (0.9)	0 (0.0)	0.00 (0.00-NA)
	Overdominant	A/A-G/G	561 (84.9)	66 (86.8)	1	0.630	487.4	501.2
		G/A	100 (15.1)	10 (13.2)	0.84 (0.42–1.70)
	Log-additive				0.76 (0.39–1.48)	0.400	486.9	500.7
rs664589	Codominant	C/C	564 (85.3)	66 (86.8)	1	0.870	489.3	507.7
		C/G	90 (13.6)	9 (11.8)	0.84 (0.40–1.74)
		G/G	7 (1.1)	1 (1.3)	1.24 (0.15–10.27)
	Dominant	C/C	564 (85.3)	66 (86.8)	1	0.680	487.4	501.2
		C/G-G/G	97 (14.7)	10 (13.2)	0.86 (0.43–1.74)
	Recessive	C/C-C/G	654 (98.9)	75 (98.7)	1	0.830	487.6	501.4
		G/G	7 (1.1)	1 (1.3)	1.27 (0.15–10.49)
	Overdominant	C/C-G/G	571 (86.4)	67 (88.2)	1	0.620	487.3	501.2
		C/G	90 (13.6)	9 (11.8)	0.83 (0.40–1.74)
	Log-additive				0.90 (0.48–1.70)	0.750	487.5	501.3
rs3200401	Codominant	C/C	484 (73.2)	47 (61.8)	1	0.035	482.9	501.3
		T/C	163 (24.7)	23 (30.3)	1.42 (0.83–2.41)
		T/T	14 (2.1)	6 (7.9)	3.89 (1.41–10.68)
	Dominant	C/C	484 (73.2)	47 (61.8)	1	0.059	484.0	497.8
		T/C-T/T	177 (26.8)	29 (38.2)	1.63 (0.99–2.68)
	Recessive	C/C-T/C	647 (97.9)	70 (92.1)	1	0.024	482.5	496.3
		T/T	14 (2.1)	6 (7.9)	3.50 (1.29–9.49)
	Overdominant	C/C-T/T	498 (75.3)	53 (69.7)	1	0.330	486.7	500.5
		T/C	163 (24.7)	23 (30.3)	1.30 (0.77–2.19)
	Log-additive				1.67 (1.11–2.51)	0.016	481.8	495.6
rs1333048	Codominant	A/A	161 (24.4)	14 (18.4)	1	0.220	486.6	505.0
		A/C	351 (53.1)	39 (51.3)	1.26 (0.66–2.39)
		C/C	149 (22.5)	23 (30.3)	1.82 (0.90–3.69)
	Dominant	A/A	161 (24.4)	14 (18.4)	1	0.240	486.2	500.0
		A/C-C/C	500 (75.6)	62 (81.6)	1.42 (0.77–2.62)
	Recessive	A/A-A/C	512 (77.5)	53 (69.7)	1	0.110	485.1	498.9
		C/C	149 (22.5)	23 (30.3)	1.55 (0.91–2.62)
	Overdominant	A/A-C/C	310 (46.9)	37 (48.7)	1	0.690	487.4	501.2
		A/C	351 (53.1)	39 (51.3)	0.91 (0.56–1.46)
	Log-additive				1.36 (0.96–1.94)	0.085	484.6	498.4
rs4977574	Codominant	A/A	184 (27.8)	15 (19.7)	1	0.070	484.3	502.7
		A/G	345 (52.2)	38 (50.0)	1.36 (0.73–2.54)
		G/G	132 (20.0)	23 (30.3)	2.21 (1.11–4.41)
	Dominant	A/A	184 (27.8)	15 (19.7)	1	0.110	485.1	498.9
		A/G-G/G	477 (72.2)	61 (80.3)	1.59 (0.88–2.88)
	Recessive	A/A-A/G	529 (80.0)	53 (69.7)	1	0.036	483.2	497.0
		G/G	132 (20.0)	23 (30.3)	1.79 (1.05–3.04)
	Overdominant	A/A-G/G	316 (47.8)	38 (50.0)	1	0.690	487.4	501.2
		A/G	345 (52.2)	38 (50.0)	0.91 (0.56–1.46)
	Log-additive				1.50 (1.06–2.13)	0.023	482.4	496.2
rs1333045	Codominant	C/C	152 (23.0)	23 (30.3)	1	0.220	486.6	505.0
		C/T	353 (53.4)	40 (52.6)	0.72 (0.42–1.26)
		T/T	156 (23.6)	13 (17.1)	0.54 (0.26–1.10)
	Dominant	C/C	152 (23.0)	23 (30.3)	1	0.140	485.4	499.2
		C/T-T/T	509 (77.0)	53 (69.7)	0.67 (0.39–1.13)
	Recessive	C/C-C/T	505 (76.4)	63 (82.9)	1	0.190	485.9	499.7
		T/T	156 (23.6)	13 (17.1)	0.67 (0.36–1.25)
	Overdominant	C/C-T/T	308 (46.6)	36 (47.4)	1	0.830	487.5	501.4
		C/T	353 (53.4)	40 (52.6)	0.95 (0.59–1.53)
	Log-additive				0.73 (0.51–1.04)	0.083	484.6	498.4
rs10757278	Codominant	A/A	164 (24.8)	15 (19.7)	1	0.180	486.2	504.6
		A/G	353 (53.4)	38 (50.0)	1.16 (0.62–2.17)
		G/G	144 (21.8)	23 (30.3)	1.81 (0.91–3.62)
	Dominant	A/A	164 (24.8)	15 (19.7)	1	0.320	486.6	500.4
		A/G-G/G	497 (75.2)	61 (80.3)	1.34 (0.74–2.43)
	Recessive	A/A-A/G	517 (78.2)	53 (69.7)	1	0.074	484.4	498.2
		G/G	144 (21.8)	23 (30.3)	1.64 (0.97–2.78)
	Overdominant	A/A-G/G	308 (46.6)	38 (50.0)	1	0.480	487.1	500.9
		A/G	353 (53.4)	38 (50.0)	0.84 (0.52–1.36)
	Log-additive			–	1.37 (0.96–1.95)	0.083	484.6	498.4

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviation: AC, adenocarcinoma.

Association Analysis of Eight SNPs in the MALAT1 and ANRIL with Different Stages of Cervical Cancer

As only 32 patients had stage III + IV cervical cancer, we analyzed the association of the eight SNPs in MALAT1 and ANRIL for patients with stage I and stage II cervical cancer. The alleles and genotypes of rs3200401 in MALAT1 were significantly different between the control and stage I groups (P=0.001 and P=0.005, respectively) (Table 9). The rs3200401C allele could be a protective factor for stage I development (OR=0.67; 95% CI=0.53–0.86) (Table 9). The allelic and genotypic frequencies of the eight SNPs were not significantly different between the control and stage II, and stage I and stage II (P>0.006) (data not shown).

Table 9

The Allele and Genotype Distribution of the SNPs in Control and CC Stage I Groups

SNPs	Alleles/Genotypes	Control n (%)	CC I Stage n (%)	OR [95% CI]	P value
rs11227209	C	1257(95.1)	657(94.1)	0.83[0.55–1.24]	0.359
	G	65(4.9)	41(5.9)
	C/C	598(90.5)	308(88.3)		0.270
	C/G	61(9.2)	41(11.7)
	G/G	2(0.3)	0(0.0)
rs619586	A	1210(91.5)	631(90.4)	0.87[0.64–1.20]	0.397
	G	112(8.5)	67(9.6)
	A/A	555(84.0)	286(81.9)		0.704
	A/G	100(15.1)	59(16.9)
	G/G	6(0.9)	4(1.1)
rs664589	C	1218(92.1)	640(91.7)	0.94[0.67–1.32]	0.728
	G	104(7.9)	58(8.3)
	C/C	564(85.3)	292(83.7)		0.256
	C/G	90(13.6)	56(16.0)
	G/G	7(1.1)	1(0.3)
rs3200401	C	1131(85.6)	558(79.9)	0.67[0.53–0.86]	0.001
	T	191(14.4)	140(20.1)
	C/C	484(73.2)	224(64.2)		0.005
	C/T	163(24.7)	110(31.5)
	T/T	14(2.1)	15(4.3)
rs1333048	A	673(50.9)	330(47.3)	0.87[0.72–1.04]	0.121
	C	649(49.1)	368(52.7)
	A/A	161(24.4)	74(21.2)		0.269
	A/C	351(53.1)	182(52.1)
	C/C	149(22.5)	93(26.6)
rs4977574	A	713(53.9)	337(48.3)	0.80[0.66–0.96]	0.016
	G	609(46.1)	361(51.7)
	A/A	184(27.8)	75(21.5)		0.044
	A/G	345(52.2)	187(53.6)
	G/G	132(20.0)	87(24.9)
rs1333045	T	665(50.3)	324(46.4)	0.86[0.71–1.03]	0.097
	C	657(49.7)	374(53.6)
	T/T	156(23.6)	69(19.8)		0.228
	T/C	353(53.4)	186(53.3)
	C/C	152(23.0)	94(26.9)
rs10757278	A	681(51.5)	334(47.9)	0.86[0.72–1.04]	0.118
	G	641(48.5)	364(52.1)
	A/A	164(24.8)	77(22.1)		0.233
	A/G	353(53.4)	180(51.6)
	G/G	144(21.8)	92(26.4)

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviation: CC, cervical cancer.

Inheritance Model Analysis of Eight SNPs in the MALAT1 and ANRIL with Different Stages of Cervical Cancer

Comparison between the control and stage I groups showed that the log-additive model was the best fit model for rs3200401 in MALAT1 (Table 10). In this model, the 2T/T+C/T genotype was associated with an increased CC (P=0.002, OR=1.47; 95% CI=1.16–1.88) (Table 10). There were no significant differences between the control and stage II, and between stage I and stage II of the eight SNPs in the inheritance model (P>0.006) (data not shown).

Table 10

The Inheritance Model Analysis Between Control and CC Stage I Groups

SNPs	Model	Genotype	Control n (%)	CC Stage I n (%)	OR (95% CI)	P value	AIC	BIC
rs11227209	Codominant	C/C	598 (90.5)	308 (88.2)	1	0.210	1295.3	1315.0
		C/G	61 (9.2)	41 (11.8)	1.32 (0.87–2.02)
		G/G	2 (0.3)	0 (0.0)	0.00 (0.00-NA)
	Dominant	C/C	598 (90.5)	308 (88.2)	1	0.240	1295.1	1309.8
		C/G-G/G	63 (9.5)	41 (11.8)	1.29 (0.85–1.96)
	Recessive	C/C-C/G	659 (99.7)	349 (100.0)	1	0.220	1295.0	1309.7
		G/G	2 (0.3)	0 (0.0)	0.00 (0.00-NA)
	Overdominant	C/C-G/G	600 (90.8)	308 (88.2)	1	0.190	1294.7	1309.5
		C/G	61 (9.2)	41 (11.8)	1.33 (0.87–2.02)
	Log-additive				1.23 (0.82–1.86)	0.320	1295.4	1310.2
rs619586	Codominant	A/A	555 (84.0)	286 (82.0)	1	0.710	1297.7	1317.4
		G/A	100 (15.1)	59 (16.9)	1.15 (0.81–1.64)
		G/G	6 (0.9)	4 (1.1)	1.25 (0.35–4.51)
	Dominant	A/A	555 (84.0)	286 (82.0)	1	0.410	1295.8	1310.5
		G/A-G/G	106 (16.0)	63 (18.1)	1.16 (0.82–1.63)
	Recessive	A/A-G/A	655 (99.1)	345 (98.8)	1	0.760	1296.3	1311.1
		G/G	6 (0.9)	4 (1.1)	1.22 (0.34–4.41)
	Overdominant	A/A-G/G	561 (84.9)	290 (83.1)	1	0.440	1295.9	1310.6
		G/A	100 (15.1)	59 (16.9)	1.15 (0.81–1.64)
	Log-additive				1.14 (0.83–1.57)	0.410	1295.7	1310.5
rs664589	Codominant	C/C	564 (85.3)	292 (83.7)	1	0.220	1295.4	1315.1
		C/G	90 (13.6)	56 (16.1)	1.21 (0.84–1.74)
		G/G	7 (1.1)	1 (0.3)	0.28 (0.03–2.27)
	Dominant	C/C	564 (85.3)	292 (83.7)	1	0.480	1295.9	1310.7
		C/G-G/G	97 (14.7)	57 (16.3)	1.14 (0.80–1.63)
	Recessive	C/C-C/G	654 (98.9)	348 (99.7)	1	0.160	1294.4	1309.2
		G/G	7 (1.1)	1 (0.3)	0.27 (0.03–2.21)
	Overdominant	C/C-G/G	571 (86.4)	293 (84.0)	1	0.290	1295.3	1310.1
		C/G	90 (13.6)	56 (16.1)	1.22 (0.85–1.75)
	Log-additive				1.06 (0.76–1.48)	0.720	1296.3	1311.1
rs3200401	Codominant	C/C	484 (73.2)	224 (64.2)	1	0.008	1288.7	1308.3
		T/C	163 (24.7)	110 (31.5)	1.44 (1.08–1.93)
		T/T	14 (2.1)	15 (4.3)	2.34 (1.10–4.95)
	Dominant	C/C	484 (73.2)	224 (64.2)	1	0.004	1288.2	1302.9
		T/C-T/T	177 (26.8)	125 (35.8)	1.51 (1.14–2.00)
	Recessive	C/C-T/C	647 (97.9)	334 (95.7)	1	0.052	1292.7	1307.4
		T/T	14 (2.1)	15 (4.3)	2.10 (1.00–4.43)
	Overdominant	C/C-T/T	498 (75.3)	239 (68.5)	1	0.026	1291.5	1306.3
		T/C	163 (24.7)	110 (31.5)	1.39 (1.04–1.85)
	Log-additive				1.47 (1.16–1.88)	0.002	1286.7	1301.5
rs1333048	Codominant	C/C	149 (22.5)	93 (26.6)	1	0.190	1295.1	1314.8
		A/C	351 (53.1)	182 (52.1)	0.79 (0.57–1.09)
		A/A	161 (24.4)	74 (21.2)	0.72 (0.49–1.05)
	Dominant	C/C	149 (22.5)	93 (26.6)	1	0.086	1293.5	1308.2
		A/C-A/A	512 (77.5)	256 (73.3)	0.77 (0.57–1.04)
	Recessive	C/C-A/C	500 (75.6)	275 (78.8)	1	0.270	1295.2	1310.0
		A/A	161 (24.4)	74 (21.2)	0.84 (0.61–1.15)
	Overdominant	C/C-A/A	310 (46.9)	167 (47.9)	1	0.590	1296.1	1310.9
		A/C	351 (53.1)	182 (52.1)	0.93 (0.72–1.21)
	Log-additive				0.84 (0.70–1.02)	0.081	1293.4	1308.1
rs4977574	Codominant	A/A	184 (27.8)	75 (21.5)	1	0.032	1291.6	1311.2
		A/G	345 (52.2)	187 (53.6)	1.32 (0.96–1.83)
		G/G	132 (20.0)	87 (24.9)	1.66 (1.13–2.44)
	Dominant	A/A	184 (27.8)	75 (21.5)	1	0.026	1291.5	1306.2
		A/G-G/G	477 (72.2)	274 (78.5)	1.42 (1.04–1.93)
	Recessive	A/A-A/G	529 (80.0)	262 (75.1)	1	0.046	1292.5	1307.2
		G/G	132 (20.0)	87 (24.9)	1.38 (1.01–1.88)
	Overdominant	A/A-G/G	316 (47.8)	162 (46.4)	1	0.780	1296.4	1311.1
		A/G	345 (52.2)	187 (53.6)	1.04 (0.80–1.35)
	Log-additive				1.29 (1.07–1.56)	0.009	1289.6	1304.3
rs1333045	Codominant	C/C	152 (23.0)	94 (26.9)	1	0.160	1294.8	1314.5
		C/T	353 (53.4)	186 (53.3)	0.81 (0.59–1.11)
		T/T	156 (23.6)	69 (19.8)	0.69 (0.47–1.02)
	Dominant	C/C	152 (23.0)	94 (26.9)	1	0.097	1293.7	1308.4
		C/T-T/T	509 (77.0)	255 (73.1)	0.77 (0.57–1.05)
	Recessive	C/C-C/T	505 (76.4)	280 (80.2)	1	0.170	1294.5	1309.3
		T/T	156 (23.6)	69 (19.8)	0.80 (0.58–1.10)
	Overdominant	C/C-T/T	308 (46.6)	163 (46.7)	1	0.780	1296.4	1311.1
		C/T	353 (53.4)	186 (53.3)	0.96 (0.74–1.25)
	Log-additive				0.83 (0.69–1.01)	0.059	1292.9	1307.6
rs10757278	Codominant	A/A	164 (24.8)	77 (22.1)	1	0.150	1294.6	1314.3
		A/G	353 (53.4)	180 (51.6)	1.07 (0.77–1.48)
		G/G	144 (21.8)	92 (26.4)	1.41 (0.97–2.06)
	Dominant	A/A	164 (24.8)	77 (22.1)	1	0.340	1295.5	1310.3
		A/G-G/G	497 (75.2)	272 (77.9)	1.16 (0.85–1.59)
	Recessive	A/A-A/G	517 (78.2)	257 (73.6)	1	0.054	1292.7	1307.5
		G/G	144 (21.8)	92 (26.4)	1.35 (1.00–1.83)
	Overdominant	A/A-G/G	308 (46.6)	169 (48.4)	1	0.420	1295.8	1310.5
		A/G	353 (53.4)	180 (51.6)	0.90 (0.69–1.17)
	Log-additive				1.19 (0.98–1.44)	0.074	1293.2	1308.0

Note: Statistically significant threshold was set at P<0.006 (0.05/8) determined by Bonferroni correction.

Abbreviation: CC, cervical cancer.

Discussion

Long non-coding RNAs (lncRNAs) and their polymorphisms play key roles in the development of different cancers. In the current study, we genotyped SNPs in ANRIL and MALAT1 in cervical cancer patients and healthy control individuals. Our results showed that rs3200401 in MALAT1 and rs4977574 in ANRIL could be associated with the development of CC.

Recently, ANRIL was shown to be upregulated in tumour tissue and function as a tumour-promoting lncRNA in different cancers. In 2018, Zhang et al reported that the expression of lncRNA ANRIL was upregulated in cervical cancer tissues and cell lines.15 In addition, the downregulation of ANRIL could suppress the proliferation and invasion of CC cell, and enhance the apoptosis, indicating that ANRIL plays an oncogenic role in CC development.15 In 2017, Khorshidi et al reported that rs4977574 in ANRIL was not associated with breast cancer in Iran.27 However, Taheri et al found that rs4977574 was associated with benign prostate hyperplasia and prostate cancer risk, and the A allele could be a protective factor for benign prostate hyperplasia (P=0.017, OR=0.66; 95% CI=0.47–0.93) and prostate cancer risk (P=0.001, OR=0.58; 95% CI=0.41–0.81).21 In 2018, Huang et al performed a meta-analysis of the association between ANRIL polymorphisms and cancer risk, and their results found that GG (GG vs AA) and the G allele were associated with an increased risk of cancer (P<0.001, OR=2.40; 95% CI =1.60–3.59, P<0.001, OR = 1.68; 95% CI = 1.35–2.08).28 In the current study, we also found that the allele and genotype frequencies of rs4977574 were significantly different between the CC and control groups (P=0.004 and P=0.014, respectively) and that the A allele might be a protective factor for CC (OR=0.80; 95% CI=0.68–0.93). The rs4977574 SNP could alter the binding sites for ETS transcription factors using the HaploReg v4 software by Ward and Kellis.20 The overexpressed ETS family transcription factor ETS-2 was found to be associated with the development of cervical cell neoplasia.29 These results indicated that rs4977574 could influence the interaction between ANRIL and ETS-2, which is associated with CC.

MALAT1 regulates alternative splicing and transcriptional regulation by modulating the levels of active serine/arginine (SR) proteins such as serine/arginine-rich splicing factor 1 (SRSF1) and SRSF2.30,31 MALAT1 is co-localised with SRSF2 splicing domains and is associated with phosphorylation of SRSF2, interacting with SR proteins as a “molecular sponge,” influencing their stability; MALAT1 regulates the alternative splicing of pre-mRNAs.30,31 In 2019, Qu et al investigated the association between MALAT1 tagSNPs (rs11227209, rs619586, rs664589, and rs3200401) and oesophageal squamous cell carcinoma (ESCC), and their results showed that rs3200401C is associated with an increased risk of ESCC.22 In the current study, we also found that the allele and genotype frequencies of rs3200401 in MALAT1 were significantly different between the control and CC groups, which indicates that the C allele is a protective factor against CC. Our results are similar to those of Qu et al rs3200401 is located in MALAT1 M region (6008-7011 nucleotides), which is one of the binding sites of SRSF2.18 Moreover, using the lncRNA SNP database, Wang et al predicted the rs3200401 potential functions, indicating that C to T leads to free energy change which may alter the structural features of MALAT1. The altered MALAT1 may result in weakened interaction between MALAT1 and SRSF2.19 SRSF2 contributed to the tumour phenotype of HPV16-positive cervical cancer cells.32 The SRSF2 depletion leads to the decreased cell proliferation and colony formation, and increased apoptosis.32 SF2/ASF, SRp20 and SRSF2 are upregulated in a cervical cancer progression model, indicating that they may have oncogenic functions.33 Thus, rs3200401 may alter the expression levels of SRSF2 and participate in CC development.

In the current study, we investigated the association of eight SNPs in ANRIL and MALAT1 between healthy control and patients with CC in a Han Chinese population. Our data showed that rs3200401 in MALAT1 and rs4977574 in ANRIL could play key roles in the development of CC. In the current study, one of the limitations should be our study only found rs3200401 in MALAT1 and rs4977574 in ANRIL were associated with the development of CC. However, the function of these two SNPs in the development of CC should be validated in vitro and in vivo. Another limitation was the relatively modest sample size. In the future, large-scale association studies, especially different populations, are required to clarify the role of these SNPs in the susceptibility and development of CC.