Harnessing the Power of Electronic Health Records and Genomics for Drug Discovery.

A long-standing recognition that information from human genetics studies has the potential to accelerate drug discovery has led to decades of research on how to leverage genetic and phenotypic information for drug discovery. Established simple and advanced statistical methods that allow the simultaneous analysis of genotype and clinical phenotype data by genome- and phenome-wide analyses, colocalization analyses with quantitative trait loci data from transcriptomics and proteomics data sets from different tissues, and Mendelian randomization are essential tools for drug development in the postgenomic era. Numerous studies have demonstrated how genomic data provide opportunities for the identification of new drug targets, the repurposing of drugs, and drug safety analyses. With an increase in the number of biobanks that enable linking in-depth omics data with rich repositories of phenotypic traits via electronic health records, more powerful ways for the evaluation and validation of drug targets will continue to expand across different disciplines of clinical research.

Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

Meta-analysis of ACE inhibitor-induced angioedema identifies novel risk locus.

BACKGROUND: Angioedema is a rare but potentially life-threatening adverse drug reaction in patients receiving angiotensin-converting enzyme inhibitors (ACEis). Research suggests that susceptibility to ACEi-induced angioedema (ACEi-AE) involves both genetic and nongenetic risk factors. Genome- and exome-wide studies of ACEi-AE have identified the first genetic risk loci. However, understanding of the underlying pathophysiology remains limited. OBJECTIVE: We sought to identify further genetic factors of ACEi-AE to eventually gain a deeper understanding of its pathophysiology. METHODS: By combining data from 8 cohorts, a genome-wide association study meta-analysis was performed in more than 1000 European patients with ACEi-AE. Secondary bioinformatic analyses were conducted to fine-map associated loci, identify relevant genes and pathways, and assess the genetic overlap between ACEi-AE and other traits. Finally, an exploratory cross-ancestry analysis was performed to assess shared genetic factors in European and African-American patients with ACEi-AE. RESULTS: Three genome-wide significant risk loci were identified. One of these, located on chromosome 20q11.22, has not been implicated previously in ACEi-AE. Integrative secondary analyses highlighted previously reported genes (BDKRB2 [bradykinin receptor B2] and F5 [coagulation factor 5]) as well as biologically plausible novel candidate genes (PROCR [protein C receptor] and EDEM2 [endoplasmic reticulum degradation enhancing alpha-mannosidase like protein 2]). Lead variants at the risk loci were found with similar effect sizes and directions in an African-American cohort. CONCLUSIONS: The present results contributed to a deeper understanding of the pathophysiology of ACEi-AE by (1) providing further evidence for the involvement of bradykinin signaling and coagulation pathways and (2) suggesting, for the first time, the involvement of the fibrinolysis pathway in this adverse drug reaction. An exploratory cross-ancestry comparison implicated the relevance of the associated risk loci across diverse ancestries.

Pharmacogenomics in Stroke and Cardiovascular Disease: State of the Art.

There is considerable interindividual variability in the response to antiplatelet and anticoagulant therapies, and this variation may be attributable to genetic variants. There has been an increased understanding of the genetic architecture of stroke and cardiovascular disease, which has been driven by advancements in genomic technologies and this has raised the possibility of more targeted pharmaceutical treatments. Pharmacogenetics promises to use a patient's genetic profile to treat those who are more likely to benefit from a particular intervention by selecting the best possible therapy. Although there are numerous studies indicating strong evidence for the effect of specific genotypes on the outcomes of vascular drugs, the adoption of pharmacogenetic testing in clinical practice has been slow. This resistance may stem from sometimes conflicting findings among pharmacogenetic studies, a lack of stroke-specific randomized controlled trials to test the effectiveness of genetically-guided therapies, and the practical and cost-effective implementation of genetic testing within the clinic. Thus, this review provides an overview of the genetic variants that influence the individual responses to aspirin, clopidogrel, warfarin and statins and the different methods for pharmacogenetic testing and guidelines for clinical implementation for stroke patients.

Genome-wide Study Identifies Association between HLA-B∗55:01 and Self-Reported Penicillin Allergy.

Hypersensitivity reactions to drugs are often unpredictable and can be life threatening, underscoring a need for understanding their underlying mechanisms and risk factors. The extent to which germline genetic variation influences the risk of commonly reported drug allergies such as penicillin allergy remains largely unknown. We extracted data from the electronic health records of more than 600,000 participants from the UK, Estonian, and Vanderbilt University Medical Center's BioVU biobanks to study the role of genetic variation in the occurrence of self-reported penicillin hypersensitivity reactions. We used imputed SNP to HLA typing data from these cohorts to further fine map the human leukocyte antigen (HLA) association and replicated our results in 23andMe's research cohort involving a total of 1.12 million individuals. Genome-wide meta-analysis of penicillin allergy revealed two loci, including one located in the HLA region on chromosome 6. This signal was further fine-mapped to the HLA-B∗55:01 allele (OR 1.41 95% CI 1.33-1.49, p value 2.04 × 10-31) and confirmed by independent replication in 23andMe's research cohort (OR 1.30 95% CI 1.25-1.34, p value 1.00 × 10-47). The lead SNP was also associated with lower lymphocyte counts and in silico follow-up suggests a potential effect on T-lymphocytes at HLA-B∗55:01. We also observed a significant hit in PTPN22 and the GWAS results correlated with the genetics of rheumatoid arthritis and psoriasis. We present robust evidence for the role of an allele of the major histocompatibility complex (MHC) I gene HLA-B in the occurrence of penicillin allergy.

Genomic Insights into Myasthenia Gravis Identify Distinct Immunological Mechanisms in Early and Late Onset Disease.

OBJECTIVE: The purpose of this study was to identify disease relevant genes and explore underlying immunological mechanisms that contribute to early and late onset forms of myasthenia gravis. METHODS: We used a novel genomic methodology to integrate genomewide association study (GWAS) findings in myasthenia gravis with cell-type specific information, such as gene expression patterns and promotor-enhancer interactions, in order to identify disease-relevant genes. Subsequently, we conducted additional genomic investigations, including an expression quantitative analysis of 313 healthy people to provide mechanistic insights. RESULTS: We identified several genes that were specifically linked to early onset myasthenia gravis including TNIP1, ORMDL3, GSDMB, and TRAF3. We showed that regulators of toll-like receptor 4 signaling were enriched among these early onset disease genes (fold enrichment = 3.85, p = 6.4 × 10-3 ). In contrast, T-cell regulators CD28 and CTLA4 were exclusively linked to late onset disease. We identified 2 causal genetic variants (rs231770 and rs231735; posterior probability = 0.98 and 0.91) near the CTLA4 gene. Subsequently, we demonstrated that these causal variants result in low expression of CTLA4 (rho = -0.66, p = 1.28 × 10-38 and rho = -0.52, p = 7.01 × 10-22 , for rs231735 and rs231770, respectively). INTERPRETATION: The disease-relevant genes identified in this study are a unique resource for many disciplines, including clinicians, scientists, and the pharmaceutical industry. The distinct immunological pathways linked to early and late onset myasthenia gravis carry important implications for drug repurposing opportunities and for future studies of drug development. ANN NEUROL 2021;90:455-463.

Interferon signature in patients with STAT1 gain-of-function mutation is epigenetically determined.

STAT1 gain-of-function (GOF) variants lead to defective Th17 cell development and chronic mucocutaneous candidiasis (CMC), but frequently also to autoimmunity. Stimulation of cells with STAT1 inducing cytokines like interferons (IFN) result in hyperphosphorylation and delayed dephosphorylation of GOF STAT1. However, the mechanism how the delayed dephosphorylation exactly causes the increased expression of STAT1-dependent genes, and how the intracellular signal transduction from cytokine receptors is affected, remains unknown. In this study we show that the circulating levels of IFN-α were not persistently elevated in STAT1 GOF patients. Nevertheless, the expression of interferon signature genes was evident even in the patient with low or undetectable serum IFN-α levels. Chromatin immunoprecipitation (ChIP) experiments revealed that the active chromatin mark trimethylation of lysine 4 of histone 3 (H3K4me3), was significantly enriched in areas associated with interferon-stimulated genes in STAT1 GOF cells in comparison to cells from healthy donors. This suggests that the chromatin binding of GOF STAT1 variant promotes epigenetic changes compatible with higher gene expression and elevated reactivity to type I interferons, and possibly predisposes for interferon-related autoimmunity. The results also suggest that epigenetic rewiring may be responsible for treatment failure of Janus kinase 1/2 (JAK1/2) inhibitors in certain patients.

Translating pharmacogenomics into clinical decisions: do not let the perfect be the enemy of the good.

The field of pharmacogenomics (PGx) is gradually shifting from the reactive testing of single genes toward the proactive testing of multiple genes to improve treatment outcomes, reduce adverse events, and decrease the burden of unnecessary costs for healthcare systems. Despite the progress in the field of pharmacogenomics, its implementation into routine care has been slow due to several barriers. However, in recent years, the number of studies on the implementation of PGx has increased, all providing a wealth of knowledge on different solutions for overcoming the obstacles that have been emphasized over the past years. This review focuses on some of the challenges faced by these initiatives, the solutions and different approaches for testing that they suggest, and the evidence that they provide regarding the benefits of preemptive PGx testing.

Variation near MTNR1A associates with early development and interacts with seasons.

Melatonin is a circadian regulatory hormone with neuroprotective properties. We have previously demonstrated the association of the genetic variant rs12506228 near the melatonin receptor 1A gene (MTNR1A) with intolerance to shift-work. Furthermore, this variant has been connected to Alzheimer's disease. Because of the previously suggested role of melatonin signalling in foetal neurocognitive and sleep development, we studied here the association of rs12506228 with early development. The study sample comprised 8-month-old infants from the Finnish CHILD-SLEEP birth cohort (n = 1,301). Parental questionnaires assessed socioemotional, communication and motor development, as well as sleep length and night awakenings. The A allele of rs12506228 showed an association with slower socioemotional (p = .025) and communication (p = .0098) development, but no direct association with sleep. However, the association of the Finnish seasons with infant sleep length interacted with rs12506228. Taken together, rs12506228 near MTNR1A, which has been previously linked to adult and elderly traits, is shown here to associate with slower early cognitive development. In addition, these results suggest that the darker seasons associate with longer infant sleep time, but only in the absence of the rs12506228 AA genotype. Because the risk allele has been connected to fewer brain MT1 melatonin receptors, these associations may reflect the influence of decreased melatonin signalling in early development.

Pathogenic implications for autoimmune mechanisms derived by comparative eQTL analysis of CD4+ versus CD8+ T cells.

Inappropriate activation or inadequate regulation of CD4+ and CD8+ T cells may contribute to the initiation and progression of multiple autoimmune and inflammatory diseases. Studies on disease-associated genetic polymorphisms have highlighted the importance of biological context for many regulatory variants, which is particularly relevant in understanding the genetic regulation of the immune system and its cellular phenotypes. Here we show cell type-specific regulation of transcript levels of genes associated with several autoimmune diseases in CD4+ and CD8+ T cells including a trans-acting regulatory locus at chr12q13.2 containing the rs1131017 SNP in the RPS26 gene. Most remarkably, we identify a common missense variant in IL27, associated with type 1 diabetes that results in decreased functional activity of the protein and reduced expression levels of downstream IRF1 and STAT1 in CD4+ T cells only. Altogether, our results indicate that eQTL mapping in purified T cells provides novel functional insights into polymorphisms and pathways associated with autoimmune diseases.

Comprehensive population-based genome sequencing provides insight into hematopoietic regulatory mechanisms.

Genetic variants affecting hematopoiesis can influence commonly measured blood cell traits. To identify factors that affect hematopoiesis, we performed association studies for blood cell traits in the population-based Estonian Biobank using high-coverage whole-genome sequencing (WGS) in 2,284 samples and SNP genotyping in an additional 14,904 samples. Using up to 7,134 samples with available phenotype data, our analyses identified 17 associations across 14 blood cell traits. Integration of WGS-based fine-mapping and complementary epigenomic datasets provided evidence for causal mechanisms at several loci, including at a previously undiscovered basophil count-associated locus near the master hematopoietic transcription factor CEBPA The fine-mapped variant at this basophil count association near CEBPA overlapped an enhancer active in common myeloid progenitors and influenced its activity. In situ perturbation of this enhancer by CRISPR/Cas9 mutagenesis in hematopoietic stem and progenitor cells demonstrated that it is necessary for and specifically regulates CEBPA expression during basophil differentiation. We additionally identified basophil count-associated variation at another more pleiotropic myeloid enhancer near GATA2, highlighting regulatory mechanisms for ordered expression of master hematopoietic regulators during lineage specification. Our study illustrates how population-based genetic studies can provide key insights into poorly understood cell differentiation processes of considerable physiologic relevance.

Translating genotype data of 44,000 biobank participants into clinical pharmacogenetic recommendations: challenges and solutions.

PURPOSE: Biomedical databases combining electronic medical records and phenotypic and genomic data constitute a powerful resource for the personalization of treatment. To leverage the wealth of information provided, algorithms are required that systematically translate the contained information into treatment recommendations based on existing genotype-phenotype associations. METHODS: We developed and tested algorithms for translation of preexisting genotype data of over 44,000 participants of the Estonian biobank into pharmacogenetic recommendations. We compared the results obtained by genome sequencing, exome sequencing, and genotyping using microarrays, and evaluated the impact of pharmacogenetic reporting based on drug prescription statistics in the Nordic countries and Estonia. RESULTS: Our most striking result was that the performance of genotyping arrays is similar to that of genome sequencing, whereas exome sequencing is not suitable for pharmacogenetic predictions. Interestingly, 99.8% of all assessed individuals had a genotype associated with increased risks to at least one medication, and thereby the implementation of pharmacogenetic recommendations based on genotyping affects at least 50 daily drug doses per 1000 inhabitants. CONCLUSION: We find that microarrays are a cost-effective solution for creating preemptive pharmacogenetic reports, and with slight modifications, existing databases can be applied for automated pharmacogenetic decision support for clinicians.

Global genetic diversity of human apolipoproteins and effects on cardiovascular disease risk.

Abnormal plasma apolipoprotein levels are consistently implicated in CVD risk. Although 30% to 60% of their interindividual variability is genetic, common genetic variants explain only 10% to 20% of these differences. Rare genetic variants may be major sources of the missing heritability, yet quantitative evaluations of their contribution to phenotypic variability are lacking. Here, we analyzed whole-genome and whole-exome sequencing data from 138,632 individuals across seven major human populations to present a systematic overview of genetic apolipoprotein variability. We provide population-specific frequencies of 38 clinically important apolipoprotein alleles and identify further 6,875 genetic variants, 33% of which are novel and 98.7% of which are rare with minor allele frequencies <1%. We predicted the functional impact of rare variants and found that their relative importance differed drastically between genes and among ethnicities. Importantly, we validated the clinical relevance of multiple variants with predicted effects by leveraging association data from the CARDIoGRAM (Coronary Artery Disease Genomewide Replication and Meta-analysis) and Global Lipids Genetics consortia. Overall, we provide a consolidated overview of population-specific apolipoprotein genetics as a valuable data resource for scientists and clinicians, estimate the importance of rare genetic variants for the missing heritability of apolipoprotein-associated disease traits, and pinpoint multiple novel apolipoprotein variants with putative population-specific impacts on serum lipid levels.

DNA breaks and chromatin structural changes enhance the transcription of autoimmune regulator target genes.

The autoimmune regulator (AIRE) protein is the key factor in thymic negative selection of autoreactive T cells by promoting the ectopic expression of tissue-specific genes in the thymic medullary epithelium. Mutations in AIRE cause a monogenic autoimmune disease called autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy. AIRE has been shown to promote DNA breaks via its interaction with topoisomerase 2 (TOP2). In this study, we investigated topoisomerase-induced DNA breaks and chromatin structural alterations in conjunction with AIRE-dependent gene expression. Using RNA sequencing, we found that inhibition of TOP2 religation activity by etoposide in AIRE-expressing cells had a synergistic effect on genes with low expression levels. AIRE-mediated transcription was not only enhanced by TOP2 inhibition but also by the TOP1 inhibitor camptothecin. The transcriptional activation was associated with structural rearrangements in chromatin, notably the accumulation of γH2AX and the exchange of histone H1 with HMGB1 at AIRE target gene promoters. In addition, we found the transcriptional up-regulation to co-occur with the chromatin structural changes within the genomic cluster of carcinoembryonic antigen-like cellular adhesion molecule genes. Overall, our results suggest that the presence of AIRE can trigger molecular events leading to an altered chromatin landscape and the enhanced transcription of low-expressed genes.

Dominance genetic variation contributes little to the missing heritability for human complex traits.

For human complex traits, non-additive genetic variation has been invoked to explain "missing heritability," but its discovery is often neglected in genome-wide association studies. Here we propose a method of using SNP data to partition and estimate the proportion of phenotypic variance attributed to additive and dominance genetic variation at all SNPs (hSNP(2) and δSNP(2)) in unrelated individuals based on an orthogonal model where the estimate of hSNP(2) is independent of that of δSNP(2). With this method, we analyzed 79 quantitative traits in 6,715 unrelated European Americans. The estimate of δSNP(2) averaged across all the 79 quantitative traits was 0.03, approximately a fifth of that for additive variation (average hSNP(2) = 0.15). There were a few traits that showed substantial estimates of δSNP(2), none of which were replicated in a larger sample of 11,965 individuals. We further performed genome-wide association analyses of the 79 quantitative traits and detected SNPs with genome-wide significant dominance effects only at the ABO locus for factor VIII and von Willebrand factor. All these results suggest that dominance variation at common SNPs explains only a small fraction of phenotypic variation for human complex traits and contributes little to the missing narrow-sense heritability problem.

IgG glycosylation and DNA methylation are interconnected with smoking.

BACKGROUND: Glycosylation is one of the most common post-translation modifications with large influences on protein structure and function. The effector function of immunoglobulin G (IgG) alters between pro- and anti-inflammatory, based on its glycosylation. IgG glycan synthesis is highly complex and dynamic. METHODS: With the use of two different analytical methods for assessing IgG glycosylation, we aim to elucidate the link between DNA methylation and glycosylation of IgG by means of epigenome-wide association studies. In total, 3000 individuals from 4 cohorts were analyzed. RESULTS: The overlap of the results from the two glycan measurement panels yielded DNA methylation of 7 CpG-sites on 5 genomic locations to be associated with IgG glycosylation: cg25189904 (chr.1, GNG12); cg05951221, cg21566642 and cg01940273 (chr.2, ALPPL2); cg05575921 (chr.5, AHRR); cg06126421 (6p21.33); and cg03636183 (chr.19, F2RL3). Mediation analyses with respect to smoking revealed that the effect of smoking on IgG glycosylation may be at least partially mediated via DNA methylation levels at these 7 CpG-sites. CONCLUSION: Our results suggest the presence of an indirect link between DNA methylation and IgG glycosylation that may in part capture environmental exposures. GENERAL SIGNIFICANCE: An epigenome-wide analysis conducted in four population-based cohorts revealed an association between DNA methylation and IgG glycosylation patterns. Presumably, DNA methylation mediates the effect of smoking on IgG glycosylation.

Single base resolution analysis of 5-hydroxymethylcytosine in 188 human genes: implications for hepatic gene expression.

To improve the epigenomic analysis of tissues rich in 5-hydroxymethylcytosine (hmC), we developed a novel protocol called TAB-Methyl-SEQ, which allows for single base resolution profiling of both hmC and 5-methylcytosine by targeted next-generation sequencing. TAB-Methyl-SEQ data were extensively validated by a set of five methodologically different protocols. Importantly, these extensive cross-comparisons revealed that protocols based on Tet1-assisted bisulfite conversion provided more precise hmC values than TrueMethyl-based methods. A total of 109 454 CpG sites were analyzed by TAB-Methyl-SEQ for mC and hmC in 188 genes from 20 different adult human livers. We describe three types of variability of hepatic hmC profiles: (i) sample-specific variability at 40.8% of CpG sites analyzed, where the local hmC values correlate to the global hmC content of livers (measured by LC-MS), (ii) gene-specific variability, where hmC levels in the coding regions positively correlate to expression of the respective gene and (iii) site-specific variability, where prominent hmC peaks span only 1 to 3 neighboring CpG sites. Our data suggest that both the gene- and site-specific components of hmC variability might contribute to the epigenetic control of hepatic genes. The protocol described here should be useful for targeted DNA analysis in a variety of applications.

Cell Specific eQTL Analysis without Sorting Cells.

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn's disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus.

A meta-analysis of gene expression signatures of blood pressure and hypertension.

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p<0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%-9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension.

Genome-wide genetic homogeneity between sexes and populations for human height and body mass index.

Sex-specific genetic effects have been proposed to be an important source of variation for human complex traits. Here we use two distinct genome-wide methods to estimate the autosomal genetic correlation (rg) between men and women for human height and body mass index (BMI), using individual-level (n = ∼44 000) and summary-level (n = ∼133 000) data from genome-wide association studies. Results are consistent and show that the between-sex genetic correlation is not significantly different from unity for both traits. In contrast, we find evidence of genetic heterogeneity between sexes for waist-hip ratio (rg = ∼0.7) and between populations for BMI (rg = ∼0.9 between Europe and the USA) but not for height. The lack of evidence for substantial genetic heterogeneity for body size is consistent with empirical findings across traits and species.