A catalog of GWAS fine-mapping efforts in autoimmune disease.

Genome-wide association studies (GWASs) have enabled unbiased identification of genetic loci contributing to common complex diseases. Because GWAS loci often harbor many variants and genes, it remains a major challenge to move from GWASs' statistical associations to the identification of causal variants and genes that underlie these association signals. Researchers have applied many statistical and functional fine-mapping strategies to prioritize genetic variants and genes as potential candidates. There is no gold standard in fine-mapping approaches, but consistent results across different approaches can improve confidence in the fine-mapping findings. Here, we combined text mining with a systematic review and formed a catalog of 85 studies with evidence of fine mapping for at least one autoimmune GWAS locus. Across all fine-mapping studies, we compiled 230 GWAS loci with allelic heterogeneity estimates and predictions of causal variants and trait-relevant genes. These 230 loci included 455 combinations of locus-by-disease association signals with 15 autoimmune diseases. Using these estimates, we assessed the probability of mediating disease risk associations across genes in GWAS loci and identified robust signals of causal disease biology. We predict that this comprehensive catalog of GWAS fine-mapping efforts in autoimmune disease will greatly help distill the plethora of information in the field and inform therapeutic strategies.

Genomic atlas of the human plasma proteome.

Although plasma proteins have important roles in biological processes and are the direct targets of many drugs, the genetic factors that control inter-individual variation in plasma protein levels are not well understood. Here we characterize the genetic architecture of the human plasma proteome in healthy blood donors from the INTERVAL study. We identify 1,927 genetic associations with 1,478 proteins, a fourfold increase on existing knowledge, including trans associations for 1,104 proteins. To understand the consequences of perturbations in plasma protein levels, we apply an integrated approach that links genetic variation with biological pathway, disease, and drug databases. We show that protein quantitative trait loci overlap with gene expression quantitative trait loci, as well as with disease-associated loci, and find evidence that protein biomarkers have causal roles in disease using Mendelian randomization analysis. By linking genetic factors to diseases via specific proteins, our analyses highlight potential therapeutic targets, opportunities for matching existing drugs with new disease indications, and potential safety concerns for drugs under development.

Genetic and Structural Analysis of a SKIV2L Mutation Causing Tricho-hepato-enteric Syndrome.

BACKGROUND: Advances in genomics have facilitated the discovery of monogenic disorders in patients with unique gastro-intestinal phenotypes. Syndromic diarrhea, also called tricho-hepato-enteric (THE) syndrome, results from deleterious mutations in SKIV2L or TTC37 genes. The main features of this disorder are intractable diarrhea, abnormal hair, facial dysmorphism, immunodeficiency and liver disease. AIM: To report on a patient with THE syndrome and present the genetic analysis that facilitated diagnosis. METHODS: Whole-exome sequencing (WES) was performed in a 4-month-old female with history of congenital diarrhea and severe failure to thrive but without hair anomalies or dysmorphism. Since the parents were first-degree cousins, the analysis focused on an autosomal recessive model. Sanger sequencing was used to validate suspected variants. Mutated protein structure was modeled to assess the effect of the mutation on protein function. RESULTS: We identified an autosomal recessive C.1891G > A missense mutation (NM_006929) in SKIV2L gene that was previously described only in a compound heterozygous state as causing THE syndrome. The mutation was determined to be deleterious in multiple prediction models. Protein modeling suggested that the mutation has the potential to cause structural destabilization of SKIV2L, either through conformational changes, interference with the protein's packing, or changes at the protein's interface. CONCLUSIONS: THE syndrome can present with a broad range of clinical features in the neonatal period. WES is an important diagnostic tool in patients with congenital diarrhea and can facilitate diagnosis of various diseases presenting with atypical features.

Genetic mapping with multiple levels of phenotypic information reveals determinants of lymphocyte glucocorticoid sensitivity.

Clinical response to glucocorticoids, steroid hormones widely used as pharmaceuticals, varies extensively in that many individuals (∼30%) show a weak response to treatment. Although little is known about the molecular basis of this variation, regulatory polymorphisms are likely to play a key role given that glucocorticoids act largely through activation of a transcription factor, the glucocorticoid receptor. In an effort to characterize the molecular basis of variation in glucocorticoid sensitivity, we measured in vitro lymphocyte glucocorticoid sensitivity and transcriptome-wide response to glucocorticoids in peripheral-blood mononuclear cells from African American healthy donors. We found that variation in lymphocyte glucocorticoid sensitivity was correlated with transcriptional response at 27 genes (false-discovery rate < 0.1). Furthermore, a genome-wide association scan revealed a quantitative trait locus (QTL) for lymphocyte glucocorticoid sensitivity (rs11129354, p = 4 × 10(-8)); it was also associated with transcriptional response at multiple genes, including many (14/27) where transcriptional response was correlated with lymphocyte glucocorticoid sensitivity. Using allelic-imbalance assays, we show that this QTL is a glucocorticoid-dependent cis-regulatory polymorphism for RBMS3, which encodes an RNA-binding protein known as a tumor suppressor. We found that siRNA-mediated knockdown of RBMS3 expression increased cellular proliferation in PBMCs, consistent with the role of the gene as a negative regulator of proliferation. We propose that differences in lymphocyte glucocorticoid sensitivity reflect variation in transcriptional response, which is influenced by a glucocorticoid-dependent regulatory polymorphism that acts in cis relative to RBMS3 and in trans to affect the transcriptional response of multiple distant genes.

Characterization of genetic loci that affect susceptibility to inflammatory bowel diseases in African Americans.

BACKGROUND & AIMS: Inflammatory bowel disease (IBD) has familial aggregation in African Americans (AAs), but little is known about the molecular genetic susceptibility. Mapping studies using the Immunochip genotyping array expand the number of susceptibility loci for IBD in Caucasians to 163, but the contribution of the 163 loci and European admixture to IBD risk in AAs is unclear. We performed a genetic mapping study using the Immunochip to determine whether IBD susceptibility loci in Caucasians also affect risk in AAs and identify new associated loci. METHODS: We recruited AAs with IBD and without IBD (controls) from 34 IBD centers in the United States; additional controls were collected from 4 other Immunochip studies. Association and admixture loci were mapped for 1088 patients with Crohn's disease, 361 with ulcerative colitis, 62 with IBD type unknown, and 1797 controls; 130,241 autosomal single-nucleotide polymorphisms (SNPs) were analyzed. RESULTS: The strongest associations were observed between ulcerative colitis and HLA rs9271366 (P = 7.5 × 10(-6)), Crohn's disease and 5p13.1 rs4286721 (P = 3.5 × 10(-6)), and IBD and KAT2A rs730086 (P = 2.3 × 10(-6)). Additional suggestive associations (P < 4.2 × 10(-5)) were observed between Crohn's disease and IBD and African-specific SNPs in STAT5A and STAT3; between IBD and SNPs in IL23R, IL12B, and C2orf43; and between ulcerative colitis and SNPs near HDAC11 and near LINC00994. The latter 3 loci have not been previously associated with IBD, but require replication. Established Caucasian associations were replicated in AAs (P < 3.1 × 10(-4)) at NOD2, IL23R, 5p15.3, and IKZF3. Significant admixture (P < 3.9 × 10(-4)) was observed for 17q12-17q21.31 (IZKF3 through STAT3), 10q11.23-10q21.2, 15q22.2-15q23, and 16p12.2-16p12.1. Network analyses showed significant enrichment (false discovery rate <1 × 10(-5)) in genes that encode members of the JAK-STAT, cytokine, and chemokine signaling pathways, as well those involved in pathogenesis of measles. CONCLUSIONS: In a genetic analysis of 3308 AA IBD cases and controls, we found that many variants associated with IBD in Caucasians also showed association evidence with these diseases in AAs; we also found evidence for variants and loci not previously associated with IBD. The complex genetic factors that determine risk for or protection against IBD in different populations require further study.

Gene expression of peripheral blood cells reveals pathways downstream of glucocorticoid receptor antagonism and nab-paclitaxel treatment.

OBJECTIVES: Whereas paclitaxel treatment is associated with leukopenia, the mechanisms that underlie this effect are not well-characterized. In addition, despite the importance of glucocorticoid signaling in cancer treatment, the genomic effects of glucocorticoid receptor antagonism by mifepristone treatment in primary human cells have never been described. METHODS: As part of a randomized phase 1 clinical trial, we used microarrays to profile gene expression in peripheral blood cells sampled from each of four patients at baseline, after placebo/nanoparticle albumin-bound paclitaxel (nab-paclitaxel) treatment (cycle 1), and after mifepristone/nab-paclitaxel treatment (cycle 2). RESULTS: We found that 63 genes were differentially expressed following treatment with nab-paclitaxel, including multiple genes in the tubulin pathway. We also found 606 genes that were differentially expressed in response to mifepristone; genes downregulated by mifepristone overlapped significantly with those previously identified as being upregulated by dexamethasone. CONCLUSION: These results provide insights into the mechanisms of paclitaxel and glucocorticoid receptor inhibition in peripheral blood cells.

Interactions between glucocorticoid treatment and cis-regulatory polymorphisms contribute to cellular response phenotypes.

Glucocorticoids (GCs) mediate physiological responses to environmental stress and are commonly used as pharmaceuticals. GCs act primarily through the GC receptor (GR, a transcription factor). Despite their clear biomedical importance, little is known about the genetic architecture of variation in GC response. Here we provide an initial assessment of variability in the cellular response to GC treatment by profiling gene expression and protein secretion in 114 EBV-transformed B lymphocytes of African and European ancestry. We found that genetic variation affects the response of nearby genes and exhibits distinctive patterns of genotype-treatment interactions, with genotypic effects evident in either only GC-treated or only control-treated conditions. Using a novel statistical framework, we identified interactions that influence the expression of 26 genes known to play central roles in GC-related pathways (e.g. NQO1, AIRE, and SGK1) and that influence the secretion of IL6.

Sequential immunotherapy: towards cures for autoimmunity.

Despite major progress in the treatment of autoimmune diseases in the past two decades, most therapies do not cure disease and can be associated with increased risk of infection through broad suppression of the immune system. However, advances in understanding the causes of autoimmune disease and clinical data from novel therapeutic modalities such as chimeric antigen receptor T cell therapies provide evidence that it may be possible to re-establish immune homeostasis and, potentially, prolong remission or even cure autoimmune diseases. Here, we propose a 'sequential immunotherapy' framework for immune system modulation to help achieve this ambitious goal. This framework encompasses three steps: controlling inflammation; resetting the immune system through elimination of pathogenic immune memory cells; and promoting and maintaining immune homeostasis via immune regulatory agents and tissue repair. We discuss existing drugs and those in development for each of the three steps. We also highlight the importance of causal human biology in identifying and prioritizing novel immunotherapeutic strategies as well as informing their application in specific patient subsets, enabling precision medicine approaches that have the potential to transform clinical care.

Proteomic Associations of Adverse Outcomes in Human Heart Failure.

BACKGROUND: Identifying novel molecular drivers of disease progression in heart failure (HF) is a high-priority goal that may provide new therapeutic targets to improve patient outcomes. The authors investigated the relationship between plasma proteins and adverse outcomes in HF and their putative causal role using Mendelian randomization. METHODS AND RESULTS: The authors measured 4776 plasma proteins among 1964 participants with HF with a reduced left ventricular ejection fraction enrolled in PHFS (Penn Heart Failure Study). Assessed were the observational relationship between plasma proteins and (1) all-cause death or (2) death or HF-related hospital admission (DHFA). The authors replicated nominally significant associations in the Washington University HF registry (N=1080). Proteins significantly associated with outcomes were the subject of 2-sample Mendelian randomization and colocalization analyses. After correction for multiple testing, 243 and 126 proteins were found to be significantly associated with death and DHFA, respectively. These included small ubiquitin-like modifier 2 (standardized hazard ratio [sHR], 1.56; P<0.0001), growth differentiation factor-15 (sHR, 1.68; P<0.0001) for death, A disintegrin and metalloproteinase with thrombospondin motifs-like protein (sHR, 1.40; P<0.0001), and pulmonary-associated surfactant protein C (sHR, 1.24; P<0.0001) for DHFA. In pathway analyses, top canonical pathways associated with death and DHFA included fibrotic, inflammatory, and coagulation pathways. Genomic analyses provided evidence of nominally significant associations between levels of 6 genetically predicted proteins with DHFA and 11 genetically predicted proteins with death. CONCLUSIONS: This study implicates multiple novel proteins in HF and provides preliminary evidence of associations between genetically predicted plasma levels of 17 candidate proteins and the risk for adverse outcomes in human HF.

Endotrophin, a Collagen VI Formation-Derived Peptide, in Heart Failure.

BACKGROUND: Endotrophin, a collagen type VI-derived peptide, mediates metabolic dysregulation, inflammation, and fibrosis in animal models, but has not been studied in human heart failure (HF). METHODS: We examined the association between circulating endotrophin and outcomes in participants suffering from HF with preserved ejection fraction (HFpEF) enrolled in the TOPCAT trial (n=205). Associations were validated in a participant-level meta-analysis (n=810) that included participants with HFpEF from the PHFS study (United States; n=174), PEOPLE cohort (New Zealand; n=168), a randomized trial of vasodilator therapy (United States; n=45), a cohort from Donostia University Hospital and University of Navarra (Spain; n=171), and the TRAINING-HF trial (Spain; n=47). We also assessed associations in HF with reduced ejection fraction in PHFS (n=1,642). RESULTS: Plasma endotrophin levels at baseline were associated with risk of future death (standardized hazard ratio [HR] = 1.74; 95% confidence interval [CI]=1.36-2.24; P<0.001) and death or HF-related hospital admission (DHFA; standardized HR=2.11; 95% CI= 1.67-2.67; P<0.001) in TOPCAT. Endotrophin improved reclassification and discrimination for these outcomes beyond the MAGGIC risk score and NT-proBNP (N-terminal pro b-type natriuretic peptide). Findings were confirmed in the participant-level meta-analysis. In participants with HF with reduced ejection fraction in PHFS, endotrophin levels were associated with death (standardized HR=1.82; 95% CI=1.66-2.00; P<0.001) and DHFA (standardized HR=1.40; 95% CI=1.31-1.50; P<0.001), but the strength of the latter association was substantially lower than for the MAGGIC risk score (standardized HR=1.93; 95% CI=1.76-2.12) and BNP (standardized HR=1.78; 95% CI=1.66-1.92). CONCLUSIONS: Circulating endotrophin levels are independently associated with future poor outcomes in patients with HF, particularly in HFpEF. (Funded by Bristol Myers Squibb; Instituto de Salud Carlos III [Spain] and European Regional Development Fund; European Commission CRUCIAL project; and the U.S. National Institutes of Health National Heart, Lung, and Blood Institute.).

Proteomic Analysis of Effects of Spironolactone in Heart Failure With Preserved Ejection Fraction.

BACKGROUND: The TOPCAT trial (Treatment of Preserved Cardiac Function Heart Failure With an Aldosterone Antagonist Trial) suggested clinical benefits of spironolactone treatment among patients with heart failure with preserved ejection fraction enrolled in the Americas. However, a comprehensive assessment of biologic pathways impacted by spironolactone therapy in heart failure with preserved ejection fraction has not been performed. METHODS: We conducted aptamer-based proteomic analysis utilizing 5284 modified aptamers to 4928 unique proteins on plasma samples from TOPCAT participants from the Americas (n=164 subjects with paired samples at baseline and 1 year) to identify proteins and pathways impacted by spironolactone therapy in heart failure with preserved ejection fraction. Mean percentage change from baseline was calculated for each protein. Additionally, we conducted pathway analysis of proteins altered by spironolactone. RESULTS: Spironolactone therapy was associated with proteome-wide significant changes in 7 proteins. Among these, CARD18 (caspase recruitment domain-containing protein 18), PKD2 (polycystin 2), and PSG2 (pregnancy-specific glycoprotein 2) were upregulated, whereas HGF (hepatic growth factor), PLTP (phospholipid transfer protein), IGF2R (insulin growth factor 2 receptor), and SWP70 (switch-associated protein 70) were downregulated. CARD18, a caspase-1 inhibitor, was the most upregulated protein by spironolactone (-0.5% with placebo versus +66.5% with spironolactone, P<0.0001). The top canonical pathways that were significantly associated with spironolactone were apelin signaling, stellate cell activation, glycoprotein 6 signaling, atherosclerosis signaling, liver X receptor activation, and farnesoid X receptor activation. Among the top pathways, collagens were a consistent theme that increased in patients receiving placebo but decreased in patients randomized to spironolactone. CONCLUSIONS: Proteomic analysis in the TOPCAT trial revealed proteins and pathways altered by spironolactone, including the caspase inhibitor CARD18 and multiple pathways that involved collagens. In addition to effects on fibrosis, our studies suggest potential antiapoptotic effects of spironolactone in heart failure with preserved ejection fraction, a hypothesis that merits further exploration.

A proteome-wide genetic investigation identifies several SARS-CoV-2-exploited host targets of clinical relevance.

Background: The virus SARS-CoV-2 can exploit biological vulnerabilities (e.g. host proteins) in susceptible hosts that predispose to the development of severe COVID-19. Methods: To identify host proteins that may contribute to the risk of severe COVID-19, we undertook proteome-wide genetic colocalisation tests, and polygenic (pan) and cis-Mendelian randomisation analyses leveraging publicly available protein and COVID-19 datasets. Results: Our analytic approach identified several known targets (e.g. ABO, OAS1), but also nominated new proteins such as soluble Fas (colocalisation probability >0.9, p=1 × 10-4), implicating Fas-mediated apoptosis as a potential target for COVID-19 risk. The polygenic (pan) and cis-Mendelian randomisation analyses showed consistent associations of genetically predicted ABO protein with several COVID-19 phenotypes. The ABO signal is highly pleiotropic, and a look-up of proteins associated with the ABO signal revealed that the strongest association was with soluble CD209. We demonstrated experimentally that CD209 directly interacts with the spike protein of SARS-CoV-2, suggesting a mechanism that could explain the ABO association with COVID-19. Conclusions: Our work provides a prioritised list of host targets potentially exploited by SARS-CoV-2 and is a precursor for further research on CD209 and FAS as therapeutically tractable targets for COVID-19. Funding: MAK, JSc, JH, AB, DO, MC, EMM, MG, ID were funded by Open Targets. J.Z. and T.R.G were funded by the UK Medical Research Council Integrative Epidemiology Unit (MC_UU_00011/4). JSh and GJW were funded by the Wellcome Trust Grant 206194. This research was funded in part by the Wellcome Trust [Grant 206194]. For the purpose of open access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission.

Individuals who become infected with the virus that causes COVID-19 can experience a wide variety of symptoms. These can range from no symptoms or minor symptoms to severe illness and death. Key demographic factors, such as age, gender and race, are known to affect how susceptible an individual is to infection. However, molecular factors, such as unique gene mutations and gene expression levels can also have a major impact on patient responses by affecting the levels of proteins in the body. Proteins that are too abundant or too scarce may mean the difference between dying from or surviving COVID-19. Identifying the molecular factors in a host that affect how viruses can infect individuals, evade immune defences or trigger severe illness, could provide new ways to treat patients with COVID-19. Such factors are likely to remain constant, even when the virus mutates into new strains. Hence, insights would likely apply across all virus strains, including current strains, such as alpha and delta, and any new strains that may emerge in the future. Using such a 'natural experiment' approach, Karim et al. compared the genetic profiles of over 30,000 COVID-19 patients and a million healthy individuals. Nine proteins were found to have an impact on COVID-19 infection and disease severity. Four proteins were ranked as top priorities for potential treatment targets. One protein, called CD209 (also known as DC-SIGN), is involved in how the virus enters the host cells, and had one of the strongest associations with COVID-19. Two proteins, called IL-6R and FAS, were involved in the immune response and could be responsible for the immune over-activation often seen in severe COVID-19. Finally, one protein, called OAS1, formed part of the body's innate antiviral defence system and appeared to reduce susceptibility to COVID-19. Knowing more about the proteins that influence the severity of COVID-19 opens up new ways to predict, protect and treat patients who may have severe or fatal reactions to infection. Indeed, one of the identified proteins (IL-6R) had already been targeted in recent clinical trials with some encouraging results. Considering CD209 as a potential receptor for the virus could provide another avenue for therapeutics, similar to previously successful approaches to block the virus' known interaction with a receptor protein. Ultimately, this research could supply an entirely new set of treatment options to help combat the COVID-19 pandemic.

Quantitative Proteomic Analysis of Diabetes Mellitus in Heart Failure With Preserved Ejection Fraction.

Diabetes mellitus (DM) is associated with a higher risk of heart failure hospitalization and mortality in patients with heart failure with preserved ejection fraction (HFpEF). Using SomaScan assays and proteomics analysis of plasma from participants in the TOPCAT (Treatment of Preserved Cardiac Function Heart Failure with an Aldosterone Antagonist) trial and the Penn Heart Failure Study, this study identified 10 proteins with significantly different expression in patients with HFpEF and DM. Of these, apolipoprotein M was found to mediate 72% (95% CI: 36% to 100%; p < 0.001) of the association between DM and the risk of cardiovascular death, aborted cardiac arrest, and heart failure hospitalization.

Plasma biomarkers associated with adverse outcomes in patients with calcific aortic stenosis.

AIMS: Enhanced risk stratification of patients with aortic stenosis (AS) is necessary to identify patients at high risk for adverse outcomes, and may allow for better management of patient subgroups at high risk of myocardial damage. The objective of this study was to identify plasma biomarkers and multimarker profiles associated with adverse outcomes in AS. METHODS AND RESULTS: We studied 708 patients with calcific AS and measured 49 biomarkers using a Luminex platform. We studied the correlation between biomarkers and the risk of (i) death and (ii) death or heart failure-related hospital admission (DHFA). We also utilized machine-learning methods (a tree-based pipeline optimizer platform) to develop multimarker models associated with the risk of death and DHFA. In this cohort with a median follow-up of 2.8 years, multiple biomarkers were significantly predictive of death in analyses adjusted for clinical confounders, including tumour necrosis factor (TNF)-α [hazard ratio (HR) 1.28, P < 0.0001], TNF receptor 1 (TNFRSF1A; HR 1.38, P < 0.0001), fibroblast growth factor (FGF)-23 (HR 1.22, P < 0.0001), N-terminal pro B-type natriuretic peptide (NT-proBNP) (HR 1.58, P < 0.0001), matrix metalloproteinase-7 (HR 1.24, P = 0.0002), syndecan-1 (HR 1.27, P = 0.0002), suppression of tumorigenicity-2 (ST2) (IL1RL1; HR 1.22, P = 0.0002), interleukin (IL)-8 (CXCL8; HR 1.22, P = 0.0005), pentraxin (PTX)-3 (HR 1.17, P = 0.001), neutrophil gelatinase-associated lipocalin (LCN2; HR 1.18, P < 0.0001), osteoprotegerin (OPG) (TNFRSF11B; HR 1.26, P = 0.0002), and endostatin (COL18A1; HR 1.28, P = 0.0012). Several biomarkers were also significantly predictive of DHFA in adjusted analyses including FGF-23 (HR 1.36, P < 0.0001), TNF-α (HR 1.26, P < 0.0001), TNFR1 (HR 1.34, P < 0.0001), angiopoietin-2 (HR 1.26, P < 0.0001), syndecan-1 (HR 1.23, P = 0.0006), ST2 (HR 1.27, P < 0.0001), IL-8 (HR 1.18, P = 0.0009), PTX-3 (HR 1.18, P = 0.0002), OPG (HR 1.20, P = 0.0013), and NT-proBNP (HR 1.63, P < 0.0001). Machine-learning multimarker models were strongly associated with adverse outcomes (mean 1-year probability of death of 0%, 2%, and 60%; mean 1-year probability of DHFA of 0%, 4%, 97%; P < 0.0001). In these models, IL-6 (a biomarker of inflammation) and FGF-23 (a biomarker of calcification) emerged as the biomarkers of highest importance. CONCLUSIONS: Plasma biomarkers are strongly associated with the risk of adverse outcomes in patients with AS. Biomarkers of inflammation and calcification were most strongly related to prognosis.

Phenome-wide Mendelian randomization mapping the influence of the plasma proteome on complex diseases.

The human proteome is a major source of therapeutic targets. Recent genetic association analyses of the plasma proteome enable systematic evaluation of the causal consequences of variation in plasma protein levels. Here we estimated the effects of 1,002 proteins on 225 phenotypes using two-sample Mendelian randomization (MR) and colocalization. Of 413 associations supported by evidence from MR, 130 (31.5%) were not supported by results of colocalization analyses, suggesting that genetic confounding due to linkage disequilibrium is widespread in naïve phenome-wide association studies of proteins. Combining MR and colocalization evidence in cis-only analyses, we identified 111 putatively causal effects between 65 proteins and 52 disease-related phenotypes ( https://www.epigraphdb.org/pqtl/ ). Evaluation of data from historic drug development programs showed that target-indication pairs with MR and colocalization support were more likely to be approved, evidencing the value of this approach in identifying and prioritizing potential therapeutic targets.

Author Correction: Genome-wide mapping of plasma protein QTLs identifies putatively causal genes and pathways for cardiovascular disease.

In the originally published version of this Article, financial support was not fully acknowledged. The sentence "KS was supported by the 'Biomedical Research Program' funds at Weill Cornell Medicine in Qatar, a program funded by the Qatar Foundation" has been added to the acknowledgement section in both the PDF and HTML versions of the Article.

Genome-wide mapping of plasma protein QTLs identifies putatively causal genes and pathways for cardiovascular disease.

Identifying genetic variants associated with circulating protein concentrations (protein quantitative trait loci; pQTLs) and integrating them with variants from genome-wide association studies (GWAS) may illuminate the proteome's causal role in disease and bridge a knowledge gap regarding SNP-disease associations. We provide the results of GWAS of 71 high-value cardiovascular disease proteins in 6861 Framingham Heart Study participants and independent external replication. We report the mapping of over 16,000 pQTL variants and their functional relevance. We provide an integrated plasma protein-QTL database. Thirteen proteins harbor pQTL variants that match coronary disease-risk variants from GWAS or test causal for coronary disease by Mendelian randomization. Eight of these proteins predict new-onset cardiovascular disease events in Framingham participants. We demonstrate that identifying pQTLs, integrating them with GWAS results, employing Mendelian randomization, and prospectively testing protein-trait associations holds potential for elucidating causal genes, proteins, and pathways for cardiovascular disease and may identify targets for its prevention and treatment.

Phenome-wide association studies across large population cohorts support drug target validation.

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

HDL-cholesterol levels and risk of age-related macular degeneration: a multiethnic genetic study using Mendelian randomization.

Background: Dyslipidemia, particularly high-density lipoprotein cholesterol (HDL-C), has recently been implicated in the pathogenesis of age-related macular degeneration (AMD), the leading cause of vision loss. However, epidemiological studies have yielded conflicting results. Methods: We investigated the causal role of plasma lipid levels in AMD in multiethnic populations comprising 16 144 advanced AMD cases and 17 832 controls of European descent, together with 2219 cases and 5275 controls of Asian descent, using Mendelian randomization in three models. Model 1 is a conventional meta-analysis which does not account for pleiotropy of instrumental variable (IV) effects. Model 2 is a univariate, inverse variance weighted regression analysis that accounts for potential unbalanced pleiotropy using MR-Egger method. Finally, Model 3 is a multivariate regression analysis that addresses pleiotropy by MR-Egger method and by adjusting for effects on other lipid traits. Results: A 1 standard deviation (SD) higher HDL-cholesterol level was associated with an odds ratio (OR) for AMD of 1.17 (95% confidence interval: 1.07-1.29) in Europeans (P = 6.88 × 10-4) and of 1.58 (1.24-2.00) in Asians (P = 2.92 × 10-4) in Model 3. The corresponding OR estimates were 1.30 (1.09-1.55) in Europeans (P = 3.18 × 10-3) and 1.42 (1.11-1.80) in Asians (P = 4.42 × 10-3) in Model 1, and 1.21 (1.11-1.31) in Europeans (P = 3.12 × 10-5) and 1.51 (1.20-1.91) in Asians (P = 7.61 × 10-4) in Model 2. Conversely, neither LDL-C (Europeans: OR = 0.96, P = 0.272; Asians: OR = 1.02, P = 0.874; Model 3) nor triglyceride levels (Europeans: OR = 0.91, P = 0.102; Asians: OR = 1.06, P = 0.613) were associated with AMD. We also assessed the association between lipid levels and polypoidal choroidal vasculopathy (PCV) in Asians, a subtype of AMD, and found a similar trend for association of PCV with HDL-C levels. Conclusions: Our study shows that high levels of plasma HDL-C are causally associated with an increased risk for advanced AMD in European and Asian populations, implying that strategies reducing HDL-C levels may be useful to prevent and treat AMD.

Mapping Variation in Cellular and Transcriptional Response to 1,25-Dihydroxyvitamin D3 in Peripheral Blood Mononuclear Cells.

The active hormonal form of vitamin D, 1,25-dihydroxyvitamin D (1,25D) is an important modulator of the immune system, inhibiting cellular proliferation and regulating transcription of immune response genes. In order to characterize the genetic basis of variation in the immunomodulatory effects of 1,25D, we mapped quantitative traits of 1,25D response at both the cellular and the transcriptional level. We carried out a genome-wide association scan of percent inhibition of cell proliferation (Imax) induced by 1,25D treatment of peripheral blood mononuclear cells from 88 healthy African-American individuals. Two genome-wide significant variants were identified: rs1893662 in a gene desert on chromosome 18 (p = 2.32 x 10-8) and rs6451692 on chromosome 5 (p = 2.55 x 10-8), which may influence the anti-proliferative activity of 1,25D by regulating the expression of nearby genes such as the chemokine gene, CCL28, and the translation initiation gene, PAIP1. We also identified 8 expression quantitative trait loci at a FDR<0.10 for transcriptional response to 1,25D treatment, which include the transcriptional regulator ets variant 3-like (ETV3L) and EH-domain containing 4 (EHD4). In addition, we identified response eQTLs in vitamin D receptor binding sites near genes differentially expressed in response to 1,25D, such as FERM Domain Containing 6 (FRMD6), which plays a critical role in regulating both cell proliferation and apoptosis. Combining information from the GWAS of Imax and the response eQTL mapping enabled identification of putative Imax-associated candidate genes such as PAIP1 and the transcriptional repressor gene ZNF649. Overall, the variants identified in this study are strong candidates for immune traits and diseases linked to vitamin D, such as multiple sclerosis.