Serum vitamin D level is inversely associated with liver fibrosis in post Kasai's portoenterostomy biliary atresia patients living with native liver.

OBJECTIVE: This study aims to investigate the association of serum vitamin D (VD) levels with the severity of liver fibrosis (LF) in post Kasai's portoenterostomy biliary atresia (PKBA) patients living with their native liver.
METHODS: In this cross-sectional study, carried out in a tertiary Children's Hospital in Taipei, Taiwan, PKBA patients living with their native liver were enrolled. Liver biochemistry data, serum 25-hydroxyvitamin D (25-OHVD), acoustic radiation force impulse (ARFI), and scores of Pediatric Quality of Life questionnaire (PedsQL) were collected.
RESULTS: All the enrolled 33 PKBA patients (36.4% males), aged 1-23 years, possessed 25-OHVD less than 30ng/ml. An inverse correlation was detected between serum 25-OHVD and ARFI (r2 = 0.175; p = 0.024). We selected a cutoff value of 23ng/mL to divide PKBA patients into two groups, as the p-value was the most significant at this point when comparing the median ARFI of two groups (p = 0.003). Ten (30.3%) had 25-OHVD≥23ng/ml (HVD group), whereas 23(69.7%) had 25-OHVD<23ng/ml (LVD group). HVD group had lower ARFI (1.13m/s vs. 1.52m/s, p = 0.003), lower aspartate transaminase (AST) (29U/L vs. 64U/L, p = 0.033), and higher scores of self-reported health-related quality of life in psychosocial functioning (86.7 vs. 77.1, p = 0.047) than LVD group.
CONCLUSION: VD levels are inversely associated with severity of LF in PKBA patients with native liver.

Introduction

Biliary atresia (BA) is the most common cause of liver cirrhosis in children. BA occurs in approximately 1.23 of 10,000 live births in Taiwan and 0.56 of 10,000 live births in European countries [1, 2]. Kasai’s portoenterostomy (KPE) is the first line therapy to reconstruct the obliterated bile ducts and restore bile flow. The introduction of KPE made long-term survival possible. In Japan, the 10-year native liver survival rate is 52.8%, while in Taiwan, the infant stool card screening program for BA started in 2002 allowed earlier KPE and yielded improvement of 5-year-survival rate with native liver (NL) from 56% to 89%[3, 4]. However, more than 60% of patients experienced complications such as ascending cholangitis (36%), portal hypertension (40%), gastrointestinal bleedings (25%) and hepatocellular carcinoma [5, 6]. Only 30% of the patients with BA living with NL had no clinical complications of chronic liver disease (CLD), and only 2–11% had absence of surgical complications and normal laboratory indices [6, 7]. Liver transplantation (LT) is considered in those who fail to clear jaundice by KPE or with progressive liver cirrhosis. With growing experience and advanced techniques of pediatric LT, the 5-year recipient survival rate is 82–98%[3].

Previously, it was believed that vitamin D (VD) deficiency is mainly caused by bile acid-related lipid malabsorption in cholestatic liver diseases, however, recent studies showed that both cholestatic and non-cholestatic liver diseases share mutual characteristic of VD deficiency [8]. In addition, VD deficiency has been linked to various adverse outcomes of chronic liver diseases in adults. In pediatric studies, serum VD levels are inversely associated with non-alcoholic steatohepatitis (NASH) and fibrosis in children with non-alcoholic fatty liver disease (NAFLD)[9]. Nevertheless, little research has been done to explore the relationship between VD status and the severity of liver fibrosis (LF) in post Kasai’s portoenterostomy BA (PKBA) patients. The aim of the present study was to investigate the association of serum VD levels with the severity of LF in PKBA patients living with their NL.

Materials and methods

Subjects

PKBA patients were consecutively recruited in this cross-sectional study between October, 2016 to March, 2018. The setting was a tertiary children’s hospital in Taipei, Taiwan. Grown-up patients who maintained regular hospital visits in MacKay Children’s Hospital were also enrolled in the study. As pure VD oral supplements were not introduced into Taiwan until 2014, the PKBA patients in our institute were encouraged to have more sun exposure, but most of them did not receive routine VD supplements at the time of study enrollment. Two patients taking vitamin and/or mineral supplements and/or medications (polyvitamin, calcitriol, calcium carbonate) known to influence serum 25-hydroxyvitamin D (25-OHVD) status were excluded from the study. Patients who had received LT were also excluded from the study. None of the enrolled patients used VD supplements upon recruitment.

The study was approved by the Institutional Review Board of the MacKay Memorial Hospital and the approved number was 17MMHIS051; the informed consents were obtained from the parents or legal guardians of the patients or the patients themselves who aged over 18 years.

Demographics and laboratory assessment

Age, gender, month-old when patients received KPE (before or after 2-month-old), and hospitalization frequency due to acute cholangitis were recorded. Biochemical parameters including aspartate transaminase (AST), alanine transaminase (ALT), gamma glutamyl transferase (rGT), direct bilirubin, total bilirubin, and 25-OHVD were obtained. The levels of 25-OHVD were measured with chemiluminescent method by Diasorin. Currently, the normal values of serum 25-OHVD are still in debate, while generally the suggested values range from 10 to 40ng/mL [10-12]. In this study, VD deficiency is defined as serum 25-OHVD levels lower than 20ng/mL and VD insufficiency is defined as serum levels between 20 and 30ng/mL.

Acoustic radiation force impulse (ARFI)

All the participants were studied with a 6-MHz probe on an Acuson S3000 US machine with a Virtual Touch Tissue Quantification (VTQ) mode for ARFI elastography (Siemens Healthcare, Erlangen, Germany). ARFI measures the shear wave speed and was regarded as the marker of elasticity in the examined tissue (expressed in meter/second). ARFI values and blood tests were obtained no more than three months apart. The ARFI values were measured independently by four pediatric gastroenterologists with over 15 years of experience doing abdominal ultrasounds and at least 1 year of experience using the ARFI modality. The operator selected a rectangular region of the right liver through the intercostal space, avoiding bile ducts and blood vessels, with the size of 1cm x 0.6cm, 4 cm to 6cm from the cutaneous surface. Patients were examined in the supine position, with their right arm in maximum abduction. Patients were asked to breathe smoothly or slightly holding their breath during the measurements of ARFI. Three measurements were made and the mean value was calculated.

Pediatric quality of life questionnaire (PedsQL)

The PedsQL Measurement Model is a reliable and valid modular approach to measuring health-related quality of life (HRQoL) in children and adolescents with acute and chronic health conditions. It is brief (contains 23 items), practical (takes less than 4 minutes to complete), developmentally appropriate (tailored for different age groups), and is translated into multiple languages. The 23 items are categorized into 4 multidimensional scales, as physical functioning (8 items), emotional functioning (5 items), social functioning (5 items), and school functioning (5 items). The 3 summary scores are total scale score (23 items), physical health summary score (8 items), and psychosocial health summary score (15 items). The interpretation of the questionnaire was made after transformation and calculation according to the PedsQL scoring procedure. The higher the score, indicates better HRQoL. Our study team was authorized to use PedsQL questionnaire to assess the quality of life of the participants. The questionnaire was distributed to the participants and their parents either by postal mails or in person during their hospital visits.

Statistical analysis

Descriptive statistics were expressed as median (interquartile range, IQR). Discrete variables were presented as frequency and percentage. Differences in continuous variables between the 2 groups were analyzed by Mann Whitney U test. Differences in categorical variables between the 2 groups were compared by Fisher’s exact test. The association of LF with serum 25-OHVD was detected using a linear regression model.

Statistical analysis was performed with the Statistical Package for the Social Sciences (SPSS) version 21.0 (SPSS, Inc., USA). All statistical tests were two-tailed, and a p-value < 0.05 was considered statistically significant.

Results

Grouping

A total of 33 patients (36.4% males) aged 1–23 years were enrolled into our study. All the patients’ 25-OHVD were less than 30ng/ml. According to the scatter diagram, there was an inverse association between serum 25-OHVD and ARFI values (r2 = 0.175; p = 0.024)(Fig 1). We set multiple cutoff values of 25-OHVD dividing the patients into two groups (higher 25-OHVD and lower 25-OHVD groups), compared the medians of ARFI of the two groups with Mann-Whitney U test and calculated the p-values (Table 1). The p-value was the most significant when the cutoff value of 25-OHVD was 23ng/mL (p = 0.003). Of the 33 patients, 10 (30.3%) had 25-OHVD≥ 23ng/ml, whereas 23(69.7%) had 25-OHVD<23ng/ml. We name the group with 25-OHVD≥ 23ng/ml as higher VD level group (HVD group), and the group with 25-OHVD<23ng/ml as lower VD group (LVD group).

Fig 1

Serum 25-OHVD had a negative correlation with ARFI (r2 = 0.175; p = 0.024).

Abbreviations: 25-OHVD: 25-hydroxyvitamin D; ARFI: Acoustic radiation force impulse.

Table 1

Different cutoff values of 25-OHVD and corresponding p-values when comparing the medians of ARFI of the two groups.

25-OHVD ng/mL as cutoff value	p-value	25-OHVD ng/mL as cutoff value	p-value
10	0.080	18	0.180
11	0.061	19	0.062
12	0.061	20	0.033
13	0.330	21	0.040
14	0.100	22	0.011
15	0.138	23	0.003
16	0.140	24	0.009
17	0.180	25	0.078

Abbreviations: 25-OHVD: 25-hydroxyvitamin D; ARFI: Acoustic radiation force impulse

Age and gender

There were 12 males and 21 females among the 33 PKBA patients. LVD group consisted of 7 males and 16 females, and HVD group consisted of 5 males and 5 females. There was no age difference between LVD and HVD groups (LVD: 12.0 (5.0–16.0), HVD: 12.0 (5.0–19.0), p = 0.814). There was no correlation between age and 25-OHVD in our patients (r2 = 0.06, p = 0.17); moreover, neither of the gender showed decline of 25-OHVD with increase of age (female: r2 = 0.093, p = 0.167; male: r2 = 0.014, p = 0.725). (Fig 2).

Fig 2

Relation between age in years and serum 25-OHVD (ng/mL) concentrations stratified by gender.

Abbreviations: 25-OHVD: 25-hydroxyvitamin D.

Serum biochemical parameters

The median AST value is higher in LVD group than HVD group (64 vs. 29, p = 0.033). The medians of other measurements including ALT, total bilirubin, direct bilirubin and rGT were also higher in LVD group comparing to those in HVD group, although there was no statistical significances (Table 2).

Table 2

Demographic and liver function tests of PKBA patients.

Parameter	All patients (n = 33)	25-OHVD < 23 (n = 23)	25-OHVD ≥ 23 (n = 10)	p-value
Age (years)	12.0 (5.5–18.5)	12.0 (5.0–16.0)	12.0 (5.0–19.0)	0.814
Male, No. (%)	12 (36.4)	7 (30.4)	5 (50)	0.433
Weight for age				0.054
5–15%	2	2 (8.7%)	0 (0)
15–50%	7	2 (8.7%)	5 (50%)
50–85%	15	11 (47.8%)	4 (40%)
85–97%	9	8 (34.8%)	1 (10%)
AST (U/L)	48 (28.5–94.5)	64 (34.5–132.75)	29 (21–46)	0.033
ALT (U/L)	78 (30–123.5)	98 (39.25–202)	32 (17–92)	0.063
Total bilirubin (mg/dL)	0.7 (0–1.45)	1.4 (1.1–1.9)	0.7 (0.7–0.7)	0.259
Direct bilirubin (mg/dL)	0.3 (0.2–0.9)	0.3 (0.2–0.9)	0.1 (0.1–0.1)	0.166
r-GT (IU/L)	103 (50–207)	121.5 (71.75–235)	64 (31–124.5)	0.139

Data are expressed as median (interquartile range) or number (percentage).

Abbreviations: 25-OHVD: 25-hydroxyvitamin D; ALT: alanine transaminase; AST: aspartate transaminase; PKBA: post Kasai’s portoenterostomy biliary atresia; rGT: gamma glutamyl transferase

ARFI values and PedsQL scores

25-OHVD had a negative correlation with ARFI (r2 = 0.175; p = 0.024) in Fig 1. Patients in LVD group had higher ARFI values (1.52 vs. 1.13, p = 0.003) than those in HVD group, indicating patients with lower 25-OHVD had more severe LF (Fig 3). The psychosocial health summary score is significantly lower in LVD group compared to HVD group (77.1 vs. 86.7, p = 0.047) (Table 3).

Fig 3

Patients in LVD group had higher ARFI values than those in HVD group (1.52m/s vs. 1.13m/s, p = 0.003), indicating patients with lower serum 25-OHVD levels had more severe liver fibrosis.

Abbreviations: LVD: lower vitamin D level; ARFI: Acoustic radiation force impulse; HVD: higher vitamin D level; 25-OHVD: 25-hydroxyvitamin D.

Table 3

Serum 25-OHVD, ARFI values, and PedsQL scores of PKBA patients.

Parameter	25-OHVD < 23 (n = 23)	25-OHVD ≥ 23 (n = 10)	p-value
25-OHVD (ng/mL)	16.7 (12.2–20.0)	27.2 (24.7–28.2)	<0.001
ARFI (m/s)	1.52 (1.23–1.67)	1.13 (0.93–1.27)	0.003
Physical score	92.2 (67.2–97.7)	90.6 (75.0–96.9)	0.877
Psychosocial score	77.1 (77.3–87.1)	86.7 (80.8–93.8)	0.047
Total score	80.6 (73.4–90.8)	87.5 (81–91.3)	0.207
Proxy physical score	90.6 (71.9–96.9)	85.9 (70.3–95.3)	0.739
Proxy psychosocial score	81.7 (68.3–91.7)	71.7 (62.9–81.3)	0.321
Proxy total score	76.0 (67.4–87.0)	83.7 (71.5–91.3)	0.482

Data are expressed as median (interquartile range).

Abbreviations: 25-OHVD: 25-hydroxyvitamin D; ARFI: Acoustic radiation force impulse; PedsQL: Pediatric Quality of Life questionnaire; PKBA: post Kasai’s portoenterostomy biliary atresia

Discussion

In this single center study, we showed that low VD status is extremely common in PKBA patients with NL in different ages. All our patients had 25-OHVD below 30ng/mL, with insufficiency in 17 patients, and deficiency in 16 patients. It is not surprising to see such universe inadequacy of serum 25-OHVD in our patient category, in which most of them have CLD. Our study observed an inverse association between serum 25-OHVD and severity of LF indicated by ARFI values. After respective analyses of multiple cutoff values dividing our patients into two groups of 25-OHVD statuses, we found that the two groups exhibited the most significant difference when the cutoff value was 23ng/mL. The median ARFI value in LVD group was 1.52m/s, which corresponded to F2 (moderate fibrosis) in METAVIR scoring system; and the median ARFI value in HVD group was 1.13m/s, which corresponded to F0 (no fibrosis).

Additionally, we found that LVD group displayed a higher AST level than HVD group (64IU/L vs. 29IU/L, p = 0.033), which probably suggested a more severe hepatic inflammation in LVD group. Self-reported HRQoL summed a lower score of psychosocial functioning in LVD group; it could be postulated that lower VD status reflects poorer liver condition and is prone to relatively impaired psychosocial function, or that lower VD status is connected to depressive disorder to some degree. However, the other HRQoL dimensions did not show concordant results as psychosocial functioning between the two groups.

VD deficiency is a global health problem that had received increasing attention recently. Normal values of serum 25-OHVD are still in debate. According to the Institute of Medicine of the National Academies in the United States, VD concentration of 20ng/mL is largely sufficient [10]. However, the Endocrine Society recommends higher thresholds, above 30ng/mL, as sufficient. Because evidence show that VD levels above 20ng/ml prevent rickets and osteomalacia; while level above 30ng/ml maximizes 25-OHVD’s effect on calcium, bone and muscle metabolism, and may have additional benefits in reducing cancers, autoimmune diseases, type 2 diabetes, cardiovascular disease, and infectious diseases. In this study, we define VD deficiency as serum 25-OHVD levels lower than 20ng/mL and insufficiency as levels between 20 and 30ng/mL [11, 12].

According to recent survey, VD deficiency accounts for 8–12% in healthy white infants and toddlers [13, 14]. In southern China, VD deficiency is 10.8% and insufficiency is 39% among children [15]. In northern Taiwan, among healthy 3–18 year-olds, 51% have VD deficiency, and 40% have VD insufficiency [16]. Beyond regulation of calcium and bone homeostasis, VD has an emerging role in proliferation, differentiation, inflammation, and immunomodulation. VD2 is synthesized from food, whereas VD3 is mainly synthesized in the epidermis after sun exposure [17]. Bile acids are requisite for gastrointestinal absorption of dietary sources of VD. VD2 and D3 are bound to VD binding protein (VDBP) in the blood circulation, undergo hydroxylation in the liver forming 25-OHVD, then hydroxylate mainly in the kidneys into 1,25(OH)2D [18]. 1,25(OH)2D is also produced by immune cells in a paracrine fashion, which participates in immunomodulation [19]. 1,25(OH)2D, the biologic active form of VD, exerts its diverse functions by binding with VD receptors (VDR), which are expressed in many different tissues [14]. It has been reported that low VD status is a common feature in different types of liver diseases [14]. Previously, it was believed that VD deficiency is attributed to bile acid-related lipid malabsorption in cholestatic liver diseases; however, recent studies show that both cholestatic and non-cholestatic liver diseases share mutual characteristic of VD deficiency [8]. In adult studies, serum 25-OHVD levels correlate inversely with liver disease severity in cirrhosis [8, 14]. Possible mechanisms include malnutrition, decreased sun exposure, reduced synthesis of VDBP, and impaired bile acid excretion. Moreover, advanced studies discover a causal association between VD deficiency and LF. VDR is abundant in hepatic stellate cells (HSCs); HSCs are the major effective cells in the wound healing process when encounters liver injury, engaging in deposition of extracellular matrix proteins. In vitro study show that activation of VDR in HSCs strongly antagonize the TGFβ signaling pathway, which is the most potent pro-fibrogenic pathway in liver [20]. In addition, preclinical studies have demonstrated that systemic administration of a VDR ligand (calcitrol or calcipotriol) effectively inhibited liver fibrogenesis in standard mouse models of liver injury [21]. Therefore, maintaining sufficient VD levels is postulated to inhibit the process of LF.

Fat-soluble vitamin (FSV) deficiency is extremely common in pre-KPE BA patients. VD deficiency alone accounts for 88–98% of pre-KPE BA patients, whereas nearly half of the patients present with two or more FSV deficiency [13, 22]. Even with successful KPE and oral supplementation of FSV, the prevalence of at least one FSV insufficiency in PKBA infants is 55% throughout six months after KPE; and the prevalence of FSV insufficiency in PKBA infants with persistent cholestasis six months after KPE is 100%, 80%, 50%, and 45%, respectively, for vitamins A, D, E, and K [23]. Bile acid secretion in the liver does not reach normal levels until 6 to 12 months after KPE is probably the main reason [22]. In addition, many studies show unsatisfactory repletion of serum 25-OHVD with enteral VD supplementations in PKBA patients, likely owing to inadequate intraluminal bile acids to solubilize orally administered VD [13, 23–25]. Many experts suggest higher doses of oral VD or parenteral VD supplements to overcome this difficulty. Interestingly, the data from our study revealed that, even years after KPE, and despite clearance of jaundice, VD deficiency remains a generalized issue in PKBA patients. This supports the current understanding that multiple factors contribute to VD deficiency instead of inadequate intraluminal bile acids alone, and highlights the necessity of persistent FSV supplementation in PKBA patients. In our institute, we now treat PKBA patients with 25-OHVD<12ng/mL with oral VD3 4000IU per day, 12ng/mL≤25-OHVD<20ng/mL with 3000IU per day, 20ng/mL≤25-OHVD<30ng/mL with 2000IU per day, and 25-OHVD≥30ng/mL with 600-800IU per day. The oral VD dosage should be tailored according to sequential serum 25-OHVD levels monitored at least once a year.

Notwithstanding the considerable literature on the association of VD status and various liver diseases, this is the first study to demonstrate the association between VD status and LF in PKBA patients living with their NL. However, a series of limitations should be considered. First, this is a cross-sectional study thus cannot provide causation. Second, the sample size was small due to disease rarity, and it did contribute to reduced power of the statistics. Third, the lack of information such as the type of diet, duration of sun exposure, and the lack of serum values of VDBP does not allow us to exclude these possible confounders. Fourth, the most accurate way to evaluate LF severity is liver histology, but we substituted with ARFI, a non-invasive method instead. Fifth, the patients’ age widely ranged from 1 to 23 years old. One study (participants aging 6 to14 year-old) showed that girls had a 4% decline in serum 25-OHVD levels for each increasing year of age [26]. However, in our study, we found no correlation between age and 25-OHVD (r2 = 0.06, p = 0.17); neither of the gender showed decline of 25-OHVD with increase of age (female: r2 = 0.093, p = 0.167; male: r2 = 0.014, p = 0.725).

Hopefully in the near future, we can determine the appropriate dosage of (enteral and parenteral) VD supplementation for PKBA patients, and elucidate a causal relationship between VD status and LF in PKBA patients, by showing decrease in ARFI values following serum increment of 25-OHVD in a larger-scale, longitudinal cohort study. As the majority of PKBA patients with NL possess CLD, it is valuable to improve their quality of life and reduce complications of CLD. If we can conclude that reaching a specific serum 25-OHVD level halts the process of LF in PKBA patients, better quality of life may be achieved, and LT may be avoided or postponed. Whether simply keeping serum 25-OHVD above 23ng/mL is of clinical benefits, or higher levels are even better for these patients, is also an issue prompting further investigations.

Conclusion

This study emphasizes the universal inadequacy of VD in PKBA patients with NL (ages range from 1 to 23 year-old and at jaundice-free phase), and the importance of monitoring their VD status. We detected a significant inverse association between 25-OHVD levels and severity of LF measured by ARFI values. Patients with serum 25-OHVD levels equal or above 23ng/mL had less severe LF than patients with serum 25-OHVD levels below 23ng/mL. Patients with higher 25-OHVD levels also had significantly lower AST levels and better self-reported HRQoL in psychosocial functioning.

Raw data used for this manuscript.

(XLSX)

Click here for additional data file.