- 1Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, United States
- 2Frederick National Laboratory for Cancer Research, Frederick, MD, United States
- 3Cancer Research Technology Program, Antibody Characterization Lab, Frederick National Laboratory for Cancer Research, Frederick, MD, United States
- 4Experimental Pathology Laboratory, Laboratory of Pathology, Center for Cancer Research, National Cancer Institute, National Institute of Health, Bethesda, MD, United States
- 5Veteran’s Administration (VA) Cooperative Studies Program, Veteran’s Administration (VA) Boston Healthcare System (151MAV), Jamaica Plain, MA, United States
- 6Department of Medicine, Boston University School of Medicine, Boston, MA, United States
- 7Pathology and Laboratory Medicine Service, Program, Veteran’s Administration (VA) Boston Healthcare System, Jamaica Plain, MA, United States
- 8Department of Pathology, Harvard Medical School, Boston, MA, United States
- 9Department of Medicine, Dana-Farber Cancer Institute, Boston, MA, United States
- 10Office of Cancer Clinical Proteomics Research, National Cancer Institute, Bethesda, MD, United States
Immunotherapies are revolutionizing cancer care, producing durable responses and potentially cures in a subset of patients. However, response rates are low for most tumors, grade 3/4 toxicities are not uncommon, and our current understanding of tumor immunobiology is incomplete. While hundreds of immunomodulatory proteins in the tumor microenvironment shape the anti-tumor response, few of them can be reliably quantified. To address this need, we developed a multiplex panel of targeted proteomic assays targeting 52 peptides representing 46 proteins using peptide immunoaffinity enrichment coupled to multiple reaction monitoring-mass spectrometry. We validated the assays in tissue and plasma matrices, where performance figures of merit showed over 3 orders of dynamic range and median inter-day CVs of 5.2% (tissue) and 21% (plasma). A feasibility study in clinical biospecimens showed detection of 48/52 peptides in frozen tissue and 38/52 peptides in plasma. The assays are publicly available as a resource for the research community.
1 Introduction
Immunotherapies, such as immune checkpoint inhibitors (1), therapeutic cancer vaccines (2), and CAR-T cell treatments (3–5), are revolutionizing cancer care. Substantial response rates are seen in subsets of patients with particularly responsive tumors (e.g., melanomas, hematologic malignancies), including durable responses and perhaps cures in some patients with advanced disease. However, for most solid tumors, response rates remain <15%, and we do not have sufficient predictive biomarkers to identify patients whose tumors are likely to respond (6). Additionally, a significant number of patients receiving immunotherapy experience immune-related adverse events (7) (irAEs). Many irAEs are manageable with systemic immunosuppression, but some can be life-threatening (5, 8) (e.g., encephalitis, fulminant myocarditis) or lead to treatment discontinuation. Thus, there is an urgent need to understand mechanisms of response and resistance to immunotherapies to design more efficacious and less toxic immunotherapies, to identify biomarkers to select patients for single agent vs. combination immunotherapies, and to develop biomarkers to predict and monitor irAEs.
Hundreds of immunomodulatory proteins in the tumor microenvironment sculpt the T cell response to cancer as part of the “cancer-immunity cycle,” (9) and it is critical that we be able to quantify these proteins in clinical and translational research settings to design and deliver improved immunotherapies. Protein expression is typically quantified by immunoassay methods [e.g., immunohistochemistry (10), MIBI-TOF (11), flow cytometry (12)] that depend upon antibodies that are often not monospecific (13). As a result, assay interferences are commonly encountered in complex biospecimens, compromising assay specificity, and limiting multiplexing of protein assays.
To circumvent the immunoassay limitations, analysis of RNA transcripts has been used as a surrogate for protein measurements (14); however, mRNA expression levels are not reliable predictors of the expression level of most proteins (15–18), nor do they correlate with protein activity (e.g., post-translational modifications). Furthermore, many technologies do not quantify soluble proteins in the tumor microenvironment, which can impact tumor immunity (19). Thus, we need new protein quantification technologies that complement immunoassays for quantifying human proteins to enable new scientific discoveries and medical insights.
Liquid chromatography-multiple reaction monitoring mass spectrometry (LC-MRM-MS) is an emerging protein quantification method (20) in which peptides released via proteolysis are quantified as stoichiometric surrogates for proteins (21, 22). In contrast to untargeted “shotgun” MS profiling-based proteomics, targeted proteomics focuses the full analytic capacity of the mass spectrometer on pre-selected peptides (and the proteins they represent) of interest. Coupling an immunoaffinity enrichment step with MRM produces immuno-MRM assays that can precisely quantify low abundance proteins (23, 24) and posttranslational modifications (25, 26). Furthermore, because the mass spectrometer is used as the detector, interferences can be readily identified and usually avoided. As a result, MRM-based assays are readily multiplexed (27, 28), and the antibodies developed for immuno-MRM need not be monospecific.
Through the incorporation of stable isotope labeled internal standards, MRM assays can be harmonized across laboratories (29, 30), even on an international stage (31). Immuno-MRM assays have been applied to make clinically relevant measurements of proteins in human cancer tissues and fluids (32), including quantifying thyroglobulin in plasma where conventional immunoassays suffer interferences (33), quantification of cardiovascular health markers in plasma (34, 35), identifying novel pharmacodynamic biomarkers (36), multiplexing quantification of inborn errors of metabolism in dried blood spots (37–39), and quantifying HER2 in tissue and bone biopsies from breast cancer patients (40–43).
In this report, we present the development and characterization of a multiplexed panel (“IO-1 panel”) of immuno-MRM assays designed to quantify immunomodulatory proteins in human tissue biopsies and biofluids. The assays target 52 peptides (46 proteins) and are part of a larger effort (44) under the Beau Biden National Cancer Moonshot (45) to accelerate scientific discovery in cancer, foster greater collaboration, and improve the sharing of data.
Fit-for-purpose bioanalytical validation was conducted for the IO-1 assay panel in tumor tissue and plasma matrices to determine performance figures of merit. The performance of the assay panel was subsequently characterized in 135 tissue biospecimens (collected from 12 different tumor types) and 45 plasma biospecimens from cancer patients. The assay panel showed robust analytical performance and the targeted peptides were widely detected in the biospecimens. Additionally, the monoclonal antibodies generated in this project were tested for use in Western blotting and protein array, and all characterization data and antibody reagents are publicly available as resources for the research community through the National Cancer Institute’s Clinical Proteomic Tumor Analysis Consortium (CPTAC) Assay Portal (46, 47) (assays.cancer.gov) and Antibody Portal (antibodies.cancer.gov).
2 Methods
2.1 Materials and Reagents
Urea (#U0631), Trizma base (#T2694), citric acid (#C0706), dimethyl sulfoxide (DMSO, #D2438), EDTA (#E7889), EGTA (#E0396), and iodoacetamide (IAM, #A3221) were obtained from Sigma (St. Louis, MO). Acetonitrile (MeCN, #A955), water (#W6, LCMS Optima® grade), trifluoroacetic acid (TFA, LC-MS grade, #85183), tris(2-carboxyethyl)phosphine (TCEP, #77720), phosphate buffered saline (PBS, #BP-399-20), ammonium bicarbonate (A643-500), xylene (#422685000), and (3-[(3-cholamidopropyl) dimethylammonio]-1-propanesulfonate) (CHAPS, #28300) detergent were obtained from Thermo Fisher Scientific (Waltham, MA). Rapigest (#186001861) was obtained from Waters (Milford, MA). Formic acid (#1.11670.1000) was obtained from EMD Millipore (Billerica, MA). Lys-C (Wako, #129-02541) and sequencing grade trypsin (#V5111, Promega, Madison, WI) were used for digestion of samples. Rabbit monoclonal antibodies were produced with Epitomics/Abcam (Cambridge, MA) and Excel Biopharm (Burlingame, CA). Mouse monoclonal antibodies were produced with Precision Antibody (Columbia, MD) and the Antibody Development Facility at the Fred Hutchinson Cancer Research Center (Seattle, WA). Light (unlabeled) synthetic peptides were obtained from Vivitide (Gardner, MA) as crude (flash purified) grade. Cleavable stable isotope-labeled (heavy) peptides from Vivitide corresponding to the tryptic analyte sequence were purified >95% by HPLC, labeled with [13C and 15N] at the tryptic C-terminal Arg or Lys, and quantified by amino acid analysis (AAA). Aliquots of peptide standards were stored in 5% acetonitrile/0.1% formic acid at −80°C until use.
2.2 Cell Lines, Culture Conditions, and Cell Lysis
HeLa (American Type Culture Collection (ATCC), #CCL-2), Jurkat (ATCC, #TIB-152), A549 (ATCC, #CCL-185), MCF7 (ATCC, #HTB-22), and NCI-H226 (ATCC, # CRL-5826) cell lines were cultured and harvested according to manufacturers’ specifications. Briefly, cells were sub-cultured and after trypsinization, were centrifuged at 1500 rpm for 6 minutes and supernatant was removed and discarded. Cells were washed once with 1X PBS, centrifuged as before, resuspended in 1X PBS, and counted. Cells were subdivided according to desired cell number in 15 mL centrifuge tubes and centrifuged as before. Cell pellets were frozen and stored at -80°C until whole cell lysis. Cell pellets were lysed using RIPA lysis and extraction buffer (Thermo Fisher Scientific, # 89900) following the manufacturer’s protocol. Mammalian Protease Inhibitor (VWR, #VWRVM250) was added according to manufacturer’s instructions. Lysate protein concentration was estimated with a BCA assay (Pierce, #23225) as described by the manufacturer’s instructions.
2.3 Human Samples
FFPE tissue and plasma samples used for fit-for-purpose validation studies were commercially acquired from BioIVT (Westbury, NY). Frozen tissue and plasma samples for determination of detectability were supplied by the Clinical Proteomics Tumor Analysis Consortium (CPTAC) as anonymized samples from consenting donors collected under IRB-approved protocols. FFPE tissue was supplied by the Veteran’s Administration from consenting donors collected under IRB-approved protocols. Tissue sub-compartment cellularity (e.g. tumor, stroma, adipocytes, lymphocytes) was calculated using the HALO Tissue Classifier (Indica Labs, NM) using a machine learning algorithm that was trained to apportion the total area of a section into the following components: epithelium, stroma, adipocytes, lymphocytes, and red blood cells; the training was based on results from two pathologists using blinded stained slides from breast cancer specimens (40).
2.4 Protein Extraction From Frozen Tissue and FFPE Samples
To produce lysates for immuno-MRM analysis, frozen tissue was cryopulverized in a cryoPREP CP-02 (Covaris, Woburn, MA) and stored frozen until analysis. 5 μL of lysis buffer (25 mM Tris, 6 M Urea, 1 mM EDTA, 1 mM EGTA, 1% (v/v) Sigma protease inhibitor (#P8340), 1% (v/v) Sigma phosphatase inhibitor cocktail 2 (#P5726), 1% (v/v) Sigma phosphatase inhibitor cocktail 3 (#P0044)) was added for each mg wet tissue weight (up to 1000 μL). The sample was vortexed for 10-15 sec and sonicated three times in a cup horn probe (filled with ice water) at 50% power for 30 seconds. The samples were stored in liquid nitrogen until the day of digestion.
FFPE samples were processed as described previously (48). Briefly, slide-mounted FFPE tissue sections were placed in a 4 slide holder. The slides were incubated three times in xylene for 3 min followed by 100% (v/v) ethanol twice for 3 min. The tissue was then hydrated twice in 85% (v/v) ethanol for 3 min, 70% (v/v) ethanol for 3 min, and distilled water for 3 min. The tissue was then blotted and scraped off the slide into a screw cap microfuge tube. To each sample tube (containing three FFPE 10 µm tissue sections), extraction buffer (0.2% RapiGest in 50 mM ammonium bicarbonate, NH4HCO3) was added and incubated at 95°C for 30 minutes with mixing at 1000 rpm (Thermomixer, Eppendorf, Enfield, CT). The samples were then cooled on ice for 5 minutes and sonicated twice in a cup horn probe (filled with ice water) at 50% power for 30 sec. The samples were then incubated at 80°C for 120 minutes with mixing at 1000 rpm and then cooled on ice for 5 min. 100 μL of 50 mM NH4HCO3, pH 8.0 was added, and the samples were sonicated twice in the cup horn probe (filled with ice water) at 50% power for 30 sec. Following processing, all samples were stored at -80°C until the day of digestion.
Protein concentrations of lysates were measured in triplicate using Micro BCA Protein Assay Kit (Pierce, #23235).
2.5 Protein Digestion
A mix of cleavable stable isotope-labeled peptide standards was added to the lysate at 200 fmol/capture. 500 μg frozen tissue lysates or 15-270 μg of protein (FFPE) was transferred to a deep-well plate for processing on an EpMotion 5057 (Eppendorf). Lysates were reduced in 30 mM TCEP for 30 minutes at 37°C with shaking, followed by alkylation with 50 mM IAM at room temperature. Lysates were then diluted with 0.8 mL 200 mM TRIS before Lys-C endopeptidase was added at a 1:50 enzyme:substrate ratio by mass and incubated for 2 hours at 37°C with mixing at 600 rpm (Thermomixer, EpMotion 5057). After 2 hours, trypsin was added at a 1:50 enzyme:substrate ratio. Digestion was carried out overnight at 37°C with mixing at 600 rpm. After 16 hours, the reaction was quenched with formic acid (final concentration 1% by volume).
Plasma (2 x 50 μL aliquots of each sample) was denatured with 150 μL of lysis buffer (25 mM Tris, 6 M Urea, 1 mM EDTA, 1 mM EGTA, 1% (v/v) Sigma protease inhibitor (#P8340), 1% (v/v) Sigma phosphatase inhibitor cocktail 2 (#P5726), 1% (v/v) Sigma phosphatase inhibitor cocktail 3 (#P0044)). The plasma was reduced with 30 mM TCEP at 37°C for 30 min, and alkylated with 50 mM IAM at room temperature for 30 min. Urea concentration was diluted 10-fold with 200 mM Tris prior to overnight digestion at 37°C with Lys-C/trypsin using a 1:50 (w/w) enzyme:substrate. Digestions were terminated with formic acid.
The mixture was desalted using Oasis HLB 96-well plates (Waters #WAT058951) and a positive pressure manifold (Waters #186005521) according to the following procedure: wash cartridge with 4 × 400 μL of 50% acetonitrile in 0.1% formic acid, equilibrate with 4 × 400 μL of 0.1% formic acid, load total volume of digest, wash with 4 × 400 μL of 0.1% formic acid, and elute with 3 × 400 μL of 50% acetonitrile in 0.1% formic acid. The eluates were lyophilized and stored at -80°C.
2.6 Peptide Immunoaffinity Enrichment
Enrichment was performed as previously described (49) with the following modifications. The final assay consisted of a mixture of 50 antibodies. Antibodies were crosslinked on protein G beads (GE Sepharose, #28-9513-79), and peptide enrichment was performed using 1 μg antibody - protein G magnetic beads for each target. Trypsin-digested lysate was resuspended in 200 μL 1X PBS + 0.01% CHAPS (pH was adjusted to 7.0 with 10 μL of 1 M Tris, pH 9). For plasma, the two digestion aliquots were combined after resuspension to a total volume of 200 μL. Beads were mixed in the incubation plate, washed twice in 1X PBS buffer + 0.01% CHAPS, washed once in 1/10X PBS + 0.01% CHAPS, and peptides were eluted in 26 μL of 5% acetic acid/3% acetonitrile/50 mM citrate. The elution plate was covered with adhesive foil and frozen at -80°C until analysis.
2.7 Liquid Chromatography Multiple Reaction Monitoring-Mass Spectrometry
LC-MS was performed with an Eksigent 425 nanoLC system with a nano autosampler and chipFLEX system (Eksigent Technologies, Dublin, CA) coupled to a 5500 QTRAP mass spectrometer (SCIEX, Foster City, CA). Peptides were loaded on a trap chip column (Reprosil C18-AQ, 0.5 mm x 200 μm, SCIEX, #804-00016) at 5 μL/min for 3 minutes using mobile phase A (0.1% formic acid in water). The LC gradient was delivered at 300 nL/minute and consisted of a linear gradient of mobile phase B (90% acetonitrile and 0.1% formic acid in water) developed from 2-14% B in 1 minute, 14-34% B in 20 minutes, 34-90% B in 2 minutes, and re-equilibration at 2% B on a 15 cm x 75 μm chip column (ChromXP 3C18-CL particles, 3 μm, SCIEX, #804-00001). The nano electrospray interface was operated in the positive ion MRM mode. Parameters for declustering potential (DP) and collision energy (CE) were taken from optimized values in Skyline. Scheduled MRM transitions used a retention time window of 210 seconds and a desired cycle time of 1.5 seconds, enabling sufficient points across a peak for quantification. A minimum of two transitions per peptide, including endogenous and spiked heavy peptides, were recorded for each light and heavy peptide.
2.8 Data Analysis
MRM data acquired on the 5500 QTRAP were analyzed by Skyline (50, 51). Peak integrations were reviewed manually, and transitions from analyte peptides were confirmed by the same retention times and relative transition areas of the light peptides and heavy stable isotope-labeled peptides. Transitions with detected interferences were not used in the data analysis. Integrated raw peak areas were exported from Skyline and total intensity was calculated using Peak Area + Background. Transitions were summed for each light/heavy pair and peak area ratios were obtained by dividing peak areas of light peptides by that of the corresponding heavy peptides (or vice versa for response curves). All measurements were filtered by the LLOQ (i.e., all measurements were required to be above the LLOQ). Peak area ratios were log (base 2) transformed for statistical analysis.
Quantitative RNA sequencing data (FPKM) was available for 74/110 frozen tissues (52). MRM peak area ratios were log (2) transformed and correlated with gene expression data using Pearson regression. Correlations were done within each tumor type and across all tumor types (n = 74).
Minimum tissue requirements were calculated by using frozen tissue results of signal-to-noise levels measured at 500 μg input. The signal-to-noise levels were compared to the LOQ for each analyte using a one-sample t-test and those analytes with 95% confidence above LOQ were further converted to input mass using linear scale dilution. The confidence interval in the minimum tissue calculation was determined based on the standard deviation of signal-to-noise distribution and degrees of freedom (equal to the number of samples for each tissue site -1).
2.9 Fit-For-Purpose Assay Validation
Four experiments (described below) were performed to characterize the analytical performance of the assays: i) digestion time course for stable recovery of peptides, ii) response curves, iii) repeatability, and iv) stability.
2.9.1 Digestion Process Time Course
Recovery of peptides was measured by a time course study using triplicate process replicates of 150 μg of protein of a pooled lysate from cell lines processed according to the trypsin digestion workflow described above. Aliquots were prepared for quenching of the digestion at four time points (2, 6, 16, and 24 hours). Digested samples were enriched and analyzed as described above.
2.9.2 Response Curve
Response curves were generated in a background matrix consisting of an equal mixture of protein lysate from 5 commercially obtained (FFPE) lung tumors, an equal mixture of lysates from 4 frozen tumors, or commercially obtained plasma. Aliquots (100 μg of tissue, 10 μL plasma) of the matrix were spiked with cleavable heavy stable isotope-labeled peptides covering the concentrations 200, 20, 2, 0.8, 0.32, 0.128, 0.0512, and 0.02048 fmol/uL for plasma lysates and 20,000, 2,000, 80, 32, 12.8, 5.12, and 2.048 fmol/mg for FFPE and frozen tissue lysates. Concentrations for spikes for heavy extended peptides corresponding to ADAM17.VDNE, CD163.LVDG, GAPDH.GALQ, TNFRSF9.NQIC, and VSIR.GHDV were ten times higher. Light peptide was also spiked into the digested lysate at 20 fmol/uL for plasma and 2,000 fmol/mg for FFPE and frozen lysates. Blanks were prepared using background matrix with light peptide (no heavy spike). All points were analyzed in four replicates. Curves were analyzed using Skyline by performing linear regression using log transformed data on all points above the lower limit of quantification. The Lower Limit of Quantifications (LLOQs) were obtained by empirically finding the lowest point on the curve that had CV <20% in the curve replicates. The upper limit of quantification (ULOQ) was determined by the highest concentration point of the response curve that was maintained in the linear range of the response. For curves that maintained linearity at the highest concentration measured, the ULOQ is a minimum estimate.
2.9.3 Repeatability
Repeatability was determined using the same matrix used to generate the response curves except 50 mM ammonium bicarbonate was used in place or Tris for the protein digestion. Cleavable heavy peptide standards were spiked at three levels (low, medium, and high) at 1.33e2, 1.33e3, 1.33e4 fmol/mg for FFPE samples and 2, 20, 200 fmol/uL for plasma sample. Spikes for heavy extended peptides corresponding to ADAM17.VDNE, CD163.LVDG, GAPDH.GALQ, TNFRSF9.NQIC, and VSIR.GHDV were ten times higher at 1.33e3, 1.33e4, 1.33e5 fmol/mg for FFPE samples and 20, 200, 2,00 fmol/uL for plasma samples. Analogous to response curve experiments, following digestion, light peptides were also spiked at 20 fmol/uL for plasma and 1,333.3 fmol/mg for FFPE lysates. Complete process triplicates (including digestion, capture, and mass spectrometry) were prepared and analyzed on five independent days. Intra-assay variation was calculated as the mean CV obtained within each day. Inter-assay variation was the CV calculated from the mean values of the five days.
2.9.4 Peptide Stability
Stability of the enriched peptides was determined by analyzing aliquots of the medium spike level sample used in repeatability studies after storage at 4°C in the autosampler for approximately 8 hours, storage at 4°C for 48 hours, after 2 freeze-thaw cycles, and after extended storage at -80°C for 5 weeks. Each condition was measured in process triplicate.
2.10 Characterization of Antibodies in Immunoblotting and Protein Array
Recombinant proteins and over-expressed lysates were obtained from Origene (Rockville, MD) or Novus Biologicals (Littleton, CO), see Supplementary Table 1 for catalog numbers. Western immunoassays were performed using traditional immunoblotting techniques according to the designed protocol described in the CPTAC Antibody Portal (antibodies.cancer.gov) SOP M-103 with the following modifications. Proper sample loading on the 4-20% Criterion TGX Stain-free precast gels (Bio-Rad, # 5678094) was verified by rapid florescence detection with ChemiDoc MP imager. Traditional immunoblotting used 10 μg/mL recombinant protein under reducing conditions (20 mL, 200 ng total protein/lane). Whole cell lysates were diluted to 2.5 mg/mL in reducing conditions (20 mL, 50 μg total protein/lane). Transfers of protein from Bio-Rad precast gels were performed by Bio-Rad Turbo-Blot at “High MW” setting for 10 minutes. Blocking of the membrane was performed using Bio-Rad Blotting Grade Blocker (Bio-Rad, #1706404) at 5% in 1X PBS/0.5% Tween-20. Primary antibodies (1 mg/mL) were diluted in 1X PBS/0.5% Tween-20 to a dilution of 1:5000 at a total volume of 25 mL. Washing of membrane was conducted using 1X PBS/0.5% Tween-20 three times. Secondary HRP-linked rabbit specific antibody (Jackson ImmunoResearch, West Grove, PA, 111-035-144) or secondary HRP-linked mouse antibody (Jackson Immunoresearch Laboratories 115-035-062) was diluted at 1:5000 in 1 X PBS/0.5% Tween-20 at a final volume of 25 mL. Immuno-detection was performed using colorimetric substrate Opti-4-CN (Bio-Rad, #1708235) at 1 mL per blot or using enhanced chemiluminescence using Clarity Western ECL Substrate (BioRad, #1705061) at 1 mL per blot. Development of immunoblot was captured using Bio-Rad ChemiDoc MP imaging system.
The Simple Western (Wes, ProteinSimple, San Jose, CA, #004-600) system was used to detect primary antibody binding to a target protein in cell lysates (MCF10A, LCL57). The Simple Westerns were performed following the procedures detailed in the CPTAC Antibody Portal (antibodies.cancer.gov) SOP M-134 with the following modifications. Cell lysates were run using the 12-230 kDa separation module, 8 x 25 capillary cartridges (ProteinSimple, #SM-W004), and detected with Anti-Rabbit Detection Module (ProteinSimple, #DM-001) or Anti-mouse Detection Module (ProteinSimple, #DM-002). Cell lysates were run at a concentration of 200 μg/mL and incubated with primary antibodies diluted 1:500.
Protein array analyses were performed using NCI-60 cell lines obtained from the Cancer Research Technology Program at NCI-Frederick (Frederick, MD). The NCI-60 cell lines were collected at the log phase growth and protein prepared by resuspending cell pellets in RIPA (Thermo Fisher Scientific, #89900) per the manufacturer’s recommendations; total protein concentration was measured by BCA Protein Assay kit (Thermo Fisher Scientific, #23225). Quantification of protein expression values was performed by well-based reverse phase protein array (RPPA) as previously reported (53, 54). Briefly, five microliters (100 ng/well) of NCI-60 cell line antigens in PBST (1X PBS, 0.1% Tween-20) were applied onto 96-well Multi-Array™ plates (96 HB SECTOR Plate, Meso Scale Discovery, Gaithersburg, MD). The plates were allowed to dry at room temperature for 2 hours. Prior to primary antibody incubation, the antigen-coated plates were blocked with 5% non-fat dry milk in PBST for 1 hour at room temperature. Target-specific antibodies were diluted (1:1000 to 1:5000) with 5% BSA in PBST. For each cell line, 25 μL of antibody were added and incubated overnight at 4 °C. The plates were washed with PBST and followed by a 90 min incubation with goat anti-mouse or anti-rabbit SULFO-TAG™ antibodies (Meso Scale Discovery) at a dilution of 1:2000 (0.5 µg/mL) containing 5% non-fat dry milk in PBST. Plates were washed, and MSD-T read buffer was added to the plate to detect binding signals using MESO QuickPlex SQ 120 reader (Meso Scale Discovery). PBST-coated wells were included on each plate as a control of non-specific binding. For each mAb, the electrochemical luminescence value of each cell line was normalized by the average value of the 60 cell lines. After normalization, levels below 0.5 were interpreted as weak or negative expression and levels above 1.5 were interpreted as strong, positive expression.
2.11 Public Availability of Data and Reagents
Targeted mass spectrometry data are available in the supplementary tables and Panorama Public (55), a public repository of targeted proteomics experiments (https://panoramaweb.org/IO1immunoMRM.url. Characterization data for assays can be found via the CPTAC Assay Portal (https://assays.cancer.gov) and antibodies are available through the CPTAC Antibody Portal (antibodies.cancer.gov) (see Supplementary Table 2 for IDs).
3 Results
3.1 Development of a Targeted, Multiplexed IO-1 Immuno-MRM Assay
Because there are no adequate predictive or prognostic biomarkers to guide delivery of immunotherapies, there is great interest in high quality assays to quantify immunomodulatory proteins for correlative studies in clinical trials in hopes of identifying single or panels of biomarkers showing clinical utility (56). Proteins and phosphorylation sites related to immunomodulatory functions were identified as targets for assay development through literature searching and expert consultation. A panel of experts in immuno-oncology was convened from academia and industry to nominate immunomodulatory proteins of interest for assay development (Table 1) to support correlative studies to identify novel biomarkers or response or immune-related adverse events. Stoichiometric surrogate peptides for the target proteins were identified using established procedures (88, 89), by mining liquid chromatography-tandem mass spectrometry (LC-MS/MS) proteomic and phosphoproteomic datasets (16, 17, 90–98) and public databases (52, 99, 100) for empirical evidence of LC-MS/MS detectability. Peptides were ranked based on the intensity and frequency of empirical observations in the datasets, chemical properties (e.g., amino acid composition, hydrophobicity), and physical properties (e.g., length, position in protein). Finally, peptides with frequent mutation sites [contained in dbSNP (101)] and PTMs [PhosphoSitePlus (102)] were avoided.
A list of 55 peptides (including multiple peptides per protein) were selected for assay development, representing 46 proteins and including 5 phosphorylation sites. Peptide sequences selected for assay development are listed in Table 1. Unlabeled (i.e., light) tryptic peptides were synthesized for assay development, and cleavable heavy stable isotope-labeled peptides were synthesized for use as internal standards. The cleavable standards incorporate 2-5 additional amino acids (corresponding to the native protein sequence) on either side of the tryptic cut sites. This strategy was adopted in order to spike the standards into the sample prior to digestion and provide a within-sample control for trypsin digestion (20). The synthetic light peptides were used to optimize the selection of specific precursor and fragment ion pairs (i.e., transitions), determining the chromatographic retention time of each peptide, and optimizing collision energy parameters in the mass spectrometer (Supplementary Table 3).
Immuno-MRM uses anti-peptide antibodies (103) for enrichment prior to LC-MRM and is applicable for quantifying expression of high and low abundance proteins as well as phosphopeptides (Figure 1A). For the immuno-MRM assay panel, we generated custom anti-peptide monoclonal antibodies to the linear peptide sequences using established approaches (25, 26). We used 50 custom monoclonal anti-peptide antibodies (Supplementary Table 2) for the multiplexed immuno-MRM assay, including 35 antibodies generated for this study and 15 antibodies from previously characterized assays (26, 40).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/765898/fimmu-12-765898-HTML-r1/image_m/fimmu-12-765898-g001.jpg)
Figure 1 Immuno-MRM enables highly multiplexed protein quantification. (A) The immuno-MRM assay workflow commences with generation of a protein lysate from the biospecimen of interest. Cleavable stable isotope labeled standards unique to each targeted peptide sequence are spiked into the sample at a known concentration. The protein mixture is converted to peptides by enzymatic digestion (Lys-C and trypsin). Custom monoclonal antibodies coupled to magnetic beads are used to enrich the endogenous peptides and labeled standards. The eluate is analyzed by multiple reaction monitoring-mass spectrometry, where analyte peptides and internal standards coelute with equivalent relative areas of monitored transitions. High sensitivity is achieved through analyte enrichment and optimization of mass spectrometer parameters for the targeted peptides. High specificity is maintained through optimal selection of fragment ion transitions. (B) Trypsin-mediated release of peptides was produced by overnight enzymatic digestion of proteins from a pool of cell lysates. The peak area ratios (light:heavy) for the 33 endogenous peptides detected were normalized to the maximum timepoint and plotted over time. Error bars are the standard deviation of three replicates. (C) Performance figures of merit for assays characterized in tissue matrix. (D) Performance figures of merit for assays characterized in plasma matrix. A representative response curve for the peptide VEIIATMK from CXCL10. Each concentration point was measured in triplicate. Distribution of R2 values from linear regression of the response curves. Distribution of lower limit of quantification (LOQ) where each point refers to concentration determined by the lowest point on the curve with less than 20% CV. Repeatability is characterized by the distribution of CV values for intra- (within day) and inter- (between day) variability at three concentrations, in addition to the measurement of endogenous peptide. Each point corresponds to the average %CV for a peptide measured at three concentrations, Low (Lo), Medium (Med), and High (Hi) in triplicate over five days (n=15 at each concentration for a peptide). Endogenous measurements refer to the intra (within day) and inter- (between day) variability of endogenous peptides detected above the LOQ in five replicates measured over 5 days (n = 25). Stability shows the distribution of %CV and relative percent difference for 4 conditions compared to immediate analysis: (i) stored at 8 hours at 4°C, (ii) 48 hours at 4°C, (iii) after two freeze-thaw cycles, and (iv) stored at -80°C for 5 weeks. For box plots, the line shows the median value, boxes show the inner quartiles, and the whiskers show 5-95% of data.
3.2 Fit-For-Purpose Method Validation of the IO-1 Immuno-MRM Assay
Trypsin digestion has been identified as a source of variability in sample processing for targeted proteomics (29, 104). Therefore, to achieve reproducible digestion, we sought to identify a timepoint where release of the tryptic peptides (light and heavy) was most stable and at a maximum. Recovery of peptides following trypsin digestion was characterized using a time course digestion study using an optimized procedure incorporating Lys-C and trypsin for digestion (105). Aliquots of pooled cell line lysates were spiked with heavy peptides and digested according to the methods described herein. To evaluate the release of endogenous tryptic peptides over time, aliquots were taken after 1, 6, 16, and 24 hours of incubation with trypsin and subjected to immunoaffinity enrichment and analysis by LC-MRM. The recovery of light peptides reached a stable maximum at the 16-hour timepoint, with two peptides showing slightly lower (~20%) signals (Figure 1B). To achieve the maximum recovery for most peptides in the assay, we used overnight (i.e., 16 hours) digestion for processing samples.
Performance of the multiplexed IO-1 immuno-MRM assay panel was characterized in tissue and plasma matrices according to published guidelines (47, 106) to establish figures of merit including linearity, limits of quantification, repeatability, and stability. The linear ranges and limits of quantification (LOQs) were determined by response curves using pooled background matrices of protein lysates from i) frozen lung tumor tissues, ii) formalin fixed and paraffin embedded (FFPE) lung tumor tissues, and iii) plasma. Aliquots (150 µg) of the pooled tissue lysates or 10 µL aliquots of plasma were spiked with synthetic cleavable heavy stable isotope-labeled standards and serially diluted to yield concentrations of 2000, 200, 20, 8, 3.2, 1.28, 0.512, 0.205, and 0 (blank) fmol/mg or fmol/µL, respectively, and each concentration level was digested in process triplicate. Synthetic light tryptic peptides were added at a constant concentration of 200 fmol/mg or 200 fmol/µL prior to desalting the digested peptides. The monoclonal antibodies were coupled to magnetic beads and used to enrich the peptides. The eluates were analyzed by LC-MRM-MS. Peptide specificity was confirmed by equivalent retention times and relative transition areas of the heavy and light peptides.
Peak area ratios (heavy:light) were plotted as a function of analyte concentration to determine the assay figures of merit in the response curves. Figures 1C, D show representative response curves measured in tissue and plasma. Lower limits of quantification (LLOQs) were determined by the lowest point with CV < 20%. Linear ranges were determined by using points on the linear regression with correlation coefficients greater than 0.9. For assays where the highest concentration point was still linear, the linear range was a minimum estimate. Figures of merit are reported for each peptide in Supplementary Table 4. Median linear dynamic range was ≥3.2 orders of magnitude in all matrices with median LLOQ 12.8 fmol/mg (range 2-2000 fmol/mg) in tissues and 0.128 fmol/µL (range 0.02-2000 fmol/µL) in plasma. The characterized LLOQs and linear ranges were the same for frozen and FFPE tissues (Supplementary Figure 1). Three peptides (LIME1.SSTC[cam]GAGT[ph]PPASSC[cam]PSLGR, LIME1.SSTC[cam]GAGTPPASSC[cam]PSLGR, and VSIR.GHDVTFYK) exhibited low signal in the curves, likely due to poor antibody activity when cross-linked to the beads, and failed to produce adequate linearity for further use.
Intra-assay (within day) and inter-assay (between day) repeatability were determined by performing complete process triplicate measurements for the multiplexed assay at three concentrations of spiked peptides over 5 days. Heavy cleavable peptide standards were spiked into 150 µg aliquots of the pooled FFPE tissue lysate or 10 µL aliquots of plasma matrix at three concentrations (10, 100, 1000 fmol/mg; low, medium, high) with addition of equal amounts of light peptides (200 fmol/mg) to each aliquot. Specificity was confirmed using the same criteria as described above. For the assays measured in tissue matrix, the median intra-assay variability was 7.5%, 5.4%, 3.5% for low to high concentration samples and the median inter-assay variability was 9.5%, 10.8% and 11.8% (low to high) (Figure 1C and Supplementary Table 5). Two outlier assays (ATM.NLS[ph]DIDQSFNK and ATM.NLSDIDQSFNK) showed high intra- or inter-assay variability in tissue samples due to signal intensities below the LOQ at the lowest concentration values. For the assays measured in plasma matrix, the median intra-assay variability was 13.6%, 5.8%, 4.5% from low to high concentration samples and the median inter-assay variability was 24.9%, 14.9% and 16.5% from low to high concentration samples (Figure 1D and Supplementary Table 5). Three assays (ANXA1.AAYLQETGKPLDETLK, PTPRC.LFLAEFQSIPR, PSMA1.ETLPAEQDLTTK) showed high variability when applied to plasma. High variability was likely due to different factors for the peptides. ANXA1.AAYL had high variability in the heavy signal, PTPRC.LFLA was likely due to instability of the peptide following processing (see below), and PSMA1.ETLP had one day with outlier values compared to other days.
To characterize expected variability in measuring endogenous proteins, heavy cleavable peptide standards were spiked into 150 µg aliquots of a pooled FFPE lung tissue lysate and 10 µL aliquots of plasma. Measurements were made using 5 complete process replicates of endogenous peptide over 5 days (n=25) (Figures 1C, D). Forty-one peptides were detected above LLOQ in the endogenous tissue sample. The median intra-assay variability for endogenous detection was 5.2% (range 2.4-16.9%) and the median inter-assay variability was 6.9% (range 0.8-21.4%). Twenty peptides were detected above LOQ in the endogenous plasma sample. The median intra-assay variability was 12.2% (range 4.1-34%) and the median inter-assay variability was 21% (range 4.5-62%).
Peptide stability was evaluated by spiking heavy peptide (200 fmol/mg) into 150 µg aliquots of the pooled FFPE tissue lysate or 200 fmol/µL into 10 µL aliquots of plasma matrix, followed by digestion and immunoaffinity capture. The samples were analyzed after storage under four conditions: (i) 4°C on the autosampler for 8 h, (ii) 4°C on the autosampler for 48 h, (iii) after two freeze-thaw cycles, and (iv) -80°C for 5 weeks; control samples were analyzed immediately. Each test case was measured in triplicate. The variability (%CV) and percent difference, comparing peak area ratio between control samples and samples with different handling conditions, were used to evaluate peptide stability. The median %CV for the test samples ranged between 2.5-8.5% in FFPE and 3.1-6.5% in plasma (Supplementary Table 6), within expectations based on repeatability measurements (Supplementary Table 5) for both matrices. The median percent difference relative to the fresh sample ranged 0.1-3.5% in FFPE and 2.0-5.6% in plasma, indicating acceptable stability for the peptides (Figures 1C, D and Supplementary Table 6). In FFPE, three assays (CCL5.EYFYTSGK, ATM.NLS[ph]DIDQSFNK, CD274.NIIQFVHGEEDLK) showed high variability and/or relative difference, indicating these peptides should not be allowed to sit on the autosampler for longer than 8 hours. Likewise, three assays (LGALS1.SFVLNLGK, CD274.NIIQFVHGEEDLK, PTPRC.LFLAEFQSIPR) showed high variability or percent difference in plasma samples after storage on the autosampler, indicating these peptides should be analyzed by LC-MRM within 8 hours.
3.3 Determination of Assay Utility in Tissue and Sample Requirements
We next applied the multiplex assay to a panel of tumor tissue specimens with two aims: (i) evaluate the utility of the assay for measuring endogenous levels of analytes in clinical biospecimens, and (ii) determine sample requirements for analyte detection in tissue where clinical material may be limited. The tissue panel included 110 frozen and 25 FFPE biospecimens collected from 12 different tumor types.
Frozen tissue biospecimens were obtained from 11 tumor types including brain, breast, colorectal, endometrium, head and neck, kidney, lung (squamous cell carcinoma and adenocarcinoma), ovarian, pancreas, and soft tissue sarcoma (Figure 2A). The distribution of cellular heterogeneity was evaluated using the H&E-stained slides, and a supervised algorithm (40) (HALO) was used to assign the fraction of tissue area pertaining to predominant cell types (e.g., adipose, lymphocytes, red blood cells, stroma, and tumor). The distribution of these cellular subtypes is plotted in Figure 2B for the frozen tumors. Overall, most of the biospecimen area can be assigned to tumor and stromal cells, with small proportions of adipose, lymphocytes, and red blood cells. Notably, ovarian tumors contained the highest fraction of adipocytes. These fractions show that the tumors were relatively good quality specimens, with median percentage of tumor cells >50%.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/765898/fimmu-12-765898-HTML-r1/image_m/fimmu-12-765898-g002.jpg)
Figure 2 The immuno-MRM assays show detection in frozen tissues. (A) Frozen tissues were obtained for 110 tumors from 11 tumor types. The number of each type is indicated in the pie chart. (B) The relative fractions of adipose, lymphocytes, red blood cells, stroma, and tumor cells were plotted for tumors with available images (108 out of 110). Cellular microheterogeneity was determined by using the HALO algorithm. Each point represents an individual tumor. Box plots show median (line), inner quartiles (box), and 5-95% range (vertical lines). (C) Distribution of peptide detection plotted as a histogram, showing the number of peptides detected above LOQ across the 110 frozen tumors. (D) Heatmap showing unsupervised clustering of analytes detected above LOQ in > 50% of tumor specimens. Peak area ratios (light:heavy) were normalized for each peptide analyte, and the z-score was used for clustering. (E) Histogram showing correlation of protein expression measured by immuno-MRM with mRNA transcript level determined by RNAseq (52). The median is indicated by a dotted line (0.559).
To determine the ability of the IO-1 immuno-MRM assay to measure endogenous proteins, we used 500 µg aliquots of the frozen tissue lysates as input. Each biospecimen was independently digested, enriched for the IO-1 analytes, and analyzed by LC-MRM. Specificity of endogenous peptide detection was assured by equal retention times and relative areas of light and heavy transitions. Overall, 48/52 peptides corresponding to 45 proteins were detected above LLOQ in the frozen tissue biospecimens. Four peptides were not detected in the biospecimen panel: phospho-ATM pS367, phospho-ATM pS2996, phospho-FAS pY291, and unmodified TNFRSF9. Most peptides were detected in all samples, with a small number of peptides detected below the LLOQ in a subset of samples (see Figure 2C). The median inter-individual CV was 77% (range 32-676%). This variation is likely due to differences between the individual tumors and not differences in the tissue type, since unsupervised clustering of the expression levels for the peptides in the IO-1 multiplex panel (Figure 2D) did not separate the samples by location. This is not surprising, since these proteins were not chosen for site-specific differences in expression, but rather for their role in immunomodulatory functions, which may vary from tumor to tumor. Finally, correlation of expression levels with mRNA (17, 52, 90–97, 107) showed a median correlation coefficient of 0.56 (Figure 2E), consistent with previous reports measuring protein-RNA correlation (15–17).
Given that clinical biospecimens often have limited available material, we next sought to determine the sample requirements for detection of endogenous protein in tissue to use as a guide for expectations in future studies. We used the results of the frozen tissue array (Figure 2) to estimate the minimum amount of tissue needed for the analyte to remain above the LOQ using the distribution of signal-to-noise values measured in tumor samples, compared to the LLOQ characterized in tissue matrix. Figure 3 shows the number of peptides predicted to be detected for tissue inputs ranging from 10-500 µg. Notably, the number of peptides detected as the amount of tissue decreases remains within ~80% of total, even with ten-fold less input. This indicates that the assays are amenable to a range of biospecimen sizes and input amounts.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/765898/fimmu-12-765898-HTML-r1/image_m/fimmu-12-765898-g003.jpg)
Figure 3 Determination of sample requirements for detection in tissue. The number of analytes predicted for decreasing amounts of tissue was determined from the signal-to-noise ratio measured using 500 μg protein digest input. Error bars show the 95% confidence interval.
To evaluate the potential for applying the assays to retrospectively collected biospecimens, we evaluated the detection of peptides in FFPE tissues representing 5 tumor types, including colorectal, head and neck, lung (squamous cell carcinoma and adenocarcinoma), and prostate (Supplementary Figure 2). Optimized sample processing conditions were used to process the samples (48), however the limited material available for these samples resulted in sample inputs ranging from 2 to over 10-fold less material compared to frozen tissue biopsies. Despite these limitations, the assay successfully quantified 33/52 peptides above LOQ in more than half of samples (Supplementary Figure 2).
3.4 The Multiplexed Assay Shows Utility for Measurements of Plasma Protein Expression
To evaluate the ability of the assay to measure protein levels in plasma, we applied the assays to a panel of 45 plasma samples from patients with 3 different tumor types (breast, colon, ovarian) (Figure 4A). We used 100 µL aliquots of plasma as input for the assay. Each specimen was independently digested, enriched for the IO-1 analytes, and analyzed by LC-MRM. Specificity of endogenous peptide detection was assured by equal retention times and relative areas of light and heavy transitions. Overall, 38/52 peptides corresponding to 37 proteins were detected above LOQ in the plasma samples. The histogram of number of peptides detected across the panel showed that over half (n=24) of peptides were detected in all samples analyzed, with a small number of peptides detected below the LOQ in a subset of samples (Figure 4B) and 14 peptides not detected in any sample, likely due to the analytes not being secreted or shed to quantifiable levels. Compared to the variation seen in the tissue panel, there were less inter-individual differences in plasma, with the median inter-individual CV being 47% (range 10-152%). Like the tissue samples, unsupervised clustering of the expression levels for the peptides in the IO-1 multiplex panel (Figure 4C) did not separate the samples by tumor type.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/765898/fimmu-12-765898-HTML-r1/image_m/fimmu-12-765898-g004.jpg)
Figure 4 Utility of the assays for measurement of protein expression in plasma. (A) Plasma samples were obtained for 45 patients with breast, colorectal, or ovarian tumors, as indicated in the pie chart. (B) Distribution of peptide detection plotted as a histogram, showing the number of peptides detected above LOQ across the 45 plasmas using 100 μL aliquots of plasma as input. (C) Heatmap showing unsupervised clustering of analytes detected above LOQ in >50% of plasma samples. Peak area ratios (light:heavy) were normalized for each peptide analyte, and the z-score was used for clustering. (D) Histogram showing the correlation of protein expression levels in patients for whom both tissue and plasma samples were available.
Finally, there were 25 plasma samples for which we also analyzed tissue specimens from the same patient. The distribution of correlation coefficients for expression levels found in tissue and plasma from the same patient are plotted in Figure 4D. Overall, 31 peptides were detected in both tissue and plasma specimens, but there was poor intra-patient correlation between the two sample types. Based on previous results (108), we expected there may be high overlap in the peptides detected; however, it was unknown how the relative quantities may be correlated. Because protein abundance may also vary depending on specimen collection conditions (e.g., blood contamination in tissue specimens and/or cell shearing during plasma collection), we also examined if blood contamination in tissue specimens affected the correlation of peptides detected in tissue and plasma. A “blood contamination” score was assigned to each tissue specimen using the intensity of gel electrophoresis bands corresponding to the molecular weights of albumin and hemoglobin in each tissue lysate; however, correlation between tissue and plasma levels did not improve using only the samples with high “blood contamination” scores (data not shown).
3.5 The Antibodies Are a Useful Resource for Additional Applications
In parallel to their development for immuno-MRM, the new antibodies generated for this study were characterized for reactivity to the target proteins in traditional immunoassay approaches. Antibodies were tested for utility in traditional Western blotting or Wes system using purified recombinant proteins (where available) and/or lysates from cancer cell lines. Blots were deemed positive if a prominent band at the expected molecular weight (given by the vendor for the recombinant proteins or based on the molecular weight reported in UniProt.org) was observed. For this set of antibodies, 47% (14/30) of the antibodies tested (14/30) were positive against recombinant proteins, and 60% (12/20) were positive against lysates (not all antibodies that were negative for the recombinant protein were tested against cell line lysates). These results were consistent with expectations based on previous work (109). Antibodies scoring positive in Western blotting were tested for application in protein array using the NCI-60 cell line array. Overall, 3 of the antibodies tested were positive, based on normalized levels above background. The data are summarized in Supplementary Table 7. The monoclonal antibodies and characterization data are available as a resource for the research community via the CPTAC antibody portal (antibodies.cancer.gov; see Supplementary Table 2).
4 Discussion
We present the development and characterization of a multiplex panel of 52 assays for quantifying 46 immuno-modulatory proteins in tissue and plasma biospecimens. The panel shows excellent quantitative characteristics in tissue and plasma, with wide dynamic range and high precision. We used a panel of 180 biospecimens (135 tissue samples, 45 plasma samples) to demonstrate feasibility of endogenous measurements using the IO-1 immuno-MRM assay panel. In addition to high detection rates in frozen biopsies, the assays were successful in FFPE tissues, opening the possibility for retrospective studies.
The immuno-MRM approach offers sufficient sensitivity for analysis of low abundance proteins due to the ability of antibodies to enrich the target peptides from complex samples. Enrichment factors are analyte-dependent and range between 102-104 (110). The enriched sample contains a much-simplified background matrix of peptides, which improves ion suppression in the mass spectrometer. While the recovery efficiencies of the peptides in the enrichment step vary, the enrichment factors greatly overcome any sample losses associated with the additional sample handling (i.e., compared to direct analysis without enrichment), making the immuno-MRM approach highly advantageous for analysis of low abundance proteins and peptides.
Quantification of protein expression offers a powerful complementarity to existing technologies. Because mRNA does not always correlate, and thus fails to predict protein levels, the direct measurement of proteins gives a more accurate link to phenotypes. Likewise, the quantitative nature of the immuno-MRM panel offers a powerful complementarity to measurements performed by IHC. IHC is semiquantitative in nature, and quantification is susceptible to technical factors as well as subjective interpretation. However, IHC does provide useful spatial context for protein expression. The immuno-MRM measurements add complementary quantitative information that furthers the potential for useful interpretation of results by pathologists.
A large advantage of LC-MRM-MS is that analytes can be readily multiplexed. We used a panel of 50 monoclonal antibodies for enriching peptides from the digested samples. The high specificity of the mass spectrometer allows for analysis of much larger multiplex groups (compared to traditional immuno-assays), since interferences can be readily detected and usually avoided, and since MS-based assays have larger dynamic ranges compared with conventional immunoassays. Indeed, immuno-MRM assays have been demonstrated for panels up to 150 antibodies (28).
To produce actionable information from tumor biopsies, the assay must be sensitive enough to quantify protein networks in core needle biopsies or 5-10 micron sections cut from FFPE blocks. Although we have demonstrated that immuno-MRM produces sufficient sensitivity and specificity in these biospecimen types, the quality of specimens collected in clinical trials and across different sites are likely to produce many specimens with much more limited material. To address this problem, future work will focus on optimizing the sensitivity of the approach by improving analyte recovery through addition-only sample processing protocols and minimizing losses due to surface adsorption in reaction vessels and lab plastics. These advances will enable a greater diversity of specimens for analysis. Another challenge in making bulk tissue measurements is the need to understand the contribution of tumor heterogeneity to the measurements. Other bulk tissue assays have addressed tissue heterogeneity by requiring a minimum percent tumor cellularity (111). The tumors in this study were associated with a minimum cellularity requirement of ≥60% (by number). We have previously demonstrated that when cellularity is low, normalization proteins can be valuable in controlling for heterogeneity. For example, in measuring the HER2 protein in samples with low tumor cellularity, we were able to improve the concordance between an immuno-MRM assay and the predicate clinical IHC assay by normalizing the MRM measurements to GAPDH (31). Therefore, based on the HER2 study, we have included GAPDH in the final validated assay to accommodate future applications where tumor cellularity may be low.
These assays are made available to the research community as part of a larger effort under the National Cancer Institute’s Beau Biden National Cancer Moonshot [APOLLO network (44)] to quantify immuno-modulatory proteins using MRM-based assays. The monoclonal antibodies, characterization data, and SOPs are freely accessible to the research community through NCI’s CPTAC Assay Portal (assays.cancer.gov) (46, 47) and CPTAC Antibody Portal (antibodies.cancer.gov). Applied to clinical biospecimens, the assays have the potential for expanding correlative studies, establishing metrics of on-target inflammation and tumor response, or studying cancer immunity mechanisms to identify new therapeutic targets.
Data Availability Statement
The names of the repository/repositories and accession number(s) can be found below: ProteomeXchange: accession number PXD029167 (doi: 10.6069/4wt1-t595) or follow the link: https://panoramaweb.org/IO1immunoMRM.url.
Ethics Statement
The studies involving human participants were reviewed and approved by Fred Hutchinson Cancer Research Center Institutional Review Office. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
AP conceived or designed the work. JW, RL, LZ, RS, AB, DH, UV, TW, JJK, RI, CL, OM, TL, SC, MT, TC, SH, SW, RR, JGK, JR, JAK, and CR acquired, analyzed, or interpreted the data. KM, ND, MB, HY, SA, EB, SG-B, WB, TH, HR, JW, and AP supervised studies. JW and AP drafted or substantively revised the work. All authors read and approved the final manuscript.
Funding
This research has been funded in whole or in part with Federal funds from the National Cancer Institute (NCI), National Institutes of Health, under the NCI Beau Biden National Cancer Moonshot (Task Order No. HHSN26100025, Applied Proteogenomics Organizational Learning and Outcomes, under Contract No. HHSN261201500003I), the NCI Clinical Proteomics Tumor Analysis Consortium (CPTAC, grant no. U01CA214114), the NCI Academic Industrial Partnership (grant no. R01CA235575), the NCI Research Specialist program (grant no. R50CA211499), in part by the NCI CPTAC Antibody Characterization Program, and a generous donation from the Aven Foundation. This work was partially supported by the VA Office of Research and Development, Cooperative Studies Program (CSP). The views expressed are those of the authors and do not necessarily reflect the position or policy of the Department of Veterans Affairs or the United States government.
Author Disclaimer
The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Karen Ketchum and the team at ESAC, Inc. for technical support in assembling data.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.765898/full#supplementary-material
References
1. Sharma P, Siddiqui BA, Anandhan S, Yadav SS, Subudhi SK, Gao J, et al. The Next Decade of Immune Checkpoint Therapy. Cancer Discov (2021) 11(4):838–57. doi: 10.1158/2159-8290.CD-20-1680
2. Saxena M, van der Burg SH, Melief CJM, Bhardwaj N. Therapeutic Cancer Vaccines. Nat Rev Cancer (2021) 21(6):360–78. doi: 10.1038/s41568-021-00346-0
3. Larson RC, Maus MV. Recent Advances and Discoveries in the Mechanisms and Functions of CAR T Cells. Nat Rev Cancer (2021) 21(3):145–61. doi: 10.1038/s41568-020-00323-z
4. Abdel-Wahab N, Shah M, Suarez-Almazor ME. Adverse Events Associated With Immune Checkpoint Blockade in Patients With Cancer: A Systematic Review of Case Reports. PloS One (2016) 11(7):e0160221. doi: 10.1371/journal.pone.0160221
5. Perdigoto AL, Kluger H, Herold KC. Adverse Events Induced by Immune Checkpoint Inhibitors. Curr Opin Immunol (2021) 69:29–38. doi: 10.1016/j.coi.2021.02.002
6. Pilard C, Ancion M, Delvenne P, Jerusalem G, Hubert P, Herfs M. Cancer Immunotherapy: It’s Time to Better Predict Patients’ Response. Br J Cancer (2021) 125(7):927–38. doi: 10.1038/s41416-021-01413-x
7. Sullivan RJ, Weber JS. Immune-Related Toxicities of Checkpoint Inhibitors: Mechanisms and Mitigation Strategies. Nat Rev Drug Discov (2021). doi: 10.1038/s41573-021-00259-5
8. Johnson DB, Balko JM, Compton ML, Chalkias S, Gorham J, Xu Y, et al. Fulminant Myocarditis With Combination Immune Checkpoint Blockade. N Engl J Med (2016) 375(18):1749–55. doi: 10.1056/NEJMoa1609214
9. Chen DS, Mellman I. Oncology Meets Immunology: The Cancer-Immunity Cycle. Immunity (2013) 39(1):1–10. doi: 10.1016/j.immuni.2013.07.012
10. Hirsch FR, McElhinny A, Stanforth D, Ranger-Moore J, Jansson M, Kulangara K, et al. PD-L1 Immunohistochemistry Assays for Lung Cancer: Results From Phase 1 of the Blueprint PD-L1 IHC Assay Comparison Project. J Thorac Oncol (2017) 12(2):208–22. doi: 10.1016/j.jtho.2016.11.2228
11. Keren L, Bosse M, Thompson S, Risom T, Vijayaragavan K, McCaffrey E, et al. MIBI-TOF: A Multiplexed Imaging Platform Relates Cellular Phenotypes and Tissue Structure. Sci Adv (2019) 5(10):eaax5851. doi: 10.1126/sciadv.aax5851
12. Gadalla R, Noamani B, MacLeod BL, Dickson RJ, Guo M, Xu W, et al. Validation of CyTOF Against Flow Cytometry for Immunological Studies and Monitoring of Human Cancer Clinical Trials. Front Oncol (2019) 9:415. doi: 10.3389/fonc.2019.00415
13. Baker M. Reproducibility Crisis: Blame it on the Antibodies. Nature (2015) 521(7552):274–6. doi: 10.1038/521274a
14. Cesano A. Ncounter((R)) PanCancer Immune Profiling Panel (NanoString Technologies, Inc., Seattle, Wa). J Immunother Cancer (2015) 3:42. doi: 10.1186/s40425-015-0088-7
15. Zhang B, Wang J, Wang X, Zhu J, Liu Q, Shi Z, et al. Proteogenomic Characterization of Human Colon and Rectal Cancer. Nature (2014) 513(7518):382–7. doi: 10.1038/nature13438
16. Mertins P, Mani DR, Ruggles KV, Gillette MA, Clauser KR, Wang P, et al. Proteogenomics Connects Somatic Mutations to Signalling in Breast Cancer. Nature (2016) 534(7605):55–62. doi: 10.1038/nature18003
17. Zhang H, Liu T, Zhang Z, Payne SH, Zhang B, McDermott JE, et al. Integrated Proteogenomic Characterization of Human High-Grade Serous Ovarian Cancer. Cell (2016) 166(3):755–65. doi: 10.1016/j.cell.2016.05.069
18. Rodriguez H, Zenklusen JC, Staudt LM, Doroshow JH, Lowy DR. The Next Horizon in Precision Oncology: Proteogenomics to Inform Cancer Diagnosis and Treatment. Cell (2021) 184(7):1661–70. doi: 10.1016/j.cell.2021.02.055
19. Gajewski TF, Meng Y, Blank C, Brown I, Kacha A, Kline J, et al. Immune Resistance Orchestrated by the Tumor Microenvironment. Immunol Rev (2006) 213:131–45. doi: 10.1111/j.1600-065X.2006.00442.x
20. CLSI. Quantitative Measurement of Proteins and Peptides by Mass Spectrometry, in CLSI Guideline C64. Wayne, PA: C.a.L.S. Institute (2021).
21. Picotti P, Bodenmiller B, Aebersold R. Proteomics Meets the Scientific Method. Nat Methods (2013) 10(1):24–7. doi: 10.1038/nmeth.2291
22. Gillette MA, Carr SA. Quantitative Analysis of Peptides and Proteins in Biomedicine by Targeted Mass Spectrometry. Nat Methods (2013) 10(1):28–34. doi: 10.1038/nmeth.2309
23. Whiteaker JR, Paulovich AG. Peptide Immunoaffinity Enrichment Coupled With Mass Spectrometry for Peptide and Protein Quantification. Clin Lab Med (2011) 31(3):385–96. doi: 10.1016/j.cll.2011.07.004
24. Whiteaker JR, Zhao L, Anderson L, Paulovich AG. An Automated and Multiplexed Method for High Throughput Peptide Immunoaffinity Enrichment and Multiple Reaction Monitoring Mass Spectrometry-Based Quantification of Protein Biomarkers. Mol Cell Proteomics (2010) 9(1):184–96. doi: 10.1074/mcp.M900254-MCP200
25. Whiteaker JR, Zhao L, Ivey RG, Sanchez-Bonilla M, Moore HD, Schoenherr RM, et al. Targeted Mass Spectrometry Enables Robust Quantification of FANCD2 Mono-Ubiquitination in Response to DNA Damage. DNA Repair (Amst) (2018) 65:47–53. doi: 10.1016/j.dnarep.2018.03.003
26. Whiteaker JR, Zhao L, Saul R, Kaczmarczyk JA, Schoenherr RM, Moore HD, et al. A Multiplexed Mass Spectrometry-Based Assay for Robust Quantification of Phosphosignaling in Response to DNA Damage. Radiat Res (2018) 189(5):505–18. doi: 10.1667/RR14963.1
27. Whiteaker JR, Zhao L, Lin C, Yan P, Wang P, Paulovich AG. Sequential Multiplexed Analyte Quantification Using Peptide Immunoaffinity Enrichment Coupled to Mass Spectrometry. Mol Cell Proteomics (2012) 11(6):M111 015347. doi: 10.1074/mcp.M111.015347
28. Ippoliti PJ, Kuhn E, Mani DR, Fagbami L, Keshishian H, Burgess MW, et al. Automated Microchromatography Enables Multiplexing of Immunoaffinity Enrichment of Peptides to Greater Than 150 for Targeted MS-Based Assays. Anal Chem (2016) 88(15):7548–55. doi: 10.1021/acs.analchem.6b00946
29. Kuhn E, Whiteaker JR, Mani DR, Jackson AM, Zhao L, Pope ME, et al. Interlaboratory Evaluation of Automated, Multiplexed Peptide Immunoaffinity Enrichment Coupled to Multiple Reaction Monitoring Mass Spectrometry for Quantifying Proteins in Plasma. Mol Cell Proteomics (2012) 11(6):M111 013854. doi: 10.1074/mcp.M111.013854
30. Netzel BC, Grant RP, Hoofnagle AN, Rockwood AL, Shuford CM, Grebe SK. First Steps Toward Harmonization of LC-MS/MS Thyroglobulin Assays. Clin Chem (2016) 62(1):297–9. doi: 10.1373/clinchem.2015.245266
31. Kennedy JJ, Abbatiello SE, Kim K, Yan P, Whiteaker JR, Lin C, et al. Demonstrating the Feasibility of Large-Scale Development of Standardized Assays to Quantify Human Proteins. Nat Methods (2014) 11(2):149–55. doi: 10.1038/nmeth.2763
32. Neubert H, Shuford CM, Olah TV, Garofolo F, Schultz GA, Jones BR, et al. Protein Biomarker Quantification by Immunoaffinity Liquid Chromatography-Tandem Mass Spectrometry: Current State and Future Vision. Clin Chem (2020) 66(2):282–301. doi: 10.1093/clinchem/hvz022
33. Kushnir MM, Rockwood AL, Roberts WL, Abraham D, Hoofnagle AN, Meikle AW. Measurement of Thyroglobulin by Liquid Chromatography-Tandem Mass Spectrometry in Serum and Plasma in the Presence of Antithyroglobulin Autoantibodies. Clin Chem (2013) 59(6):982–90. doi: 10.1373/clinchem.2012.195594
34. van den Broek I, Romijn FP, Nouta J, van der Laarse A, Drijfhout JW, Smit NP, et al. Automated Multiplex LC-MS/MS Assay for Quantifying Serum Apolipoproteins A-I, B, C-I, C-II, C-III, and E With Qualitative Apolipoprotein E Phenotyping. Clin Chem (2016) 62(1):188–97. doi: 10.1373/clinchem.2015.246702
35. Kushnir MM, Rockwood AL, Strathmann FG, Frank EL, Straseski JA, Meikle AW. LC-MS/MS Measurement of Parathyroid Hormone-Related Peptide. Clin Chem (2016) 62(1):218–26. doi: 10.1373/clinchem.2015.244012
36. Jones GN, Rooney C, Griffin N, Roudier M, Young LA, Garcia-Trinidad A, et al. Prad50: A Novel and Clinically Applicable Pharmacodynamic Biomarker of Both ATM and ATR Inhibition Identified Using Mass Spectrometry and Immunohistochemistry. Br J Cancer (2018) 119(10):1233–43. doi: 10.1038/s41416-018-0286-4
37. Collins CJ, Chang IJ, Jung S, Dayuha R, Whiteaker JR, Segundo GRS, et al. Rapid Multiplexed Proteomic Screening for Primary Immunodeficiency Disorders From Dried Blood Spots. Front Immunol (2018) 9:2756. doi: 10.3389/fimmu.2018.02756
38. Collins CJ, Yi F, Dayuha R, Whiteaker JR, Ochs HD, Freeman A, et al. Multiplexed Proteomic Analysis for Diagnosis and Screening of Five Primary Immunodeficiency Disorders From Dried Blood Spots. Front Immunol (2020) 11:464. doi: 10.3389/fimmu.2020.00464
39. Jung S, Whiteaker JR, Zhao L, Yoo HW, Paulovich AG, Hahn SH. Quantification of ATP7B Protein in Dried Blood Spots by Peptide Immuno-SRM as a Potential Screen for Wilson’s Disease. J Proteome Res (2017) 16(2):862–71. doi: 10.1021/acs.jproteome.6b00828
40. Kennedy JJ, Whiteaker JR, Kennedy LC, Bosch DE, Lerch ML, Schoenherr RM, et al. Quantification of Human Epidermal Growth Factor Receptor 2 by Immunopeptide Enrichment and Targeted Mass Spectrometry in Formalin-Fixed Paraffin-Embedded and Frozen Breast Cancer Tissues. Clin Chem (2021) 67(7):1008–18. doi: 10.1093/clinchem/hvab047
41. Schoenherr RM, Voytovich UJ, Whiteaker JR, Gadi VK, Gralow JR, Paulovich AG. Quantification of Estrogen Receptor, Progesterone Receptor, and Human Epidermal Growth Factor Receptor 2 Protein Expression in Bone Biopsies by Targeted Mass Spectrometry Without Acid Decalcification. Clin Chem (2020) 66(11):1459–61. doi: 10.1093/clinchem/hvaa220
42. Hembrough T, Thyparambil S, Liao WL, Darfler MM, Abdo J, Bengali KM, et al. Application of Selected Reaction Monitoring for Multiplex Quantification of Clinically Validated Biomarkers in Formalin-Fixed, Paraffin-Embedded Tumor Tissue. J Mol Diagn (2013) 15(4):454–65. doi: 10.1016/j.jmoldx.2013.03.002
43. Do M, Kim H, Yeo I, Lee J, Park IA, Ryu HS, et al. Clinical Application of Multiple Reaction Monitoring-Mass Spectrometry to Human Epidermal Growth Factor Receptor 2 Measurements as a Potential Diagnostic Tool for Breast Cancer Therapy. Clin Chem (2020) 66(10):1339–48. doi: 10.1093/clinchem/hvaa178
44. Fiore LD, Rodriguez H, Shriver CD. Collaboration to Accelerate Proteogenomics Cancer Care: The Department of Veterans Affairs, Department of Defense, and the National Cancer Institute’s Applied Proteogenomics OrganizationaL Learning and Outcomes (APOLLO) Network. Clin Pharmacol Ther (2017) 101(5):619–21. doi: 10.1002/cpt.658
45. Sharpless NE, Singer DS. Progress and Potential: The Cancer Moonshot. Cancer Cell (2021) 39(7):889–94. doi: 10.1016/j.ccell.2021.04.015
46. Whiteaker JR, Halusa GN, Hoofnagle AN, Sharma V, MacLean B, Yan P, et al. Using the CPTAC Assay Portal to Identify and Implement Highly Characterized Targeted Proteomics Assays. Methods Mol Biol (2016) 1410:223–36. doi: 10.1007/978-1-4939-3524-6_13
47. Whiteaker JR, Halusa GN, Hoofnagle AN, Sharma V, MacLean B, Yan P, et al. CPTAC Assay Portal: A Repository of Targeted Proteomic Assays. Nat Methods (2014) 11(7):703–4. doi: 10.1038/nmeth.3002
48. Kennedy JJ, Whiteaker JR, Schoenherr RM, Yan P, Allison K, Shipley M, et al. Optimized Protocol for Quantitative Multiple Reaction Monitoring-Based Proteomic Analysis of Formalin-Fixed, Paraffin-Embedded Tissues. J Proteome Res (2016) 15(8):2717–28. doi: 10.1021/acs.jproteome.6b00245
49. Zhao L, Whiteaker JR, Pope ME, Kuhn E, Jackson A, Anderson NL, et al. Quantification of Proteins Using Peptide Immunoaffinity Enrichment Coupled With Mass Spectrometry. J Vis Exp (2011) 53):2812. doi: 10.3791/2812
50. MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, et al. Skyline: An Open Source Document Editor for Creating and Analyzing Targeted Proteomics Experiments. Bioinformatics (2010) 26(7):966–8. doi: 10.1093/bioinformatics/btq054
51. Pino LK, Searle BC, Bollinger JG, Nunn B, MacLean B, MacCoss MJ. The Skyline Ecosystem: Informatics for Quantitative Mass Spectrometry Proteomics. Mass Spectrom Rev (2020) 39(3):229–44. doi: 10.1002/mas.21540
52. Edwards NJ, Oberti M, Thangudu RR, Cai S, McGarvey PB, Jacob S, et al. The CPTAC Data Portal: A Resource for Cancer Proteomics Research. J Proteome Res (2015) 14(6):2707–13. doi: 10.1021/pr501254j
53. Williams R, Chung JY, Ylaya K, Whiteley G, Hewitt SM. Characterizations and Validations of Novel Antibodies Toward Translational Research. Proteomics Clin Appl (2010) 4(6-7):618–25. doi: 10.1002/prca.200900186
54. Chung JY, Lee SJ, Kris Y, Braunschweig T, Traicoff JL, Hewitt SM. A Well-Based Reverse-Phase Protein Array Applicable to Extracts From Formalin-Fixed Paraffin-Embedded Tissue. Proteomics Clin Appl (2008) 2(10-11):1539–47. doi: 10.1002/prca.200800005
55. Sharma V, Eckels J, Schilling B, Ludwig C, Jaffe JD, MacCoss MJ, et al. Panorama Public: A Public Repository for Quantitative Data Sets Processed in Skyline. Mol Cell Proteomics (2018) 17(6):1239–44. doi: 10.1074/mcp.RA117.000543
56. Hegde PS, Chen DS. Top 10 Challenges in Cancer Immunotherapy. Immunity (2020) 52(1):17–35. doi: 10.1016/j.immuni.2019.12.011
57. Moss ML, Jin SL, Milla ME, Bickett DM, Burkhart W, Carter HL, et al. Cloning of a Disintegrin Metalloproteinase That Processes Precursor Tumour-Necrosis Factor-Alpha. Nature (1997) 385(6618):733–6. doi: 10.1038/385733a0
58. Arcone R, Arpaia G, Ruoppolo M, Malorni A, Pucci P, Marino G, et al. Structural Characterization of a Biologically Active Human Lipocortin 1 Expressed in Escherichia Coli. Eur J Biochem (1993) 211(1-2):347–55. doi: 10.1111/j.1432-1033.1993.tb19904.x
59. Munder M, Schneider H, Luckner C, Giese T, Langhans CD, Fuentes JM, et al. Suppression of T-Cell Functions by Human Granulocyte Arginase. Blood (2006) 108(5):1627–34. doi: 10.1182/blood-2006-11-010389
60. Lumsden JM, McCarty T, Petiniot LK, Shen R, Barlow C, Wynn TA, et al. Immunoglobulin Class Switch Recombination is Impaired in Atm-Deficient Mice. J Exp Med (2004) 200(9):1111–21. doi: 10.1084/jem.20041074
61. Cheah MT, Chen JY, Sahoo D, Contreras-Trujillo H, Volkmer AK, Scheeren FA, et al. CD14-Expressing Cancer Cells Establish the Inflammatory and Proliferative Tumor Microenvironment in Bladder Cancer. Proc Natl Acad Sci USA (2015) 112(15):4725–30. doi: 10.1073/pnas.1424795112
62. Imam R, Chang Q, Black M, Yu C, Cao W. CD47 Expression and CD163(+) Macrophages Correlated With Prognosis of Pancreatic Neuroendocrine Tumor. BMC Cancer (2021) 21(1):320. doi: 10.1186/s12885-021-08045-7
63. Lajaunias F, Dayer JM, Chizzolini C. Constitutive Repressor Activity of CD33 on Human Monocytes Requires Sialic Acid Recognition and Phosphoinositide 3-Kinase-Mediated Intracellular Signaling. Eur J Immunol (2005) 35(1):243–51. doi: 10.1002/eji.200425273
64. Stein R, Mattes MJ, Cardillo TM, Hansen HJ, Chang CH, Burton J, et al. CD74: A New Candidate Target for the Immunotherapy of B-Cell Neoplasms. Clin Cancer Res (2007) 13(18 Pt 2):5556s–63s. doi: 10.1158/1078-0432.CCR-07-1167
65. Dankner M, Gray-Owen SD, Huang YH, Blumberg RS, Beauchemin N. CEACAM1 as a Multi-Purpose Target for Cancer Immunotherapy. Oncoimmunology (2017) 6(7):e1328336. doi: 10.1080/2162402X.2017.1328336
66. Yang M, Lu J, Zhang G, Wang Y, He M, Xu Q, et al. CXCL13 Shapes Immunoactive Tumor Microenvironment and Enhances the Efficacy of PD-1 Checkpoint Blockade in High-Grade Serous Ovarian Cancer. J Immunother Cancer (2021) 9(1):e001136. doi: 10.1136/jitc-2020-001136
67. Bastid J, Cottalorda-Regairaz A, Alberici G, Bonnefoy N, Eliaou JF, Bensussan A. ENTPD1/CD39 is a Promising Therapeutic Target in Oncology. Oncogene (2013) 32(14):1743–51. doi: 10.1038/onc.2012.269
68. Dostert C, Grusdat M, Letellier E, Brenner D. The TNF Family of Ligands and Receptors: Communication Modules in the Immune System and Beyond. Physiol Rev (2019) 99(1):115–60. doi: 10.1152/physrev.00045.2017
69. Nicholls C, Li H, Liu JP. GAPDH: A Common Enzyme With Uncommon Functions. Clin Exp Pharmacol Physiol (2012) 39(8):674–9. doi: 10.1111/j.1440-1681.2011.05599.x
70. Wang B, Jiang H, Zhou T, Ma N, Liu W, Wang Y, et al. Expression of ICOSL is Associated With Decreased Survival in Invasive Breast Cancer. PeerJ (2019) 7:e6903. doi: 10.7717/peerj.6903
71. Ma Z, Li W, Yoshiya S, Xu Y, Hata M, El-Darawish Y, et al. Augmentation of Immune Checkpoint Cancer Immunotherapy With IL18. Clin Cancer Res (2016) 22(12):2969–80. doi: 10.1158/1078-0432.CCR-15-1655
72. Kim Y, Shin Y, Kang GH. Prognostic Significance of CD103+ Immune Cells in Solid Tumor: A Systemic Review and Meta-Analysis. Sci Rep (2019) 9(1):3808. doi: 10.1038/s41598-019-40527-4
73. Li Y, Yang S, Yue H, Yuan D, Li L, Zhao J, et al. Unraveling LGALS1 as a Potential Immune Checkpoint and a Predictor of the Response to Anti-PD1 Therapy in Clear Cell Renal Carcinoma. Pathol Oncol Res (2020) 26(3):1451–8. doi: 10.1007/s12253-019-00710-4
74. Bommhardt U, Schraven B, Simeoni L. Beyond TCR Signaling: Emerging Functions of Lck in Cancer and Immunotherapy. Int J Mol Sci (2019) 20(14):3500. doi: 10.3390/ijms20143500
75. Betzler AC, Theodoraki MN, Schuler PJ, Doscher J, Laban S, Hoffmann TK, et al. NF-kappaB and Its Role in Checkpoint Control. Int J Mol Sci (2020) 21(11):3949. doi: 10.3390/ijms21113949
76. Jiang T, Xu X, Qiao M, Li X, Zhao C, Zhou F, et al. Comprehensive Evaluation of NT5E/CD73 Expression and its Prognostic Significance in Distinct Types of Cancers. BMC Cancer (2018) 18(1):267. doi: 10.1186/s12885-018-4073-7
77. Huang KC, Chiang SF, Chen TW, Chen WT, Yang PC, Ke TW, et al. Prognostic Relevance of Programmed Cell Death 1 Ligand 2 (PDCD1LG2/PD-L2) in Patients With Advanced Stage Colon Carcinoma Treated With Chemotherapy. Sci Rep (2020) 10(1):22330. doi: 10.1038/s41598-020-79419-3
78. Bergom C, Gao C, Newman PJ. Mechanisms of PECAM-1-Mediated Cytoprotection and Implications for Cancer Cell Survival. Leuk Lymphoma (2005) 46(10):1409–21. doi: 10.1080/10428190500126091
79. Li Y, Huang J, Sun J, Xiang S, Yang D, Ying X, et al. The Transcription Levels and Prognostic Values of Seven Proteasome Alpha Subunits in Human Cancers. Oncotarget (2017) 8(3):4501–19. doi: 10.18632/oncotarget.13885
80. Markosyan N, Li J, Sun YH, Richman LP, Lin JH, Yan F, et al. Tumor Cell-Intrinsic EPHA2 Suppresses Anti-Tumor Immunity by Regulating PTGS2 (COX-2). J Clin Invest (2019) 129(9):3594–609. doi: 10.1172/JCI127755
81. Raiter A, Zlotnik O, Lipovetsky J, Mugami S, Dar S, Lubin I, et al. A Novel Role for an Old Target: CD45 for Breast Cancer Immunotherapy. Oncoimmunology (2021) 10(1):1929725. doi: 10.1080/2162402X.2021.1929725
82. Samstein RM, Riaz N. The DNA Damage Response in Immunotherapy and Radiation. Adv Radiat Oncol (2018) 3(4):527–33. doi: 10.1016/j.adro.2018.08.017
83. Haura EB, Turkson J, Jove R. Mechanisms of Disease: Insights Into the Emerging Role of Signal Transducers and Activators of Transcription in Cancer. Nat Clin Pract Oncol (2005) 2(6):315–24. doi: 10.1038/ncponc0195
84. Henle AM, Nassar A, Puglisi-Knutson D, Youssef B, Knutson KL. Downregulation of TAP1 and TAP2 in Early Stage Breast Cancer. PloS One (2017) 12(11):e0187323. doi: 10.1371/journal.pone.0187323
85. Zhang D, Bi J, Liang Q, Wang S, Zhang L, Han F, et al. VCAM1 Promotes Tumor Cell Invasion and Metastasis by Inducing EMT and Transendothelial Migration in Colorectal Cancer. Front Oncol (2020) 10:1066. doi: 10.3389/fonc.2020.01066
86. Lines JL, Pantazi E, Mak J, Sempere LF, Wang L, O’Connell S, et al. VISTA is an Immune Checkpoint Molecule for Human T Cells. Cancer Res (2014) 74(7):1924–32. doi: 10.1158/0008-5472.CAN-13-1504
87. Altan M, Kidwell KM, Pelekanou V, Carvajal-Hausdorf DE, Schalper KA, Toki MI, et al. Association of B7-H4, PD-L1, and Tumor Infiltrating Lymphocytes With Outcomes in Breast Cancer. NPJ Breast Cancer (2018) 4:40. doi: 10.1038/s41523-018-0095-1
88. Hoofnagle AN, Whiteaker JR, Carr SA, Kuhn E, Liu T, Massoni SA, et al. Recommendations for the Generation, Quantification, Storage, and Handling of Peptides Used for Mass Spectrometry-Based Assays. Clin Chem (2016) 62(1):48–69. doi: 10.1373/clinchem.2015.250563
89. Lange V, Picotti P, Domon B, Aebersold R. Selected Reaction Monitoring for Quantitative Proteomics: A Tutorial. Mol Syst Biol (2008) 4:222. doi: 10.1038/msb.2008.61
90. Wang LB, Karpova A, Gritsenko MA, Kyle JE, Cao S, Li Y, et al. Proteogenomic and Metabolomic Characterization of Human Glioblastoma. Cancer Cell (2021) 39(4):509–528 e20. doi: 10.1016/j.ccell.2021.01.006
91. Krug K, Jaehnig EJ, Satpathy S, Blumenberg L, Karpova A, Anurag M, et al. Proteogenomic Landscape of Breast Cancer Tumorigenesis and Targeted Therapy. Cell (2020) 183(5):1436–1456 e31. doi: 10.1016/j.cell.2020.10.036
92. Vasaikar S, Huang C, Wang X, Petyuk VA, Savage SR, Wen B, et al. Proteogenomic Analysis of Human Colon Cancer Reveals New Therapeutic Opportunities. Cell (2019) 177(4):1035–1049 e19. doi: 10.1016/j.cell.2019.03.030
93. Huang C, Chen L, Savage SR, Eguez RV, Dou Y, Li Y, et al. Proteogenomic Insights Into the Biology and Treatment of HPV-Negative Head and Neck Squamous Cell Carcinoma. Cancer Cell (2021) 39(3):361–379 e16. doi: 10.1016/j.ccell.2020.12.007
94. Clark DJ, Dhanasekaran SM, Petralia F, Pan J, Song X, Hu Y, et al. Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma. Cell (2019) 179(4):964–983 e31. doi: 10.1016/j.cell.2019.10.007
95. Clark DJ, Dhanasekaran SM, Petralia F, Pan J, Song X, Hu Y, et al. Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma. Cell (2020) 180(1):207. doi: 10.1016/j.cell.2019.12.026
96. Gillette MA, Satpathy S, Cao S, Dhanasekaran SM, Vasaikar SV, Krug K, et al. Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma. Cell (2020) 182(1):200–225.e35. doi: 10.1016/j.cell.2020.06.013
97. Hu Y, Pan J, Shah P, Ao M, Thomas SN, Liu Y, et al. Integrated Proteomic and Glycoproteomic Characterization of Human High-Grade Serous Ovarian Carcinoma. Cell Rep (2020) 33(3):108276. doi: 10.1016/j.celrep.2020.108276
98. Kennedy JJ, Yan P, Zhao L, Ivey RG, Voytovich UJ, Moore HD, et al. Immobilized Metal Affinity Chromatography Coupled to Multiple Reaction Monitoring Enables Reproducible Quantification of Phospho-Signaling. Mol Cell Proteomics (2016) 15(2):726–39. doi: 10.1074/mcp.O115.054940
99. Deutsch EW, Bandeira N, Sharma V, Perez-Riverol Y, Carver JJ, Kundu DJ, et al. The ProteomeXchange Consortium in 2020: Enabling ’Big Data’ Approaches in Proteomics. Nucleic Acids Res (2020) 48(D1):D1145–52. doi: 10.1093/nar/gkz984
100. Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE Database and Related Tools and Resources in 2019: Improving Support for Quantification Data. Nucleic Acids Res (2019) 47(D1):D442–50. doi: 10.1093/nar/gky1106
101. Sherry ST, Ward M, Sirotkin K. dbSNP-Database for Single Nucleotide Polymorphisms and Other Classes of Minor Genetic Variation. Genome Res (1999) 9(8):677–9. doi: 10.1101/gr.9.8.677
102. Hornbeck PV, Kornhauser JM, Tkachev S, Zhang B, Skrzypek E, Murray B, et al. PhosphoSitePlus: A Comprehensive Resource for Investigating the Structure and Function of Experimentally Determined Post-Translational Modifications in Man and Mouse. Nucleic Acids Res (2012) 40(Database issue):D261–70. doi: 10.1093/nar/gkr1122
103. Schoenherr RM, Huang D, Voytovich UJ, Ivey RG, Kennedy JJ, Saul RG, et al. A Dataset Describing a Suite of Novel Antibody Reagents for the RAS Signaling Network. Sci Data (2019) 6(1):160. doi: 10.1038/s41597-019-0166-7
104. Addona TA, Abbatiello SE, Schilling B, Skates SJ, Mani DR, Bunk DM, et al. Multi-Site Assessment of the Precision and Reproducibility of Multiple Reaction Monitoring-Based Measurements of Proteins in Plasma. Nat Biotechnol (2009) 27(7):633–41. doi: 10.1038/nbt.1546
105. Glatter T, Ludwig C, Ahrne E, Aebersold R, Heck AJ, Schmidt A. Large-Scale Quantitative Assessment of Different in-Solution Protein Digestion Protocols Reveals Superior Cleavage Efficiency of Tandem Lys-C/trypsin Proteolysis Over Trypsin Digestion. J Proteome Res (2012) 11(11):5145–56. doi: 10.1021/pr300273g
106. Carr SA, Abbatiello SE, Ackermann BL, Borchers C, Domon B, Deutsch EW, et al. Targeted Peptide Measurements in Biology and Medicine: Best Practices for Mass Spectrometry-Based Assay Development Using a Fit-for-Purpose Approach. Mol Cell Proteomics (2014) 13(3):907–17. doi: 10.1074/mcp.M113.036095
107. Dou Y, Kawaler EA, Cui Zhou D, Gritsenko MA, Huang C, Blumenberg L, et al. Proteogenomic Characterization of Endometrial Carcinoma. Cell (2020) 180(4):729–748 e26. doi: 10.1016/j.cell.2020.01.026
108. Zhang H, Liu AY, Loriaux P, Wollscheid B, Zhou Y, Watts JD, et al. Mass Spectrometric Detection of Tissue Proteins in Plasma. Mol Cell Proteomics (2007) 6(1):64–71. doi: 10.1074/mcp.M600160-MCP200
109. Schoenherr RM, Saul RG, Whiteaker JR, Yan P, Whiteley GR, Paulovich AG. Anti-Peptide Monoclonal Antibodies Generated for Immuno-Multiple Reaction Monitoring-Mass Spectrometry Assays Have a High Probability of Supporting Western Blot and ELISA. Mol Cell Proteomics (2015) 14(2):382–98. doi: 10.1074/mcp.O114.043133
110. Whiteaker JR, Sharma K, Hoffman MA, Kuhn E, Zhao L, Cocco AR, et al. Targeted Mass Spectrometry-Based Assays Enable Multiplex Quantification of Receptor Tyrosine Kinase, MAP Kinase, and AKT Signaling. Cell Rep Methods (2021) 1(3):1–11. doi: 10.1016/j.crmeth.2021.100015
Keywords: immunotherapy, cancer, correlative biomarkers, mass spectrometry, immuno-MRM
Citation: Whiteaker JR, Lundeen RA, Zhao L, Schoenherr RM, Burian A, Huang D, Voytovich U, Wang T, Kennedy JJ, Ivey RG, Lin C, Murillo OD, Lorentzen TD, Thiagarajan M, Colantonio S, Caceres TW, Roberts RR, Knotts JG, Reading JJ, Kaczmarczyk JA, Richardson CW, Garcia-Buntley SS, Bocik W, Hewitt SM, Murray KE, Do N, Brophy M, Wilz SW, Yu H, Ajjarapu S, Boja E, Hiltke T, Rodriguez H and Paulovich AG (2021) Targeted Mass Spectrometry Enables Multiplexed Quantification of Immunomodulatory Proteins in Clinical Biospecimens. Front. Immunol. 12:765898. doi: 10.3389/fimmu.2021.765898
Received: 27 August 2021; Accepted: 22 October 2021;
Published: 11 November 2021.
Edited by:
Subramaniam Malarkannan, Medical College of Wisconsin, United StatesReviewed by:
Christoph Borchers, McGill University, CanadaArthur Salomon, Brown University, United States
Copyright © 2021 Whiteaker, Lundeen, Zhao, Schoenherr, Burian, Huang, Voytovich, Wang, Kennedy, Ivey, Lin, Murillo, Lorentzen, Thiagarajan, Colantonio, Caceres, Roberts, Knotts, Reading, Kaczmarczyk, Richardson, Garcia-Buntley, Bocik, Hewitt, Murray, Do, Brophy, Wilz, Yu, Ajjarapu, Boja, Hiltke, Rodriguez and Paulovich. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Amanda G. Paulovich, apaulovi@fredhutch.org