Identifying therapeutic drug targets using bidirectional effect genes.

Prioritizing genes for translation to therapeutics for common diseases has been challenging. Here, we propose an approach to identify drug targets with high probability of success by focusing on genes with both gain of function (GoF) and loss of function (LoF) mutations associated with opposing effects on phenotype (Bidirectional Effect Selected Targets, BEST). We find 98 BEST genes for a variety of indications. Drugs targeting those genes are 3.8-fold more likely to be approved than non-BEST genes. We focus on five genes (IGF1R, NPPC, NPR2, FGFR3, and SHOX) with evidence for bidirectional effects on stature. Rare protein-altering variants in those genes result in significantly increased risk for idiopathic short stature (ISS) (OR = 2.75, p = 3.99 × 10-8). Finally, using functional experiments, we demonstrate that adding an exogenous CNP analog (encoded by NPPC) rescues the phenotype, thus validating its potential as a therapeutic treatment for ISS. Our results show the value of looking for bidirectional effects to identify and validate drug targets.

Exome sequencing and characterization of 49,960 individuals in the UK Biobank.

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

Interleukin-18 as a drug repositioning opportunity for inflammatory bowel disease: A Mendelian randomization study.

Support from human genetics increases the probability of success in drug development. However, few examples exist of successful genomically-driven drug repositioning. Given that a Mendelian form of severe enterocolitis is due to up-regulation of the interleukin-18 (IL18) signaling pathway, and pharmacologic inhibition of IL18 has been shown to reverse this enterocolitis, we undertook a Mendelian randomization study to test the causal effect of elevated IL18 levels on inflammatory bowel disease susceptibility (IBD) in 12,882 cases and 21,770 controls. Mendelian randomization is an established method to assess the role of biomarkers in disease etiology in a manner that minimizes confounding and prevents reverse causation. Using three SNPs that explained almost 7% of the variance in IL18 level, we found that each genetically predicted standard deviation increase in IL18 was associated with an increase in IBD susceptibility (odds ratio = 1.22, 95% CI = 1.11-1.34, P-value = 6 × 10-5). This association was further validated in 25,042 IBD cases and 34,915 controls (odds ratio = 1.13, 95% CI = 1.05-1.20). Recently, an anti-IL18 monoclonal antibody, which decreased free IL18 levels, was found to be safe, yet ineffective in a phase II trial for type 2 diabetes. Taken together, these genomic findings implicated IBD as an alternative indication for anti-IL18 therapy, which should be tested in randomized controlled trials.

Characterization of ADME gene variation in 21 populations by exome sequencing.

OBJECTIVE: Proteins involving absorption, distribution, metabolism, and excretion (ADME) play a critical role in drug pharmacokinetics. The type and frequency of genetic variation in the ADME genes differ among populations. The aim of this study was to systematically investigate common and rare ADME coding variation in diverse ethnic populations by exome sequencing. MATERIALS AND METHODS: Data derived from commercial exome capture arrays and next-generation sequencing were used to characterize coding variation in 298 ADME genes in 251 Northeast Asians and 1181 individuals from the 1000 Genomes Project. RESULTS: Approximately 75% of the ADME coding sequence was captured at high quality across the joint samples harboring more than 8000 variants, with 49% of individuals carrying at least one 'knockout' allele. ADME genes carried 50% more nonsynonymous variation than non-ADME genes (P=8.2×10) and showed significantly greater levels of population differentiation (P=7.6×10). Out of the 2135 variants identified that were predicted to be deleterious, 633 were not on commercially available ADME or general-purpose genotyping arrays. Forty deleterious variants within important ADME genes, with frequencies of at least 2% in at least one population, were identified as candidates for future pharmacogenetic studies. CONCLUSION: Exome sequencing was effective in accurately genotyping most ADME variants important for pharmacogenetic research, in addition to identifying rare or potentially de novo coding variants that may be clinically meaningful. Furthermore, as a class, ADME genes are more variable and less sensitive to purifying selection than non-ADME genes.

Precision medicine, genomics and drug discovery.

The hope for precision medicine has long been on the drug discovery horizon, well before the Human Genome Project gave it promise at the turn of the 21st century. In oncology, the concept has finally been realized and is now firmly embedded in ongoing drug discovery programs, and with many recent therapies involving some level of patient/disease stratification, including some highly personalized treatments. In addition, several drugs for rare diseases have been recently approved or are in late-stage clinical development, and new delivery modalities in cell and gene therapy and oligonucleotide approaches are yielding exciting new medicines for rare diseases of unmet need. For common complex diseases, however, the GWAS-driven advances in annotation of the genetic architecture over the past decade have not led to a concomitant shift in refined treatments. Similarly, attempts to disentangle treatment responders from non-responders via genetic predictors in pharmacogenetics studies have not met their anticipated success. It is possible that common diseases are simply lagging behind due to the inherent time lag with drug discovery, but it is also possible that their inherent multifactorial nature and their etiological and clinical heterogeneity will prove more resistant to refined treatment paradigms. The emergence of population-based resources in electronic health records, coupled with the rapid expansion of mobile devices and digital health may help to refine the measurement of phenotypic outcomes to match the exquisite detail emerging at the molecular level.

A phenome-wide association study of a lipoprotein-associated phospholipase A2 loss-of-function variant in 90 000 Chinese adults.

BACKGROUND: Lipoprotein-associated phospholipase A2 (Lp-PLA2) has been implicated in development of atherosclerosis; however, recent randomized trials of Lp-PLA2 inhibition reported no beneficial effects on vascular diseases. In East Asians, a loss-of-function variant in the PLA2G7 gene can be used to assess the effects of genetically determined lower Lp-PLA2 METHODS: PLA2G7 V279F (rs76863441) was genotyped in 91 428 individuals randomly selected from the China Kadoorie Biobank of 0.5 M participants recruited in 2004-08 from 10 regions of China, with 7 years' follow-up. Linear regression was used to assess effects of V279F on baseline traits. Logistic regression was conducted for a range of vascular and non-vascular diseases, including 41 ICD-10 coded disease categories. RESULTS: PLA2G7 V279F frequency was 5% overall (range 3-7% by region), and 9691 (11%) participants had at least one loss-of-function variant. V279F was not associated with baseline blood pressure, adiposity, blood glucose or lung function. V279F was not associated with major vascular events [7141 events; odds ratio (OR) = 0.98 per F variant, 95% confidence interval (CI) 0.90-1.06] or other vascular outcomes, including major coronary events (922 events; 0.96, 0.79-1.18) and stroke (5967 events; 1.00, 0.92-1.09). Individuals with V279F had lower risks of diabetes (7031 events; 0.91, 0.84-0.98) and asthma (182 events; 0.53, 0.28-0.98), but there was no association after adjustment for multiple testing. CONCLUSIONS: Lifelong lower Lp-PLA2 activity was not associated with major risks of vascular or non-vascular diseases in Chinese adults. Using functional genetic variants in large-scale prospective studies with linkage to a range of health outcomes is a valuable approach to inform drug development and repositioning.

Impact of genetically supported target selection on R&D productivity.

This corrects the article DOI: 10.1038/nrd.2016.164.

The support of human genetic evidence for approved drug indications.

Over a quarter of drugs that enter clinical development fail because they are ineffective. Growing insight into genes that influence human disease may affect how drug targets and indications are selected. However, there is little guidance about how much weight should be given to genetic evidence in making these key decisions. To answer this question, we investigated how well the current archive of genetic evidence predicts drug mechanisms. We found that, among well-studied indications, the proportion of drug mechanisms with direct genetic support increases significantly across the drug development pipeline, from 2.0% at the preclinical stage to 8.2% among mechanisms for approved drugs, and varies dramatically among disease areas. We estimate that selecting genetically supported targets could double the success rate in clinical development. Therefore, using the growing wealth of human genetic data to select the best targets and indications should have a measurable impact on the successful development of new drugs.

Identification of a variant in KDR associated with serum VEGFR2 and pharmacodynamics of Pazopanib.

PURPOSE: VEGF receptor (VEGFR) kinases are important drug targets in oncology that affect function of systemic endothelial cells. To discover genetic markers that affect VEGFR inhibitor pharmacodynamics, we performed a genome-wide association study of serum soluble vascular VEGFR2 concentrations [sVEGFR2], a pharmacodynamic biomarker for VEGFR2 inhibitors. EXPERIMENTAL DESIGN: We conducted a genome-wide association study (GWAS) of [sVEGFR2] in 736 healthy Old Order Amish volunteers. Gene variants identified from the GWAS were genotyped serially in a cohort of 128 patients with advanced solid tumor with baseline [sVEGFR2] measurements, and in 121 patients with renal carcinoma with [sVEGFR2] measured before and during pazopanib therapy. RESULTS: rs34231037 (C482R) in KDR, the gene encoding sVEGFR2 was found to be highly associated with [sVEGFR2], explaining 23% of the variance (P = 2.7 × 10(-37)). Association of rs34231037 with [sVEGFR2] was replicated in 128 patients with cancer with comparable effect size (P = 0.025). Furthermore, rs34231037 was a significant predictor of changes in [sVEGFR2] in response to pazopanib (P = 0.01). CONCLUSION: Our findings suggest that genome-wide analysis of phenotypes in healthy populations can expedite identification of candidate pharmacogenetic markers. Genotyping for germline variants in KDR may have clinical utility in identifying patients with cancer with unusual sensitivity to effects of VEGFR2 kinase inhibitors.

Prospective validation of HLA-DRB1*07:01 allele carriage as a predictive risk factor for lapatinib-induced liver injury.

PURPOSE: Liver injury is a serious adverse event leading to permanent discontinuation of lapatinib in affected patients. This study aimed to validate previously associated major histocompatibility complex (MHC) variants as predictors of risk of liver injury by using a large, randomized, placebo-controlled trial of lapatinib in human epidermal growth factor receptor 2-positive, early-stage breast cancer (Tykerb Evaluation After Chemotherapy [TEACH]: Lapatinib Versus Placebo In Women With Early-Stage Breast Cancer). PATIENTS AND METHODS: The frequency of ALT elevation cases was compared among four MHC variants in 1,194 patients randomly assigned to lapatinib. Cumulative ALT elevation time courses during treatment were also compared between carriers and noncarriers of specified MHC variants. RESULTS: In lapatinib-treated patients, there was a significant difference in ALT case incidence between HLA carriers and noncarriers. The highly correlated alleles HLA-DRB1*07:01 and HLA-DQA1*02:01 (study frequency, 22.4%) were associated with ALT elevation (odds ratio, 14) between cases (n = 37) and controls (n = 1,071). These associations strengthened at higher ALT elevation thresholds and in Hy's Law cases. In lapatinib-treated patients, the overall risk for National Cancer Institute-Common Terminology Criteria for Adverse Events grade 3 ALT elevation (> 5× upper limit of normal) was 2.1%; HLA allele carriers had an increased risk of 7.7%; in noncarriers, risk was reduced to 0.5%, comparable to ALT elevation for all patients receiving placebo. The increase in ALT case incidence in the lapatinib arm showed no evidence of plateau during 1 year of lapatinib treatment. CONCLUSION: These results validate HLA-DRB1*07:01 allele carriage as a predictor of increased risk of lapatinib-induced liver injury and implicate an immune pathology. The HLA association could support clinical management of patients experiencing hepatotoxicity during lapatinib treatment.

Truncated robust distance for clinical laboratory safety data monitoring and assessment.

Laboratory safety data are routinely collected in clinical studies for safety monitoring and assessment. We have developed a truncated robust multivariate outlier detection method for identifying subjects with clinically relevant abnormal laboratory measurements. The proposed method can be applied to historical clinical data to establish a multivariate decision boundary that can then be used for future clinical trial laboratory safety data monitoring and assessment. Simulations demonstrate that the proposed method has the ability to detect relevant outliers while automatically excluding irrelevant outliers. Two examples from actual clinical studies are used to illustrate the use of this method for identifying clinically relevant outliers.

An abundance of rare functional variants in 202 drug target genes sequenced in 14,002 people.

Rare genetic variants contribute to complex disease risk; however, the abundance of rare variants in human populations remains unknown. We explored this spectrum of variation by sequencing 202 genes encoding drug targets in 14,002 individuals. We find rare variants are abundant (1 every 17 bases) and geographically localized, so that even with large sample sizes, rare variant catalogs will be largely incomplete. We used the observed patterns of variation to estimate population growth parameters, the proportion of variants in a given frequency class that are putatively deleterious, and mutation rates for each gene. We conclude that because of rapid population growth and weak purifying selection, human populations harbor an abundance of rare variants, many of which are deleterious and have relevance to understanding disease risk.

Deep sequencing of the LRRK2 gene in 14,002 individuals reveals evidence of purifying selection and independent origin of the p.Arg1628Pro mutation in Europe.

Genetic variation in LRRK2 predisposes to Parkinson disease (PD), which underpins its development as a therapeutic target. Here, we aimed to identify novel genotype-phenotype associations that might support developing LRRK2 therapies for other conditions. We sequenced the 51 exons of LRRK2 in cases comprising 12 common diseases (n = 9,582), and in 4,420 population controls. We identified 739 single-nucleotide variants, 62% of which were observed in only one person, including 316 novel exonic variants. We found evidence of purifying selection for the LRRK2 gene and a trend suggesting that this is more pronounced in the central (ROC-COR-kinase) core protein domains of LRRK2 than the flanking domains. Population genetic analyses revealed that LRRK2 is not especially polymorphic or differentiated in comparison to 201 other drug target genes. Among Europeans, we identified 17 carriers (0.13%) of pathogenic LRRK2 mutations that were not significantly enriched within any disease or in those reporting a family history of PD. Analysis of pathogenic mutations within Europe reveals that the p.Arg1628Pro (c4883G>C) mutation arose independently in Europe and Asia. Taken together, these findings demonstrate how targeted deep sequencing can help to reveal fundamental characteristics of clinically important loci.

Deep resequencing unveils genetic architecture of ADIPOQ and identifies a novel low-frequency variant strongly associated with adiponectin variation.

Increased adiponectin levels have been shown to be associated with a lower risk of type 2 diabetes. To understand the relations between genetic variation at the adiponectin-encoding gene, ADIPOQ, and adiponectin levels, and subsequently its role in disease, we conducted a deep resequencing experiment of ADIPOQ in 14,002 subjects, including 12,514 Europeans, 594 African Americans, and 567 Indian Asians. We identified 296 single nucleotide polymorphisms (SNPs), including 30 amino acid changes, and carried out association analyses in a subset of 3,665 subjects from two independent studies. We confirmed multiple genome-wide association study findings and identified a novel association between a low-frequency SNP (rs17366653) and adiponectin levels (P = 2.2E-17). We show that seven SNPs exert independent effects on adiponectin levels. Together, they explained 6% of adiponectin variation in our samples. We subsequently assessed association between these SNPs and type 2 diabetes in the Genetics of Diabetes Audit and Research in Tayside Scotland (GO-DARTS) study, comprised of 5,145 case and 6,374 control subjects. No evidence of association with type 2 diabetes was found, but we were also unable to exclude the possibility of substantial effects (e.g., odds ratio 95% CI for rs7366653 [0.91-1.58]). Further investigation by large-scale and well-powered Mendelian randomization studies is warranted.