IFN Stimulated Gene Expression in the Liver is a Better Predictor of Treatment Response in Chronic Hepatitis C than the IL28B (IFN lambda 3) Genotype

BACKGROUND: Therapy of chronic hepatitis C (CHC) with pegIFNα/ribavirin achieves a sustained virologic response (SVR) in ∼55%. Pre-activation of the endogenous interferon system in the liver is associated with non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. METHODS: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. RESULTS: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. CONCLUSIONS: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

A comparison of Bayesian and frequentist approaches to incorporating external information for the prediction of prostate cancer risk.

We present the most comprehensive comparison to date of the predictive benefit of genetics in addition to currently used clinical variables, using genotype data for 33 single-nucleotide polymorphisms (SNPs) in 1,547 Caucasian men from the placebo arm of the REduction by DUtasteride of prostate Cancer Events (REDUCE®) trial. Moreover, we conducted a detailed comparison of three techniques for incorporating genetics into clinical risk prediction. The first method was a standard logistic regression model, which included separate terms for the clinical covariates and for each of the genetic markers. This approach ignores a substantial amount of external information concerning effect sizes for these Genome Wide Association Study (GWAS)-replicated SNPs. The second and third methods investigated two possible approaches to incorporating meta-analysed external SNP effect estimates - one via a weighted PCa 'risk' score based solely on the meta analysis estimates, and the other incorporating both the current and prior data via informative priors in a Bayesian logistic regression model. All methods demonstrated a slight improvement in predictive performance upon incorporation of genetics. The two methods that incorporated external information showed the greatest receiver-operating-characteristic AUCs increase from 0.61 to 0.64. The value of our methods comparison is likely to lie in observations of performance similarities, rather than difference, between three approaches of very different resource requirements. The two methods that included external information performed best, but only marginally despite substantial differences in complexity.

ADME pharmacogenetics: investigation of the pharmacokinetics of the antiretroviral agent lopinavir coformulated with ritonavir.

BACKGROUND: An ADME (absorption, distribution, metabolism and excretion)-pharmacogenetics association study may identify functional variants relevant to the pharmacokinetics of lopinavir co-formulated with ritonavir (LPV/r), a first-line anti-HIV agent. METHODS: An extensive search of literature and web resources helped select ADME genes and single nucleotide polymorphisms (SNPs, functional and HapMap tagging SNPs) with a proven or potentially relevant role in LPV/r pharmacokinetics. The study followed a two-stage design. Stage 1 (discovery) considered a Caucasian population (n=638) receiving LPV/r, where we selected 117 individuals with low LPV clearance (cases) and 90 individuals with high clearance (controls). Genotyping was performed by a 1536-SNP customized GoldenGate Illumina BeadArray. Stage 2 (confirmation) represented a replication study of candidate SNPs from the stage 1 in 148 individuals receiving LPV/r. The analysis led to formal population pharmacokinetic-pharmacogenetic modeling of demographic, environmental and candidate SNP effects. RESULTS: One thousand three hundred and eighty SNPs were successfully genotyped. Nine SNPs prioritized by the stage 1 analysis were brought to replication. Stage 2 confirmed the contribution of two functional SNPs in SLCO1B1, one functional SNP in ABCC2 and a tag SNP of the CYP3A locus in addition to body weight effect and ritonavir coadministration. According to the population pharmacokinetic-pharmacogenetic model, genetic variants explained 5% of LPV variability. Individuals homozygous rs11045819 (SLCO1B1*4) had a clearance of 12.6 l/h, compared with 5.4 l/h in the reference group, and 3.9 l/h in individuals with two or more variant alleles of rs4149056 (SLCO1B1*5), rs717620 (ABCC2) or rs6945984 (CYP3A). A subanalysis confirmed that although a significant part of the variance in LPV clearance was attributed to fluctuation in ritonavir levels, genetic variants had an additional effect on LPV clearance. CONCLUSION: The two-stage strategy successfully identified genetic variants affecting LPV/r pharmacokinetics. Such a general approach of ADME pharmacogenetics should be generalized to other drugs.

Computational analysis of the genetic and environmental contributions to disease-related human phenotypes: From predicting adverse lipid response of HIV patients receiving antiretroviral therapy to genome-wide association studies of cardiovascular and lipid related disorders

Genetic variants influence the risk to develop certain diseases or give rise to differences in drug response. Recent progresses in cost-effective, high-throughput genome-wide techniques, such as microarrays measuring Single Nucleotide Polymorphisms (SNPs), have facilitated genotyping of large clinical and population cohorts. Combining the massive genotypic data with measurements of phenotypic traits allows for the determination of genetic differences that explain, at least in part, the phenotypic variations within a population. So far, models combining the most significant variants can only explain a small fraction of the variance, indicating the limitations of current models. In particular, researchers have only begun to address the possibility of interactions between genotypes and the environment. Elucidating the contributions of such interactions is a difficult task because of the large number of genetic as well as possible environmental factors.In this thesis, I worked on several projects within this context. My first and main project was the identification of possible SNP-environment interactions, where the phenotypes were serum lipid levels of patients from the Swiss HIV Cohort Study (SHCS) treated with antiretroviral therapy. Here the genotypes consisted of a limited set of SNPs in candidate genes relevant for lipid transport and metabolism. The environmental variables were the specific combinations of drugs given to each patient over the treatment period. My work explored bioinformatic and statistical approaches to relate patients' lipid responses to these SNPs, drugs and, importantly, their interactions. The goal of this project was to improve our understanding and to explore the possibility of predicting dyslipidemia, a well-known adverse drug reaction of antiretroviral therapy. Specifically, I quantified how much of the variance in lipid profiles could be explained by the host genetic variants, the administered drugs and SNP-drug interactions and assessed the predictive power of these features on lipid responses. Using cross-validation stratified by patients, we could not validate our hypothesis that models that select a subset of SNP-drug interactions in a principled way have better predictive power than the control models using "random" subsets. Nevertheless, all models tested containing SNP and/or drug terms, exhibited significant predictive power (as compared to a random predictor) and explained a sizable proportion of variance, in the patient stratified cross-validation context. Importantly, the model containing stepwise selected SNP terms showed higher capacity to predict triglyceride levels than a model containing randomly selected SNPs. Dyslipidemia is a complex trait for which many factors remain to be discovered, thus missing from the data, and possibly explaining the limitations of our analysis. In particular, the interactions of drugs with SNPs selected from the set of candidate genes likely have small effect sizes which we were unable to detect in a sample of the present size (<800 patients).In the second part of my thesis, I performed genome-wide association studies within the Cohorte Lausannoise (CoLaus). I have been involved in several international projects to identify SNPs that are associated with various traits, such as serum calcium, body mass index, two-hour glucose levels, as well as metabolic syndrome and its components. These phenotypes are all related to major human health issues, such as cardiovascular disease. I applied statistical methods to detect new variants associated with these phenotypes, contributing to the identification of new genetic loci that may lead to new insights into the genetic basis of these traits. This kind of research will lead to a better understanding of the mechanisms underlying these pathologies, a better evaluation of disease risk, the identification of new therapeutic leads and may ultimately lead to the realization of "personalized" medicine.

Genetic association analyses reveal a role for vitamin D insufficiency in hepatitis C virus-associated hepatocellular carcinoma development

Background: V itamin D insufficiency has been associated with the occurrence of various types of cancer, but causal relationships remain elusive. Methods: Associations between t he r isk o f HCV-related HCC development and CYP2R1 , GC, and DHCR7 genotypes, which are genetic determinants of reduced 25-OH-vitamin D3 (25[OH]D3) serum levels, were determined. Results: A t otal of 5604 HCV-infected patients, 1279 with a nd 4325 without progression to HCC, w ere identified. The well-known association between 25(OH)D3 s erum levels and variations in CYP2R1 ( rs1993116, rs10741657), GC ( rs2282679), a nd DHCR7 ( rs7944926, rs12785878) g enotypes was also apparent in patients w ith chronic hepatitis C. The same genotypes of t hese single nucleotide polymorphisms (SNPs), w hich are associated with reduced 25(OH)D3 s erum levels, were significantly associated with HCV-associated HCC (P=0.07 [OR=1.13] for CYP2R1 , P=0.007 [OR=1.56] for GC, P=0.003 [OR=1.42] for DHCR7; ORs for risk genotypes). In contrast, no association between t hese genetic variations and the o utcome of antiviral therapy with pegylated interferon-α and ribavirin ( P>0.2 for e ach SNP) or liver fibrosis progression rate (P>0.2 for each SNP) was observed, s uggesting a specific influence o f the genetic d eterminants of 25(OH)D3 s erum levels o n hepatocarcinogenesis. Conclusions: Our data suggest a relatively weak but functionally relevant role for vitamin D in the prevention of HCV-related HCC development. Controlled clinical trials to assess the benefit of vitamin D supplementation in HCVinfected patients with advanced liver fibrosis or cirrhosis are warranted.

Contribution of genetic background, traditional risk factors, and HIV-related factors to coronary artery disease events in HIV-positive persons.

BACKGROUND: Persons infected with human immunodeficiency virus (HIV) have increased rates of coronary artery disease (CAD). The relative contribution of genetic background, HIV-related factors, antiretroviral medications, and traditional risk factors to CAD has not been fully evaluated in the setting of HIV infection. METHODS: In the general population, 23 common single-nucleotide polymorphisms (SNPs) were shown to be associated with CAD through genome-wide association analysis. Using the Metabochip, we genotyped 1875 HIV-positive, white individuals enrolled in 24 HIV observational studies, including 571 participants with a first CAD event during the 9-year study period and 1304 controls matched on sex and cohort. RESULTS: A genetic risk score built from 23 CAD-associated SNPs contributed significantly to CAD (P = 2.9 × 10(-4)). In the final multivariable model, participants with an unfavorable genetic background (top genetic score quartile) had a CAD odds ratio (OR) of 1.47 (95% confidence interval [CI], 1.05-2.04). This effect was similar to hypertension (OR = 1.36; 95% CI, 1.06-1.73), hypercholesterolemia (OR = 1.51; 95% CI, 1.16-1.96), diabetes (OR = 1.66; 95% CI, 1.10-2.49), ≥ 1 year lopinavir exposure (OR = 1.36; 95% CI, 1.06-1.73), and current abacavir treatment (OR = 1.56; 95% CI, 1.17-2.07). The effect of the genetic risk score was additive to the effect of nongenetic CAD risk factors, and did not change after adjustment for family history of CAD. CONCLUSIONS: In the setting of HIV infection, the effect of an unfavorable genetic background was similar to traditional CAD risk factors and certain adverse antiretroviral exposures. Genetic testing may provide prognostic information complementary to family history of CAD.

Genome-wide association study of kidney function decline in individuals of European descent.

Genome-wide association studies (GWASs) have identified multiple loci associated with cross-sectional eGFR, but a systematic genetic analysis of kidney function decline over time is missing. Here we conducted a GWAS meta-analysis among 63,558 participants of European descent, initially from 16 cohorts with serial kidney function measurements within the CKDGen Consortium, followed by independent replication among additional participants from 13 cohorts. In stage 1 GWAS meta-analysis, single-nucleotide polymorphisms (SNPs) at MEOX2, GALNT11, IL1RAP, NPPA, HPCAL1, and CDH23 showed the strongest associations for at least one trait, in addition to the known UMOD locus, which showed genome-wide significance with an annual change in eGFR. In stage 2 meta-analysis, the significant association at UMOD was replicated. Associations at GALNT11 with Rapid Decline (annual eGFR decline of 3 ml/min per 1.73 m(2) or more), and CDH23 with eGFR change among those with CKD showed significant suggestive evidence of replication. Combined stage 1 and 2 meta-analyses showed significance for UMOD, GALNT11, and CDH23. Morpholino knockdowns of galnt11 and cdh23 in zebrafish embryos each had signs of severe edema 72 h after gentamicin treatment compared with controls, but no gross morphological renal abnormalities before gentamicin administration. Thus, our results suggest a role in the deterioration of kidney function for the loci GALNT11 and CDH23, and show that the UMOD locus is significantly associated with kidney function decline.Kidney International advance online publication, 10 December 2014; doi:10.1038/ki.2014.361.

Genome-wide meta-analysis identifies six novel loci associated with habitual coffee consumption

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91 462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10-8).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee.

Mapping Bias Overestimates Reference Allele Frequencies at the HLA Genes in the 1000 Genomes Project Phase I Data.

Next-generation sequencing (NGS) technologies have become the standard for data generation in studies of population genomics, as the 1000 Genomes Project (1000G). However, these techniques are known to be problematic when applied to highly polymorphic genomic regions, such as the human leukocyte antigen (HLA) genes. Because accurate genotype calls and allele frequency estimations are crucial to population genomics analyses, it is important to assess the reliability of NGS data. Here, we evaluate the reliability of genotype calls and allele frequency estimates of the single-nucleotide polymorphisms (SNPs) reported by 1000G (phase I) at five HLA genes (HLA-A, -B, -C, -DRB1, and -DQB1). We take advantage of the availability of HLA Sanger sequencing of 930 of the 1092 1000G samples and use this as a gold standard to benchmark the 1000G data. We document that 18.6% of SNP genotype calls in HLA genes are incorrect and that allele frequencies are estimated with an error greater than ±0.1 at approximately 25% of the SNPs in HLA genes. We found a bias toward overestimation of reference allele frequency for the 1000G data, indicating mapping bias is an important cause of error in frequency estimation in this dataset. We provide a list of sites that have poor allele frequency estimates and discuss the outcomes of including those sites in different kinds of analyses. Because the HLA region is the most polymorphic in the human genome, our results provide insights into the challenges of using of NGS data at other genomic regions of high diversity.

Local effects drive heterozygosity-fitness correlations in an outcrossing long-lived tree.

Heterozygosity-fitness correlations (HFCs) have been used to understand the complex interactions between inbreeding, genetic diversity and evolution. Although frequently reported for decades, evidence for HFCs was often based on underpowered studies or inappropriate methods, and hence their underlying mechanisms are still under debate. Here, we used 6100 genome-wide single nucleotide polymorphisms (SNPs) to test for general and local effect HFCs in maritime pine (Pinus pinaster Ait.), an iconic Mediterranean forest tree. Survival was used as a fitness proxy, and HFCs were assessed at a four-site common garden under contrasting environmental conditions (total of 16 288 trees). We found no significant correlations between genome-wide heterozygosity and fitness at any location, despite variation in inbreeding explaining a substantial proportion of the total variance for survival. However, four SNPs (including two non-synonymous mutations) were involved in significant associations with survival, in particular in the common gardens with higher environmental stress, as shown by a novel heterozygosity-fitness association test at the species-wide level. Fitness effects of SNPs involved in significant HFCs were stable across maritime pine gene pools naturally growing in distinct environments. These results led us to dismiss the general effect hypothesis and suggested a significant role of heterozygosity in specific candidate genes for increasing fitness in maritime pine. Our study highlights the importance of considering the species evolutionary and demographic history and different spatial scales and testing environments when assessing and interpreting HFCs.

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance.

BackgroundBipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain.AimsWe sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL.MethodTo detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals.ResultsIntegrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85×10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54×10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals.ConclusionsOur findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Prevalence and determinants of periodic limb movements in the general population.

OBJECTIVE: Periodic limb movements during sleep (PLMS) are sleep phenomena characterized by periodic episodes of repetitive stereotyped limb movements. The aim of this study was to describe the prevalence and determinants of PLMS in a middle to older aged general population. METHODS: Data from 2,162 subjects (51.2% women, mean age = 58.4 ± 11.1 years) participating in a population-based study (HypnoLaus, Lausanne, Switzerland) were collected. Assessments included laboratory tests, sociodemographic data, personal and treatment history, and full polysomnography at home. PLMS index (PLMSI) was determined, and PLMSI > 15/h was considered as significant. RESULTS: Prevalence of PLMSI > 15/h was 28.6% (31.3% in men, 26% in women). Compared to subjects with PLMSI ≤ 15/h, subjects with PLMSI > 15/h were older (p < 0.001), were predominantly males (p = 0.007), had a higher proportion of restless legs syndrome (RLS; p < 0.001), had a higher body mass index (p = 0.001), and had a lower mean glomerular filtration rate (p < 0.001). Subjects with PLMSI > 15/h also had a higher prevalence of diabetes, hypertension, and beta-blocker or hypnotic treatments. The prevalence of antidepressant use was higher, but not statistically significant (p = 0.07). Single nucleotide polymorphisms (SNPs) within BTBD9 (rs3923809), TOX3 (rs3104788), and MEIS1 (rs2300478) genes were significantly associated with PLSMI > 15/h. Conversely, mean hemoglobin and ferritin levels were similar in both groups. In the multivariate analysis, age, male gender, antidepressant intake, RLS, and rs3923809, rs3104788, and rs2300478 SNPs were independently associated with PLMSI > 15/h. INTERPRETATION: PLMS are highly prevalent in our middle-aged European population. Age, male gender, RLS, antidepressant treatment, and specific BTBD9, TOX3, and MEIS1 SNP distribution are independent predictors of PLMSI > 15/h. ANN NEUROL 2016;79:464-474.

Host Genetics of Response to Hepatitis B Vaccine: A Systematic Review

Background. Hepatitis B virus (HBV) is an important cause of chronic viral disease worldwide and can be life threatening. While a safe and effective vaccine is widely available, 5 to 10% of healthy vaccinees fail to achieve a protective anti-hepatitis B surface antigen antibody (anti-HBs) titer (>10mIU/ml). A limited number of studies investigated host genetics of the response to HBV vaccine. To our knowledge, no comprehensive overview of genetic polymorphisms both within and outside the HLA system has been done so far. Aim. The aim of this study was to perform a systematic review of the literature of human genetics influencing immune response after hepatitis B vaccination. Methods. Literature searches using keywords were conducted in the electronic databases Medline, Embase and ISI Web of Science the cut-off date being March 2014. After selection of papers according to stringent inclusion criteria, relevant information was systematically collected from the remaining articles, including demographic data, number of patients, schedule and type of vaccine, phenotypes, genes and single nucleotide polymorphisms (SNPs) genotyping results and their association with immune response to hepatitis B vaccine. Results. The literature search produced a total of 1968 articles from which 46 studies were kept for further analyses. From these studies, data was extracted for 19 alleles from the human leukocyte antigen (HLA) region that were reported as significant at least twice. Among those alleles, 9 were firmly associated with vaccine response outcome (DQ2 [DQB1*02 and DQB1*0201], DR3 [DRB1*03 and DRB1*0301], DR7 [DRB1*07 and DRB1*0701], C4AQ0, DPB1*0401, DQ3, DQB1*06, DRB1*01 and DRB1*13 [DRB1*1301]). In addition, data was extracted for 55 different genes from which 13 extra-HLA genes had polymorphisms that were studied by different group of investigators or by the same group with a replication study. Among the 13 genes allowing comparison, 4 genes (IL-1B, IL-2, IL-4R and IL- 6) revealed no significant data, 6 genes (IL-4, IL-10, IL-12B, IL-13, TNFA, IFNG and TLR2) were explored with inconsistent results and 2 genes (CD3Z and ITGAL) yielded promising results as their association with vaccine response was confirmed by a replication approach. Furthermore, this review produced a list of 46 SNPs from 26 genes that were associated with immune response to vaccine only once, providing novel candidates to be tested in datasets from existing genome-wide association studies (GWAS). Conclusion. To the best of our knowledge, this is the first systematic review of immunogenetic studies of response to hepatitis B vaccine. While this work reassesses the role of several HLA alleles on vaccine response outcome, the associations with polymorphisms in genes outside the HLA region were rather inconsistent. Moreover, this work produced a list of 46 significant SNPs that were reported by a single group of investigators, opening up some interesting possibilities for further research.

Genome-wide association study of metabolic traits reveals novel gene-metabolite-disease links.

Metabolic traits are molecular phenotypes that can drive clinical phenotypes and may predict disease progression. Here, we report results from a metabolome- and genome-wide association study on (1)H-NMR urine metabolic profiles. The study was conducted within an untargeted approach, employing a novel method for compound identification. From our discovery cohort of 835 Caucasian individuals who participated in the CoLaus study, we identified 139 suggestively significant (P<5×10(-8)) and independent associations between single nucleotide polymorphisms (SNP) and metabolome features. Fifty-six of these associations replicated in the TasteSensomics cohort, comprising 601 individuals from São Paulo of vastly diverse ethnic background. They correspond to eleven gene-metabolite associations, six of which had been previously identified in the urine metabolome and three in the serum metabolome. Our key novel findings are the associations of two SNPs with NMR spectral signatures pointing to fucose (rs492602, P = 6.9×10(-44)) and lysine (rs8101881, P = 1.2×10(-33)), respectively. Fine-mapping of the first locus pinpointed the FUT2 gene, which encodes a fucosyltransferase enzyme and has previously been associated with Crohn's disease. This implicates fucose as a potential prognostic disease marker, for which there is already published evidence from a mouse model. The second SNP lies within the SLC7A9 gene, rare mutations of which have been linked to severe kidney damage. The replication of previous associations and our new discoveries demonstrate the potential of untargeted metabolomics GWAS to robustly identify molecular disease markers.

Validated SNPs for eGFR and their associations with albuminuria.

Albuminuria and reduced glomerular filtration rate are manifestations of chronic kidney disease (CKD) that predict end-stage renal disease, acute kidney injury, cardiovascular disease and death. We hypothesized that SNPs identified in association with the estimated glomerular filtration rate (eGFR) would also be associated with albuminuria. Within the CKDGen Consortium cohort (n= 31 580, European ancestry), we tested 16 eGFR-associated SNPs for association with the urinary albumin-to-creatinine ratio (UACR) and albuminuria [UACR >25 mg/g (women); 17 mg/g (men)]. In parallel, within the CARe Renal Consortium (n= 5569, African ancestry), we tested seven eGFR-associated SNPs for association with the UACR. We used a Bonferroni-corrected P-value of 0.003 (0.05/16) in CKDGen and 0.007 (0.05/7) in CARe. We also assessed whether the 16 eGFR SNPs were associated with the UACR in aggregate using a beta-weighted genotype score. In the CKDGen Consortium, the minor A allele of rs17319721 in the SHROOM3 gene, known to be associated with a lower eGFR, was associated with lower ln(UACR) levels (beta = -0.034, P-value = 0.0002). No additional eGFR-associated SNPs met the Bonferroni-corrected P-value threshold of 0.003 for either UACR or albuminuria. In the CARe Renal Consortium, there were no associations between SNPs and UACR with a P< 0.007. Although we found the genotype score to be associated with albuminuria (P= 0.0006), this result was driven almost entirely by the known SHROOM3 variant, rs17319721. Removal of rs17319721 resulted in a P-value 0.03, indicating a weak residual aggregate signal. No alleles, previously demonstrated to be associated with a lower eGFR, were associated with the UACR or albuminuria, suggesting that there may be distinct genetic components for these traits.