Genome-Wide Association Study (GWAS) on Bilirubin Concentrations in Subjects with Metabolic Syndrome: Sex-Specific GWAS Analysis and Gene-Diet Interactions in a Mediterranean Population
Authors: Coltell Oscar , Asensio Eva M. , Sorlí José V. , Barragán Rocio , Fernández-Carrión Rebeca , Portolés Olga , Ortega-Azorín Carolina , Martínez-LaCruz Raul , González José I. , Zanón-Moreno Vicente , Gimenez-Alba Ignacio , Fitó Montserrat , Ros Emilio , Ordovas Jose M. , Corella Dolores 
Published: January 2019

Although, for decades, increased serum bilirubin concentrations were considered a threatening sign of underlying liver disease and had been associated with neonatal jaundice, data from recent years show that bilirubin is a powerful antioxidant and suggest that slightly increased serum bilirubin concentrations are protective against oxidative stress-related diseases, such as cardiovascular diseases. Therefore, a better understanding of the gene-diet interactions in determining serum bilirubin concentrations is needed. None of the previous genome-wide association studies (GWAS) on bilirubin concentrations has been stratified by sex. Therefore, considering the increasing interest in incorporating the gender perspective into nutritional genomics, our main aim was to carry out a GWAS on total serum bilirubin concentrations in a Mediterranean population with metabolic syndrome, stratified by sex. Our secondary aim was to explore, as a pilot study, the presence of gene-diet interactions at the GWAS level. We included 430 participants (188 men and 242 women, aged 55–75 years, and with metabolic syndrome) in the PREDIMED Plus-Valencia study. Global and sex-specific GWAS were undertaken to analyze associations and gene-diet interaction on total serum bilirubin. Adherence (low and high) to the Mediterranean diet (MedDiet) was analyzed as the dietary modulator. In the GWAS, we detected more than 55 SNPs associated with serum bilirubin at p < 5 × 10−8 (GWAS level). The top-ranked were four SNPs (rs4148325 (p = 9.25 × 10−24), rs4148324 (p = 9.48 × 10−24), rs6742078 (p = 1.29 × 10−23), rs887829 (p = 1.39 × 10−23), and the rs4148324 (p = 9.48 × 10−24)) in the UGT1A1 (UDP glucuronosyltransferase family 1 member A1) gene, which replicated previous findings revealing the UGT1A1 as the major locus. In the sex-specific GWAS, the top-ranked SNPs at the GWAS level were similar in men and women (the lead SNP was the rs4148324-UGT1A1 in both men (p = 4.77 × 10−11) and women (p = 2.15 × 10−14), which shows homogeneous genetic results for the major locus. There was more sex-specific heterogeneity for other minor genes associated at the suggestive level of GWAS significance (p < 1 × 10−5). We did not detect any gene-MedDiet interaction at p < 1 × 10−5 for the major genetic locus, but we detected some gene-MedDiet interactions with other genes at p < 1 × 10−5, and even at the GWAS level for the IL17B gene (p = 3.14 × 10−8). These interaction results, however, should be interpreted with caution due to our small sample size. In conclusion, our study provides new data, with a gender perspective, on genes associated with total serum bilirubin concentrations in men and women, and suggests possible additional modulations by adherence to MedDiet.