Human SRMAtlas: A Resource of Targeted Assays to Quantify the Complete Human Proteome.

The ability to reliably and reproducibly measure any protein of the human proteome in any tissue or cell type would be transformative for understanding systems-level properties as well as specific pathways in physiology and disease. Here, we describe the generation and verification of a compendium of highly specific assays that enable quantification of 99.7% of the 20,277 annotated human proteins by the widely accessible, sensitive, and robust targeted mass spectrometric method selected reaction monitoring, SRM. This human SRMAtlas provides definitive coordinates that conclusively identify the respective peptide in biological samples. We report data on 166,174 proteotypic peptides providing multiple, independent assays to quantify any human protein and numerous spliced variants, non-synonymous mutations, and post-translational modifications. The data are freely accessible as a resource at http://www.srmatlas.org/, and we demonstrate its utility by examining the network response to inhibition of cholesterol synthesis in liver cells and to docetaxel in prostate cancer lines.

A computational framework for the inference of protein complex remodeling from whole-proteome measurements.

Protein complexes are responsible for the enactment of most cellular functions. For the protein complex to form and function, its subunits often need to be present at defined quantitative ratios. Typically, global changes in protein complex composition are assessed with experimental approaches that tend to be time consuming. Here, we have developed a computational algorithm for the detection of altered protein complexes based on the systematic assessment of subunit ratios from quantitative proteomic measurements. We applied it to measurements from breast cancer cell lines and patient biopsies and were able to identify strong remodeling of HDAC2 epigenetic complexes in more aggressive forms of cancer. The presented algorithm is available as an R package and enables the inference of changes in protein complex states by extracting functionally relevant information from bottom-up proteomic datasets.

Systematic characterization of pan-cancer mutation clusters.

Cancer genome sequencing has shown that driver genes can often be distinguished not only by the elevated mutation frequency but also by specific nucleotide positions that accumulate changes at a high rate. However, properties associated with a residue's potential to drive tumorigenesis when mutated have not yet been systematically investigated. Here, using a novel methodological approach, we identify and characterize a compendium of 180 hotspot residues within 160 human proteins which occur with a significant frequency and are likely to have functionally relevant impact. We find that such mutations (i) are more prominent in proteins that can exist in the on and off state, (ii) reflect the identity of a tumor of origin, and (iii) often localize within interfaces which mediate interactions with other proteins or ligands. Following, we further examine structural data for human protein complexes and identify a number of additional protein interfaces that accumulate cancer mutations at a high rate. Jointly, these analyses suggest that disruption and dysregulation of protein interactions can be instrumental in switching functions of cancer proteins and activating downstream changes.

Mass Spectrometric Exploration of the Biochemical Basis of Living Systems.

Predicting how a system behaves under changing conditions is an essential component of science and engineering. The ability to make accurate predictions about the system indicates that it is well understood and provides the opportunity to simulate the response to conditions that would be empirically difficult or impossible to test. In the life sciences, the term systems biology was introduced to articulate the notion that the molecular and phenotypic response of a cell or organism to perturbations is the result of interplay of a multitude of molecules. The ability to predict the behavior of such complex molecular systems remains challenging and inevitably requires the involvement of different types of models and data that support them. In this article, we discuss a range of data-driven models that have proven particularly useful for predicting the behavior of biological systems at different levels of complexity and the matching data generation methods that support them. We specifically focus on predictions based on protein or proteome data generated by mass spectrometry. We describe three case studies that represent frequently encountered situations in systems biology.

Systematic cell-based phenotyping of missense alleles empowers rare variant association studies: a case for LDLR and myocardial infarction.

A fundamental challenge to contemporary genetics is to distinguish rare missense alleles that disrupt protein functions from the majority of alleles neutral on protein activities. High-throughput experimental tools to securely discriminate between disruptive and non-disruptive missense alleles are currently missing. Here we establish a scalable cell-based strategy to profile the biological effects and likely disease relevance of rare missense variants in vitro. We apply this strategy to systematically characterize missense alleles in the low-density lipoprotein receptor (LDLR) gene identified through exome sequencing of 3,235 individuals and exome-chip profiling of 39,186 individuals. Our strategy reliably identifies disruptive missense alleles, and disruptive-allele carriers have higher plasma LDL-cholesterol (LDL-C). Importantly, considering experimental data refined the risk of rare LDLR allele carriers from 4.5- to 25.3-fold for high LDL-C, and from 2.1- to 20-fold for early-onset myocardial infarction. Our study generates proof-of-concept that systematic functional variant profiling may empower rare variant-association studies by orders of magnitude.

RNAi-based functional profiling of loci from blood lipid genome-wide association studies identifies genes with cholesterol-regulatory function.

Genome-wide association studies (GWAS) are powerful tools to unravel genomic loci associated with common traits and complex human disease. However, GWAS only rarely reveal information on the exact genetic elements and pathogenic events underlying an association. In order to extract functional information from genomic data, strategies for systematic follow-up studies on a phenotypic level are required. Here we address these limitations by applying RNA interference (RNAi) to analyze 133 candidate genes within 56 loci identified by GWAS as associated with blood lipid levels, coronary artery disease, and/or myocardial infarction for a function in regulating cholesterol levels in cells. Knockdown of a surprisingly high number (41%) of trait-associated genes affected low-density lipoprotein (LDL) internalization and/or cellular levels of free cholesterol. Our data further show that individual GWAS loci may contain more than one gene with cholesterol-regulatory functions. Using a set of secondary assays we demonstrate for a number of genes without previously known lipid-regulatory roles (e.g. CXCL12, FAM174A, PAFAH1B1, SEZ6L, TBL2, WDR12) that knockdown correlates with altered LDL-receptor levels and/or that overexpression as GFP-tagged fusion proteins inversely modifies cellular cholesterol levels. By providing strong evidence for disease-relevant functions of lipid trait-associated genes, our study demonstrates that quantitative, cell-based RNAi is a scalable strategy for a systematic, unbiased detection of functional effectors within GWAS loci.

Influenza A virus infection of human primary dendritic cells impairs their ability to cross-present antigen to CD8 T cells.

Influenza A virus (IAV) infection is normally controlled by adaptive immune responses initiated by dendritic cells (DCs). We investigated the consequences of IAV infection of human primary DCs on their ability to function as antigen-presenting cells. IAV was internalized by both myeloid DCs (mDCs) and plasmacytoid DCs but only mDCs supported viral replication. Although infected mDCs efficiently presented endogenous IAV antigens on MHC class II, this was not the case for presentation on MHC class I. Indeed, cross-presentation by uninfected cells of minute amounts of endocytosed, exogenous IAV was -300-fold more efficient than presentation of IAV antigens synthesized by infected cells and resulted in a statistically significant increase in expansion of IAV-specific CD8 T cells. Furthermore, IAV infection also impaired cross-presentation of other exogenous antigens, indicating that IAV infection broadly attenuates presentation on MHC class I molecules. Our results suggest that cross-presentation by uninfected mDCs is a preferred mechanism of antigen-presentation for the activation and expansion of CD8 T cells during IAV infection.

Entry of human papillomavirus type 16 by actin-dependent, clathrin- and lipid raft-independent endocytosis.

Infectious endocytosis of incoming human papillomavirus type 16 (HPV-16), the main etiological agent of cervical cancer, is poorly characterized in terms of cellular requirements and pathways. Conflicting reports attribute HPV-16 entry to clathrin-dependent and -independent mechanisms. To comprehensively describe the cell biological features of HPV-16 entry into human epithelial cells, we compared HPV-16 pseudovirion (PsV) infection in the context of cell perturbations (drug inhibition, siRNA silencing, overexpression of dominant mutants) to five other viruses (influenza A virus, Semliki Forest virus, simian virus 40, vesicular stomatitis virus, and vaccinia virus) with defined endocytic requirements. Our analysis included infection data, i.e. GFP expression after plasmid delivery by HPV-16 PsV, and endocytosis assays in combination with electron, immunofluorescence, and video microscopy. The results indicated that HPV-16 entry into HeLa and HaCaT cells was clathrin-, caveolin-, cholesterol- and dynamin-independent. The virus made use of a potentially novel ligand-induced endocytic pathway related to macropinocytosis. This pathway was distinct from classical macropinocytosis in regards to vesicle size, cholesterol-sensitivity, and GTPase requirements, but similar in respect to the need for tyrosine kinase signaling, actin dynamics, Na⁺/H⁺ exchangers, PAK-1 and PKC. After internalization the virus was transported to late endosomes and/or endolysosomes, and activated through exposure to low pH.

Characterization of drug-induced transcriptional modules: towards drug repositioning and functional understanding.

In pharmacology, it is crucial to understand the complex biological responses that drugs elicit in the human organism and how well they can be inferred from model organisms. We therefore identified a large set of drug-induced transcriptional modules from genome-wide microarray data of drug-treated human cell lines and rat liver, and first characterized their conservation. Over 70% of these modules were common for multiple cell lines and 15% were conserved between the human in vitro and the rat in vivo system. We then illustrate the utility of conserved and cell-type-specific drug-induced modules by predicting and experimentally validating (i) gene functions, e.g., 10 novel regulators of cellular cholesterol homeostasis and (ii) new mechanisms of action for existing drugs, thereby providing a starting point for drug repositioning, e.g., novel cell cycle inhibitors and new modulators of α-adrenergic receptor, peroxisome proliferator-activated receptor and estrogen receptor. Taken together, the identified modules reveal the conservation of transcriptional responses towards drugs across cell types and organisms, and improve our understanding of both the molecular basis of drug action and human biology.

A uniquely efficacious type of CFTR corrector with complementary mode of action.

Three distinct pharmacological corrector types (I, II, III) with different binding sites and additive behavior only partially rescue the F508del-cystic fibrosis transmembrane conductance regulator (CFTR) folding and trafficking defect observed in cystic fibrosis. We describe uniquely effective, macrocyclic CFTR correctors that were additive to the known corrector types, exerting a complementary "type IV" corrector mechanism. Macrocycles achieved wild-type-like folding efficiency of F508del-CFTR at the endoplasmic reticulum and normalized CFTR currents in reconstituted patient-derived bronchial epithelium. Using photo-activatable macrocycles, docking studies and site-directed mutagenesis a highly probable binding site and pose for type IV correctors was identified in a cavity between lasso helix-1 (Lh1) and transmembrane helix-1 of membrane spanning domain (MSD)-1, distinct from the known corrector binding sites. Since only F508del-CFTR fragments spanning from Lh1 until MSD2 responded to type IV correctors, these likely promote cotranslational assembly of Lh1, MSD1, and MSD2. Previously corrector-resistant CFTR folding mutants were also robustly rescued, suggesting substantial therapeutic potential for type IV correctors.

Immunocapture LC-MS(/MS) assays for biotherapeutic and biomarker proteins: the European Bioanalysis Forum continuing discussions on scientific and regulatory challenges.

The use of LC-MS(/MS) assays to quantify (biotherapeutic or biomarker) proteins is commonplace and well accepted across industry. There is a good understanding on the added value over conventional analytical technologies (i.e., ligand-binding assays). In fact, the impact of combining small- and large-molecule technologies for large-molecule analysis has played a significant part in bringing the bioanalytical communities closer together and building a mutual respect and understanding between scientists. This paper from the European Bioanalysis Forum presents a history of the journey and future perspectives for hybrid assays, with focus on the unanswered scientific questions, including regulatory discussions to be had. Hybrid assays are essentially a combination of ligand-binding assays and MS, and the ICH M10 guideline does not address this approach directly. Decision-based acceptance criteria are still being discussed, and the industry should continue to do so.

Quantitative proteomic landscapes of primary and recurrent glioblastoma reveal a protumorigeneic role for FBXO2-dependent glioma-microenvironment interactions.

BACKGROUND: Recent efforts have described the evolution of glioblastoma from initial diagnosis to post-treatment recurrence on a genomic and transcriptomic level. However, the evolution of the proteomic landscape is largely unknown. METHODS: Sequential window acquisition of all theoretical fragment ion spectra mass spectrometry (SWATH-MS) was used to characterize the quantitative proteomes of two independent cohorts of paired newly diagnosed and recurrent glioblastomas. Recurrence-associated proteins were validated using immunohistochemistry and further studied in human glioma cell lines, orthotopic xenograft models, and human organotypic brain slice cultures. External spatial transcriptomic, single-cell, and bulk RNA sequencing data were analyzed to gain mechanistic insights. RESULTS: Although overall proteomic changes were heterogeneous across patients, we identified BCAS1, INF2, and FBXO2 as consistently upregulated proteins at recurrence and validated these using immunohistochemistry. Knockout of FBXO2 in human glioma cells conferred a strong survival benefit in orthotopic xenograft mouse models and reduced invasive growth in organotypic brain slice cultures. In glioblastoma patient samples, FBXO2 expression was enriched in the tumor infiltration zone and FBXO2-positive cancer cells were associated with synaptic signaling processes. CONCLUSIONS: These findings demonstrate a potential role of FBXO2-dependent glioma-microenvironment interactions to promote tumor growth. Furthermore, the published datasets provide a valuable resource for further studies.

NF-κB drives epithelial-mesenchymal mechanisms of lung fibrosis in a translational lung cell model.

In the progression phase of idiopathic pulmonary fibrosis (IPF), the normal alveolar structure of the lung is lost and replaced by remodeled fibrotic tissue and by bronchiolized cystic airspaces. Although these are characteristic features of IPF, knowledge of specific interactions between these pathological processes is limited. Here, the interaction of lung epithelial and lung mesenchymal cells was investigated in a coculture model of human primary airway epithelial cells (EC) and lung fibroblasts (FB). Single-cell RNA sequencing revealed that the starting EC population was heterogenous and enriched for cells with a basal cell signature. Furthermore, fractions of the initial EC and FB populations adopted distinct pro-fibrotic cell differentiation states upon cocultivation, resembling specific cell populations that were previously identified in lungs of patients with IPF. Transcriptomic analysis revealed active NF-κB signaling early in the cocultured EC and FB, and the identified NF-κB expression signatures were found in "HAS1 High FB" and "PLIN2+ FB" populations from IPF patient lungs. Pharmacological blockade of NF-κB signaling attenuated specific phenotypic changes of EC and prevented FB-mediated interleukin-6, interleukin-8, and CXC chemokine ligand 6 cytokine secretion, as well as collagen α-1(I) chain and α-smooth muscle actin accumulation. Thus, we identified NF-κB as a potential mediator, linking epithelial pathobiology with fibrogenesis.

Functional and spatial proteomics profiling reveals intra- and intercellular signaling crosstalk in colorectal cancer.

Precision oncology approaches for patients with colorectal cancer (CRC) continue to lag behind other solid cancers. Functional precision oncology-a strategy that is based on perturbing primary tumor cells from cancer patients-could provide a road forward to personalize treatment. We extend this paradigm to measuring proteome activity landscapes by acquiring quantitative phosphoproteomic data from patient-derived organoids (PDOs). We show that kinase inhibitors induce inhibitor- and patient-specific off-target effects and pathway crosstalk. Reconstruction of the kinase networks revealed that the signaling rewiring is modestly affected by mutations. We show non-genetic heterogeneity of the PDOs and upregulation of stemness and differentiation genes by kinase inhibitors. Using imaging mass-cytometry-based profiling of the primary tumors, we characterize the tumor microenvironment (TME) and determine spatial heterocellular crosstalk and tumor-immune cell interactions. Collectively, we provide a framework for inferring tumor cell intrinsic signaling and external signaling from the TME to inform precision (immuno-) oncology in CRC.

Vascular functioning and the water balance of ripening kiwifruit (Actinidia chinensis) berries.

Indirect evidence suggests that water supply to fleshy fruits during the final stages of development occurs through the phloem, with the xylem providing little water, or acting as a pathway for water loss back to the plant. This inference was tested by examining the water balance and vascular functioning of ripening kiwifruit berries (Actinidia chinensis var. chinensis 'Hort16A') exhibiting a pre-harvest 'shrivel' disorder in California, and normal development in New Zealand. Dye labelling and mass balance experiments indicated that the xylem and phloem were both functional and contributed approximately equally to the fruit water supply during this stage of development. The modelled fruit water balance was dominated by transpiration, with net water loss under high vapour pressure deficit (D(a)) conditions in California, but a net gain under cooler New Zealand conditions. Direct measurement of pedicel sap flow under controlled conditions confirmed inward flows in both the phloem and xylem under conditions of both low and high D(a). Phloem flows were required for growth, with gradual recovery after a step increase in D(a). Xylem flows alone were unable to support growth, but did supply transpiration and were responsive to D(a)-induced pressure fluctuations. The results suggest that the shrivel disorder was a consequence of a high fruit transpiration rate, and that the perception of complete loss or reversal of inward xylem flows in ripening fruits should be re-examined.

Pairwise effects between lipid GWAS genes modulate lipid plasma levels and cellular uptake.

Complex traits are characterized by multiple genes and variants acting simultaneously on a phenotype. However, studying the contribution of individual pairs of genes to complex traits has been challenging since human genetics necessitates very large population sizes, while findings from model systems do not always translate to humans. Here, we combine genetics with combinatorial RNAi (coRNAi) to systematically test for pairwise additive effects (AEs) and genetic interactions (GIs) between 30 lipid genome-wide association studies (GWAS) genes. Gene-based burden tests from 240,970 exomes show that in carriers with truncating mutations in both, APOB and either PCSK9 or LPL ("human double knock-outs") plasma lipid levels change additively. Genetics and coRNAi identify overlapping AEs for 12 additional gene pairs. Overlapping GIs are observed for TOMM40/APOE with SORT1 and NCAN. Our study identifies distinct gene pairs that modulate plasma and cellular lipid levels primarily via AEs and nominates putative drug target pairs for improved lipid-lowering combination therapies.

Generation of a zebrafish SWATH-MS spectral library to quantify 10,000 proteins.

Sequential window acquisition of all theoretical mass spectra (SWATH-MS) requires a spectral library to extract quantitative measurements from the mass spectrometry data acquired in data-independent acquisition mode (DIA). Large combined spectral libraries containing SWATH assays have been generated for humans and several other organisms, but so far no publicly available library exists for measuring the proteome of zebrafish, a rapidly emerging model system in biomedical research. Here, we present a large zebrafish SWATH spectral library to measure the abundance of 104,185 proteotypic peptides from 10,405 proteins. The library includes proteins expressed in 9 different zebrafish tissues (brain, eye, heart, intestine, liver, muscle, ovary, spleen, and testis) and provides an important new resource to quantify 40% of the protein-coding zebrafish genes. We employ this resource to quantify the proteome across brain, muscle, and liver and characterize divergent expression levels of paralogous proteins in different tissues. Data are available via ProteomeXchange (PXD010876, PXD010869) and SWATHAtlas (PASS01237).

Quantitative Interactomics in Primary T Cells Provides a Rationale for Concomitant PD-1 and BTLA Coinhibitor Blockade in Cancer Immunotherapy.

Deciphering how TCR signals are modulated by coinhibitory receptors is of fundamental and clinical interest. Using quantitative interactomics, we define the composition and dynamics of the PD-1 and BTLA coinhibitory signalosomes in primary effector T cells and at the T cell-antigen-presenting cell interface. We also solve the existing controversy regarding the role of the SHP-1 and SHP-2 protein-tyrosine phosphatases in mediating PD-1 coinhibition. PD-1 predominantly recruits SHP-2, but when absent, it recruits SHP-1 and remains functional. In contrast, BTLA predominantly recruits SHP-1 and to a lesser extent SHP-2. By separately analyzing the PD-1-SHP-1 and PD-1-SHP-2 complexes, we show that both dampen the TCR and CD28 signaling pathways equally. Therefore, our study illustrates how comparison of coinhibitory receptor signaling via quantitative interactomics in primary T cells unveils their extent of redundancy and provides a rationale for designing combinations of blocking antibodies in cancer immunotherapy on the basis of undisputed modes of action.

High-throughput proteomic analysis of FFPE tissue samples facilitates tumor stratification.

Formalin-fixed, paraffin-embedded (FFPE), biobanked tissue samples offer an invaluable resource for clinical and biomarker research. Here, we developed a pressure cycling technology (PCT)-SWATH mass spectrometry workflow to analyze FFPE tissue proteomes and applied it to the stratification of prostate cancer (PCa) and diffuse large B-cell lymphoma (DLBCL) samples. We show that the proteome patterns of FFPE PCa tissue samples and their analogous fresh-frozen (FF) counterparts have a high degree of similarity and we confirmed multiple proteins consistently regulated in PCa tissues in an independent sample cohort. We further demonstrate temporal stability of proteome patterns from FFPE samples that were stored between 1 and 15 years in a biobank and show a high degree of the proteome pattern similarity between two types of histological regions in small FFPE samples, that is, punched tissue biopsies and thin tissue sections of micrometer thickness, despite the existence of a certain degree of biological variations. Applying the method to two independent DLBCL cohorts, we identified myeloperoxidase, a peroxidase enzyme, as a novel prognostic marker. In summary, this study presents a robust proteomic method to analyze bulk and biopsy FFPE tissues and reports the first systematic comparison of proteome maps generated from FFPE and FF samples. Our data demonstrate the practicality and superiority of FFPE over FF samples for proteome in biomarker discovery. Promising biomarker candidates for PCa and DLBCL have been discovered.

Similarities and Differences of Blood N-Glycoproteins in Five Solid Carcinomas at Localized Clinical Stage Analyzed by SWATH-MS.

Cancer is mostly incurable when diagnosed at a metastatic stage, making its early detection via blood proteins of immense clinical interest. Proteomic changes in tumor tissue may lead to changes detectable in the protein composition of circulating blood plasma. Using a proteomic workflow combining N-glycosite enrichment and SWATH mass spectrometry, we generate a data resource of 284 blood samples derived from patients with different types of localized-stage carcinomas and from matched controls. We observe whether the changes in the patient's plasma are specific to a particular carcinoma or represent a generic signature of proteins modified uniformly in a common, systemic response to many cancers. A quantitative comparison of the resulting N-glycosite profiles discovers that proteins related to blood platelets are common to several cancers (e.g., THBS1), whereas others are highly cancer-type specific. Available proteomics data, including a SWATH library to study N-glycoproteins, will facilitate follow-up biomarker research into early cancer detection.