Biochemical and Bioinformatic Characterization of Patients with a Sodium Taurocholate Cotransporting Polypeptide Mutation

authors:

avatar Xuan Li 1 , avatar XueQian Zhou 1 , avatar Hong Yang 1 , avatar LiangJun Zhang 1 , avatar XiaoXun Zhang 1 , avatar Jin Chai 1 , *

Department of Gastroenterology; Institute of Digestive Diseases of PLA; Cholestatic Liver Diseases Center and Center for Metabolic Associated Fatty Liver Disease, The First Affiliated Hospital (Southwest Hospital), Third Military Medical University (Army Medical University), Chongqing 400038, China

how to cite: Li X, Zhou X, Yang H, Zhang L, Zhang X, et al. Biochemical and Bioinformatic Characterization of Patients with a Sodium Taurocholate Cotransporting Polypeptide Mutation. Hepat Mon. 2022;22(1):e121842. https://doi.org/10.5812/hepatmon-121842.

Abstract

Background:

SLC10A1 codes for the sodium taurocholate cotransporting polypeptide (NTCP). The SLC10A1S267F mutation is associated with loss of function of bile acid (BA) uptake and defined as a new type of hypercholanemia. This kind of hypercholanemia is characterized by high levels of serum BA. However, limited studies have been conducted on this topic.

Objectives:

This study aimed to describe the biochemical and bioinformatic characterization of patients with an SLC10A1S267F mutation, as well as to dissect pathogenesis in hypercholanemia.

Methods:

In this study, a total of 12 individuals (including 5 homozygous, 3 heterozygous, and 4 wild-type individuals) were recruited. Whole-genome sequencing (WGS) and Sanger sequencing were used to confirm the genotype. Tests of liver function, renal function, and serum lipid level, in addition to routine blood tests, were performed to evaluate the clinical consequences of patients with an SLC10A1S267F mutation. The ClinVar website and protein prediction tools were used to analyze other cholesterol and BAs related gene mutations in SLC10A1S267F patients, as well as to evaluate their possible effects on serum BA levels of patients.

Results:

All SLC10A1S267F homozygous patients displayed high levels of BAs. Liver and renal functions were generally normal. According to previous reports, homozygous patients are prone to vitamin D deficiency and deviated blood lipids. However, all homozygous individuals had normal levels of blood lipids, thyroid hormones, and vitamin D (25(OH)D). Moreover, except for the SLC10A1S267F mutation, according to the WGS results, multiple gene mutations were found in 5 homozygous and might affect the level of BAs, but the SLC10A1S267F mutation still is the most important reason resulting in a high level of BAs.

Conclusions:

This study provided a more detailed description of the SLC10A1S267F mutation-induced hypercholanemia, delivering a new idea that there might be some mutations in SLC10A1S267F homozygotes, probably influencing BA metabolism.

1. Background

Bile is mainly composed of bile salts, phospholipids, and cholesterol. Bile salt is a kind of sodium salt of the conjugate of bile acid (BA) with either glycine or taurine. The conversion of cholesterol to BAs is closely related to lipid metabolism, fat solute vitamins, and hormones (1). BA levels differ in different liver diseases, and the level of BA in the liver reflects and influences liver function to a certain extent. However, the long-term clinical consequences of high serum BA levels remain unknown.

Hepatic uptake of BA is mainly mediated by sodium-dependent or independent transporters. The sodium taurocholate cotransporting polypeptide (NTCP; gene symbol SLC10A1 in humans) is a major transporter of conjugated BA from the serum compartment into hepatocytes (2, 3). Besides, NTCP also transports additional substrates, such as thyroid hormones and drugs (4, 5).

Mutations in the SLC10A1 gene, whose cytogenetic location is on chromosome 14 (14q24.1), may alter the expression and/or function of the NTCP protein, thus influencing bile salts transportation (6-10). The single nucleotide variant (p.Ser267Phe/c.800C>T) of NTCP/SLC10A1 is associated with high levels of BA. The minor allele frequency (MAF) of SLC10A1S267F in East Asia was reported to be between 8% and 12% (9, 11, 12). The SLC10A1S267F variant exhibits a near-complete loss of function for bile salt uptake; yet, it possesses a fully normal transport function for the non-BA substrates such as estrone sulfate (8). This indicates that Ser267 may be placed in a protein region specific for BA transport. Therefore, any clinically impactful changes in the SLC10A1S267F variant may affect liver detoxification or increase susceptibility to hepatotoxicity in patients.

According to recent references in which only a total of 48 individuals, including 38 children and 10 adults (13-19), have been reported as SLC10A1S267F homozygous patients, the mutation of SLC10A1S267F is probably related to elevated serum BAs (a condition often referred to as hypercholanemia). However, the biochemical description of this mutation and whether it is the only reason influencing BAs in patients remain unclear. In this study, we identified 5 SLC10A1S267F homozygous patients and investigated the clinical manifestations and long-term consequences. Here, we also analyze if there exists a combination of multiple gene mutations to influence serum BA levels in SLC10A1S267F patients.

2. Methods

2.1. Human Subject

A total of 12 Han Chinese individuals were recruited at the First Affiliated Hospital of Army Medical University. We identified that 5 individuals with high serum levels of BAs were SLC10A1S267F homozygous, while the other 7 individuals were heterozygous (n = 3) and wild-type (n = 4). We collected basic clinical information and peripheral blood samples from the individuals. Lab examinations, including liver and renal function tests, serum lipids level tests, and routine blood tests, were performed for each subject. This study was approved by the Ethics Committee of the First Affiliated Hospital of Army Medical University, PLA (People's Liberation Army of China) (code: KY2021116). All the study subjects or their legal guardians signed written informed consent.

2.2. Whole-Genome Sequencing

We performed a whole-genome sequencing (WGS) analysis of 12 individuals in this study. Using TruSeq Nano DNA LT Sample Preparation Kits (Illumina Inc, San Diego, CA, USA), the libraries were constructed. Briefly, using an S220 Focused-Ultrasonicator (Covaris, USA), the genomic DNA was sheared into fragments with a length of ~350 base pairs (bp). Using an Illumina sequencing platform, a HiSeq X Ten platform (Illumina Inc, San Diego, CA, USA), the final libraries were sequenced after polymerase chain reaction (PCR) amplification and purification, and 150-bp paired-end reads were generated. The raw sequence data were processed using the NGS-QC toolkit (20) for quality control of sequence data. The clean reads were aligned to the reference genome (GRCh38/hg19) using the Burrows-Wheeler Aligner (21) (BWA), and the average coverage is 29.34X. After alignment, Picard was employed to mark duplicate reads, and SAMtools (22) was used to transform the reads into compressed binary alignment map (BAM) files for sorting and indexing. The genome analysis toolkit (GATK) (23) were used to verify the accuracy of single nucleotide polymorphisms (SNP) and insertion and deletion (InDel). Finally, ANNOVAR (24) was used to perform annotation for the variants.

2.3. Gene Variation Screening

After completing the annotation and filtering steps, we focused on genes known to code for cholesterol and BA metabolism. Firstly, we used “bile acid” and “cholesterol” as keywords and selected all genes coding for human cholesterol metabolism, especially for BA metabolism, using the UniProt database and obtained a total of 560 genes. We further retrieved the SNP file from the WGS results and obtained 709 SNPs (except for SLC10A1S267F) from the recruited subjects. Because homozygotes have greater phenotype and clinical importance, we decided to narrow down the number of SNPs by excluding the heterozygous mutation, meaning that those SNP genotypes should be homozygous in individuals 1 to 5. A total of 560 genes were filtered and reduced to 62 genes covering 76 missense variants. Secondly, to evaluate benign or deleterious of those variants, we searched for their clinical significance in ClinVar or from previous studies and classified them as benign or deleterious variants. Meanwhile, for mutations not previously reported, we used MAF < 0.05 in the public database of 1000 genomes project for the East Asian population and protein prediction tools, including sorting tolerant from intolerant (SIFT; dbNSFP version 3.3a) and polymorphism phenotyping version 2 (Polyphen 2 version 2.2), to forecast the underlying functional impacts. The protein function prediction of SIFT is based on the SIFT score calculated by sequence homology and the physical/chemical similarity between amino acids (deleterious means sift ≤ 0.05; tolerated means sift > 0.05). On the other hand, we used the polyphen 2 HDIV score for protein function prediction, which is based on the HumanDiv database for complex diseases (probably damaging: score ≥ 0.957; possibly damaging: Score = 0.453 ~ 0.956; benign: score ≤ 0.452). Besides, we examined the pathogenesis of the variants reported in all related studies. There were 31 possibly damaging variations that might be related to hypercholanemia (Table 1).

Table 1.

Variations in Bile Acids and Cholesterol Metabolism Genes, Detected by Whole - Genome Sequencing

GenersVariationClinical SignificanceMAFFunctional Effects aGenotype of Individual 1/2/3/4/5/6/7/8/9/10/11/12
ABCA12rs7560008c.T1375Ap.S459TBenign0.0016ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
ABCB11rs2287622c.T1331Cp.V444AProgressive familial intrahepatic cholestasis> 0.01ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HE/HO
ABCC2rs927344c.A116Tp.Y39FDubin - Johnson syndrome0.0054ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
ADGRG6rs1262686c.A3296Gp.Q1099RBenign0.0004DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
AKR1C4rs4880718c.A749Gp.Q250R/0.0002ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
APOBrs679899c.C1853Tp.A618VFamilial hypercholesterolemia, Familial hypobetalipoproteinemia> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
APOBrs584542c.A6937Gp.I2313V/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HE/HO
APOL1rs2239785c.G394Ap.E132K/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HE/HO/HE
APOL1rs136175c.G630Ap.M210I/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HE/HO/HE
ATP11Crs2491014c.T342Gp.C114W/0.0079ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
ATP8B1rs222581c.G3454Ap.A1152T/0.0002DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
CD44rs1467558c.T689Cp.I230T/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
CTAGE5rs7140561c.T17Cp.V6A/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HE/HE/WT/HE/HE/HE
CTAGE5rs1950952c.G991Cp.E331Q/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
CUBNrs2796835c.C8150Gp.S2717WBenign0.0000ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
CUBNrs1276712c.G6485Ap.C2162YBenign0.0066ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
DBTrs12021720c.A1150Gp.S384GBenign> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
EGFrs4698803c.A2633Tp.E878VBenign> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
FCGBPrs782538403c.G3650Ap.R1217Q/0.0020ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
FCGBPrs782416774c.A3700Gp.T1234A/0.0048ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
FCGBPrs782342257c.G3707Ap.R1236Q/NoneDeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
GCrs9016c.A1334Gp.H445R/0.0072DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
LIPCrs6083c.A644Gp.N215SHepatic lipase deficiency> 0.01ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
LIPCrs3829462c.C1068Ap.F356LHepatic lipase deficiency> 0.01ToleratedHO/HO/HO/HO/HO/HO/HO/HE/HE/HO/HE/HO
PCSK9rs562556c.G1420Ap.V474IFamilial hypercholesterolemia, Familial hypobetalipoproteinemia> 0.01ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
PEX16rs10742772c.G346Ap.V116IBenign0.0000ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
PROM2rs12992066c.A1523Gp.Q508R/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HE/HO/HO/HO/HE/HE
SLC14A1rs11877062c.C10Tp.R4W/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HE/HE/HO
TFrs2692696c.A1342Gp.I448V/0.0068ToleratedHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
TTC39Brs1407977c.A1051Gp.I351V/> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO
VPS33Brs11073964c.G1459Ap.G487SArthrogryposis with renal dysfunction and cholestasis syndrome> 0.01DeleteriousHO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO/HO

2.4. Sanger Sequencing Validation

Sanger sequencing was used to confirm the SLC10A1S267F mutation in all recruited individuals. PCR primers were designed and used to amplify the 309/419-bp fragment covering the SLC10A1 gene (forward primer, 5’-GCTAGAAACTTGCTTGTTG-3'; reverse primer, 5’-CTCTGAGTGTATGTGGGGT-3’ or forward primer, 5’-GAGTGCAGTTGATGGAAGT-3’; reverse primer, 5’-CTGAGTGTATGTGGGGTTT-3’). The SLC10A1 gene information was obtained from the National Center for Biotechnology Information database.

2.5. Statistical Analysis

Data were analyzed using SPSS version 23 (SPSS Inc, Chicago, Ill, USA). Continuous variables with normal distribution were presented as mean ± SEM; non-normal variables were reported as median (interquartile range). A comparison of the clinical laboratory results among the 3 study groups (wild-type, heterozygote, and homozygote) was analyzed by either parametric or non-parametric tests (1-way analysis of variance [ANOVA] or Kruskal-Wallis omnibus test) as appropriate. P values less than 0.05 were considered statistically significant.

3. Results

3.1. Description of the Phenotype

We recruited a total of 12 subjects and divided them into 3 groups; homozygote, heterozygote, and wild-type (Figure 1; Table 2). Telephone follow-ups revealed good cognition competence and a normal social ability for both adults and children. Adults of homozygotes and heterozygotes had well-developed secondary sexual characteristics and were fertile. Consistent with previous studies (13-19), the level of serum BA was much higher in homozygotes than in wild-type individuals (P = 0.027), meaning that all 5 homozygotes had hypercholanemia, whereas heterozygotes and wild-type individuals had normal serum BA levels. Many factors could affect serum BA levels, needing to be fully excluded before diagnosing SLC10A1S267F mutation-induced hypercholanemia. Thus, a comprehensive medical history was obtained from each homozygote to determine any other diseases that could lead to high levels of serum BA. However, we did not find any individuals with symptoms of cholestatic liver diseases/malabsorptive diseases or long-term medication history.

Bioinformatic analysis of the SLC10A1 gene mutation. A, schematic view of NTCP exons and multiple sequence alignment information of NTCP; B, three - dimensional structures, prediction of NTCP wild - type and mutation; C, the Sanger sequencing results of the SLC10A1 gene in 9 individuals. Individuals 1 to 5 were homozygous for the SLC10A1 variants c.800c> T (p.Ser267Phe), individuals 6 to 8 were heterozygous for the SLC10A1 variants c.800c> T (p.Ser267Phe), and individual 9 was wild - type for the SLC10A1 variants c.800c> T (p.Ser267Phe). Individuals 10 to 12 were not included in this analysis due to the limited peripheral blood sample.
Bioinformatic analysis of the SLC10A1 gene mutation. A, schematic view of NTCP exons and multiple sequence alignment information of NTCP; B, three - dimensional structures, prediction of NTCP wild - type and mutation; C, the Sanger sequencing results of the SLC10A1 gene in 9 individuals. Individuals 1 to 5 were homozygous for the SLC10A1 variants c.800c> T (p.Ser267Phe), individuals 6 to 8 were heterozygous for the SLC10A1 variants c.800c> T (p.Ser267Phe), and individual 9 was wild - type for the SLC10A1 variants c.800c> T (p.Ser267Phe). Individuals 10 to 12 were not included in this analysis due to the limited peripheral blood sample.
Table 2.

Clinical Characteristics of Individuals

VariablesHomozygotesHeterozygotesWild-Type
Individual123456789101112
SexMaleMaleFemaleFemaleFemaleFemaleMaleFemaleFemaleMaleFemaleMale
Age (y)36211753633319 months434417
TBA (umol/L)33.173435.2471.94.9287.72.83.61.79.1

We compared liver and renal functions between the 3 groups (Table 3). The levels of biochemical markers of liver injury were normal, except for high alkaline phosphatase activity in children because of the accelerated growth. The globulin levels decreased but not significantly. Also, the A/G ratio increased in the homozygotes (P = 0.036) and heterozygotes (P = 0.039) but with no clinical phenotype. The level of total bilirubin was lower in homozygotes than in wild-type individuals (P = 0.014) and heterozygotes but still within the normal range; the decreased total bilirubin might influence iron absorption in the long term in homozygotes.

Table 3.

Biochemical Parameters of the Study Cohort Grouped by Genotype a

Genotype (rs2296651)Wild-Type (GG) (n = 4)Heterozygous (GA) (n = 3)Homozygous (AA) (n = 5)
TestRef.Unit
Liver/Renal Function
ALT0 - 42(IU/L)14.87 ± 2.7621.00 ± 18.2520.00 ± 11.77
AST0 - 42(IU/L)17.70 (17.65, 23.35)16.00 (15.50, 21.50)30.00 (26.00, 42.00)
ALP34 - 114(IU/L)192.33 ± 199.9559.33 ± 32.32190.20 ± 90.73
GGT5 - 50(IU/L)13.67 ± 3.2120.67 ± 9.0720.67 ± 7.20
TBA0 - 10(umol/L)2.80 (2.25, 3.20)4.90 (3.40, 5.95)35.20 (34.00, 47.00) b
TBIL6 - 21(umol/L)12.70 ± 3.9311.03 ± 1.376.82 ± 2.17 b
DBIL0 - 6(umol/L)2.57 ± 0.674.50 ± 1.842.98 ± 0.89
TP66 - 83(g/L)72.83 ± 0.9072.30 ± 1.9370.40 ± 5.69
Alb38 - 51(g/L)46.60 (45.40, 47.20)45.20 (45.0, 46.6)49.00 (48.80, 49.00)
G25 - 38(g/L)26.63 ± 2.7224.00 ± 5.1920.84 ± 4.09
A/G1.2 - 2.51.76 ± 0.251.77 ± 0.17 b2.45 ± 0.49 b
UA155 - 428(umol/L)374.67 ± 124.64286.33 ± 138.60277.80 ± 87.42
Urea1.7 - 8.3(mmol/L)3.97 ± 1.073.81 ± 0.633.12 ± 1.26
Cre45 - 84(umol/L)61.40 (46.20, 61.6)43.00 (42.50, 57.35)36.00 (21.00, 38.00)
Blood Fat
TCH3.1 - 5.7(mmol/L)4.63 ± 0.304.00 ± 1.144.16 ± 0.90
TG0.4 - 1.73(mmol/L)1.28 (1.02, 2.185)0.76 (0.53, 1.56)1.07 (0.78, 1.43)
HDL - C0.9 - 2(mmol/L)1.49 ± 0.521.44 ± 0.391.52 ± 0.31
LDL - C2.07 - 3.1(mmol/L)2.73 ± 0.282.36 ± 0.812.64 ± 0.91
Serum Thyroid Hormone
T31.3 - 3.1(nmol/L)1.35 (1.30, 1.58)1.70 (1.65, 1.76)2.09 (1.86, 2.16)
T466 - 181(nmol/L)98.70 (94.41, 121.22)91.43 (89.87, 94.35)110.05 (104.74, 144.10)
FT33.1 - 6.8(pmol/L)2.97 ± 0.702.57 ± 0.093.28 ± 0.52
FT412 - 22(pmol/L)15.82 (15.69, 16.55)16.49 (14.96, 17.09)21.49 (18.22, 23.09)
TSH0.27 - 4.20(uIU/mL)0.78 ± 0.331.01 ± 0.390.56 ± 0.19
Vitamin D
Vitamin D (25(OH)D)25 - 200(pmol/L)50.60 (47.18, 55.02)45.06 (42.86, 47.55)32.26 (30.60, 34.57)
Routine Blood Tests
WBC3.5 - 9.5(109/L)7.78 ± 1.005.69 ± 0.647.48 ± 2.56
NEU1.8 - 6.3(109/L)4.47 ± 0.813.33 ± 0.712.62 ± 0.88 b
LYM1.1 - 3.2(109/L)2.79 ± 0.301.79 ± 0.10 b3.97 ± 2.74
MON0.1 - 0.6(109/L)0.35 ± 0.160.39 ± 0.010.51 ± 0.21
BAS0 - 0.06(109/L)0.00 (0.00.0.11)0.01 (0.005, 0.01)0.02 (0.01, 0.03)
EOS0.02 - 0.52(109/L)0.13 (0.09.0.21)0.20 (0.12, 0.23)0.15 (0.05, 0.17)
RBC4.3 - 5.8(1012/L)4.88 ± 0.334.95 ± 0.504.79 ± 0.56
HGB115 - 150(g/L)141.00 ± 9.64141.00 ± 24.02127.80 ± 17.58
MCV82 - 100(%)86.57 ± 5.2688.13 ± 12.6386.04 ± 13.52
MCH27 - 34(pg)28.93 ± 1.8728.43 ± 4.8726.92 ± 3.71
MCHC316 - 354(g/L)333.00 (333.00, 335.00)318.00 (313.00, 331.00)311.00 (309.00, 317.00)
RDW - SD37 - 54(%)41.20 ± 1.0146.07 ± 3.56 b43.58 ± 1.82
RDW - CV11 - 16(%)13.00 (12.75, 13.50)13.70 (13.25, 14.95)14.10 (14.00, 14.40)
PLT125 - 350(109/L)287.00 ± 51.47195.67 ± 17.67 b310.40 ± 25.15 c
MPV9 - 13(fl)9.70 ± 0.269.47 ± 1.029.66 ± 0.79
PDW9 - 17(fl)11.50 ± 0.1016.13 ± 0.90 b15.66 ± 0.52 b

Apart from the liver and renal functions, we also determined the levels of blood lipids (including total cholesterol, triglyceride, high-density lipoprotein cholesterol [HDL-C], and low-density lipoprotein cholesterol [LDL-C]), thyroid hormones (T3, T4, FT3, FT4, TSH), and vitamin D (25(OH)D). Previous studies (13, 14, 16) have reported that patients with a homozygous SLC10A1S267F mutation were prone to vitamin D deficiency and abnormal blood lipid levels. However, in our study, all homozygotes had normal blood lipids, thyroid hormones, and vitamin D levels (Table 3). Compared to heterozygotes and wild-type subjects, in homozygotes, some indices (such as T3, T4, FT3, and FT4) increased while TSH decreased; however, all parameters were within the reference range.

The clinical chemistry parameters are presented in Table 1. The neutrophil count in homozygotes decreased (P = 0.015), whereas the lymphocyte and monocyte counts were slightly elevated. Some biochemical markers, including platelet (PLT) and platelet distribution width (PDW) in homozygotes, significantly differed compared with the heterozygotes (P = 0.001) and wild-type subjects (P < 0.001), respectively–but with a small clinical difference (Table 3). Interestingly, in heterozygotes, the lymphocyte count and PLT significantly decreased (P = 0.046 and P = 0.09), and PDW increased (P < 0.01). Taken together, these findings suggest that individuals with or without an SLC10A1S267F mutation are healthy. However, the small differences in results should be studied further.

3.2. Validation of the SLC10A1S267F Mutation Using WGS and Sanger Sequencing

WGS and Sanger sequencing were performed to confirm the presence of the SLC10A1S267F mutation in the subjects. Individuals 1 to 5 were homozygotes, individuals 2 to 7 were heterozygotes at the variation allele, and individuals 9 to 12 were wild-type (Figure 1B). The Sanger sequencing results were consistent with WGS results.

3.3. Relationship of the SLC10A1S267F Mutation and Other Gene Mutations in Homozygotes

We noted that in some cases, SLC10A1S267F homozygotes with a combination of other gene mutations or diseases showed different symptoms, such as jaundice and/or liver dysfunction, late fetal loss, or urticarial rashes (13, 15, 17, 25, 26). We confirmed that serum BA levels increased asymptotically in SLC10A1S267F homozygotes. Approximately 90% of cholesterol would be converted into BA by a series of enzymes and reactions in the liver. Therefore, cholesterol metabolism-related gene mutations could affect BA metabolism and might influence serum BA levels. To confirm if other variants could influence serum BA levels in homozygotes, we performed further gene analysis on these individuals. A total of 62 genes (coding for cholesterol metabolism and/or BA metabolism) covering 76 SNP were selected. The ClinVar public database and MAF were used to evaluate those variants. We also used protein prediction tools to predict underlying functional impacts in metabolism and find out deleterious variants. After filtering, 31 possibly damaging SNP were identified (see Table 1). According to references and prediction tools, except for the SLC10A1S267F mutation, there are 5 SNPs that have been thought to be the most possible factors in influencing serum BA levels in recruited homozygotes (ABCB11V444A, ABCC2Y39F, APOBA618V, PSCK9V474I, and VPS33BG487S). None of the previous studies about SLC10A1S267F have mentioned these SNPs, and SLC10A1S267F in homozygotes was thought to be the most likely cause of high BA levels. Although we did not figure out the relationship between the 5 SNPs and hypercholesterolemia in this study, these variations will be studied in our future research.

4. Discussion

Non-synonymous variants of SLC10A1, including c.668T>C(p.Ile223Thr), c.800C>T(p.Ser267Phe), c.836T>C(p.Ile279Thr), c.940A>G(p.Lys314Glu), and c.755G>A(p.Arg252His), have been reported previously. The p.Ser267Phe variant exhibits a loss of function for bile salt uptake, but it keeps a fully normal transport function for the non-BA substrate.

In this study, a biochemical assessment of recruited individuals was done to enrich the characterization of the SLC10A1S267F mutation, which is more comprehensive compared with previous studies. Bioinformatics and clinical evidence in this study suggested that this variant was none-pathogenic, which is consistent with previously reported cases of NTCP deficiency. However, the neutrophil count in 2 children (individuals 3 and 8) was lower than normal, while their lymphocyte count was above normal. Two children with a homozygous SLC10A1S267F mutation also presented decreased globulin, suggesting that researchers and doctors should focus on the immunologic function of homozygotes and their susceptibility to pathogens.

Bile and bile salts are vital for endogenous compounds and metabolic products, such as the digestion of dietary fat and fat-soluble vitamins (27-30). According to previous studies, fat-soluble compounds (such as vitamin D) deviate from the normal range in SLC10A1S267F homozygotes. However, our patients displayed slight changes in vitamin D and serum lipid levels. Thus, the specific molecular role of SLC10A1S267F in fat solubilization and absorption needs to be further elucidated.

Some studies have mentioned that patients with SLC10A1S267F display worse clinical phenotype (such as liver dysfunction), always accompanied by other damage variants. Also, our previous results indicated that Slc10a1S267F homozygous mice displayed normal levels of serum BA and did not cause hypercholanemia (26). Thus, to our knowledge, there are probably other variants in SLC10A1S267F homozygotes; in this regard, 5 candidate SNPs are described below.

ABCB11 and ABCC2 are genes coding for the bile salt export pump (BSEP) and multidrug resistance-associated protein 2 (MRP2), respectively. BSEP and MRP2 are bile salt transporters located at the canalicular membrane of hepatocytes (6). BSEP mainly participates in hepatic BA and lipid homeostasis. In most cases, mutations in BSEP are closely related to inherited progressive cholestatic liver disorders in infancy, known as progressive familial intrahepatic cholestasis (PFIC) (31, 32). Missense variant ABCB11V444A is related not only to PFIC but also to a risk factor of intrahepatic cholestasis in pregnancy or drug-induced liver injury (33-35). Therefore, we speculate that 444Ala constitutes a susceptibility factor for the development of hypercholanemia under certain conditions, such as a combination mutation of SLC10A1S267F. MRP2 can transport divalent bile salts. The ABCC2Y39F mutation has found in SLC10A1S267F homozygotes. Some mutations of single amino acid in ABCC2 could confer MRP2 new transport capacity for bile salts, indicating that we should pay attention to this SNP and its underlying effect on serum BA levels.

Mutations of APOBA618V and PSCK9V474I were reported as related factors of familial hypercholesterolemia in the ClinVar database. Apolipoprotein B (APOB) is a major protein constituent of chylomicrons (apo B-48), LDL (apo B-100), and VLDL (apo B-100) (36). The APOBA618V mutation is associated with higher LDL-C levels (odds ratio [OR] = 1.41) (37), meaning that cholesterol and BA levels in SLC10A1S267F homozygotes could be influenced to some extent. The coexistence of APOBA618V and APOBA618V should be studied. The main function of PCSK9 is the degradation of the LDL receptor. In some cases, PCSK9 increases apoB100 production, reduces APOB degradation to increase circulating LDL levels and contributes to hypercholesterolemia and atherosclerosis (38). However, PSCK9V474I seems to be a protection variant associated with lower LDL and cholesterol (39). This might be the reason why the serum lipids in SLC10A1S267F homozygotes were normal in this study.

VPS33B is a class C vacuolar protein sorting protein (40). Approximately 75% of arthrogryposis, renal dysfunction, and cholestasis syndrome (ARC) patients harbor germline mutations in VPS33B. ARC is characterized by congenital joint contractures, renal tubular acidosis, and neonatal jaundice, leading to death in infancy. However, in milder cases, patients survive into adulthood and increase serum BA values but normal bilirubin levels (41). The VPS33BG487S mutation in SLC10A1S267F homozygotes is predicted as a non-specific risk factor in ARC. However, the relationship between VPS33BG487S and SLC10A1S267F mutation still needs to be studied.

From the above discussion, we suggest that although the SLC10A1S267F mutation is the most important reason for high serum BA levels, we detected some mutations in homozygotes that might influence BA metabolism. Further studies on these variants could provide a deeper understanding of the etiology of hypercholesterolemia.

Although we provided a comprehensive clinical, biochemical explanation of the SLC10A1S267F mutation and discussed how to explore the etiology of the SLC10A1S267F mutation-related hypercholesterolemia, there are some limitations. Due to experimental conditions and equipment limitations, we only used medical history data to exclude cholestasis/ malabsorptive diseases or long-term medication history (such as anti-hyperlipidemic agents). Also, we did not evaluate the insulin resistance status of the subjects.

Apart from the function of transporting molecules, NTCP has other functions. Convincing human genetic evidence supports the role of NTCP as a cellular receptor for hepatic B virus (HBV). Recently, several epidemiology studies have reported SLC10A1S267F homozygotes with lower HBV entry and hepatocyte infection compared to wild-type individuals. SLC10A1S267F homozygotes may escape from HBV infection and present decreased cholesterol levels as a consequence of impaired BA uptake (42). According to the marked higher allele frequency of p.S267F in Asia, this SNP might be a positive selection of Asian for its advantage in conferring resistance to chronic HBV. We aim to design more in vivo and in vitro experiments to study and analyze these variants of SLC10A1 in the next phase of our investigation. We have already constructed Slc10a1S267F mutation mice and prepared for further study to determine whether there are deleterious effects on the population and how 267Phe influences HBV infection. To understand the etiology of hypercholanemia, a mouse model of a Slc10a1S267F mutation combined with other related variants is also needed. The effect of p.Ser267Phe on other BA transport proteins should also be explored.

4.1. Conclusions

This study presents an in-depth characterization of NTCP p.Ser267Phe mutations in 5 homozygotes by clinical examination, genetic analysis, and functional prediction. The results supported that the p.Ser267Phe mutation has no deleterious clinical effects in both children and adults.

Acknowledgements

References

  • 1.

    Russell DW. Fifty years of advances in bile acid synthesis and metabolism. J Lipid Res. 2009;50 Suppl:S120-5. [PubMed ID: 18815433]. [PubMed Central ID: PMC2674696]. https://doi.org/10.1194/jlr.R800026-JLR200.

  • 2.

    Hofmann AF. Bile acids: trying to understand their chemistry and biology with the hope of helping patients. Hepatology. 2009;49(5):1403-18. [PubMed ID: 19296471]. https://doi.org/10.1002/hep.22789.

  • 3.

    Hofmann AF, Hagey LR. Bile acids: chemistry, pathochemistry, biology, pathobiology, and therapeutics. Cell Mol Life Sci. 2008;65(16):2461-83. [PubMed ID: 18488143]. https://doi.org/10.1007/s00018-008-7568-6.

  • 4.

    Meier PJ, Stieger B. Bile salt transporters. Annu Rev Physiol. 2002;64:635-61. [PubMed ID: 11826283]. https://doi.org/10.1146/annurev.physiol.64.082201.100300.

  • 5.

    Suzuki H, Sugiyama Y. Transport of drugs across the hepatic sinusoidal membrane: sinusoidal drug influx and efflux in the liver. Semin Liver Dis. 2000;20(3):251-63. [PubMed ID: 11076394]. https://doi.org/10.1055/s-2000-8408.

  • 6.

    Trauner M, Boyer JL. Bile salt transporters: molecular characterization, function, and regulation. Physiol Rev. 2003;83(2):633-71. [PubMed ID: 12663868]. https://doi.org/10.1152/physrev.00027.2002.

  • 7.

    Ho RH, Tirona RG, Leake BF, Glaeser H, Lee W, Lemke CJ, et al. Drug and bile acid transporters in rosuvastatin hepatic uptake: function, expression, and pharmacogenetics. Gastroenterology. 2006;130(6):1793-806. [PubMed ID: 16697742]. https://doi.org/10.1053/j.gastro.2006.02.034.

  • 8.

    Ho RH, Leake BF, Roberts RL, Lee W, Kim RB. Ethnicity-dependent polymorphism in Na+-taurocholate cotransporting polypeptide (SLC10A1) reveals a domain critical for bile acid substrate recognition. J Biol Chem. 2004;279(8):7213-22. [PubMed ID: 14660639]. https://doi.org/10.1074/jbc.M305782200.

  • 9.

    Claro da Silva T, Polli JE, Swaan PW. The solute carrier family 10 (SLC10): beyond bile acid transport. Mol Aspects Med. 2013;34(2-3):252-69. [PubMed ID: 23506869]. [PubMed Central ID: PMC3602841]. https://doi.org/10.1016/j.mam.2012.07.004.

  • 10.

    Vaz FM, Paulusma CC, Huidekoper H, de Ru M, Lim C, Koster J, et al. Sodium taurocholate cotransporting polypeptide (SLC10A1) deficiency: conjugated hypercholanemia without a clear clinical phenotype. Hepatology. 2015;61(1):260-7. [PubMed ID: 24867799]. https://doi.org/10.1002/hep.27240.

  • 11.

    Yan H, Peng B, Liu Y, Xu G, He W, Ren B, et al. Viral entry of hepatitis B and D viruses and bile salts transportation share common molecular determinants on sodium taurocholate cotransporting polypeptide. J Virol. 2014;88(6):3273-84. [PubMed ID: 24390325]. [PubMed Central ID: PMC3957944]. https://doi.org/10.1128/JVI.03478-13.

  • 12.

    An P, Zeng Z, Winkler CA. The Loss-of-Function S267F Variant in HBV Receptor NTCP Reduces Human Risk for HBV Infection and Disease Progression. J Infect Dis. 2018;218(9):1404-10. [PubMed ID: 29905807]. [PubMed Central ID: PMC6151084]. https://doi.org/10.1093/infdis/jiy355.

  • 13.

    Deng M, Mao M, Guo L, Chen FP, Wen WR, Song YZ. Clinical and molecular study of a pediatric patient with sodium taurocholate cotransporting polypeptide deficiency. Exp Ther Med. 2016;12(5):3294-300. [PubMed ID: 27882152]. [PubMed Central ID: PMC5103782]. https://doi.org/10.3892/etm.2016.3752.

  • 14.

    Liu R, Chen C, Xia X, Liao Q, Wang Q, Newcombe PJ, et al. Homozygous p.Ser267Phe in SLC10A1 is associated with a new type of hypercholanemia and implications for personalized medicine. Sci Rep. 2017;7(1):9214. [PubMed ID: 28835676]. [PubMed Central ID: PMC5569087]. https://doi.org/10.1038/s41598-017-07012-2.

  • 15.

    Qiu JW, Deng M, Cheng Y, Atif RM, Lin WX, Guo L, et al. Sodium taurocholate cotransporting polypeptide (NTCP) deficiency: Identification of a novel SLC10A1 mutation in two unrelated infants presenting with neonatal indirect hyperbilirubinemia and remarkable hypercholanemia. Oncotarget. 2017;8(63):106598-607. [PubMed ID: 29290974]. [PubMed Central ID: PMC5739759]. https://doi.org/10.18632/oncotarget.22503.

  • 16.

    Song YZ, Deng M. [Sodium taurocholate cotransporting polypeptide deficiency manifesting as cholestatic jaundice in early infancy: a complicated case study]. Zhongguo Dang Dai Er Ke Za Zhi. 2017;19(3):350-4. Chinese. [PubMed ID: 28302211]. [PubMed Central ID: PMC7390148].

  • 17.

    Van Herpe F, Waterham HR, Adams CJ, Mannens M, Bikker H, Vaz FM, et al. NTCP deficiency and persistently raised bile salts: an adult case. J Inherit Metab Dis. 2017;40(3):313-5. [PubMed ID: 28283843]. https://doi.org/10.1007/s10545-017-0031-9.

  • 18.

    Lin H, Qiu JW, Rauf YM, Lin GZ, Liu R, Deng LJ, et al. Sodium Taurocholate Cotransporting Polypeptide (NTCP) Deficiency Hidden Behind Citrin Deficiency in Early Infancy: A Report of Three Cases. Front Genet. 2019;10:1108. [PubMed ID: 31788003]. [PubMed Central ID: PMC6856633]. https://doi.org/10.3389/fgene.2019.01108.

  • 19.

    Yan YY, Wang MX, Gong JY, Liu LL, Setchell KDR, Xie XB, et al. Abnormal Bilirubin Metabolism in Patients With Sodium Taurocholate Cotransporting Polypeptide Deficiency. J Pediatr Gastroenterol Nutr. 2020;71(5):e138-41. [PubMed ID: 33093374]. https://doi.org/10.1097/MPG.0000000000002862.

  • 20.

    Dai M, Thompson RC, Maher C, Contreras-Galindo R, Kaplan MH, Markovitz DM, et al. NGSQC: cross-platform quality analysis pipeline for deep sequencing data. BMC Genomics. 2010;11 Suppl 4. S7. [PubMed ID: 21143816]. [PubMed Central ID: PMC3005923]. https://doi.org/10.1186/1471-2164-11-S4-S7.

  • 21.

    Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010;26(5):589-95. [PubMed ID: 20080505]. [PubMed Central ID: PMC2828108]. https://doi.org/10.1093/bioinformatics/btp698.

  • 22.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078-9. [PubMed ID: 19505943]. [PubMed Central ID: PMC2723002]. https://doi.org/10.1093/bioinformatics/btp352.

  • 23.

    McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297-303. [PubMed ID: 20644199]. [PubMed Central ID: PMC2928508]. https://doi.org/10.1101/gr.107524.110.

  • 24.

    Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16). e164. [PubMed ID: 20601685]. [PubMed Central ID: PMC2938201]. https://doi.org/10.1093/nar/gkq603.

  • 25.

    Li H, Deng M, Guo L, Qiu JW, Lin GZ, Long XL, et al. Clinical and molecular characterization of four patients with NTCP deficiency from two unrelated families harboring the novel SLC10A1 variant c.595A>C (p.Ser199Arg). Mol Med Rep. 2019;20(6):4915-24. [PubMed ID: 31661128]. [PubMed Central ID: PMC6854589]. https://doi.org/10.3892/mmr.2019.10763.

  • 26.

    Pan Q, Luo G, Qu J, Chen S, Zhang X, Zhao N, et al. A homozygous R148W mutation in Semaphorin 7A causes progressive familial intrahepatic cholestasis. EMBO Mol Med. 2021;13(11). e14563. [PubMed ID: 34585848]. [PubMed Central ID: PMC8573601]. https://doi.org/10.15252/emmm.202114563.

  • 27.

    Boyer JL. Bile formation and secretion. Compr Physiol. 2013;3(3):1035-78. [PubMed ID: 23897680]. [PubMed Central ID: PMC4091928]. https://doi.org/10.1002/cphy.c120027.

  • 28.

    Nathanson MH, Boyer JL. Mechanisms and regulation of bile secretion. Hepatology. 1991;14(3):551-66. [PubMed ID: 1874500].

  • 29.

    Shneider BL, Fox VL, Schwarz KB, Watson CL, Ananthanarayanan M, Thevananther S, et al. Hepatic basolateral sodium-dependent-bile acid transporter expression in two unusual cases of hypercholanemia and in extrahepatic biliary atresia. Hepatology. 1997;25(5):1176-83. [PubMed ID: 9141436]. https://doi.org/10.1002/hep.510250521.

  • 30.

    Stieger B. Recent insights into the function and regulation of the bile salt export pump (ABCB11). Curr Opin Lipidol. 2009;20(3):176-81. [PubMed ID: 19684528]. https://doi.org/10.1097/MOL.0b013e32832b677c.

  • 31.

    Jansen PL, Muller M. The molecular genetics of familial intrahepatic cholestasis. Gut. 2000;47(1):1-5. [PubMed ID: 10861251]. [PubMed Central ID: PMC1727973]. https://doi.org/10.1136/gut.47.1.1.

  • 32.

    Jansen PL, Muller M, Sturm E. Genes and cholestasis. Hepatology. 2001;34(6):1067-74. [PubMed ID: 11731993]. https://doi.org/10.1053/jhep.2001.29625.

  • 33.

    Ali I, Khalid S, Stieger B, Brouwer KLR. Effect of a Common Genetic Variant (p.V444A) in the Bile Salt Export Pump on the Inhibition of Bile Acid Transport by Cholestatic Medications. Mol Pharm. 2019;16(3):1406-11. [PubMed ID: 30608704]. [PubMed Central ID: PMC6467649]. https://doi.org/10.1021/acs.molpharmaceut.8b01124.

  • 34.

    Dixon PH, van Mil SW, Chambers J, Strautnieks S, Thompson RJ, Lammert F, et al. Contribution of variant alleles of ABCB11 to susceptibility to intrahepatic cholestasis of pregnancy. Gut. 2009;58(4):537-44. [PubMed ID: 18987030]. https://doi.org/10.1136/gut.2008.159541.

  • 35.

    Lang C, Meier Y, Stieger B, Beuers U, Lang T, Kerb R, et al. Mutations and polymorphisms in the bile salt export pump and the multidrug resistance protein 3 associated with drug-induced liver injury. Pharmacogenet Genomics. 2007;17(1):47-60. [PubMed ID: 17264802]. https://doi.org/10.1097/01.fpc.0000230418.28091.76.

  • 36.

    Sniderman AD, Thanassoulis G, Glavinovic T, Navar AM, Pencina M, Catapano A, et al. Apolipoprotein B Particles and Cardiovascular Disease: A Narrative Review. JAMA Cardiol. 2019;4(12):1287-95. [PubMed ID: 31642874]. [PubMed Central ID: PMC7369156]. https://doi.org/10.1001/jamacardio.2019.3780.

  • 37.

    Hayat M, Kerr R, Bentley AR, Rotimi CN, Raal FJ, Ramsay M. Genetic associations between serum low LDL-cholesterol levels and variants in LDLR, APOB, PCSK9 and LDLRAP1 in African populations. PLoS One. 2020;15(2). e0229098. [PubMed ID: 32084179]. [PubMed Central ID: PMC7034850]. https://doi.org/10.1371/journal.pone.0229098.

  • 38.

    Reiss AB, Shah N, Muhieddine D, Zhen J, Yudkevich J, Kasselman LJ, et al. PCSK9 in cholesterol metabolism: from bench to bedside. Clin Sci (Lond). 2018;132(11):1135-53. [PubMed ID: 29895529]. https://doi.org/10.1042/CS20180190.

  • 39.

    Chuan J, Qian Z, Zhang Y, Tong R, Peng M. The association of the PCSK9 rs562556 polymorphism with serum lipids level: a meta-analysis. Lipids Health Dis. 2019;18(1):105. [PubMed ID: 31036026]. [PubMed Central ID: PMC6489332]. https://doi.org/10.1186/s12944-019-1036-1.

  • 40.

    Hanley J, Dhar DK, Mazzacuva F, Fiadeiro R, Burden JJ, Lyne AM, et al. Vps33b is crucial for structural and functional hepatocyte polarity. J Hepatol. 2017;66(5):1001-11. [PubMed ID: 28082148]. [PubMed Central ID: PMC5387182]. https://doi.org/10.1016/j.jhep.2017.01.001.

  • 41.

    Smith H, Galmes R, Gogolina E, Straatman-Iwanowska A, Reay K, Banushi B, et al. Associations among genotype, clinical phenotype, and intracellular localization of trafficking proteins in ARC syndrome. Hum Mutat. 2012;33(12):1656-64. [PubMed ID: 22753090]. [PubMed Central ID: PMC3746110]. https://doi.org/10.1002/humu.22155.

  • 42.

    Cheng X, Wang Y, Tian J, Zhou L, Chen X, Guo H, et al. SLC10A1 S267F variant influences susceptibility to HBV infection and reduces cholesterol level by impairing bile acid uptake. J Viral Hepat. 2019;26(10):1178-85. [PubMed ID: 31177598]. https://doi.org/10.1111/jvh.13157.