Jorgen Engmann

Background: Despite the growing interest in the use of human genomic data for drug target identification and validation, the extent to which the spectrum of human disease has been addressed by genome-wide association studies (GWAS), or by drug development, and the degree to which these efforts overlap remain unclear.Methods: In this study we harmonize and integrate different data sources to create a sample space of all the human drug targets and diseases and identify points of convergence or divergence of GWAS and drug development efforts.Results: We show that only 612 of 11,158 diseases listed in Human Disease Ontology have an approved drug treatment in at least one region of the world. Of the 1414 diseases that are the subject of preclinical or clinical phase drug development, only 666 have been investigated in GWAS. Conversely, of the 1914 human diseases that have been the subject of GWAS, 1121 have yet to be investigated in drug development.Conclusions: We produce target-disease indication lists to help the pharmaceutical industry to prioritize future drug development efforts based on genetic evidence, academia to prioritize future GWAS for diseases without effective treatments, and both sectors to harness genetic evidence to expand the indications for licensed drugs or to identify repurposing opportunities for clinical candidates that failed in their originally intended indication.

Drug target Mendelian randomization (MR) studies use DNA sequence variants in or near a gene encoding a drug target, that alter the target’s expression or function, as a tool to anticipate the effect of drug action on the same target. Here we apply MR to prioritize drug targets for their causal relevance for coronary heart disease (CHD). The targets are further prioritized using independent replication, co-localization, protein expression profiles and data from the British National Formulary and clinicaltrials.gov. Out of the 341 drug targets identified through their association with blood lipids (HDL-C, LDL-C and triglycerides), we robustly prioritize 30 targets that might elicit beneficial effects in the prevention or treatment of CHD, including NPC1L1 and PCSK9, the targets of drugs used in CHD prevention. We discuss how this approach can be generalized to other targets, disease biomarkers and endpoints to help prioritize and validate targets during the drug development process.

Background: Nuclear magnetic resonance (NMR) spectroscopy allows triglycerides to be subclassified into 14 different classes based on particle size and lipid content. We recently showed that these subfractions have differential associations with cardiovascular disease events. Here we report the distributions and define reference interval ranges for 14 triglyceride-containing lipoprotein subfraction metabolites. Methods: Lipoprotein subfractions using the Nightingale NMR platform were measured in 9073 participants from four cohort studies contributing to the UCL-Edinburgh-Bristol consortium. The distribution of each metabolite was assessed, and reference interval ranges were calculated for a disease-free population, by sex and age group (65 years), and in a subgroup population of participants with cardiovascular disease or type 2 diabetes. We also determined the distribution across body mass index and smoking status. Results: The largest reference interval range was observed in the medium very-low density lipoprotein subclass (2.5th 97.5th percentile; 0.08 to 0.68 mmol/L). The reference intervals were comparable among male and female participants, with the exception of triglyceride in high-density lipoprotein. Triglyceride subfraction concentrations in very-low density lipoprotein, intermediate-density lipoprotein, low-density lipoprotein and high-density lipoprotein subclasses increased with increasing age and increasing body mass index. Triglyceride subfraction concentrations were significantly higher in ever smokers compared to never smokers, among those with clinical chemistry measured total triglyceride greater than 1.7 mmol/L, and in those with cardiovascular disease, and type 2 diabetes as compared to disease-free subjects. Conclusion: This is the first study to establish reference interval ranges for 14 triglyceride-containing lipoprotein subfractions in samples from the general population measured using the nuclear magnetic resonance platform. The utility of nuclear magnetic resonance lipid measures may lead to greater insights for the role of triglyceride in cardiovascular disease, emphasizing the importance of appropriate reference interval ranges for future clinical decision making.

Carotid artery intima media thickness (cIMT) and carotid plaque are measures of subclinical atherosclerosis associated with ischemic stroke and coronary heart disease (CHD). Here, we undertake meta-analyses of genome-wide association studies (GWAS) in 71,128 individuals for cIMT, and 48,434 individuals for carotid plaque traits. We identify eight novel susceptibility loci for cIMT, one independent association at the previously-identified PINX1 locus, and one novel locus for carotid plaque. Colocalization analysis with nearby vascular expression quantitative loci (cis-eQTLs) derived from arterial wall and metabolic tissues obtained from patients with CHD identifies candidate genes at two potentially additional loci, ADAMTS9 and LOXL4. LD score regression reveals significant genetic correlations between cIMT and plaque traits, and both cIMT and plaque with CHD, any stroke subtype and ischemic stroke. Our study provides insights into genes and tissue-specific regulatory mechanisms linking atherosclerosis both to its functional genomic origins and its clinical consequences in humans.

Sleep is an essential human function but its regulation is poorly understood. Identifying genetic variants associated with quality, quantity and timing of sleep will provide biological insights into the regulation of sleep and potential links with disease. Using accelerometer data from 85,670 individuals in the UK Biobank, we performed a genome-wide association study of 8 accelerometer-derived sleep traits, 5 of which are not accessible through self-report alone. We identified 47 genetic associations across the sleep traits ( P -8 ) and replicated our findings in 5,819 individuals from 3 independent studies. These included 26 novel associations for sleep quality and 10 for nocturnal sleep duration. The majority of newly identified variants were associated with a single sleep trait, except for variants previously associated with restless legs syndrome that were associated with multiple sleep traits. Of the new associated and replicated sleep duration loci, we were able to fine-map a missense variant (p.Tyr727Cys) in PDE11A , a dual-specificity 3’,5’-cyclic nucleotide phosphodiesterase expressed in the hippocampus, as the likely causal variant. As a group, sleep quality loci were enriched for serotonin processing genes and all sleep traits were enriched for cerebellar-expressed genes. These findings provide new biological insights into sleep characteristics.

The causal nature of the association between hypovitaminosis D and poor cognitive function in mid- to later-life is uncertain. Using a Mendelian randomisation(MR) approach, we examined the causal relationship between 25(OH)D and cognitive function. Data came from 172,349 participants from 17 cohorts. DHCR7(rs12785878), CYP2R1 rs12794714) and their combined synthesis score were chosen to proxy 25(OH)D. Cognitive tests were standardised into global and memory scores. Analyses were stratified by 25(OH)D tertiles, sex and age. Random effects meta-analyses assessed associations between 25(OH)D and cognitive function. Associations of serum 25(OH)D with global and memory-related cognitive function were non-linear (lower cognitive scores for both low and high 25(OH)D, p_curvature ≤ 0.006), with much of the curvature attributed to a single study. DHCR7, CYP2R1, and the synthesis score were associated with small reductions in 25(OH)D per vitamin D-decreasing allele. However, coefficients for associations with global or memory-related cognitive function were non-significant and in opposing directions for DHCR7 and CYP2R1, with no overall association observed for the synthesis score. Coefficients for the synthesis score and global and memory cognition were similar when stratified by 25(OH)D tertiles, sex and age. We found no evidence for serum 25(OH)D concentration as a causal factor for cognitive performance in mid- to later life.

Background: The implications of different adiposity measures on cardiovascular disease etiology remain unclear. In this article, we quantify and contrast causal associations of central adiposity (waist-to-hip ratio adjusted for body mass index [WHRadjBMI]) and general adiposity (body mass index [BMI]) with cardiometabolic disease. Methods: Ninety-seven independent single-nucleotide polymorphisms for BMI and 49 single-nucleotide polymorphisms for WHRadjBMI were used to conduct Mendelian randomization analyses in 14 prospective studies supplemented with coronary heart disease (CHD) data from CARDIoGRAMplusC4D (Coronary Artery Disease Genome-wide Replication and Meta-analysis [CARDIoGRAM] plus The Coronary Artery Disease [C4D] Genetics; combined total 66 842 cases), stroke from METASTROKE (12 389 ischemic stroke cases), type 2 diabetes mellitus from DIAGRAM (Diabetes Genetics Replication and Meta-analysis; 34 840 cases), and lipids from GLGC (Global Lipids Genetic Consortium; 213 500 participants) consortia. Primary outcomes were CHD, type 2 diabetes mellitus, and major stroke subtypes; secondary analyses included 18 cardiometabolic traits. Results: Each one standard deviation (SD) higher WHRadjBMI (1 SD≈0.08 U) associated with a 48% excess risk of CHD (odds ratio [OR] for CHD, 1.48; 95% confidence interval [CI], 1.28–1.71), similar to findings for BMI (1 SD≈4.6 kg/m ² ; OR for CHD, 1.36; 95% CI, 1.22–1.52). Only WHRadjBMI increased risk of ischemic stroke (OR, 1.32; 95% CI, 1.03–1.70). For type 2 diabetes mellitus, both measures had large effects: OR, 1.82 (95% CI, 1.38–2.42) and OR, 1.98 (95% CI, 1.41–2.78) per 1 SD higher WHRadjBMI and BMI, respectively. Both WHRadjBMI and BMI were associated with higher left ventricular hypertrophy, glycemic traits, interleukin 6, and circulating lipids. WHRadjBMI was also associated with higher carotid intima-media thickness (39%; 95% CI, 9%–77% per 1 SD). Conclusions: Both general and central adiposity have causal effects on CHD and type 2 diabetes mellitus. Central adiposity may have a stronger effect on stroke risk. Future estimates of the burden of adiposity on health should include measures of central and general adiposity.

The druggable genome and genome-wide association study data reveal new drug development and repurposing opportunities.

Objective We investigated discrimination and calibration of cardiovascular disease (CVD) risk scores when genotypic was added to phenotypic information. The potential of genetic information for those at intermediate risk by a phenotype-based risk score was assessed. Methods Data were from seven prospective studies including 11 851 individuals initially free of CVD or diabetes, with 1444 incident CVD events over 10 years' follow-up. We calculated a score from 53 CVD-related single nucleotide polymorphisms and an established CVD risk equation ‘QRISK-2’ comprising phenotypic measures. The area under the receiver operating characteristic curve (AUROC), detection rate for given false-positive rate (FPR) and net reclassification improvement (NRI) index were estimated for gene scores alone and in addition to the QRISK-2 CVD risk score. We also evaluated use of genetic information only for those at intermediate risk according to QRISK-2. Results The AUROC was 0.635 for QRISK-2 alone and 0.623 with addition of the gene score. The detection rate for 5% FPR improved from 11.9% to 12.0% when the gene score was added. For a 10-year CVD risk cut-off point of 10%, the NRI was 0.25% when the gene score was added to QRISK-2. Applying the genetic risk score only to those with QRISK-2 risk of 10%–Conclusion The gene score produced minimal incremental population-wide utility over phenotypic risk prediction of CVD. Tailored prediction using genetic information for those at intermediate risk may have clinical utility.

Impaired kidney function, as measured by reduced estimated glomerular filtration rate (eGFR), has been associated with increased risk of coronary heart disease (CHD) in observational studies, but it is unclear whether this association is causal or the result of confounding or reverse causation. In this study we applied Mendelian randomisation analysis using 17 genetic variants previously associated with eGFR to investigate the causal role of kidney function on CHD. We used 13,145 participants from the UCL-LSHTM-Edinburgh-Bristol (UCLEB) Consortium and 194,427 participants from the Coronary ARtery DIsease Genome-wide Replication and Meta-analysis plus Coronary Artery Disease (CARDIoGRAMplusC4D) consortium. We observed significant association of an unweighted gene score with CHD risk (odds ratio = 0.983 per additional eGFR-increasing allele, 95% CI = 0.970–0.996, p = 0.008). However, using weights calculated from UCLEB, the gene score was not associated with disease risk (p = 0.11). These conflicting results could be explained by a single SNP, rs653178, which was not associated with eGFR in the UCLEB sample, but has known pleiotropic effects that prevent us from drawing a causal conclusion. The observational association between low eGFR and increased CHD risk was not explained by potential confounders and there was no evidence of reverse causation, therefore leaving the remaining unexplained association as an open question.