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 (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, MRM-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.65e-08), ZCWPW1-PILRB (rs1476679, P=2.73e-32), CTSH-CTSH (rs3784539, P=2.88e-24) and HESX1-RETN (rs186108507, P=8.39e-08), 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.90e-7). DLB risk locus GBA showed positive trans effects on seven inter-related CSF levels in DLB patients only. No pQTLs were identified for frontotemporal dementia, either for the total sample as for analyses performed within FTD only. pQTL 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.

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.

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 proteomic data identified 23 risk genes of Alzheimer's disease.

Characterizing the genetic mechanisms underlying Alzheimer's disease (AD) dementia is crucial for developing new therapeutics. Proteome-wide association study (PWAS) integrating proteomics data with genome-wide association study (GWAS) summary data was shown as a powerful tool for detecting risk genes. The identified PWAS risk genes can be interpretated as having genetic effects mediated through the genetically regulated protein abundances. Existing PWAS analyses of AD often rely on the availability of individual-level proteomics and genetics data of a reference cohort. Leveraging summary-level protein quantitative trait loci (pQTL) reference data of multiple relevant tissues is expected to improve PWAS findings for studying AD. Here, we applied our recently developed OTTERS tool to conduct PWAS of AD dementia, by leveraging summary-level pQTL data of brain, cerebrospinal fluid (CSF), and plasma tissues, and multiple statistical methods. For each target protein, imputation models of the protein abundance with genetic predictors were trained from summary-level pQTL data, estimating a set of pQTL weights for considered genetic predictors. PWAS p-values were obtained by integrating GWAS summary data of AD dementia with estimated pQTL weights. PWAS p-values from multiple statistical methods were combined by the aggregated Cauchy association test to yield one omnibus PWAS p-value for the target protein. We identified significant PWAS risk genes through omnibus PWAS p-values and analyzed their protein-protein interactions using STRING. Their potential causal effects were assessed by the probabilistic Mendelian randomization (PMR-Egger). As a result, we identified a total of 23 significant PWAS risk genes for AD dementia in brain, CSF, and plasma tissues, including 7 novel findings. We showed that 15 of these risk genes were interconnected within a protein-protein interaction network involving the well-known AD risk gene of APOE and 5 novel findings, and enriched in immune functions and lipids pathways including positive regulation of immune system process, positive regulation of macrophage proliferation, humoral immune response, and high-density lipoprotein particle clearance. Existing biological evidence was found to relate our novel findings with AD. We validated the mediated causal effects of 14 risk genes (60.8%). In conclusion, we identified both known and novel PWAS risk genes, providing novel insights into the genetic mechanisms in brain, CSF, and plasma tissues, and targeted therapeutics development of AD dementia. Our study also demonstrated the effectiveness of integrating public available summary-level pQTL data with GWAS summary data for mapping risk genes of complex human diseases.

Abnormal Brain Protein Abundance and Cross-tissue mRNA Expression in Amyotrophic Lateral Sclerosis.

Due to the limitations of the present risk genes in understanding the etiology of amyotrophic lateral sclerosis (ALS), it is necessary to find additional causative genes utilizing novel approaches. In this study, we conducted a two-stage proteome-wide association study (PWAS) using ALS genome-wide association study (GWAS) data (N = 152,268) and two distinct human brain protein quantitative trait loci (pQTL) datasets (ROSMAP N = 376 and Banner N = 152) to identify ALS risk genes and prioritized candidate genes with Mendelian randomization (MR) and Bayesian colocalization analysis. Next, we verified the aberrant expression of risk genes in multiple tissues, including lower motor neurons, skeletal muscle, and whole blood. Six ALS risk genes (SCFD1, SARM1, TMEM175, BCS1L, WIPI2, and DHRS11) were found during the PWAS discovery phase, and SARM1 and BCS1L were confirmed during the validation phase. The following MR (p = 2.10 × 10-7) and Bayesian colocalization analysis (ROSMAP PP4 = 0.999, Banner PP4 = 0.999) confirmed the causal association between SARM1 and ALS. Further differential expression analysis revealed that SARM1 was markedly downregulated in lower motor neurons (p = 7.64 × 10-3), skeletal muscle (p = 9.34 × 10-3), and whole blood (p = 1.94 × 10-3). Our findings identified some promising protein candidates for future investigation as therapeutic targets. The dysregulation of SARM1 in multiple tissues provides a new way to explain ALS pathology.

CSF N-acylethanolamine acid amidase level and Parkinson's disease risk: A mendelian randomization study.

BACKGROUND: Neuroinflammation is involved in the progression of Parkinson's disease (PD), and N-acylethanolamine acid amidase (NAAA) is involved in regulating inflammation by hydrolyzing bioactive lipid mediators called N-acylethanolamines (NAEs). However, the causal relationship between cerebrospinal fluid (CSF) NAAA protein levels and the risk of PD remains unclear. This study aimed to explore the causal effect of CSF NAAA levels on PD risk through Mendelian randomization (MR) analysis. METHOD: Genome-wide association study (GWAS) summary statistics for CSF NAAA protein quantitative trait loci (pQTL) and GWAS summary statistics for PD were obtained from publicly available databases. Inverse-variance weighted (IVW) was the main causal estimation method for MR analysis. In addition, the maximum likelihood, MR Egger regression, and weighted median were used to supplement the IVW results. Finally, various sensitivity tests were performed to verify the reliability of the MR findings. RESULTS: In the initial MR analysis, the IVW showed that CSF NAAA protein levels significantly increased PD risk (odds ratio [OR] = 1.17, 95% confidence interval [CI]: 1.01-1.35, P = 0.031). This finding was further validated in a replicate MR analysis (OR = 1.20, 95% CI: 1.02-1.41, P = 0.027). Sensitivity analysis showed that MR results were stable and not affected by heterogeneity and horizontal pleiotropy. CONCLUSION: The present MR study supports a causal relationship between elevated CSF NAAA protein levels and increased PD risk.

Proteome-Wide Multicenter Mendelian Randomization Analysis to Identify Novel Therapeutic Targets for Lung Cancer.

INTRODUCTION: Lung cancer (LC) remains a leading cause of cancer mortality worldwide, underscoring the urgent need for novel therapeutic targets. The integration of Mendelian randomization (MR) with proteomic data presents a novel approach to identifying potential targets for LC treatment. METHODS: This study utilized a proteome-wide MR analysis, leveraging publicly available data from genome-wide association studies (GWAS) and protein quantitative trait loci (pQTL) studies. We analyzed genetic association data for LC from the TRICL-ILCCO Consortium and proteomic data from the Decode cohort. The MR framework was employed to estimate the causal effects of specific proteins on LC risk, supplemented by external validation, co-localization analyses, and exploration of protein-protein interaction (PPI) networks. RESULTS: Our analysis identified five proteins (TFPI, ICAM5, SFTPB, COL6A3, EPHB1) with significant associations to LC risk. External validation confirmed the potential therapeutic relevance of ICAM5 and SFTPB. Co-localization analyses and PPI network exploration provided further insights into the biological pathways involved and their potential mechanistic roles in LC pathogenesis. CONCLUSION: The study highlights the power of integrating genomic and proteomic data through MR analysis to uncover novel therapeutic targets for lung cancer. The identified proteins, particularly ICAM5 and SFTPB, offer promising directions for future research and development of targeted therapies, demonstrating the potential to advance personalized medicine in lung cancer treatment.

The genetic regulation of protein expression in cerebrospinal fluid.

Studies of the genetic regulation of cerebrospinal fluid (CSF) proteins may reveal pathways for treatment of neurological diseases. 398 proteins in CSF were measured in 1,591 participants from the BioFINDER study. Protein quantitative trait loci (pQTL) were identified as associations between genetic variants and proteins, with 176 pQTLs for 145 CSF proteins (P < 1.25 × 10-10 , 117 cis-pQTLs and 59 trans-pQTLs). Ventricular volume (measured with brain magnetic resonance imaging) was a confounder for several pQTLs. pQTLs for CSF and plasma proteins were overall correlated, but CSF-specific pQTLs were also observed. Mendelian randomization analyses suggested causal roles for several proteins, for example, ApoE, CD33, and GRN in Alzheimer's disease, MMP-10 in preclinical Alzheimer's disease, SIGLEC9 in amyotrophic lateral sclerosis, and CD38, GPNMB, and ADAM15 in Parkinson's disease. CSF levels of GRN, MMP-10, and GPNMB were altered in Alzheimer's disease, preclinical Alzheimer's disease, and Parkinson's disease, respectively. These findings point to pathways to be explored for novel therapies. The novel finding that ventricular volume confounded pQTLs has implications for design of future studies of the genetic regulation of the CSF proteome.

Integration of Candida albicans-induced single-cell gene expression data and secretory protein concentrations reveal genetic regulators of inflammation.

Both gene expression and protein concentrations are regulated by genetic variants. Exploring the regulation of both eQTLs and pQTLs simultaneously in a context- and cell-type dependent manner may help to unravel mechanistic basis for genetic regulation of pQTLs. Here, we performed meta-analysis of Candida albicans-induced pQTLs from two population-based cohorts and intersected the results with Candida-induced cell-type specific expression association data (eQTL). This revealed systematic differences between the pQTLs and eQTL, where only 35% of the pQTLs significantly correlated with mRNA expressions at single cell level, indicating the limitation of eQTLs use as a proxy for pQTLs. By taking advantage of the tightly co-regulated pattern of the proteins, we also identified SNPs affecting protein network upon Candida stimulations. Colocalization of pQTLs and eQTLs signals implicated several genomic loci including MMP-1 and AMZ1. Analysis of Candida-induced single cell gene expression data implicated specific cell types that exhibit significant expression QTLs upon stimulation. By highlighting the role of trans-regulatory networks in determining the abundance of secretory proteins, our study serve as a framework to gain insights into the mechanisms of genetic regulation of protein levels in a context-dependent manner.