Rs6757 in microRNA-3976 binding site of CD147 confers risk of hepatocellular carcinoma in South Chinese population
World Journal of Surgical Oncology volume 20, Article number: 260 (2022)
Cluster of differentiation 147 (CD147) overexpression plays a key role in the proliferation, differentiation, invasion, metastasis, and prognosis of hepatocellular carcinoma (HCC). The aim of this study was to explore the relationship between rs6757 and the HCC risk in the South Chinese population, and the functional significance of rs6757 by affecting the efficacy of microRNA-3976 (miR-3976) binding to the CD147 3′-UTR.
We performed a retrospective case-control study to analyze the association between rs6757 and the risk of HCC. We chose candidate microRNAs with the potential of interacting with rs6757 through a series of silico analyses. A luciferase reporter gene assay was implemented to detect the binding extent of microRNAs to each polymorphic allele of rs6757.
An obvious association between rs6757 and the risk of HCC was detected in C vs. T (OR = 1.826, 95% CI [1.263–2.642]), CC vs. TT (OR = 4.513, 95% CI [1.510–13.489]), dominant genetic model (OR = 1.824, 95% CI [1.120–2.965]), and recessive genetic model (OR = 3.765, 95% CI [1.286–11.020]). Bioinformatics analysis indicated that miR-3976 binding sites containing the rs6757-T allele had lower free energies than those with the C allele, the lower free energies, the higher affinities. Luciferase activity was remarkably decreased by miR-3976 binding to the CD147 3′-UTR bearing rs6757 T allele, which could be reversed by miR-3976 inhibitors. Furthermore, miR-3976 reduced the luciferase expression in a manner of dose-dependent when cotransfected with constructs with the CD147-TT-pSICHECK2.
The research we have done suggests that rs6757 confers the CD147 allele-specific translational suppression by miR-3976, which provides a theoretical basis for antineoplastic therapy targeting CD147.
HCC is one of the most common types of all cancers worldwide, and is the fifth most common cancer among men and the ninth most common cancer among women [1, 2]. Globally, HCC is the third leading cause of cancer-related deaths . Although it is commonly recognized that HCC is caused by a combination of genetic, immunologic, and environmental factors, its complete etiology has not been fully understood.
CD147, as well as basigin (BSG), is a transmembrane glycoprotein which regulates the level of matrix metalloproteinase (MMPs) through cell–matrix and cell–cell interaction. Previous work on CD147 have indicated that CD147 is involved in a variety of pathophysiological functions . CD147 is highly expressed in a variety of tumors and metastatic cancer cells [5, 6]. The expression level of CD147 is not only related to the age of the patient and tumor type, but also related to the clinical stage and histopathologic type of the tumor . Disruption of Arf6-mediated CD147 trafficking markedly decreased the migration and invasion of liver cancer cells, while high expression of the Arf6-CD147 signaling components in HCC was closely correlated with poor clinical outcome of patients . Hypophosphorylated CD147 promotes the migration and invasion of HCC cells and correlates with an unfavorable prognosis in HCC patients . TGF-β1 stimulation upregulated CD147 expression and mediated the differentiation of HCC cells. The overexpression of CD147 induced the dedifferentiation and enhanced the malignancy of HCC cells, and also increased the transcriptional expression of TGF-β1 by activating β-catenin . Considering the important role of CD147 in tumors, we speculated that CD147 genetic polymorphisms might influence the occurrence and subsequent development of HCC.
MicroRNAs (miRNAs) are endogenous non-coding RNAs about 22 nt which may regulate gene expression by targeting mRNAs for degradation or translational repression at the post-transcriptional level . The previous work on miRNAs have indicated that miRNAs are associated with many biological processes such as cell proliferation, apoptosis, differentiation, and tumorigenesis . Recently, more and more studies focus on single nucleotide polymorphisms (SNPs) in miRNAs target sites (miRSNPs) which can enhance or weaken miRNA–mRNA interaction s[13,14,15]. MiRSNPs of susceptibility genes could affect complex traits/diseases through changing the function of miRNAs [16, 17]. Any changes in the recognition site could destroy binding sites or modify binding affinity, which could result in the evasion of miRNA negative regulation . However, until now there have been no studies reported about the association between the miRSNPs of CD147 and the risk of HCC.
Therefore, the purpose of this study was to investigate the association between the rs6757 and the risk of HCC in South Chinese population through a new mechanism of precise gene regulation by miRNA binding ability.
Materials and methods
Patients and controls
We performed a case-control study including 155 HCC patients (135 males and 20 females with an average age of 50.7) and 151 healthy controls (98 males and 53 females with an average age of 45.1) on a random basis admitted to the hospital for treatment at the First Affiliated Hospital of Guangxi Medical University from January 2017 to June 2017. All HCC patients were positive for hepatitis B surface antigen (HBsAg) and had never received radiotherapy or chemotherapy. The diagnostic criteria for primary liver cancer are the guidelines for the diagnosis and treatment of primary liver cancer formulated by the Ministry of Health of the People's Republic of China in the 2017 edition. A total of 151 healthy controls were recruited from the Physical Examination Center of the First Affiliated Hospital of Guangxi Medical University during the same period. Inclusion criteria were HBsAg positive; liver function, serum tumor markers, blood lipids, blood glucose, and other test results were normal; and X-ray or B-ultrasound inspection showed no tumors.
DNA isolation and genotyping
Whole blood samples were collected by vacuum blood collection with EDTA-K2 anticoagulant. Genomic DNA was isolated according to the instructions of a Blood Genomic DNA Extraction Kit (Axygen Scientific Inc., USA), and stored at – 80 °C. Primers for the rs6757 were designed using Primer Premier software (version 6.0, Palo Alto, CA). The upstream and downstream primers used were 5′-GCCTGTCTGAAGCCAATGC-3′ and 5′-AACATGCCGAGAGGAGTAACA-3′. The PCR conditions for amplifying targets were predenature 1 cycle of 5 min at 94 °C, then 35 cycles of 30 s at 94 °C, 30 s at 65 °C, and 30 s at 72 °C, followed by 5 min at 72°C. PCR products were directly sequenced to determine the CD147 rs6757 by Genomic Company (Shenggong Biological Engineering Co, Shanghai, China).
In silico analyses
MirSNP (http://cmbi.bjmu.edu.cn/mirsnp), a database of polymorphisms predicting miRNA–mRNA binding sites, identifies miRNA-related SNPs in GWAS SNPs and eQTLs . Apart from this, we analyzed the miRNA binding sites in CD147 3′-UTR using miRBase (http://www.mirbase.org/index.shtml) and TargetScan (http://www.targetscan.org/). With these bioinformatics tools, we identified hsa-miR-3976 potentially binding to a stretch of sequence harboring the rs6757. Subsequently, we utilized RNAhybrid (https://bibiserv.cebitec.uni-bielefeld.de/rnahybird) to predict the miRNA-target binding structures and energies.
Plasmid construction and luciferase assay
The TT genotype of rs6757 was amplified by PCR using CD147-TT-F: 5′-ccgctcgagGGCAGGTGGCCCGAGGACGC-3′ and CD147-TT-R: 5′-ataagaatgcggccgcGAGGGTGGAGGTGGGGGCGA-3′, which were tagged with the XhoI and NotI restriction enzyme sites (underlined) respectively. The rs6757 TT genotype inserts were digested by XhoI/NotI, gel-purified, extracted, and then subcloned into the XhoI/NotI site of pSICHECK2 vector (Promega). The prepared plasmid was designated as CD147-TT-pSICHECK2. The inserts were sequenced to verify the validity and orientation relative to the luciferase gene at Genomic Company (Shenggong Biological Engineering Co, Shanghai, China). Using the site-directed mutagenesis, we constructed CD147-CC-pSICHEC-K2 that contains rs6757 CC genotype sequence by using CD147-CC-F: 5′-GAGGACGGCCGGCTCTCCATAGCACCAGGGCTCACGTGG-3′, and CD147-CC-R: 5′-CCACGTGAGCCCTGGTGCTATGGAGAGCCGGCCGTCCTC-3′ primers. The inserts of CD147-CC-pSICHECK2 were also confirmed by sequencing.
HEK293T cells (from Shenggong Biological Engineering Co.) were cultured with high-sugar DMEM medium containing 10% fetal bovine serum at 37 °C in humidified 5% CO2.
Transfection and luciferase assay
2 × 104 cells per well were seeded in 96-multiwell plates. After the cells reached 50–60% confluence, we used Lipofectamine 2000 reagent (Invitrogen) to transfect the prepared CD147-TT-pSICHECK2 or CD147-CC-pSICHECK2 (500 ng), 50 nmol/L hsa-miR-3976 mimics, hsa-miR-3976 control (as a control for hsa-miR-3976, the specific sequence is: UCACAACCUCCUAGAAAGAGUAGA), 100 nmol/L hsa-miR-3976 inhibitor (reverse sequence of hsa-miR-3976), and 100 nmol/L hsa-miR-3976 negative control (as a control for hsa-miR-3976 inhibitor, the specific sequence is UUUGUACUACACACAAAAGUACUG) into the cells. Over and above these, various concentrations (1.0, 2.5, 5.0, and 10 nM) of hsa-miR-3976 mimics were added to the transfection system to observe the effect of hsa-miR-3976 on luciferase activity at different concentrations. Each transfection was carried out in triplicate. The data were analyzed by normalizing firefly luciferase activity with that of the Renilla luciferase.
Continuous variables were compared using an unpaired t-test. The differences of alleles and genotypes frequency distribution between the groups were statistically analyzed using the χ2 test. The Hardy–Weinberg equilibrium (HWE) was tested by a Q test with one degree of freedom [14, 15]. The association between CD147 rs6757 and the risk of HCC was estimated using the odds ratio (OR) and the 95% confidence interval (CI). After adjustment for age, gender, smoking, and drinking, adjusted ORs and 95% CIs were calculated by logistic regression. The value of relative luciferase activity were presented as means ± SD and analyzed by Student’s t test. All statistics were performed by SPSS software version 19.0 and Stata software version 11.0. It was considered to be statistically significant when P < 0.05.
A total of 155 HCC patients and 151 controls were enrolled in this study. There were significant differences in gender, age, and drinking status between HCC patients and controls. However, no significant difference was found regarding smoking (Table 1).
Association of the CD147 rs6757:T>C with HCC susceptibility
The T/C variant (rs6757) was genotyped by sequencing in all 306 participants (Table 2). As shown in Table 2, the frequencies of the rs6757 genotypes and alleles between the HCC patients and controls were statistically significant (χ2 = 10.721, p = 0.005; χ2 = 10.364, p = 0.001). The distributions of genotypes in controls and HCC patients were compatible with HWE (PControls = 0.558, PHCC = 0.351). The detailed results are shown in Table 3. After adjusting for confounding factors including age, gender, smoking, and drinking, significant associations between the T/C variant (rs6757) and the risk of HCC were detected in C vs. T (OR = 1.826, 95% CI [1.263–2.642]), CC vs. TT (OR = 4.513, 95% CI [1.510–13.489]), dominant genetic model (OR = 1.824, 95% CI [1.122–2.965]), and recessive genetic model (OR = 3.765, 95% CI [1.286–11.020]). C variant (rs6757) was considered to be an independent risk factor for HCC. The detailed results are shown in Table 4.
Rs6757 is located in the predicted binding sites of hsa-miR-3976 in the CD147 3′-UTR
With the development of bioinformatics, there are more and more websites or software tools for predicting miRNA target genes online. Through integrating and analyzing the prediction results of different websites or software such as MirSNP (http://cmbi.bjmu.edu.cn/mirsnp), miRBase (http://www.mirbase.org/index.shtml), and TargetScan, we found that there were multiple potential miRNAs in the CD147 3′-UTR region. After intersecting these predicted results and comparing them in the NCBI database, we found that the rs6757 was exactly located in the seed region of has-miR-3976 (Fig. 1).
The forecasted minimal folding energies of the hsa-miR-3976 target duplexes for the rs6757-T and -C alleles were − 24.7 kcal/mol and − 20.2 kcal/mol respectively (Fig. 2). Obviously, the minimal folding energies of the duplexes rs6757-T-containing (− 24.7 kcal/mol) is lower than that of the duplexes rs6757-C-containing (− 20.2 kcal/mol). The lower minimal folding energies, the more stable binding of hsa-miR-3976 to the CD147 3′-UTR, which indicated more efficient translational repression for the CD147 transcript containing rs6757 T variant. The predicted sum of the ΔΔG value was 4.5 kcal/mol, which suggested that the rs6757 was more likely to be a functional SNP.
Functional analyses of rs6757
In order to detect whether the T/C variant (rs6757) affected the CD147 translational regulation by affecting the binding of miR-3976 to CD147 3′-UTR, we constructed the luciferase reporter vectors containing rs6757 different genotypes (Figs. 3, 4, 5, and 6). Luciferase expression was significantly reduced following transfection with CD147-TT-pSICHECK2 (P < 0.05); this translational suppression of miR-3976 can be reversed by its inhibitor (Fig. 7A). However, there was no significant statistical difference between the groups following transfection with CD147-CC-pSICHECK2 (Fig. 7B). Not only so, miR-3976 dose-dependently reduced the luciferase expression when cotransfected with constructs with the CD147-TT-pSICHECK2 (Fig. 7C). However, miR-3976 showed no effect on the luciferase expression when cotransfected with constructs with the CD147-CC-pSICHECK2 (Fig. 7C).
In accordance with the principle of genetic diversity, differences in evolution and environment result in differences in genes and phenotypes between different populations [20, 21]. Therefore, it is important to study the relationship between these genetic variants and genetic susceptibility of different people. To our knowledge, no study has been conducted on rs6757:T>C and HCC risk in the South Chinese population.
In this research, we calculated the associations between the rs6757 and risk of HCC in 155 HCC patients and 151 controls. We find an obvious relationship between C/T variant (rs6757) and the genetic susceptibility of HCC, participants carrying allele C had a distinctly increased risk of HCC (OR = 1.826, 95% CI [1.263–2.642]; P = 0.001).
As shown in Fig. 1, the T/C variant (rs6757) is located in the seed region, nucleotides 2−7, of miRNA-3976, which means that the C/T variant (rs6757) has the potential to disrupt the binding of miR-3976 to CD147 3′-UTR. We predicted and calculated the minimal free energies of the secondary structures of miR-3976 bounding to the CD147 3′-UTR containing T/C variant (rs6757) by computational modeling. The minimal folding energy and the ΔΔG value were used to assess the impact of SNPs on the binding ability of the miRNA with its target mRNA . The minimal folding energies of the duplexes rs6757-T-containing (− 24.7 kcal/mol) is lower than that of the duplexes rs6757-C-containing (− 20.2 kcal/mol). The lower minimal folding energies, the more stable binding of miR-3976 to the CD147 mRNA. The sum of the ΔΔG value was 4.5 kcal/mol, the higher sum of the ΔΔG value of the rs6757, the more likely to be a functional SNP.
As expected, luciferase reporter gene assay showed that miRNA-3976 suppressed the luciferase construct bearing the TT, but not the CC, genotype of the rs6757 (Fig. 7). Our investigation confirm that the C/T variant (rs6757) has an obvious effect on miR-3976 binding affinity, which causes greater suppression when miR-3976 bounding to CD147 3′-UTR containing the T allele rather than the C allele.
Despite many new discoveries mentioned above, there are some limitations in this study. First, although we selected health individuals without signs or symptoms of HCC as the control group, it should be noticed that the controls did not undergo a CT or MRI plain scan. Second, the sample size in the present study was still not large enough to reveal genes with allele frequencies between 5 and 10% . Third, taking into account geographic variations, the prevalence of C/T variant (rs6757) in the South Chinese population may bias the results based on the single-center case-control studies. Fourth, which confounding factors should be selected for multivariate analysis has always been a controversial problem. Generally, most researchers choose the following two methods to solve this problem. One is that only dissimilar confounding factors should be corrected. The other suggests all confounding factors should be corrected. In this paper, we adopted the second approach. Actually, no matter which method you choose, the statistical results have no significant change. Finally, We simply verified that the rs6757:T>C would affect on the binding of miR-3976 to CD147 3′-UTR using the luciferase reporter assay in HEK293T cells. Next, we will study the relationship between the rs6757 and HCC in vitro at least two liver cell lines.
To sum up, we have revealed that CD147 is a direct target of miR-3976, and the C/T variant (rs6757) could affect on the binding of miR-3976 to CD147 3′-UTR, which result in a statistically significant association between the rs6757 and risk of HCC in the South Chinese population. However, further studies are needed to explore the new mechanism of subtle gene regulation of CD147 through miR-3976 binding capacity affected by rs6757.
Availability of data and materials
El-Serag HB. Epidemiology of viral hepatitis and hepatocellular carcinoma. Gastroenterology. 2012;142:1264–73.
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136:E359–86.
Grandhi MS, Kim AK, Ronnekleiv-Kelly SM, Kamel IR, Ghasebeh MA, Pawlik TM. Hepatocellular carcinoma: from diagnosis to treatment. Surg Oncol. 2016;25:74–85.
Xiong L, Edwards CK 3rd, Zhou L. The biological function and clinical utilization of CD147 in human diseases: a review of the current scientific literature. Int J Mol Sci. 2014;15:17411–41.
Papoutsi M, Kurz H, Schächtele C, Marmé D, Christ B, Pröls F, et al. Induction of the blood–brain barrier marker neurothelin/HT7 in endothelial cells by a variety of tumors in chick embryos. Histochem Cell Biol. 2000;113:105–13.
Ran XY, Huang J, Zhang HZ, Jiang ZM, Chen J. Factors related to biologic behavior in giant cell tumor of bone. Zhonghua Bing Li Xue Za Zhi. 2013;42:669–74.
Walters DK, Arendt BK, Jelinek DF. CD147 regulates the expression of MCT1 and lactate export in multiple myeloma cells. Cell Cycle. 2013;12:3175–83.
Qi S, Su L, Li J, Zhang C, Ma Z, Liu G, et al. Arf6-driven endocytic recycling of CD147 determines HCC malignant phenotypes. J Exp Clin Cancer Res. 2019;38:471.
Jin J, Wang SJ, Cui J, Li L, Li JY, Liu FL, et al. Hypo-phosphorylated CD147 promotes migration and invasion of hepatocellular carcinoma cells and predicts a poor prognosis. Cell Oncol (Dordr). 2019;42:537–54.
Wu J, Lu M, Li Y, Shang YK, Wang SJ, Meng Y, et al. Regulation of a TGF-β1-CD147 self-sustaining network in the differentiation plasticity of hepatocellular carcinoma cells. Oncogene. 2016;35:5468–79.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Tian XP, Wang CY, Jin XH, Li M, Wang FW, Huang WJ, et al. Acidic microenvironment up-regulates exosomal miR-21 and miR-10b in early-stage hepatocellular carcinoma to promote cancer cell proliferation and metastasis. Theranostics. 2019;9:1965–79.
Teng MS, Hsu LA, Juan SH, Lin WC, Lee MC, Su CW, et al. A GDF15 3' UTR variant, rs1054564, results in allele-specific translational repression of GDF15 by hsa-miR-1233-3p. PLoS One. 2017;12:e0183187.
Shasttiri A, Rostamian Delavar M, Baghi M, Dehghani Ashkezari M, Ghaedi K. SNP rs10800708 within the KIF14 miRNA binding site is linked with breast cancer. Br J Biomed Sci. 2019;76:46–8.
Pirooz HJ, Jafari N, Rastegari M, Fathi-Roudsari M, Tasharrofi N, Shokri G, et al. Functional SNP in microRNA-491-5p binding site of MMP9 3'-UTR affects cancer susceptibility. J Cell Biochem. 2018;119:5126–34.
Sun G, Yan J, Noltner K, Feng J, Li H, Sarkis DA, et al. SNPs in human miRNA genes affect biogenesis and function. RNA. 2009;15:1640–51.
Lin J, Horikawa Y, Tamboli P, Clague J, Wood CG, Wu X. Genetic variations in microRNA-related genes are associated with survival and recurrence in patients with renal cell carcinoma. Carcinogenesis. 2010;31:1805–12.
Yuan Y, Weidhaas JB. Functional microRNA binding site variants. Mol Oncol. 2019;13:4–8.
Liu C, Zhang F, Li T, Lu M, Wang L, Yue W, et al. MirSNP, a database of polymorphisms altering miRNA target sites, identifies miRNA-related SNPs in GWAS SNPs and eQTLs. BMC Genomics. 2012;13:661.
Leimar O. The evolution of phenotypic polymorphism: randomized strategies versus evolutionary branching. Am Nat. 2005;165:669–81.
Messer PW, Ellner SP, Hairston NG Jr. Can Population Genetics Adapt to Rapid Evolution? Trends Genet. 2016;32:408–18.
Moraghebi M, Maleki R, Ahmadi M, Negahi AA, Abbasi H, Mousavi P. In silico analysis of polymorphisms in microRNAs deregulated in Alzheimer disease. Front Neurosci. 2021;15:631852.
Stranger BE, Nica AC, Forrest MS, Dimas A, Bird CP, Beazley C, et al. Population genomics of human gene expression. Nat Genet. 2007;39:1217–24.
My deepest gratitude goes first and foremost to Professor Jinbao Zong, my supervisor, whose instructive suggestions and expert advice in the field of hepatocellular carcinoma have helped me a lot in writing this thesis. I hereby extend my grateful thanks to the other teachers as they generously helped me collect materials and gave me many invaluable suggestions.
No funding approved.
Ethics approval and consent to participate
This study protocol was performed in accordance with the ethical standards as laid down in the 1964 Declaration of Helsinki and was approved by the Ethics Committee of Guangxi Medical University. Most participants sign informed consent, due to some participants' illiteracy, physical disability or other limiting condition, verbal consent is also an alternative.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Guo, F., Li, H., Wang, L. et al. Rs6757 in microRNA-3976 binding site of CD147 confers risk of hepatocellular carcinoma in South Chinese population. World J Surg Onc 20, 260 (2022). https://doi.org/10.1186/s12957-022-02724-w