SLCO1B1 c.388A > G variant incidence and the severity of hyperbilirubinemia in Indonesian neonates
BMC Pediatrics volume 19, Article number: 212 (2019)
It has been established that genetic factors play a substantial role in the development of neonatal hyperbilirubinemia. The population of Indonesia and other Southeast Asian countries has similar, yet different genetic makeup compared to the rest of Asia. Aside from UGT1A1, variants of SLCO1B1 have also been known to contribute to the severity of neonatal hyperbilirubinemia in Asian populations. Since there has been no report on SLCO1B1 polymorphism in relation with hyperbilirubinemia in Indonesia, this study aims to explore incidence of SLCO1B1*1B polymorphism in Indonesia based on 3 hospitals from different provinces and population, and their association with hyperbilirubinemia severity.
Our study included 88 neonates with mild and moderate-severe hyperbilirubinemia from 3 NICU in hospitals representing homogenous and heterogenous populations: Biak General Hospital Papua, Cipto Mangunkusumo Hospital (Jakarta), and M Yunus Hospital (Bengkulu). We collected samples between November 2016 and September 2017. DNA was obtained from existing samples of the patients from previous studies and were subjected to Polymerase Chain Reaction – Restriction Fragment Length Polymorphism (PCR-RFLP). We analyzed the *1B variant located in exon 5 of SLCO1B1 with TaqI restriction endonuclease. Clinical, demographic, and laboratory data was also collected from medical records and parents’ interviews.
The most dominant variant of SLCO1B1*1B in our population is the homozygous G/G (68.18%), followed by heterozygous A/G (26.14%), and wild type A/A (5.68%). The heterozygous A/G had an Odds Ratio (OR) of 0.73 (95% CI 0.10–5.2) and homozygous G/G with OR of 0.51 (95%CI 0.08–3.27), both were not significant. Genotypic distribution across the different centers were also similar and not significant. The significant risk factors for moderate-severe hyperbilirubinemia were the population the neonate originated from (p = < 0.001) and the delivery location (p = 0.001), while SLCO1B1*1B was not associated with the different severity of hyperbilirubinemia.
SLCO1B1*1B is not associated with higher bilirubin levels among neonates with hyperbilirubinemia in Indonesia. Further study is needed to find other potentially important genetic polymorphisms in the development of severe hyperbilirubinemia in Indonesia.
Neonatal hyperbilirubinemia (jaundice) occurs in about 60% of all term infants . In general, most of the increase in bilirubin levels is physiological and usually will improve by the end of the first week of life. An increase in total serum bilirubin level above 5 mg/dL can cause jaundice that is visible to observers . This condition can be severe if the infant was born prematurely, where the incidence increases to 80% . In these infants, hyperbilirubinemia can develop into severe neonatal jaundice (SNNJ) that can cause bilirubin-induced neurological damage (BIND) and ultimately lead to irreversible neurodevelopmental impairment or even death .
Aside from prematurity, genetic polymorphisms have been reported as an important factor in the development of SNNJ [4,5,6]. The most commonly-studied gene is the uridine diphosphate glucuronosyltransferase 1A1 (UGT1A1) that is closely associated with Crigler-Najjar’s syndrome type I and II, and Gilberts syndrome. These conditions are frequently found in Caucasian and African populations [7, 8]. Another gene that has been gaining attention, one that is mostly found in Asian populations, is OATP2/SLCO1B1 . SLCO1B1 (solute carrier organic anion transporter family member 1b1) encodes the solute carrier organic anion transporter family member 1B1 and is responsible for the absorption of unconjugated and conjugated bilirubin in the liver. This gene, located in chromosome 12p12, consists of 15 exons (1 coding exon and 1 non-coding exon) and 14 introns, as well as containing 2073 nucleotides in the gene-coding area . Neonates carrying this gene polymorphism are associated with reduced bilirubin clearance which induces hyperbilirubinemia . Rotor Syndrome, an autosomal recessive disease associated with elevated conjugated bilirubin, is an example of a disease frequently associated with mutations in SLCO1B1 and SLCO1B3. The most frequent polymorphisms found in this gene are SLCO1B1*1B (388A > G), SLCO1B1*4 (463C > A), and SLCO1B1*5 (521 T > C), and the most commonly associated with hyperbilirubinemia is the SLCO1B1*1B (388A > G) [9, 11,12,13].
Based on previous studies, there is an association between SLCO1B1 with racial variety in regard to the severity of unconjugated neonatal hyperbilirubinemia. In Indian, Japanese, Taiwanese, and Chinese neonates, the incidence of polymorphism is higher in infants with higher bilirubin levels compared to those with lower bilirubin levels [9, 11, 12, 14]. However, in studies of Malaysian and Thai hyperbilirubinemic neonates, this SNP was not frequently found [15, 16]. In Indonesia, polymorphism studies on neonatal hyperbilirubinemia has yet to be conducted extensively. Therefore, this study provides valuable insight into the association of genetic factors especially SLCO1B1 with the severity of neonatal hyperbilirubinemia in Indonesia.
Materials and methods
We conducted a cross-sectional study of patients in neonatal intensive care unit across three hospitals; M. Yunus General Hospital Bengkulu, Biak General Hospital Papua, and Cipto Mangunkusumo Hospital (CMH) Jakarta based on total population sampling. We collected samples and clinical data from November 2016 to September 2017, followed by molecular analysis in 2018. Diagnosis of hyperbilirubinemia was based on clinical jaundice as observed by a neonatologist based on Kramer’s index between day 3–7 since birth before receiving phototherapy, and peak Total Serum Bilirubin (TSB) was then taken prior to initiation of phototherapy. Mild hyperbilirubinemia was considered below 13 mg/dL, and moderate to severe was considered above and equal to 13 mg/dL . Both inborn and outborn, preterm and term infants were included. Patients with hemolytic anemia, cephalhematoma, neonatal sepsis, ABO and rhesus incompatibility, maternal diabetes, and other comorbidities that would affect the total serum bilirubin level were excluded through clinical and laboratory approach. Clinical and sociodemographic variables were obtained through interviews and medical records. Informed consent was obtained from both parents prior to inclusion into the study. This study was approved by the Ethics Committee of Faculty of Medicine, Universitas Indonesia.
Capillary blood samples were obtained from 88 infants and were used to measure their peak total serum bilirubin using the ADVIA Chemistry Total Bilirubin 2 device (diazo method)  taken between day 3 and 7 before phototherapy. A blood sample was also collected from each infant and stored at − 20 °C prior to DNA extraction.
DNA extraction and manipulation
DNA extraction from blood samples was conducted using the QIAmp Blood DNA Mini Kit (QIAGEN, Germany) according to the manufacturer’s protocol as previously described by Amandito et al. . Following extraction, the DNA concentration was then measured by Nano-spectrophotometer (Thermo-Scientific) for use in polymerase chain reaction - restriction fragment length polymorphism (PCR-RFLP) analysis for determining the SNP in the exon 5 of SLCO1B1.
As much as 2uL of DNA was amplified using KOD FX Neo Kit (Toyobo) with each 25uL mix containing 60 ng extracted DNA, 12.5 uL KOD FX Neo Buffer, 5uL dNTPs (2 mM each of dATP, dCTP, dGTP, dTTP), 0.75uL of each primer (final concentration 8 uM), and 0.7 uL KOD FX Neo DNA polymerase (Novagen, Germany). Amplification was carried out with the following parameters: initial denaturation for 2 min at 95 °C, followed by 35 cycles of 1 min denaturation at 94 °C, annealing for 30 s at 57 °C, extension for 30 s at 68 °C, and final polymerisation for 5 min at 68 °C, and hold at 8 °C for 60 min (Biometra, Germany). All PCR products were visualized by agarose gel electrophoresis in 2% (wt/vol) agarose gels with TBE buffer.
To detect the 388A > G SNP (rs2306283), as according to a previous study, we used the primers 388F 5′-ATA ATG GTG CAA ATA AAG GGG-3′ (IDTDNA) for Forward Primer and 388R 5′-ACT ATC TCA GGT GAT GCT CTA-3′ (IDTDNA) for Reverse Primer . For digestion, we used the restriction endonuclease TaqI (NEB) and CutSmart buffer (NEB). The digestion procedure consisted of incubation at 65 °C, followed by enzyme inactivation in 80 °C in a Thermal Cycler PCR MJ Mini (Bio-Rad). The sequenced products that were digested were analyzed by electrophoresis on 3% agarose gels with TBE buffer (Fig. 1).
We used SPSS version 24 for data management and statistical analysis. We conducted descriptive analysis for our study prior to inclusion for analytical study. We conducted bivariate analysis between independent variables and dependent variables using parametric and non-parametric methods, as appropriate. Variables with p value < 0.25 were included in the multivariate analysis using multiple logistic regression. We used two-sided p values in our analysis. Gender, population, exclusive breastfeeding, delivery method, delivery location, sibling requiring phototherapy, and SNPs were categorical data, while gestational age, birth weight, and mother’s age were numerical data.
We enrolled 88 hyperbilirubinemic neonates from three different hospitals; Bengkulu (a west coast city on Sumatra Island), Jakarta (capital of Indonesia located on the island of Java), and Biak (a far-eastern remote town on Papua Island), consisting of 47 in the control group (mild hyperbilirubinemia), and 41 in the case group (moderate-severe hyperbilirubinemia). Clinical data of all the patient groups based on the severity of hyperbilirubinemia are shown in Table 1. The mean bilirubin level of all neonates was 13.96 mg/dL with a SD of 6.57, mean gestational age of 33.23 weeks for control and 35.76 weeks for case group, median birth weight of 1580 g for the control group and 2500 g for the case group, and comparable gender between the two groups. Population, delivery location, gestational age, and birth weight were statistically-significant risk factors for moderate-severe hyperbilirubinemia. On the other hand, gender, exclusive breastfeeding, delivery method, sibling requiring phototherapy, and age of the mother were not significant risk factors, and therefore these factors were not included in the multivariate analysis.
For the multiple logistic regression analysis (Table 2), the most significant factor was population with OR of 4.62 (95% CI 1.79–11.97) and p value of < 0.001, while delivery location, gestational age, and birth weight were no longer significant.
Detection of SLCO1B1*1B polymorphism
Table 3 shows the genotypic frequencies of SLCO1B1*1B variants due to the SNP in exon 5. In the control group, most of the SNPs were homozygote G/G with 34 neonates (72.3%), and wildtype A/A was the least common with only 2 neonates (4.3%). The same is in the case group with 26 neonates (63.4%) homozygote G/G and 3 neonates (7.3%) were wild-type A/A. No statistical significance was achieved when analyzing the two groups based on SLCO1B1 SNP variant. Table 4 shows the analysis between SNPs and bilirubin level which shows no statistical significance.
Table 5 highlights the distribution of SLCO1B1*1B between the three centers which represents three different ethnicities and environments, with CMH, Jakarta being the metropolitan capital city with heterogenous ethnicity, M. Yunus, Bengkulu an isolated west coast city with a homogenous ethnicity, and Biak, a remote east town on Papua, also with a homogenous ethnicity. In all three centers, homozygote variant of SLCO1B*1B was the most dominant (75, 57.7, and 50%, respectively), with the heterozygotic variant being the second-most dominant and the wild-type being the least dominant. There was no significant association between the population and the pattern of SLCO1B1*1B polymorphic variance. Table 6 shows the demographic data across the three different hospitals. Between them, only gender (p = 0.004) and ethnicity (p < 0.001) were significantly different.
We conducted additional analysis from both groups and within each group to find any association between SLCO1B1*1B and bilirubin level and found there was no significant association.
Physiological conditions such as the short neonatal red blood cell’s lifespan, absorption of bilirubin by the liver and a non-optimal conjugation process of bilirubin can cause hyperbilirubinemia [1, 2]. Premature birth carries a higher risk of severe hyperbilirubinemia because of the immaturity of the red blood cells and an impaired hepatic cell function responsible for neonatal the bilirubin metabolism in neonates [2, 20]. The cause for this condition also varies between populations due to suspected contributions of genetic factors towards the pathogenesis of hyperbilirubinemia.
SLCO1B1 is a sinusoidal membrane protein expressed at the basolateral membrane of hepatocytes and is putatively involved in facilitating the hepatic uptake of unconjugated bilirubin [9, 13]. The exon 5 mutation of 388G > A is a nonsynonymous SNP with a missense mutation that causes a single amino acid substitution (asparagine to aspartic acid) encoded by codons 130 of SLCO1B1. This change disrupts the function of the SLCO1B1 transporter in transporting bilirubin from blood to the liver, leading to bilirubin concentration remaining high in the blood, which, in turn, causes hyperbilirubinemia .
Aside from the well-established risk factors of low gestational age and low birth weight, we found a significant association between the severity of hyperbilirubinemia and difference of populations across Jakarta, Bengkulu, and Papua, as well as with delivery location.
SLCO1B1*1B in both the heterozygous and homozygous variants was predominantly found (94.32%) in the neonates that we studied. Compared with the study of UGT1A1*60 and UGT1A1*6 from our previous study, SCLO1B1*1B was similar to *60 which showed a high incidence albeit with no significant association with the higher bilirubin level (severity) in neonatal hyperbilirubinemia . Other studies on SCLO1B1 in different countries also showed similar findings to our study concerning the incidence of polymorphisms. In a study in India, there were 77.7% hyperbilirubinemic neonates having the variant *1B in hyperbilirubinemic neonates, and 87.7% in Malaysia. All these studies demonstrate that more than half of the studied population has the *1B variant [14, 16, 21]. Both Asian countries share a lot in common in terms of genetic makeup with Indonesia.
Regarding the association between SLCO1B1*1B and higher bilirubin level in infants, there are still inconclusive evidence from different studies of different populations. In a study by Huang et al. in 2004, they reported the *1B SNP to be more frequent among hyperbilirubinemic infants with an OR of 2.01. In the study by Liu et al., the Chinese variant of *1B was also found to be associated with hyperbilirubinemia in the Guangdong population, but not in the Yunnan population . Whereas other studies in Fujian and Taiwan also failed to prove the association of *1B with neonatal hyperbilirubinemia, one study in Taiwan showed that Taiwanese neonates with the minor allele of *1B were at high risk to develop severe hyperbilirubinemia [22,23,24].
Some ethnicities are shown to be more at risk for severe hyperbilirubinemia than others such as African-American and Asian . This is most likely caused by the different prevalence of genetic polymorphisms responsible for metabolism or transport pathways . Variations of specific genes can also differ between the individuals in a smaller population group. An example of this is the intra-ethnic differences of UGT1A1 polymorphism in the Chinese populations of Han, She, and Dong origin . Global analysis of SLCO1B1 also suggests that SLCO1B1 diversity is greater within populations than between populations . Conversely, we did not find any intra-ethnic differences in our three populations for SLCO1B1 polymorphism.
From the comparison of genotypic distribution and clinical factor among the three different centers comprising population of different sub-ethnic origins, we observed a significant association with the severity of hyperbilirubinemia independent from the history of exclusive breastfeeding and other established clinical risk factors. This suggests the presence of a potentially undiscovered important risk factor for the development severe hyperbilirubinemia in the Indonesian population. Further research effort in Indonesia is warranted to investigate other potential genetic factors in different pathways leading to hyperbilirubinemia, including BLVRA and HMOX1 [29,30,31,32].
Our study is limited by the fact that we used a bank of DNA from a larger project, in which some subjects’ DNA samples were not sufficient for genetic analysis. Therefore, we were unable to study the whole population that we originally collected from the three centers. We also did not use consistently homogenous populations in the comparisons, and this could have obscured certain ethnicities where a different polymorphism variant could be dominant instead of *1B. The method of SNP analysis using PCR-RFLP is also less efficient and effective compared to other techniques such as DNA sequencing. However, we believe it is still sufficient in order to obtain national preliminary data that would have been otherwise severely lacking.
SLCO1B1*1B was found to be a common occurrence in Indonesian hyperbilirubinemic neonates. There was no statistically significant differences between occurrence of the SLCO1B1*1B variant and severity of neonatal hyperbilirubinemia. There was, however, a significant difference in the severity depending on the center where the neonate was enrolled. Future studies are required that includes a more diverse population to represent the entirety of the ethnically-diverse Indonesian population and exploring additional genetic risk factors is highly recommended.
Availability of data and materials
The datasets generated and/or analysed during the current study are not publicly available due to regulation by the hospital but are available from the corresponding author on reasonable request.
Cipto Mangunkusumo Hospital
polymerase chain reaction–restriction fragment length polymorphism
Solute Carrier Organic Anion Transporter Family Member 1b1
Statistical package for the social sciences
Total serum bilirubin
Porter ML, Dennis BL. Hyperbilirubinemia in the term newborn. Am Fam Physician. 2002;65(344):599–605.
Lauer BJ, Spector ND. Hyperbilirubinemia in the newborn. Pediatr Rev. 2011;32(8):341–9.
Watchko JF, Tiribelli C. Bilirubin-induced neurologic damage—mechanisms and management approaches. N Engl J Med. 2013;369(21):2021–30.
Rohsiswatmo R, Amandito R. Hiperbilirubinemia pada neonatus> 35 minggu di Indonesia; pemeriksaan dan tatalaksana terkini. Sari Pediatri. 2018;20(2):115–22.
Maruo Y, Nishizawa K, Sato H, Doida Y, Shimada M. Association of neonatal hyperbilirubinemia with bilirubin UDP-glucuronosyltransferase polymorphism. Pediatrics. 1999;103(6):1224–7.
Travan L, Lega S, Crovella S, Montico M, Panontin E, Demarini S. Severe neonatal hyperbilirubinemia and UGT1A1 promoter polymorphism. J Pediatr. 2014;165(1):42–5.
Kaplan M, Slusher T, Renbaum P, et al. n UDP-glucuronosyltransferase 1A1 promoter polymorphism in Nigerian neonates. Pediatr Res. 2008;63(1):109–11.
Beutler E, Gelbart T, Demina A. Racial variability in the UDP-glucuronosyltransferase 1 (UGT1A1) promoter: a balanced polymorphism for regulation of bilirubin metabolism? Proc Natl Acad Sci. 1998;95(14):8170–4.
Liu J, Long J, Zhang S, Fang X, Luo Y. The impact of SLCO1B1 genetic polymorphisms on neonatal hyperbilirubinemia: a systematic review with meta-analysis. J Pediatr. 2013;89(5):434–43.
Johnson AD, Kavousi M, Smith AV, et al. Genome-wide association meta-analysis for total serum bilirubin levels. Hum Mol Genet. 2009;18(14):2700–10.
Kim S-R, Saito Y, Sai K, et al. Genetic variations and frequencies of major haplotypes in SLCO1B1 encoding the transporter OATP1B1 in Japanese subjects: SLCO1B1* 17 is more prevalent than* 15. Drug Metab Pharmacokinet. 2007;22(6):456–61.
Liu J, Long J, Zhang S, Fang X, Luo Y. Polymorphic variants of SLCO1B1 in neonatal hyperbilirubinemia in China. Ital J Pediatr. 2013;39(1):49.
Lu A, Zhong D, Xie X. Relationship between SLCO1B1 polymorphism and neonatal hyperbilirubinemia. Guangxi Med J. 2015;37:592–6.
D'Silva S, Colah RB, Ghosh K, Mukherjee MB. Combined effects of the UGT1A1 and OATP2 gene polymorphisms as major risk factor for unconjugated hyperbilirubinemia in Indian neonates. Gene. 2014;547(1):18–22.
Prachukthum S, Nunnarumit P, Pienvichit P, et al. Genetic polymorphisms in Thai neonates with hyperbilirubinemia. Acta Paediatr. 2009;98(7):1106–10.
Wong F, Boo N, Othman A. Risk factors associated with unconjugated neonatal hyperbilirubinemia in Malaysian neonates. J Trop Pediatr. 2013;59(4):280–5.
Paludetto R, Mansi G, Raimondi F, et al. Moderate hyperbilirubinemia induces a transient alteration of neonatal behavior. Pediatrics. 2002;110(4):e50.
Farrell C-JL, Carter AC. Serum indices: managing assay interference. Ann Clin Biochem. 2016;53(5):527–38.
Amandito R, Putradista R, Jikesya C, et al. UGT1A1 gene and neonatal hyperbilirubinemia: a preliminary study from Bengkulu, Indonesia. BMC Res Notes. 2018;11(1):172.
Mishra S, Agarwal R, Deorari AK, Paul VK. Jaundice in the newborns. Indian J Pediatr. 2008;75(2):157.
Kohlrausch FB, de Cássia Estrela R, Barroso PF, Suarez-Kurtz G. The impact of SLCO1B1 polymorphisms on the plasma concentration of lopinavir and ritonavir in HIV-infected men. Br J Clin Pharmacol. 2010;69(1):95–8.
Weng Y-H, Chiu Y-W, Cheng S-W, Yang C-Y. Risk assessment of gene variants for neonatal hyperbilirubinemia in Taiwan. BMC Pediatr. 2016;16(1):144.
Zhou J, Yang C, Zhu W, et al. Identification of genetic risk factors for neonatal hyperbilirubinemia in Fujian Province, southeastern China: a case-control study. Biomed Res Int. 2018;2018.
Lin R, Wang X, Wang Y, et al. Association of polymorphisms in four bilirubin metabolism genes with serum bilirubin in three Asian populations. Hum Mutat. 2009;30(4):609–15.
Setia S, Villaveces A, Dhillon P, Mueller BA. Neonatal jaundice in Asian, white, and mixed-race infants. Arch Pediatr Adolesc Med. 2002;156(3):276–9.
Ramamoorthy A, Pacanowski M, Bull J, Zhang L. Racial/ethnic differences in drug disposition and response: review of recently approved drugs. Clin Pharmacol Ther. 2015;97(3):263–73.
Zhang A, Xing Q, Qin S, et al. Intra-ethnic differences in genetic variants of the UGT-glucuronosyltransferase 1A1 gene in Chinese populations. Pharmacogenomics J. 2007;7(5):333.
Pasanen MK, Neuvonen PJ, Niemi M. Global analysis of genetic variation in SLCO1B1; 2008.
Huang CS, Huang MJ, Chen ES, Chang PF, Chen WC. Glucose-6-phosphate dehydrogenase deficiency, the UDP-glucuronosyl transferase 1A1 gene, and neonatal hyperbilirubinemia. Gastroenterology. 2002;123(1):127–33.
Tiwari PK, Sethi A, Basu S, Raman R, Kumar A. Heme oxygenase-1 gene variants and hyperbilirubinemia risk in north Indian newborns. Eur J Pediatr. 2013;172(12):1627–32.
Chou H-C, Chen M-H, Yang H-I, et al. 211 G to a variation of UDP-glucuronosyl transferase 1A1 gene and neonatal breastfeeding jaundice. Pediatr Res. 2011;69(2):170–4.
Li Y, Wu T, Chen L, Zhu Y. Associations between G6PD, OATP1B1 and BLVRA variants and susceptibility to neonatal hyperbilirubinaemia in a Chinese Han population. J Paediatr Child Health. 2019;55(1):1–7.
We would like to thank the laboratory staff Bude Narti and NICU nursing staff of CMH, M. Yunus hospital, and Biak hospital for their contribution to this study. We would extend our utmost gratitude to Prof. Claudio Tiribelli, MD, PhD of Liver Research Center, Italian Liver Foundation for his scientific expertise. We would also give our appreciation to Alex Qinyang Liu, MBBS of Department of Surgery, Prince of Wales Hospital, Hong Kong and Nick Heng, PhD of Faculty of Dentistry, University of Otago for their invaluable English polishing assistance of our manuscript.
Funding was supported by Universitas Indonesia Research Grant for Q1 and Q2 International Journal Publication 2019 (No. NKB-0194/UN2.R3.1/HKP.05.00/2019) to AM and RR, and partly by Ministry of Research, Technology and Higher Education Republic of Indonesia Research Grant Hibah PUPT 2017 (No. 2601/UN2.R3.1/HKP05.00/2017) to AM. The funding body has no role in the design of the study and collection, analysis, nor interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All authors hereby declare that all experiments have been examined and approved by the Ethics Committee of the Faculty of Medicine, Universitas Indonesia (No. 995/UN2.F1/ETIK/2016, protocol number 16–11-432) and have therefore been performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki. A written informed consent was obtained from the parents before inclusion into the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Amandito, R., Rohsiswatmo, R., Halim, M. et al. SLCO1B1 c.388A > G variant incidence and the severity of hyperbilirubinemia in Indonesian neonates. BMC Pediatr 19, 212 (2019). https://doi.org/10.1186/s12887-019-1589-1
- Indonesia, neonatal hyperbilirubinemia