Role of S-Palmitoylation by ZDHHC13 in Mitochondrial function and Metabolism in Liver

Palmitoyltransferase (PAT) catalyses protein S-palmitoylation which adds 16-carbon palmitate to specific cysteines and contributes to various biological functions. We previously reported that in mice, deficiency of Zdhhc13, a member of the PAT family, causes severe phenotypes including amyloidosis, alopecia, and osteoporosis. Here, we show that Zdhhc13 deficiency results in abnormal liver function, lipid abnormalities, and hypermetabolism. To elucidate the molecular mechanisms underlying these disease phenotypes, we applied a site-specific quantitative approach integrating an alkylating resin-assisted capture and mass spectrometry-based label-free strategy for studying the liver S-palmitoylome. We identified 2,190 S-palmitoylated peptides corresponding to 883 S-palmitoylated proteins. After normalization using the membrane proteome with TMT10-plex labelling, 400 (31%) of S-palmitoylation sites on 254 proteins were down-regulated in Zdhhc13-deficient mice, representing potential ZDHHC13 substrates. Among these, lipid metabolism and mitochondrial dysfunction proteins were overrepresented. MCAT and CTNND1 were confirmed to be specific ZDHHC13 substrates. Furthermore, we found impaired mitochondrial function in hepatocytes of Zdhhc13-deficient mice and Zdhhc13-knockdown Hep1–6 cells. These results indicate that ZDHHC13 is an important regulator of mitochondrial activity. Collectively, our study allows for a systematic view of S-palmitoylation for identification of ZDHHC13 substrates and demonstrates the role of ZDHHC13 in mitochondrial function and metabolism in liver.

thiopropyl Sepharose. Compared to ABE, this strategy is simpler and the captured peptides are more acceptable to MS owing to the absence of biotin labelling; hence, this method has increased the sensitivity for the identification of palmitoylated peptides. However, palmitoylation sites are difficult to determine unambiguously when a peptide contains two or more cysteine residues. Metabolic labelling with a palmitic acid analog followed by click chemistry (MLCC) has been applied in palmitoylated protein research 16,17 . Nevertheless, MLCC has been used in mammalian cells rather than in live animals, and only capture proteins in the metabolic labelling time course.
To date, palmitoylation has been investigated partly because PATs have been linked to several human conditions, including Huntington's disease (ZDHHC13 and ZDHHC17) 18,19 , schizophrenia (ZDHHC8) 20 , type 1 diabetes (ZDHHC17) 21 , intellectual disability (ZDHHC9 and ZDHHC15) 8,22 , and multiple types of cancers (ZDHHC2, ZDHHC9, ZDHHC11, ZDHHC14, and ZDHHC20) 23 . Furthermore, mice with an F233 deletion in ZDHHC21 show abnormalities in skin homeostasis and hair defects 24 . To understand the function of disease-associated PATs, it is critical to link enzymes to their corresponding substrates. To date, only six proteomics studies have focused on identifying substrates of specific PATs, such as ZDHHC2, ZDHHC5, and ZDHHC17 [25][26][27][28] . We previously reported that ZDHHC13 is a novel regulator of postnatal skeletal development and bone mass acquisition via palmitoylation of MT1-MMP 29 . Moreover, ZDHHC13 is important for hair anchoring and skin barrier functions, likely mediated by the ZDHHC13 substrate, cornifelin 30 . Previously, we generated a Zdhhc13-deficient mouse model through N-ethyl-N-nitrosourea (ENU) mutagenesis 31 . The mutation creates a truncated protein (R425X) that lacks the zinc-finger DHHC-CRD domain (positions 426-476) and the active site (C456) for the formation of an S-palmitoyl cysteine, resulting in nonsense-mediated mRNA decay of Zdhhc13 mRNA. Compared with other PAT family member-defective mouse models, Zdhhc13-deficient mice exhibit severe phenotypes, including alopecia, osteoporosis, and amyloidosis. In this study, we aimed to investigate the metabolic abnormalities, including hypermetabolism, and defective lipid metabolism of Zdhhc13-deficient mice. Furthermore, given that few ZDHHC13 substrates have yet been identified, we studied the liver S-palmitoylome using a quantitative proteomics approach that combined the alkylating RAC method and an MS-based label-free strategy to identify candidate ZDHHC13 substrates.

Results
Abnormal liver function, lipid abnormalities, and hypermetabolism in Zdhhc13-deficient mice. In addition to our previous report that Zdhhc13-deficient mice presented with osteoporosis, alopecia, skin abnormalities, and amyloidosis 31 , the mice also exhibited elevated liver transaminase. In this work, we found that these mice had low blood triglyceride and cholesterol levels, as well as decreased lipid accumulation in the liver (Fig. 1a,b). We further observed that these mutant mice displayed excessive sweating and poor weight gain despite adequate water and food intake. Therefore, we examined the metabolic state and efficiency of Zdhhc13-deficient mice and found that VO 2 consumption, VCO 2 production, and heat production were significantly increased in the mutant mice compared to those of wild-type (WT) mice (Fig. 1c). However, serum levels of T4 and T3 were not significantly different in WT and mutant mice (T3: WT: 1.05 ± 0.24 ng/ml, mutant mice: 1.27 ± 0.32 ng/ml; T4: WT: 4.6 ± 0.43 ng/ml, mutant mice: 4.3 ± 0.17 ng/ml), suggesting that this hypermetabolism was not due to hyperthyroidism.
Identification and characterization of the liver S-palmitoylome. To explore the pathological mechanisms underlying hypermetabolism and defective lipid metabolism in Zdhhc13-deficient mice, we analysed tissue samples from the liver, a key metabolic organ. For site-specific identification and quantitation of the S-palmitoylome and systematic delineation of S-palmitoylated substrates of ZDHHC13 from WT and Zdhhc13-deficient mice, we used a simple multiplex strategy that combines an alkylating RAC method and an informatics-assisted label-free strategy 32 (Fig. 2a). Using alkylating agents for blocking free and disulfide cysteines, the false positive signals that originate from protein or reagent disulfide exchange can be minimized. With the covalent tagging of these cysteines, S-palmitoyl groups on proteins were reduced by hydroxylamine, generating a reducing form of the cysteine residue. Followed by thiol-resin capture and eluted by TCEP, unambiguous palmitoylation sites can be specifically identified as reduced cysteine residues. To further increase identification coverage and decrease quantitation of heterogeneity from individual mice, we also utilized peptide alignment to discover signals of un-identified S-palmitoylated peptides in different LC-MS/MS replicates based on the corresponding retention time and m/z of identified peptides. Finally, all S-palmitoylated peptides were relatively quantified with extracted ion chromatogram (XIC). In addition, expression levels of membrane proteins were quantified by TMT10-plex labelling for normalization.
In total, 2,190 unique S-palmitoylated peptides (1,905 sites) corresponding to 883 proteins were identified from three pairs of WT and Zdhhc13-deficient mice (mascot score ≥25, p < 0.05, Table S1), including well-known palmitoylated proteins such as CD82, CPS1, Gnas, Rras, Gphn, TMX, Ctnnd2 and M6pr. Among these, 311 (35%) were novel S-palmitoylated proteins that had not been previously identified (Table S2). These results demonstrate that alkylating RAC combined with peptide alignment of the liver S-palmitoylome provides greater sensitivity for identification of palmitoylated proteins than that of the previously used ABE system. To explore the potential functional roles of these 883 S-palmitoylated proteins in mice, subcellular localizations and functions were annotated using the Ingenuity Pathway Analysis (IPA) Knowledge Base and the Gene Ontology (GO) consortium. As shown in Fig. 2b, 50% of identified S-palmitoylated proteins are located in the membrane fraction, including the plasma membrane (17%), mitochondrial membrane (15%), endoplasmic reticulum (10%), Golgi apparatus (1%), and other membrane (7%). In addition, 31% of S-palmitoylated proteins are located in the cytosol, 11% in the nucleus, and 3% in the extracellular space. The functional categorization of these S-palmitoylated proteins revealed that most proteins are enzymes (54.4%) and many are involved in metabolism (Fig. 2c). In addition, some S-palmitoylated proteins were categorized as transporters (12.5%), translation regulators (2%), transcription regulators (1.5%), and transmembrane receptors (1%).

Identification and characterization of candidate ZDHHC13 substrates. Candidate ZDHHC13
substrates were identified as S-palmitoylated peptides that were significantly down-regulated (≤0.5-fold) in triplicate samples of Zdhhc13-deficient mice vs. WT mice. In order to remove the effects of changes in protein expression on S-palmitoylation levels, we used the corresponding protein level of each protein to normalize the S-palmitoylation levels. Following quantitation of protein levels with TMT10-plex labelling, expression levels of approximately 96% of proteins (3,117 of 3,246) were unchanged, with fold-changes between 0.71 and 1.41 (average ratio = 1.06 ± 0.34 with 95% confidence). Based on the biological functions and subcellular localizations, some outlier proteins as well as general contaminant proteins from immunoprecipitation or affinity capture (such as ribosomal proteins, refer to CRAPome database: http://www.crapome.org/) were filtered. The remaining 254 proteins corresponding to 400 S-palmitoylation sites (369 peptides) were annotated as potential ZDHHC13 substrates (Table S4).
We further studied the canonical pathway and functional annotations by IPA and found that many of these proteins are associated with liver-related disorders and metabolism, especially fatty acid metabolism and mitochondrial dysfunction ( Table 1, Fig. 3a,b). Of these proteins, 46% are involved in metabolism, including MCAT, SDHA, NDUFV1, NDUFS1, FASN, ACAA1, and ACAA2 (Table S5). A systematic network of protein-protein interactions showed that these candidate ZDHHC13 substrates are associated with energy production, metabolic diseases, and the maintenance of cellular function (Fig. 3b). In addition, we identified many enzymes involved in lipid metabolism. In order to eliminate the possibility that these enzymes form covalent acyl-enzyme intermediates during catalysis, we confirmed the palmitoylation of three proteins, ACAA2, MCAT, and CTNND1, by MLCC and verified that they are indeed palmitoylated proteins (Fig. S1). Taken together, the quantitative study of the S-palmitoylome indicated that ZDHHC13 and S-palmitoylated proteins that are down-regulated in Zdhhc13-deficient mice play important roles in metabolism and mitochondrial dysfunction.
Decreased palmitoylation of candidate ZDHHC13 substrates in Zdhhc13-deficient mice. We chose nine candidate ZDHHC13 substrates based on their functions in mitochondria and lipid metabolism, as well as their potential to cause metabolic abnormalities, to further validate S-palmitoylation levels using western blot analysis. SDHA, NDUFV1, NDUFS1, and UQCRC1, are involved in regulating the electron transport chain in mitochondria, while MCAT is involved in mitochondrial fatty acid synthesis, and FASN, ACAA1, ACAA2, and CTNND1 are involved in fatty acid synthesis and β-oxidation. As shown in Fig. 4, these proteins were significantly less palmitoylated in Zdhhc13-deficient mice compared with those of WT mice. These results are consistent
with the proteomics data and suggest that reduced palmitoylation of candidate ZDHHC13 substrates may cause abnormal phenotypes in Zdhhc13-deficient mice.
Confirmation of MCAT and CTNND1 as ZDHHC13 palmitoylation substrates. We further chose two of the above proteins to confirm that they are indeed ZDHHC13 substrates. Previous studies have shown that Mcat-knockout mice exhibited similar phenotypes to those of Zdhhc13-deficient mice, including kyphosis and alopecia 33 . Here, we showed that palmitoylation levels of MCAT were significantly (2-fold) higher in ZDWT cells than basal palmitoylation levels in ZDK cells (Fig. 5a). CTNND1, a known anti-adipogenic protein 34 , has been identified as a palmitoylated protein in the lipid rafts of human prostate cancer cells in a proteomics study 35 . The S-palmitoylation level of CTNND1 was 1.25-fold higher in ZDWT co-expressing cells than in ZDK cells (Fig. 5b). These results confirm that MCAT and CTNND1 are ZDHHC13 substrates.
Our proteomics results identified C104 and C618 as palmitoylation sites in MCAT and CTNND1, respectively (Table 1). To further examine the fidelity of MS-based identification of palmitoylation sites, we generated two mutants, MCAT C104G and CTNND1 C618G, in which the cysteine residues were replaced with glycine. Reduced palmitoylation levels were observed for both mutant proteins compared with those of WT proteins  Table 1. Partial list of candidate ZDHHC13 substrates involved in mitochondrial dysfunction and lipid metabolism. Candidate ZDHHC13 substrates are list in line included the uniport number, protein name, gene name, mass to charge ratio (m/z), peptide charge, peptide sequence which the modification cysteine is shown in lower case letter, S-palmitoylated modification site, Palm/proteome, Palm level ("+" indicated the sample treated with Hydroxylamine), Proteome level. (K indicates Zdhhc13-deficient mice).
( Fig. 5c). These findings not only verify that C104 of MCAT and C618 of CTNND1 are ZDHHC13 palmitoylation sites, but also confirm the fidelity of the alkylating RAC method in detecting novel palmitoylation sites.

Impaired mitochondrial function in hepatocytes of Zdhhc13-deficient mice and
Zdhhc13-knockdown Hep1-6 cells. The systematic proteomics study of ZDHHC13 identified several candidates for S-palmitoylated substrates involved in the regulation of mitochondrial function (Fig. 3); therefore, these proteins may result in impaired mitochondrial function in Zdhhc13-deficient mice. To evaluate the effect of ZDHHC13 on mitochondrial function, primary hepatocytes isolated from WT and Zdhhc13-deficient mice were assay by the Seahorse XF96 extracellular flux analysers. The results showed that the proton-leak OCR was slightly higher in Zdhhc13-deficient hepatocytes than in WT hepatocytes, and the spare respiratory capacity in the presence of 4 μM FCCP was higher in Zdhhc13-deficient hepatocytes than in WT hepatocytes. The basal mitochondrial OCR and ATP-link OCR did not show significant difference (Fig. 6a). We further examined the mitochondrial membrane potential, mitochondrial ROS level, and oxidative damage biomarkers. The membrane potential was assessed by measuring changes in the JC-1 aggregate to JC-1 monomer ratio. The results showed that the membrane potential were slightly lower in Zdhhc13-deficient hepatocytes than in WT hepatocytes (Fig. 6b). The mitochondrial ROS level was assessed by MitoSOX Red staining. The data showed that mitochondrial ROS were obviously increased in Zdhhc13-deficient hepatocytes than in WT hepatocytes (Fig. 6c). For the observation of hypermetabolism and heat production of the Zdhhc13-deficient mice, we examined the protein expression of the uncoupling protein (UCP). The data show that the UCP1 expression in the hepatocyte of the Zdhhc13-deficient mice was higher than in the WT. And the oxidative damage biomarkers assessed by SOD2 and Carboxymethyl Lysine (CML) were all higher in Zdhhc13-deficient hepatocytes than in WT hepatocytes (Fig. 6d). Moreover, based on IPA annotation, 18 candidate ZDHHC13 substrates, including LDHA, PRDX1, and NDUFS1, were identified as inhibitors of ROS production [36][37][38] , and CAT was identified as a ROS scavenger 39 (Table S6). Collectedly, our data support the mitochondrial dysfunction and oxidative stress in the hepatocytes of Zdhhc13-deficient mice. To further investigate the effect of ZDHHC13 on mitochondrial function, we knocked-down Zdhhc13 expression in Hep1-6 cells using shRNA. Zdhhc13 mRNA expression was reduced by 10-30% in four Zdhhc13 shRNA clones (A1, B1, C1 and A5) as compared to expression levels of scrambled control shRNA (A6) (Fig. 7a). The spare respiratory capacity was lower in Zdhhc13-knockdown cells than in the scramble control in the presence of 0.5 μM FCCP (Fig. 7b). Silencing of Zdhhc13 led to lower mitochondrial membrane potential compared with that of the scrambled control (Fig. 7c). As expected, knockdown of Zdhhc13 also led to the production of more mitochondrial ROS (Fig. 7d). In summary, these data support that knockdown of Zdhhc13 disrupts mitochondrial function in Hep 1-6 cells.

Discussion
In this study, we observed Zdhhc13-deficient mice to exhibit hypermetabolism and defective lipid metabolism and attempted to elucidate the molecular mechanisms underlying these disease phenotypes. We applied an alkylating RAC-based label-free quantitative proteomic strategy to systematically identify S-palmitoylated candidate ZDHHC13 substrates in the liver. With good blocking efficiency and an irreversible chemical reaction to reduce S-palmitoylation false positives, a total of 2,190 unique S-palmitoylated peptides corresponding to 883 proteins were identified and 1905 unambiguous S-palmitoylation sites were determined. Among the identified S-palmitoylated proteins, 311 were novel palmitoylated proteins with unambiguous palmitoylation sites that had not been observed previously by other proteomics studies. We not only successfully identified more than twice the number of palmitoylated proteins as had been previously determined in acyl-RAC-based palmitoylome study using MMTS blocking 15 , but we also effectively determined precise S-palmitoylation sites. This technical advance helps overcome the challenge presented by the fact that some palmitoylated proteins exhibit relatively low levels of palmitoylation or have multiple cysteine residues in their peptide sequences. In addition, complete blocking also decreases the false positive signal from free cysteines and from protein-protein or protein-reagent disulfide linkages 40 . When compared with the results of two palmitoyl-proteomics studies involving ZDHHC5 and ZDHHC17, our method identified 254 and 105 of the S-palmitoylated proteins found in these earlier studies, respectively. This indicates that different PATs may share the same substrates. Compensatory effects between differentially expressed PATs may regulate the extensive palmitoylation of proteins. In addition, diverse palmitoylation profiles may be due to tissue-specific expression of PATs. Quantitative proteomics analysis revealed potential candidate substrates of ZDHHC13, as defined by proteins with unambiguous palmitoylation sites where the palmitoylation level was ≥0.5-fold lower in Zdhhc13-deficient mutants. Utilizing the peptide alignment strategy for label-free quantitation, we further increased the identification coverage and decreased heterogeneity from individual mice. Integrating multiplex quantitation with TMT10-plex isotopic labelling, the palmitoylation levels of these potential candidate substrates were further examined by normalizing with each protein's expression level. However, in some cases, decreased palmitoylation can affect protein abundance, as palmitoylation is strongly related to protein stability 3 . Thus, further validation of the palmitoylation-mediated biological functions of candidate substrates should be performed in the future.
Zdhhc13-deficient mice present with metabolic abnormalities, including non-thyroidal hypermetabolism and hyperhidrosis. Compared to WT mice, Zdhhc13-deficient mice consume more food and water but exhibit poor weight gain. Our MS-based systematic analysis demonstrated that 111 of 307 candidate ZDHHC13 substrates were associated with metabolic disease, including mitochondrial dysfunction. Of these, SDHA, NDUFA7, NDUFV1, NDUFS1, and UQCRC1 are part of the mitochondrial respiratory chain, suggesting that defective S-palmitoylation of such ZDHHC13 substrates may contribute to mitochondrial dysfunction. While the ZDHHC13 enzyme is known to be Golgi-localized, we found that many of its substrates (examined in Fig. 4) are nuclear encoded and localized in the mitochondria. We suspect that many of these proteins are palmitoylated in the cytoplasm by ZDHHC13 prior to being imported into the mitochondria. Besides, a sequence motif recognized by Ankyrin repeat domain of ZDHHC13/17 was identified 41 . This may support that ZDHHC13 recognize its substrates by protein sequence rather than a three dimensional protein structure. Notably, palmitoylation was shown to influence the translocation of some proteins, including IRGM1 and BAX, to mitochondria 42,43 . Moreover, a confirmed mitochondrial palmitoylated protein, CPS1 44 , was also identified as a ZDHHC13 substrate in our study. This suggests that the localization of palmitoylated proteins to mitochondria is a naturally occurring cellular phenomenon. Additionally, high mitochondrial ROS levels in Zdhhc13-deficient hepatocyte and Zdhhc13-knockdown cells may result from inefficient transfer of electrons across the electron transport chain or abnormal removal of ROS. Indeed, 18 candidate ZDHHC13 substrates function involved in ROS production (Table S6). In humans, mitochondrial dysfunction underlies a clinical phenotype of non-thyroidal hypermetabolism called Luft's syndrome 45 . Patients with this syndrome exhibit profuse perspiration and lose weight despite adequate food and water intake; at the cellular level, their mitochondria respire profusely and give off excess energy as heat. The causative gene for Luft's syndrome is still unknown. It would be interesting to investigate whether Luft's syndrome in humans is caused by a defect in one of these ZDHHC13 substrates.
Zdhhc13-deficient mice in this study exhibited mitochondrial dysfunction similar to that in mice carrying a proofreading-deficient mitochondrial DNA polymerase 46 or Mcat-knockout mice 33 . MCAT, a ZDHHC13 substrate identified in this study (Fig. 5), is functionally involved in the de novo synthesis of fatty acyl and lipoyl moieties in mitochondria 47 . The non-palmitoylated mutant form of MCAT (C104G) did not exhibit an altered cellular localization (Fig. S2). This suggests that eliminating the palmitoylation of MCAT may disturb its enzyme activity rather than influence its cellular localization.
Zdhhc13-deficient mice also exhibited lipid abnormalities with low blood triglyceride and cholesterol levels and decreased lipid accumulation in the liver (Fig. 1a,b). Our proteomics results identified several ZDHHC13 candidate substrates that were functionally involved in lipid metabolism, such as FASN, which is involved in lipogenesis, and ACSL1, which is involved in β-oxidation. In these proteins, palmitoylation may provide a membrane anchor for preferred localizations, such as the inner mitochondrial membrane, or may influence their binding affinity with substrates. Since we identified many enzymes involved in lipid metabolism, it is possible that these enzymes form covalent acyl-enzyme intermediates during catalysis. To eliminate this possibility, we confirmed the palmitoylation of three proteins, ACAA2, MCAT, and CTNND1, by MLCC and demonstrated that they are indeed palmitoylated proteins (Fig. S1). In addition, CTNND1, a specific ZDHHC13 substrate confirmed in this work, is a well-known anti-adipogenic factor 48 . It has been reported that cytoplasmic CTNND1 inhibits GLUT-4 membrane trafficking by regulation of RHO and RAC 34 . As the palmitoylated site C618 is located in the RHO binding domain of CTNND1 49 , it is possible that loss of palmitoylation at C618 may disturb the binding of CTNND1 to RHO, thereby inhibiting lipid accumulation. It would be interesting to know whether palmitoylation of CTNND1 plays a role in inhibiting lipid accumulation, and this should be investigated in the future. In summary, this is the first study to demonstrate that defects in ZDHHC13, a PAT, leads to phenotypes in mice, including non-thyroidal hypermetabolism, mitochondrial dysfunction, and abnormal lipid metabolism. Here, we used a label-free quantitative proteomics strategy to comprehensively study the liver S-palmitoylome, allowing identification of candidate ZDHHC13 substrates and their functions. Silencing of Zdhhc13 confirmed that ZDHHC13 is important for mitochondrial function. This study not only revealed the S-palmitoylome for discovery of ZDHHC13 substrates but also shed light on the roles of ZDHHC13 and S-palmitoylation in mitochondrial dysfunction and hypermetabolism. We anticipate that this study will be useful in developing therapeutic strategies for diseases associated with mitochondrial dysfunction and hypermetabolism in the future.

Materials and Methods
Mice and genotyping and cell lines. The Zdhhc13-deficient mice were generated by ENU mutagenesis as described previously 31 . Genotype was analysed by sequencing tail genomic DNA. All the animals and protocols (IACUC number: 11-05-187) used in this study were approved by the Institutional Animal Care and Utilization Committee of Academia Sinica. All methods were performed in accordance with the relevant guidelines and regulations. HEK-293T, Hep1-6, 3T3-L1 and MCF-7 cells were purchased and cultured as described by the American Tissue Culture Collection (ATCC).
Metabolic rate measurement. The metabolic performance was measured using a PhenoMaster System (TSE Lab Master System, Homburg, Germany). Six-week-old mice (n = 8) were placed in the PhenoMaster open-circuit indirect calorimeter at room temperature (22 °C-24 °C). Mice were allowed to acclimatize in the chambers for at least 24 h. Food and water were provided ad libitum and measured by the built-in automated instruments. VO 2 , VCO 2 , and heat production were measured for the following 48 h.
Extraction of membrane proteins from mice liver. Liver tissue samples from 6-week-old mice were homogenized in LB buffer (150 mM NaCl, 50 mM Tris, 5 mM EDTA; pH 7.2) with 10 mM N-ethylmaleimide (NEM, Sigma-Aldrich), 2X protease inhibitor (PI, Roche), and 2 mM PMSF for 2 min at 24,000 rpm in a Polytron PT3000 homogenizer. Tissues were further sonicated using a Branson Sonifier 150 with 10 cycles of 1 s ON and 2 s OFF. The total membrane fraction was collected by high speed centrifugation (200,000 × g, 30 min, 4 °C) and was resuspended in 3 mL LB buffer with 1.7% Triton X-100, 10 mM NEM, 1X PI, and 1 mM PMSF for 1 h at 4 °C with end to end rotation. Particles were removed by low-speed centrifugation (250 × g, 5 min, 4 °C). The membrane pellet was precipitated by chloroform-methanol precipitation and incubated with 4SB buffer (4% SDS, 50 mM NEM, 50 mM Tris, 5 min EDTA; pH 7.2) at 37 °C until it dissolved. The membrane proteins were applied in MS analysis and western blotting.
Alkylating Resin-assisted capture (RAC), on-resin digestion, and peptide elution. Detection of S-palmitoylation was performed using the alkylating RAC method 50 with minor modifications as our previous strategy 40 . In brief, the concentration of NEM-blocked membrane proteins from liver tissues of WT or Zdhhc13-deficient mice was adjusted to 2 mg/mL and reacted with final 5 mM TCEP (Sigma-Aldrich) at 37 °C for 30 min to reduce all disulfide bonds. Then the samples were reacted with final 200 mM iodoacetamide (IAM, Sigma-Aldrich) in blocking buffer (20 mM HEPES buffer, pH 7.7, plus 5% (w/v) SDS) at 37 °C for 2 h in the dark to block all of the reduced cysteine residues. After acetone precipitation to remove excess IAM, the mixture was centrifuged at 16000 × g for 20 min. The protein pellet was then washed three times with 95% ice-cold ethanol and resuspended in 20 mM HEPES (pH 7.2)/1 mM EDTA buffer/1% SDS. The protein concentration was determined using the BCA Protein Assay Kit (Pierce) and adjusted to 2 mg/ml. The S-palmitoyl groups on proteins were reduced with or without final 0.7 M hydroxylamine (HA, pH 7.2, Sigma, St. Louis, MO) at 37 °C for 2 h to reduce the thioester linkages. After removing excess HA by acetone precipitation and washed three times with 95% ice cold ethanol, the protein pellets were resuspended in 20 mM HEPES (pH 7.2)/1 mM EDTA buffer/1% SDS. In particular, equal amounts of proteins (1.5 mg for mass spectrometry and 150 μg for immunoblot experiments) were harvested and the concentration was adjusted to 1-2 mg/mL. For quantitation by MS analysis, 1 μg BSA protein was added into each sample as the internal standard protein. Then, 300 μl (5:1, proteins: sepharose, w/v) of activated thiol sepharose 4B (GE Healthcare) pre-washing three times with degas PBS was reacted with each sample at 30 °C for 1 h with frequent vortexing. After washing three times by PBS and removing the buffer, the capture proteins on sepharoses were resuspended by 25 mM TEABC (Sigma-Aldrich) and digested with 25 μg trypsin at 37 °C for 16 h with frequent vortexing. The supernatant of digested samples was removed and the captured sepharoses were washed three times by PBS. The captured peptides on sepharoses were resuspended in 25 mM TEABC, eluted by final 5 mM TCEP at 37 °C for 30 min, and dried completely under vacuum. The eluted peptides were desalted by C18 Zip-tip TM (Millipore) and subjected to downstream MS analysis. Triplicate biological experiments and LC-MS/MS analyses were performed on each sample to obtain confident identification and quantitation.
TMT10-plex labelling of membrane proteins. Membrane proteins from each sample after reduction and alkylating were precipitated by ice acetone and resuspended by 6 M urea in 50 mM HEPES (pH8.2). Protein extracts (10 μg for each) were then diluted to 1 M urea with 25 mM TEABC, and 0.5 μg of trypsin was added and incubated overnight at 37 °C. After stopped by acidification with 0.1% TFA (v/v), peptides were dried and resuspended by 20 μl of 100 mM TEABC. Ten micrograms of peptides from each sample was labeled with TMT reagent. TMT10-plex reagents (0.8 mg) were dissolved in anhydrous acetonitrile (80 μL), of which 8 μL was added into the peptides making final acetonitrile concentration of approximately 30% (v/v). The labelling reaction proceeded for 1 h at room temperature and then was quenched with 1.6 μl of 5% HA to a final concentration of 0.3% (v/v). The TMT-labeled samples were mixed equally and vacuum centrifuged to near dryness. Finally, 20 μg of TMT10-labeled peptide mixtures were fractionated by high pH reverse-phase stage-tip to 8 fractions and dried.
Mass Spectrometry and data analysis. Eluted palmitoylated peptides mixtures from alkylating RAC and TMT10-labeled peptide mixtures from membrane subfraction were analysed by LTQ-Orbitrap XL and LTQ-Orbitrap Fusion MS, respectively. All MS raw data were processed with Proteome Discover version 2.1 (Thermo Fisher Scientific). For label-free quantification of palmitoylated peptides, data analysis was performed using IDEAL-Q software [34]. Membrane proteomic quantitation was performed by Proteome Discover version 2.1. Full experimental procedures are available in the supplementary data.
Protein annotation. For subcellular localization, molecular function annotations, and information linked to diseases, all identified proteins were analysed using the Ingenuity Pathway Analysis (IPA, http://www.ingenuity.com) Knowledge Base and the Gene Ontology (GO) consortium. The potential protein-protein interactions and networks of the quantified S-palmitoylated proteins were further annotated by literature mining and IPA analysis.

Measurement of O 2 consumption rate.
Cells was seeded onto an XF 96-well microplate at a density of 5000 cells/well in hepatocyte or 10000 cells/well in Hep1-6 cells. Cells were incubated in DMEM for 1 h and in DMEM without bicarbonate for 30 min. The microplate was loaded into a XF-96 Extracellular Flux Analyzer (Seahorse Bioscience) following manufacturer's instructions. Oligomycin, FCCP and Rotenone (Sigma-Aldrich) were injected into the medium sequentially to detect O 2 consumption rate of cells. All experiments were performed at 37 °C.
The informations about Blood Chemistry, Detection of lipid accumulation by Oil red O stain, Liquid Chromatography and Tandem Mass Spectrometry, Mass Spectrometry Raw data analysis, Label-free quantitation of palmitoylated peptides, Western blotting, Zdhhc13-knockdown in Hep1-6 and 3T3-L1 cells, JC1 staining, Measurement of mitochondrial reactive oxygen species (ROS), Metabolic labelling combined with click chemistry and Immunofluorescence assay were addressed in supplementary information.