Omnibus proteome-wide association study identifies 43 risk genes for Alzheimer disease dementia.

Transcriptome-wide association study (TWAS) tools have been applied to conduct proteome-wide association studies (PWASs) by integrating proteomics data with genome-wide association study (GWAS) summary data. The genetic effects of PWAS-identified significant genes are potentially mediated through genetically regulated protein abundance, thus informing the underlying disease mechanisms better than GWAS loci. However, existing TWAS/PWAS tools are limited by considering only one statistical model. We propose an omnibus PWAS pipeline to account for multiple statistical models and demonstrate improved performance by simulation and application studies of Alzheimer disease (AD) dementia. We employ the Aggregated Cauchy Association Test to derive omnibus PWAS (PWAS-O) p values from PWAS p values obtained by three existing tools assuming complementary statistical models-TIGAR, PrediXcan, and FUSION. Our simulation studies demonstrated improved power, with well-calibrated type I error, for PWAS-O over all three individual tools. We applied PWAS-O to studying AD dementia with reference proteomic data profiled from dorsolateral prefrontal cortex of postmortem brains from individuals of European ancestry. We identified 43 risk genes, including 5 not identified by previous studies, which are interconnected through a protein-protein interaction network that includes the well-known AD risk genes TOMM40, APOC1, and APOC2. We also validated causal genetic effects mediated through the proteome for 27 (63%) PWAS-O risk genes, providing insights into the underlying biological mechanisms of AD dementia and highlighting promising targets for therapeutic development. PWAS-O can be easily applied to studying other complex diseases.

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast.

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein coexpression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

Integration of genetic fine-mapping and multi-omics data reveals candidate effector genes for hypertension.

Genome-wide association studies of blood pressure (BP) have identified >1,000 loci, but the effector genes and biological pathways at these loci are mostly unknown. Using published association summary statistics, we conducted annotation-informed fine-mapping incorporating tissue-specific chromatin segmentation and colocalization to identify causal variants and candidate effector genes for systolic BP, diastolic BP, and pulse pressure. We observed 532 distinct signals associated with ≥2 BP traits and 84 with all three. For >20% of signals, a single variant accounted for >75% posterior probability, 65 were missense variants in known (SLC39A8, ADRB2, and DBH) and previously unreported BP candidate genes (NRIP1 and MMP14). In disease-relevant tissues, we colocalized >80 and >400 distinct signals for each BP trait with cis-eQTLs and regulatory regions from promoter capture Hi-C, respectively. Integrating mouse, human disorder, gene expression and tissue abundance data, and literature review, we provide consolidated evidence for 436 BP candidate genes for future functional validation and discover several potential drug targets.

Promises and Challenges of populational Proteomics in Health and Disease.

Advances in proteomic assay technologies have significantly increased coverage and throughput, enabling recent increases in the number of large-scale population-based proteomic studies of human plasma and serum. Improvements in multiplexed protein assays have facilitated the quantification of thousands of proteins over a large dynamic range, a key requirement for detecting the lowest-ranging, and potentially the most disease-relevant, blood-circulating proteins. In this perspective, we examine how populational proteomic datasets in conjunction with other concurrent omic measures can be leveraged to better understand the genomic and non-genomic correlates of the soluble proteome, constructing biomarker panels for disease prediction, among others. Mass spectrometry workflows are discussed as they are becoming increasingly competitive with affinity-based array platforms in terms of speed, cost, and proteome coverage due to advances in both instrumentation and workflows. Despite much success, there remain considerable challenges such as orthogonal validation and absolute quantification. We also highlight emergent challenges associated with study design, analytical considerations, and data integration as population-scale studies are run in batches and may involve longitudinal samples collated over many years. Lastly, we take a look at the future of what the nascent next-generation proteomic technologies might provide to the analysis of large sets of blood samples, as well as the difficulties in designing large-scale studies that will likely require participation from multiple and complex funding sources and where data sharing, study designs, and financing must be solved.

Connecting dementia risk loci to the CSF proteome identifies pathophysiological leads for dementia.

Genome-wide association studies have successfully identified many genetic risk loci for dementia, but exact biological mechanisms through which genetic risk factors contribute to dementia remains unclear. Integrating CSF proteomic data with dementia risk loci could reveal intermediate molecular pathways connecting genetic variance to the development of dementia. We tested to what extent effects of known dementia risk loci can be observed in CSF levels of 665 proteins [proximity extension-based (PEA) immunoassays] in a deeply-phenotyped mixed memory clinic cohort [n = 502, mean age (standard deviation, SD) = 64.1 (8.7) years, 181 female (35.4%)], including patients with Alzheimer's disease (AD, n = 213), dementia with Lewy bodies (DLB, n = 50) and frontotemporal dementia (FTD, n = 93), and controls (n = 146). Validation was assessed in independent cohorts (n = 99 PEA platform, n = 198, mass reaction monitoring-targeted mass spectroscopy and multiplex assay). We performed additional analyses stratified according to diagnostic status (AD, DLB, FTD and controls separately), to explore whether associations between CSF proteins and genetic variants were specific to disease or not. We identified four AD risk loci as protein quantitative trait loci (pQTL): CR1-CR2 (rs3818361, P = 1.65 × 10-8), ZCWPW1-PILRB (rs1476679, P = 2.73 × 10-32), CTSH-CTSH (rs3784539, P = 2.88 × 10-24) and HESX1-RETN (rs186108507, P = 8.39 × 10-8), of which the first three pQTLs showed direct replication in the independent cohorts. We identified one AD-specific association between a rare genetic variant of TREM2 and CSF IL6 levels (rs75932628, P = 3.90 × 10-7). DLB risk locus GBA showed positive trans effects on seven inter-related CSF levels in DLB patients only. No pQTLs were identified for FTD loci, either for the total sample as for analyses performed within FTD only. Protein QTL variants were involved in the immune system, highlighting the importance of this system in the pathophysiology of dementia. We further identified pQTLs in stratified analyses for AD and DLB, hinting at disease-specific pQTLs in dementia. Dissecting the contribution of risk loci to neurobiological processes aids in understanding disease mechanisms underlying dementia.

Proteomics: Progress and Promise of High-Throughput Proteomics in Chronic Kidney Disease.

Current proteomic tools permit the high-throughput analysis of the blood proteome in large cohorts, including those enriched for chronic kidney disease (CKD) or its risk factors. To date, these studies have identified numerous proteins associated with cross-sectional measures of kidney function, as well as with the longitudinal risk of CKD progression. Representative signals that have emerged from the literature include an association between levels of testican-2 and favorable kidney prognosis and an association between levels of TNFRSF1A and TNFRSF1B and worse kidney prognosis. For these and other associations, however, understanding whether the proteins play a causal role in kidney disease pathogenesis remains a fundamental challenge, especially given the strong impact that kidney function can have on blood protein levels. Prior to investing in dedicated animal models or randomized trials, methods that leverage the availability of genotyping in epidemiologic cohorts-including Mendelian randomization, colocalization analyses, and proteome-wide association studies-can add evidence for causal inference in CKD proteomics research. In addition, integration of large-scale blood proteome analyses with urine and tissue proteomics, as well as improved assessment of posttranslational protein modifications (e.g., carbamylation), represent important future directions. Taken together, these approaches seek to translate progress in large-scale proteomic profiling into the promise of improved diagnostic tools and therapeutic target identification in kidney disease.

Sex-biased genetic regulation of inflammatory proteins in the Dutch population.

BACKGROUND: Significant differences in immune responses, prevalence or susceptibility of diseases and treatment responses have been described between males and females. Despite this, sex-differentiation analysis of the genetic architecture of inflammatory proteins is largely unexplored. We performed sex-stratified meta-analysis after protein quantitative trait loci (pQTL) mapping using inflammatory biomarkers profiled using targeted proteomics (Olink inflammatory panel) of two population-based cohorts of Europeans. RESULTS: Even though, around 67% of the pQTLs demonstrated shared effect between sexes, colocalization analysis identified two loci in the males (LINC01135 and ITGAV) and three loci (CNOT10, SRD5A2, and LILRB5) in the females with evidence of sex-dependent modulation by pQTL variants. Furthermore, we identified pathways with relevant functions in the sex-biased pQTL variants. We also showed through cross-validation that the sex-specific pQTLs are linked with sex-specific phenotypic traits. CONCLUSION: Our study demonstrates the relevance of genetic sex-stratified analysis in the context of genetic dissection of protein abundances among individuals and reveals that, sex-specific pQTLs might mediate sex-linked phenotypes. Identification of sex-specific pQTLs associated with sex-biased diseases can help realize the promise of individualized treatment.

Genetic control of the human brain proteome.

We generated an online brain pQTL resource for 7,376 proteins through the analysis of genetic and proteomic data derived from post-mortem samples of the dorsolateral prefrontal cortex of 330 older adults. The identified pQTLs tend to be non-synonymous variation, are over-represented among variants associated with brain diseases, and replicate well (77%) in an independent brain dataset. Comparison to a large study of brain eQTLs revealed that about 75% of pQTLs are also eQTLs. In contrast, about 40% of eQTLs were identified as pQTLs. These results are consistent with lower pQTL mapping power and greater evolutionary constraint on protein abundance. The latter is additionally supported by observations of pQTLs with large effects' tending to be rare, deleterious, and associated with proteins that have evidence for fewer protein-protein interactions. Mediation analyses using matched transcriptomic and proteomic data provided additional evidence that pQTL effects are often, but not always, mediated by mRNA. Specifically, we identified roughly 1.6 times more mRNA-mediated pQTLs than mRNA-independent pQTLs (550 versus 341). Our pQTL resource provides insight into the functional consequences of genetic variation in the human brain and a basis for novel investigations of genetics and disease.

Identification of potential therapeutic targets for nonischemic cardiomyopathy in European ancestry: an integrated multiomics analysis.

BACKGROUND: Nonischemic cardiomyopathy (NISCM) is a clinical challenge with limited therapeutic targets. This study aims to identify promising drug targets for NISCM. METHODS: We utilized cis-pQTLs from the deCODE study, which includes data from 35,559 Icelanders, and SNPs from the FinnGen study, which includes data from 1,754 NISCM cases and 340,815 controls of Finnish ancestry. Mendelian randomization (MR) analysis was performed to estimate the causal relationship between circulating plasma protein levels and NISCM risk. Proteins with significant associations underwent false discovery rate (FDR) correction, followed by Bayesian colocalization analysis. The expression of top two proteins, LILRA5 and NELL1, was further analyzed using various NISCM datasets. Descriptions from the Human Protein Atlas (HPA) validated protein expression. The impact of environmental exposures on LILRA5 was assessed using the Comparative Toxicogenomics Database (CTD), and molecular docking identified the potential small molecule interactions. RESULTS: MR analysis identified 255 circulating plasma proteins associated with NISCM, with 16 remaining significant after FDR correction. Bayesian colocalization analysis identified LILRA5 and NELL1 as significant, with PP.H4 > 0.8. LILRA5 has a protective effect (OR = 0.758, 95% CI, 0.670-0.857) while NELL1 displays the risk effect (OR = 1.290, 95% CI, 1.199-1.387) in NISCM. Decreased LILRA5 expression was found in NISCM such as diabetic, hypertrophic, dilated, and inflammatory cardiomyopathy, while NELL1 expression increased in hypertrophic cardiomyopathy. HPA data indicated high LILRA5 expression in neutrophils, macrophages and endothelial cells within normal heart and limited NELL1 expression. Immune infiltration analysis revealed decreased neutrophil in diabetic cardiomyopathy. CTD analysis identified several small molecules that affect LILRA5 mRNA expression. Among these, Estradiol, Estradiol-3-benzoate, Gadodiamide, Topotecan, and Testosterone were found to stably bind to the LILRA5 protein at the conserved VAL-15 or THR-133 residues in the Ig-like C2 domain. CONCLUSION: Based on European Ancestry Cohort, this study reveals that LILRA5 and NELL1 are potential therapeutic targets for NISCM, with LILRA5 showing particularly promising prospects in diabetic cardiomyopathy. Several small molecules interact with LILRA5, implying potential clinical implication.

Exploration of protein and genetic targets causing atrioventricular block: mendelian-randomization analyses based on eQTL data and pQTL data.

BACKGROUND: Atrioventricular block (AVB) is a heterogeneous group of arrhythmias. AVB can lead to sudden arrest of the heart and subsequent syncope or sudden cardiac death. Few scholars have investigated the underlying molecular mechanisms of AVB. Finding molecular markers can facilitate understanding of AVB and exploration of therapeutic targets. METHODS: Two-sample Mendelian randomization (MR) analysis was undertaken with inverse variance weighted (IVW) model and Wald ratio as the primary approach. Reverse MR analysis was undertaken to identify the associated protein targets and gene targets. Expression quantitative trait loci (eQTL) data from the eQTLGen database and protein quantitative trait loci (pQTL) data from three previous large-scale proteomic studies on plasma were retrieved as exposure data. Genome-wide association study (GWAS) summary data (586 cases and 379,215 controls) for AVB were retrieved from the UK Biobank database. Colocalization analyses were undertaken to identify the effect of filtered markers on outcome data. Databases (DrugBank, Therapeutic Target, PubChem) were used to identify drugs that interacted with targets. RESULTS: We discovered that 692 genes and 42 proteins showed a significant correlation with the AVB phenotype. Proteins (cadherin-5, sTie-1, Notch 1) and genes (DNAJC30, ABO) were putative molecules to AVB. Drug-interaction analyses revealed anticancer drugs such as tyrosine-kinase inhibitors and TIMD3 inhibitors could cause AVB. Other substances (e.g. toxins, neurological drugs) could also cause AVB. CONCLUSIONS: We identified the proteins (cadherin-5, sTie-1, Notch 1) and gene (DNAJC30, ABO) targets associated with AVB pathogenesis. Anticancer drugs (tyrosine-kinase inhibitors, TIMD3 inhibitors), toxins, or neurological drugs could also cause AVB.

Blood DNA methylomic signatures associated with CSF biomarkers of Alzheimer's disease in the EMIF-AD study.

INTRODUCTION: We investigated blood DNA methylation patterns associated with 15 well-established cerebrospinal fluid (CSF) biomarkers of Alzheimer's disease (AD) pathophysiology, neuroinflammation, and neurodegeneration. METHODS: We assessed DNA methylation in 885 blood samples from the European Medical Information Framework for Alzheimer's Disease (EMIF-AD) study using the EPIC array. RESULTS: We identified Bonferroni-significant differential methylation associated with CSF YKL-40 (five loci) and neurofilament light chain (NfL; seven loci) levels, with two of the loci associated with CSF YKL-40 levels correlating with plasma YKL-40 levels. A co-localization analysis showed shared genetic variants underlying YKL-40 DNA methylation and CSF protein levels, with evidence that DNA methylation mediates the association between genotype and protein levels. Weighted gene correlation network analysis identified two modules of co-methylated loci correlated with several amyloid measures and enriched in pathways associated with lipoproteins and development. DISCUSSION: We conducted the most comprehensive epigenome-wide association study (EWAS) of AD-relevant CSF biomarkers to date. Future work should explore the relationship between YKL-40 genotype, DNA methylation, and protein levels in the brain. HIGHLIGHTS: Blood DNA methylation was assessed in the EMIF-AD MBD study. Epigenome-wide association studies (EWASs) were performed for 15 Alzheimer's disease (AD)-relevant cerebrospinal fluid (CSF) biomarker measures. Five Bonferroni-significant loci were associated with YKL-40 levels and seven with neurofilament light chain (NfL). DNA methylation in YKL-40 co-localized with previously reported genetic variation. DNA methylation potentially mediates the effect of single-nucleotide polymorphisms (SNPs) in YKL-40 on CSF protein levels.

An optimal variant to gene distance window derived from an empirical definition of cis and trans protein QTLs.

BACKGROUND: A genome-wide association study (GWAS) correlates variation in the genotype with variation in the phenotype across a cohort, but the causal gene mediating that impact is often unclear. When the phenotype is protein abundance, a reasonable hypothesis is that the gene encoding that protein is the causal gene. However, as variants impacting protein levels can occur thousands or even millions of base pairs from the gene encoding the protein, it is unclear at what distance this simple hypothesis breaks down. RESULTS: By making the simple assumption that cis-pQTLs should be distance dependent while trans-pQTLs are distance independent, we arrive at a simple and empirical distance cutoff separating cis- and trans-pQTLs. Analyzing a recent large-scale pQTL study (Pietzner in Science 374:eabj1541, 2021) we arrive at an estimated distance cutoff of 944 kilobasepairs (95% confidence interval: 767-1,161) separating the cis and trans regimes. CONCLUSIONS: We demonstrate that this simple model can be applied to other molecular GWAS traits. Since much of biology is built on molecular traits like protein, transcript and metabolite abundance, we posit that the mathematical models for cis and trans distance distributions derived here will also apply to more complex phenotypes and traits.

Identification of 969 protein quantitative trait loci in an African American population with kidney disease attributed to hypertension.

Investigations into the causal underpinnings of disease processes can be aided by the incorporation of genetic information. Genetic studies require populations varied in both ancestry and prevalent disease in order to optimize discovery and ensure generalizability of findings to the global population. Here, we report the genetic determinants of the serum proteome in 466 African Americans with chronic kidney disease attributed to hypertension from the richly phenotyped African American Study of Kidney Disease and Hypertension (AASK) study. Using the largest aptamer-based protein profiling platform to date (6,790 proteins or protein complexes), we identified 969 genetic associations with 900 unique proteins; including 52 novel cis (local) associations and 379 novel trans (distant) associations. The genetic effects of previously published cis-protein quantitative trait loci (pQTLs) were found to be highly reproducible, and we found evidence that our novel genetic signals colocalize with gene expression and disease processes. Many trans- pQTLs were found to reflect associations mediated by the circulating cis protein, and the common trans-pQTLs are enriched for processes involving extracellular vesicles, highlighting a plausible mechanism for distal regulation of the levels of secreted proteins. Thus, our study generates a valuable resource of genetic associations linking variants to protein levels and disease in an understudied patient population to inform future studies of drug targets and physiology.

Characterization of Proteome Variation During Modern Maize Breeding.

The success of modern maize breeding has been demonstrated by remarkable increases in productivity with tremendous modification of agricultural phenotypes over the last century. Although the underlying genetic changes of the maize adaptation from tropical to temperate regions have been extensively studied, our knowledge is limited regarding the accordance of protein and mRNA expression levels accompanying such adaptation. Here we conducted an integrative analysis of proteomic and transcriptomic changes in a maize association panel. The minimum extent of correlation between protein and RNA levels suggests that variation in mRNA expression is often not indicative of protein expression at a population scale. This is corroborated by the observation that mRNA- and protein-based coexpression networks are relatively independent of each other, and many pQTLs arise without the presence of corresponding eQTLs. Importantly, compared with transcriptome, the subtypes categorized by the proteome show a markedly high accuracy to resemble the genomic subpopulation. These findings suggest that proteome evolved under a greater evolutionary constraint than transcriptome during maize adaptation from tropical to temperate regions. Overall, the integrated multi-omics analysis provides a functional context to interpret gene expression variation during modern maize breeding.

Identification of drought responsive proteins and related proteomic QTLs in barley.

Drought is a major abiotic stress that negatively influences crop yield. Breeding strategies for improved drought resistance require an improved knowledge of plant drought responses. We therefore applied drought to barley recombinant inbred lines and their parental genotypes shortly before tillering. A large-scale proteomic analysis of leaf and root tissue revealed proteins that respond to drought in a genotype-specific manner. Of these, Rubisco activase in chloroplast, luminal binding protein in endoplasmic reticulum, phosphoglycerate mutase, glutathione S-transferase, heat shock proteins and enzymes involved in phenylpropanoid biosynthesis showed strong genotype×environment interactions. These data were subjected to genetic linkage analysis and the identification of proteomic QTLs that have potential value in marker-assisted breeding programs.

Extensive co-regulation of neighboring genes complicates the use of eQTLs in target gene prioritization.

Identifying causal genes underlying genome-wide association studies (GWASs) is a fundamental problem in human genetics. Although colocalization with gene expression quantitative trait loci (eQTLs) is often used to prioritize GWAS target genes, systematic benchmarking has been limited due to unavailability of large ground truth datasets. Here, we re-analyzed plasma protein QTL data from 3,301 individuals of the INTERVAL cohort together with 131 eQTL Catalog datasets. Focusing on variants located within or close to the affected protein identified 793 proteins with at least one cis-pQTL where we could assume that the most likely causal gene was the gene coding for the protein. We then benchmarked the ability of cis-eQTLs to recover these causal genes by comparing three Bayesian colocalization methods (coloc.susie, coloc.abf, and CLPP) and five Mendelian randomization (MR) approaches (three varieties of inverse-variance weighted MR, MR-RAPS, and MRLocus). We found that assigning fine-mapped pQTLs to their closest protein coding genes outperformed all colocalization methods regarding both precision (71.9%) and recall (76.9%). Furthermore, the colocalization method with the highest recall (coloc.susie - 46.3%) also had the lowest precision (45.1%). Combining evidence from multiple conditionally distinct colocalizing QTLs with MR increased precision to 81%, but this was accompanied by a large reduction in recall to 7.1%. Furthermore, the choice of the MR method greatly affected performance, with the standard inverse-variance-weighted MR often producing many false positives. Our results highlight that linking GWAS variants to target genes remains challenging with eQTL evidence alone, and prioritizing novel targets requires triangulation of evidence from multiple sources.

Identification of novel pQTL-SNPs associated with lung adenocarcinoma risk: A multi-stage study.

BACKGROUND AND OBJECTIVE: To explore the association between protein quantitative trait loci (pQTL-SNPs) and the risk of LUAD. METHODS: "Blood +" high depth blood proteomics analysis was performed on plasma from female LUAD patients and female healthy controls, and combined with proteomics data from tumors and adjacent non-tumor tissues of female LUAD patients to screen proteins uniformly expressed in plasma and tissues. pQTL-SNPs were then screened through multiple databases and subjected to multilevel screening. The associations between selected pQTL-SNPs and LUAD risk were evaluated by Female Lung Cancer Consortium in Asia GWAS (FLCCA GWAS). Enzyme linked immunosorbent assay (ELISA) is used to determine the levels of candidate protein. RESULTS: A total of 7 pQTL-SNPs were significantly associated with altered LUAD risk (p < 0.05). Meanwhile, the expression of their corresponding target proteins were all decreased in both plasma and tumor tissues of LUAD cases, which may play a role of tumor suppressor proteins. After mutation of 3 pQTL-SNPs (rs7683000, rs73224660, and rs2776937), the expression of corresponding target proteins BST1 and NRP1 decreased, and as potential tumor suppressor proteins, which may promote tumorigenesis and further increasing the risk of developing LUAD (OR >1, p < 0.05); while after mutation the other pQTL-SNP rs62069916, the corresponding target protein APOH expression was increased, while as a potential tumor suppressor protein, which may inhibit tumorigenesis and further reduced the risk of developing LUAD (OR <1, p < 0.05). In addition, the expression of NRP1 and APOH were significant decreased in LUAD cell lines and validated in plasma of LUAD patients. CONCLUSION: A total of 4 pQTL-SNPs (rs7683000, rs73224660, rs2776937, and rs62069916) may associate with altered LUAD risk by regulating the expression of target proteins (BST1, NRP1, and APOH) after mutation.

Polygenic scores and Mendelian randomization identify plasma proteins causally implicated in Alzheimer's disease.

Background: An increasing body of evidence suggests that neuroinflammation is one of the key drivers of late-onset Alzheimer's disease (LOAD) pathology. Due to the increased permeability of the blood-brain barrier (BBB) in older adults, peripheral plasma proteins can infiltrate the central nervous system (CNS) and drive neuroinflammation through interactions with neurons and glial cells. Because these inflammatory factors are heritable, a greater understanding of their genetic relationship with LOAD could identify new biomarkers that contribute to LOAD pathology or offer protection against it. Methods: We used a genome-wide association study (GWAS) of 90 different plasma proteins (n = 17,747) to create polygenic scores (PGSs) in an independent discovery (cases = 1,852 and controls = 1,990) and replication (cases = 799 and controls = 778) cohort. Multivariate logistic regression was used to associate the plasma protein PGSs with LOAD diagnosis while controlling for age, sex, principal components 1-2, and the number of APOE-e4 alleles as covariates. After meta-analyzing the PGS-LOAD associations between the two cohorts, we then performed a two-sample Mendelian randomization (MR) analysis using the summary statistics of significant plasma protein level PGSs in the meta-analysis as an exposure, and a GWAS of clinically diagnosed LOAD (cases = 21,982, controls = 41,944) as an outcome to explore possible causal relationships between the two. Results: We identified four plasma protein level PGSs that were significantly associated (FDR-adjusted p < 0.05) with LOAD in a meta-analysis of the discovery and replication cohorts: CX3CL1, hepatocyte growth factor (HGF), TIE2, and matrix metalloproteinase-3 (MMP-3). When these four plasma proteins were used as exposures in MR with LOAD liability as the outcome, plasma levels of HGF were inferred to have a negative causal relationship with the disease when single-nucleotide polymorphisms (SNPs) used as instrumental variables were not restricted to cis-variants (OR/95%CI = 0.945/0.906-0.984, p = 0.005). Conclusion: Our results show that plasma HGF has a negative causal relationship with LOAD liability that is driven by pleiotropic SNPs possibly involved in other pathways. These findings suggest a low transferability between PGS and MR approaches, and future research should explore ways in which LOAD and the plasma proteome may interact.

Proteome-wide association studies using summary pQTL data of three tissues identified 30 risk genes of Alzheimer's disease dementia.

Background: Proteome-wide association study (PWAS) integrating proteomic data with genome-wide association study (GWAS) summary data is a powerful tool for studying Alzheimer's disease (AD) dementia. Existing PWAS analyses of AD often rely on the availability of individual-level proteomic and genetic data of a reference panel. Leveraging summary protein quantitative trait loci (pQTL) reference data of multiple AD-relevant tissues is expected to improve PWAS findings of AD dementia. Methods: We conducted PWAS of AD dementia by integrating publicly available summary pQTL data of brain, cerebrospinal fluid (CSF), and plasma tissues, with the latest GWAS summary data of AD dementia. For each target protein per tissue, we employed our recently published OTTERS tool to obtain omnibus PWAS p-value, to test whether the genetically regulated protein abundance in the corresponding tissue is associated with AD dementia. Protein-protein interactions and enriched pathways of identified significant PWAS risk genes were analyzed by STRING. The potential causal effects of these PWAS risk genes were assessed by probabilistic Mendelian randomization analyses. Results: We identified 30 unique significant PWAS risk genes for AD dementia, including 11 for brain, 9 for CSF, and 16 for plasma tissues. Four of these were shared by at least two tissues, and gene MAPK3 was found in all three tissues. We found that 11 of these PWAS risk genes were associated with AD or AD pathological hall marks as shown in GWAS Catalog; 18 of these were detected by transcriptome-wide association studies (TWAS); and 25 of these, including 8 out of 9 novel genes, were interconnected within a protein-protein interaction network involving the well-known AD risk gene APOE. Especially, these PWAS risk genes were enriched in immune response, glial cell proliferation, and high-density lipoprotein particle clearance pathways. Mediated causal effects were validated for 13 PWAS risk genes (43.3%). Conclusions: Our findings provide novel insights into the genetic mechanisms of AD dementia in brain, CSF, and plasma tissues, and targets for developing therapeutic interventions. We also demonstrated the effectiveness of integrating summary pQTL and GWAS data for mapping risk genes of complex human diseases.