Standardized Profiling of The Membrane-Enriched Proteome of Mouse Dorsal Root Ganglia (DRG) Provides Novel Insights Into Chronic Pain*

Chronic pain is a complex disease with limited treatment options. Several profiling efforts have been employed with the aim to dissect its molecular underpinnings. However, generated results are often inconsistent and nonoverlapping, which is largely because of inherent technical constraints. Emerging data-independent acquisition (DIA)-mass spectrometry (MS) has the potential to provide unbiased, reproducible and quantitative proteome maps - a prerequisite for standardization among experiments. Here, we designed a DIA-based proteomics workflow to profile changes in the abundance of dorsal root ganglia (DRG) proteins in two mouse models of chronic pain, inflammatory and neuropathic. We generated a DRG-specific spectral library containing 3067 DRG proteins, which enables their standardized quantification by means of DIA-MS in any laboratory. Using this resource, we profiled 2526 DRG proteins in each biological replicate of both chronic pain models and respective controls with unprecedented reproducibility. We detected numerous differentially regulated proteins, the majority of which exhibited pain model-specificity. Our approach recapitulates known biology and discovers dozens of proteins that have not been characterized in the somatosensory system before. Functional validation experiments and analysis of mouse pain behaviors demonstrate that indeed meaningful protein alterations were discovered. These results illustrate how the application of DIA-MS can open new avenues to achieve the long-awaited standardization in the molecular dissection of pathologies of the somatosensory system. Therefore, our findings provide a valuable framework to qualitatively extend our understanding of chronic pain and somatosensation.

Although acute pain is an evolutionary adaptive response, chronic pain is considered a pathology of the nervous system (1,2). Treatment of chronic pain presents a big challenge as current therapies target molecules with key physiological functions throughout the body resulting in minimal safety profiles. Hence, the identification of proteins specifically involved in chronic pain states is of paramount importance for designing optimized interventions.
In vertebrates, noxious insults are detected by specialized peripheral sensory neurons (so-called nociceptors) such as those of dorsal root ganglia (DRG) 1 . These DRG neurons innervate the skin, muscles and all organs where they serve as primary cellular detectors of noxious stimuli of various qualities (chemical, thermal or mechanical) including tissue injury and inflammation (3,4). Ultimately, their activity shapes the afferent message transmitted to the brain where the sensation of pain is generated. The function of nociceptors relies on dedicated protein networks (3,4). During the last decade enormous progress has been made regarding the identification of ion channels, receptors and signaling proteins involved in nociception and pain (4 -6). However, a systematic and quantitative comparison of protein dynamics across different pain states remains a major challenge.
Several -omics approaches have been developed to study pain-related changes in gene expression (7,8) or of proteins (9 -14). Deep-sequencing transcriptomics (RNA-Seq) is able to quantify the majority of the transcriptome (7,8), but only half of mRNA changes account for protein variability (15) limiting the ability to relate results to protein function. On the other hand, proteome profiling schemes using data-dependent acquisition (DDA; commonly referred to as shotgun mass spectrometry) are semi-stochastic and biased toward detecting the most abundant proteins in complex samples, all of which leads to well-documented effects of undersampling (16 -19). Consequently, an average shotgun mass spectrometry (MS) experiment fails to detect 84% of peptides present in a given sample (16). Furthermore, as many as 30 -60% of sampled peptides may differ between technical replicates because of limited reproducibility (20) precluding reliable proteome profiling. Besides other factors related to experimental design, these technical limitations are largely responsible for the fact that the pain research community has accumulated inconsistent and nonoverlapping lists of regulated proteins during pain (9 -14).
A potential solution might be offered by emerging alternative mass spectrometry-based methods such as data-independent acquisition mass spectrometry (DIA-MS) (21)(22)(23). In contrast to random selection and fragmentation of only the most abundant precursor ions in shotgun approaches, one major hallmark of DIA-MS lies in the concurrent selection and fragmentation of nearly all precursor ions irrespective of their abundance (as long as their abundance is above the detection limit of the mass spectrometer). Precursor ion selection is solely dependent on user-defined mass windows spanning several hundreds of mass-to-charge ratios (m/z). Consequently, one implementation of this technique is called Sequential Windowed Acquisition of all THeoretical mass spectra (SWATH) (21)(22)(23). However, recorded fragment ion spectra are of chimeric nature as they may derive from several coeluting precursor ions, which does not permit the application of classical database searches to decipher protein identities. Instead, a prevalent approach for DIA-MS data analysis is the use of reference spectral libraries containing information on retention time and fragmentation spectra for all included peptides (21)(22)(23). Thus, DIA-MS is capable to quantify thousands of proteins present in a reference spectral library across multiple samples in a highly reproducible and experimenter-independent fashion (21)(22)(23). Ergo, this technique is especially suited to globally interrogate pathology-related proteome alterations (23)(24)(25); however, DIA-MS has so far not been applied to decipher the protein signature of chronic pain.
Here, we present a multi-layer workflow consisting of mouse behavior, biochemistry and DIA-MS for standardized detection and quantification of the changes in the DRG proteome in two mouse models of chronic pain. We reproducibly quantified the same set of 2526 proteins in all biological and technical replicates, and revealed 141 differentially regulated candidates, of which dozens have not been characterized in the somatosensory system. The biological significance of our data is reflected by our functional validation experiments suggesting the involvement of selected candidates in nociceptive signaling and pain. In summary, our extensive data set establishes a valuable framework to qualitatively extend our understanding of chronic pain and somatosensation.

EXPERIMENTAL PROCEDURES
Animals-Adult C57Bl/6J mice (6 -8 weeks old at the time of surgery or injection) were in-house bred and kept in a temperatureand humidity-controlled environment under a 12 h (h) light/dark cycle with food and water provided ad libitum. All animal experiments were approved and carried out in compliance with institutional guidelines and guidelines of the Landesamt fü r Verbraucherschutz und Lebensmittelsicherheit of Lower Saxony, Germany (33.9-42502-04-14/1638 and 3392-42502-04-13/1077).
Pain Paradigms and Measurement of Mechanical Hypersensitivity-Chronic Inflammatory Pain-Ten microliters of Complete Freund's Adjuvant (CFA, 1.0 mg/ml, Sigma-Aldrich, Carlsbad, CA) were injected into the plantar surface of one hind paw. Control mice were injected with 10 l of phosphate-buffered saline (PBS; Veh; Life Technologies). Mechanical hypersensitivity was measured with a dynamic aesthesiometer (Ugo Basile) 24 h after injection according to standard procedures (26) and the manufacturerЈs manual (Ugo Basile, Monvalle VA, Italy). In brief, mice were acclimatized for 2 h and paw withdrawal latencies were determined upon application of a graded force (0 -10 g in 40 s) to the plantar surface of each hind paw. Paw withdrawal latencies were averaged from at least four readings per paw. In order to minimize stress (and associated potential alterations of the proteome) because of behavioral testing (27,28), we intended to separate behavioral testing as much as possible from tissue collection in time. Therefore, we did not test mechanical hypersensitivity in CFA-and Veh-injected mice destined for tissue collection for subsequent experiments (calcium imaging, biochemistry and mass spectrometry). Instead, we used a separate cohort of mice to assess mechanical hypersensitivity. This separate mouse cohort was CFAinjected in parallel (control: Veh injection) and subsequently tested for mechanical hypersensitivity. Separate mouse cohorts consisted of at least 7 animals per DIA-MS replicate. All of them exhibited pronounced mechanical hypersensitivity as described previously (26,29) and shown in Fig. 1B. Moreover, all CFA-injected animals in this study showed pronounced paw edema (26), which was absent in Vehinjected mice.
Neuropathic Pain-The spared nerve injury (SNI) paradigm was induced according to standard protocols (30). Mice were anesthetized with isoflurane (4% induction, 1.8% maintenance) and the sciatic nerve and its branches, the common peroneal, tibial and sural nerves, were exposed. Common peroneal and tibial nerves were ligated with a 6.0 silk suture (Braun) and transected distal to the ligation, removing 2-3 mm of each nerve. Care was taken to avoid damage to the sural nerve. Sham surgery was performed in a similar manner without ligating and transecting the nerves. Mechanical sensitivity was assessed 2 days before and 7, 14, 21, and 26 days after surgery, as described above, but applying the force to the lateral side of the plantar surface of both hind paws. Only mice that exhibited evident mechanical hypersensitivity (manifested as reduced withdrawal latencies) at all tested time points (90% of all mice receiving SNI) were included in the study and used for tissue collection 28 days after surgery. As explained above, tissue collection (at day 28 post surgery) was separated in time from behavioral testing (at day 26 post surgery) for the sake of stress-reduction.
Injection of Rotenone and PACMA 31-Ten microliters of CFA was unilaterally injected into the plantar surface of a hindpaw of mice. After 24 h, mice were acclimatized for 1-2 h, and injected into the plantar surface of the same hind paw with either 5 l of Veh (0.63% DMSO/PBS) or a suspension of 5 l of 2.5 mM Rotenone (Sigma-Aldrich) in 0.63% DMSO/PBS. Mice were acclimatized for another 15 min (min) and mechanical sensitivity determined by a plantar aesthesiometer (Ugo Basile) as described above up to 1 h after injection of Rotenone/Veh. Maximally 200 l of PACMA 31 (Tocris, Bristol, UK; 20 mg/kg body weight) or Veh (6.9% DMSO in sterile PBS) were injected intraperitoneally. Mice were then acclimatized for 2 h and mechanical sensitivity determined by a plantar aesthesiometer (Ugo Basile) as described above. Withdrawal latencies were measured up to 6 h after injection of PACMA 31/Veh. At least two independent mouse cohorts were tested for each agent. In line with other studies (31)(32)(33), neither Rotenone nor PACMA 31 affected general health or motor coordination.
Membrane-enriched Fractions of DRG Lysates-Chemicals were obtained from Roth, unless stated otherwise. 24 h after induction of inflammatory pain and 28 days after induction of neuropathic pain mice were euthanized with CO 2 , decapitated, ipsilateral lumbar DRG were dissected, immediately snap-frozen in liquid nitrogen and stored at Ϫ80°C until further use. To prepare membrane-enriched fractions we used a protocol modified from Lu and colleagues (34). Lumbar DRG of 10 -13 mice per condition were homogenized in high salt buffer (2 M NaCl, 10 mM Tris/HCl, pH 7.5, complete protease inhibitor mixture, Roche), and centrifuged at 600 ϫ g for 10 min, 4°C. Supernatants were centrifuged at 16,000 ϫ g for 15 min, 4°C, and the pellets resuspended twice in 100 l carbonate buffer (0.1 mM Na 2 CO 3 , pH 11.3), followed by centrifugation at 16,000 ϫ g for 15 min, 4°C. Pellets were then resuspended in 100 l of 1% RapiGest (Waters, Eschborn, Germany) in ammonium bicarbonate buffer (50 mM (NH 4 )HCO 3 , pH 8.5), incubated overnight, 4°C, and centrifuged at 16,000 ϫ g for 15 min, 4°C. For acetone-precipitation, 500 l of ice-cold acetone were added to the supernatants and samples were incubated for 2 h, Ϫ20°C, and centrifuged at 14,000 ϫ g for 30 min. Pellets were washed with ice-cold 80% EtOH and centrifuged at 14,000 ϫ g for 30 min. Supernatants were taken off and pellets air-dried. For the in-solution digest, 10 l of 1% RapiGest in ammonium bicarbonate buffer were added followed by 10 l of 50 mM dithiothreitol (Calbiochem) in 25 mM ammonium bicarbonate buffer, and samples were incubated at 850 rpm for 1 h, 56°C on a thermoshaker. Next, 10 l of 100 mM iodoacetamide (Sigma-Aldrich, Munich, Germany) in 25 mM ammonium bicarbonate buffer were added and samples were incubated at 850 rpm for 1 h, 37°C. Seventy microliters of trypsin (Promega, Mannheim, Germany) solution in 25 mM ammonium bicarbonate buffer were added (1:20 enzyme to protein ratio) and samples were incubated overnight, 37°C. Twenty microliters of 5% trifluoroacetic acid (Roth, Mannheim, Germany) were added, samples were incubated for 2 h, 37°C, and centrifuged at 16000 ϫ g for 30 min, room temperature. Supernatants were lyophilized and stored at Ϫ20°C until further analysis.
Sample Preparation, DIA-MS and Data Analysis-All solvents were HPLC-grade from Sigma-Aldrich and all chemicals were obtained from Sigma-Aldrich if not stated otherwise. All steps of DIA-MS and its analysis were performed by Biognosys AG (Zuerich, Switzerland) essentially as described (23); in short: Sample Preparation-The 12 samples (3 biological replicates per experimental group) were prepared for mass spectrometric analysis. The peptides were re-suspended in 8 M urea/0.1 M ammonium bicarbonate at room temperature (RT). Peptides were cleaned up for mass spectrometric analysis using a C18 Ultra MicroSpin Column (The Nest Group, Southborough, MA) according to the manufacturer's instructions. The peptide concentration was determined using a Spectro-starNano spectrometer and was in the expected range of 1 g/l for all samples. Peptides were dried down to complete dryness and resuspended in liquid chromatography (LC) solvent A (40 l of 1% acetonitrile in water with 0.1% formic acid).
Following, samples were analyzed for spectral library generation and Hyper Reaction Monitoring (HRM)-MS profiling. The LC gradient was 5-35% solvent B in 120 min followed by 35-100% B (3% water in acetonitrile with 0.1% formic acid) in 2 min and 100% B for 8 min (total gradient length was 240 min). For shotgun analyses on the Q Exactive a modified TOP12 method from Olsen et al. was used (35). Prior to mass spectrometric acquisition, the HRM Calibration Kit (Biognosys) was spiked into the samples at the manufacturer's instructions (36).
Spectral Library Generation-To generate the spectral library, we used one sample for each of the four experimental groups and pooled the results. Thus, 4 h gradient LC-MS/MS acquisitions were performed in shotgun proteomics mode using a Q Exactive mass spectrometer. The LC-MS/MS mass spectrometric data were analyzed using Maxquant 1.4 software (37) with the following settings: Trypsin/P as protease; maximal two missed cleavages; carbamidomethyl (C) as fixed modification; oxidation (M) and acetylation (at protein N-term) as variable modifications; false discovery rate on peptide and protein level set to 1%; mass tolerance for precursor ions 20 ppm and fragment ions (MS/MS) 20 ppm. A mouse Uniprot FASTA database (state 05.06.2014) was used for the search engine. The final spectral library contained 16,971 peptide assays including charge states and modifications (of which 15,850 were proteotypic) corresponding to 3067 proteins (of which 2530 were defined by at least one proteotypic peptide). On average, there were 5.5 peptides per protein.
DIA-MS and Data Analysis-The 12 samples of the four conditions and three replicates were acquired in HRM mode in a replicate blocking manner to avoid performance bias in the acquisition (supplemental Fig. S3). The HRM method consisted of loops of 1 survey scan and 19 SWATHs (i.e. isolation windows) spanning 400 to 1200 m/z. The LC gradient was 2 h. This resulted in on average 8 data points per peak, sufficient for quantitative analysis. The 12 HRM measurements were analyzed using the Spectronaut™ 5.0 software from Biognosys AG (Zuerich, Switzerland). The spectral library generated in this project was used for the targeted analysis approach. The false discovery rate (FDR) on peptide level was set to 1%. The resulting data were filtered using row based extraction, for a peptide identified in one measurement; all quantitative values were extracted in all measurements resulting in full profiles. The resulting data tables were normalized using local normalization controlling for loading differences and spray bias (supplemental Fig. S4).
Quality control was performed by a control sample (human blood plasma) measured before and after profiling MS measurements in shotgun proteomics. The analysis of the two control measurements revealed that the performance stayed constant over the time of acquisition (data not shown). To identify significantly differential expressed proteins, a state comparison analysis was performed comparing either CFA versus Veh or SNI versus Sham using a pairwise t test. The resulting p values were corrected for multiple testing using the Benjamini-Hochberg and q-value method to control the overall false discovery rates (38). State comparison analysis and the analysis of the detection of fold change were performed using MSstats 2.3.2 (39) essentially as described (23).
All data have been deposited to Peptide Atlas (http://www.peptideatlas. org, No. PASS00826, the username is PASS00826 and the password is ZE5945at).
The metadata analysis to compare all 141 regulated candidates with proteins and genes previously shown to be associated with pain states in humans and rodents was done manually by extensive literature search using Pubmed (NCBI) and the web-resource Pain Networks (painnetworks.org) (40). For the purpose of determining the amount of novel proteins identified in our approach (i.e. proteins that have not been investigated in the context of pain before), we disregarded all genes/proteins that (1) were shown to be involved in nociception and pain elsewhere, (2) are included in Pain Networks as pain-related genes or (3) are considered members of the mitochondrial electron transport chain (ETC) as inhibitors of mitochondrial ETC have been demonstrated to be involved in nociception and pain (32,41). Gene ontology (GO) analysis was performed by uploading mouse gene IDs into the search GO data engine (42). To determine enriched GO terms we used the DAVID functional annotation tool to assign genes with their affiliate terms and to order them by enrichment (43,44). Only significantly enriched GO terms (p Ͻ 0.05 after Benjamini-Hochberg correction) are reported. To visualize molecular interaction networks and predict direct and indirect associations among all significantly regulated proteins in each pain model we uploaded accession numbers (Uniprot; mus musculus) of differentially regulated proteins into STRING (Search Tool for the Retrieval of Interacting Genes/ Proteins) (45). STRING settings: Confidence view; confidence level 0.4; clustering algorithm MCL set to 2. Of note, we only display "connected proteins," i.e. associations among significantly regulated proteins (Supplemental Table 3). As disconnected nodes are hidden for the purpose of visualizing molecular interaction networks, not all significantly regulated proteins are included in the presented protein networks (Fig. 4).
Western Blots-Membrane-enriched fractions of lumbar DRG of 12 mice per condition were obtained as described above. Instead of resuspending pellets in RapiGest (Waters), pellets were resuspended in 30 l PBS/lithium dodecyl sulfate sample buffer (Life Technologies) and incubated for 1.5 h, 4°C, while shaking. 1ϫ NuPAGE sample reducing agent (Life Technologies) was added and equal amounts of protein separated by 1D NuPAGE. Protein gels of membrane-enriched fractions were blotted on PVDF membranes using the iBlot System (Life Technologies) and blocked for 30 min in 5% milk/PBS. Western blots were incubated for 1 h, RT, with primary antibodies to ␤-Actin (mouse, 1:500, Sigma-Aldrich #A1978), to Annexin A2 (rabbit; 1:100; Abcam #ab41803, Cambridge, UK), to sodium-potassium ATPase alpha2 (rabbit, 1:100, Millipore, Darmstadt, Germany #07-674) and Serca1 (rabbit, 1:100, Alomone, Jerusalem, Israel #ACP011) in 1% milk/PBS, and additionally overnight, 4°C, at a three times lower concentration. Blots were incubated for 2 h at RT with secondary antibodies conjugated to Alexa Fluor 680 or 790 (1:8000, Life Technologies, #A10038/A10043) in 1% milk/PBS, washed extensively and imaged using infrared imaging (Odyssey, LI-COR, Bad Homburg, Germany). For our orthogonal validation experiments (Fig. 3) Western blots were cut to be able to separately probe for high (Serca1) and low molecular weight proteins (␤-Actin) (please see supplemental Fig. S7 for full blots). Western blots were repeated two times. Only for presentation purposes size, brightness, and contrast levels of paired blots (i.e. experiments done in parallel and probed with the same antibody) were adjusted using Photoshop CS6 (Adobe).
Image Acquisition and Analysis of Immunohistochemistry-Stained cryosections were imaged at a Zeiss Axio Observer Z1 inverted epifluorescence microscope. Images for all experimental groups were taken using identical acquisition parameters. For each experimental group, at least six pictures were randomly analyzed. Raw images were analyzed using NIH ImageJ (46). Signal was considered positive when higher than mean background fluorescence plus three times the standard deviation, measured from at least 10 random unstained tissue spots. Mean background fluorescence intensity was subtracted from the mean surface fluorescence intensity. Signal particle size threshold was set at 15 square pixels to account for nonspecific staining. Values were averaged among experiments. As the SNI model involves transection/injury of two major branches of the sciatic nerve (30) we separately analyzed big fascicles (corresponding to branches of the tibial and common peroneal nerves) and small fascicles (corresponding to the sural nerve and smaller branches of the tibial and common peroneal nerves, e.g. articular branch) in cross sections of sciatic nerves (47). To be consistent, we applied the same procedure to cross sections obtained from Sham animals as well as CFA-and Veh-injected animals. The analysis of fluorescence intensity was only performed on cross sections. Only for presentation purposes brightness and contrast of paired images were adjusted using Photoshop CS6 (Adobe).
Experiments were performed 24 h after plating. Assay buffer consisted of regular Tyrodes's (in mM): NaCl 140, KCl 4, CaCl 2 2, Glucose 10, MgCl 2 2, HEPES 10, with an osmolality of 295-300 mOsm and pH 7.4. Neuronal cultures were incubated with 2.5 M Fura2-AM (Life Technologies) and 0.04% Pluronic F-127 (Life Technologies) in Tyrode's, for 30 min, 37°C, washed and mounted on a recording chamber constantly superfused by a gravity-driven system. To block PMCA activity, HEPES was replaced by 10 mM Trizma and pH adjusted to 8.8 (49). Block of Serca activity was obtained by further addition of 150 nM Thapsigargin (TG; Sigma-Aldrich). At this concentration TG has been shown to be specific for Serca while not affecting voltage-gated calcium currents in mouse DRG neurons (50). The fluorophore was excited alternately at 340 nm and 380 nm, and emission acquired at 510 nm, using the Zeiss Axio Observer Z1. Intracellular calcium concentration is expressed as the 340/380 ratio. All experimental groups to be compared were processed in parallel using the same culture preparation. At least two coverslips from three independent culture preparations were analyzed as described (48). The mean intensity value was assessed with ImageJ following background subtraction and particle analysis as described above.
Experimental Design and Statistical Rationale-All DIA-MS experiments were performed in triplicates, i.e. 3 biological replicates were analyzed. Experimental conditions were CFA-treated animals and animals with spared nerve injury (SNI). Control conditions involved Veh-treated animals (control for CFA-treated animals) and animals receiving Sham surgery (control for SNI). DIA-MS experiments involved pooled tissue from 10 -13 animals/replicate/condition. To identify significantly differential expressed proteins, a state comparison analysis was performed comparing either CFA versus Veh or SNI versus Sham using a pairwise t test. The resulting p values were corrected for multiple testing using the Benjamini-Hochberg and qvalue method to control the overall false discovery rates. State comparison analysis and the analysis of the detection of fold change were performed using MSstats 2.3.2 (39) essentially as described (23).
For statistical analysis of all other experiments, two-tailed unpaired t test was used for single comparisons and one-way ANOVA was used for multiple comparisons. Statistical analysis was performed with GraphPad Prism if not stated otherwise. No statistical methods were used to predetermine sample sizes for any experiment, but our sample sizes are similar to those generally employed in the field. Compiled data are expressed as mean Ϯ sem; in Fig. Legends n indicates the sample number, N indicates the number of independent experiments, P denotes the significance (*, p Ͻ 0.05, **, p Ͻ 0.01, ***, p Ͻ 0.001) and refers to the respective control in each experimental group if not noted otherwise, and NS indicates "not significant."

Overview of Our Integrated Workflow: From Murine Chronic
Pain Models to DIA-MS Proteome Profiling-In order to study proteins relevant for somatosensory signaling by DIA-MS, we generated a customized spectral library (22,51,52) from DRG isolated under pain and control conditions, respectively. To this end, we employed two well-established mouse models of chronic pain with different etiology (26,53) (Fig. 1A): inflammatory pain induced by injection of Complete Freund's Adjuvant (CFA; Vehicle (Veh) injections served as controls) and neuropathic pain using the spared nerve injury model, in which the tibial and common peroneal branch of the sciatic nerve are transected (SNI; control groups received Sham surgery; please refer to Experimental Procedures for further details). Both models mimic clinically relevant features of chronic pain patients including mechanical hypersensitivity (54,55). Therefore, we assessed hypersensitivity to innocuous punctuate mechanical stimulation (allodynia) by the plantar aesthesiometer test (26) in treated (ipsilateral) hindpaws compared with the contralateral side as well as to Veh and Sham surgery controls, respectively ( Fig. 1B and supplemental Fig.  S1). Only mouse cohorts showing prominent mechanical allodynia were included in our study ( Fig. 1B; please refer to Experimental Procedures for further details). Animals were FIG. 1. Experimental workflow. A, Representation of the four behavioral groups used in this study: (from left to right) Vehicle (Veh) and CFA injection (CFA) into one hind paw, exposed but not transected (Sham) and exposed as well as transected branches of the sciatic nerve (SNI). B, Representative assessment of hypersensitivity to mechanical stimulation reveals pronounced mechanical allodynia upon induction of inflammatory and neuropathic pain. Mechanical allodynia is manifested by decreased paw withdrawal latencies of the treated side (ipsilateral; Ipsi) in CFA and SNI animals (24 h after CFA-injection and 28 days after SNI). Paw withdrawal latencies were unaltered in Veh-injected (control for CFA-injection) and Sham animals (control for SNI) and contralateral paws (Contra) of all groups (***, p Ͻ 0.001 by one-way ANOVA followed by BonferroniЈs multiple comparison test; n ϭ 21-24 animals/condition from three independent cohorts). C, Representative Western blot demonstrates enrichment for the membrane marker Na ϩ /K ϩ -ATPase and decrease of the cytosolic markers Actin and AnnexinA2 (AnxA2) in membrane-enriched fractions (S4). As expected, opposing results were obtained in cytosolic-enriched fractions (S1 and S2; n ϭ 3 independent preparations). D, Membrane-enriched samples (S4) were analyzed by DIA-MS. sacrificed for dissection of ipsilateral lumbar DRG (lDRG) 24 h after induction of inflammatory pain and 28 days after induction of neuropathic pain (and respective controls). The time points were chosen based on previous publications reporting robust mechanical allodynia (26,29,30) and our assessment of mechanical hypersensitivity (Fig. 1B and supplemental Fig. S1).
The function and interplay of proteins expressed in sensory neurons and other cell types of DRG (e.g. peripheral glial cells and various immune cells) shape nociceptive signaling (56,57). Especially transmembrane and membrane-associated proteins (hereafter summarized as "membrane proteins") are critically involved in stimulus detection, neuronal activity and consequently, the sensation of pain (4 -6, 58 -61). Hence, they represent attractive targets for pharmacotherapy (62,63). Yet, the experimental analysis of membrane proteins is hampered by their low abundance and pronounced hydrophobicity (64). As a corollary to these difficulties, the membrane proteome is traditionally understudied despite its paramount importance for diverse (patho-) physiological processes including pain (4-6, 58-61, 63). For these reasons we aimed at assessing the potential of DIA-MS to profile a membrane-enriched proteome of DRG. To this end we performed biochemical sub-fractionation and obtained membrane-enriched lDRG lysates of each experimental and control condition (Fig. 1C). These were subsequently used for DIA-MS-based profiling (Fig. 1D).
A DRG-specific Spectral Library-We prepared three independent biological replicates per pain paradigm (CFA versus Veh and SNI versus Sham) each consisting of membraneenriched lDRG lysates from 10 -13 mice. One replicate per condition was used to generate the spectral library on a Q Exactive mass spectrometer. Data were analyzed using Maxquant 1.4 software (37), setting the false discovery rate on peptide and protein level to 1%. A mouse Uniprot FASTA database (state 05.06.2014) was used for the search engine. The final spectral library contained 16,971 peptide assays (of which 15,850 were proteotypic) corresponding to 3067 proteins (of which 2530 were defined by at least one proteotypic peptide) (supplemental Table S1). To examine the representation of cellular compartments and functions in our spectral library, we systematically analyzed gene ontology (GO) terms with the Database for Annotation, Visualization and Integrated Discovery (DAVID; david.ncifcrf.gov). Spectral library components broadly covered diverse cellular compartments and molecular functions (supplemental Fig. S2). As expected, this analysis confirmed the enrichment of membrane-associated proteins (including membrane-passing, i.e. transmembrane proteins) in our spectral library (supplemental Fig. S2). Among the latter were several ion channels with known roles in somatosensation and pain, most of which have not been identified in a previous proteomic study of DRG membranes (65), e.g. Piezo2 and TRPA1 channels (supplemental Table S2). Further comparison with published proteomes of neuromas, sciatic nerve, DRG, and spinal cord (harboring presynaptic terminals of DRG neurons) demonstrates, that our spectral library constitutes the largest compendium of DRG proteins reported to date (9,10,(12)(13)(14)65). The generated resource is publically available in PeptideAtlas (66) and constitutes an optimized resort for unbiased, consistent and highly accurate protein quantification in future DIA-MS experiments of any mouse sample of interest.
Differential Regulation of Proteins During Inflammatory and Neuropathic Pain-We next used the generated spectral library for the targeted search of DIA-MS runs conducted with three biological replicates of each of the four conditions: induction of inflammatory pain by CFA versus Veh-injection and induction of neuropathic pain by SNI versus Sham (Fig. 1). The 12 samples were acquired in Hyper Reaction Monitoring (HRM) mode in a replicate blocking manner to avoid performance bias in the acquisition (23), which resulted in high quality raw data acquisition (supplemental Fig. S3). The spectral library was applied setting the false discovery rate (FDR) of detection at 1% on a peptide level. This analysis enabled the profiling of more than 14,500 peptides representing the same set of 2526 proteins across all replicates and samples (2581 proteins in all CFA and Veh replicates; 2600 proteins in all SNI and Sham replicates; supplemental Table S3), which reflects the high reproducibility of the presented workflow ( Fig. 2A). The resulting data tables were normalized using local normalization in order to account for loading differences and spray bias (67) (supplemental Fig. S4). The analysis of the coefficient of variation of the quantified signals (which determines the variance from the biological pool, sample preparation and instrument performance) was in the expected range below 20% (68) indicating conserved and reliable quantification performance (supplemental Fig. S5).
Upon statistical comparison of these data (corrected for multiple testing using the q-value (38); for complete data set please see supplemental Table S3), we observed little overall proteome changes as visualized in the heatmap and volcano plots (Fig. 2). Only around 3% of all quantified proteins showed significant up-and downregulations in each pain paradigm (64 and 77 proteins in inflammatory and neuropathic pain, respectively; Fig. 2A) and 12 proteins were commonly altered in both pain models, suggesting a high degree of pain model-specificity ( Fig. 2A; supplemental Table S4). Gene Ontology (GO) analysis revealed that regulated proteins were expressed in different cellular compartments (Fig. 2B,  2C) and covered distinct as well as overlapping molecular functions (supplemental Fig. S6; supplemental Table S4). Intriguingly, several significantly regulated candidates of our data set have previously been associated with diverse pathologies involving the somatosensory system (including painful ones) in humans and rodents (examples are given in Table I), demonstrating the ability of our workflow to uncover meaningful candidates. For instance, Dipeptidylpeptidase 4 (Dpp4) was significantly up-regulated in both pain models (Table I  and supplemental Tables S3 and S4). Dpp4 is a multifunc- tional cell surface glycoprotein with peptidase activity controlling the activity of numerous peptides and proteins (69). Interestingly, Dpp4 Ϫ/Ϫ mice show enhanced nociceptive pain (70), whereas in vivo inhibition of Dpp4 produces analgesia in chronic diabetic neuropathy in rats (71). Similarly, Heat shock protein 90b1 (HSP90b1), commonly regulated in both pain models of our study (Table I), is implicated in diabetic neuropathy in rodents (72,73). Other candidates in our study are involved in painful human pathologies such as familial hemiplegic migraine type II caused by mutations of the ATPase ATP1A2 (ATP1A2 is commonly regulated in our study, Table I) (74) or Charcot-Marie-Tooth 4A neuropathy linked to mutations of Ganglioside-induced differentiation associated-protein 1 (Gdap1) (Gdap1 is regulated during inflammatory pain in our study, Table I) (75).
Conversely, a detailed literature search revealed that roughly 50% of significantly altered proteins identified by our approach have not been investigated in the context of nociception or pain before (please see Experimental Procedures for details on literature search). Hence, they represent potentially promising candidates for future studies on their role in painful conditions.
Volcano plot for all proteins profiled in inflammatory (CFA/Veh) pain. Dotted horizontal line represents the significance level, i.e. q-value ϭ 0.05. Hence, all proteins above this dotted line exhibit a q-value Ͻ 0.05 and are considered to be significantly regulated. Right, Assigned Gene Ontology (GO) Cellular Component (CC) categories and the relative abundance of significantly regulated proteins during inflammatory pain. C, Left, Volcano plot for all proteins profiled in neuropathic (SNI/Sham) pain. Significance level q-value Ͻ 0.05 as in B. Right, Assigned Gene Ontology (GO) Cellular Component (CC) categories and the relative abundance of significantly regulated proteins during neuropathic pain.

TABLE I Examples of significantly regulated proteins in our study, which have been reported to be implicated in rodent pain models or vertebrate pathologies involving the somatosensory system. Of note, this table contains only examples of such proteins for which functional significance was demonstrated in vivo. In our study, blue-labeled proteins were significantly regulated during inflammatory pain, red-labeled proteins during neuropathic pain and green-labeled proteins in both pain models. dHMNs: distal hereditary motor neuropathies
Correlation of proteomic changes with altered protein function and expression-Next, we intended to validate our results by orthogonal methods. To this end, we focused on sarco/ endoplasmic reticulum calcium-ATPase (Serca) as it was differentially regulated in the two pain models we used for this study, i.e. Serca was up-regulated in CFA and downregulated in SNI animals (supplemental Table S4). Using independent mouse cohorts, we induced inflammatory and neuropathic pain and compared the abundance of Serca in DRG by Western blotting (Fig. 3A) as well as by immunohistochemistry (Fig.  3B-3D). Of note, animals were sacrificed for tissue dissection at the same time points as for the DIA-MS screen, i.e. 24 h after CFA-injection (controls: Veh-injection) and 28 days after SNI (controls: Sham surgery). In total accordance with our MS data, we recapitulated the differential and pain model-specific regulation of Serca not only in DRG (Fig. 3A, 3B, 3D; please see supplemental Fig. S7 for full blots), but also in distinct fascicles of the sciatic nerve (Fig. 3C, 3D; supplemental Fig. S8).
We next sought to investigate whether these changes in Serca abundance are indicative of altered Serca function. To this end, we cultured DRG neurons from CFA-and Veh-injected animals and performed ratiometric calcium imaging (Fig. 3E,  3F). Serca is a calcium-ATPase in the endoplasmic reticulum (ER) membrane and as such controls cytoplasmic calcium levels by active calcium transport from the cytosol to the ER. Hence, monitoring cytoplasmic calcium levels upon application of the Serca inhibitor thapsigargin can indirectly assess Serca activity. In order to do so without interference by the plasma membrane calcium ATPase (PMCA), we used an established protocol consisting of initial block of PMCA followed by coapplication of thapsigargin at a concentration shown to be specific for Serca (please see Experimental Procedures for details) (49,50). Interestingly, application of both blockers revealed increased cytosolic calcium levels in DRG neurons from CFAinjected animals suggesting higher Serca activity compared with Veh-injected controls (Fig. 3E, 3F). In contrast, intracellular calcium levels were comparable when only PMCA was inhibited indicating no difference of PMCA activity among conditions (Fig.  3E, 3F). Importantly, and in line with our functional data, a retrospective in silico interrogation of our proteomics results found comparable PMCA levels in CFA or SNI animals and their respective controls (q-values for PMCA: CFA: 0.357; SNI: 0.645; supplemental Table S3). These data highlight the capability of DIA-MS profiling to predict functional alterations of distinct proteins within the same cellular pathway.
Major Cellular Signaling Networks are Altered in Both Chronic Pain Models-The high statistical power and reproducibility of our DIA-MS screen provided us with the opportunity to meaningfully interrogate the data set in respect to changes in protein networks. To this end, we tested for evidence of associations among all significantly regulated proteins in both pain models (q-value Ͻ 0.05, supplemental Table  S3) using the web-based application STRING (45). Indeed, we found associations among 40 out of 64 proteins regulated during inflammatory pain and 56 out of 77 regulated proteins during neuropathic pain, respectively (Fig. 4A, 4B). Despite the certain degree of ambiguity intrinsic to functional categorization based on previously annotated information, this anal-ysis underlined mutual interconnections among regulated proteins and associated functional processes (Fig. 4). For example, in the neuropathic pain model we observed the down-regulation of structural proteins involved in cytoskeleton organization and myelin maintenance (supplemental Table  S3 and Fig. 4B; pink-colored protein nodes). Alterations in these proteins have been reported in another neuropathic pain model and are consistent with nerve injury and repair processes (11).
Moreover, our network analysis indicated differential alterations of multiple proteins implicated in mitochondrial function (Fig. 4A, 4B: represented by orange-brown-colored protein nodes) and protein biosynthesis/maturation (Fig. 4A, 4B: represented by green-colored protein nodes). Both cellular pathways are crucially interconnected determinants of apoptosis, energy metabolism as well as oxidative and ER stress (76 -79). Furthermore, the notion is evolving that mitochondrial abnormality (33, 80 -82) and ER stress (77,79,83) contribute to several forms of neuropathic pain. Mitochondrial dysfunction in painful conditions appears to be multifactorial as it may arise from changes in the activity of the electron transport chain (ETC) (33, 80 -82), a concomitant increase of reactive oxygen species leading to oxidative stress (84) or mitochondrial calcium handling (85). Similarly, mitochondrial fission represents a crucial step for cellular homeostasis and has recently been implicated in painful neuropathies brought on by cancer chemotherapy and anti-HIV medication (41). ER stress represents another significant driver of pain in diverse neuropathic models such as diabetic neuropathy (77,79,83). ER stress may be caused by a disequilibrium in ER homeostasis regarding protein demand, synthesis and folding. Its pronounced impact on pain phenotypes has recently been demonstrated in healthy animals where initiation of ER stress signaling cascade generated lasting pain, which could be reversed by ER stress blockers (77).
In addition to their impact on neuropathic pain, we observed a previously unreported regulation of proteins involved in mitochondrial function (orange-brown-colored nodes) and protein biosynthesis/maturation (green-colored nodes) upon CFA-induced inflammatory pain (Fig. 4A). Prompted by this finding, we set out to test whether manipulation of selected candidates in vivo may affect mouse pain behaviors upon CFA-injection. For this purpose, we looked for agents neither affecting general health nor motor function in naïve animals and at the same time being suitable to generally manipulate  Table S3), which exhibit reported associations among each other are displayed in the depicted networks (number of observed interactions versus expected interactions in a random protein set: inflammatory pain: 111 versus 22; neuropathic pain: 162 versus 27; assessed by STRING). Association networks were generated by STRING and depict the "confidence view," i.e. thicker lines represent stronger interactions (for details please see Experimental Procedures). Application of the MCL algorithm (STRING) highlights protein clusters by connecting protein nodes with continuous lines. Colors of protein nodes denote protein clusters involved in mitochondrial function, cytoskeleton and myelin maintenance, protein biosynthesis and maturation and others, respectively. Intercluster edges are displayed by dashed-lines. mitochondria and protein maturation in the ER in vivo (31,32). For the first, we used the mitochondrial complex I inhibitor Rotenone (32,33) as many of the regulated proteins are members of the respiratory electron chain complex I (e.g. Ndufs4, Ndufv2; Fig 4A). For the latter, we used PACMA 31, a nonselective inhibitor of protein disulfide isomerases (PDI) (31), of which three isoforms (PDIA1, PDIA3, and PDIA6) were strongly up-regulated during inflammatory pain ( Fig. 4A and supplemental Table S3). Indeed unilateral injection of Rotenone into the ipsilateral paw of CFA-treated mice resulted in diminished hypersensitivity to innocuous mechanical stimulation delivered with a plantar aesthesiometer compared with Veh-injected controls (Fig. 5A). This finding considerably extends previous studies that reported the preventive action of modulating mitochondrial function by application of inhibitors before induction and chronification of inflammatory pain (32,41). Likewise, we observed significantly less CFA-mediated mechanical allodynia upon intraperitoneal injection of PACMA 31 compared with Veh-injected mice (Fig. 5B).
Even though we used similar concentrations of Rotenone and PACMA 31 as described previously (31)(32)(33), we cannot unambiguously determine whether the observed analgesic action of Rotenone and PACMA 31 is solely mediated by inhibition of complex I and PDI in somatosensory neurons, respectively. Sensory neuron-specific conditional knockout mice would be required to critically interrogate potential offtarget effects. These may derive from nonspecific action either on other proteins in sensory neurons and/or effects on nonneuronal cells in vivo (e.g. different cell types in the skin or DRG). The latter hypothesis is especially relevant to systemic injection of PACMA 31, yet neither general motor coordination nor withdrawal latencies of contralateral (noninjected) paws were affected (Fig. 5) arguing against major side-effects elicited by compound injection. These observations are in line with other studies using these compounds (31)(32)(33).
In summary, our results suggest an acute analgesic potential of mitochondrial complex I and PDI inhibitors in fully established inflammatory pain and provide a proof-of-principle that pathologically relevant proteins and protein networks can be uncovered with our profiling approach.

DISCUSSION
The unbiased characterization of protein signaling networks is key for understanding the molecular signature of health and disease. Here we present an integrated workflow combining mouse behavior and DIA-MS-based protein profiling to standardize the reproducible detection and quantification of several thousand proteins in two mouse models of chronic pain. We found substantial and pain model-specific alterations of several dozens of proteins, among which 50% have not been linked to nociception and pain before. Strikingly, many of those represent novel and functionally uncharacterized candidates. Comparison with known biology paired with our validation experiments strongly suggests that indeed proteins relevant for pain and somatosensory signaling were identified. Taken together, we present (1) the so far largest compendium of the DRG proteome and (2) a list of proteins specifically regulated during chronic pain states. These constitute valuable resources enabling further investigations of pathologies of the somatosensory system.
We focused our analysis on DRG isolated from male mice subjected to two different chronic pain states. As DRG harbor primary sensory neurons, which are the interface of the nervous system with external and internal environments, they are crucially implicated in the beginning of somatosensory pathways including nociception (3,4,86). This, together with their accessibility, renders DRG key cellular targets for analgesic therapies (2,4,87).
Because of the cellular complexity of DRG, we cannot assign the detected changes in protein abundance to specific cell types, at least at the level of mass spectrometry. However, we partially solved this limitation employing MS-independent and cell type-specific techniques for validation of our results. For example, we show that the regulation of Serca occurs in sensory neurons as revealed by calcium-imaging in primary DRG cultures as well as by immunohistochemistry on DRG and sciatic nerves.
The current profiling study recapitulates known biology regarding the identification of proteins previously linked to pain. FIG. 5. Mitochondrial and PDI activity play a role for chronic inflammatory pain in mice. A, Rotenone-mediated analgesia of established chronic inflammatory pain. Illustration of paw withdrawal latencies of CFA-injected animals, which received ipsilaterally (ipsi) either Veh or Rotenone (Rot) 24 h after CFA-application (***, p Ͻ 0.001 and *, p Ͻ 0.05 by one-way ANOVA followed by BonferroniЈs multiple comparison test; 17 animals per condition from 3 independent mouse cohorts). B, Inhibition of PDI by PACMA 31 attenuated established chronic inflammatory pain. Illustration of paw withdrawal latencies of CFA-injected animals, which received either Veh or PACMA 31 (PACMA) via intraperitoneal injection 24 h after CFA-application (***, p Ͻ 0.001 and *, p Ͻ 0.05 by one-way ANOVA followed by Bonfer-roniЈs multiple comparison test; 10 animals per condition from 2 independent mouse cohorts). Importantly, withdrawal responses of contralateral (i.e. untreated) paws were unaffected arguing against unspecific motor effects upon intraperitoneal injection of PACMA 31.
In addition, we provide a proof-of-principle with orthogonal validation experiments that our DIA-MS workflow is capable to uncover functionally meaningful changes. Using various methods, we correlated our DIA-MS data with altered Serca function and expression in DRG neurons. Serca is crucially involved in neuronal calcium and ER homeostasis by governing ER calcium uptake and sequestration. Consequently, the activity of Serca requires precise control to enable adequate neuronal responses to different conditions. In fact, Sercaactivity in DRG neurons is modulated by noxious substances (88), protein kinase C (89) and upon painful neuropathy (90). In turn, injection of Serca-inhibitors has recently been shown to elicit ER stress-mediated nociceptive responses in vivo (77). Our profiling data suggest diminished abundance of Serca in neuropathic mice that is in complete agreement with previous reports (77,90). Strikingly, we observed the opposite regulation, i.e. augmented Serca levels and function specifically during inflammatory pain. Although the underlying mechanisms require further investigations, which are beyond the scope of this study, these results propose the existence of signaling mechanisms converging on Serca in a pain modelspecific manner.
Integration of DIA-MS profiling together with protein network analysis offers another level of understanding the molecular signature of diseases. Using the web-based application STRING, we visualized associations among significantly regulated proteins (Fig. 4), which revealed that components of major signaling pathways are altered during chronic pain, e.g. mitochondrial function as well as protein synthesis and maturation. Although this notion is not new and many transcriptional profiling techniques provide a systems-wide view on cellular networks (7,8,29,91), proteome profiling has so far made only minor contributions (92)(93)(94). Undersampling and limited reproducibility inherent to commonly used shotgun proteomics approaches (16 -19) represent major obstacles for proteome-based systems biology. This is in stark contrast to the emerging importance of "network medicine" as a powerful strategy to correct functional misalignments of complex cellular processes (2,76,95). Consequently, it might provide promising tools for achieving analgesia as exemplified by studies linking mitochondrial dysfunction (33, 80 -82) and ER stress to neuropathic pain (77,79,83). However, the underlying molecular changes and altered proteins are only scarcely understood preventing the design of targeted interventions with therapeutic potential. Here, we report highly reproducible DIA-proteome profiling of 2526 proteins, which granted unique insights into significant alterations of candidate proteins belonging to distinct cellular signaling networks. These data are not only in agreement with previous studies implicating mitochondrial (33, 80 -82) and ER (77,79,83) abnormalities in neuropathic pain, but importantly (i) identify affected members of these protein networks and (ii) provide evidence for their differential involvement also in inflammatory pain. Indeed, we could demonstrate attenuation of fully es-tablished chronic inflammatory pain by application of mitochondrial complex I or protein disulfide isomerase inhibitors in vivo. Taken together, we expect the current study to spark further research to decipher the exact mechanisms how the here detected protein changes contribute to pain pathology.
Although we affirm selected candidates by examples of orthogonal validation as well as correlation with known pain biology, we did not detect several previously reported regulations of proteins during chronic pain. For example, we did not observe a regulation of TRPA1 in our study whereas other reports demonstrated its up-regulation in DRG of CFA-injected mice (96). Several technical considerations might underlie these discrepancies. On one hand, we pooled ipsilateral lumbar DRG from different levels in order to obtain sufficient tissue per DIA-MS replicate. Although this is common practice in most -omics studies, it is noteworthy that even neurons within one DRG are highly heterogeneous, as affected and nonaffected neurons are intermingled upon induction of chronic pain in mice (97). Furthermore, the relative proportions of affected and nonaffected neurons in one DRG critically depend on the injection site of agents in respect to the projection sites of individual nerve branches innervating the hind paw (97). Indeed, differences between affected versus nonaffected nerve branches might potentially explain our results showing disparate regulation of Serca in distinct sciatic nerve fascicles. On the other hand, our study differs from previous proteomics studies in respect to biochemical sample preparation. With the goal of deciphering the (patho-) physiologically relevant (4 -6, 58 -61, 63) but traditionally understudied membrane proteome (64) we used membrane-enriched lysates of DRG. Consequently, the depletion of soluble proteins in our samples might explain the marginal overlap of our results with published proteomics studies using wholecell lysates (9 -14).
In addition, species, tissue and pain model differences have to be considered when interpreting and comparing data among studies (97). For instance, Gemes and colleagues reported elevated activity of plasma membrane calcium-ATPase (PMCA) in the neuropathic pain model of spinal nerve ligation in rats (98). In contrast, we employed the CFA-induced inflammatory pain model and the neuropathic pain model of spared nerve injury in mice and could not detect any alterations of PMCA, neither by DIA-MS profiling nor on the functional level using calcium-imaging in sensory neurons harvested from CFA-injected mice.
The here presented spectral library is enriched with fingerprints for transmembrane and membrane-associated proteins found in diverse cellular compartments like the plasma membrane, membranes of the endoplasmic reticulum as well as mitochondrial membranes. For example, TRPA1 and TRPV1 channels are molecular noxious stimulus detectors and critically implicated in nociceptive signaling and pain (4, 99 -101). Therefore, they are considered as potential targets for pharmacological interventions (98,(102)(103)(104). Moreover, Piezo2 channels confer mechanosensitivity to sensory neurons (105,106) as well as Merkel cells of the skin (107,108), and articular cartilage cells (109) highlighting the applicability of our spectral library to different tissues beyond DRG. However, despite successful enrichment of membrane proteins, our spectral library (and therefore our DIA-MS results) is not complete. Although we provide fingerprints for the sodium channels Nav1.8 and Nav1.9, the spectral library lacks information on the sodium channel Nav1.7 known to be involved in several pain states (5). The reasons can be manifold, ranging from relative low expression levels to insufficient solubilization or localization in detergent-resistant microdomains -all factors known to render membrane protein analysis challenging (64). Nonetheless, and compared with published proteomics studies on pain (9 -14) or DRG membranes (65), our spectral library constitutes a unique resource for many hundreds of DRG membrane proteins, which has the potential to enormously facilitate research on these proteins in any tissue. Moreover, the nature of DIA-MS profiling allows any experimenter to continuously extend this library, and use it for in-silico analysis of generated data sets without requiring new experiments. Therefore, the depth of protein profiling will steadily improve once experiments in DIA-mode are commonly used by the research community.
Another powerful indication of the value of our data set is the high rate of reproducibility we achieved in our DIA-MS experiments. We reproducibly profiled the same set of 2,526 proteins in all 12 DIA-MS experiments performed, i.e. in three replicates of each of the four conditions. Taking into account that each replicate consisted of independent rounds of surgeries/injections, tissue extraction and biochemical preparation, all being subject to high variability, the observed reproducibility is remarkable and qualitatively superior to shotgun proteomics studies (16 -19).
Thus, the (1) enrichment of membrane proteins in our spectral library paired with the (2) high rate of reproducibility, and (3) the possibility for future in-silico examination of generated data sets implement a standardized toolbox for profiling these proteins in any tissue of interest.
Collectively, our findings support the notion that the here generated DRG protein compendium represents a rich and valuable resource for the research community. It provides a framework for comprehensive characterization of particular molecular signatures underlying diverse pathologies, well beyond the here described application in mouse models of chronic pain.