Bilirubin is independently associated with oxidized LDL levels in young obese patients
- Henrique Nascimento1, 2, 3,
- Ana Inês Alves1, 2,
- Susana Coimbra1, 2, 4,
- Cristina Catarino1, 2, 3,
- Diana Gomes2,
- Elsa Bronze-da-Rocha1, 2, 3,
- Elísio Costa1, 2, 3,
- Petronila Rocha-Pereira5,
- Luísa Aires6, 7,
- Jorge Mota6,
- Helena Ferreira Mansilha8,
- Carla Rêgo9,
- Alice Santos-Silva1, 2, 3 and
- Luís Belo1, 2, 3, 10Email author
© Nascimento et al.; licensee BioMed Central. 2015
Received: 11 November 2014
Accepted: 14 January 2015
Published: 23 January 2015
Bilirubin can prevent lipid oxidation in vitro, but the association in vivo with oxidized low-density lipoprotein (Ox-LDL) levels has been poorly explored. Our aim is to the association of Ox-LDL with total bilirubin (TB) levels and with variables related with metabolic syndrome and inflammation, in young obese individuals.
125 obese patients (13.4 years; 53.6% females) were studied. TB, lipid profile including Ox-LDL, markers of glucose metabolism, and levels of C-reactive protein (CRP) and adiponectin were determined. Anthropometric data was also collected.
In all patients, Ox-LDL correlated positively with BMI, total cholesterol, LDLc, triglycerides (TG), CRP, glucose, insulin and HOMAIR; while inversely with TB and HDLc/Total cholesterol ratio (P < 0.05 for all). In multiple linear regression analysis, LDLc, TG, HDLc and TB levels were significantly associated with Ox-LDL (standardized Beta: 0.656, 0.293, −0.283, −0.164, respectively; P < 0.01 for all). After removing TG and HDLc from the analysis, HOMAIR was included in the regression model. In this new model, LDLc remained the best predictor of Ox-LDL levels (β = 0.665, P < 0.001), followed by TB (β = −0.202, P = 0.002) and HOMAIR (β = 0.163, P = 0.010).
Lower bilirubin levels may contribute to increased LDL oxidation in obese children and adolescents, predisposing to increased cardiovascular risk.
KeywordsBilirubin Oxidized-LDL Oxidative stress Atherosclerosis Pediatric obesity
Bilirubin has important antioxidant and anti-inflammatory properties. Actually, bilirubin can protect low-density lipoprotein (LDL) from oxidation  and persons expressing Gilbert syndrome (i.e., unconjugated hyperbilirubinemia) present lower circulating levels of oxidized-LDL (Ox-LDL) . Stojanov et al.  reported that healthy young men (18 to 22 years old) with high total bilirubin (TB) have lower concentration of Ox-LDL. Furthermore, a protective association between serum bilirubin level and cardiovascular disease (CVD) and death exists in both men and women . Besides, serum bilirubin levels were inversely associated with body mass index (BMI) values in adults , and in a study involving American children and adolescents, bilirubin concentrations were found to be inversely associated with metabolic syndrome, which seemed to be related to the insulin resistance (IR) status .
Obese children and adolescents already present with altered CVD risk markers, namely higher IR, triglycerides (TG), LDL cholesterol (LDLc) and C-reactive protein (CRP) concentrations and lower high-density lipoprotein cholesterol (HDLc) and adiponectin levels compared to controls . It was also reported that Ox-LDL level is higher in overweight/obese adolescents (compared to lean adolescents) and that obesity and IR are main independent predictors of higher Ox-LDL [8, 9]. However, the association of bilirubin with Ox-LDL levels in young obese patients is unknown.
The aim of our study is to evaluate the association of Ox-LDL with clinical data of young obese patients and several biochemical variables, including total bilirubin, lipid profile, and markers of glucose metabolism.
Methods and procedures
Obese children and adolescents, aged 11–18 years, were identified from medical records, at the outpatient clinics of pediatric obesity in two hospitals in Oporto - Portugal. A group of children from 5 primary and 2 middle and high public schools from Oporto suburban setting were also recruited to this study, enlarging the obese group.
The study protocol was approved by the Committee on Ethics of both Hospitals and by the Review Committee of the Scientific Board of the Faculty of Sport of the University of Oporto.
A total of 125 obese children and adolescents participated in the study after informed and written consent of their parents. Smokers, subjects with endocrine disorders, inflammatory, hepatic or infectious diseases or under any therapy that could interfere with our results were excluded from the study.
Procedures and assays
Anthropometric characterization and clinical evaluation
All participants were subjected to clinical examination. Height and weight were measured. Obesity was defined as body mass index (BMI) z-score greater than +1.65 for age and gender, according to 2000 Centre for Disease Control and Prevention (CDC) growth charts. Waist circumference was measured at the superior border of the iliac crest, according to the protocol of the NHANES . Waist-to-Height ratio (WtHr) was calculated and a value greater than 0.5 was considered as predictor of cardiometabolic risk.
Blood was collected by venipuncture in EDTA containing tubes, after overnight fasting (10-12 h) and processed within 2 h of collection. Aliquots of plasma were made, and immediately stored at −80 °C until assayed.
Plasma levels of CRP were determined by immunoturbidimetry [CRP (latex) High-Sensitivity, Roche Diagnostics] and quantification of TB was performed by a colorimetric test (diazotized sulfanilic acid reaction, Roche Diagnostics). Adiponectin plasmatic concentration was evaluated by using a standard commercial enzyme-linked immunoassay (Bender MedSystems, San Diego, CA, USA).
The determination of circulating levels of glucose and insulin was performed by using routine automated technology (ABX Diagnostics). Homeostasis model assessment of insulin resistance (HOMAIR) was calculated according to Matthews .
The lipid profile was performed in an autoanalyzer (Cobas Mira S, Roche, Basel, Switzerland), using commercially available kits; total cholesterol and TG concentrations were determined by enzymatic colorimetric tests (cholesterol oxidase-phenol aminophenazone and glycerol-3-phosphate oxidase-phenol aminophenazone methods, Roche, resp.); HDLc was measured using enzymatic colorimetric test, after selective separation of HDLc fraction (direct HDL cholesterol, Roche); LDLc was calculated using Friedewald formula ; Ox-LDL was measured directly in plasma by using a two-site enzyme immunoassay (oxidized LDL ELISA, Mercodia, Uppsala, Sweden).
The distributions of continuous variables were analyzed using Kolmogorov-Smirnov tests to assess significant departures from Normality. Normally distributed variables are presented as mean ± SD (standard deviation) and those non-normally distributed are presented as median (interquartile range). Comparisons between two groups were performed using Student’s unpaired t-test or Mann–Whitney U test. The association between categorical variables was analyzed using chi-squared test and Fisher’s exact test.
The strength of the association between the variables was estimated by Pearson correlation coefficient, after log transformation of the variables (whenever necessary). To evaluate the contribution of the different variables (clinical and biochemical) to Ox-LDL levels, multiple regression analysis was performed, using stepwise selection, with an entry criteria of P <0.05. Age, sex, Tanner stage, weight, BMI, BMI z-score, waist circumference, WtHr, TG, total cholesterol, LDLc, HDLc, glucose, insulin, HOMAIR, CRP, adiponectin and TB were included in the analysis as independent variables.
Statistical analysis was performed using the Statistical Package for Social Sciences (SPSS, version 21.0) for Windows (SPSS Inc., Armonk, NY, USA). Statistical significance was accepted at P less than 0.05.
Anthropometric, nutritional and biochemical variables for total obese patients, and according to oxidized-LDL (Ox-LDL) levels (lower or higher than the median value observed for total patients, 39.2 U/l)
All obese patients
Ox-LDL < 39.2 U/l
Ox-LDL ≥ 39.2 U/l
(n = 125)
(n = 62)
(n = 63)
13.4 ± 1.9
13.3 ± 1.8
13.6 ± 1.9
Sex (females (%))
Tanner (pre-pub (%))
161.8 ± 9.0
162.4 ± 9.2
161.3 ± 8.8
89.3 ± 22.4
86.7 ± 18.9
91.9 ± 25.3
BMI (kg/m 2 )
33.7 ± 6.2
32.6 ± 4.9
34.9 ± 7.1
2.28 ± 0.32
2.24 ± 0.29
2.32 ± 0.36
Waist circumference (cm)
104.9 ± 15.2
103.7 ± 13.3
105.9 ± 16.8
0.65 ± 0.08
0.64 ± 0.07
0.66 ± 0.09
Total cholesterol (mg/dl)
Total bilirubin (μmol/l)
No statistically significant differences were found for metabolic health markers between clinical and school-based samples (data not shown).
In all patients, Ox-LDL correlated positively and significantly with BMI (r = 0.197, P = 0.027), total cholesterol (r = 0.630, P < 0.001), LDLc (r = 0.678, P < 0.001), TG (r = 0.510, P < 0.001), CRP (r = 0.224, P = 0.012), glucose (r = 0.196, P = 0.029), insulin (r = 0.186, P = 0.038) and HOMAIR (r = 0.209, P = 0.019); moreover, Ox-LDL levels correlated inversely and significantly with TB (r = −0.252, P = 0.005) and HDLc/Total cholesterol ratio (r = −0.732, P < 0.001).
Main determinants of oxidized LDL (Ox-LDL) levels by multiple linear regression analysis
Model 1 (R2 = 0.672)
Model 2 (R2 = 0.524)*
As far as we know this is the first study demonstrating an association of bilirubin with Ox-LDL levels in young obese patients.
The strong positive correlations found between Ox-LDL and both LDLc and TG (with their inclusion in the regression model of Ox-LDL) are in accordance with literature. Ox-LDL level is likely to increase with the number of LDL particles, and high levels of TG lead to the formation of smaller, denser LDL species, which are known to be more prone to oxidative modification [13, 14]. On the other hand, HDL is known as an atheroprotective protein and is capable to protect LDL from oxidative modification in vitro[15, 16]. This is in agreement with the inclusion of HDL in our regression analysis of Ox-LDL (β = −0.283, P < 0.001). Bilirubin was also an independent predictor of Ox-LDL levels (Table 2, model1). Actually, it was previously demonstrated that bilirubin is also capable to prevent LDL from oxidation in vitro and two recent studies involving young males  and patients under hemodialysis  reported a similar association in vivo between bilirubin and Ox-LDL, suggesting a potential cardiovascular risk protection of bilirubin.
In children, Ox-LDL was found to associate with adiposity and, independently of body fatness, with IR . Indeed, previous reports have demonstrated that obesity and IR are main independent predictors of higher Ox-LDL [8, 9]. Ryder et al. reported that IR is an independent predictor of Ox-LDL when whole-body insulin sensitivity index is used, but not when HOMAIR is used. Another study found that Ox-LDL is significantly correlated with fasting insulin and HOMAIR, but these relationships disappeared after adjustment for BMI . In the present study, BMI and markers of glucose metabolism were positively and significantly correlated with Ox-LDL levels. However, only after removing TG and HDLc from the analysis, HOMAIR was included in the regression model (model 2; Table 2). In this new model, bilirubin was the second best predictor of Ox-LDL levels (β = −0.202, P = 0.002), after LDLc (kept in the analysis) and followed by HOMAIR (β = 0.163, P = 0.010), explaining approximately 52% of the variance. It is noteworthy that none adiposity marker entered this new model, probably due to adjustment for glucose metabolism markers. Indeed, both insulin and HOMAIR were highly statistically positively correlated with adiposity markers, namely with BMI, BMI z-score, waist circumference and WtHr (r > 0.420, P < 0.001 for all).
In the present study, bilirubin was included in both regression models, independently of confounders. These results suggest that the antioxidant and/or anti-inflammatory properties of bilirubin may explain per se its negative association with Ox-LDL. However, it is also possible that other antioxidant compounds may contribute to this association. That is, bilirubin levels may actually be a marker of the status of the antioxidant system. Besides bilirubin, the endogenous antioxidant system includes other compounds that have a role in the protection from the deleterious effects of oxygen free radicals.
A previous work from our group demonstrated that young obese patients with higher body fat percentage present lower circulating levels of bilirubin and a higher degree of inflammation . In the present study we did not assess body composition (due to logistic reasons) and the only anthropometric marker related with Ox-LDL was BMI. However, it would be interesting to investigate if body fat mass in obese patients relates concomitantly with bilirubin and Ox-LDL. Furthermore, in the present study, the lack of a significant association between waist circumference and Ox-LDL levels contrast with previous findings in lean and overweight/obese adolescents  and may indicate that this is not a consistent or sensitive predictor of Ox-LDL amongst different youth populations. Small inter-individual differences in this measurement (performed at different institutions and only in obese patients) may have contributed to the absence of the association but are unlikely to explain this result by itself. Indeed, the correlation between waist circumference and Ox-LDL was far from being statistically significant. In addition, waist circumference presented well known correlations with other variables, such as insulin, HOMAIR, TG and HDLc (all statistically significant; data not shown).
Adiponectin is an important antidiabetic and antiatherogenic adipokine. In obese children and adolescents, adiponectin levels are inversely and significantly correlated with BMI, TG, and HDLc/Total cholesterol, and correlated positively and significantly with HDLc . In the present study we found no significant association between adiponectin and Ox-LDL levels, not suggesting a major role of this adipokine in regulating Ox-LDL level in young obese patients.
In conclusion, TB is independently and negatively associated with OX-LDL levels in young obese patients. Thus, patients presenting lower bilirubin may be at increased atherogenic risk and strategies to reduce Ox-LDL levels should be seek early in life.
Body mass index
High density lipoprotein
- HOMAIR :
Homeostasis model assessment of insulin resistance
Low density lipoprotein
The authors wish to thank: 1) the technician team from Porto Hospital Centre and Hospital São João for their expert assistance on blood collection, 2) all the participants that were involved in this project and their parents.
This work was funded by European Regional Development Fund (FEDER, Europe) funds through the Operational Competitiveness Programme (COMPETE, Portugal) and by National Funds through Fundação para a Ciência e a Tecnologia (FCT, Portugal) under the project FCOMP-01-0124-FEDER-028613 (PTDC/DTP-DES/0393/2012).
- Wu TW, Fung KP, Wu J, Yang CC, Weisel RD: Antioxidation of human low density lipoprotein by unconjugated and conjugated bilirubins. Biochem Pharmacol 1996, 51:859–862. 10.1016/0006-2952(95)02395-XView ArticlePubMedGoogle Scholar
- Boon AC, Hawkins CL, Bisht K, Coombes JS, Bakrania B, Wagner KH, et al.: Reduced circulating oxidized LDL is associated with hypocholesterolemia and enhanced thiol status in Gilbert syndrome. Free Radic Biol Med 2012, 52:2120–2127. 10.1016/j.freeradbiomed.2012.03.002View ArticlePubMed CentralPubMedGoogle Scholar
- Stojanov M, Stefanovic A, Dzingalasevic G, Ivanisevic J, Miljkovic M, Mandic-Radic S, et al.: Total bilirubin in young men and women: association with risk markers for cardiovascular diseases. Clin Biochem 2013, 46:1516–1519. 10.1016/j.clinbiochem.2013.06.020View ArticlePubMedGoogle Scholar
- Horsfall LJ, Nazareth I, Petersen I: Cardiovascular events as a function of serum bilirubin levels in a large, statin-treated cohort. Circulation 2012, 126:2556–2564. 10.1161/CIRCULATIONAHA.112.114066View ArticlePubMedGoogle Scholar
- McArdle PF, Whitcomb BW, Tanner K, Mitchell BD, Shuldiner AR, Parsa A: Association between bilirubin and cardiovascular disease risk factors: using Mendelian randomization to assess causal inference. BMC Cardiovasc Disord 2012, 12:16. 10.1186/1471-2261-12-16View ArticlePubMed CentralPubMedGoogle Scholar
- Lin LY, Kuo HK, Hwang JJ, Lai LP, Chiang FT, Tseng CD, et al.: Serum bilirubin is inversely associated with insulin resistance and metabolic syndrome among children and adolescents. Atherosclerosis 2009, 203:563–568. 10.1016/j.atherosclerosis.2008.07.021View ArticlePubMedGoogle Scholar
- Nascimento H, Costa E, Rocha S, Lucena C, Rocha-Pereira P, Rego C, et al.: Adiponectin and markers of metabolic syndrome in obese children and adolescents: impact of 8-mo regular physical exercise program. Pediatr Res 2014, 76:159–165. 10.1038/pr.2014.73View ArticlePubMedGoogle Scholar
- Norris AL, Steinberger J, Steffen LM, Metzig AM, Schwarzenberg SJ, Kelly AS: Circulating oxidized LDL and inflammation in extreme pediatric obesity. Obesity (Silver Spring) 2011, 19:1415–1419. 10.1038/oby.2011.21View ArticleGoogle Scholar
- Ryder JR, Vega-Lopez S, Djedjos CS, Shaibi GQ: Abdominal adiposity, insulin resistance, and oxidized low-density lipoproteins in Latino adolescents. Diabetol Metab Syndr 2013, 5:72. 10.1186/1758-5996-5-72View ArticlePubMed CentralPubMedGoogle Scholar
- National Center for Health Statistics Centers for Disease Control and Prevention Hyattsville M: ANALYTIC AND REPORTING GUIDELINES: The Third National Health and Nutrition Examination Survey, NHANES III (1988–94). Book ANALYTIC AND REPORTING GUIDELINES: The Third National Health and Nutrition Examination Survey, NHANES III (1988–94) 1996.Google Scholar
- Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC: Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 1985, 28:412–419. 10.1007/BF00280883View ArticlePubMedGoogle Scholar
- Friedewald WT, Levy RI, Fredrickson DS: Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem 1972, 18:499–502.PubMedGoogle Scholar
- Chait A, Brazg RL, Tribble DL, Krauss RM: Susceptibility of small, dense, low-density lipoproteins to oxidative modification in subjects with the atherogenic lipoprotein phenotype, pattern B. Am J Med 1993, 94:350–356. 10.1016/0002-9343(93)90144-EView ArticlePubMedGoogle Scholar
- de Graaf J, Hak-Lemmers HL, Hectors MP, Demacker PN, Hendriks JC, Stalenhoef AF: Enhanced susceptibility to in vitro oxidation of the dense low density lipoprotein subfraction in healthy subjects. Arterioscler Thromb 1991, 11:298–306. 10.1161/01.ATV.11.2.298View ArticlePubMedGoogle Scholar
- Kontush A: HDL-mediated mechanisms of protection in cardiovascular disease. Cardiovasc Res 2014, 103:341–349. 10.1093/cvr/cvu147View ArticlePubMedGoogle Scholar
- Mertens A, Holvoet P: Oxidized LDL and HDL: antagonists in atherothrombosis. Faseb J 2001, 15:2073–2084. 10.1096/fj.01-0273revView ArticlePubMedGoogle Scholar
- do Sameiro-Faria M, Kohlova M, Ribeiro S, Rocha-Pereira P, Teixeira L, Nascimento H, et al.: Potential cardiovascular risk protection of bilirubin in end-stage renal disease patients under hemodialysis. Biomed Res Int 2014, 2014:175286.PubMed CentralPubMedGoogle Scholar
- Kelly AS, Jacobs DR Jr, Sinaiko AR, Moran A, Steffen LM, Steinberger J: Relation of circulating oxidized LDL to obesity and insulin resistance in children. Pediatr Diabetes 2010, 11:552–555. 10.1111/j.1399-5448.2009.00640.xView ArticlePubMed CentralPubMedGoogle Scholar
- Belo L, Nascimento H, Kohlova M, Bronze-da-Rocha E, Fernandes J, Costa E, et al.: Body fat percentage is a major determinant of total bilirubin independently of UGT1A1*28 polymorphism in young obese. PLoS One 2014, 9:e98467. 10.1371/journal.pone.0098467View ArticlePubMed CentralPubMedGoogle Scholar
- Nascimento H, Silva L, Lourenco P, Weinfurterova R, Castro E, Rego C, et al.: Lipid profile in Portuguese obese children and adolescents: interaction of apolipoprotein E polymorphism with adiponectin levels. Arch Pediatr Adolesc Med 2009, 163:1030–1036.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.