Identification of quantitative trait loci underlying proteome variation in human lymphoblastoid cells

Proteome-wide Systems Genetics to Identify Functional Regulators of Complex Traits

Proteomic technologies now enable the rapid quantification of thousands of proteins across genetically diverse samples. Integration of these data with systems-genetics analyses is a powerful approach to identify new regulators of economically important or disease-relevant phenotypes in various populations. In this review, we summarize the latest proteomic technologies and discuss technical challenges for their use in population studies. We demonstrate how the analysis of correlation structure and loci mapping can be used to identify genetic factors regulating functional protein networks and complex traits. Finally, we provide an extensive summary of the use of proteome-wide systems genetics throughout fungi, plant, and animal kingdoms and discuss the power of this approach to identify candidate regulators and drug targets in large human consortium studies.

QTL Analysis Beyond eQTLs

Expression quantitative trait locus (eQTL) analysis is a powerful method to understand the association between genetic variant and gene expression; it also has potential impact for the study of transcription medicine for human complex disease. In the past two decades, the researchers focus on studying the eQTL, while more and more evidence shows that the regulatory genetic variants locating noncoding region have strong effect for the gene expression. More and more researchers working on eQTL analysis realize the importance of other types of QTLs beyond eQTL. In this chapter, we will explore some QTLs beyond eQTLs that show the regulatory association with eQTLs and explain the underlying link among these types of QTLs.

The MS4A gene cluster is a key modulator of soluble TREM2 and Alzheimer's disease risk

Soluble triggering receptor expressed on myeloid cells 2 (sTREM2) in cerebrospinal fluid (CSF) has been associated with Alzheimer's disease (AD). TREM2 plays a critical role in microglial activation, survival, and phagocytosis; however, the pathophysiological role of sTREM2 in AD is not well understood. Understanding the role of sTREM2 in AD may reveal new pathological mechanisms and lead to the identification of therapeutic targets. We performed a genome-wide association study (GWAS) to identify genetic modifiers of CSF sTREM2 obtained from the Alzheimer's Disease Neuroimaging Initiative. Common variants in the membrane-spanning 4-domains subfamily A (MS4A) gene region were associated with CSF sTREM2 concentrations (rs1582763; P = 1.15 × 10-15); this was replicated in independent datasets. The variants associated with increased CSF sTREM2 concentrations were associated with reduced AD risk and delayed age at onset of disease. The single-nucleotide polymorphism rs1582763 modified expression of the MS4A4A and MS4A6A genes in multiple tissues, suggesting that one or both of these genes are important for modulating sTREM2 production. Using human macrophages as a proxy for microglia, we found that MS4A4A and TREM2 colocalized on lipid rafts at the plasma membrane, that sTREM2 increased with MS4A4A overexpression, and that silencing of MS4A4A reduced sTREM2 production. These genetic, molecular, and cellular findings suggest that MS4A4A modulates sTREM2. These findings also provide a mechanistic explanation for the original GWAS signal in the MS4A locus for AD risk and indicate that TREM2 may be involved in AD pathogenesis not only in TREM2 risk-variant carriers but also in those with sporadic disease.

Genomic and Phenomic Research in the 21st Century

The field of human genomics has changed dramatically over time. Initial genomic studies were predominantly restricted to rare disorders in small families. Over the past decade, researchers changed course from family-based studies and instead focused on common diseases and traits in populations of unrelated individuals. With further advancements in biobanking, computer science, electronic health record (EHR) data, and more affordable high-throughput genomics, we are experiencing a new paradigm in human genomic research. Rapidly changing technologies and resources now make it possible to study thousands of diseases simultaneously at the genomic level. This review will focus on these advancements as scientists begin to incorporate phenome-wide strategies in human genomic research to understand the etiology of human diseases and develop new drugs to treat them.

The human brainome: network analysis identifies HSPA2 as a novel Alzheimer’s disease target

Our hypothesis is that changes in gene and protein expression are crucial to the development of late-onset Alzheimer’s disease. Previously we examined how DNA alleles control downstream expression of RNA transcripts and how those relationships are changed in late-onset Alzheimer’s disease. We have now examined how proteins are incorporated into networks in two separate series and evaluated our outputs in two different cell lines. Our pipeline included the following steps: (i) predicting expression quantitative trait loci; (ii) determining differential expression; (iii) analysing networks of transcript and peptide relationships; and (iv) validating effects in two separate cell lines. We performed all our analysis in two separate brain series to validate effects. Our two series included 345 samples in the first set (177 controls, 168 cases; age range 65–105; 58% female; KRONOSII cohort) and 409 samples in the replicate set (153 controls, 141 cases, 115 mild cognitive impairment; age range 66–107; 63% female; RUSH cohort). Our top target is heat shock protein family A member 2 (HSPA2), which was identified as a key driver in our two datasets. HSPA2 was validated in two cell lines, with overexpression driving further elevation of amyloid-β40 and amyloid-β42 levels in APP mutant cells, as well as significant elevation of microtubule associated protein tau and phosphorylated-tau in a modified neuroglioma line. This work further demonstrates that studying changes in gene and protein expression is crucial to understanding late onset disease and further nominates HSPA2 as a specific key regulator of late-onset Alzheimer’s disease processes.10.1093/brain/awy215_video1awy215media15824729224001.

Allele-specific SHAPE-MaP assessment of the effects of somatic variation and protein binding on mRNA structure

The impact of inherited and somatic mutations on messenger RNA (mRNA) structure remains poorly understood. Recent technological advances that leverage next-generation sequencing to obtain experimental structure data, such as SHAPE-MaP, can reveal structural effects of mutations, especially when these data are incorporated into structure modeling. Here, we analyze the ability of SHAPE-MaP to detect the relatively subtle structural changes caused by single-nucleotide mutations. We find that allele-specific sorting greatly improved our detection ability. Thus, we used SHAPE-MaP with a novel combination of clone-free robotic mutagenesis and allele-specific sorting to perform a rapid, comprehensive survey of noncoding somatic and inherited riboSNitches in two cancer-associated mRNAs, TPT1 and LCP1 Using rigorous thermodynamic modeling of the Boltzmann suboptimal ensemble, we identified a subset of mutations that change TPT1 and LCP1 RNA structure, with approximately 14% of all variants identified as riboSNitches. To confirm that these in vitro structures were biologically relevant, we tested how dependent TPT1 and LCP1 mRNA structures were on their environments. We performed SHAPE-MaP on TPT1 and LCP1 mRNAs in the presence or absence of cellular proteins and found that both mRNAs have similar overall folds in all conditions. RiboSNitches identified within these mRNAs in vitro likely exist under biological conditions. Overall, these data reveal a robust mRNA structural landscape where differences in environmental conditions and most sequence variants do not significantly alter RNA structural ensembles. Finally, predicting riboSNitches in mRNAs from sequence alone remains particularly challenging; these data will provide the community with benchmarks for further algorithmic development.

Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis

Background

There is a pressing need in rheumatoid arthritis (RA) to identify patients who will not respond to first-line disease-modifying anti-rheumatic drugs (DMARD). We explored whether differences in genomic architecture represented by a chromosome conformation signature (CCS) in blood taken from early RA patients before methotrexate (MTX) treatment could assist in identifying non-response to DMARD and, whether there is an association between such a signature and RA specific expression quantitative trait loci (eQTL).

Methods

We looked for the presence of a CCS in blood from early RA patients commencing MTX using chromosome conformation capture by EpiSwitch™. Using blood samples from MTX responders, non-responders and healthy controls, a custom designed biomarker discovery array was refined to a 5-marker CCS that could discriminate between responders and non-responders to MTX. We cross-validated the predictive power of the CCS by generating 150 randomized groups of 59 early RA patients (30 responders and 29 non-responders) before MTX treatment. The CCS was validated using a blinded, independent cohort of 19 early RA patients (9 responders and 10 non-responders). Last, the loci of the CCS markers were mapped to RA-specific eQTL.

Results

We identified a 5-marker CCS that could identify, at baseline, responders and non-responders to MTX. The CCS consisted of binary chromosome conformations in the genomic regions of IFNAR1, IL-21R, IL-23, CXCL13 and IL-17A. When tested on a cohort of 59 RA patients, the CCS provided a negative predictive value of 90.0% for MTX response. When tested on a blinded independent validation cohort of 19 early RA patients, the signature demonstrated a true negative response rate of 86 and a 90% sensitivity for detection of non-responders to MTX. Only conformations in responders mapped to RA-specific eQTL.

Conclusions

Here we demonstrate that detection of a CCS in blood in early RA is able to predict inadequate response to MTX with a high degree of accuracy. Our results provide a proof of principle that a priori stratification of response to MTX is possible, offering a mechanism to provide alternative treatments for non-responders to MTX earlier in the course of the disease.

Electronic supplementary material

The online version of this article (10.1186/s12967-018-1387-9) contains supplementary material, which is available to authorized users.

Human iPSC-derived neurons and lymphoblastoid cells for personalized medicine research in neuropsychiatric disorders

The development and clinical implementation of personalized medicine crucially depends on the availability of high-quality human biosamples; animal models, although capable of modeling complex human diseases, cannot reflect the large variation in the human genome, epigenome, transcriptome, proteome, and metabolome. Although the biosamples available from public biobanks that store human tissues and cells may represent the large human diversity for most diseases, these samples are not always sufficient for developing biomarkers for patient-tailored therapies for neuropsychiatric disorders. Postmortem human tissues are available from many biobanks; nevertheless, collections of neuronal human cells from large patient cohorts representing the human diversity remain scarce. Two tools are gaining popularity for personalized medicine research on neuropsychiatric disorders: human induced pluripotent stem cell-derived neurons and human lymphoblastoid cell lines. This review examines and contrasts the advantages and limitations of each tool for personalized medicine research.

Emerging Affinity-Based Proteomic Technologies for Large-Scale Plasma Profiling in Cardiovascular Disease

Plasma biomarkers that reflect molecular states of the cardiovascular system are central for clinical decision making. Routinely used plasma biomarkers include troponins, natriuretic peptides, and lipoprotein particles, yet interrogate only a modest subset of pathways relevant to cardiovascular disease. Systematic profiling of a larger portion of circulating plasma proteins (the plasma proteome) will provide opportunities for unbiased discovery of novel markers to improve diagnostic or predictive accuracy. In addition, proteomic profiling may inform pathophysiological understanding and point to novel therapeutic targets. Obstacles for comprehensive proteomic profiling include the immense size and structural heterogeneity of the proteome, and the broad range of abundance levels, as well. Proteome-wide, untargeted profiling can be performed in tissues and cells with tandem mass spectrometry. However, applications to plasma are limited by the need for complex preanalytical sample preparation stages limiting sample throughput. Multiplexing of targeted methods based on capture and detection of specific proteins are therefore receiving increasing attention in plasma proteomics. Immunoaffinity assays are the workhorse for measuring individual proteins but have been limited for proteomic applications by long development times, cross-reactivity preventing multiplexing, specificity issues, and incomplete sensitivity to detect proteins in the lower range of the abundance spectrum (below picograms per milliliter). Emerging technologies to address these issues include nucleotide-labeled immunoassays and aptamer reagents that can be automated for efficient multiplexing of thousands of proteins at high sample throughput, coupling of affinity capture methods to mass spectrometry for improved specificity, and ultrasensitive detection systems to measure low-abundance proteins. In addition, proteomics can now be integrated with modern genomics tools to comprehensively relate proteomic profiles to genetic variants, which may both influence binding of affinity reagents and serve to validate the target specificity of affinity assays. The application of deep quantitative proteomic profiling to large cohorts has thus become increasingly feasible with emerging affinity methods. The aims of this article are to provide the broad readership of Circulation with a timely overview of emerging methods for affinity proteomics and recent progress in cardiovascular medicine based on such methods.

Connecting genetic risk to disease end points through the human blood plasma proteome

Genome-wide association studies (GWAS) with intermediate phenotypes, like changes in metabolite and protein levels, provide functional evidence to map disease associations and translate them into clinical applications. However, although hundreds of genetic variants have been associated with complex disorders, the underlying molecular pathways often remain elusive. Associations with intermediate traits are key in establishing functional links between GWAS-identified risk-variants and disease end points. Here we describe a GWAS using a highly multiplexed aptamer-based affinity proteomics platform. We quantify 539 associations between protein levels and gene variants (pQTLs) in a German cohort and replicate over half of them in an Arab and Asian cohort. Fifty-five of the replicated pQTLs are located in trans. Our associations overlap with 57 genetic risk loci for 42 unique disease end points. We integrate this information into a genome-proteome network and provide an interactive web-tool for interrogations. Our results provide a basis for novel approaches to pharmaceutical and diagnostic applications.

Individual genetic variation can affect the levels of protein in blood, but detailed data sets linking these two types of data are rare. Here, the authors carry out a genome-wide association study of levels of over a thousand different proteins, and describe many new SNP-protein interactions.

Associations Between Common and Rare Exonic Genetic Variants and Serum Levels of 20 Cardiovascular-Related Proteins: The Tromsø Study

Supplemental Digital Content is available in the text.

Background—

Genetic variation can be used to study causal relationships between biomarkers and diseases. Here, we identify new common and rare genetic variants associated with cardiovascular-related protein levels (protein quantitative trait loci [pQTLs]). We functionally annotate these pQTLs, predict and experimentally confirm a novel molecular interaction, and determine which pQTLs are associated with diseases and physiological phenotypes.

Methods and Results—

As part of a larger case–control study of venous thromboembolism, serum levels of 51 proteins implicated in cardiovascular diseases were measured in 330 individuals from the Tromsø Study. Exonic genetic variation near each protein’s respective gene (cis) was identified using sequencing and arrays. Using single site and gene-based tests, we identified 27 genetic associations between pQTLs and the serum levels of 20 proteins: 14 associated with common variation in cis, of which 6 are novel (ie, not previously reported); 7 associations with rare variants in cis, of which 4 are novel; and 6 associations in trans. Of the 20 proteins, 15 were associated with single sites and 7 with rare variants. cis-pQTLs for kallikrein and F12 also show trans associations for proteins (uPAR, kininogen) known to be cleaved by kallikrein and with NTproBNP. We experimentally demonstrate that kallikrein can cleave proBNP (NTproBNP precursor) in vitro. Nine of the pQTLs have previously identified associations with 17 disease and physiological phenotypes.

Conclusions—

We have identified cis and trans genetic variation associated with the serum levels of 20 proteins and utilized these pQTLs to study molecular mechanisms underlying disease and physiological phenotypes.

Autoimmune diseases - connecting risk alleles with molecular traits of the immune system

Genome-wide strategies have driven the discovery of more than 300 susceptibility loci for autoimmune diseases. However, for almost all loci, understanding of the mechanisms leading to autoimmunity remains limited, and most variants that are likely to be causal are in non-coding regions of the genome. A critical next step will be to identify the in vivo and ex vivo immunophenotypes that are affected by risk variants. To do this, key cell types and cell states that are implicated in autoimmune diseases will need to be defined. Functional genomic annotations from these cell types and states can then be used to resolve candidate genes and causal variants. Together with longitudinal studies, this approach may yield pivotal insights into how autoimmunity is triggered.

Genetic studies of plasma analytes identify novel potential biomarkers for several complex traits

Genome-wide association studies of 146 plasma protein levels in 818 individuals revealed 56 genome-wide significant associations (28 novel) with 47 analytes. Loci associated with plasma levels of 39 proteins tested have been previously associated with various complex traits such as heart disease, inflammatory bowel disease, Type 2 diabetes and multiple sclerosis. These data suggest that these plasma protein levels may constitute informative endophenotypes for these complex traits. We found three potential pleiotropic genes: ABO for plasma SELE and ACE levels, FUT2 for CA19-9 and CEA plasma levels and APOE for ApoE and CRP levels. We also found multiple independent signals in loci associated with plasma levels of ApoH, CA19-9, FetuinA, IL6r and LPa. Our study highlights the power of biological traits for genetic studies to identify genetic variants influencing clinically relevant traits, potential pleiotropic effects and complex disease associations in the same locus.

Allele-Selective Transcriptome Recruitment to Polysomes Primed for Translation: Protein-Coding and Noncoding RNAs, and RNA Isoforms

mRNA translation into proteins is highly regulated, but the role of mRNA isoforms, noncoding RNAs (ncRNAs), and genetic variants remains poorly understood. mRNA levels on polysomes have been shown to correlate well with expressed protein levels, pointing to polysomal loading as a critical factor. To study regulation and genetic factors of protein translation we measured levels and allelic ratios of mRNAs and ncRNAs (including microRNAs) in lymphoblast cell lines (LCL) and in polysomal fractions. We first used targeted assays to measure polysomal loading of mRNA alleles, confirming reported genetic effects on translation of OPRM1 and NAT1, and detecting no effect of rs1045642 (3435C>T) in ABCB1 (MDR1) on polysomal loading while supporting previous results showing increased mRNA turnover of the 3435T allele. Use of high-throughput sequencing of complete transcript profiles (RNA-Seq) in three LCLs revealed significant differences in polysomal loading of individual RNA classes and isoforms. Correlated polysomal distribution between protein-coding and non-coding RNAs suggests interactions between them. Allele-selective polysome recruitment revealed strong genetic influence for multiple RNAs, attributable either to differential expression of RNA isoforms or to differential loading onto polysomes, the latter defining a direct genetic effect on translation. Genes identified by different allelic RNA ratios between cytosol and polysomes were enriched with published expression quantitative trait loci (eQTLs) affecting RNA functions, and associations with clinical phenotypes. Polysomal RNA-Seq combined with allelic ratio analysis provides a powerful approach to study polysomal RNA recruitment and regulatory variants affecting protein translation.

Extensive allele-specific translational regulation in hybrid mice

Translational regulation is mediated through the interaction between diffusible trans-factors and cis-elements residing within mRNA transcripts. In contrast to extensively studied transcriptional regulation, cis-regulation on translation remains underexplored. Using deep sequencing-based transcriptome and polysome profiling, we globally profiled allele-specific translational efficiency for the first time in an F1 hybrid mouse. Out of 7,156 genes with reliable quantification of both alleles, we found 1,008 (14.1%) exhibiting significant allelic divergence in translational efficiency. Systematic analysis of sequence features of the genes with biased allelic translation revealed that local RNA secondary structure surrounding the start codon and proximal out-of-frame upstream AUGs could affect translational efficiency. Finally, we observed that the cis-effect was quantitatively comparable between transcriptional and translational regulation. Such effects in the two regulatory processes were more frequently compensatory, suggesting that the regulation at the two levels could be coordinated in maintaining robustness of protein expression.

Genomic variation. Impact of regulatory variation from RNA to protein

The phenotypic consequences of expression quantitative trait loci (eQTLs) are presumably due to their effects on protein expression levels. Yet the impact of genetic variation, including eQTLs, on protein levels remains poorly understood. To address this, we mapped genetic variants that are associated with eQTLs, ribosome occupancy (rQTLs), or protein abundance (pQTLs). We found that most QTLs are associated with transcript expression levels, with consequent effects on ribosome and protein levels. However, eQTLs tend to have significantly reduced effect sizes on protein levels, which suggests that their potential impact on downstream phenotypes is often attenuated or buffered. Additionally, we identified a class of cis QTLs that affect protein abundance with little or no effect on messenger RNA or ribosome levels, which suggests that they may arise from differences in posttranslational regulation.

Identification and validation of genetic variants that influence transcription factor and cell signaling protein levels

Many genetic variants associated with human disease have been found to be associated with alterations in mRNA expression. Although it is commonly assumed that mRNA expression changes will lead to consequent changes in protein levels, methodological challenges have limited our ability to test the degree to which this assumption holds true. Here, we further developed the micro-western array approach and globally examined relationships between human genetic variation and cellular protein levels. We collected more than 250,000 protein level measurements comprising 441 transcription factor and signaling protein isoforms across 68 Yoruba (YRI) HapMap lymphoblastoid cell lines (LCLs) and identified 12 cis and 160 trans protein level QTLs (pQTLs) at a false discovery rate (FDR) of 20%. Whereas up to two thirds of cis mRNA expression QTLs (eQTLs) were also pQTLs, many pQTLs were not associated with mRNA expression. Notably, we replicated and functionally validated a trans pQTL relationship between the KARS lysyl-tRNA synthetase locus and levels of the DIDO1 protein. This study demonstrates proof of concept in applying an antibody-based microarray approach to iteratively measure the levels of human proteins and relate these levels to human genome variation and other genomic data sets. Our results suggest that protein-based mechanisms might functionally buffer genetic alterations that influence mRNA expression levels and that pQTLs might contribute phenotypic diversity to a human population independently of influences on mRNA expression.

Proteomic studies related to genetic determinants of variability in protein concentrations

Genetic variation has multiple effects on the proteome. It may influence the expression level of proteins, modify their sequences through single nucleotide polymorphisms, the occurrence of allelic variants, or alternative splicing (ASP) events. This perspective paper summarizes the major effects of genetic variability on protein expression and isoforms and provides an overview of proteomics techniques and methods that allow studying the effects of genetic variability at different levels of the proteome. The paper provides an overview of recent quantitative trait loci studies performed to explore the effect of genetic variation on protein expression (pQTL). Finally it gives a perspective view on advances in proteomics technology and the role of the Chromosome-Centric Human Proteome Project (C-HPP) by creating large-scale resources that may facilitate performing more comprehensive pQTL experiments in the future.

Influence of genetic variation on plasma protein levels in older adults using a multi-analyte panel

Proteins, widely studied as potential biomarkers, play important roles in numerous physiological functions and diseases. Genetic variation may modulate corresponding protein levels and point to the role of these variants in disease pathophysiology. Effects of individual single nucleotide polymorphisms (SNPs) within a gene were analyzed for corresponding plasma protein levels using genome-wide association study (GWAS) genotype data and proteomic panel data with 132 quality-controlled analytes from 521 Caucasian participants in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort. Linear regression analysis detected 112 significant (Bonferroni threshold p = 2.44×10⁻⁵) associations between 27 analytes and 112 SNPs. 107 out of these 112 associations were tested in the Indiana Memory and Aging Study (IMAS) cohort for replication and 50 associations were replicated at uncorrected p<0.05 in the same direction of effect as those in the ADNI. We identified multiple novel associations including the association of rs7517126 with plasma complement factor H-related protein 1 (CFHR1) level at p<1.46×10⁻⁶⁰, accounting for 40 percent of total variation of the protein level. We serendipitously found the association of rs6677604 with the same protein at p<9.29×10⁻¹¹². Although these two SNPs were not in the strong linkage disequilibrium, 61 percent of total variation of CFHR1 was accounted for by rs6677604 without additional variation by rs7517126 when both SNPs were tested together. 78 other SNP-protein associations in the ADNI sample exceeded genome-wide significance (5×10⁻⁸). Our results confirmed previously identified gene-protein associations for interleukin-6 receptor, chemokine CC-4, angiotensin-converting enzyme, and angiotensinogen, although the direction of effect was reversed in some cases. This study is among the first analyses of gene-protein product relationships integrating multiplex-panel proteomics and targeted genes extracted from a GWAS array. With intensive searches taking place for proteomic biomarkers for many diseases, the role of genetic variation takes on new importance and should be considered in interpretation of proteomic results.

Integrative genomics: quantifying significance of phenotype-genotype relationships from multiple sources of high-throughput data

Given recent advances in the generation of high-throughput data such as whole-genome genetic variation and transcriptome expression, it is critical to come up with novel methods to integrate these heterogeneous datasets and to assess the significance of identified phenotype-genotype relationships. Recent studies show that genome-wide association findings are likely to fall in loci with gene regulatory effects such as expression quantitative trait loci (eQTLs), demonstrating the utility of such integrative approaches. When genotype and gene expression data are available on the same individuals, we and others developed methods wherein top phenotype-associated genetic variants are prioritized if they are associated, as eQTLs, with gene expression traits that are themselves associated with the phenotype. Yet there has been no method to determine an overall p-value for the findings that arise specifically from the integrative nature of the approach. We propose a computationally feasible permutation method that accounts for the assimilative nature of the method and the correlation structure among gene expression traits and among genotypes. We apply the method to data from a study of cellular sensitivity to etoposide, one of the most widely used chemotherapeutic drugs. To our knowledge, this study is the first statistically sound quantification of the overall significance of the genotype-phenotype relationships resulting from applying an integrative approach. This method can be easily extended to cases in which gene expression data are replaced by other molecular phenotypes of interest, e.g., microRNA or proteomic data. This study has important implications for studies seeking to expand on genetic association studies by the use of omics data. Finally, we provide an R code to compute the empirical false discovery rate when p-values for the observed and simulated phenotypes are available.

Identification of genetic variants influencing the human plasma proteome

Genetic variants influencing the transcriptome have been extensively studied. However, the impact of the genetic factors on the human proteome is largely unexplored, mainly due to lack of suitable high-throughput methods. Here we present unique and comprehensive identification of genetic variants affecting the human plasma protein profile by combining high-throughput and high-resolution mass spectrometry (MS) with genome-wide SNP data. We identified and quantified the abundance of 1,056 tryptic-digested peptides, representing 163 proteins in the plasma of 1,060 individuals from two population-based cohorts. The abundance level of almost one-fifth (19%) of the peptides was found to be heritable, with heritability ranging from 0.08 to 0.43. The levels of 60 peptides from 25 proteins, 15% of the proteins studied, were influenced by cis-acting SNPs. We identified and replicated individual cis-acting SNPs (combined P value ranging from 3.1 × 10(-52) to 2.9 × 10(-12)) influencing 11 peptides from 5 individual proteins. These SNPs represent both regulatory SNPs and nonsynonymous changes defining well-studied disease alleles such as the ε4 allele of apolipoprotein E (APOE), which has been shown to increase risk of Alzheimer's disease. Our results show that high-throughput mass spectrometry represents a promising method for large-scale characterization of the human proteome, allowing for both quantification and sequencing of individual proteins. Abundance and peptide composition of a protein plays an important role in the etiology, diagnosis, and treatment of a number of diseases. A better understanding of the genetic impact on the plasma proteome is therefore important for evaluating potential biomarkers and therapeutic agents for common diseases.

Identification of cis-regulatory variation influencing protein abundance levels in human plasma

Proteins are central to almost all cellular processes, and dysregulation of expression and function is associated with a range of disorders. A number of studies in human have recently shown that genetic factors significantly contribute gene expression variation. In contrast, very little is known about the genetic basis of variation in protein abundance in man. Here, we assayed the abundance levels of proteins in plasma from 96 elderly Europeans using a new aptamer-based proteomic technology and performed genome-wide local (cis-) regulatory association analysis to identify protein quantitative trait loci (pQTL). We detected robust cis-associations for 60 proteins at a false discovery rate of 5%. The most highly significant single nucleotide polymorphism detected was rs7021589 (false discovery rate, 2.5 × 10(-12)), mapped within the gene coding sequence of Tenascin C (TNC). Importantly, we identified evidence of cis-regulatory variation for 20 previously disease-associated genes encoding protein, including variants with strong evidence of disease association show significant association with protein abundance levels. These results demonstrate that common genetic variants contribute to the differences in protein abundance levels in human plasma. Identification of pQTLs will significantly enhance our ability to discover and comprehend the biological and functional consequences of loci identified from genome-wide association study of complex traits. This is the first large-scale genetic association study of proteins in plasma measured using a novel, highly multiplexed slow off-rate modified aptamer (SOMAmer) proteomic platform.

Targeted protein-omic methods are bridging the gap between proteomic and hypothesis-driven protein analysis approaches

While proteomic methods have illuminated many areas of biological protein space, many fundamental questions remain with regard to systems-level relationships between mRNAs, proteins and cell behaviors. While mass spectrometric methods offer a panoramic picture of the relative expression and modification of large numbers of proteins, they are neither optimal for the analysis of predefined targets across large numbers of samples nor for assessing differences in proteins between individual cells or cell compartments. Conversely, traditional antibody-based methods are effective at sensitively analyzing small numbers of proteins across small numbers of conditions, and can be used to analyze relative differences in protein abundance and modification between cells and cell compartments. However, traditional antibody-based approaches are not optimal for analyzing large numbers of protein abundances and modifications across many samples. In this article, we will review recent advances in methodologies and philosophies behind several microarray-based, intermediate-level, 'protein-omic' methods, including a focus on reverse-phase lysate arrays and micro-western arrays, which have been helpful for bridging gaps between large- and small-scale protein analysis approaches and have provided insight into the roles that protein systems play in several biological processes.

Proteomics in evolutionary ecology: linking the genotype with the phenotype

The study of the proteome (proteomics), which includes the dynamics of protein expression, regulation, interactions and its function, has played a less prominent role in evolutionary and ecological investigations in comparison with the study of the genome and transcriptome. There are, however, a number of arguments suggesting that this situation should change. First, the proteome is closer to the phenotype than the genome or the transcriptome, and as such may be more directly responsive to natural selection, and thus closely linked to adaptation. Second, there is evidence of a low correlation between protein and transcript expression levels across genes in many different organisms. Finally, there have been some recent important technological improvements in proteomics methods that make them feasible, practical and useful to address a wide range of evolutionary questions even in nonmodel organisms. The different proteomic methods, their limitations and problems when interpreting empirical data are described and discussed. In addition, the proteomic literature pertaining to evolutionary ecology is reviewed with examples, and potential applications of proteomics in a variety of evolutionary contexts are outlined. New proteomic research trends such as the study of posttranslational modifications and protein-protein interactions, as well as the combined use of the different -omics approaches, are discussed in relation to the development of a more functional and integrated perspective, needed for achieving a more comprehensive knowledge of evolutionary change.

Biomarker discovery from the top down: Protein biomarkers for efficient virus transmission by insects (Homoptera: Aphididae) discovered by coupling genetics and 2-D DIGE

Yellow dwarf viruses cause the most economically important virus diseases of cereal crops worldwide and are vectored by aphids. The identification of vector proteins mediating virus transmission is critical to develop sustainable virus management practices and to understand viral strategies for circulative movement in all insect vectors. Previously, we applied 2-D DIGE to an aphid filial generation 2 population to identify proteins correlated with the transmission phenotype that were stably inherited and expressed in the absence of the virus. In the present study, we examined the expression of the DIGE candidates in previously unstudied, field-collected aphid populations. We hypothesized that the expression of proteins involved in virus transmission could be clinically validated in unrelated, virus transmission-competent, field-collected aphid populations. All putative biomarkers were expressed in the field-collected biotypes, and the expression of nine of these aligned with the virus transmission-competent phenotype. The strong conservation of the expression of the biomarkers in multiple field-collected populations facilitates new and testable hypotheses concerning the genetics and biochemistry of virus transmission. Integration of these biomarkers into current aphid-scouting methodologies will enable rational strategies for vector control aimed at judicious use and development of precision pest control methods that reduce plant virus infection.

Genomewide analysis of inherited variation associated with phosphorylation of PI3K/AKT/mTOR signaling proteins

While there exists a wealth of information about genetic influences on gene expression, less is known about how inherited variation influences the expression and post-translational modifications of proteins, especially those involved in intracellular signaling. The PI3K/AKT/mTOR signaling pathway contains several such proteins that have been implicated in a number of diseases, including a variety of cancers and some psychiatric disorders. To assess whether the activation of this pathway is influenced by genetic factors, we measured phosphorylated and total levels of three key proteins in the pathway (AKT1, p70S6K, 4E-BP1) by ELISA in 122 lymphoblastoid cell lines from 14 families. Interestingly, the phenotypes with the highest proportion of genetic influence were the ratios of phosphorylated to total protein for two of the pathway members: AKT1 and p70S6K. Genomewide linkage analysis suggested several loci of interest for these phenotypes, including a linkage peak for the AKT1 phenotype that contained the AKT1 gene on chromosome 14. Linkage peaks for the phosphorylated:total protein ratios of AKT1 and p70S6K also overlapped on chromosome 3. We selected and genotyped candidate genes from under the linkage peaks, and several statistically significant associations were found. One polymorphism in HSP90AA1 was associated with the ratio of phosphorylated to total AKT1, and polymorphisms in RAF1 and GRM7 were associated with the ratio of phosphorylated to total p70S6K. These findings, representing the first genomewide search for variants influencing human protein phosphorylation, provide useful information about the PI3K/AKT/mTOR pathway and serve as a valuable proof of concept for studies integrating human genomics and proteomics.

Progress in mining the human proteome for disease applications

Currently available technologies allow in-depth analysis of multiple facets of the proteome that have clinical relevance and that complement current genomics-based approaches. Although some progress has been made in our knowledge of the human proteome in health and in disease, there is an urgent need to chart a coherent road map with clearly defined milestones to guide proteomics efforts. Areas of emphasis include: (1) building resources, (2) filling gaps in our understanding of biological variation, and (3) systematically characterizing proteome alterations that occur in well-defined disease states, all of which require an organized and collaborative effort.

Tangible benefits of the aphid Acyrthosiphon pisum genome sequencing for aphid proteomics: Enhancements in protein identification and data validation for homology-based proteomics

Homology-driven proteomics promises to reveal functional biology in insects with sparse genome sequence information. A proteomics study comparing plant virus transmission competent and refractive genotypes of the aphid Schizaphis graminum isolated numerous candidate proteins involved in virus transmission, but limited genome sequence information hampered their identification. The complete genome of the pea aphid, Acyrthosiphon pisum, released in 2008, enabled us to double the number of protein identifications beyond what was possible using available EST libraries and other insect sequences. This was concomitant with a dramatic increase of the number of MS and MS/MS peptide spectra matching the genome-derived protein sequence. LC-MS/MS proved to be the most robust method of peptide detection. Cross-matching spectral data to multiple EST sequences and error tolerant searching to identify amino acid substitutions enhanced the percent coverage of the Schizaphis graminum proteins. 2-D electrophoresis provided the protein pI and MW which enabled the refinement of the candidate protein selection and provided a measure of protein abundance when coupled to the spectral data. Thus, the homology-based proteomics pipeline for insects should include efforts to maximize the number of peptide matches to the protein to increase certainty in protein identification and relative protein abundance.

Genetics coupled to quantitative intact proteomics links heritable aphid and endosymbiont protein expression to circulative polerovirus transmission

Yellow dwarf viruses in the family Luteoviridae, which are the causal agents of yellow dwarf disease in cereal crops, are each transmitted most efficiently by different species of aphids in a circulative manner that requires the virus to interact with a multitude of aphid proteins. Aphid proteins differentially expressed in F2 Schizaphis graminum genotypes segregating for the ability to transmit Cereal yellow dwarf virus-RPV (CYDV-RPV) were identified using two-dimensional difference gel electrophoresis (DIGE) coupled to either matrix-assisted laser desorption ionization-tandem mass spectrometry or online nanoscale liquid chromatography coupled to electrospray tandem mass spectrometry. A total of 50 protein spots, containing aphid proteins and proteins from the aphid's obligate and maternally inherited bacterial endosymbiont, Buchnera, were identified as differentially expressed between transmission-competent and refractive aphids. Surprisingly, in virus transmission-competent F2 genotypes, the isoelectric points of the Buchnera proteins did not match those in the maternal Buchnera proteome as expected, but instead they aligned with the Buchnera proteome of the transmission-competent paternal parent. Among the aphid proteins identified, many were involved in energy metabolism, membrane trafficking, lipid signaling, and the cytoskeleton. At least eight aphid proteins were expressed as heritable, isoelectric point isoform pairs, one derived from each parental lineage. In the F2 genotypes, the expression of aphid protein isoforms derived from the competent parental lineage aligned with the virus transmission phenotype with high precision. Thus, these isoforms are candidate biomarkers for CYDV-RPV transmission in S. graminum. Our combined genetic and DIGE approach also made it possible to predict where several of the proteins may be expressed in refractive aphids with different barriers to transmission. Twelve proteins were predicted to act in the hindgut of the aphid, while six proteins were predicted to be associated with the accessory salivary glands or hemolymph. Knowledge of the proteins that regulate virus transmission and their predicted locations will aid in understanding the biochemical mechanisms regulating circulative virus transmission in aphids, as well as in identifying new targets to block transmission.