The Molecular Analysis of rs11614913 Polymorphism from miRNA196a Gene and Its Relationship with TNF-α Gene Expression in Cervical Cancer

authors:

avatar Ahmad Hamta ORCID 1 , * , avatar Fatemeh Hajihassani 1

Biology Department, Faculty of Sciences, Arak University, Arak, Iran

how to cite: Hamta A, Hajihassani F. The Molecular Analysis of rs11614913 Polymorphism from miRNA196a Gene and Its Relationship with TNF-α Gene Expression in Cervical Cancer. Jentashapir J Cell Mol Biol. 2020;11(1):e101796. https://doi.org/10.5812/jjcmb.101796.

Abstract

Background:

Cervical cancer (CC) is one of the most common malignant tumors in women, which has been diagnosed as fourth cancer in females worldwide. In addition to human papillomavirus (HPV), genetic factors, including altered expression of some microRNAs and mutations in tumor necrosis factor α (TNF-α) gene, are involved in this cancer.

Objectives:

This study aimed to investigate the rs11614913 polymorphism from the miRNA196a gene and its association with the expression of the TNF-α gene in cervical cancer for early diagnosis and treatment.

Methods:

In this study, 52 samples of pre-cancerous and cancerous lesions, and 50 tissue samples were collected from healthy subjects in an Iranian population. DNA was extracted from the samples, and rs11614913 polymorphism of the miRNA196a gene was investigated by PCR. RNA was extracted from the samples, and the expression of the miRNA196a and TNF-α genes were evaluated. Finally, for data analysis, Epi Info software version 7.1.3.10 and MedCalc Version 19.2.0 were used.

Results:

The frequency of CC, TC, and TT genotypes from rs11614913 polymorphism of miRNA196a gene was 0.58, 0.34, and 0.08, respectively, but in the healthy group it was 0.36, 0.46, and 0.18, respectively. The results also showed that the expression of miRNA196a and TNF-α genes in the patient group was higher than the control group.

Conclusions:

Based on the results of this study, a significant correlation was found between CC genotype and rs11614913 polymorphism of miRNA196a gene and TNF-α gene expression in the cervical cancer sample. Therefore, investigating these factors in patients with cervical cancer may be helpful.

1. Background

One of the most common malignant tumors in women is cervical cancer (CC), which has been diagnosed as fourth cancer in females worldwide (1, 2). Cervical cancer is a component of preventable cancers, and effective screening and appropriate diagnosis and treatment programs have reduced cervical cancer death rates in developed countries (3, 4). Biological, social, and health factors contribute to the formation of this cancer, among which the human papillomavirus (HPV) is known as the most important risk factor for cervical cancer (3-5).

In addition to the HPV virus, genetic factors such as altering the expression of some MicroRNAs (miRNA) also contribute to the development of cervical cancer (6). MicroRNAs can act as oncogenes or tumor suppressors by inhibiting the expression of cancer-related target genes. Besides, functional differences between different types of tumors and various stages of cancer are associated with the expression of microRNAs (7-9).

Moreover, in some studies, the aberrant expression of “mir196a” has been reported in cervical cancer (10). On the other hand, single nucleotide polymorphisms (SNPs) can also be effective in changing the expression of miRNA targets (11). It has also been shown in various studies that miR-196a2 C.T (rs11614913) is involved in various cancers, so we also explored it in this study (12). On the other hand, tumor necrosis factor α (TNF-α) is one of the pro-inflammatory cytokines, which directly contributes to oncogene activation, DNA damage, and, development of cervical lesions (13-15).

2. Objectives

In the present study, the molecular analysis of rs11614913 polymorphism from the miRNA196a gene and its relationship with TNF-α gene expression in cervical cancer was investigated in the Iranian population.

3. Methods

In this study, 52 formalin-fixed paraffin-embedded (FFPE) tissue samples from patients with pre-cancerous and cancerous lesions of cervical tissue and 50 healthy tissue samples were obtained from Khatam Hospital in Tehran. The specimens were examined histopathologically and confirmed by a pathologist.

3.1. DNA Extraction

After deparaffinization of blocks, DNA was extracted by DNA extraction CinnaGen Inc. according to the manufacturer’s protocol and stored at -20ºC.

3.2. Tetra-ARMS PCR Reaction for Proliferation of rs11614913 Polymorphism

Genotype for SNP miR-196a rs11614913 C/T was amplified by T-ARMS PCR in a thermal cycler (Primus 25, Peqlab, Germany) by CinnaGen kits, as well as synthetic primers synthesized by CinnaGen Co (Table 1). The material was poured into PCR-specific tubes, and one sample was considered a negative control that Instead of DNA, 3.5 μL of distilled water was added. The reaction mixture was prepared in a 25 μL volume and after preparation was transferred to the thermocycler apparatus with the program listed in Table 2. Also, the formation of the desired parts in PCR products was investigated by electrophoresis on 2.5% agarose gel.

Table 1.

Sequences of Primers Used for rs11614913 Polymorphism from miRNA196a-2 Gene

PrimerPrimer SequenceReference
Outer-F5’ TCTCTAATCCTTAGGGAGGTTGTGGG 3’(16)
Outer-R5’AAATAAGGGTTCTCCAGACTTGTTCTGC 3’(16)
Inner-F5’AATTTTAAACTCGGCAACAAGAAACGGT 3’(16)
Inner-R5’ GACATAAACCGACTGATGTAACTCCGG 3’(16)
Table 2.

PCR Program for Polymorphism rs11614913 from miRNA196a-2 Gene

SegmentTemperatureIncubation TimeNumber of Cycle
First denaturation94ºC5 minx 1
Denaturation94ºC45 secx 35
Annealing60ºC30 secx 35
Extension72ºC1 minx 35
Final extension72ºC10 minx 1

3.3. RNA Isolation and Quantitative Real-Time PCR (qRT-PCR)

Total RNA was extracted from the tissues with a CinnaPure RNA kit according to the manufacturer protocols. RNA was quantified by NanoDrop spectrophotometer (BioTek Epoch, America).

Before the reverse transcription reaction was performed, to remove any DNA-contamination, the RNA sample was admixed with a DNA-removing mixture for 2 minutes at 42°C. After DNA removal, the QuantaScript RT enzyme was transcribed, and the RNA was transcribed to the cDNA. Synthesis of cDNA from these genes was performed using the QuantiEect® Reverse Transcription Kit (QIAGEN) and specific stem-loop primers for all types of genes according to the protocol.

3.4. Sequence of Specific Stem-Loop Primers for Genes Includes

SNORD-47 (17): 5’-GTCGTATGCAGAGCAGGGTCCGAGG-TATTCGCACTGCATACGACAACCTC-3’

miRNA196a (18):

5’-CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGA-GCCCAACAA-3’

The sequences of primers used are given in Table 3.

Table 3.

Primer Sequences for Real-Time PCR

GenePrimer SequenceReference
miRNA196aF- 5’GAGGCGTGGCAGACTATGC-3’(19)
R- 5’CTTGTACTCCGTCAGCGTGA-3’
SNORD-47F- 5’ATCACTGTAAAACCG TTCCA-3’(20)
R- 5’GAGCAGGGTCCGAGGT-3’
TNF-αF- 5’CCCAGGCAGTCAGATCATCTTC-3’(21)
R- 5’AGCTGCCCCTCAGCTTGA-3’

The expression level of TNF-α and miR-196 in tissues was measured by quantitative Real-time PCR (qRT-PCR) with Snord-47 as the internal reference. Real-time PCR was performed using the ABI 7500 Real-time PCR system (ABI, America) and SYBR GREEN PCR Master Mix (Takara) kit. The amplification conditions for were qRT-PCR were: initial denaturation-95ºc for 5 min, denaturation- 95ºc for 50 sec for 35 cycles, annealing- 60ºc for 40 sec for 35 cycles, extension-75ºc for 1min for 35 cycles and final extension-72ºc for 8 min. The relative expression of miRNA was calculated using Ct values. The difference between the threshold of the desired gene and the internal control gene (housekeeping) can be achieved by obtaining the relative expression of the gene by the 2-ΔΔCt method. Finally, for data analysis, Epi Info software (version 7.1.3.10) and MedCalc (Version 19.2.0) were used.

4. Results

4.1. Tetra-ARMS PCR Results

The frequency of genotypes from rs11614913 polymorphism of the miRNA196a-2 gene was evaluated by using the Tetra-ARMS PCR method. First, the rs11614913 polymorphism of the miRNA196a-2 gene was amplified in all specimens using the PCR method. The desired fragment is 362 bp for miRNA196a-2, which is performed using foreign primers and acts as an internal control in the PCR method, with 362 bp bands present in all samples. According to the statistical analysis of the rs11614913 polymorphism of the miRNA196a-2 gene, the chance of suffering from those with a CC genotype is more than those with TT genotype. We concluded that the CC genotype was a polymorphism of the miRNA196a-2 gene that predisposed cervical cancer (Table 4).

Table 4.

Frequency of Types of CC, TC, and TT Genotypes from rs11614913 Polymorphism in Healthy and Diseased Patients

Genotype/AlleleTotal Number of Patients (N = 52), No. (%)The Total Number of Healthy People (N = 50), No. (%)P ValueaOR 95% CI
TT4 (7.69)9 (18)0.128
CT18 (34.61)23 (46)0.2421.609
CC30 (57.69)18 (36)0.0290.412
CC30 (57.69)18 (36)
TT+CT22 (42.3)33 (64)0.0020.414
CC+CT48 (92.3)41 (82)
TT4 (7.69)9 (18)0.020.347
T26 (25)41 (41)
C78 (75)59 (59)0.0162.085

4.2. Comparison of TNF-α Gene Expression and miRNA196a Expression

In this method, the normalized CT value is measured relative to an untreated sample, and we also need internal standards, their CT values should be deducted from the CT specimen value (normalization).

The relative difference of the test sample versus the control is calculated by 2-∆∆CT formula:

Fold change=2-CT=2-(Ct case-Ct control )=2-Ctcase /2-Ct control 
Fold change miRNA196a=2-Ctcase /2-Ct control =4.4/1.7=2.59
Fold change TNFa=2-Ctcase /2-Ct control =4.2/1.8=2.33

The results of the comparison of miRNA196a and TNF-α gene expression in both healthy and patient groups showed that miRNA196a gene expression was about 2.59 fold, and TNF-α was about 2.33 fold higher in the patients than in the controls (Figure 1).

Comparison of TNF-α and miRNA gene expression (relative to the U47 gene after normalization) in patients and controls
Comparison of TNF-α and miRNA gene expression (relative to the U47 gene after normalization) in patients and controls

4.3. The Relationship Between Types of miRNA196a Gene Polymorphisms by TNF-α Gene Expression

The results showed that there is a significant relationship between the frequency of TT, TC, and CC genotypes and the level of TNF-α expression in the affected individuals (Table 5 and Figure 2). In other words, the CC genotype increases the expression of the TNF-α gene in comparison to the TT genotype of miRNA196a gene polymorphism.

Comparison of miRNA196a gene polymorphism genotypes with the expression of the TNF-α gene
Comparison of miRNA196a gene polymorphism genotypes with the expression of the TNF-α gene
Table 5.

Relationship Between miRNA196a Polymorphism Genotypes and the Expression of the TNF-α Gene

Genotype/AlleleFrequency of Genotypes (%)TNF-α Gene Expression (%)OR 95% CIP Valuea
TT860.2820.287
TC35200.1360.003
CC58749.350.0006
T2511
C75890.3710.01

5. Discussion

This study involves several research aspects. First, the molecular analysis of rs11614913 polymorphism from miRNA196a gene in cervical cancer, and second, the study of the expression of miRNA196a and TNF-α genes in cervical cancer, and finally, the relationship between TNF-α gene expression and genotype polymorphisms rs11614913 from miRNA196a gene in cervical cancer.

Numerous studies have shown the regulatory role of miRNAs in various types of cancer, in which the role of miRNA196a in various types of malignancies has been proven, and rs11614913 as a polymorphism of miRNA196a gene in some cancers is effective (22, 23). The role of TNF-α in various types of cancer has also been shown in numerous studies (13-15). Therefore, a multifaceted study of these factors can help identify biomarkers that are effective in cervical cancer.

Our results in this study showed that the C allele was associated with an increased risk of cervical cancer, and, on the other hand, a significant correlation was found between CC genotype and rs11614913 polymorphism from miRNA196a gene in the cervical cancer sample. Previous studies have also demonstrated the relationship between this SNP and other cancers. For example, one study demonstrated that carriers of the variant homozygote CC of miR-196a-2 rs11614913 were more likely to develop epithelial ovarian cancer (EOC) compared with wild‑type homozygote TT and heterozygote CT carriers (23).

Bodal et al. (24) showed that heterozygous genotype of miR196a2 and combined polymorphism of miR-146 and miR196a2 genes were associated with increased risk of breast cancer in North Indian Women. Guan et al. (7) found that statistically significant associations with HPV16-positive SCCOP and survival for Hsa-mir-196a2 rs11614913. Also, for the first time Tian et al. found that variant genotype CC of miR-196a2 rs11614913 was associated with a significantly increased risk of lung cancer in the Chinese population (25).

However, some studies, unlike our study, did not find a significant relationship between this polymorphism and cancer. For example, Hashemi et al. found no statistically significant association between miR-196a2 rs11614913 variant and PCa risk in the Iranian population (26). Very limited studies have been conducted on the rs11614913 polymorphism of the mir196a gene in cervical cancer. In 2017, the results of Srivastava et al. in India showed that rs3746444 T/C gene polymorphism of the Hsa-miRNA499 gene had a significant relationship with the risk of cervical cancer, while polymorphisms of miRNA146a and miRNA196a2 did not show any association with cervical cancer. Although in the same study of Srivastava, miRNA196a2 gene polymorphism in smokers has a high risk of cervical cancer, no significant changes in the risk of cervical cancer are observed for other polymorphisms (11).

Genetic polymorphism involves a person’s susceptibility to cancer and heavily dependent on environmental mutagens. The reason for this difference in research findings can be in different types of races and geographical conditions, as well as differences in the type of life in each region (27). On the other hand, in the present study, the increase in the expression of the miRNA196a gene in people with cervical cancer has been observed in comparison with healthy subjects, and various studies have also shown similar results to this study. In 2016, a study by Yang et al. In China showed increased a196miRNA gene expression in cancerous ovarian tissue and benign ovarian tissue compared with normal ovarian epithelial tissue (28).

Another study by Liu et al. in China revealed that the serum level of miRNA196a was elevated in patients with cervical cancer compared to healthy subjects. Moreover, the serum level of miRNA196a was associated with CIN grade and various clinical parameters important for cervical cancer (29). Hou et al. also observed a significant increase in miRNA196a in people with cervical cancer compared to healthy subjects (30). Therefore, the expression of the mir196a gene can be considered an effective biomarker for cervical cancer.

In the present study, like mir196a, the expression of TNF-α has also increased in cervical cancer, and other studies have also shown the importance of TNF-α in cancers. For example, a study by Govan et al. in South Africa found that TNF-α plays an important role in all phases of cervical cancer (31). Also, in a study by Al Obeed et al., an increase in TNF-α expression in cancerous tissues was observed in patients with colorectal cancer compared with healthy subjects (15). Furthermore, in another study conducted by Piura et al. in Israel, there was an increase in the expression of TNF-α in cancerous tissues of ovarian cancer compared with healthy tissues (32).

5.1. Conclusions

Considering the results of this study, the expression of TNF-α gene in patients was higher than the control group, and further studies showed that CC genotype increases the expression of TNF-α gene by miRNA196a gene polymorphism than TT genotype. It can be concluded that TNF-α can be a good biomarker for cervical cancer, and more studies with larger sample sizes and different ethnicities can help confirm this acclaim.

Acknowledgements

References

  • 1.

    Berman TA, Schiller JT. Human papillomavirus in cervical cancer and oropharyngeal cancer: one cause, two diseases. Cancer. 2017;123(12):2219-29.

  • 2.

    Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: a cancer journal for clinicians. 2018;68(6):394-424.

  • 3.

    Kuguyo O, Matimba A, Tsikai N, Magwali T, Madziyire M, Gidiri M, et al. Cervical cancer in Zimbabwe: a situation analysis. The Pan African medical journal. 2017;27.

  • 4.

    Momenimovahed Z, Salehiniya H. Cervical cancer in Iran: integrative insights of epidemiological analysis. BioMedicine. 2018;8(3).

  • 5.

    Hemmat N, Baghi HB. Human papillomavirus E5 protein, the undercover culprit of tumorigenesis. Infectious agents and cancer. 2018;13(1):31.

  • 6.

    Pardini B, De Maria D, Francavilla A, Di Gaetano C, Ronco G, Naccarati A. MicroRNAs as markers of progression in cervical cancer: a systematic review. BMC cancer. 2018;18(1):696.

  • 7.

    Guan X, Sturgis EM, Song X, Liu Z, El-Naggar AK, Wei Q, et al. Pre-microRNA variants predict HPV16-positive tumors and survival in patients with squamous cell carcinoma of the oropharynx. Cancer letters. 2013;330(2):233-40.

  • 8.

    He B, Pan Y, Cho WC, Xu Y, Gu L, Nie Z, et al. The association between four genetic variants in microRNAs (rs11614913, rs2910164, rs3746444, rs2292832) and cancer risk: evidence from published studies. PloS one. 2012;7(11). e49032.

  • 9.

    Ma XP, Zhang T, Peng B, Yu L, Jiang DK. Association between microRNA polymorphisms and cancer risk based on the findings of 66 case-control studies. PloS one. 2013;8(11). e79584.

  • 10.

    Liu C, Lin J, Li L, Zhang Y, Chen W, Cao Z, et al. HPV16 early gene E5 specifically reduces miRNA-196a in cervical cancer cells. Scientific reports. 2015;5:7653.

  • 11.

    Srivastava S, Singh S, Fatima N, Mittal B, Srivastava AN. Pre-microRNA Gene Polymorphisms and Risk of Cervical Squamous Cell Carcinoma. Journal of clinical and diagnostic research: JCDR. 2017;11(9):GC01.

  • 12.

    Xu Y, Gu L, Pan Y, Li R, Gao T, Song G, et al. Different effects of three polymorphisms in MicroRNAs on cancer risk in Asian population: evidence from published literatures. PLoS One. 2013;8(6). e65123.

  • 13.

    Babapour N, Mehramiz M, Rastgar Moghadam A, Behboodi N, Yousefi Z, Maftouh M, et al. Association of TNF‐308 G> A polymorphism located in tumor necrosis factor a with the risk of developing cervical cancer and results of pap smear. Journal of cellular biochemistry. 2018.

  • 14.

    Li X, Yin G, Li J, Wu A, Yuan Z, Liang J, et al. The Correlation Between TNF-α Promoter Gene Polymorphism and Genetic Susceptibility to Cervical Cancer. Technology in cancer research & treatment. 2018;17:1533033818782800.

  • 15.

    Al Obeed OA, Alkhayal KA, Al Sheikh A, Zubaidi AM, Vaali-Mohammed M, Boushey R, et al. Increased expression of tumor necrosis factor-α is associated with advanced colorectal cancer stages. World Journal of Gastroenterology: WJG. 2014;20(48):18390.

  • 16.

    Amin-Beidokhti M, Mirfakhraie R, Zare-Karizi S, Karamoddin F. The role of parental microRNA alleles in recurrent pregnancy loss: an association study. Reproductive biomedicine online. 2017;34(3):325-30.

  • 17.

    Jahanafrooz Z, Motamed N, Bakhshandeh B. Effects of miR-21 downregulation and silibinin treatment in breast cancer cell lines. Cytotechnology. 2017;69(4):667-80.

  • 18.

    Xu M, Qiang F, Gao Y, Kang M, Wang M, Tao G, et al. Evaluation of a novel functional single-nucleotide polymorphism (rs35010275 G> C) in MIR196A2 promoter region as a risk factor of gastric cancer in a Chinese population. Medicine. 2014;93(26).

  • 19.

    Qi P, Wang L, Zhou B, Yao WJ, Xu S, Zhou Y, et al. Associations of miRNA polymorphisms and expression levels with breast cancer risk in the Chinese population. Genet Mol Res. 2015;14(2):6289-96.

  • 20.

    Kazemzadeh S, Farsinejad AR, Sabour Takanlu J, Kaviani S, Atashi A, Soleimani M, et al. miR-146a: A possible tumor suppressor in multiple myeloma. Scientific Journal of Iran Blood Transfus Organ. 2016;13(3):224-32.

  • 21.

    Rasmi Y, Bagheri M, Faramarz-Gaznagh S, Nemati M, Khadem-Ansari MH, Saboory E, et al. Transcriptional activity of tumor necrosis factor-alpha gene in peripheral blood mononuclear cells in patients with coronary slow flow. ARYA atherosclerosis. 2017;13(4):196.

  • 22.

    Hong Y, Kang H, Kwak J, Park BL, You C, Kim Y, et al. Association between microRNA196a2 rs11614913 genotypes and the risk of non-small cell lung cancer in Korean population. Journal of Preventive Medicine and Public Health. 2011;44(3):125.

  • 23.

    Song Z, Wu Y, Zhao H, Liu C, Cai H, Guo B, et al. Association between the rs11614913 variant of miRNA-196a-2 and the risk of epithelial ovarian cancer. Oncology letters. 2016;11(1):194-200.

  • 24.

    Bodal VK, Sangwan S, Bal MS, Kaur M, Sharma S, Kaur B. Association between Microrna 146a and Microrna 196a2 Genes Polymorphism and Breast Cancer Risk in North Indian Women. Asian Pacific journal of cancer prevention: APJCP. 2017;18(9):2345.

  • 25.

    Tian T, Shu Y, Chen J, Hu Z, Xu L, Jin G, et al. A functional genetic variant in microRNA-196a2 is associated with increased susceptibility of lung cancer in Chinese. Cancer Epidemiology and Prevention Biomarkers. 2009;18(4):1183-7.

  • 26.

    Hashemi M, Moradi N, Ziaee SAM, Narouie B, Soltani MH, Rezaei M, et al. Association between single nucleotide polymorphism in miR-499, miR-196a2, miR-146a and miR-149 and prostate cancer risk in a sample of Iranian population. Journal of advanced research. 2016;7(3):491-8.

  • 27.

    Wünsch Filho V, Zago MA. Modern cancer epidemiological research: genetic polymorphisms and environment. Revista de saude publica. 2005;39:490-7.

  • 28.

    Yang B, Li S, Ma L, Liu H, Liu J, Shao J. Expression and mechanism of action of miR-196a in epithelial ovarian cancer. Asian Pacific journal of tropical medicine. 2016;9(11):1105-10.

  • 29.

    Liu P, Xin F, Ma CF. Clinical significance of serum miR-196a in cervical intraepithelial neoplasia and cervical cancer. Genet Mol Res. 2015;14(4):17995-8002.

  • 30.

    Hou T, Ou J, Zhao X, Huang X, Huang Y, Zhang Y. MicroRNA-196a promotes cervical cancer proliferation through the regulation of FOXO1 and p27 Kip1. British journal of cancer. 2014;110(5):1260.

  • 31.

    Govan VA, Carrara HR, Sachs JA, Hoffman M, Stanczuk GA, Williamson A. Ethnic differences in allelic distribution of IFN-g in South African women but no link with cervical cancer. Journal of Carcinogenesis. 2003;2:3.

  • 32.

    Piura B, Medina L, Rabinovich A, Dyomin V, Levy RS, Huleihel M. Distinct expression and localization of TNF system in ovarian carcinoma tissues: possible involvement of TNF-α in morphological changes of ovarian cancerous cells. Anticancer research. 2014;34(2):745-52.