1. Context
Hepatitis B virus (HBV) infection is one of the major public health concerns worldwide, with a total of 248 million chronically infected patients (1); it is the main cause of cirrhosis and hepatocellular carcinoma. The outcome of HBV infection is varied: spontaneous clearance, HBV carriers, patients with chronic HBV infection, cirrhosis, or hepatocellular carcinoma, which is influenced by environmental, viral and host genetic factors. HBV infection exhibits a significant family clustering, and relatives with blood relationship are more susceptible to HBV infection than ones without blood relationship (2). In addition, studies by Lin et al. (3) showed the coincidence rate of identical twin HBV infection is higher than that of fraternal twins, which may mean that genetic factors were the important causes of HBV infection.
In recent years, the discovery of sodium taurocholate cotransporting polypeptide (NTCP) was a major breakthrough in the researches of HBV infection mechanism. NTCP is the coding product of SLC10A1, a member of the 10th solute carrier gene family. Early studies suggested that the main function of NTCP is to absorb and transported back to the stem cells at a 2:1 ratio of Na+ reabsorbed by the portal vein to conjugated bile acid, which plays an important role in the intestinal-hepatic circulation of bile acid (4). Recent studies found that NTCP is a functional receptor of hepatocytes infected with HBV, which mediates HBV invasion and infection through specific binding with the pre-S1 antigen of HBV envelope protein (5). The study by Yan et al. (6) revealed that mutation of p.S267F (rs2296651) in SLC10A1 gene can decrease the expression level of NTCP or the affinity between NTCP and HBV; thereby, reducing the infection of HBV to hepatocytes. To date, a series of studies show a correlation between NTCP gene variants and HBV infection, but the results of some studies are contradictory. For instance, the study by Su et al. (7) focusing on Tibetans and Uygurs from West of China showed no association between SLC10A1 (rs2296651) polymorphism and HBV chronicity. In contrast, the investigation by Peng et al. (8) indicated that gene variant of NTCP was significantly correlated with resistance to chronic HBV infection.
2. Objectives
The current study mainly aimed at evaluating the association between NTCP gene polymorphisms and spontaneous clearance of HBV in Asian population by summarizing the results of previous studies.
3. Data Sources
Two independent investigators systematically searched for articles studying the association between NTCP gene polymorphisms and HBV in the following databases: PubMed, Embase, Web of Science, Scopus, Chinese National Knowledge Infrastructure, and Chinese WanFang Data. The latest search update was 09 August 2019. The search strategy was relied on a combination of “NTCP OR sodium taurocholate cotransporting polypeptide OR SLC10A1”, “hepatitis B virus OR HBV OR CHB” and “polymorphism OR variant OR SNP OR gene mutation OR genetic variants”, without any limitation of language or publishing year. In addition, the reference list of eligible articles and relevant reviews were manually searched to find other potential sources of studies.
4. Study Selection
4.1. Inclusion Criteria
Each study was independently assessed by two investigators, and the third party opinions were sought in case of disagreements. The main inclusion criteria in the current meta-analysis were: (1) evaluation of the association between NTCP gene polymorphisms and spontaneous clearance of HBV; (2) reporting the number of genotypes in case and control groups, and providing sufficient information to calculate odds ratio (OR) and 95% confidence interval (CI); (3) selection of HBV persistent infection as the case and spontaneous clearance of infection as the control groups, while the diagnostic criteria were clearly described; and (4) case-control study. The exclusion criteria were: (1) overlapping data, and (2) enrolling patients co-infected with HCV, HDV, or HIV. If there were multiple literature reports on the same population, only the largest study was included in the meta-analysis. The quality of literature was assessed by the Newcastle-Ottawa scale (NOS), which consists of three parts, as the selection of subjects (0 - 4 points), the comparability between groups (0 - 2 points), and the measurement of exposure (0 - 3 points) (Supplementary File Appendix 1).
4.2. Definitions
HBV persistent infection is defined as patients with positive HBsAg for more than six months. HBV spontaneous clearance refers to the ones negative HBsAg, but positive for HBsAb, without a history of vaccination and hepatitis.
5. Data Extraction
For each article included in the analysis, the following data were extracted: the first author, year of publication, country and ethnicity of participants, sample size, age, gender, allele and genotype frequencies of SNP, and genotyping methods. The two researchers extracted data independently and agreed on inconsistent contents through discussions with the third party.
5.1. Statistical Analysis
In the current meta-analysis, pooled OR and 95%CI were performed for dominant (GG vs. GA+AA), recessive (GG+GA vs. AA), codominant (GG vs. AA, GA vs. AA, GG vs. GA) and alleles (G vs. A allele) genetic model to analyze the association between NTCP gene polymorphisms and spontaneous clearance of HBV infection. Since the number of AA genotypes of rs2296651 in some studies was zero in the case and control groups, the association of rs2296651 was only analyzed in dominant, codominant, and alleles models. First, Q-test and I2 statistic were used to assess heterogeneity between studies (P < 0.1 and I2 > 50% indicating significant heterogeneity). Subsequently, based on the results of the heterogeneity, the pooled OR value and the 95% CI were calculated using the random effects model or the fixed effect model. The subgroup analysis was performed based on Hardy-Weinberg equilibrium (HWE) when significant heterogeneity was found or the frequency distribution of the control group gene in the study did not satisfy the HWE (P < 0.1). Next, the sensitivity analysis was used as a method of recalculating the combined OR values after omitting the studies one by one. Publication bias was measured by the funnel plots and Egger regression test. All of the above statistical analysis methods were performed by Stata V. 15.0 software (StataCorp, College Station, TX, USA). Unless otherwise stated, a two-tailed P < 0.05 was considered the level of significance.
6. Results
6.1. Literature Retrieval and Literature Quality Evaluation
A total of 201 literatures were retrieved. According to inclusion-exclusion criteria, nine case-control studies (7-15) published from 2014 to 2019 were finally included in the meta-analysis. The literature selection chart is shown in Figure 1. The meta-analysis included 14165 subjects, including 6400 spontaneous clearance controls and 7765 patients with persistent HBV infection (Supplementary File Appendix 2). In the association study of rs2296651 by Zhang et al. (13) (P = 0.004) and Cheng et al. (14) (P = 0.060), the frequency distribution of the control group gene did not satisfy the HWE, neither did the studies of Wu (11) (P < 0.001) and Su et al. (7, 9) (P = 0.024, P = 0.035) in the rs7154439. The basic characteristics and quality scores of the included studies are shown in Table 1.
Figure 1.
Flow chart of literature search and screening
Table 1.Characteristic of the Included Studies
| Study ID | Country | Ethnicity | SNP | Control | Case | Gender, Male (%, Control/Case) | Age (Control/Case) | Method | NOS |
|---|---|---|---|---|---|---|---|---|---|
| Su, Z. (2014) (9) | Southwestern China | Han | rs4646287 rs7154439 | 352 | 581 | 54.8/70.4 | (50 ± 12)/(47 ± 13)a | TD PCR | 6 |
| Peng, L. (2015) (8) | Southern China | Han | rs2296651 | 90 | 99 | NA | NA | Sanger | 7 |
| Chen, X. (2016) (10) | Central China | Han | rs7154439 | 713 | 1023 | 60.1/64.1 | 37.3[2 - 82]/36.2[3 - 81]b | TaqMan | 7 |
| Su, Z. (2016) (7) | Western China | Tibetans | rs2296651 rs4646287 rs7154439 | 432 | 439 | 53.9/67.4 | 40[33, 55]/39[30, 49]c | TD PCR | 6 |
| Uygurs | 234 | 197 | 59.8/61.9 | 48[39, 58]/29[23, 42]c | |||||
| Wu, W. (2017) (11) | Southern China | Han | rs2296651 rs7154439 | 595 | 586 | 65.0/73.7 | (45.53 ± 11.95)/(45.68 ± 11.53)a | iMLDR | 8 |
| Zhang, Y. (2017) (13) | Southeastern China | Han | rs2296651 rs4646287 rs7154439 | 2116 | 2535 | 74.1/73.1 | (41.3 ± 13.8)/(41.2 ± 12.1)a | Mass Array | 7 |
| Su, M. (2017) (12) | Southern China | Han | rs2296651 rs4646287 | 199 | 305 | 47.2/70.8 | (45.4 ± 11.3)/(42.2 ± 13.9)a | iMLDR | 6 |
| Cheng, X. (2019) (14) | Central China | Han | rs2296651 | 1463 | 1390 | 41.2/39.7 | (40.4 ± 13.9)/(39.5 ± 14.3)a | TaqMan | 7 |
| Chuaypen, N. (2019) (15) | Thailand | Thai | rs2296651 rs4646287 | 206 | 610 | 66.5/69.7 | (51.3 ± 5.9)/(52.6 ± 12.2)a | TaqMan | 7 |
Abbreviations: iMLDR, improved multiplex ligation detection reaction; NA, not available; NOS, Newcastle-Ottawa scale; SNP, single nucleotide polymorphism; TD PCR, touchdown polymerase chain reaction.
aValues are expressed as mean ± SD.
bValues are expressed as mean (range).
cValues are expressed as median (interquartile range).
6.2. Association Between the rs2296651 Polymorphism and Spontaneous Clearance of HBV
Considering the heterogeneity, the association between rs2296651 polymorphism and self-limited clearance of HBV was analyzed using the random effects model. Pooled results indicated that rs2296651 A allele was a protection factor against persistent HBV infection (G vs. A allele: OR = 0.74, 95%CI = 0.58 - 0.94). The results were similar under the dominant model (GG vs. GA+AA: OR = 0.73, 95%CI = 0.58 - 0.93) and codominant model (GG vs. GA: OR = 0.74, 95%CI = 0.59 - 0.93) (Table 2).
Table 2.Association Between SLC10A1 and Spontaneous Clearance of HBV
| SNP | Number of Studies | Genetic Model | P Value | OR (95%CI) | I2 (%) | Q-test | Egger’s P | |
|---|---|---|---|---|---|---|---|---|
| rs2296651 | 7 | Dominant | GG vs. GA+AA | 0.010a | 0.73 (0.58, 0.93) | 51.1 | 0.046 | 0.277 |
| Codominant | GG vs. GA | 0.008a | 0.74 (0.59, 0.93) | 45.2 | 0.078 | 0.216 | ||
| Allele | G vs. A allele | 0.014a | 0.74 (0.58, 0.94) | 56.0 | 0.026 | 0.332 | ||
| rs4646287 | 5 | Dominant | GG vs. GA+AA | 0.150 | 1.10 (0.97, 1.24) | 0.0 | 0.420 | 0.849 |
| Recessive | GG+GA vs. AA | 0.332 | 1.26 (0.79, 2.02) | 14.1 | 0.324 | 0.862 | ||
| Codominant1 | GG vs. AA | 0.304 | 1.28 (0.80, 2.05) | 14.1 | 0.325 | 0.874 | ||
| Codominant2 | GA vs. AA | 0.530 | 1.17 (0.72, 1.89) | 17.0 | 0.304 | 0.755 | ||
| Codominant3 | GG vs. GA | 0.206 | 1.09 (0.96, 1.24) | 0.0 | 0.810 | 0.239 | ||
| Allele | G vs. A allele | 0.116 | 1.10 (0.98, 1.23) | 0.0 | 0.765 | 0.719 | ||
| rs7154439 | 5 | Dominant | GG vs. GA+AA | 0.323 | 1.04 (0.96, 1.14) | 0.0 | 0.940 | 0.795 |
| Recessive | GG+GA vs. AA | 0.109 | 0.72 (0.48, 1.07) | 60.3 | 0.028 | 0.076 | ||
| Codominant1 | GG vs. AA | 0.136 | 0.75 (0.51, 1.10) | 57.1 | 0.040 | 0.080 | ||
| Codominant2 | GA vs. AA | 0.069 | 0.67 (0.43, 1.03) | 63.9 | 0.016 | 0.066 | ||
| Codominant3 | GG vs. GA | 0.134 | 1.07 (0.98, 1.17) | 0.0 | 0.774 | 0.140 | ||
| Allele | G vs. A allele | 0.739 | 1.01 (0.94, 1.09) | 0.0 | 0.718 | 0.293 |
Abbreviation: SNP, single nucleotide polymorphism.
aP < 0.05.
Next, subgroup analysis was performed according to whether the control population met HWE. In the allele model, the analysis results of both subgroups suggested the protective effect of A allele (met HWE: OR = 0.62, 95%CI = 0.39 - 0.97; out of HWE: OR = 0.86, 95%CI = 0.75 - 0.98). Yet, only the subgroups of the studies satisfying HWE showed significant association in the dominant model (GG vs. GA+AA: OR = 0.60, 95%CI = 0.40 - 0.91) and codominant model (GG vs. GA: OR = 0.61, 95%CI = 0.43 - 0.85). No significant heterogeneity existed among studies that were consistent with HWE under the codominant model. In the dominant model, after excluding studies by Zhang et al. (13) and Cheng et al. (14) due to deviation from HWE, the heterogeneity among the studies decreased (I2 decreased from 51.1% to 46.3%) (Figure 2).
Figure 2.
The forest plots for subgroup analysis of the association between rs2296651 polymorphism and spontaneous clearance of HBV. A, B and C, represent the dominant (GG vs. GA+AA), codominant (GG vs. GA) and allele genetic model (G vs. A), respectively. Odds ratios and respective 95% confidence intervals for the different studies included were shown.
6.3. Association Between the rs4646287 Polymorphism and Spontaneous Clearance of HBV
Under all six genetic models, there was no significant heterogeneity among the studies. Thus, the fixed effect model was used for analysis. Pooled results indicated no association between rs4646287 polymorphism and spontaneous clearance of HBV infection (Table 2).
6.4. Association Between the rs7154439 Polymorphism and Spontaneous Clearance of HBV
The correlation between rs7154439 polymorphism and self-limited clearance of HBV was analyzed by the fixed effects model (GG vs. GA+AA, GG vs. GA, G vs. A allele) and the random-effect model (GG+GA vs. AA, GG vs. AA, GA vs. AA) according to the degree of inter-study heterogeneity. No evidence indicated the association between rs7154439 polymorphism and self-limited clearance of HBV under all six genetic models (Table 2).
Subgroup analysis was conducted by recessive model and codominant model 1 and 2 according to whether the HWE was satisfied. In the subgroups of the studies satisfying HWE, there was also no significant association between rs7154439 polymorphism and self-limited clearance of HBV. Results showed no significant heterogeneity within each subgroup (Figure 3).
Figure 3.
The forest plots for subgroup analysis of the association between rs7154439 polymorphism and spontaneous clearance of HBV. A, B and C, represent the recessive (GG+GA vs. AA), codominant1 (GG vs. AA) and codominant2 (GA vs. AA), respectively. Odds ratios and respective 95% confidence intervals for the different studies included were shown.
6.5. Sensitivity Analysis
Sensitivity analysis was performed to assess the impact of each study on the pooled results by omitting each study separately in all genetic models. The results of the analysis of SNP of rs2296651 and rs4646287 indicated that no study had a significant impact on overall results. However, as for recessive model, codominant models 1 and 2 in rs7154439, after excluding studies of Zhang et al. (13), the results were reversed and showed the protective effect of A allele against chronic HBV infection.
6.6. Publication Bias
The potential publication bias was evaluated by funnel plots and Egger regression test. The shape of all funnel plots did not reveal obvious asymmetry and Egger regression also indicated that there is no evidence of potential publication bias in all genetic models (Table 2, Figure 4).
Figure 4.
The funnel plots of the association between NTCP gene polymorphisms and spontaneous clearance of HBV under the dominant genetic model. A, B and C, represent rs2296651, rs4646287 and rs7154439, respectively.
7. Conclusions
A study by Yan et al. (6) found that the NTCP molecule with the S267F (rs2296651) mutation loses its ability to transport bile acids and bind to HBV-PreS1 and HBV. Yet, the effects of this mutation on HBV infection and subsequent liver disease in patients with chronic HBV are not fully understood. According to the results of the current meta-analysis, SLC10A1 rs2296651 A allele was a protective factor against HBV persistent infection in Asian populations. The frequency of genotype GA and AA in patients with persistent HBV infection were higher than those of the spontaneous recovery group, and the difference was statistically significant. Subgroup analysis and sensitivity analysis showed similar results. Furthermore, after omitting the study by Zhang et al. (13) and Cheng et al. (14) due to deviation from HWE, the heterogeneity among the studies decreased under the dominant and codominant model, indicating that the conclusion was robust and deviation from HWE may be one of the source of heterogeneity in the overall analysis. Interestingly, the current mutation of the rs2296651 was only observed in Asian populations with significant racial specificity. For instance, the S267F variant was absent among Moroccans both in healthy controls and patients with chronic HBV patients (16). Perhaps this is one of the reasons for the species-specificity of HBV infection.
The current study also indicated that the frequency of rs4646287 genotype AA in the persistent infection group was slightly higher than that of the control group, but the difference was not statistically significant, which may indicate that rs4646287 gene polymorphism is not associated with susceptibility to HBV. In terms of disease progression, study by Su et al. (9) showed that rs4646287 gene polymorphism was associated with HBV-related cirrhosis and liver cancer.
As for SLC10A1 rs7154439, the current study results for all six models indicated no association between the polymorphism and HBV susceptibility. Furthermore, the combination of recessive, codominant model 1 and 2 in subgroup analysis showed the same results after limiting the study which satisfied HWE, and there was no significant heterogeneity among the studies. Although the deviation from HWE can partly explain the heterogeneity in the overall analysis, sensitivity analysis still suggested that the results were not robust. Given the limited number of studies available for the present meta-analysis, further studies are needed to assess the role of rs7154439 variation in the spontaneous clearance of HBV.
To the best of authors’ knowledge, the current study is the first meta-analysis on the association of NTCP gene polymorphisms with HBV susceptibility. The control population in the current study was spontaneous clearance of HBV with a clear history of HBV virus infection. Yet, the meta-analysis also has the following limitations. First, local databases in other regions were not searched, which may lead to the inclusion of studies mostly focused on Asian populations (especially Chinese). Thus, subgroup analysis cannot be carried out to explore the association between NTCP gene polymorphisms and persistent HBV infection in different ethnical populations. Second, the lack of frequency data combined with genotypes and other environmental factors limited the study ability to further assess gene-environment interactions such as smoking and drinking.
In conclusion, the present meta-analysis results suggested that SLC10A1 rs2296651 polymorphism may play a significant role against persistent HBV infection in Asian populations, while the variation of NTCP gene rs4646287 was not associated with persistent HBV susceptibility. As for rs7154439, more researches are needed to confirm the current study results indicating no association between rs7154439 variant and spontaneous clearance. Further studies investigating the NTCP gene variant in other ethnic populations are also essential.
