Metabolites downstream of predicted loss-of-function variants inform relationship to disease

A small minority (< 3%) of protein-coding genetic variants are predicted to lead to loss of protein function. However, these predicted loss-of-function (pLOF) variants can provide insight into mode of transcriptional effect. To examine how these changes are propagated to phenotype, we determined associations with downstream metabolites. We performed association analyses of 37 pLOF variants - previously reported to be significantly associated with disease in >400,000 subjects in UK Biobank - with metabolites. We conducted these analyses in three community-based cohorts: the Framingham Heart Study (FHS) Offspring Cohort, FHS Generation 3, and the KORA F4 cohort. We identified 19 new low-frequency or rare (minor allele frequency (MAF) <5%) pLOF variant-metabolite associations, and 12 new common (MAF > 5%) pLOF variant-metabolite associations. Rare pLOF variants in the genes BTN3A2, ENPEP, and GEM that have been associated with blood pressure in UK Biobank, were associated with vasoactive metabolites indoxyl sulfate, asymmetric dimethylarginine (ADMA), and with niacinamide, respectively. A common pLOF variant in gene CCHCR1, associated with asthma in UK Biobank, was associated with histamine and niacinamide in FHS Generation 3, both reported to play a role in this disease. Common variants in olfactory receptor gene OX4C11 that associated with blood pressure in UK Biobank were associated with the nicotine metabolite cotinine, suggesting an interaction between altered olfaction, smoking behaviour, and blood pressure. These findings provide biological validity for pLOF variant-disease associations, and point to the effector roles of common metabolites. Such an approach may provide novel disease markers and therapeutic targets.