Diverse Effects of the NTCP p.Ser267Phe Variant on Disease Progression During Chronic HBV Infection and on HBV preS1 Variability

The sodium taurocholate co-transporting polypeptide (NTCP) acts as a cellular receptor for the hepatitis B virus (HBV) infection on host hepatocytes. We aim to investigate how the NTCP p.Ser267Phe variant affects HBV-related disease progression and analyze viral genomic variability under a host genetic background carrying the p.Ser267Phe variant. A total of 3187 chronic hepatitis B (CHB) patients were enrolled and genotyped for the p.Ser267Phe variant. The variant's association with disease progression was evaluated by logistic regression analysis. We also enrolled 83 treatment-naive CHB patients to analyze the variability of the HBV preS1 region. The frequency of the NTCP p.Ser267Phe variant was significantly lower in patients diagnosed with acute-on-chronic liver failure [OR (95% CI) = 0.33 (0.18–0.58), P = 1.34 × 10⁻⁴], cirrhosis [OR (95% CI) = 0.47 (0.31–0.72), P = 4.04 × 10⁻⁴], and hepatocellular carcinoma [OR (95% CI) = 0.54 (0.34–0.86), P = 9.83 × 10⁻³] as compared with CHB controls under the additive model after adjustment. Furthermore, the percentage of amino acid mutations in HBV preS1 region was significantly higher in the NTCP p.Ser267Phe heterozygote group than in the NTCP wild type homozygote group (P < 0.05). We herein demonstrate that the NTCP p.Ser267Phe variant is a protective factor reducing CHB patient risk for liver failure, cirrhosis, and hepatocellular carcinoma. A host genetic background carrying NTCP p.Ser267Phe exerts selective pressure on the virus, leading to more variability.

Introduction

Hepatitis B virus (HBV) infection constitutes a public health challenge, which affects more than 250 million people worldwide (Trepo et al., 2014). Chronic HBV infection confers a higher risk of hepatitis, liver cirrhosis, and hepatocellular carcinoma (HCC) (Lok and McMahon, 2007). Varying clinical outcomes of HBV infection are affected by a complex combination of virus and host factors, such as viral genomic variability, host immunological state, and host genetic background (Hu et al., 2013; Al-Qahtani et al., 2018; Tan et al., 2018b). The mechanism of HBV-related disease progression has not yet been fully elucidated.

The HBV virus comprises an external envelope composed of surface glycoproteins, an icosahedral nucleocapsid, and a 3.2 kb partially double-stranded DNA genome. The three surface proteins [small (S), middle (M), and large (L)] have a common C-terminal of S domain, but the M and L proteins have different N-terminal domain lengths (Heermann et al., 1984; Seeger and Mason, 2000). HBV manifests significant hepatotropism, and only human and Tupaia hepatocytes are susceptible to HBV infection. Recently, it was discovered that HBV entry into human hepatocytes is mediated by the receptor sodium taurocholate co-transporting polypeptide (NTCP) expressed by the host (Yan et al., 2012; Ni et al., 2014). The preS1 domain of large envelope proteins is responsible for its binding with NTCP and involved in virus–host receptor interaction (Barrera et al., 2005; Glebe et al., 2005; Yan et al., 2012).

NTCP, encoded by SLC10A1, is a sodium–bile acid cotransporter that facilitates their transfer into hepatocytes from portal blood as part of the enterohepatic circulation (Stieger, 2011). Previous studies demonstrated that specific single nucleotide polymorphisms (SNPs) alter the physiological function of NTCP including bile salt homoeostasis, HBV entry, and clinical outcomes of HBV infection (Ho et al., 2004; Peng et al., 2015; Liu et al., 2017). Most SLC10A1 SNPs have distributions related to ethnicity, and the non-synonymous mutation that encodes the p.Ser267Phe variant (S267F, c.800 G>A, rs2296651) is specific to Asian patient populations (Pan et al., 2011; Lee et al., 2017). Our previous work found an association between the chronic hepatitis B (CHB) resistance and the p.Ser267Phe variant (Peng et al., 2015), and several other studies have come to the same conclusion (Hu et al., 2016; Wang et al., 2017; An et al., 2018). Furthermore, some studies found that the Ser267Phe variant is inversely associated with HBV-related disease progression to cirrhosis and HCC (Hu et al., 2016; An et al., 2018). However, several studies did not observe a protection effect of the minor allele of the p.Ser267Phe variant against HBV infection or HBV-related disease progression (Su et al., 2014; Zhang et al., 2017). We, therefore, performed a genetic association study, enrolling a large cohort of 3,187 chronic HBV infected patients to comprehensively assess the relationship between the NTCP p.Ser267Phe variant and the HBV-related clinical outcomes including cirrhosis, liver failure and HCC.

HBV has been evolving together with Homo sapiens over millennia and many studies have reported the evolutionary “arms race” between the HBV and the host immune system (Kramvis et al., 2018; Lumley et al., 2018). Recently, two studies reported that HBV-related hepadnaviruses drive the evolutionary history of their mammalian host receptor and induce signatures of adaptive selection in the cellular receptor NTCP (Jacquet et al., 2018; Takeuchi et al., 2018). Host genes also influence the evolution of the virus. For example, the genetic polymorphisms of the TRIM5 and APOBEC3G genes in primates exert selective pressure on the evolution of simian immunodeficiency virus (Kirmaier et al., 2010; Compton et al., 2012). Previous studies report that the p.Ser267Phe variant of NTCP abolishes or reduces the infection of HBV in vivo and in vitro (Yan et al., 2014; Peng et al., 2015). Therefore, we hypothesize whether the p.Ser267Phe variant of NTCP exerts selective pressure on HBV to drive virus evolution. Considering that preS1 region of HBV is responsible for its binding with NTCP and that this interaction mediates the early viral entry process (Yan et al., 2012), we, therefore, analyzed the variability and conservation of preS1 region in the presence of the host's p.Ser267Phe variant.

Patients and Methods

Cohort and CHB Progression Group Definition

From May 2008 to January 2012, we collected data from 1,899 Han Chinese patients with chronic HBV infection from the Third Affiliated Hospital of Sun Yat-Sen University in Guangzhou, Guangdong Province, China. A total of 1,817 patients from the previous cohort who met the criteria for the present study were included in this study as cohort I, and their clinical data were last updated on April 30, 2018. From January 2016 to April 2018, 1,370 patients with chronic HBV infection were enrolled from the same center as cohort II. A total of 3,187 individuals with chronic HBV infection were enrolled in the project.

Cases in cohorts I and II were consistent with the inclusion and exclusion criteria for chronic HBV infection as in our previous studies (Peng et al., 2015). The diagnosis of chronic HBV infection was based on seropositivity of hepatitis B surface antigen (HBsAg) over 6 months in accordance with the Chinese guideline of prevention and treatment for CHB (Chinese Society of Hepatology and Chinese Society of Infectious Diseases, 2011). Cases who were pregnant, were coinfected with hepatitis A, C, D, E, or human immunodeficiency virus, and were suffering from autoimmune diseases or from alcohol or drug-induced hepatitis or fatty liver disease were excluded. Patients in the CHB control group were screened of laboratory tests (platelet count, international normalized ratio of blood clotting, and alpha-fetoprotein) and medical imaging examinations [hepatic ultrasound, computed Tomography (CT), or magnetic resonance imaging (MRI)] to exclude the presence of cirrhosis, HCC, and liver failure. Cases in the cirrhosis group should be the evidence of decompensation such as image diagnosis of splenomegaly, portal hypertension ascites, hypoproteinemia (<36 g/L), or thrombocytopenia (<1 × 10⁹/L). Patients with liver failure or HCC were excluded from the cirrhosis group. Inclusion criteria for the acute-on-chronic liver failure (ACLF) group were according to consensus recommendations of the Asian Pacific Association for the Study of the Liver (APASL) (Sarin et al., 2009). Hepatic encephalopathy (HE) is defined by the scale provided in the AASLD 2014 Hepatic Encephalopathy Guidelines (AASLD and EASL, 2014), and the HE group only enrolled the overt HE cases. Patients in the HCC group were diagnosed based on confirmed histopathological evidence by two independent pathologists, or at least two imaging examinations (hepatic ultrasound together with CT or MRI) with elevated serum α-fetoprotein level.

This study was carried out in accordance with the guidelines of the 1975 Declaration of Helsinki. Before implementation of the project, approval from the Ethics Committee of the Third Affiliated Hospital of Sun Yat-Sen University was obtained. All subjects were informed and provided written informed consent.

Clinical Parameters

We collected clinical parameters including serum alanine aminotransferase (ALT), aspartate aminotransferase (AST), HBsAg, hepatitis B e-antigen (HBeAg), and HBV DNA loads from the hospital information system, and the test method was described as previously reported (Zhang et al., 2016). The HBV DNA load was monitored prior to antiviral therapy. Other baseline information was collected during each patient's first clinical examination after disease onset.

Genomic DNA Extraction, Sequencing, and Genotyping

Genomic DNA was extracted from whole blood using the Blood Genomic DNA Extraction Kit (Tiangen, Beijing, China) according to the manufacturer's instruction. The genomic DNA was stored at −80°C in a freezer until used. Five target SNPs (rs201339654, rs202018997, rs2296651, rs61745930, and rs759531965) in SLC10A1, which were reported to affect HBV infection or bile acid transport in vitro, were selected for analysis (Ho et al., 2004; Pan et al., 2011; Lou et al., 2014; Fu et al., 2017; Muller et al., 2018). Polymerase chain reactions (PCR) were carried out strictly following a standardized protocol. Information on these five SNPs and their sequencing primers are provided in Tables S1, S2.

We first performed a SNaPshot multiplex assay with the five SNPs identified in 574 DNA samples from cohort II. Apart from rs2296651 G>A, no heterozygous or mutant homozygous individuals of the other four SNPs were found. Therefore, we followed the Sanger sequencing method, consistent with previous studies, to genotype rs2296651 in the remaining 796 samples from cohort II. SNPs were visualized using GeneMapper 4.1 (Applied Biosystems, Foster City, CA, USA).

HBV Viral DNA Isolation and HBV Sequencing

We extracted HBV DNA from 200 ul of patient serum samples using a QIAamp MinElute Virus Spin Kit (QIAGEN, Hilden, Germany), following procedures provided by the manufacturer. A DNA fragment containing the preS1 region was specifically amplified using Q5 Hot Start High-Fidelity Master Mix (New England Biolabs, MA, USA) with the following primers: forward, 5′-AAGGTGGGAAACTTTACGGG-3′; reverse, 5′-TGACAWACTTTCCAATCAATAGG-3′. PCR conditions are as follows: 98°C for 30 s; 98°C for 10 s, 55°C for 20 s, and 72°C for 50 s for 35 cycles, and 72°C for 2 min (final extension). PCR products were sequenced using the BigDye Terminator v3.1 Cycle Sequencing Kit (Invitrogen, CA, USA), and sequencing was performed on an ABI 3730XL DNA analyzer. The HBV sub-genomic fragments were submitted to the NCBI website to identify the genotype (http://www.ncbi.nlm.nih.gov/projects/genotyping/formpage.cgi). We aligned and compared the amplified HBV sequences with the HBV genotype B (GenBank accession no. D00330) or genotype C (GenBank accession no. AB033556) reference sequence to determine preS1 fragment variants.

Data Analysis

Statistical analysis of SNP associations was performed using PLINK software (Purcell et al., 2007). Variation associated with diagnosis was tested using logistic modeling with correction for gender, age, lg HBV DNA, AST, ALT, and HBeAg status. The codominant, dominant, recessive, and additive genetic models were applied to calculate the genotype frequencies. A Bonferroni- adjusted P value < 0.05 was considered statistically significant. Odds ratios (OR) and 95% confidence intervals (95% CI) were also calculated. Linear regression modeling of SNP genotype and lg HBV DNA, AST, ALT, and HBeAg status was applied to estimate the correlation, with correction for diagnosis. Deviation from the Hardy–Weinberg equilibrium (HWE), expectation for any genetic variant was analyzed using the chi-square test. The distributions of point variations at the amino acid level between groups were analyzed using Fisher's exact test. Amino acid mutation sites appearing in more than two sequences were included in the analysis. Between-group differences in qualitative data (gender, HBeAg status) were assessed using chi-square tests, and the quantitative data (age and lgHBV) were assessed using Student's t-tests. A two-tailed P value < 0.05 was considered statistically significant.

Results

Demographic and Clinical Characteristics

A total of 3,187 participants, including CHB controls (n = 1,426), individuals with ACLF (n = 535), individuals with cirrhosis (n = 753), individuals with HCC (n = 473), were enrolled in the study. The demographic and clinical characteristics of the recruited population are summarized in Table 1. Participants had a mean age of 44.22 years and comprised 2,632 men and 555 women. Mean age of each group increased in the order of CHB (35.56 ± 10), ACLF (41.61 ± 10), Cirrhosis (CIR) (46.67 ± 10), and HCC (47.82 ± 10). The genotypes of the p.Ser267Phe variant in the CHB group, CIR group, HCC group, and all participants combined were distributed in line with HWE (P > 0.05). In ACLF group, the allele frequency of the p.Ser267Phe variant was lower than expected and did not conform to HWE (P = 0.017).

TABLE 1

Table 1. Demographic and clinical features in the study groups.

The NTCP p.Ser267Phe Variant and Risk for HBV- Related ACLF

To assess the association between the p.Ser267Phe variant and the HBV-related ACLF, patients with ACLF were compared with CHB controls (Table 2). The minor allele frequency (MAF) of the p.Ser267Phe variant in the ACLF group was significantly lower than that in the CHB controls [OR (95% CI) = 0.31 (0.20–0.49); P = 6.59 × 10⁻⁸]. Under the dominant model, the frequency of GA+AA genotypes in ACLF patients was significantly lower than that in the CHB controls [OR (95% CI) = 0.30 (0.16–0.54); P = 6.05 × 10⁻⁵]. In the additive model, the estimated allelic OR is 0.44 (P = 7.67 × 10⁻⁷). Genetic analysis demonstrated that the p.Ser267Phe variant is inversely correlated with HBV-related ACLF, a result that is consistent with our previous work.

TABLE 2

Table 2. NTCP p.Ser267Phe variant and risk for acute-on-chronic liver failure (ACLF).

HE is a complex neuropsychiatric syndrome that occurs during liver failure (Swaminathan et al., 2018). Recent study has shown that HE is associated with increased blood levels of ammonia and bile acids, especially the conjugated bile acids (Xie et al., 2018). NTCP plays an important role in maintaining homeostasis enterohepatic circulation of bile acids (Hagenbuch and Meier, 1994; Kullak-Ublick et al., 2004). In order to determine whether HE is related to the p.Ser267Phe variant genotype, we further compared 204 patients with overt HE with 332 ACLF patients without HE. In the additive model, the estimated allelic OR is 2.50 with a P value of 0.038, suggesting that NTCP p.Ser267Phe is associated with increased risk of HE in liver failure patients (Table S3).

The Association Between NTCP p.Ser267Phe Variant and Development of Cirrhosis and HCC

Next, we examined whether the presence of the p.Ser267Phe variant affects the development of HBV-related cirrhosis and HCC. Genetic analysis demonstrated that the p.Ser267Phe variant protects CHB patients from HBV-induced cirrhosis progression regardless of whether the model was codominant (OR = 0.46; P = 1.79 × 10⁻³), dominant (OR = 0.46; P = 3.81 × 10⁻⁴), or additive models (OR = 0.47; P = 4.04 × 10⁻⁴; Table 3).

TABLE 3

Table 3. NTCP p.Ser267Phe variant and risk for cirrhosis (CIR).

Since HCC always occurs on the foundation of cirrhosis (Affo et al., 2017), we assessed the association between the NTCP p.Ser267Phe variant and the HCC. Firstly, we compared HCC patients with CHB controls. The results of the dominant (OR = 0.55; P = 0.01) and additive models (OR = 0.54; P = 9.83 × 10⁻³) suggested a 45% reduction in risk for developing HCC (Table 4). Secondly, to determine whether the p.Ser267Phe variant affects risk of developing HCC in individuals with cirrhosis, we compared cirrhosis patients with and without HCC. The results of allele frequency (P = 0.93), dominant model (P = 0.77), and additive model (P = 0.86) did not show statistically significant differences (Table S4). This suggests that the p.Ser267Phe variant does not reduce the risk of cirrhosis developing into HCC.

TABLE 4

Table 4. NTCP p.Ser267Phe variant and development of hepatocellular carcinoma (HCC).

The NTCP p.Ser267Phe Variant and Overall Disease Progression

We sought to determine whether the presence of the p.Ser267Phe mutation might affect the overall disease progression in CHB patients. We combined the ACLF, cirrhosis, and HCC groups together as the overall disease progression group and then compared it with the CHB group. The frequency of the minor A allele of the p.Ser267Phe variant in the overall disease progression group was significantly lower than that in CHB in control group [OR (95% CI) = 0.46 (0.36–0.60); P = 2.56 × 10⁻¹⁰] (Table 5). The estimated ORs of individuals with HBV-related disease progression and CHB controls who are heterozygous (GA) and minor allele homozygous (AA) relative to that of wild type homozygotes (GG) were 0.43 under the dominant model, indicating a protective effect of the A allele against disease progression. However, the recessive model failed to produce a significant difference, probably because only six homozygous AA carriers were enrolled in the study. Overall, we found that the minor allele (A) plays a protective effect in CHB patients, conferring an approximately 70% reduction in risk for developing liver failure (Table 2), a 50% reduction in risk for developing cirrhosis (Table 3), and a 40% reduction in risk for developing HCC (Table 4).

TABLE 5

Table 5. NTCP p.Ser267Phe variant and overall disease progression.

NTCP p.Ser267Phe Variant and Clinical Parameters

Additionally, we tested whether the p.Ser267Phe variant had an impact on clinical parameters including HBeAg state and levels of HBV DNA, ALT, and AST. Under an additive model, the p.Ser267Phe variant is inversely associated with ALT and AST levels (β = −114.6, P = 1.55 × 10⁻⁴ and β = −129.4, P = 3.14 × 10⁻⁴, respectively). However, HBV DNA levels and HBeAg status showed no significant differences (Table S5).

Host Genetic Background of NTCP p.Ser267Phe Variant and HBV preS1 Variability

To determine whether the SLC10A1 genotype exerts selective pressure on HBV, we analyzed the NTCP-interacting domain of the HBV preS1 region. First, we aligned and compared preS1 amino acids (AAs) 1-59 of 2520 genotype B and of 2600 genotype C HBV sequences from an online HBV database (https://hbvdb.ibcp.fr/HBVdb/HBVdbIndex) (Hayer et al., 2013). The results indicate that the preS1 regions of AAs 1-59 in HBV genotype B and C are conserved (Figure 1). Overall preS1 AA composition was visualized using the WebLogo tool (http://weblogo.threeplusone.com/) (Crooks et al., 2004). Next, 83 treatment-naive CHB patients with a viral load >4 lgIU/mL were enrolled in HBV sequence analysis, of which 50 were NTCP wild type homozygotes and 33 were NTCP S267F heterozygotes. The demographic and clinical features of these 83 patients are provided in the Table S6. There are 6 individuals with homozygous p.Ser267Phe mutation in the cohort, and the demographic and clinical features are provided in the Table S7. We can access only the serum from one individual whose HBV DNA viral load is 112 IU/mL, since the other 3 homozygotes are negative for HBV DNA viral load while another 2 homozygotes are lost follow up. The HBV sequences from individuals with homozygous p.Ser267Phe mutation are not included in this study. PreS1 AA 1-59 sequences were aligned and visualized according to the HBV genotype and the NTCP polymorphism status (Figure 1). The view conservation of amino acid sequences was displayed in Figure S1. Four (16.67%) out of 24 genotype B patients with heterozygous S267F mutation harbor the S17A mutation in the preS1 domain, whereas the S17A mutation was not observed in the 39 genotype B patients with wild type NTCP. The percentages of amino acid mutations at positions 17 (which change Ser to Ala) and 56 (which change Asn to Trp) of HBV genotype B were significantly different between the NTCP wild type homozygote and the NTCP S267F heterozygote groups (P = 0.0178 and 0.0178, respectively; Table 6). Two (22.22%) out of 9 genotype C patients with heterozygous S267F mutation harbor the S17A and D47N mutations in the preS1 domain, whereas the S17A and D47N mutations were not observed in the 11 genotype C patients with wild type NTCP. The percentages of amino acid mutations at positions 17 and 42 of HBV genotype C were higher in the NTCP S267F heterozygote group than in the NTCP wild type homozygote group, although these differences were not statistically significant (Table 6). The pooled mutation frequencies of AAs 1-59 in both HBV genotypes B and C are statistically significant higher in the NTCP S267F heterozygote group than in the NTCP wild type homozygote group (1.48 vs. 0.22%, P < 0.001 and 1.89 vs. 0.62%, P = 0.045, respectively). These data indicate that host genetic status at NTCP p.Ser267Phe may exert selective pressure on the HBV virus to drive mutation of the preS1 region.

FIGURE 1

Figure 1. Sequence logos of preS1 amino acids 1–59 from multiple HBV sequences. HBV preS1 amino acids 1–59 sequences were aligned, and the alignment was visualized using the WebLogo tool (http://weblogo.threeplusone.com/). The sequence logos represent calculation of the information content of each amino acid position. 2520 genotype B HBV sequences (A) and 2600 genotype C HBV sequences (D) were obtained from an online HBV database (https://hbvdb.ibcp.fr/HBVdb/). Thirty-nine sequences (B) from the NTCP WT group and 24 sequences (C) from the NTCP S267F group were HBV genotype B. 11 sequences (E) from the NTCP WT group and 9 sequences (F) from the NTCP S267F group were HBV genotype C.

TABLE 6

Table 6. Frequency of amino acid changes in the NTCP-interacting domain of HBV preS1 region.

Discussion

NTCP has been well-established for playing an important role in the enterohepatic circulation of bile acids (Hagenbuch and Meier, 1994; Kullak-Ublick et al., 2004). In 2012, it was identified as a functional cell surface receptor enabling HBV and HDV entries into human hepatocytes in in vitro studies (Yan et al., 2012). Genetic polymorphisms in SLC10A1 have been found to affect NTCP protein function, and some mutations are believed to be specific to certain ethnic backgrounds (Pan et al., 2011). In this study, we screened five SNPs (rs201339654, rs202018997, rs2296651, rs61745930, and rs759531965) in the NTCP gene that have been reported to alter its function in vitro (Ho et al., 2004; Pan et al., 2011; Lou et al., 2014; Fu et al., 2017; Muller et al., 2018). This represents the first demonstration that four of these SNPs are all wild type in 574 Han Chinese patients with CHB. The MAF rate of NTCP p.Ser267Phe variant is 4.52% (288/6374) in Han Chinese patients with CHB in our cohort, and the result is consistent with a previous study showing that the p.Ser267Phe mutation is specific to individuals of Asian descent (Pan et al., 2011). Our previous study has found that this mutation site is associated with resistance to CHB in a Han Chinese cohort of 1,899 CHB patients and 1,828 healthy controls (Peng et al., 2015), which is consistent with studies from Taiwan by Hu et al. (2016) and from Korea by Lee et al. (2017).

In this large cohort study enrolling 3,187 CHB patients, we confirmed the NCTP p.Ser267Phe variant is inversely associated with HBV-related disease progression. We found that the minor allele of p.Ser267Phe variant plays a protective effect in CHB patients conferring an approximately 70% reduction in risk for developing liver failure, 50% reduction in risk for developing cirrhosis, and 40% reduction in risk for developing HCC. Our conclusions are generally consistent with those of previous large cohort studies. Hu et al. reported that p.Ser267Phe mutation is inversely associated with progression to cirrhosis and HCC by comparing patients with cirrhosis or cirrhosis-related HCC to those with neither cirrhosis nor HCC in a community-based CHB cohort (Hu et al., 2016). An et al. reported that the NTCP p.Ser267Phe variant is associated with decreased risk for developing HBV-related cirrhosis and HCC in 1,117 Han Chinese patients with HBV infection (An et al., 2018). Admittedly, we found some inconsistent results from other studies, which may be due to differences in numbers of cases enrolled. Lee et al. did not find that the p.Ser267Phe variant is associated with HCC development in a cohort of 1,376 CHB Korean patients (Lee et al., 2017). Wang et al. did not observe any protective effect of the p.Ser267Phe variant when comparing 117 HCC patients with 866 CHB patients but reported an inverse association between this variant and the HCC when performing a systematic review and meta-analysis including 956 HCC patients and 3,759 CHB patients (Wang et al., 2017). In summary, our findings strongly support a beneficial role for the p.Ser267Phe variant.

The mechanism of the protective effect conferred by the p.Ser267Phe variant has not yet been elucidated. We hypothesize that it may be through reduced HBV infection of hepatocytes and bile acid circulation. Yan et al. (2014) and Liu et al. (2018) reported that the p.Ser267Phe mutation renders substantial loss of function or inefficiency in the ability of NTCP to support HBV infection in vitro. Our previous work demonstrated that the p.Ser267Phe mutation lies within the channel domain for ligand translocation using a structural model (Peng et al., 2015). This mutation may, therefore, reduce HBV transmission and decrease liver inflammation by interfering with NTCP binding to HBV. Regarding bile acids, these are thought to be cytotoxic, and previous studies have suggested that excessive bile acid concentrations can induce hepatocyte injury by activating the death receptor pathway (Faubion et al., 1999). Ho et al. found that the p.Ser267Phe variant of NTCP produces a near-complete loss of function for bile acid uptake from portal circulation into hepatocytes (Ho et al., 2004). This indicates the possibility of reduced intrahepatic cytotoxic bile acid accumulation and a lower likelihood of hepatic inflammation or hepatocyte injury, resulting in a lower risk of HBV-related cirrhosis, HCC, and ACLF disease progression. In our previous work, we found that homozygosity for p.Ser267Phe in SLC10A1 is associated with increased serum bile acid levels and reduced vitamin D levels, which indicates that the bile acid enterohepatic circulation is disturbed and that bile acids in the intestine are insufficient (Liu et al., 2017). Recently, Ma et al. found that bile acid metabolism is mediated by gut microbiome and regulates liver cancer via antitumor immunosurveillance by natural killer T cells38 (Ma et al., 2018). In the present study, we also found that the p.Ser267Phe variant is inversely associated with HE among liver failure patients, which indicates that this mutation might affect the enterohepatic circulation of bile acids or the gut microbiome. However, the association between the NTCP p.Ser267Phe variant and the gut microbiome is not yet clear and requires further study.

Since the p.Ser267Phe mutation is associated with HBV infection and CHB disease progression, whether this mutation exerts selection pressure for viral adaption to interact with the cellular receptor of NTCP needs to be examined. HBV is considered to bind NTCP through the N-terminal pre-S domain of the HBV large envelope protein, and the preS1 domain from Gly-2 to Val-47 (based on HBV genotype D, which correspond to Gly-13 to Val-58 for HBV genotypes B and C) of the N-terminus is regarded as the NTCP binding site (Yan et al., 2012). Two previous studies that enrolled 18 patients from Spain and 12 from Korea reported that the NTCP-interacting domain of the HBV preS1 region shows a high degree of conservation, especially in essential residues for NTCP binding (Lee et al., 2017; Casillas et al., 2018). However, our study found more AA mutation positions, such as S17A, D42N, and N56W (based on HBV genotypes B and C) in the preS1 region among NTCP S267F heterozygotes. It appears that a host genotype encoding NTCP p.Ser267Phe exerts selective pressure on the HBV virus to drive mutation of the preS1 region. The possible mechanism underlying the association between the NTCP S267F mutation and higher preS1 (aa 1-59) mutation rate is still unknown. We hypothesized that the possible pathway may be through bile acid transportation. Verrier et al. reported that bile acid, transported through NTCP, modulates the expression of interferon-stimulated genes to affect HCV Infection (Verrier et al., 2016). When the bile acid transportation by NTCP is inhibited by PreS1 peptides, both the Huh7.5.1-NTCP cells and primary human hepatocytes stimulate the expression of genes involved in the IFN-a response. Previous studies have shown that the p.Ser267Phe mutation abolishes most of its function of bile acid uptake in vitro experiments (Ho et al., 2004; Yan et al., 2014). Several mutation homozygous cases also demonstrate that the presence of NTCP p.Ser267Phe variant is associated with specific changes in circulating bile acids (Deng et al., 2016; Liu et al., 2017; Tan et al., 2018a). The NTCP S267F mutation may regulate HBV infection by modulating the interferon signaling, which is dependent on its function of bile acid transportation and independent of its receptor function. Unlike the in vitro experiment, our data showed that 6 HBV-infected patients (3 for CHB, 2 for ACLF, and 1 for cirrhosis) were homozygotes for NTCP S267F mutation. The NTCP p.Ser267Phe mutation cannot completely block HBV infection in vivo for 8.57% (15/175) of p.Ser267Phe mutation homozygotes were reported to be HBV-infected. In detail, 6 of 32 homozygotes in our cohort, 5 of 42 homozygotes in Hu's study (Hu et al., 2016), 2 of 2 homozygotes in Wang's study (Wang et al., 2017), 1 of 2 homozygotes in An's study (An et al., 2018), 1 of 12 homozygotes in Zhang's study (Zhang et al., 2017), and 0 of 85 homozygotes in Nfor's study (Nfor et al., 2018) were reported to be HBV infection. The HBV-infected p.Ser267Phe mutation homozygotes all had low or negative viral loads. A recent study reported that the HBV strains isolated from a CHB patient homozygous for NTCP S267F exhibit two mutations, S17A and D27E, in the preS1 region (Liu et al., 2018). The authors also demonstrate that the HBV with S17A and D27E mutation uses wild type NTCP as a cellular receptor, and the NTCP p.Ser267Phe variant supports HBV infection with a low efficiency. From these data, it seems likely that mutations in the NTCP contact site in the preS1 domain of the large envelope protein may compensate for the NTCP S267F mutation. Whether mutation of preS1 caused by NTCP polymorphisms leads to higher binding affinity with the NTCP protein needs to be further studied, which may be important to optimize hepatitis B and D virus entry inhibitor drugs.

In the present study, we conducted a comprehensive evaluation of genetic associations of the NTCP p.Ser267Phe variant with CHB disease progression. Meanwhile, we evaluated whether the mutation inflicts evolutionary pressure on the NTCP-interacting domain of the HBV preS1 region. However, there are some limitations to our study. First, the distribution of the minor allele was not in accordance with the Hardy–Weinberg equilibrium (p < 0.001) in the ACLF group, and we did not perform additional population stratification. Second, we found the protective association between the minor allele GA or AA genotypes with HBV-related disease progression. However, our data also indicated that the minor allele was not associated with the HBV DNA levels and HBeAg status aside from the ALT and AST levels. Our study did not provide additional evidence to address the question how does the p.Ser267Phe variation in NTCP alter the HBV-related clinical outcomes. More experiments and evidence are needed to elucidate the role of NTCP and the p.Ser267Phe mutation in HBV-related disease progression. Third, we only analyzed the preS1 region rather than the whole genome of the HBV virus by performing Sanger sequencing. As previously discussed, the preS1 region of HBV is responsible for binding with NTCP, and if there is a selective pressure exerted by the p.Ser267Phe variant, then the preS1 domain is most likely to evolve adaptive selection to compensate for the mutated NTCP. Our data indicated there are more mutations in the HBV preS1 domain in the NTCP p.Ser267Phe heterozygote. Whether these HBV mutations caused by the p.Ser267Phe variant are specifically observed in the preS1 region or are also seen in other fragments need to be further studied. The percentage of wild type and mutation in HBV sequences from “one patient” and the low-frequency mutations were also inaccessible in our data because Sanger sequencing is not precisely quantifiable and is unable to detect alleles below a threshold of 15–20% (Rohlin et al., 2009). In order to address these question, deep sequencing of the HBV preS1 domain is currently underway. After comprehensive analysis of the adaptive evolution of HBV under host genetic background carrying NTCP p.Ser267Phe, it is necessary to demonstrate whether the adaptive variant affects HBV infection in vitro cell culture experiments using HepG2-NTCP cells. It may be interesting to compare the blocking capacities of synthetic preS1 peptide with mutation sites and wild type peptide in inhibiting HBV or HDV entry. At the same time, experiments are needed to assess whether the attachment of preS1 peptide with mutation sites is more robust than the wild type peptide in binding with NTCP protein.

In summary, by assessing all the HBV-related disease progressionas in a large cohort of Chinese CHB patients, the NTCP p.Ser267Phe variant was demonstrated to be protective, reducing the risk of liver cirrhosis, hepatocellular carcinoma, and liver failure in CHB patients. The host genotype encoding NTCP p.Ser267Phe exerts selective pressure on the HBV virus to drive further mutations in the HBV preS1 region. Together with other important risk factors, genetic variation in NCTP and HBV preS1 variability can likely be incorporated into cirrhosis, HCC, and ACLF prediction models to help identify individuals with higher or lower risk of CHB-related advanced liver diseases in Asian areas of endemic HBV infection.

Author Contributions

LP and CX conceived and designed the study. FY, LW, WX, YL, LZ, and GN recruited the cohort and data. JS, QJ, YZ, and TC extracted the DNA and sequenced the genotype. FY and LW analyzed the data and wrote the manuscript. LP and CX critically revised the manuscript. All authors reviewed the manuscript.

Funding

This study was supported by grants from the National Natural Science Foundation of China (No.81873572, 81570539, 81472259, 81101864, 81202319), National Science and Technology Major Project (2018ZX10302204), Guangzhou Science and Technology Plan Projects (201510010292), Science and technology project of Guangdong Province (2015A020212007), and Sun Yat-Sen University Clinical Research 5010 Program (2015004).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We gratefully thank Professor Yiming Wang, Dr. Zhao Qiang and Dr. Ruihong Liu of Sun Yat-sen University for their help in sequencing and genotyping. We also thank Professor Yuehua Huang of the Third Affiliated Hospital of Sun Yat-sen University for help in recruiting the cohort. We thank all the participants of the study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2019.00018/full#supplementary-material

References

Affo, S., Yu, L. X., and Schwabe, R. F. (2017). The role of cancer-associated fibroblasts and fibrosis in liver cancer. Annu. Rev. Pathol. 12, 153–186. doi: 10.1146/annurev-pathol-052016-100322

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Qahtani, A. A., Al-Anazi, M. R., Nazir, N., Abdo, A. A., Sanai, F. M., Al-Hamoudi, W. K., et al. (2018). The correlation between hepatitis B virus precore/core mutations and the progression of severe liver disease. Front. Cell. Infect. Microbiol. 8:355. doi: 10.3389/fcimb.2018.00355

PubMed Abstract | CrossRef Full Text | Google Scholar

American Association for the Study of Liver Diseases (AASLD), European Association for the Study of the Liver (EASL). (2014). Hepatic encephalopathy in chronic liver disease: 2014 practice guideline by the european association for the study of the liver and the american association for the study of liver diseases. J. Hepatol. 61, 642–659. doi: 10.1016/j.jhep.2014.05.042

CrossRef Full Text | Google Scholar

An, P., Zeng, Z., and Winkler, C. A. (2018). The Loss-of-Function S267F variant in HBV receptor NTCP reduces human risk for HBV infection and disease progression. J. Infect. Dis. 218, 1404–1410. doi: 10.1093/infdis/jiy355

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrera, A., Guerra, B., Notvall, L., and Lanford, R. E. (2005). Mapping of the hepatitis B virus pre-S1 domain involved in receptor recognition. J. Virol. 79, 9786–9798. doi: 10.1128/jvi.79.15.9786-9798.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Casillas, R., Tabernero, D., Gregori, J., Belmonte, I., Cortese, M. F., Gonzalez, C., et al. (2018). Analysis of hepatitis B virus preS1 variability and prevalence of the rs2296651 polymorphism in a Spanish population. World J. Gastroenterol. 24, 680–692. doi: 10.3748/wjg.v24.i6.680.

PubMed Abstract | CrossRef Full Text | Google Scholar

Chinese Society of Hepatology and Chinese Society of Infectious Diseases, CMA. (2011) [The guideline of prevention and treatment for chronic hepatitis B (2010 version)]. Zhonghua Gan Zang Bing Za Zhi 19, 13–24. doi: 10.3760/cma.j.issn.1007-3418.2011.01.007.

CrossRef Full Text | Google Scholar

Compton, A. A., Hirsch, V. M., and Emerman, M. (2012). The host restriction factor APOBEC3G and retroviral Vif protein coevolve due to ongoing genetic conflict. Cell Host Microbe 11, 91–98. doi: 10.1016/j.chom.2011.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Crooks, G. E., Hon, G., Chandonia, J. M., and Brenner, S. E. (2004). WebLogo: a sequence logo generator. Genome Res. 14, 1188–1190. doi: 10.1101/gr.849004

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, M., Mao, M., Guo, L., Chen, F. P., Wen, W. R., and Song, Y. Z. (2016). Clinical and molecular study of a pediatric patient with sodium taurocholate cotransporting polypeptide deficiency. Exp. Ther. Med. 12, 3294–3300. doi: 10.3892/etm.2016.3752

PubMed Abstract | CrossRef Full Text | Google Scholar

Faubion, W. A., Guicciardi, M. E., Miyoshi, H., Bronk, S. F., Roberts, P. J., Svingen, P. A., et al. (1999). Toxic bile salts induce rodent hepatocyte apoptosis via direct activation of Fas. J. Clin. Invest. 103, 137–145. doi: 10.1172/jci4765

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, L., Hu, H., Liu, Y., Jing, Z., and Li, W. (2017). Woodchuck sodium taurocholate cotransporting polypeptide supports low-level hepatitis B and D virus entry. Virology 505, 1–11. doi: 10.1016/j.virol.2017.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Glebe, D., Urban, S., Knoop, E. V., Cag, N., Krass, P., Grun, S., et al. (2005). Mapping of the hepatitis B virus attachment site by use of infection-inhibiting preS1 lipopeptides and tupaia hepatocytes. Gastroenterology 129, 234–245. doi: 10.1053/j.gastro.2005.03.090

PubMed Abstract | CrossRef Full Text | Google Scholar

Hagenbuch, B., and Meier, P. J. (1994). Molecular cloning, chromosomal localization, and functional characterization of a human liver Na+;/bile acid cotransporter. J. Clin. Invest. 93, 1326–1331. doi: 10.1172/jci117091

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayer, J., Jadeau, F., Deleage, G., Kay, A., Zoulim, F., and Combet, C. (2013). HBVdb: a knowledge database for Hepatitis B Virus. Nucleic Acids Res. 41(Database issue), D566–D570. doi: 10.1093/nar/gks1022

PubMed Abstract | CrossRef Full Text | Google Scholar

Heermann, K. H., Goldmann, U., Schwartz, W., Seyffarth, T., Baumgarten, H., and Gerlich, W. H. (1984). Large surface proteins of hepatitis B virus containing the pre-s sequence. J. Virol. 52, 396–402.

PubMed Abstract | Google Scholar

Ho, R. H., Leake, B. F., Roberts, R. L., Lee, W., and Kim, R. B. (2004). Ethnicity-dependent polymorphism in Na+;-taurocholate cotransporting polypeptide (SLC10A1) reveals a domain critical for bile acid substrate recognition. J. Biol. Chem. 279, 7213–7222. doi: 10.1074/jbc.M305782200

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, H. H., Liu, J., Lin, Y. L., Luo, W. S., Chu, Y. J., Chang, C. L., et al. (2016). The rs2296651 (S267F) variant on NTCP (SLC10A1) is inversely associated with chronic hepatitis B and progression to cirrhosis and hepatocellular carcinoma in patients with chronic hepatitis B. Gut 65, 1514–1521. doi: 10.1136/gutjnl-2015-310686

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, Z., Liu, Y., Zhai, X., Dai, J., Jin, G., Wang, L., et al. (2013). New loci associated with chronic hepatitis B virus infection in Han Chinese. Nat. Genet. 45, 1499–1503. doi: 10.1038/ng.2809

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacquet, S., Pons, J. B., De Bernardo, A., Ngoubangoye, B., Cosset, F. L., Regis, C., et al. (2018). Evolution of Hepatitis B Virus receptor NTCP reveals differential pathogenicity and species-specificities of hepadnaviruses in primates, rodents and bats. J. Virol. doi: 10.1128/jvi.01738-18. [Epub ahead of print].

CrossRef Full Text | Google Scholar

Kirmaier, A., Wu, F., Newman, R. M., Hall, L. R., Morgan, J. S., O'Connor, S., et al. (2010). TRIM5 suppresses cross-species transmission of a primate immunodeficiency virus and selects for emergence of resistant variants in the new species. PLoS Biol. 8:e1000462. doi: 10.1371/journal.pbio.1000462

PubMed Abstract | CrossRef Full Text | Google Scholar

Kramvis, A., Kostaki, E. G., Hatzakis, A., and Paraskevis, D. (2018). Immunomodulatory function of HBeAg related to short-sighted evolution, transmissibility, and clinical manifestation of hepatitis B virus. Front. Microbiol. 9:2521. doi: 10.3389/fmicb.2018.02521

PubMed Abstract | CrossRef Full Text | Google Scholar

Kullak-Ublick, G. A., Stieger, B., and Meier, P. J. (2004). Enterohepatic bile salt transporters in normal physiology and liver disease. Gastroenterology 126, 322–342. doi: 10.1053/j.gastro.2003.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H. W., Park, H. J., Jin, B., Dezhbord, M., Kim, D. Y., Han, K. H., et al. (2017). Effect of S267F variant of NTCP on the patients with chronic hepatitis B. Sci. Rep. 7:17634. doi: 10.1038/s41598-017-17959-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Xu, G., Gao, Z., Zhou, Z., Guo, G., Li, D., et al. (2018). The p.Ser267Phe variant of sodium taurocholate cotransporting polypeptide (NTCP) supports HBV infection with a low efficiency. Virology 522, 168–176. doi: 10.1016/j.virol.2018.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, R., Chen, C., Xia, X., Liao, Q., Wang, Q., Newcombe, P. J., et al. (2017). Homozygous p.Ser267Phe in SLC10A1 is associated with a new type of hypercholanemia and implications for personalized medicine. Sci. Rep. 7:9214. doi: 10.1038/s41598-017-07012-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Lok, A. S., and McMahon, B. J. (2007). Chronic hepatitis B. Hepatology 45, 507–539. doi: 10.1002/hep.21513

PubMed Abstract | CrossRef Full Text | Google Scholar

Lou, X. Y., Zhang, W., Wang, G., Hu, D. L., Guo, D., Tan, Z. R., et al. (2014). The effect of Na+;/taurocholate cotransporting polypeptide (NTCP) c.800C > T polymorphism on rosuvastatin pharmacokinetics in Chinese healthy males. Pharmazie 69, 775–779. doi: 10.1691/ph.2014.4576

PubMed Abstract | CrossRef Full Text | Google Scholar

Lumley, S. F., McNaughton, A. L., Klenerman, P., Lythgoe, K. A., and Matthews, P. C. (2018). Hepatitis B Virus adaptation to the CD8+; T cell response: consequences for host and pathogen. Front. Immunol. 9:1561. doi: 10.3389/fimmu.2018.01561

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, C., Han, M., Heinrich, B., Fu, Q., Zhang, Q., Sandhu, M., et al. (2018). Gut microbiome-mediated bile acid metabolism regulates liver cancer via NKT cells. Science 360:eaan5931. doi: 10.1126/science.aan5931

PubMed Abstract | CrossRef Full Text | Google Scholar

Muller, S. F., Konig, A., Doring, B., Glebe, D., and Geyer, J. (2018). Characterisation of the hepatitis B virus cross-species transmission pattern via Na+/taurocholate co-transporting polypeptides from 11 New World and Old World primate species. PLoS ONE 13:e0199200. doi: 10.1371/journal.pone.0199200

PubMed Abstract | CrossRef Full Text | Google Scholar

Nfor, O. N., Wu, M. F., Debnath, T., Lee, C. T., Lee, W., Liu, W. H., et al. (2018). Hepatitis B virus infection in Taiwan: the role of NTCP rs2296651 variant in relation to sex. J. Viral Hepat. 25, 1116–1120. doi: 10.1111/jvh.12912.

PubMed Abstract | CrossRef Full Text | Google Scholar

Ni, Y., Lempp, F. A., Mehrle, S., Nkongolo, S., Kaufman, C., Fälth, M., et al. (2014). Hepatitis B and D viruses exploit sodium taurocholate co-transporting polypeptide for species-specific entry into hepatocytes. Gastroenterology 146, 1070–1083. doi: 10.1053/j.gastro.2013.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, W., Song, I. S., Shin, H. J., Kim, M. H., Choi, Y. L., Lim, S. J., et al. (2011). Genetic polymorphisms in Na+;-taurocholate co-transporting polypeptide (NTCP) and ileal apical sodium-dependent bile acid transporter (ASBT) and ethnic comparisons of functional variants of NTCP among Asian populations. Xenobiotica 41, 501–510. doi: 10.3109/00498254.2011.555567

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, L., Zhao, Q., Li, Q., Li, M., Li, C., Xu, T., et al. (2015). The p.Ser267Phe variant in SLC10A1 is associated with resistance to chronic hepatitis B. Hepatology 61, 1251–1260. doi: 10.1002/hep.27608

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

Rohlin, A., Wernersson, J., Engwall, Y., Wiklund, L., Bjork, J., and Nordling, M. (2009). Parallel sequencing used in detection of mosaic mutations: comparison with four diagnostic DNA screening techniques. Hum. Mutat. 30, 1012–1020. doi: 10.1002/humu.20980

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarin, S. K., Kumar, A., Almeida, J. A., Chawla, Y. K., Fan, S. T., Garg, H., et al. (2009). Acute-on-chronic liver failure: consensus recommendations of the Asian Pacific Association for the study of the liver (APASL). Hepatol. Int. 3, 269–282. doi: 10.1007/s12072-008-9106-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Seeger, C., and Mason, W. S. (2000). Hepatitis B Virus biology. Microbiol. Mol. Biol. Rev. 64, 51–68. doi: 10.1128/MMBR.64.1.51-68.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Stieger, B. (2011). The role of the sodium-taurocholate cotransporting polypeptide (NTCP) and of the bile salt export pump (BSEP) in physiology and pathophysiology of bile formation. Handb. Exp. Pharmacol. 201, 205–259. doi: 10.1007/978-3-642-14541-4_5

CrossRef Full Text | Google Scholar

Su, Z., Li, Y., Liao, Y., Cai, B., Chen, J., Zhang, J., et al. (2014). Association of the gene polymorphisms in sodium taurocholate cotransporting polypeptide with the outcomes of hepatitis B infection in Chinese Han population. Infect. Genet. Evol. 27, 77–82. doi: 10.1016/j.meegid.2014.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Swaminathan, M., Ellul, M. A., and Cross, T. J. (2018). Hepatic encephalopathy: current challenges and future prospects. Hepat. Med. 10, 1–11. doi: 10.2147/hmer.s118964

PubMed Abstract | CrossRef Full Text | Google Scholar

Takeuchi, J. S., Fukano, K., Iwamoto, M., Tsukuda, S., Suzuki, R., Aizaki, H., et al. (2018). A single adaptive mutation in sodium taurocholate cotransporting polypeptide induced by hepadnaviruses determines virus species specificity. J. Virol. doi: 10.1128/jvi.01432-18. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, H. J., Deng, M., Qiu, J. W., Wu, J. F., and Song, Y. Z. (2018a). Monozygotic twins suffering from sodium taurocholate cotransporting polypeptide deficiency: a case report. Front. Pediatr. 6:354. doi: 10.3389/fped.2018.00354

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, W., Xia, J., Dan, Y., Li, M., Lin, S., Pan, X., et al. (2018b). Genome-wide association study identi fi es HLA-DR variants conferring risk of HBV-related acute-on-chronic liver failure. Gut 67, 757–766. doi: 10.1136/gutjnl-2016-313035

PubMed Abstract | CrossRef Full Text | Google Scholar

Trepo, C., Chan, H. L., and Lok, A. (2014). Hepatitis B virus infection. Lancet 384, 2053–2063. doi: 10.1016/s0140-6736(14)60220-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Verrier, E. R., Colpitts, C. C., Bach, C., Heydmann, L., Zona, L., Xiao, F., et al. (2016). Solute carrier NTCP regulates innate antiviral immune responses targeting hepatitis C virus infection of hepatocytes. Cell Rep. 17, 1357–1368. doi: 10.1016/j.celrep.2016.09.084

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, P., Mo, R., Lai, R., Xu, Y., Lu, J., Zhao, G., et al. (2017). Genetic variations of NTCP are associated with susceptibility to HBV infection and related hepatocellular carcinoma. Oncotarget 8, 105407–105424. doi: 10.18632/oncotarget.22211

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, G., Wang, X., Jiang, R., Zhao, A., Yan, J., Zheng, X., et al. (2018). Dysregulated bile acid signaling contributes to the neurological impairment in murine models of acute and chronic liver failure. EBioMedicine 37, 294–306. doi: 10.1016/j.ebiom.2018.10.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, H., Peng, B., Liu, Y., Xu, G., He, W., Ren, B., et al. (2014). Viral entry of hepatitis B and D viruses and bile salts transportation share common molecular determinants on sodium taurocholate cotransporting polypeptide. J. Virol. 88, 3273–3284. doi: 10.1128/JVI.03478-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, H., Zhong, G., Xu, G., He, W., Jing, Z., Gao, Z., et al. (2012). Sodium taurocholate cotransporting polypeptide is a functional receptor for human hepatitis B and D virus. Elife 3:e00049. doi: 10.7554/eLife.00049

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, G. L., Zhang, T., Ye, Y. N., Liu, J., Zhang, X. H., Xie, C., et al. (2016). Complement factor 3 could be an independent risk factor for mortality in patients with HBV related acute-on-chronic liver failure. Biomed Res. Int. 2016:3524842. doi: 10.1155/2016/3524842

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Li, Y., Wu, M., Cao, P., Liu, X., Ren, Q., et al. (2017). Comprehensive assessment showed no associations of variants at the SLC10A1 locus with susceptibility to persistent HBV infection among Southern Chinese. Sci. Rep. 7:46490. doi: 10.1038/srep46490

PubMed Abstract | CrossRef Full Text | Google Scholar