Profiling of UGT1A1*6, UGT1A1*60, UGT1A1*93, and UGT1A1*28 Polymorphisms in Indonesian Neonates With Hyperbilirubinemia Using Multiplex PCR Sequencing
Radhian Amandito1,2* 
Rinawati Rohsiswatmo2 
Erica Carolina3 
Rizka Maulida4 
Windhi Kresnawati5 
Amarila Malik3*- 1Neonatal Intensive Care Unit, Pondok Indah General Hospital, Jakarta, Indonesia
- 2Division of Perinatology, Department of Pediatrics, Faculty of Medicine, Cipto Mangunkusumo General Hospital, Universitas Indonesia, Jakarta, Indonesia
- 3Division of Pharmaceutical Microbiology and Biotechnology, Faculty of Pharmacy, Universitas Indonesia, Depok, Indonesia
- 4Department of Epidemiology, Faculty of Public Health, Universitas Indonesia, Depok, Indonesia
- 5Neonatology Unit, Biak General Hospital, Biak Numfor, Indonesia
Background: Single nucleotide polymorphism (SNP) variants of the uridine diphosphate glucuronosyltransferase 1A1 (UGT1A1) gene have been studied as an important factor in neonatal hyperbilirubinemia (jaundice) severity. Specific ethnicities, including Asians, have certain SNPs that appear more frequently than others.
Aim: To identify the most common SNPs in Indonesian neonates and their association with the severity of neonatal hyperbilirubinemia.
Methods: Eighty-eight inborn and outborn jaundiced infants from three different hospitals (Bengkulu, Jakarta, Biak Papua) across Indonesia were enrolled in this cross-sectional study and their peak total serum bilirubin (TSB) levels assessed. SNP variant analyses of the TATAA box, promoter, and exon 1 regions of UGT1A1 gene from 78 of the 88 infants were carried out using the SNaPshotR Multiplex Polymerase Chain Reaction (PCR) System followed by DNA sequencing.
Results: We detected SNP variants UGT1A1*28, UGT1A1*60, UGT1A1*93, and UGT1A1*6 in our population. Mean total serum bilirubin (TSB) was 14.59 ± 5.57 mg/dL. Bivariate analyses using delivery location, gestational age, birth weight, mother's age, and ethnicity were shown to be associated with moderate-to-severe hyperbilirubinemia (p < 0.05). None of the four SNPs appeared to be associated with moderate-to-severe hyperbilirubinemia. In multivariate analysis, however, only the “other ethnic group” (e.g., Chinese, Bengkulu, Papua, Bima) category showed an association with moderate-to-severe hyperbilirubinemia, with an odds ratio of 6.49 (95% CI 1.01–41.67; p < 0.05).
Conclusions: We found that the UGT1A1*60 is the most common SNP detected in neonates with hyperbilirubinemia in the Indonesian population. Interestingly, in Indonesia, UGT1A1 polymorphisms do not appear to be associated with differences in the severity of hyperbilirubinemia.
Introduction
Neonatal hyperbilirubinemia (jaundice) is commonly-encountered in daily medical practice. Availability of early detection of hyperbilirubinemia enables clinicians to treat these patients promptly with phototherapy. This improvement in healthcare reduces the neonatal mortality and morbidity rates of jaundice-related conditions. In most cases, the neonatal jaundice is physiological. It is a benign transitional phenomenon related to immaturity and hemolysis of red blood cells, liver, and gastrointestinal system, in which the problem resolves with little to no therapy (1–3). However, in some cases, jaundice is pathological and with rapidly-increasing levels of TSB (total serum bilirubin) that requires urgent therapy and scrutinized observation (3). On studies have shown that the risk factors for hyperbilirubinemia consist of hemolytic diseases, as well as environmental and genetic origins, especially when multiple factors are co-expressed (1, 3, 4). The most common gene variants reported to be associated are those encoding for the enzyme G6PD (glucose-6-phosphate dehydrogenase), OATP/SLCO1B1 (hepatic solute carrier organic anion transporter 1B1), and UGT1A1 (uridine diphosphate glucuronosyltransferase 1A1) (1, 5). UGT1A1 has been widely-studied in different populations and the findings appear to differ between ethnicities (6–9).
The human UGT1A1 is the enzyme responsible for bilirubin conjugation with glucuronic acid. UGT1A1 genetic variants that lead to reduced activity and decreased expression of the enzyme have been shown to be associated with non-hemolytic hyperbilirubinemia syndromes such as Gilbert's syndrome (GS) and Crigler-Najjar's syndrome types I and II (CN I and CN II, respectively) (1, 10–15). A variant of the TATAA box that has an additional TA insertion (A(TA)7TAA) (UGT1A1*28) was found to be associated with GS and CNs, with the mutation being found most frequently in Caucasian and African populations (1, 13, 16–19). In the East Asian population, however, UGT1A1*6 (211G > A) is the most prevalent single nucleotide polymorphism (SNP) (16, 17, 20, 21). This variant leads to an amino acid substitution of glycine-71 to arginine. In the promoter region of the gene, variant UGT1A1*60 (3279T > G) was found to be a common SNP in Egyptian neonates with neonatal hyperbilirubinemia (22). UGT1A1*60 is a SNP located in the phenobarbital-responsive enhancer module (gtPBREM) associated with decreased transcriptional activity of the UGT1A1 promoter (11, 12, 18). More recently, UGT1A1*93 (3156G > A), a SNP located between the TATAA box and gtPBREM, was found to show a high extent of linkage disequilibrium between UGT1A1*93 and UGT1A1*28 (11).
We have previously found that the UGT1A1*6 and UGT1A1*60 SNPs were prevalent in the ethnic Bengkulu population in Indonesia (23). However, the prevalence of the UGT1A1*28 and UGT1A1*93 SNPs have not been previously investigated in Indonesia, with the only other studies performed with the Malaysian population (11, 13, 14). This study aimed to identify the most common SNPs associated with neonatal jaundice in the Indonesian population.
Materials and Methods
Study Population
This was a cross-sectional, multi-center study conducted in three hospitals across Indonesia: (i) M. Yunus General Hospital in Bengkulu province; (ii) Biak General Hospital in West Papua province; and (iii) Cipto Mangunkusumo Hospital in Jakarta. The data collection process lasted 12 months (September 2016 to September 2017). Eighty-eight neonates who were diagnosed with clinical jaundice (based on Kramer's index) (24). by a neonatologist between days 3–7 post-natal were enrolled. Both pre-term and full-term, in-born and out-born neonates were included. We applied exclusion criteria in order to obtain unbiased data to prove the extent of genetic factors. Exclusion criteria includes neonates with signs of hemolytic anemia including those caused by ABO incompatibility (assessed through blood grouping of mother and neonate, direct Coomb's test of neonate, and fall in Hemoglobin of neonate), cephalhematoma, sepsis, major congenital malformations, neonates with conjugated bilirubin >20% of the TSB, and G6PD deficiency. These conditions were evaluated through clinical and laboratory examinations. Out of the 88 neonates that were originally enrolled, several had very low levels of DNA present in the sample, thus complicating the SNP analysis. Peak TSB was measured from venous blood (taken between day 3 and 7 before phototherapy was administered) using the ADVIA Chemistry Total Bilirubin 2 device (25). Peak TSB levels <13 mg/dL indicated mild hyperbilirubinemia, and moderate-to-severe hyperbilirubinemia was indicated by levels ≥13 mg/dL (26). We used this classification based on Paludetto et al. (26) due to the slight variation period when the TSB was taken. Other variables obtained included gestational age, birth weight, sex, mother's age at delivery, delivery method, delivery location, exclusive breastfeeding, history of older sibling requiring phototherapy, and the parents' ethnicity.
The study was approved by the Health Research Ethics Committee—Universitas Indonesia and Cipto Mangunkusumo Hospital. We obtained prior consent from the parents of all included neonates in the study.
Detection of SNPs
Peripheral blood samples were stored at −20°C in EDTA-containing vacutainer tubes until required. Genomic DNA (gDNA) was purified from blood samples using the QIAamp DNA Blood Mini Kit (QIAGEN, Germany) as described previously (23). The gDNAs were subsequently used for multiplex PCR experiments and SNP genotyping.
Analysis of DNA Sample by PCR
All oligonucleotide primers used in this study are listed in Supplementary Table 1. Genotyping of three SNPs in the promoter and coding regions of UGT1A1, namely UGT1A1*60 (promoter, 3279T>G), UGT1A1*93 (promoter, 3156G>A), and UGT1A1*6 (Exon 1, 211G>A) was carried out using the SNaPshotR Multiplex kit. Genotyping of the TATAA box sequence repeat (UGT1A1*28) in the UGT1A1 promoter region was determined by cloning and sequencing of PCR amplicons.
A total volume 25 μl of PCR Mix was prepared which consisted of 12.5 μl PCR Master Mix MyFi MixTM (Bioline, USA), 0.5 μl forward primer (10 μM), 0.5 μl reverse primer (10 μM), 1.0 μl DNA template (50 ng/μl), and 10.5 μl nuclease-free water. PCR was carried out using the following parameters: initial denaturation of 5 min at 95°C followed by 30 cycles of denaturation of 95°C for 30 s, annealing at 55°C for 30 s, extension at 72°C for 1 min, and a final polymerization at 72°C for 10 min. The presence of amplicons was confirmed by agarose gel electrophoresis (Supplementary Figure 1) prior to DNA sequencing. DNA sequence data were analyzed using the GeneMapper 4.0 (Applied Biosystems) system (1st BASE Sdn. Bhd., Malaysia).
Amplicons for cloning were amplified in 25 μl of PCR mix as follows; 12.5 μl PCR Master Mix (KOD FX Neo, Toyobo), 1.0 μl Forward primer (10 μM), 1.0 μl Reverse primer (10 μM), 1.0 μl DNA template (50 ng/μl), and 9.5 μl nuclease free water. PCR amplification was carried out as follows: initial denaturation of 3 min at 95°C followed by 35 cycles of denaturation of 95°C for 20 s, annealing at 55°C for 30 s, extension at 68°C for 30 min, and final extension at 68°C for 5 min. One microliter aliquots of individual unpurified 207-bp amplicons (Supplementary Figure 2) were then inserted into the pJET1.2 cloning vector (ThermoFisher, USA) and transformed into competent Escherichia coli DH5α cells according to the manufacturer's protocol (27). Colony PCR was used to screen for positive clones and recombinant plasmids purified from 18-h cultures were sequenced with both UGT1A1*28_F and UGT1A1*28_R primers (1st BASE, Malaysia).
Statistical Analysis
STATA version 12 (Macintosh version) was used for data management and statistical analyses. The variables were first presented descriptively with 88 subjects. However, only 78 subjects were analyzed further for bivariate and multivariate analyses as 10 participants had to be excluded due to incomplete data. Bivariate analysis was conducted between independent and dependent variables using chi-square/Fisher's exact test, Student's t-test, and Kruskal-Wallis. Variables with p < 0.25 were included in the multivariate analysis using logistic regression. We used two-sided p-values in our analysis with a p < 0.05 level of significance. Gender, hospital, ethnic group, exclusive breastfeeding, delivery method, delivery location, sibling requiring phototherapy, ethnic groups, and SNPs were categorical data, while gestational age, birth weight, and mother's age were expressed as numerical data. Ethnicity was based on ethnic group majority.
Results
Clinical characteristics of 88 neonates from all three study centers (Bengkulu, Jakarta, and Biak) are shown in Table 1. There was an almost equal proportion of female (54.5%) and male neonatal subjects (45.5%). More than half of the neonates were breastfed exclusively (60.2%), about half of the newborns were born via Cesarean-section (56.8%), and almost two-thirds were delivered in hospital (65.9%). Most of the neonates (96.6%) had 335 no sibling who required phototherapy. The neonates originated from many different ethnic groups with no dominant ethnic group. The average gestational age was 34.5 ± 3.7 weeks, average birth weight was 2189.1 ± 849.8 g, and average mother's age was 29.3 ± 6.2 years old. The mean peak TSB was 14.59 ± 5.57 mg/dL, and the number of neonates with severe hyperbilirubinemia (TSB > 25 mg/dL) is 4. Our bivariate analysis results of the clinical characteristic are shown in Table 2. Delivery location, gestational age, birth weight, mother's age, and ethnic groups are shown to be associated with moderate-to-severe hyperbilirubinemia (p < 0.05).
Table 1. Clinical characteristics of the neonates (n = 88).
Table 2. Variables associated with the severity of hyperbilirubinemia.
The peak TSB measurements in neonates in the wild-type group and all UGT1A1 SNP variants are shown in Table 3. UGT1A1*6 only yielded 84 results with most neonates were wild-type (90.48%), UGT1A1*60 yielded 86 results with half of the neonates being wild-type (50%), and UGT1A1*93 only yielded 82 results with majority of neonates being wild-type (80.49%). In UGT1A1*28, almost all neonates were wild-type (97.73%). Bivariate analysis results of SNPs and hyperbilirubinemia are shown in Table 4. None of the SNPs showed any association with moderate-to-severe hyperbilirubinemia.
Table 3. Frequency of SNPs detected.
Table 4. SNP distribution in mild and moderate-severe hyperbilirubinemia groups.
As shown in Table 5, variables with p > 0.25 from bivariate analyses (Tables 2, 4) were put into multivariate analysis using logistic regression. Only the minority, i.e. others (Chinese, Bengkulu, Papua, Bima) ethnic group showed association with moderate-to-severe bilirubinemia (p = 0.049) with a POR of 6.49 compared to Javanese neonates.
Table 5. Variables used in logistic regression analysis.
UGT1A1*60 (promoter, 3279T > G), UGT1A1*93 (promoter, 3156G > A), and UGT1A1*6 (Exon 1, 211G > A), as well as TATAA box sequence repeat of promoter, UGT1A1*28, were considered the most frequently reported SNPs as previously studied in several Asian countries (5, 10, 28). We also previously studied UGT1A1*60 and UGT1A1*6 in Bengkulu, Sumatra population performing PCR-RFLP (23). Other three SNPs, i.e., c.686C4A, c.1091C4T, and c.1456T4G, were not detected.
Discussion
The present investigation has confirmed and expanded on our previous findings on the association of UGT1A1*60 and UGT1A1*6 with neonatal hyperbilirubinemia in Indonesia (23). In addition, the present SNP analysis has also included two further SNPs, UGT1A1*28, and UGT1A1*93 to determine their prevalence, if any, in the condition. Furthermore, we have increased the ethnic diversity of patients to ensure a more accurate representation of Indonesia's multi-ethnic population.
Certain characteristics of our population, e.g., preterm mean gestational age, low mean birth weight, and more exclusive breastfeeding, were as expected (29). Previous studies reported that premature birth is associated with a defective UGT1A1 protein causing immaturity of the conjugating enzyme, therefore increasing the risk of severe hyperbilirubinemia (30, 31). In addition, prematurity is also associated with higher red blood cell, hepatocyte, and gastrointestinal cell immaturity leading to hemolysis and, consequently, a build-up of bilirubin (24, 32). Exclusive breastfeeding is also a very well-known risk factor for hyperbilirubinemia in the form of breastmilk jaundice (2, 10). Breastmilk jaundice may occur due to certain factors present in human breastmilk (i.e., IL 1ß, IL6, ß-glucuronidase, epidermal growth factor, alpha-fetoprotein), which inhibits the conjugation of bilirubin that facilitates its excretion (33). Chou et al. (10) reported an association between UGT1A1 polymorphisms and breastfeeding jaundice (10). They concluded that both factors have to be present in order for higher development of hyperbilirubinemia to occur. This might be why not all neonates harboring the UGT1A1*28 polymorphism and with exclusive breastfeeding suffer from hyperbilirubinemia, and vice versa. Our current findings did not support this hypothesis. In addition, our term and healthy infants start breastfeeding in the delivery room when considered clinically stable. The pathophysiology of breastmilk jaundice involves the occurrence of weight loss after birth subsequent to fasting (34). Through immediate breastfeeding, we suspect that the expected increase in TSB due to breastmilk jaundice does not occur in our population. This might explain why breastfeeding is not significant in our study. Bivariate analysis also revealed that the delivery location [OR 5.56 (1.38–31.91)] and the mother's age were significant. Infants born in a hospital are more routinely examined for comorbidities including hyperbilirubinemia during the early days post-natal when they are still in hospital. It is therefore safer for mothers who have known risk factors, or history of a child with severe hyperbilirubinemia, to deliver their child in a hospital as they can identified early and treated appropriately. In North America, hyperbilirubinemia is the most important reason for readmission of neonates to the hospital (35). With regard to the mother's age, this has always been considered as a risk factor and predictive of extreme neonatal hyperbilirubinemia (36). In contrary, our study showed that infants of an older mother appears to be more protected against moderate-to-severe hyperbilirubinemia. We have to consider that the maternal health and antenatal examinations are inadequate in terms of quality in Indonesia (37). There may also be more younger mothers who do not attend routine antenatal visits due to lack of knowledge (38). In the Indonesian population, it is likely that a combination of different risk factors eventually leads to a higher risk of moderate-to-severe neonatal hyperbilirubinemia.
In previous in vitro studies, the UGT1A1*6 mutation has been proven to have decreased isozyme activity (5). Heterozygous forms display 60% lower activity, while the activity in homozygous forms is even less at <35%. Whereas, the presence of UGT1A1*6 is known to be a risk factor for mild adult hyperbilirubinemia (GS) in Caucasian and East Asian populations (16, 19, 28, 39), the prevalence of this SNP is low in the Southeast Asian population (17, 23, 28). Our study confirms this with higher occurrence of wild-type UGT1A1 than the UGT1A1*6 SNP (90.5 vs. 9.5%) with an OR of 0.72 (0.11–5.36). When comparing between the mild and moderate-to-severe hyperbilirubinemia groups, we also did not find a significant association (p = 0.69). This shows that UGT1A1*6 is not prevalent, and that it is not a contributing factor toward the severity of hyperbilirubinemia in Southeast Asian populations.
We found a similar frequency of UGT1A1*60 SNPs in our population compared to wild-type. Those with this SNP are widely-known to have decreased transcriptional activity of UGT1A1 (about 60%) via an impaired gtPBREM (12, 18). This finding, however, contrasts what was reported for the Malaysian population, in which a higher frequency of UGT1A1*60 SNPs was observed (11). Interestingly, whilst the peak TSB measurements were higher in wild-type compared to the *60 variants, the differences were not statistically significant (p = 0.652). SNPs of *60 was found to co-exist with *28 (40). It has been hypothesized that co-existence of both *60 and *28 SNPs would decrease the transcriptional activity of UGT1A1 by 70% (18). However, our study shows that no association could be found between individual nor multiple SNPs with developing moderate-to-severe hyperbilirubinemia compared to mild hyperbilirubinemia [OR of 1.05 (0.36–3.05) for heterozygotes and 2.4 (0.21–125.93) for homozygotes].
Another commonly found SNP combination that has been reported is UGT1A1*93 and UGT1A1*60. The *93 SNP can be found not only in hyperbilirubinemic neonates, but also healthy neonates (4). The co-existence with *60 SNPs can also explain this marked increase in TSB level. In our study, we found only a few number of *93 SNPs and their existence was not predictive of moderate-to-severe hyperbilirubinemia in Indonesian neonates with an OR of 1.01 (0.26–4.31) (p = 0.98). This is consistent with findings from a Chinese study and a study comparing Uzbekistan and Japanese populations, in which wild-type *93 were more frequent than their SNPs (41, 42).
In previous studies, the enzyme activity in UGT1A1 with TA insertions in the normal sequence A(TA)6TAA of the TATAA box promoter region significantly reduces enzyme activity (60–80% of normal) (2, 11, 39). However, in another functional study, longer TA repeats appear to be a protective factor for hyperbilirubinemia. According to the American Academy of Pediatrics, whereas Caucasian and Asian populations are more susceptible to hyperbilirubinemia due to the presence of longer TA repeats, African populations are less susceptible (43).
On the other hand, fewer TA repeats, which causes an overall increase in UGT1A1 expression, has been reported to reduce TSB levels in African neonates (4, 44). However, data on Southeast Asian populations remain inconclusive. For example, Yusoff et al. (14) showed a low incidence of this SNP, i.e., no significant association, with severe hyperbilirubinemia in Malaysian neonates (14). This is supported by our data for the Indonesian population which showed only one sample exhibiting this SNP (having moderate-to-severe hyperbilirubinemia). Considering our relatively low sample size, this result should be addressed, not only with a much larger population size, but also a more ethnically-diverse one.
UGT1A1*28 and UGT1A1*60 have been demonstrated to show strong linkage disequilibrium (40), and co-expression in these two promoter regions decreases the transcriptional activity to 37% of normal, more so when the homozygous alleles are present (30% of normal) (18). In our study, very few *28 SNPs were found and therefore we could not prove this effect of this SNP combination. We should always keep in mind however, that the clinical manifestation of hyperbilirubinemia, especially moderate-to-severe hyperbilirubinemia, would not be present unless active hemolysis is occurring. Therefore, it is of vital importance that assessment of the causes of hyperbilirubinemia in infants also always includes factors related with hemolysis.
From the bivariate analysis, we included several variables that are eligible for multivariate analysis. We found that the only significant factor is ethnicity, specifically the groups comprising other ethnicities (Chinese, Bengkulu, Papua, and Bima; p = 0.049). Due to the small number of samples from individual ethnic groups, we were unable to identify which ethnic group has a higher risk factor toward moderate-to-severe hyperbilirubinemia in Indonesia. What we can infer from this is that there are intra-ethnic differences of UGT1A1 polymorphisms in the Indonesian population, not unlike that observed in the Han, She, and Dong populations in China (45). However, we found that there is no significant association between ethnicity and all of the SNPs that we studied (p > 0.05). This prompts further investigation on what could differentiate between the ethnicities, both from their clinical and genetic factors that could be involved in the development of hyperbilirubinemia.
UGT1A1 is not only involved in neonatal hyperbilirubinemia, but certain variants are also important in drug delivery. One of them, UGT1A1*60, is located in the gtPBREM. The assessment of this genetic variation may be useful to consider phenobarbital treatment depending on their respective ethnic group populations. Considering the varied genetic makeup of Indonesia and Southeast Asian populations, even compared to Indian and East Asian populations, further studies in not only UGT1A1, but other genes, e.g., HMOX1, BLVRA, and SLCO1B1, especially those closely-related with hyperbilirubinemia or drug delivery are warranted (46–48). Approach and treatment of infants with hyperbilirubinemia and other UGT1A1-related drug therapies should be based on the considerations of SNPs specific for Southeast Asian populations.
There were several limitations to the present study. Firstly, bilirubin measurement timing was not possible to be taken at all exactly the same time post-natal. However, we ensured that the TSB level that we analyzed was at its peak for each neonate. Secondly, the results of this study could not be generalized to all neonates in Indonesia due to the limited number of cases and very limited number of participating hospitals/provinces. Additionally, discrepancy between numbers of samples that were available for analysis in different groups might also cause some experimental bias. The study design that we prefer is case control study, but due to limitations in the execution of our research we only were able to conduct a cross-sectional study. We also excluded most of the risk factors of hyperbilirubinemia, however this becomes a limitation since many other studies concluded that other risk factors aside from genetic polymorphism need to be present to result in a severe hyperbilirubinemia. A more comprehensive future investigation with a much larger sample size from other regions of Indonesia, is needed to define a nationwide correlation of genomic variation and hyperbilirubinemia severity in Indonesia.
In conclusion, our current study found that UGT1A1*6, UGT1A1*60, UGT1A1*93, and UGT1A1*28 are present in Indonesian population. The distribution and frequency of SNPs in our Indonesian population is, however, different when compared to other Asian populations, as only UGT1A1*60 was found to be prevalent, i.e., in at least half of the neonates we studied. No SNP in our study was found to be associated with the severity of hyperbilirubinemia in Indonesian neonates. Out of all the clinical factors, only ethnic grouping was considered significant and therefore future studies of other genes relevant to hyperbilirubinemia in Indonesia should include a more ethnically-diverse population.
Data Availability
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.
Author Contributions
This work was carried out in collaboration between all authors. RA, RR, and AM designed the study. RA and WK acquired the data. RA, EC, and AM performed experiments. RA, RM, EC, RR, and WK, analyzed the data. RA, RM, and AM wrote the manuscript. EC, WK, and RR gave technical support and conceptual advice. All authors read and approved the final manuscript.
Funding
Funding was supported by a Hibah PUPT 2017 Research Grant (No. 2601/UN2.R3.1/HKP05.00/2017, Directorate for Higher Education, Republic of Indonesia) to AM. The funding body had no role in the design of the study and collection, analysis, interpretation of data, or in preparation of the manuscript.
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
The authors would like to thank the staff at our neonatal wards and NICUs especially Bude Narti for their helpful co-operation in this study. We would extend our utmost appreciation to Prof. Claudio Tiribelli, of the Italian Liver Foundation for scientific advice and support. We would also give our deepest gratitude to Nick Heng, University of Otago for invaluable English polishing assistance of our manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2019.00328/full#supplementary-material
Abbreviations
CN, Crigler-Najjar Syndrome; G6PD, glucose-6-phosphate dehydrogenase enzyme; gDNA, Genomic DNA; GS, Gilbert's syndrome; gtPBREM, phenobarbital-responsive enhancer module; OATP/SLCO1B1, hepatic solute carrier organic anion transporter 1B1; PCR-RFLP, Polymerase Chain Reaction-Restriction Fragment Length Polymorphism; POR, prevalence odds ratio; TSB, total serum bilirubin; UGT1A1, UDP-glucoronosyltransferase 1A1.
References
1. Watchko JF, Lin Z, Clark RH, Kelleher AS, Walker MW, Spitzer AR, et al. Complex multifactorial nature of significant hyperbilirubinemia in neonates. Pediatrics. (2009) 124:e868–77. doi: 10.1542/peds.2009-0460
2. Zhou Y, Wang SN, Li H, Zha W, Wang X, Liu Y, et al. Association of UGT1A1 variants and hyperbilirubinemia in breast-fed full-term Chinese infants. PLoS ONE. (2014) 9:e104251. doi: 10.1371/journal.pone.0104251
3. Lin Z, Fontaine J, Watchko JF. Coexpression of gene polymorphisms involved in bilirubin production and metabolism. Pediatrics. (2008) 122:e156–62. doi: 10.1542/peds.2007-3249
4. 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:18–22. doi: 10.1016/j.gene.2014.05.047
5. Huang MJ, Kua KE, Teng HC, Tang KS, Weng HW, Huang CS. Risk factors for severe hyperbilirubinemia in neonates. Pediatric Res. (2004) 56:682–9. doi: 10.1203/01.PDR.0000141846.37253.AF
6. Watchko JF, Daood MJ, Biniwale M. Understanding neonatal hyperbilirubinaemia in the era of genomics. In: Paper presented at: Seminars in Neonatology. Philadelphia, PA (2002). doi: 10.1053/siny.2002.0102
7. Yu Z, Zhu K, Wang L, Liu Y, Sun J. Association of neonatal hyperbilirubinemia with UGT1A1 gene polymorphisms: a meta-analysis. Med. Sci. Monitor. (2015) 21:3104. doi: 10.12659/MSM.894043
8. Kaniwa N, Kurose K, Jinno H, Tanaka-Kagawa T, Saito Y, Saeki M, et al. Racial variability in haplotype frequencies of UGT1A1 and glucuronidation activity of a novel single nucleotide polymorphism 686C> T (P229L) found in an African-American. Drug Metab. Disposit. (2005) 33:458–65. doi: 10.1124/dmd.104.001800
9. Tiwari PK, Bhutada A, Agarwal R, Basu S, Raman R, Kumar A. UGT1A1 gene variants and clinical risk factors modulate hyperbilirubinemia risk in newborns. J Perinatol. (2014) 34:120–4. doi: 10.1038/jp.2013.140
10. Chou HC, Chen MH, Yang HI, Su YN, Hsieh WS, Chen CY, et al. 211 G to a variation of UDP-glucuronosyl transferase 1A1 gene and neonatal breastfeeding jaundice. Pediatr Res. (2011) 69:170–4. doi: 10.1203/PDR.0b013e31820263d2
11. Yusoff S, Takeuchi A, Ashi C, Tsukada M, Ma'amor NH, Zilfalil BA, et al. A polymorphic mutation, c.-3279T> G, in the UGT1A1 promoter is a risk factor for neonatal jaundice in the Malay population. Pediatr Res. (2010) 67:401–406. doi: 10.1203/PDR.0b013e3181d22f78
12. Sugatani J, Mizushima K, Osabe M, Yamakawa K, Kakizaki S, Takagi H, et al. Transcriptional regulation of human UGT1A1 gene expression through distal and proximal promoter motifs: implication of defects in the UGT1A1 gene promoter. Naunyn-Schmiedeberg's Arch Pharmacol. (2008) 377:597–605. doi: 10.1007/s00210-007-0226-y
13. Sutomo R, Talib NA, Yusoff NM, Van Rostenberghe H, Sadewa AH, Sunarti, et al. Screening for G71R mutation of the UGT1A1 gene in the Javanese-Indonesian and Malay-Malaysian populations. Pediatr Int. (2004) 46:565–9. doi: 10.1111/j.1442-200x.2004.01959.x
14. Yusoff S, Van Rostenberghe H, Yusoff NM, Talib NA, Ramli N, Ismail NZ, et al. Frequencies of A (TA) 7TAA, G71R, and G493R mutations of the UGT1A1 gene in the Malaysian population. Neonatology. (2006) 89:171–6. doi: 10.1159/000088844
15. Huang MJ, Chen YC, Huang YY, Yang SS, Chen PL, Huang CS. Effect of UDP-glucuronosyltransferase 1A1 activity on risk for developing Gilbert's syndrome. Kaohsiung J Med Sci. 2019. doi: 10.1002/kjm2.12077
16. Liu W, Chang LW, Xie M, Li WB, Rong ZH, Wu L, et al. Correlation between UGT1A1 polymorphism and neonatal hyperbilirubinemia of neonates in Wuhan. J Huazhong Univer Sci Technol. (2017) 37:740–3. doi: 10.1007/s11596-017-1797-6
17. Agrawal SK, Kumar P, Rathi R, Sharma N, DAS R, Prasad R, et al. UGT1A1 gene polymorphisms in North Indian neonates presenting with unconjugated hyperbilirubinemia. Pediatr Res. (2009) 65:675–80. doi: 10.1203/PDR.0b013e31819ed5de
18. Mazur-Kominek K, Romanowski T, Bielawski K, Kiełbratowska B, Preis K, Domzalska – Popadiuk I, et al. Association between uridin diphosphate glucuronosylotranserase 1A1 (UGT1A1) gene polymorphism and neonatal hyperbilirubinemia. Acta Biochimica Polonica. (2017) 64:351–6. doi: 10.18388/abp.2016_1450
19. Travan L, Lega S, Crovella S, Montico M, Panontin E, Demarini S. Severe neonatal hyperbilirubinemia and UGT1A1 promoter polymorphism. J Pediatr. (2014) 165:42–5. doi: 10.1016/j.jpeds.2014.03.013
20. Mokhtar WA, Allam RM, Zidan N, Hamed MS. Association between UGT1A1 gene GLY71Arg polymorphism and severe unexplained indirect hyperbilirubinemia among neonates. Alex J Pediatr. (2018) 31:76. doi: 10.4103/AJOP.AJOP_18_18
21. Mehrad-Majd H, Haerian MS, Akhtari J, Ravanshad Y, Azarfar A, Mamouri G. Effects of Gly71Arg mutation in UGT1A1 gene on neonatal hyperbilirubinemia: a systematic review and meta-analysis. J Maternal-Fetal Neonatal Med. (2019) 32:1575–85. doi: 10.1080/14767058.2017.1410789
22. Tomerak RH, Helal NF, Shaker OG, Yousef MA. Association between the specific UGT1A1 promoter sequence variant (c-3279T>G) and unconjugated neonatal hyperbilirubinemia. J Trop Pediatr. (2016) 62:457–63. doi: 10.1093/tropej/fmw031
23. Amandito R, Putradista R, Jikesya C, Utaminingsih D, Rusin J, Rohsiswatmo R, et al. UGT1A1 gene and neonatal hyperbilirubinemia: a preliminary study from Bengkulu, Indonesia. BMC Res Notes. (2018) 11:172. doi: 10.1186/s13104-018-3284-y
24. Bhutani VK, Vilms RJ, Hamerman-Johnson L. Universal bilirubin screening for severe neonatal hyperbilirubinemia. J Perinatol. (2010) 30:S6. doi: 10.1038/jp.2010.98
25. Farrell CJ, Carter AC. Serum indices: managing assay interference. Ann Clin Biochem. (2016) 53:527–38. doi: 10.1177/0004563216643557
26. Paludetto R, Mansi G, Raimondi F, Romano A, Crivaro V, Bussi M, et al. Moderate hyperbilirubinemia induces a transient alteration of neonatal behavior. Pediatrics. (2002) 110:e50. doi: 10.1542/peds.110.4.e50
27. Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press (1989).
28. Wu XJ, Zhong DN, Xie XZ, Ye DZ, Gao ZY. UGT1A1 gene mutations and neonatal hyperbilirubinemia in Guangxi Heiyi Zhuang and Han populations. Pediatr Res. (2015) 78:585. doi: 10.1038/pr.2015.134
29. Khoury MJ, Calle EE, Joesoef RM. Recurrence risk of neonatal hyperbilirubinemia in siblings. Am J Dis Child. (1988) 142:1065–9. doi: 10.1001/archpedi.1988.02150100059026
30. Kawade N, Onishi S. The prenatal and postnatal development of UDP-glucuronyltransferase activity towards bilirubin and the effect of premature birth on this activity in the human liver. Biochem J. (1981) 196:257–60. doi: 10.1042/bj1960257
31. Kaplan M, Renbaum P, Levy-Lahad E, Hammerman C, Lahad A, Beutler E. Gilbert syndrome and glucose-6-phosphate dehydrogenase deficiency: a dose-dependent genetic interaction crucial to neonatal hyperbilirubinemia. Proc Natl Acad Sci USA. (1997) 94:12128–32. doi: 10.1073/pnas.94.22.12128
32. Engle WA, Tomashek KM, Wallman C. “Late-preterm” infants: a population at risk. Pediatrics. (2007) 120:1390–401. doi: 10.1542/peds.2007-2952
33. Bratton S, Stern M. Breast milk Jaundice. In: StatPearls [Internet]. Treasure Island, FL: StatPearls Publishing (2019).
34. Bertini G, Dani C, Tronchin M, Rubaltelli FF. Is breastfeeding really favoring early neonatal jaundice. Pediatrics. (2001) 107:E41. doi: 10.1542/peds.107.3.e41
35. Janssen PA, Saxell L, Page LA, Klein MC, Liston RM, Lee SK. Outcomes of planned home birth with registered midwife versus planned hospital birth with midwife or physician. CMAJ. (2009) 181:377–83. doi: 10.1503/cmaj.081869
36. Newman TB, Xiong B, Gonzales VM, Escobar GJ. Prediction and prevention of extreme neonatal hyperbilirubinemia in a mature health maintenance organization. Arch Pediatr Adol Med. (2000) 154:1140–7. doi: 10.1001/archpedi.154.11.1140
37. Dasar RK. Hasil Utama Riskesdas 2018. Kemenkes Balitbangkes. (2018). Available online at: http://dinus.ac.id/repository/docs/ajar/hasil-riskesdas-2018.pdf
38. Acuin CS, Khor GL, Liabsuetrakul T, Achadi EL, Htay TT, Firestone R, et al. Maternal, neonatal, and child health in southeast Asia: towards greater regional collaboration. Lancet. (2011) 377:516–25. doi: 10.1016/S0140-6736(10)62049-1
39. Long J, Zhang S, Fang X, Luo Y, Liu J. Neonatal hyperbilirubinemia and Gly71Arg mutation of UGT1A1 gene: a Chinese case–control study followed by systematic review of existing evidence. Acta Paediatr. (2011) 100:966–71. doi: 10.1111/j.1651-2227.2011.02176.x
40. Moyer AM, Skierka JM, Kotzer KE, Kluge ML, Black JL, Baudhuin LM. Clinical UGT1A1 genetic analysis in pediatric patients: experience of a reference laboratory. Mol Diagn Ther. (2017) 21:327–35. doi: 10.1007/s40291-017-0265-0
41. Maeda H, Hazama S, Shavkat A, Okamoto K, Oba K, Sakamoto J, et al. Differences in UGT1A1, UGT1A7, and UGT1A9 polymorphisms between Uzbek and Japanese populations. Mol Diag Ther. (2014) 18:333–42. doi: 10.1007/s40291-014-0083-6
42. Ma L, Chen Y, Yang C, Jiang H, Zhu J, Cheng Y. Association of UGT1A1 (* 28,* 60 and* 93) polymorphism with the adverse reactions of irinotecan chemotherapy in extensive stage small cell lung cancer. Zhonghua zhong liu za zhi. (2015) 37:29–32. doi: 10.3760/cma.j.issn.0253-3766.2015.01.006
43. Hyperbilirubinemia AAoPSo. Management of hyperbilirubinemia in the newborn infant 35 or more weeks of gestation. Pediatrics. (2004) 114:297. doi: 10.1542/peds.114.1.297
44. 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 USA. (1998) 95:8170–4. doi: 10.1073/pnas.95.14.8170
45. Zhang A, Xing Q, Qin S, Du J, Wang L, Yu L, et al. Intra-ethnic differences in genetic variants of the UGT-glucuronosyltransferase 1A1 gene in Chinese populations. Pharm J. (2007) 7:333. doi: 10.1038/sj.tpj.6500424
46. Kaplan M, Wong R, Stevenson D. Heme oxygenase-1 promoter polymorphisms: do they modulate neonatal hyperbilirubinemia? J Perinatol. (2017) 37:901. doi: 10.1038/jp.2017.6
47. 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) 1:1–7. doi: 10.1111/jpc.14346
Keywords: Indonesia, snapshot, polymorphism, UGT1A1, unconjugated hyperbilirubinemia
Citation: Amandito R, Rohsiswatmo R, Carolina E, Maulida R, Kresnawati W and Malik A (2019) Profiling of UGT1A1*6, UGT1A1*60, UGT1A1*93, and UGT1A1*28 Polymorphisms in Indonesian Neonates With Hyperbilirubinemia Using Multiplex PCR Sequencing. Front. Pediatr. 7:328. doi: 10.3389/fped.2019.00328
Received: 21 May 2019; Accepted: 22 July 2019;
Published: 07 August 2019.
Edited by:
Arjan Te Pas, Leiden University, NetherlandsReviewed by:
Ronald James Wong, Stanford University, United StatesChristian Victor Hulzebos, University Medical Center Groningen, Netherlands
Copyright © 2019 Amandito, Rohsiswatmo, Carolina, Maulida, Kresnawati and Malik. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Radhian Amandito, amandito.radhian@gmail.com; Amarila Malik, amarila.malik@ui.ac.id; amarila.malik@gmail.com