Accounting for Protein Subcellular Localization: A Compartmental Map of the Rat Liver Proteome*

Accurate knowledge of the intracellular location of proteins is important for numerous areas of biomedical research including assessing fidelity of putative protein-protein interactions, modeling cellular processes at a system-wide level and investigating metabolic and disease pathways. Many proteins have not been localized, or have been incompletely localized, partly because most studies do not account for entire subcellular distribution. Thus, proteins are frequently assigned to one organelle whereas a significant fraction may reside elsewhere. As a step toward a comprehensive cellular map, we used subcellular fractionation with classic balance sheet analysis and isobaric labeling/quantitative mass spectrometry to assign locations to >6000 rat liver proteins. We provide quantitative data and error estimates describing the distribution of each protein among the eight major cellular compartments: nucleus, mitochondria, lysosomes, peroxisomes, endoplasmic reticulum, Golgi, plasma membrane and cytosol. Accounting for total intracellular distribution improves quality of organelle assignments and assigns proteins with multiple locations. Protein assignments and supporting data are available online through the Prolocate website (http://prolocate.cabm.rutgers.edu). As an example of the utility of this data set, we have used organelle assignments to help analyze whole exome sequencing data from an infant dying at 6 months of age from a suspected neurodegenerative lysosomal storage disorder of unknown etiology. Sequencing data was prioritized using lists of lysosomal proteins comprising well-established residents of this organelle as well as novel candidates identified in this study. The latter included copper transporter 1, encoded by SLC31A1, which we localized to both the plasma membrane and lysosome. The patient harbors two predicted loss of function mutations in SLC31A1, suggesting that this may represent a heretofore undescribed recessive lysosomal storage disease gene.

Accurate knowledge of the intracellular location of proteins is important for numerous areas of biomedical research including assessing fidelity of putative proteinprotein interactions, modeling cellular processes at a system-wide level and investigating metabolic and disease pathways. Many proteins have not been localized, or have been incompletely localized, partly because most studies do not account for entire subcellular distribution. Thus, proteins are frequently assigned to one organelle whereas a significant fraction may reside elsewhere. As a step toward a comprehensive cellular map, we used subcellular fractionation with classic balance sheet analysis and isobaric labeling/quantitative mass spectrometry to assign locations to >6000 rat liver proteins. We provide quantitative data and error estimates describing the distribution of each protein among the eight major cellular compartments: nucleus, mitochondria, lysosomes, peroxisomes, endoplasmic reticulum, Golgi, plasma membrane and cytosol. Accounting for total intracellular distribution improves quality of organelle assignments and assigns proteins with multiple locations. Protein assignments and supporting data are available online through the Prolocate website (http://prolocate.cabm.rutgers. edu). As an example of the utility of this data set, we have used organelle assignments to help analyze whole exome sequencing data from an infant dying at 6 months of age from a suspected neurodegenerative lysosomal storage disorder of unknown etiology. Sequencing data was prioritized using lists of lysosomal proteins comprising wellestablished residents of this organelle as well as novel candidates identified in this study. The latter included copper transporter 1, encoded by SLC31A1, which we localized to both the plasma membrane and lysosome. The patient harbors two predicted loss of function mutations in SLC31A1, suggesting that this may represent a heretofore undescribed recessive lysosomal storage disease gene. An outstanding challenge in cell biology is to understand how proteins are organized and interact in biological pathways. Analogous to the organs found in multicellular organisms, the cell contains organelles, which are macromolecular assemblies whose components orchestrate specialized functions. Knowledge of the organellar residence of a given protein can provide valuable clues to its physiological role whereas conversely, understanding the protein composition of a given organelle helps describe its functional capabilities. In addition, accurate knowledge of the location of proteins is critical in assessing the biological significance of a variety of big data initiatives including protein network and pathway analysis (1,2) as well as facilitating studies on the genetic basis of disease (3).
Microscopy and subcellular fractionation are typically used to ascertain protein localization. Fluorescence and electronmicroscopy-based approaches (4) can localize proteins to morphologically recognizable cellular structures and to each other, potentially with extremely high spatial and temporal resolution but experimental parameters may result in incorrect assignments. For example, when expressing exogenous proteins, the presence of a tag and/or nonphysiological steady-state levels can disrupt their normal trafficking and localization. For immunolocalization of endogenous proteins, there is increasing awareness that antibody specificity is a serious area of concern (5). In addition, even with high quality reagents, it may be difficult to determine whether a signal represents all or a subset of the protein of interest because of a variety of factors including destruction of a tag after targeting of a fusion protein to a particular organelle, quenching of fluorescent labels under certain cellular environments or relative accessibility of antibodies to different cellular structures. It is also worth noting that unless conducted with rigorous quantitation and appropriate controls, there is a tendency with morphological approaches to focus on fluorescence that is high intensity and punctate, versus lower intensity but diffuse.
Subcellular fractionation is a biochemical approach where cells are disrupted and cellular structures separated, typically by centrifugation based on their sedimentation coefficients (6). Proteins in fractions are measured using immunoassays, activity assays, or MS, the latter being amenable to highthroughput analysis. Analytical fractionation can be extremely quantitative, but there are limitations. While it is generally assumed that the protein composition of the vesicles represents the parent structures, some degree of mixing or relocalization during the fractionation process is possible. In addition, different types of structures may comigrate in a given fractionation scheme, potentially limiting resolution. Nonetheless, if one remains aware of these caveats when evaluating data, fractionation is a powerful tool for understanding the organization of the cell (7)(8)(9).
Many cellular localization studies employing subcellular fractionation have focused on cataloging proteins associated with individual organelles and therefore cannot employ balance sheet analysis ("bookkeeping") to account for the distribution of each protein throughout the cell. A failure to account for total distribution is a significant problem because, even if the bulk of a particular protein is distributed elsewhere within the cell, identification of some portion of this protein within a given organelle frequently can become dogma as its sole site of residence. Not only is this incorrect per se, this can lead to erroneous conclusions regarding function of a protein if it is considered only in the context of its organelle of minor residence. In contrast, bookkeeping is a central tenet of analytical subcellular fractionation and is used to tally recoveries, providing confidence that the distribution reflects the total cellular population of each protein.
An early study by Mann and coworkers used "protein correlation profiling" to follow distribution of multiple proteins in mouse liver subcellular fractions (10), and, more recently, high-quality maps have been published for mouse pluripotent stem cells (11) and HeLa cells (12). While these studies do list multiple residencies for some proteins, with the exception of estimating relative distribution among nuclear, cytosol, and organellar fractions (12), relative distributions among different organelles were not estimated.
We have combined multiple orthogonal subcellular fractionation schemes with high-resolution mass-spectrometry as a platform for an approach that applies the principles of bookkeeping to estimate the distribution of rat liver proteins associated with the eight major compartments that together account for most of the cell: nucleus, mitochondria, lysosomes, peroxisomes, endoplasmic reticulum (ER) 1 , Golgi, plasma membrane (PM) and cytosol. This has allowed us to compile a draft map of the mammalian liver proteome, providing information regarding cellular locations for approximately a third of the genes predicted to be expressed in liver based on transcriptome analysis (13). This database (Prolocate) is accessible through a web portal (prolocate.cabm.rutgers.edu) with the ability to query the subcellular distribution of individual proteins and examine underlying data, compile lists of candidate residents of different organelles, and identify proteins with similar multicompartmental distributions.

EXPERIMENTAL PROCEDURES
Subcellular Fractionation-All experiments and procedures involving live animals (adult male Wistar rats) were conducted in compliance with approved Institutional Animal Care and Use Committee protocols. In Experiments (Expts) A and C, animals were fasted overnight in solid bottom cages to reduce glycogen stores prior to fractionation, whereas animals used in Expts B and D were fed ad libitum.
Step-by-step subcellular fractionation protocols are provided in Supplemental Methods and the overall scheme is shown in Fig. 1A. To summarize, liver homogenization and classical differential centrifugation were performed essentially as described to produce N, M, L, P, and S fractions (14). Note that in this protocol, pellets are washed by a resuspension-recentrifugation procedure, with the combined supernatants being used for subsequent centrifugation steps. The L fraction was subfractionated to yield an L1 pellet and L2 supernatant as shown in Fig. 1A and as described in supplemental Methods. Further fractionation of the L1 fraction in a Nycodenz step gradient was conducted essentially as described (15) to prepare a fraction highly enriched in lysosomes over peroxisomes ("Nyc2") ( Fig. 1A and supplemental Methods). In Expts C and D, L1 fractions prepared from control and Triton WR 1339-("Triton") injected rats (16) were adjusted to a density of 1.18 g/ml, centrifuged, and the tube divided into top and bottom fractions ( Fig. 1A and supplemental Methods). Marker enzyme activity and protein assays were conducted as described previously (17).
Experimental Design and Statistical Rationale-This study consists of four experiments using rat liver, with livers from two animals used for each subcellular fractionation. Since feeding status does not appreciably affect the overall distributions of the vast majority of liver proteins (see Results), there were two biological replicates for each type of experiment. Expts A and B were biological replicates that included analysis of all differential centrifugation fractions. In addition, Expt A contained a single fraction from a Nycodenz gradient performed on the L1 differential centrifugation fraction to help distinguish lysosomal and peroxisomal proteins. Expts C and D represent biological replicates of isopycnic centrifugation experiments on lysosomal/peroxisome-enriched fractions from control and Triton-treated an-imals. Given the overall concordance from each set of experiments involving biological replicates, sample sizes were considered acceptable.
Mass Spectrometry-Samples indicated in Table I were reduced, alkylated and processed for MS by in-gel trypsin digestion as described (18). Note that the homogenate (H), used for bookkeeping purposes, was reconstituted by mixing appropriate ratios of the initial low speed pellet (N fraction) and the pooled post nuclear supernatant (E) shortly before sample processing to minimize agglutination and sampling error. Two purified histidine-tagged bacterial proteins (Northeast Structural Genomics Consortium target identification numbers DrR57 and GmR40, kindly provided by Dr. Guy Montelione) were added before digestion as internal standards as described (17). GmR40 was added at a constant protein weight/weight ratio of 1 part per 300 parts experimental sample. DrR57 was added at different amounts to yield weight/weight ratios of 0:1, 1:600, 1:300, or 1:100 of DrR57:experimental sample. Tryptic peptides were extracted and 100 g of each sample labeled with isobaric reagents (Expts A, B, and D, iTRAQ-8plex, Expt C, TMT-6plex) using manufacturer's protocols. All individual labeled samples were combined into one and desalted. To reduce sample complexity, the resulting pooled sample was further fractionnated by strong cation ion exchange (SCX) and/or alkaline reverse phase ("high pH," or HP) chromatography prior to LC-MS. Fractions resulting from off-line chromatography fractionation were then each analyzed by LC-MS using a Thermo Orbitrap Velos mass spectrometer. The number of individual LC-MS runs conducted for each experiment were as follows: Expt A, 288 (10 different SCX fractions and 254 different 2D (SCX x HP) fractions); Expt B, 252 (unfractionated, 16 different SCX fractions and 189 different 2D fractions); Expt C, 272 (unfractionated, 26 different SCX fractions, 40 different HP fractions, and 163 different 2D fractions); and Expt D, 199 (14 different SCX fractions and 153 different 2D fractions). (Note that some fractions were run in replicate, thus the number of fractions analyzed is less than the total number of LC-MS/MS runs). Chromatography details for prefractionation and LC-MS have been deposited in MassIVE in the methods folder for each experiment, as have LC-MS raw files, peak lists, search results, and processed data. GPM .xml search result files containing supporting information for spectral assignments can be viewed using the public installation of the GPM (http://h.thegpm.org/tandem/thegpm_upview.html) as well as using local installations (ftp://ftp.thegpm.org/projects/gpm/gpmxe-installer/). MS peak list data files (.mgf) were created using Proteome Discoverer version 1.4 (Thermo Fisher Scientific). Individual reporter ion intensities were abstracted from mgf files by summing intensities of ions within Ϯ20 ppm of each reporter ion theoretical m/z. The mgf files were further preprocessed to include only spectra where at least two of the reporter ions had intensities Ն1000. Spectra were matched to peptide sequences using a local implementation of the Global Proteome Machine (GPM) (14,15) with X!Tandem version SLEDGE-HAMMER (2013.09.01). All mgf files within a given experiment were searched together to produce a merged "MudPIT" output file. Target proteins consisted of 24,989 unique rat sequences (ENSEMBL Rnor_5.0.75 processed to retain only the first entry for any sequence that was assigned to two or more protein accession numbers), the two bacterial internal standards and 41 potential contaminants (dust/ contact proteins and trypsin from the GPM cRAP database). All search parameters are listed in the GPM xml files. Briefly, searches allowed for up to one missed trypsin cleavage and used precursor and product mass tolerances of 10 and 20 ppm, respectively. Initial searching was done using complete carboxamidomethylation at cysteine and selenocysteine, complete isobaric labeling at amino-terminal peptide residues and lysine, and potential methionine oxidation. During model refinement, the following potential modifications were allowed: round 1, methionine and tryptophan oxidation, asparagine and glutamine deamidation; round 2, methionine and tryptophan dioxidation; round 3, isobaric labeling at tyrosine and loss of label at the amino terminus and lysine. The spectrum synthesis option was set to "no." Peptide false positive rates (FPR) were calculated by GPM and are as follows: Expt A, 0.58%; Expt B, 0.63%; Expt C, 0.56% and Expt D, 0.85%. Peptide assignments for all spectra with peptide expectation values of Յ 0.1 were exported to excel files and merged with processed reporter ion intensities (supplemental Workbooks S2A-S2D, GPM… worksheets). For the latter, fragment ion m/z were first linearly recalibrated so that the most intense reporter ion had zero mass error. Intensities associated with each reporter ion (m/z Ϯ10 ppm) were adjusted for spillover into neighboring channels using vendor-supplied correction factors.
After mapping spectra to peptides, we conducted a series of steps to improve data quality. All steps are documented in spreadsheets for each experiment (supplemental Workbooks S2A-S2D), downloadable through the Prolocate website (prolocate.cabm.rutgers.edu) and also deposited in MassIVE.
1. The first step was to filter spectra based on peptide assignment. Peptides were not used for quantification unless they met all the following criteria: fully tryptic with no missed cleavages; complete isobaric labeling of amino-terminal residues and lysines; and absence of adventitious modifications deemed to increase variability (asparagine or glutamine deamidation, methionine di-oxidation, tryptophan mono-and di-oxidation, and isobaric labeling of tyrosine at positions other than the amino terminus). Spectra meeting these criteria are flagged "Y" in the "acceptable peptide" columns in supplemental Workbooks S2A-S2D.
2. We used the internal standard GmR40 added in equal amounts to all samples to help correct for potential differences in efficiencies of digestion, labeling and recovery prior to pooling the samples in each experiment. Here, the intensity of each reporter ion from a given spectrum was divided by the corresponding GmR40 reporter ion intensity (average of all spectra from acceptable GmR40 peptides). Corrected normalized reporter ion ratios were then calculated.
3. We conducted a balance sheet analysis for each spectra. Here, we first scaled the normalized reporter ion ratio (based on labeling of same amount of protein in each sample) to that which would be present in the entire fraction derived from one gram of liver. We then estimated recoveries by summing fractions of interest and dividing by the starting material (i.e. Expts A and B, (NϩMϩL1ϩL2ϩPϩS)/H; Expts C and D, (top sucrose ϩ bottom sucrose)/L1 for both the control and Triton-treated samples). In theory, the amount of peptide in fractionated material should be the same as that of the starting material, and this ratio can be used to check data quality. Spectra with recoveries less than 2/3 or greater than 3/2 were eliminated from further analysis.
Gene Product Assignments-A single gene can be associated with multiple protein accession numbers because of alternative splicing, protein processing, and other considerations such as database annotation. Protein accession numbers are mapped to the gene name (obtained using the ENSEMBL Biomart tool) which is used as the primary identifier in this study. In some cases, peptides could be equally well assigned to multiple, genetically distinct proteins. Here, we assigned the peptide to the protein with the highest GPM expectation value which was generally the one with the greatest sequence coverage in the search. In addition to this "best" protein assignment, supplemental Workbooks S2A-S2D lists all alternative protein accession numbers, gene accession numbers, and gene names associated with a given peptide. Organelle assignments were generated for gene products represented by at least one rat protein with a GPM expectation score of Յ 10 Ϫ3 .
Multiple Assignments Using Constrained Proportionate Assignment (CPA)-We developed this procedure for analysis of Expts A and B. CPA uses enrichment factors (EF) (levels of a protein of interest in different centrifugation fractions compared with its level in starting material), expressed as profiles, to assign proteins to different subcellular compartments. To do this, we select several proteins that have well-established residencies for each of the compartments (single-compartment marker or reference proteins, supplemental Table  S2) and use these to produce an average profile for each compartment (eight compartments for Expt A, seven compartments for Expt B). We then fit the profile of each individual protein in our data set to an optimal combination of the reference profiles, thereby obtaining estimates of its distribution among all compartments.
The detailed CPA procedure conducted on Expt A is as follows: 1. We first calculate seven EFs (N, M, L1, L2, P, S and Nyc2) for each spectrum. EFs are equivalent to relative specific activities (RSAs) used in classical subcellular fractionation (14) and the steps for these calculations are detailed in supplemental Workbook S2A. EFs are log2 transformed after setting values Ͼ16 and Ͻ1/16 to 16 and 1/16, respectively. Spectra with unacceptable recoveries or assigned to unacceptable peptides are dropped from further consideration (see above).
2. We then perform an outlier screen for proteins with at least three spectra. For each EF for a given protein, we use all spectra to compute "normal scores," given by z ϭ (x -mean (x))/std (x), where x is the log2 transformed EF. We discard any spectrum for which any z lies in the outer 1% tails of the normal distribution. This process is repeated for all proteins.
3. Once outlier spectra have been removed, the next step is to determine a mean profile for each protein. If a protein has at least 4 spectra and at least 3 different sequences (peptides), we ordinarily have enough data to fit a random effects model. This model computes a weighted average and standard error of the log2 EFs from all spectra associated with a given protein, accounting for the fact that spectra are nested within sequences. This computation is carried out using the "lmer" function in the "lme4" R package (19). The result is essentially the mean and standard error of the observations, with an adjustment for the nested structure of the data. This procedure prevents a sequence with a very large number of spectra from dominating the estimates of the mean and standard error. Occasionally the "lmer" program will fail to converge for a particular EF. If that happens or if the condition of having at least 4 spectra and 3 different sequences is not met, we compute the simple mean and standard error of the log2 transformed data for that EF. If there are fewer than 3 spectra, we cannot compute a standard error, and thus only report the mean. When this process has been completed, every protein has a mean profile consisting of mean EFs for each of the seven fractions, denotedx ϭ (x 1 , x 2 , …, x 7 ).
4. Using the reference proteins listed in supplemental Table S2, we average all spectra containing acceptable log2 transformed EFs to form profiles for each of the eight compartments (Cytosol, ER, Golgi, Lysosome, Mitochondria, Nucleus, Peroxisome, and PM). The eight compartmental profiles are denoted ass 1 ,s 2 , …,s 8 , where each compartmental profiles j is a list of its seven mean EFs.
5. For each protein, we find estimated proportions p 1 , p 2 , …, p 8 so that weighted profileỹ ϭ p 1s 1 ϩ … ϩ p 8s 8 is as close as possible to the observed profilex, subject to the constraints 0 Յ p j Յ 1 for all j, and jϭ1 8 p j ϭ 1 where "close" is defined by the minimizing the sum of squares of the differences Q ϭ iϭ1 In other words, we may view the p j as proportional allocations of the eight standard profiles to formỹ, which is as close as possible to the observedx for this particular protein. This constrained optimization is carried out using the "spg" function in the R package "BB" (20) (also manuscript in preparation, Varadhan R "Projection of infeasible iterates in constrained nonlinear optimization using BB::spg"). These estimated proportions are referred to as classification coefficients and are used to estimate the relative amount of a given protein that is associated with each compartment.
6) For proteins with Ն3 acceptable spectra, we also compute 95% confidence limits for the p j by creating 1000 bootstrap samples of the original data, and fitting the CPA method to each of them. The 95% confidence limits were obtained by ordering the estimates, and selecting the 2.5% and 97.5% percentiles of the estimates.
A similar procedure is used for Expt B (supplemental Workbook S2B), except there is no Nyc2 fraction, and we estimate distribution among seven compartments, merging the lysosomal and peroxisomal compartments into a Lyso/Perox compartment.
Classification of Lysosomal and Nonlysosomal Proteins Using Gaussian Clustering-This procedure used for analysis of Expts C and D is a variant of classical Fischer discriminant analysis. The purpose of the Triton density shift experiments is to use the differential separation of lysosomal proteins from nonlysosomal proteins when animals are injected with Triton WR-1339, as described in Results. The ratio of a particular protein in the top (t) to bottom (b) sucrose fractions from the Triton-treated animals is denoted Tt/Tb, and the corresponding ratio for uninjected control animals is denoted Ct/Cb. For nonlysosomal proteins, Triton has little or no effect on these ratios, whereas for lysosomal proteins, Tt/Tb Ͼ Ct/Cb. Our classification procedure is as follows: 1. The ratios are log2 transformed after setting values of Ͼ16 and Ͻ1/16 to 16 and 1/16, respectively. Spectra with unacceptable recoveries or assigned to unacceptable peptides are dropped from further consideration (see above).
2. Next we remove outliers similarly to that described in Step 2 for the CPA procedure (see above) i.e. for each protein we remove any spectrum for which either Log2(Tt/Tb) or Log2(Ct/Cb) lies in the outer 1% tail of a normal distribution fitted to all the spectra of that protein.
3. Having obtained a clean set of spectra, we compute the mean values of Log2(Tt/Tb) and Log2(Ct/Cb) for each protein and also the covariance matrix for these quantities. We use reference protein means and covariances for Log2(Tt/Tb) and Log2(Ct/Cb) as a "training" set for our classification procedure, which is a generalization of Fisher quadratic discriminant analysis. Our procedure uses the training set to define a mixture of two bivariate normal distributions with the aforementioned means and covariances, one for the lysosomal reference proteins and one for the nonlysosomal ("other") reference proteins. We initially evaluated the proteins used as markers in Expts A and B as well as those whose assignments were consistent in Expts A and B with various thresholds for classification coefficient point estimates and lower confidence limits. All showed a similar pattern, but given the complete concordance of proteins with a lower confidence limit of Ն 0.7 in both Expts A and B (Table II), we used these for purposes of classification, excluding GZMA, EPX, and MPO from the lysosomal classification set as these are found in specialized lysosome-like granules but do not undergo the Triton shift. We use the R procedure "Mclust" from the "Mclust" package (21,22) to define these two normal distributions and to predict, for each protein, the probability that it belongs to the lysosomal cluster or the other cluster. The predictions consist of probability assignments to these groups.
We generate confidence limits for each protein using a parametric bootstrap procedure. Here, for each protein, using the covariance matrix computed as described above, we generate a bivariate normal parametric bootstrap sample centered at the mean [Log2(Tt/Tb), Log2(Ct/Cb)]. We use these estimates to generate probabilities of the two subcellular locations (lysosomal or other) for all proteins using the procedure described in step 3 above. We repeat this bootstrap pro-cedure 1000 times to obtain 1000 predictions of subcellular locations for each protein. We obtain 95% confidence intervals by selecting the 25th and 750th values of the sorted bootstrap estimates.
Other Classification Analyses-We performed support vector machine clustering using the "svm" function in the "e1071" R package with the marker proteins (supplemental Table S2) as a training set. We performed hierarchical clustering using the "hclust" function in the R package. Clustering was based on Euclidean distances between points using the complete linkage method.
Distances Between Proteins-The distance calculator on the Prolocate website can be used to find proteins with similar distributions, regardless of organelle assignments. Consider two proteins i and j from Expt A with profiles s i and s j (each of which is a list of the seven log2-transformed EFs). The distance between them is the Euclidean distance d ij is defined analogously, with the profiles containing six EFs (because there is no Nyc2 column in this experiment). For Expt C, the "profiles" are simply Log2(Ct/Cb) and Log2(Tt/Tb), so that d ij To utilize data from Expts A and B (for proteins contained in both), we simply sum the distances, d ij Distances for other combinations are similarly obtained by summing the appropriate individual distances.
Mapping ENSEMBL Rat Proteins to Other Databases and/or to Mouse and Human Orthologs-Mapping results are in supplemental Workbook S3. A list of 9,871 ENSEMBL protein sequences identified from MS searches with associated accession numbers and gene names were rooted in the ENSEMBL Rnor_5.0.75 database. Corresponding rat UniProt entries were obtained from the top BLASTP hit against the UniProt database (Major release-2015_10). Matches were considered acceptable if they covered Ն75% of the ENSEMBL rat protein sequence and had an expectation score Յ0.001. Mouse and human orthologs of rat proteins were identified via a BLASTP Reciprocal Best Hits (RBH) approach where each rat protein sequence was BLASTed against the Mouse and Human proteomes (UniProt, major release 2015_10, ENSEMBL GRCm38.p4 (Mus musculus) and GRCh38.p5 (Homo sapiens)). The top three BLAST hits in each ortholog proteome were BLASTed back against the 9871 rat proteins identified by MS. The data set in supplemental Workbook S3, worksheet "Blasts with statistics" utilizes a y4x nomenclature to denote the results of these analyses where x denotes the BLAST rank of an ortholog candidate protein whereas y denotes the BLAST rank of the original rat query protein when the ortholog candidate protein is BLASTed back against the rat proteome. For example, the nomenclature 141 is used for a RBH identified ortholog where the first BLAST identified protein in the ortholog proteome (x ϭ 1), when BLASTed back against the rat proteome, provides a BLAST result where the original query rat protein is the first identified protein (y ϭ 1). BLAST hits were only considered where e-values of pair-wise sequence alignments were Յ0.1 and the pair-wise sequence alignments covered Ն75% of the length of the query rat protein. Although true RBH identified orthologs (141) where identified for the majority of rat proteins, there were instances where the analyses returned proteins not considered true RBH orthologs. For example, when the first ortholog candidate was BLASTed back against the rat proteome, the original rat query protein was the second returned BLAST hit (241). Analysis of such results suggests that they primarily arise from multiple isoforms in the ortholog proteome.
supplemental Workbook S3 also contains accession numbers and assignments from a variety of other studies listing organelle assignments (worksheets "CDbBS" and "AllOtherSudies"). These are mapped to the identifiers in our study in worksheet "Blasts_with_ lookups." This is reproduced in a more streamlined form in supplemental Workbook S1 worksheet "OtherStudies" and assignments compared with our study in supplemental Workbook S1 worksheet "ProteinDistributionExptsA-D." For the mouse liver global subcellular fractionation study of Foster et al. (10), only proteins denoted as having a single "Yes" when considering all compartments listed in their supplemental Table S5 (refined location worksheet) were used for the analysis, allowing "Probably," "Co-migrating," and "No" in other fields (supplemental Workbook S3, worksheet AllOtherStudies).
Curation of a Benchmark Set-For the Compartments database human benchmark set of 12,892 proteins (23), only 7955 were assigned to a single compartment ("ϩ" for one compartment and either "-" or not listed for all other compartments) used in our classification scheme (supplemental Workbook S3, Worksheet CDbBS). Of these, 2256 orthologs were present in our entire data set, and 809 were assigned to our high stringency data set. Creation of the "curated" CDbBS was conducted in two rounds. We initially reviewed 197 proteins, 34 of which the primary compartment identified in our study and in the CDbBS agreed. Here, two investigators each with years of experience in evaluating subcellular location (MB and DES), blinded to the Prolocate and CDbBS assignments and to whether these agreed or disagreed, independently evaluated evidence for subcellular localization based on the validity of existing experimental evidence in published papers, evidence in the Human Protein Atlas (HPA) (24), and, occasionally, common knowledge (e.g. ribosomal proteins are present both in the cytosol and on ER). A third investigator (PL) then reviewed/arbitrated these assignments in a nonblinded manner. There was strong evidence supporting 31 of the 34 proteins that agreed, with inconclusive (but not contradictory) information for the remaining three. Given that the common assignments had excellent agreement we accepted the remaining common assignments. The second round of analysis was conducted on the remaining 118 proteins where the primary compartment identified in our study and in the CDbBS disagreed. Here, the list was split between the first two investigators for initial blind evaluation, and the third investigator reviewed and arbitrated these assignments. Individual investigator assignments, rationale and PubMed references supporting these assignments are documented in supplemental Workbook S1, Worksheet OtherStudies).
Exome Sequencing Analysis-Whole exome sequencing at ϳ100ϫ coverage was conducted using the Life technologies (ABI) TargetSeq enrichment kit and the SOLiD 5500XL sequencing system at the Center for Targeted Therapy (University of Texas MD Anderson Cancer Center). Variants were called using LifeScope Genomic Analysis Software (ThermoFisher), the CLC Genomics Workbench (Qiagene), and the Genome Analysis Toolkit (25). Based on the ANNOVAR annotation (26), intergenic, intronic, synonymous variants, and variants with Ͼ 1% allele frequency in the ExAC data set and/or Ͼ 5% allele frequency in the 1000 Genome Project data set were removed. Candidate disease causal variants were selected for further analysis by Sanger sequencing.

Enrichment of Organelles in Subcellular Fractions-We
conducted four independent experiments to assign protein location. These represent biological replicates of two different experiments, the first of which was to generally assign organellar location whereas the second was to provide additional data to distinguish between proteins residing in peroxisomes and lysosomes.
The initial fractionation steps were based on the classical differential centrifugation method developed by de Duve and coworkers (6). Here, a rat liver homogenate (H) is separated into five fractions denoted as N (nuclear fraction), M (heavy mitochondrial fraction), L (light mitochondrial fraction), P (microsomal fraction) and S (final supernatant) (14). We modified this procedure by subjecting the L fraction to a second centrifugation step to produce an L1 and L2 fraction (Fig.   1A). Marker enzyme analysis by colorimetric or fluorometric in vitro assays revealed that mitochondria, ER, PM and cytosol have distinct distributions from each other and from lysosomes and peroxisomes, whereas the latter two organ- elles have similar distributions (Fig. 1B). Note that it is axiomatic that in subcellular fractionation, organelles are depleted or enriched in the various fractions but are never purified to homogeneity.
Lysosomes and peroxisomes are difficult to distinguish by differential centrifugation thus we used two different approaches to separate these organelles. One exploits intrinsic differences in the densities of organelles in a Nycodenz gradient (15). Here, L1 (the fraction in which peroxisomes and lysosomes are the most enriched) was further fractionated by centrifugation on a bottom-loaded Nycodenz step gradient (Fig. 1A). Marker enzyme analysis revealed that lysosomes and peroxisomes were enriched and depleted, respectively, at the interface between the 1.105 and 1.135 g/cm 3 density solutions (Fig. 1B, Nycodenz Centrifugation, fraction Nyc2). The other approach relies on the selective decrease in the density of liver lysosomes when rats are administered Triton WR 1339 ("Triton") (16). Triton is a lipoprotein lipase inhibitor, inducing a marked increase in circulating triglyceride-rich lipoproteins which are subsequently endocytosed and accumulate in liver cells (27). L1 fractions from control or Tritontreated rats were dispersed into a sucrose solution (final density 1.18 g/cm 3 ), centrifuged, and the tube divided into top and bottom fractions that were collected separately (Fig. 1A). Marker enzyme analysis revealed that lysosomes are preferentially shifted into the top fraction by the Triton treatment (Fig. 1B, Sucrose Centrifugation).
Protein localization was then assigned based on isobaric labeling/quantitative MS analysis of selected fractions from the four experiments (Table I). There are intrinsic properties of isobaric labeling experiments that may introduce stochastic and systematic errors and we addressed these in the following manner: 1. Coisolation of labeled background peptides prior to MS2 introduces background signal leading to reporter ion compression and reduction of dynamic range (28). To address this, for most LC-MS runs, we prefractionated peptide mixtures by strong cation exchange and/or alkaline reverse phase chromatography to reduce sample complexity and thus background peptides. Monitoring of an internal standard added in different amounts to each sample revealed that there was some ratio compression but the response was linear (supple-mental Fig. S1), with r 2 ranging from 0.95 to 0.99 for the four experiments.
2. Samples are subjected to multiple experimental procedures before pooling and each step can potentially result in differential digestion and labeling efficiencies as well as different degrees of deamination and oxidation. To address this, we only include spectra assigned to fully tryptic peptides with no missed cleavages, fully labeled at expected residues (amino termini and lysines), and lack of adventitious modifications other than methionine monooxidation.
3. To help eliminate noise and quantification errors, we used classic balance sheet analysis to compare the amount of material represented by each spectrum that was present in all fractions with respect to starting material, and only included those with recoveries ranging from 2/3 to 3/2 (Experimental Procedures and supplemental Fig. S2).
Of the ϳ1.8 million spectra assigned to rat peptides in the four experiments, ϳ900,000 were considered quantifiable based on our inclusion criteria (supplemental Table S1). These mapped to a total of 9871 different proteins, which was reduced to 6920 when only considering those represented by two or more peptides and at least three spectra within a given experiment (Table I).
Protein Assignment to Subcellular Organelles-Expts A and B consisted of proteomic analyses of all differential centrifugation fractions and the initial rat liver homogenates and were used to establish primary subcellular localization of different proteins. A single lysosome-enriched fraction from a Nycodenz gradient ("Nyc2") conducted on the L1 fraction was included in Expt A. Note that the Nycodenz fractionation was not performed in Expt B and thus lysosomes and peroxisomes are not distinguishable in this experiment. Expts C and D consisted of a proteomic analysis of a density-based separation of the L1 fraction from control and Triton-treated rats, and were used to help properly distinguish between lysosomal and peroxisomal localizations. It is worth elaborating on the methods used for assignment of protein location to allow researchers to understand and make best use of this information.
We used a custom approach that we term "constrained proportionate assignment" (CPA) adapted from (20) in Expts A and B to estimate the probability that a given protein is located within a given compartment as well as its distribution among all compartments. CPA is described in detail under Experimental Procedures but given its importance in our analysis, it is worth briefly reviewing the overall method. The basis of CPA is quite simple. If a protein was present in two or more different compartments, its distribution among the centrifugation fractions would reflect the sedimentation properties of the compartments, proportionally weighted to the residencies within these compartments. For our analysis of Expt A, we postulated that each protein resides in one or more of the eight major subcellular compartments: mitochon-dria, lysosomes, peroxisomes, ER, Golgi, PM, nucleus, and cytosol. For each quantifiable spectrum, analogous to the relative specific activity term used in classical subcellular fractionation marker enzyme analysis, we calculated enrichment factors for N, M, L1, L2, P, S and Nyc2. These data were used to create a profile for each protein consisting of its enrichment factor in each fraction. Based on published studies, we selected sets of 4 -6 proteins with rigorously established cellular localization as archetypal markers for each of the eight cellular compartments (supplemental Table S2) and used these to create reference profiles ( Fig. 2A, Table 2) determined from MS data. Black-yellow dashed lines, consensus profiles for indicated compartments. Other lines, data from marker enzyme activity assays. A: black solid lines (mitochondria, cytochrome oxidase; lysosome, ␤-galactosidase; peroxisome, catalase; ER, alkaline ␣-glucosidase; PM, alkaline phosphodiesterase; cytosol, dipeptidylpeptidase III) or in black dashed lines (lysosome, TPP1; peroxisome, urate oxidase; ER, NADH cytochrome c reductase). B: as above, except in the lysosome/peroxisome panel, the peroxisomal maker catalase is denoted by a blue solid line and urate oxidase was not quantified. dashed lines). In general, the profiles for each compartment were distinct from each other and consistent with respective marker enzyme analyses (Fig. 2). Each protein was then fit to an optimal linear combination of the eight reference profiles, with optimality defined by minimizing the sum of squares of differences between the profiles of the protein and of the references. The eight "classification" coefficients (each constrained to values ranging from 0 to 1 and summing to 1) correspond to a proportionate assignment of residencies of each protein to the eight compartments, with single or multiple residencies possible. Similar analysis was conducted with Expt B except that because we used a seven-compartment distribution, not distinguishing between lysosomal and peroxisomal locations (Fig. 2B). The coefficient of this "Lyso/Perox" classifier would be comparable to the sum of the lysosomal and peroxisomal classification coefficients from Expt A.
It is important to note that accuracy of classification depends both on the quality of the quantitative MS data and the degree of separation of different compartments by subcellular fractionation. To estimate data quality, for proteins represented by at least three spectra, we used variance in the data and a bootstrap procedure to generate confidence intervals for each classification coefficient. For degree of separation by subcellular fractionation, inspection of the marker protein profiles indicates good resolution of most organelles from each other with some exceptions. In particular, ER, Golgi, and PM are only subtly different, which may lead to some uncertainty on whether a protein is assigned properly among these three compartments, especially when considering proteins with multi-compartment distributions.
Comparison of the classification coefficients for proteins quantified in both Expts A and B indicates a strong positive correlation (Fig. 3, all, r 2 ϭ 0.79). Not unexpectedly, Golgi and PM show the greatest scatter. It is likely that this reflects intrinsic difficulties in distinguishing ER, Golgi, and PM by differential centrifugation alone rather than differences in underlying protein distributions in the two experiments. If one sums the values of the ER, Golgi, and PM classification coefficients to create a microsomal classification, the overall correlation improves (Fig. 3, All*, r 2 ϭ 0.88).
There is good agreement when comparing the predominant compartmental classification for proteins quantified in both Expts A and B (Table II). There is 82.4% concordance when considering all overlapping 2952 proteins. Overall agreement increases to 94.9% when comparing the 2133 proteins where the classification coefficient point estimate is Ն0.5. Concordance improves to 97.3% for the 1877 proteins with classification coefficient lower confidence limits Ն0.5, rises to 99.8% for the 1277 proteins with point estimates Ն0.7, and reaches 100% for the 1094 proteins with lower confidence limits Ն0.7. Thus, the greater the stringency, the more confident the assignment of predominant compartmental localization of a given protein, at the expense of overall proteome coverage. Fig. 4 illustrates separations achieved in Expts A and B for individual proteins using principal component plots, with symbols indicating predominant compartment assignment for those with point estimates Ն0.7. This also can be illustrated using heatmaps (supplemental Workbook S1, worksheets HeatA and HeatB). Note that plots A and B appear different  (Table I) are plotted against each other and fit using linear regression. For LysoϩPerox, the sum of the two classification coefficients in Expt A is plotted against the combined lysosome/peroxisome classification coefficient from Expt B. For ERϩGolgiϩPM (microsomal distribution), the sum of the three classification coefficients is plotted against each other. For All*, the sum of the ER, Golgi and PM coefficients are used instead of the individual coefficients.
because Expt A contains an additional dimension that contributes to the principal components (the Nyc2 fraction), which also allows lysosomes and peroxisomes to be classified independently. However, in both Expts A and B, proteins with the same predominant classification coefficient from the CPA analysis cluster, although not unexpectedly, ER, Golgi, and PM proteins are not well resolved. Support vector machine clustering yielded single compartmental assignments that agreed well with the predominant compartment assigned in our CPA procedure (data not shown). Similar agreement was found using hierarchical clustering. This underscores the utility of the CPA method, which, in addition to agreeing with more conventional procedures in assigning proteins to a major compartment, also provides estimates for proportional residence in multiple compartments.
In Expts C and D (supplemental Workbooks S2C and S2D), the Triton shift approach was conducted specifically to distinguish lysosomal from other proteins present in the L1 frac- tion. Fig. 5A shows the relative distribution of generally accepted lysosomal, peroxisomal and mitochondrial proteins after centrifugation in 1.18/cm 3 sucrose. All the lysosomal proteins show a dramatic shift in the Triton-treated samples. Triton-treatment had little effect on peroxisomal proteins. In contrast, there was a small effect on mitochondrial proteins. The basis for this is unclear but one possibility is that this is because of the normal process of mitochondrial turnover and degradation within the lysosome. It also is important to note that for multicompartmental proteins with lysosomal component and a nonlysosomal, nonperoxisomal component, the Triton shift experiments will overestimate the lysosomal proportion of the total population.
We explored different methods to assign proteins based on Triton shift experiments. One possibility was to fit the data to the same set of rigorously established reference proteins used for Expts A and B (Fig. 5B). We also considered using proteins that had consistent classification in Expts A and B with various thresholds for inclusion. The proteins assigned with a lower confidence limit of Ն 0.7 in both Expts A and B (Fig. 5C) have a pattern quite similar to that of the established markers (Fig. 5B). Given that the former are much more numerous, we used these proteins for purposes of classification. Here, we divided the proteins into a lysosomal marker set and a marker set encompassing all other compartments. We then used finite Gaussian mixture clustering (21) (Experimental Procedures) to calculate "Lyso" and "Other" classification coefficients for all proteins in each experiment, generating confidence limits for proteins with at least three spectra.
We also compared the fractionation of all 1923 proteins common to both Expts C and D (Fig. 5D). There was good correlation in the distribution found in fasted and fed animals regardless of Triton treatment, indicating minimal overall effects of nutritional status. Because of the greater protein coverage, we rely primarily on Expt C for corroborating lysosomal assignment in Expt A and for distinguishing between lysosomal and peroxisomal residence in Expt B.
Benchmarking protein assignments-While our assignments were largely internally consistent, it was critical to Fig. 5. Analysis of Triton-shift experiments. A, Distribution of select lysosomal, peroxisomal and mitochondrial proteins following sucrose density centrifugation of L1 fractions from control (C) and Tritontreated (T) rats. Proteins were selected as having Ն4 peptides in Expts C and D that have been localized to mitochondria (from MitoCarta v2.0 (29) with a FDR Յ0.001); lysosomes (from compendium reviewed in (53) omitting MPO (54)); and peroxisomes (PeroxisomeDB2.0 (32) omitting SOD1(17)). B, Distribution of proteins with at least 2 peptides and 3 spectra (gray points) with other symbols indicating marker proteins from supplemental Table S2. C, As in Panel B, with other symbols indicating assignments from Expt A for proteins with a predominant compartment lower confidence interval Ն 0.7 in Expts A and B. Lines indicate boundaries for proportional assignment between lysosomal and other locations (dashed line 1:1; dotted lines, 2:1 and 1:2). D, Comparison of Expts C and D. benchmark them against proteins with established locations. For this purpose, we selected a "high stringency set" of 3034 proteins whose major compartment was consistent in Expts A, B and C, with a minimum of 2 peptides and 3 spectra for Expt A and at least one peptide and spectrum from Expts B and C (supplemental Workbook S1). We compared these to the Compartments Database Benchmark Set (23) (CDbBS), which lists locations for 12,892 human proteins. Of these, we found 809 rat orthologs in our high stringency set that were annotated in CDbBS to be associated with only one of our eight compartments. As shown in Fig. 6A, there was moderately good concordance for these proteins. However, given that there is no "gold standard" for compartmental localization, we conducted a blinded analysis to evaluate the published experimental evidence supporting the CDbBS categorizations (see Experimental Procedures). This allowed us to discard questionable entries and create a "curated" database containing 588 proteins. Our major assignments now agree extremely well with those in this high stringency curated CDbBS (Fig. 6B). This not only provides additional confidence in results from analytical subcellular fractionation but also underscores the inaccuracies present in existing subcellular localization databases, of which further examples are described below.
Comparing Prolocate Protein Assignments with Organelle Databases-We have compared our data against several commonly used protein localization data sets (see below). However, investigators can readily compare our assignments with any data set of interest using simple software tools (see Supplemental Materials, "Comparing custom data sets with Prolocate". The MitoCarta study utilizes a combination of MS, GFPtagging and/or bioinformatics criteria to estimate the likelihood in terms of a false discovery rate (FDR) for mitochondrial localization of Ͼ20,000 different mammalian gene products (29,30), 2657 of which are present in our high stringency data set. For these, our Mito classification coefficient decreased with increasing MitoCarta FDR (Fig. 7A). The MitoCarta study also contains a list of proteins that are designated as having strong support of mitochondrial localization, 621 of which are in our high stringency set. Most of these have low MitoCarta FDRs, but some do not, and it is likely that some represent false positives. Our mitochondrial classification coefficients are considerably higher for the designated mitochondrial Mi-toCarta proteins with low FDRs compared with those with higher FDRs (Fig. 7B and 7C, see also supplemental Workbook S1, worksheet heatMITOCARTA). Use of both databases in tandem should provide increased confidence in mitochondrial assignments, as well as identifying possibly incorrectly assigned proteins.
The MitoMiner database is a searchable database that has compiled an integrated mitochondrial protein index (IMPI) listing putative rat, mouse and human mitochondrial proteins (31). Our high stringency data set contains 998 of the proteins on the rat IMPI list. Only half of these are classified as predominantly mitochondrial in Expt A, whereas for the others, many have Mito classification coefficients near zero (Fig. 7D, see also supplemental Workbook S1, worksheet heatIMPI). These data suggest that there are a number of putative mitochondrial proteins on the IMPI list that are worthy of further evaluation.
The Peroxisome Database (32) (peroxisomeDB) is a catalogue of the peroxisomal proteome of multiple species, including rat, compiled from experimental literature review, homology and bioinformatic annotation. Fifty-eight rat proteins in the peroxisomeDB were present in the Prolocate high stringency data set and a majority were assigned a primary peroxisomal localization. However, for a significant proportion of these (19 proteins), we assigned primary localization to compartments other than the peroxisome (Fig. 7E, see also supplemental Workbook S1, worksheet heatPEROXDB). Most of these were assigned a minor peroxisomal location in Prolocate and this is consistent with published evidence for residence in multiple cellular locations. However, seven had a peroxisomal classification coefficient of Ͻ0.01 in Expt A and these warrant further scrutiny regarding intracellular location. Similar results were obtained when our data were compared with a single study (33) focusing on proteomic analysis of peroxisomes (data not shown).
The CLEAR (Coordinated Lysosomal Expression and Regulation) network is a group of genes that are predicted to regulate lysosomal biogenesis and function via the transcription factor EB (TFEB) (34,35). Members of the CLEAR network are characterized by the presence of binding sites for TFEB within their respective promoters. The Prolocate high stringency data set contains 159 proteins encoded by genes containing CLEAR elements. For 50 of these, we report a primary lysosomal residence but the majority (109 proteins) are assigned to other organelles (Fig. 7F, see also supplemental Workbook S1, worksheet heatCLEAR) although a minor lysosomal residence was detected for some. A primarily nonlysosomal location for proteins involved in lysosome biogenesis or function is quite conceivable and consistent with this, our data indicate that the presence of a CLEAR element should not be interpreted as an indicator of lysosomal residence.
Other studies have combined MS with subcellular fractionation to determine global protein localization. While these did not estimate proportionate assignments among all the compartments, and also had some differences in definition of compartments, we were able to conduct some comparisons. Note that we compare the Prolocate high confidence data set to all protein orthologs in the other studies (supplemental Tables S3-S5), but in the summary below, we only report final numbers for proteins found in analogous compartments. An early study used rate-zonal centrifugation and protein profile correlation analysis to classify proteins in a mouse liver homogenate (10). Of the 379 mouse liver proteins assigned to a single compartment, 93.6% agreed with our predominant classification (supplemental Table S3). However, lysosomes and peroxisomes were not part of the mouse liver classification scheme, and the protein correlation profile classifications were heavily biased toward mitochondrial and cytosolic proteins (344 assignments). Two additional studies were published while this manuscript was under initial review. One primarily used isopycnic centrifugation and a ten-channel iso- baric labeling MS3 approach (termed "Hyperlopit") to map proteins in mouse pluripotent stem cells (11). Of the 435 overlapping proteins, 93.3% agreed with our assignments (supplemental Table S4). The latest study analyzed the human HeLa cell line by differential centrifugation and SILAC labeling to map proteins present in an organelle-enriched fraction (essentially corresponding to a combined MϩLϩP fraction), and also used a label-free approach to estimate proportional residence among the "N", "MϩLϩP", and "S" equivalents (12). When considering the 687 HeLa proteins assigned to individual organelles with very high, high and medium confidence, 95.6% agreed with the Prolocate assignments (supplemental Table S5). This increased to 97.3% when using 548 HeLa proteins assigned with very high and high confidence (supplemental Table S5). This remarkable degree of agreement lends confidence to these resources, not only for the proteins identified in multiple studies, but also for the many proteins unique to each.
Applications of Prolocate-Proteins with Complex Distributions-One advantage of our approach is that we estimate the entire distribution of a protein among eight different cellular compartments. We have also generated matrices listing distances between all protein pairs identified in each experiment (supplemental Workbooks S4A-S4D), and the Prolocate interface allows integration of the data to find proteins with similar distributions to any protein of interest measured in one or more experiment. Two select examples are as follows: COG Complex-The COG complex comprises 8 individual subunits that together mediate retrograde vesicle transport within the Golgi. We identified all 8 subunits (COG1-8) and, consistent with the known properties of this complex, they exhibited an extremely similar pattern of localization, with ϳ30% in the Golgi and ϳ70% in the cytosolic fraction (Fig.  8A). This agrees with morphological (36) and subcellular fractionation (37) studies that localize members of the COG complex to membrane-associated Golgi and cytosol, respectively. Remarkably, when we search for proteins that colocalize with COG1 using the Prolocate distance calculator (Expts A and C), COGs 3, 4, 5, 7, and 8 represent 5 out of the top 6 highest scoring hits.
AP Complexes-We also investigated the distribution of individual proteins in AP complexes, which are multi subunit complexes involved in different vesicular transport pathways (38). Based on our data, intracellular distribution of the individual proteins of AP complexes 1, 2, 3, and 5 were remarkably consistent (Fig. 8B). Proteins of the AP1 complex were primarily localized to the Golgi and the cytosol, which is consistent with the role of AP1 in trans-Golgi to endosome trafficking. Proteins of the AP2 complex were restricted to the PM and cytosol, with no lysosomal localization, consistent with the role of AP2 in clathrin-dependent endocytosis. Proteins of the AP3 complex were ϳ70% localized to the lysosome, with the remainder being cytosolic. AP3 plays a role in the biogenesis of lysosome and related organelles e.g. endosomes. Distribution of the proteins of AP5 was more heterogeneous than that of the other AP complexes but the lysosome was the primary location. The function of AP5 is unclear but it may play a role in endosomal sorting. As with the COG complex, using the Prolocate distance calculator to identify proteins that colocate with individual AP complex components primarily returns the other complex members as the highest scoring hits.
Candidate Disease Genes-Knowledge of protein location is useful in investigating disease pathways and we present one example here that arises from our long-term interest in lysosomal diseases of unknown etiology (39 -41). In the course of this research, we performed whole exome sequencing on a number of unsolved cases accumulated in the course of our studies but found potential mutations in an unwieldy number (thousands) of genes. We thus used Prolocate to prioritize candidate disease genes based on lysosomal location of the respective encoded proteins.
One of the proteins on our lysosomal candidate list, SLC31A1, also called CTR1, is involved in copper transport (42). The recent HeLa cell study (12) localizes SLC31A1 to the plasma membrane whereas previous reports localized it to the PM with some evidence for an intracellular pool in different cell lines (43). More recently, evidence has been reported that may or may not support some lysosomal residence (44,45). We assign SLC31A1 to both the lysosome and to the plasma membrane (Fig. 8C).
Initial whole exome sequencing data highlighted SLC31A1 as the primary candidate in one of our potential neurodegenerative lysosomal disease subjects, an infant who died in 1988 at age 6 months. Sanger sequencing revealed compound heterozygosity for two potential mutations in SLC31A1 (Fig.  9). A nonsense mutation at Arg90 results in a severely truncated protein and it is presumably a null allele. A missense Val181Leu mutation is of unclear significance but predicted to be deleterious by Sorting Intolerant From Tolerant (SIFT) analysis software (46).
Several observations support the possibility that SLC31A1 deficits may be the cause of disease in this case. First, neither of these changes was found in the large scale Exome Aggregation Consortium sequence of 60,706 individuals (47) indicating they are not common polymorphisms. Second, a homozygous mouse knockout for SLC31A1 is embryonic lethal whereas the heterozygote displays tissue-specific abnormalities including reduced copper levels in brain (48). Thus, the nonsense mutation, in combination with a hypomorphic allele, could be responsible for disease. Third, a missense mutation (Arg90Gly, coincidentally at the site of the nonsense mutation reported here) in SLC31A1 has been suggested to be associated with recessive cognitive disease (49). Fourth, defects in other proteins involved in copper transport e.g. a coppertransporting ATPase in Wilson disease, results in neurologic and other abnormalities (50). Thus, although we do not have access to pedigree and detailed case history information, there is intriguing evidence for SLC31A1 as a potential human disease gene and it warrants further investigation. Use of the Prolocate database to filter the whole exome sequencing data was a key step in identifying this candidate. DISCUSSION We report here our progress to date in establishing a global atlas of protein location. Our approach differs from other quantitative MS efforts in the field in several important respects. First, we use balance sheet analysis to filter data for acceptable recovery. Second, we have implemented statistical approaches that allow estimation of the distribution of each protein among multiple compartments. Third, we estimate confidence limits for protein distributions. This is important as there are experimental errors inherent in any quantitative measurement and this is particularly true for datadependent MS, where there is oversampling of abundant proteins and occasional incorrect assignments of spectra to peptides and of peptides to proteins. Fourth, evidence for location of individual proteins of interest should be scrutinized closely and to facilitate this, the Prolocate website allows inspection of the fractionation profiles of all component peptides assigned to each protein. Profiles of individual spectra, peptides and proteins have been archived on the Prolocate Website and in supplemental Workbooks 2A-2D.
It is worthwhile noting that our analysis depends on assigning proteins to reference compartments but they do not cover all cellular structures (e.g. different types of endosomes, stress granules, lipid droplets). If a different cellular structure fractionated similarly to one of the reference compartments, proteins in that structure would be assigned to that reference compartment. Alternatively, if a protein resided in a single structure that fractionated distinctly from the eight reference compartments, it would be assigned to multiple reference compartments. In such cases, the organelle assignments may be misleading but identification of proteins with similar fractionation properties (e.g. using the distance calculator) may still reveal potential associations with a common compartment or sets of compartments. Expanding the number of fractionation procedures used as well as use of additional compartment and subcompartment markers should allow incorporation of additional cellular structures into future classification schemes. Nonetheless, concordant results from a series of complementary approaches including analytical fractionation and microscopy are required for the lowest level of uncertainty in assignment.
We use rat liver as a source for analytical centrifugation because of effective fractionation procedures tailored to this tissue (6), and inability to achieve equivalent high-quality fractionation using mouse liver (51). In addition, fresh tissue is highly relevant to in vivo physiology. However, although hepatocytes make up the bulk of the liver, there are additional cell types present, and it is possible that this introduces additional heterogeneity in the centrifugation properties of organelles. Cultured cells may provide a more homogeneous source, with the tradeoff that cells undergo changes in culture with accompanying perturbations in protein trafficking, even after short time periods (52). Nonetheless, the remarkable concordance between our assignments on rat liver proteins with those from studies on mouse pluripotent stem cells (11) and HeLa cell (12) proteomes, suggests that such concerns are not a serious obstacle in assigning locations to most proteins. It is worth noting that the use of different cell sources (e.g. rat and mouse liver, cultured mouse and human cells) and different approaches to fractionation (differential, velocity gradient, isopycnic centrifugation) will results in data sets with some overlap and, thus provide independent validation, but many proteins will be unique to each study. Consolidating these indi-vidual studies will be an important step toward establishing a high-quality map of the mammalian subcellular proteome.
In summary, we provide subcellular localization data for over 6000 rat liver proteins, representing approximately a third of the predicted liver proteome (13). This Prolocate database is readily accessible through a web portal. Prolocate can be used to query the location of individual proteins and access the supporting data, and to generate lists of candidate proteins for different cellular compartments. A key feature of Prolocate is the ability to identify sets of proteins that have similar single-or multi-compartmental distributions that may reflect functional interactions within the cell. While future studies will extend coverage of the proteome, explore the relationships between subcellular location and post-translational modification, and address additional cellular compartments, the current map and analytical tools should provide a valuable resource for the biomedical research community. To determine if these were on same or different copy of the SLC31A1 gene, long range PCR was conducted to amplify 2.1 kb fragments (forward primer, CAAGCAGTCTGACCAAAAGGT; reverse primer, CTTCAACAACTTCCCACTGCA) containing both regions of interest. After subcloning, sequencing of individual plasmids indicated compound heterozygosity.