Abstract
microRNAs’ (miRNAs) loci may influence hepatocellular carcinoma (HCC) development. Many recent studies have assessed the relationship between miRNA-499, miRNA-146a, and miRNA-196a2 loci and HCC risk. However, the observed results are conflicting. A total of 584 HCC patients and 923 age- and sex-matched controls were recruited. The correlation of miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 rs11614913 with HCC development was assessed. In the <53-year-old subgroup, a correlation of the rs2910164 locus with HCC risk was found (GG/CG vs. CC: adjusted p = 0.011, GG vs. CC: adjusted p = 0.021 and CG vs. CC: adjusted p = 0.027). The association between miRNA-146a rs2910164 and the risk of HCC was also found in the never smoking (GG/CG vs. CC: adjusted p = 0.011 and CG vs. CC: adjusted p = 0.018). Using false-positive report probability method and power value, we identified that miRNA-146a rs2910164 conferred a risk to HCC in the <53-year-old and never-smoking subgroups. In conclusion, this study indicates rs2910164 may be a risk factor for HCC, especially in the <53-year-old and never-smoking subgroups.
Introduction
In China, an estimated 466,100 new liver carcinoma (LC) cases and 422,100 LC-related deaths occurred in 2015, which might rank as the fourth most frequent cancer and the third leading cause of their related death (Chen et al., 2016). Hepatocellular carcinoma (HCC) is regarded as one of the major subtypes of LC. The etiology of HCC is thought to be associated with a prevalence of hepatitis B virus (HBV) in China. However, HCC is a complex disease and a number of hereditary and environmental factors may contribute to the development of HCC.
microRNAs (miRNAs) contain about 22 nucleotides. They are common noncoding RNA molecules. In eukaryotes, miRNAs are of vital importance for regulation of gene expression (Bartel, 2004). Many investigations have suggested that variations in normal functions of cells are associated with miRNAs. Deficiency of miRNAs can result in a number of diseases. It is believed that miRNAs are important for many vital roles (e.g., apoptosis, development, proliferation, differentiation, and carcinogenesis) (Drakaki and Iliopoulos, 2009; Palma et al., 2014; Wang et al., 2016b; Zhao et al., 2017; Kim et al., 2018; Ou et al., 2018). In addition, accumulating evidence highlights that miRNAs may contribute to infection, immune, and inflammatory response (Jovanovic and Hengartner, 2006). As well, the association of miRNAs with the development of cancer was also studied (He et al., 2018; Zhou et al., 2019). It is reported that HCC is associated with infection of HBV or hepatitis C virus (Ghouri et al., 2017; Stewart et al., 2019). Thus, miRNAs may be implicated in the etiology of HCC (Zhang et al., 2012; Fiorino et al., 2016; Hayes and Chayama, 2016; Sadri Nahand et al., 2019).
There are current studies focusing on the correlation of the rs3746444 locus in miRNA-499 with HCC development (Akkiz et al., 2011a; Xiang et al., 2012). Some reports have indicated that the rs3746444 A>G single nucleotide polymorphism (SNP) in miRNA-499 increased the risk of HCC (Xiang et al., 2012; Chu et al., 2014; Ma et al., 2014; Qi et al., 2014). Several publications showed that rs3746444 in miRNA-499 might decrease the risk of HCC (Zou and Zhao, 2013). However, more case–control studies suggested that rs3746444 variants were not associated with HCC risk (Akkiz et al., 2011a; Shan et al., 2013; Yan et al., 2015; Farokhizadeh et al., 2019). Recently, some meta-analyses explored the relationship between this locus and HCC risk and the observations were also controversial. Some meta-analyses suggested that the rs3746444 A>G polymorphism in miRNA-499 contributed to HCC development (Liu et al., 2014b; Ye et al., 2016; Yu et al., 2016, 2017; Qiu et al., 2018; Yang et al., 2018). However, other meta-analyses failed to find an association of the rs3746444 SNP in miRNA-499 with HCC risk (Hu et al., 2013; Wang et al., 2013, 2016a; Xu et al., 2013a, b; Liu et al., 2014a; Zhu et al., 2016b; Zheng et al., 2017). Thus, the association between this SNP and the development of HCC was conflicting.
Khanizadeh et al. (2019) and Al-Qahtani et al. (2017) reported that rs2910164 C>G SNP in miRNA-146a conferred an increased susceptibility to chronic infection by HBV. Thus, this polymorphism may facilitate the development of HCC. Another study also showed that this SNP might influence the development of HCC (Xu et al., 2008). However, other studies suggested that rs2910164 SNP might not play a vital role in HCC development (Akkiz et al., 2011b; Zhang et al., 2011, 2013; Shan et al., 2013; Chu et al., 2014; Qi et al., 2014; Li et al., 2015; Riazalhosseini et al., 2016; Farokhizadeh et al., 2019). Several meta-analyses also observed conflicting findings (He et al., 2012; Wang et al., 2012; Hu et al., 2013; Xu et al., 2013a; Yin et al., 2013; Zhang et al., 2017; Hao et al., 2018). The inconsistent results may be due to multiple gene, social, and environmental factors.
MiRNA-196a2 rs11614913 T>C, another common miRNA SNP, was also explored for a relationship to the risk of HCC. Rs11614913 is located in the 3′-UTR of the miRNA-196a2 precursor. This SNP alters the transcription of mature miRNA-196a (Hoffman et al., 2009) and rs11614913 might be involved in HBV infection and its related liver complications (Al-Qahtani et al., 2017). In a meta-analysis, the miRNA-196a2 rs11614913 C allele was found to be associated with an increased risk of HCC (Liu et al., 2018). However, most of the previous studies included had small sample sizes and the observation may not be accurate.
In this study, we included 1507 participants and assessed the relationship of miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 rs11614913 with HCC risk.
Materials and Methods
Subjects
In total, 525 male and 59 female patients with HCC were recruited from the Fujian Medical University Union Hospital (Fuzhou City, China) and No. 900 Hospital of Chinese People's Liberation Army (Fuzhou City, China) between 2002 and 2016. In this study, each patient had surgery and was pathologically confirmed with HCC. The control group constituted 923 subjects (835 male and 88 female) without any cancer history. Each participant signed a written informed consent. By using a questionnaire, we collected the following information: drinking, age, sex, and tobacco consumption. By using biochemical examination, we evaluated the status of chronic HBV infection in each subject. This investigation was approved by the Ethics Committee of Changzhou No. 3 People's Hospital (Approval No. 2018-04).
Polymorphisms selection
To determine the role of miRNA polymorphisms in the development of HCC, we selected miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 rs11614913 SNPs according to some published literatures. In these studies, rs3746444, rs2910164, and rs11614913 were found to be correlated with cancer (Chen et al., 2018; Ahmad et al., 2019; Tang et al., 2019; Yan et al., 2019).
Genotying
We collected 1507 blood samples with an EDTA anticoagulant and they were stored at −80°C. Extracting DNA was carefully carried out from leucocytes. The quality and concentration was measured by NanoDrop ND-1000 microspectrophotometer (Thermo Fisher Scientific, Inc., Waltham, MA). The concentration of genomic DNA was 20–50 ng/μL. The 260/280 nm ratio was 1.8–2.0. Using an ABI 3730XL sequencer, the SNPscan™ method was used to genotype the miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 rs11614913 loci (Chen et al., 2018). A second technician selected 4% DNA samples randomly, and the results indicated that genotypes of the first assay were repeatable.
Statistical analyses
The mean ± standard deviation was used to express the distribution of age. The difference of age was measured with Student's t-test. The categorical variable (e.g., the number of genotypes, alcohol consumption, age, gender, and smoking) between the two groups were compared with a χ
2-test. Odds ratios (ORs) and their confidence intervals (95% CIs) were used to determine the relationship between miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 rs11614913 and HCC risk. An online Hardy–Weinberg equilibrium (HWE) calculation program (
Results
Case–control study
A total of 1507 sex/age-matched participants (584 HCC cases/923 controls) were recruited. No significant difference was found in terms of age (HCC cases vs. controls [mean age]: 53.17 ± 11.76 years vs. 53.72 ± 9.97 years, p = 0.327), and gender (number of cases vs. controls [male/female]: 525/59 vs. 835/88, p = 0.717). The ratio of tobacco consumption was similar between two groups (number of cases vs. controls [never/ever]: 374/210 vs. 596/327, p = 0.834). Differences in the distribution of chronic HBV infection and alcohol drinking were found (chronic HBV infection of HCC cases vs. controls [yes/no]: 412/172 vs. 85/838, p < 0.001; drinking of cases vs. controls [never/ever]: 414/170 vs. 775/148, p < 0.001). In controls, the genotype of the miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 rs11614913 conformed to the HWE. The minor allele frequency (MAF) of miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 was 0.15, 0.38, and 0.43, respectively (Table 1).
Primary Information for miR-146a rs2910164 C>G, miR-196a2 rs11614913 T > C and miR-499 rs3746444 A>G Polymorphisms
HWE, Hardy–Weinberg equilibrium; MAF, minor allele frequency; SNP, single nucleotide polymorphism.
Relationship of miRNA SNPs with HCC
The genotype of the miRNA SNPs are summarized in the Supplementary Data. Table 2 lists the frequency of rs3746444 variants in miRNA-499. MiRNA-499 rs3746444 AA, AG, and GG variant frequencies were 71.13% (409), 26.78% (154), and 2.09% (12) in HCC patients and 72.64% (669), 24.97% (230), and 2.39% (22) in controls. In the included controls, MAF of rs3746444 in miRNA-499 was similar to the HapMap-CHB database (MAF = 0.15). We identified a null association between rs3746444 in miRNA-499 and the development of HCC (AG vs. AA: p = 0.454; GG vs. AA: p = 0.754, AG/GG vs. AA: p = 0.527; GG vs. AA/AG: p = 0.704). With adjustments for age, sex, HBsAg, smoking, and drinking, we also found that there was no association between rs3746444 in miRNA-499 and HCC risk (AG vs. AA: p = 0.331; GG vs. AA: p = 0.845, AG/GG vs. AA: p = 0.332; GG vs. AA/AG: p = 0.912).
The Relationship of miRNA Polymorphisms with Hepatocellular Carcinoma
Bold values are statistically significant (p < 0.05).
Adjusted for age, sex, smoking, status of chronic HBV infection, and drinking.
CI, confidence interval; HBV, hepatitis B virus; miRNAs, microRNAs; OR, odds ratio.
Table 2 summarizes the frequency of rs2910164 variants in miRNA-146a. The CC, CG, and GG variant frequencies were 33.74% (194), 51.65% (297), and 14.61% (84) in HCC patients and 38.98% (359), 46.91% (432), and 14.12% (130) in controls. In controls, the MAF of miRNA-146a rs2910164 was similar to the HapMap-CHB database (MAF = 0.35). We found that miRNA-146a CG or CG/GG increased the risk of HCC (CG vs. CC: p = 0.039 and CG/GG vs. CC: p = 0.041). With adjustments for age, sex, HBsAg, smoking, and drinking, we also found that there was an association between miRNA-146a rs2910164 and HCC susceptibility (CG/GG vs. CC: p = 0.043).
miRNA-196a2 rs11614913 TT, TC, and CC variant frequencies were 31.48% (181), 48.87% (281), and 19.65% (113) in HCC patients and 31.38% (289), 51.47% (474), and 17.16% (158) in controls. In controls, the MAF of miRNA-196a2 rs11614913 was 0.43. We did not identify any association between rs11614913 in miRNA-196a2 and development of HCC (TC vs. TT: p = 0.650; CC vs. TT: p = 0.393, CT/CC vs. TT: p = 0.968; CC vs. TT/TC: p = 0.223). With adjustments for age, sex, HBsAg, smoking, and drinking, we also found that there was no correlation between rs11614913 in miRNA-196a2 and HCC (TC vs. TT: p = 0.772; CC vs. TT: p = 0.343, CT/CC vs. TT: p = 0.896; CC vs. TT/TC: p = 0.214).
Stratification analysis
Table 3 lists the number of miRNA-146a rs2910164 variants in different subgroups. In the <53-year-old subgroup, a correlation of the rs2910164 locus with HCC risk was found (GG/CG vs. CC: adjusted p = 0.011, GG vs. CC: adjusted p = 0.021 and CG vs. CC: adjusted p = 0.027). The association between miRNA-146a rs2910164 and the risk of HCC was also found in the never-smoking and never-drinking subgroups (never-smoking subgroup: GG/CG vs. CC: adjusted p = 0.011 and CG vs. CC: adjusted p = 0.018; never-drinking subgroup: GG/CG vs. CC: adjusted p = 0.042 and CG vs. CC: adjusted p = 0.039). When we conducted stratification analysis by the status of chronic HBV infection, we identified that miRNA-146a rs2910164 conferred a risk to HCC in individuals with chronic HBV infection (GG/CG vs. CC: adjusted p = 0.035 and CG vs. CC: adjusted p = 0.028). However, a null association was found in other subgroups.
Stratified Analyses Between miRNA-146a rs2910164 C > G Polymorphism and Hepatocellular Carcinoma Risk by Sex, Age, Chronic Hepatitis B Virus Infection, Smoking Status, and Alcohol Consumption
Bold values are statistically significant (p < 0.05).
For miRNA-146a rs2910164 C>G, the genotyping was successful in 575 (98.46%) HCC cases, and 921 (99.78%) controls.
Adjusted for age, sex, smoking, drinking, chronic HBV infection (besides stratified factors accordingly) in a multiple logistic regression model.
HCC, hepatocellular carcinoma.
For miRNA-499 rs3746444 and miRNA-196a2 rs11614913, we did not find any association between these SNPs and HCC risk (Tables 4 and 5).
Stratified Analyses Between miRNA-196a2 rs11614913 T > C Polymorphism and Hepatocellular Carcinoma Risk by Sex, Age, Chronic Hepatitis B Virus Infection, Smoking Status, and Alcohol Consumption
For miR-196a2 rs11614913 T>C, the genotyping was successful in 575 (98.46%) HCC cases, and 921 (99.78%) controls.
Adjusted for age, sex, smoking, drinking, chronic HBV infection (besides stratified factors accordingly) in a multiple logistic regression model.
Stratified Analyses Between miRNA-499 rs3746444 A>G Polymorphism and Hepatocellular Carcinoma Risk by Sex, Age, Chronic Hepatitis B Virus Infection, Smoking Status, and Alcohol Consumption
For miRNA-499 rs3746444 A>G, the genotyping was successful in 575 (98.46%) HCC cases, and 921 (99.78%) controls.
Adjusted for age, sex, smoking, drinking, chronic HBV infection (besides stratified factors accordingly) in a multiple logistic regression model.
Combination analysis of miRNA polymorphisms
Table 6 summarized the combined analyses of miRNA polymorphisms. Rs2910164/rs3746444, rs11614913/rs3746444, rs11614913/rs2910164, and rs11614913/rs2910164/rs3746444 combinations were used to analyze the interaction of SNP-SNP with HCC risk (Qiu et al., 2019; Thakur et al., 2019). Compared with rs2910164 CC/rs3746444 AA, we found that rs2910164 CG/rs3746444 AA may confer an increased risk to HCC (OR 1.43, 95% CI 1.09–1.87, p = 0.010).
Combination Analysis of miRNA Loci (rs3746444, rs2910164, and rs11614913) in Hepatocellular Carcinoma Cases and Controls
Bold values are statistically significant (p < 0.05).
Power calculation and FPRP method
For miRNA-146a rs2910164 SNP, the power value was 0.768 in CG/GG versus CC genetic model. In the <53-year-old subgroup, the power values were 0.837 in CG versus CC, 0.875 in GG versus CC, and 0.928 in CG/GG versus CC. In the never-smoking subgroup, the power values were 0.861 in CG versus CC and 0.842 in CG/GG versus CC. In the never-drinking subgroup, the power values were 0.759 in CG versus CC and 0.730 in CG/GG versus CC. In individuals with chronic HBV infection subgroup, the power values were 0.646 in CG versus CC and 0.616 in CG/GG versus CC. In combination analysis of miRNA polymorphisms, compared with rs2910164CC/3746444AA, rs2910164CG/3746444AA could not promote the development of HCC (power = 0.735). Using FPRP method and power value, we found that miRNA-146a rs2910164 conferred a risk to HCC in the <53-year-old and never-smoking subgroups.
Discussion
Accumulating evidence suggests that miRNAs are implicated in a new mechanism in the development of HCC. Recently, several studies have been performed to identify the relationship of miRNA variants with susceptibility to HCC. However, these studies reported conflicting findings. In this study, we included 584 HCC cases and 923 controls and explored the association of miRNAs variants with HCC susceptibility. We found that rs2910164 increased the susceptibility to HCC in the never-smoking and the <53-year-old subgroups.
Currently, more gene-related investigations have focused on the relationship of rs2910164 C>G polymorphism in miRNA-146a with cancer development. However, most of the polymorphisms only make a low penetrance for cancer risk. In this study, we performed a study involving 1507 subjects to conduct an evaluation between rs2910164 variants and HCC susceptibility. Recently, the correlation of the variants in rs2910164 with HCC risk has been extensively studied (Xu et al., 2008; Akkiz et al., 2011b; Zhang et al., 2011, 2013; Shan et al., 2013; Chu et al., 2014; Qi et al., 2014; Li et al., 2015; Riazalhosseini et al., 2016; Farokhizadeh et al., 2019). In addition, two meta-analyses have identified an association between rs2910164 variants in miRNA-146a and the overall susceptibility to cancer (Ma et al., 2013; Nikolic et al., 2015). However, the potential correlation between the variants in rs2910164 and HCC development is conflicting. Xu et al. (2008) reported that rs2910164 SNP might influence a risk of HCC. However, other investigations suggested that rs2910164 C>G polymorphism in miRNA-146a was not associated with the risk of HCC (Akkiz et al., 2011b; Zhang et al., 2011, 2013; Shan et al., 2013; Chu et al., 2014; Qi et al., 2014; Li et al., 2015; Riazalhosseini et al., 2016; Farokhizadeh et al., 2019).
Considering most of these studies had a small sample size, the association of rs2910164 C > G polymorphism to HCC development was still unclear. In this study, 1507 participants were recruited to evaluate a correlation of rs2910164 SNP with HCC risk. In this study, we identified that the rs2910164 C > G locus might increase the susceptibility of HCC, especially in the never-smoking and the <53-year-old subgroups. It is notable that rs2910164 SNP conferred a risk to chronic infection of HBV (Al-Qahtani et al., 2017; Khanizadeh et al., 2019), and the rs2910164 CC variant altered spontaneous viral clearance of HBV (Khanizadeh et al., 2019). It could be explained that the rs2910164 C > G polymorphism might influence the status of chronic HBV infection and subsequently increased the susceptibility of HCC. In addition, the association could be influenced by interaction of gene and environmental factors. Wang et al. (2018) reported miR-146a rs2910164 G allele could increase the risk of HCC in a more recent meta-analysis, which was similar to our findings.
MiR-146a is a negative regulator of immune response. It was suggested that mRNAs might regulate Toll-like receptors (TLRs) by combining with 3′-untranslated regions (O'Neill et al., 2011). What is more, some bioinformatic investigations indicated an interaction of miR-146a and TLR4 (Yang et al., 2011; Li and Shi, 2013). A previous study found that an elevated miR-146a expression might lead to the drug resistance and milieu of immunosuppression (Khorrami et al., 2017). Compared with rs2910164 CC/GC genotype, the mir-146a expression was higher in rs2910164 GG genotype patients with inflammatory bowel disease (Zhu et al., 2016a). Additionally, in never-smoking lung cancer patients, the expression of miR-146a also increased in patients with rs2910164 GG variant than those with rs2910164 CC/CG variants (Jeon et al., 2014). Thus, miR-146a rs2910164 C > G locus may be implicated in the development of HCC by promoting the expression miR-146a.
Combined influence of three miRNA polymorphisms was explored in this study. More and more publications reported that the interaction of gene–gene may confer the susceptibility to complex disorders. Recently, some case–control studies suggested the influence of combined genotypes of these miRNA SNPs to risk of diseases (Rah et al., 2013; Thakur et al., 2019). In the current investigation, rs2910164 CG/rs3746444 AA was found to be associated with a tendency risk to HCC. In a previous study, Rah et al. (2013) reported that rs2910164 CG/rs3746444 AA might decrease the risk of premature ovarian failure. However, in the other study, the association of rs2910164 CG/rs3746444 AA with risk of cervical cancer was not found (Thakur et al., 2019), which was not consistent with our findings. Therefore, the roles of rs2910164 CG/rs3746444 AA might be different in different diseases. The potential association should be verified in future studies.
In this study, some merits should be acknowledged. First, our investigation included larger sample size to detect the correlation between miRNA variants and the development of HCC. Second, the MAFs of miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 rs11614913 were similar to the database for a Chinese population. Third, the distribution of miRNA variants was consistent with the HWE in controls. Finally, we observed that rs2910164 increased a risk of HCC in the chronic HBV infection subgroup. It is assumed that the rs2910164 G allele may influence spontaneous viral clearance of HBV and is implicated in the susceptibility to HCC subsequently.
Additionally, several limitations in the present investigation also should be addressed. First, due to a hospital-based study, bias might not have been avoided. Second, we only included five major risk factors (e.g., age, sex, HBsAg, smoking, and drinking) in this study. Other potential environmental factors should not have been ignored. Third, in some subgroups, the moderate sample size might limit the power of stratification analysis. Fourth, the association of SNPs in miRNAs with the expression and function was not studied. Fifth, the participants were relatively small in some comparisons, and these results might be underpowered. Finally, we only studied the relationship of miRNA-499 rs3746444, miRNA-146a rs2910164, and miRNA-196a2 with HCC risk. Additionally, we did not explore the association of other SNPs in miRNAs with HCC development.
In summary, our study suggests that miRNA-146a rs2910164 C > G SNP increases the risk of HCC. In the future, more studies with detailed environmental factors and large sample sizes should be carried out.
Footnotes
Acknowledgments
The authors appreciate all subjects who participated in this study. They wish to thank Dr. Yan Liu (Genesky Biotechnologies, Inc., Shanghai, China) for technical support.
Author Contributions
Conceived and designed the experiments: Y.C., L.L.; Performed the experiments: S.Z., L.C., Y.W.; Analyzed the data: Y.W., W.T.; Contributed reagents/materials/analysis tools: Y.C., L.L.; Wrote the article: S.Z., L.C.
Disclosure Statement
No competing financial interests exist.
Funding Information
The project was supported by the Application and Basic Research Project of Changzhou City (CJ20180068).
Supplementary Material
Supplementary Data
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
