Proteomic analysis of protein lysine 2-hydroxyisobutyrylation (Khib) in soybean leaves
BMC Plant Biology volume 23, Article number: 23 (2023)
Protein lysine 2-hydroxyisobutyrylation (Khib) is a novel post-translational modification (PTM) discovered in cells or tissues of animals, microorganisms and plants in recent years. Proteome-wide identification of Khib-modified proteins has been performed in several plant species, suggesting that Khib-modified proteins are involved in a variety of biological processes and metabolic pathways. However, the protein Khib modification in soybean, a globally important legume crop that provides the rich source of plant protein and oil, remains unclear.
In this study, the Khib-modified proteins in soybean leaves were identified for the first time using affinity enrichment and high-resolution mass spectrometry-based proteomic techniques, and a systematic bioinformatics analysis of these Khib-modified proteins was performed. Our results showed that a total of 4251 Khib sites in 1532 proteins were identified as overlapping in three replicates (the raw mass spectrometry data are available via ProteomeXchange with the identifier of PXD03650). These Khib-modified proteins are involved in a wide range of cellular processes, particularly enriched in biosynthesis, central carbon metabolism and photosynthesis, and are widely distributed in subcellular locations, mainly in chloroplasts, cytoplasm and nucleus. In addition, a total of 12 sequence motifs were extracted from all identified Khib peptides, and a basic amino acid residue (K), an acidic amino acid residue (E) and three aliphatic amino acid residues with small side chains (G/A/V) were found to be more preferred around the Khib site. Furthermore, 16 highly-connected clusters of Khib proteins were retrieved from the global PPI network, which suggest that Khib modifications tend to occur in proteins associated with specific functional clusters.
These findings suggest that Khib modification is an abundant and conserved PTM in soybean and that this modification may play an important role in regulating physiological processes in soybean leaves. The Khib proteomic data obtained in this study will help to further elucidate the regulatory mechanisms of Khib modification in soybean in the future.
Protein post-translational modification (PTM) is known to be the most effective strategy for altering protein properties, such as stability, activity, cellular localization, and protein–protein interactions [1, 2], and play key roles in a variety of biological processes and metabolic pathways [3, 4]. Currently, more than 400 different types of PTMs have been discovered in various organisms . Among these PTMs, lysine residue of proteins has been found to be one of the most common modification sites, where acylation modification of lysine is a type of modification that has been extensively studied in recent years, including methylation (Kme), acetylation (Kac), crotonylation (Kcr), butyrylation (Kbu), glutarylation (Kglu), succinylation (Ksucc), malonylation (Kmal), propionylation (Kpr), ubiquitination (Kub) and 2-hydroxyisobutyrylation (Khib) [6, 7]. In recent years, due to the development of high-resolution mass spectrometry-based proteomic techniques, proteome-wide identification and analysis of lysine acylation has been performed in a variety of organisms, including eukaryotes and prokaryotes, and these types of lysine acylation modifications have been found to be a family of PTMs that are physiologically relevant [2, 8].
2-Hydroxyisobutyrate is a short-chain fatty acid present in a large number of biofluids in humans, including blood, urine, and feces  and is a precursor for the synthesis of 2-hydroxyisobutyryl coenzyme A . The lysine 2-hydroxyisobutyrylation (Khib) modification is a new protein PTM, which was first identified on histones from mammalian cells by Dai et al. in 2014 . However, recent studies have shown that this modification also occurs in nonhistone proteins in various compartments of cells. To date, Khib modification has been identified in human or animal [9,10,11,12,13], fungus (Aspergillus flavus, Fusarium graminearum, Candida albicans, bacteria (Proteus mirabilis and Toxoplasma gondii) and plants. All these studies indicate that Khib modification is not only widely present in various organisms, but also involved in a wide variety of physiological functions.
In plants, the global identification of Khib-modified proteins has been carried out in several different species, such as rice (Oryza sativa L.) [19, 20], wheat (Triticum aestivum L.) [5, 21, 22], rhubarb (Rheum palmatum L.) , Physcomitrella patens and Arabidopsis siliques. In wheat, 6328 Khib sites in 2186 proteins were identified in three replicates from the roots of cultivar Jimai 44 . However, in the leaves of this wheat cultivar, only 3004 Khib sites in 1104 proteins were identified . While a similar number of Khib sites and Khib-modified proteins were identified in the leaves of wheat cultivar Ak58, which identified 3348 Khib sites in 1074 proteins . In addition, in the leaves of Chinese herb rhubarb, 4333 overlapping Khib sites in 1525 proteins were identified in three independent tests . These results indicate that the number of Khib sites and Khib-modified proteins in leaves is much lower than that in roots. Moreover, the Khib-modified proteins whether identified in wheat roots and leaves, rhubarb leaves, or rice flowers and seedlings are all distributed in various compartments of the cell, with most Khib-modified proteins located in the chloroplast, cytoplasm and nucleus. Meanwhile, the Khib-modified proteins identified in these plants show a variety of functions, such as binding, catalytic activity, and sensitivity of structural molecules [20, 23]. Based on the accumulated research evidence of Khib modification, this modification is considered to be widely present in plant proteins and is an evolutionarily conserved PTM [19, 23].
Soybean [Glycine max (L.) Merr.] is an important legume crop that is widely grown worldwide . As one of the most important economic crops in the world, soybean provides an abundant source of plant protein and oil for human and livestock, accounting for approximately 56% of global oilseed production . The leaf is an important organ of the plant and is considered the energy factory of the plant, which efficiently absorbs energy from sunlight and converts it into bioenergy through photosynthesis, thus, the leaf largely determines the yield of soybean . Therefore, a better understanding of the function of leaf proteins will provide the basis for research to improve soybean yields.
A previous study on the global analysis of lysine acetylation (Kac) modification in soybean leaves has shown that this modification is involved in the regulation of proteins related to ribosome activity, protein biosynthesis, carbohydrate and energy metabolism, photosynthesis and fatty acid metabolism . However, the Khib modification in soybean leaves has not been investigated. To better understand the physiological role of soybean leaves and the potential regulatory function of Khib modification on its physiological activity, we performed the first systematic analysis of protein Khib modification in soybean leaves. Our data suggest that Khib modification is widely present in proteins of soybean leaves and is involved in various biological processes. The Khib sites and proteins identified in this study not only expand our understanding of the functional role of Khib modification, but also provide a platform for further exploration of Khib modification in soybean physiology and biology.
Results and discussion
Global identification of Khib modification in soybean leaves
The Khib modification is a newly discovered modification first found on histone proteins in mammalian cells . Since then, this protein modification has been found in various organisms, including animals, microorganisms and plants. In plants, proteome-wide Khib modifications have been profiled in Arabidopsis siliques, rhubarb , wheat [5, 22], rice [20, 29], etc. Soybean is an important economic crop grown worldwide . Previous studies have reported global analysis of protein Kacmodifiaction in soybean seeds  and leaves  and found this modification participate in various biological processes. However, Khib modification has not been identified in soybean. To gain a comprehensive understanding of the potential role of Khib in the growth, development, and physiological activities of soybean, a systematic analysis of Khib modification in soybean leaves was performed in this study. Total proteins were extracted from soybean leaves at different growth stages and the Khib-modified proteins were then identified by three independent tests, including affinity enrichment using Khib pan antibody and nano-HPLC–MS/MS analysis.
The identification results of the above mentioned three replicates are shown in Fig. 1 and Table S1. In total, 7287, 7439 and 6932 Khib-modified peptides were identified, respectively, of which 4251 peptides overlapped in three replicates (Fig. 1B), with peptide scores above 40 and mass errors within 5 ppm (Fig. 1A), indicating that highly accurate and high-quality MS data were obtained and that the data met the requirements for further analysis. The length distribution of Khib-modified peptides was analyzed and found to vary from 7 to 34 amino acid residues, but most of the peptides were between 7 and 20 residues in length (Fig. 1D).
Furthermore, the Khib-modified peptides identified in the three replicates matched 2195, 2227 and 2152 proteins, respectively, with 1532 proteins overlapping in the three replicates (Fig. 1C). The 1532 Khib-modified proteins contained a range of 1–30 Khib sites per protein, of which 42.1% contained 1 Khib site, 43.5% contained 2–5 Khib sites, 11.1% contained 6–10 Khib sites, and the remaining 3.3% contained more than10 Khib sites (Fig. 1E), indicating that most of the proteins contained multiple Khib sites. Surprisingly, several proteins contained more than 28 Khib sites, such as lipoxygenase (accession number: A0A0R0G554, 28 sites), Glutamine amidotransferase type-2 domain-containing protein (accession number: A0A368UIE9, 30 sites) and an uncharacterized protein (accession number: I1MRV6, 29 sites) (Table S1). The distribution of Khib-modified peptide lengths and the number of Khib sites contained in individual proteins in soybean were similar to those found in wheat [21, 22] and rice . The identification of such a large number of Khib-modified sites and proteins in soybean not only indicates that Khib modification is an abundant and complex PTM, but also indicates that this modification may have a very important role in regulating protein functions in soybean. In a word, we have characterized the lysine 2-hydroxyisobutyrylome in soybean for the first time, and the obtained data on Khib modification will help to elucidate the potential regulatory role of Khib in the physiological functions of soybean.
Functional characterization, subcellular localization and secondary structure analysis of Khib-modified proteins in soybean
Considering that Khib-modified proteins may have important physiological roles in soybean, we performed Gene Ontology (GO) classification analysis to facilitate the understanding of the potential functions of Khib-modified proteins that identified in soybean. GO annotation is commonly used to determine the role of genes or proteins based on three independent ontologies, including biological process, molecular function and cellular component .
The classification results were shown in Fig. 2A and Table S2, from which it is clear that Khib-modified proteins are involved in a variety of biological processes, molecular functions and cellular components. In the classification of “biological process”, most Khib-modified proteins were classified into three major categories: metabolic processes, cellular processes and individual organ processes, accounting for 37.8%, 29.2% and 22.0% of the total number of Khib-modified proteins, respectively (Fig. 2A, cyan bars). In the classification of “molecular function”, most Khib-modified proteins were classified into two categories, namely catalytic activity (43.7%) and binding (36.7%) (Fig. 2A, purple bars), which is in agreement with the classification of biological processes, and both results suggest that enzymatic proteins related to metabolism are more preferred to be Khib modified in soybean. This finding suggests that Khib modification may play an important role in metabolism and cellular process. As for the classification of “cellular component”, 37.3% of Khib proteins were in the cell, 26.2% in the macromolecular complex, 21.4% in the organelle and 14.9% in the membrane (Fig. 2A, orange bars). The subcellular localization of Khib-modified proteins was also predicted using WoLF PSORT software, which showed that most of the proteins were distributed in chloroplasts (46.0%), cytoplasm (28.1%) and nucleus (10.0%) (Fig. 2B and Table S3).
In addition, we performed secondary structure analysis of the identified proteins to further assess the effect of Khib modification on protein function. The results showed that 27.2% of all Khib sites were located in the α-helix, 6.4% in the β-strand, and up to 66.4% in the coil (Fig. 2C and Table S4), indicating that Khib has a preference for secondary structures. However, the distribution trend of Khib sites was not significantly different from the distribution of all lysine sites (Fig. 2C). Moreover, we further assessed the surface accessibility of the identified Khib sites and found that 38.9% of the Khib sites were exposed on the protein surface, slightly lower than all lysine sites (40.5%) (Fig. 2D and Table S4). The lower surface accessibility of the Khib sites indicates that Khib modification may occur in a selective process, which is similar to the percentage of surface accessibility of Khib-modified proteins found in wheat . Furthermore, the distribution trends of the functional classification, subcellular localization and secondary structure of Khib proteins found in soybean are very similar to those of Kac-modified proteins in soybean  and also to those of Khib-modified proteins in other plants, such as wheat [5, 20, 22] and rhubarb . Taken together, these findings suggest that protein Khib modification may be a conserved mechanism for regulating protein function in plants.
Enrichment analysis of Khib-modified proteins in soybean
To further characterize the nature of Khib-modified proteins in soybean, we performed enrichment analysis both of GO annotation and KEGG pathways for the Khib-modified proteins identified in this study. As shown in Fig. 3 and Table S5, the Khib-modified proteins were found to be significantly enriched in both the GO annotation and KEGG pathway.
In the enrichment of GO molecular function, structural constituent of ribosome, structural molecule activity, NAD binding, hydrogen ion transmembrane transporter activity, and ligase activity, coenzyme binding, cofactor binding and oxidoreductase activity, etc., were found to be significantly enriched (Fig. 3A, purple bars, and Fig. S1). The enrichment of these molecule functions indicates that those proteins whose functions are related to ribosome and metabolism may prefer to be Khib modified. The enrichment result of GO biological process showed that glycerol ether metabolic process, glycosyl compound biosynthetic process hydrogen transport, ribose phosphate biosynthetic process, purine-containing compound biosynthetic process were the top 5 highly enriched biological processes (Fig. 3A, cyan bars, and Fig. S1). These five enriched biological processes and the remaining significantly enriched biological processes shown in Fig. 3A are mainly related to metabolism and biosynthesis, which is consistent with the enrichment of GO molecule functions, suggesting that Khib modification is more important for the regulation of biosynthesis and metabolism than on other physiological functions. For the enrichment analysis of GO cellular component, the Khib-modified proteins were significantly enriched in photosynthetic membrane, cytoplasmic part, cytoplasm, and macromolecular complex, etc. (Fig. 3A, orange bars, and Fig. S1).
Consistent with the GO enrichments, a number of pathways related to metabolism and biosynthesis were found to be significantly enriched in KEGG pathway analysis (Fig. 3B, Fig. S2 and Table S6). Among these enriched pathways, citrate cycle (TCA cycle), glyoxylate and dicarboxylate metabolism, oxidative phosphorylation, glycine, serine and threonine metabolism, pyruvate metabolism, glycolysis/gluconeogenesis, pentose phosphate pathway, and fructose and mannose metabolism are involved in central carbon metabolism (the Khib-modified proteins in glycolysis/gluconeogenesis and TCA is shown in Fig. 4A, Khib-modified proteins in pentose phosphate pathway and oxidative phosphorylation shown in Fig. S3A), and pathways of ribosome, selenocompound metabolism, monobactam biosynthesis, glycine, serine and threonine metabolism, lysine biosynthesis, glutathione metabolism, aminoacyl-tRNA biosynthesis, and arginine biosynthesis are involved in biosynthesis (the Khib proteins in ribosome shown in Fig. 4B). These findings also indicated that proteins related to metabolism and biosynthesis tended to be Khib modified. Moreover, several pathways involved in photosynthesis and fatty acid metabolism were also found to be markedly enriched (the Khib proteins in photosynthesis shown in Fig. 4C, and the Khib proteins in carbon fixation in photosynthetic organisms and fatty acid biosynthesis shown in Fig. S3B and C), suggesting that photosynthesis and fatty acid metabolism in soybean may also be regulated by Khib modification. In addition, the pathways associated with proteasome, 2-oxocarboxylic acid metabolism, ascorbate and aldarate metabolism, and porphyrin and chlorophyll metabolism were also found to be significantly enriched (Fig. 3B).
Taken together, the enrichment analysis results suggest that Khib-modified proteins related to biosynthesis and metabolism are extremely important for the physiological functions of soybean leaves. Notably, Khib-modified proteins involved in metabolism and biosynthesis were also found to be significantly enriched in other plants, such as rice [19, 20], wheat [5, 22], and rhubarb . The enrichment of Khib-modified proteins associated with central carbon metabolism and biosynthesis found in plants for which global analyses of Khib modifications have been performed to date suggests that Khib modifications of these proteins may be the most important regulation pattern in plants. Furthermore, Khib-modified proteins enriched in photosynthesis were found not only in soybean in this study, but also in wheat [21, 22] and rice . As we all known that photosynthesis is a fundamental biological process involving the conversion of solar energy into chemical energy, as well as a primary producer of oxygen, which is crucial for plants and even all life on Earth . In fact, photosynthesis-related pathways including photosynthesis-antenna proteins, photosynthesis and carbon fixation in photosynthesis organisms were the top 3 enriched KEGG pathways found in this study (Fig. 3B), indicating that Khib modification is also important for photosynthesis in plants. Therefore, the discovery of Khib modification of photosynthesis-related proteins will help to fully elucidate the mechanism of photosynthesis of leaves.
Analysis of Khib site properties reveals conserved motifs in soybean
To further analyze the natural properties of Khib-modified proteins in soybean, we investigated the sequence motifs around the Khib sites of all identified Khib-modified peptides using the Motif-X software, a tool for extracting overrepresented patterns . In total, 12 conserved motifs were obtained with amino acid residues ranging from -10 to + 10 surrounding the Khib site, such as K*Khib, Khib*K, V*Khib, A*Khib, Khib*A, EKhib, GKhib and Khib*G (“*” indicates one or more random amino acid residues, Fig. 5A and Table S7). Investigation of these motifs found that only five residues are preferred flanking the Khib site, including a basic amino acid residue (K), an acidic amino acid residue (E), and three aliphatic amino acid residues with small side chains (G, A and V). Thus, based on the different properties of these five amino acid residues, these 12 motifs could be divided into three patterns: + 5 to + 8 (except + 6) or -7 position for a basic amino acid residue (K……Khib….K.KK), -1 position for an acidic amino acid residue (EKhib), and -1 to -5 (except -4) or + 3 position for an aliphatic amino acid residue with a small side chain (V.AAAKhib.A/G) (“.” indicates a random amino acid residue at that position).
Consistent with the results of motif analysis, the heat map of the frequency of occurrence of amino acid residues flanking the Khib site showed that A, G, K and V at several positions were overrepresented around the Khib site (Fig. 5B). However, C, M and S were significantly underrepresented at all positions from -10 to + 10 flanking the Khib site. Among these motifs, the most distributed motifs are Khib…..K and Khib….K, matching to 472 (15.3%) and 426 (13.8%) Khib peptides, respectively (Fig. 5C). However, due to the other two motifs flanked by K (K…..Khib and Khib……K) match fewer peptides, the third type of motif pattern, namely Khib flanked by G, A and V, matched the largest number of Khib peptides among the three types of motif patterns mentioned above, accounting for 53.2% of all matched peptides (Fig. 5C).
The motif pattern of Khib site flanking by K at several positions is highly similar to the sequence motifs obtained in other plants, such as wheat  and rhubarb . Additionally, although motifs with K flanking the Khib site were also found in rice, K only occurs at the position of + 9 or -1 . Unlike the motif pattern of Khib site flanking by K, the motif pattern of Khib site flanking by E, namely EKhib, was found in rhubarb , rice  and Arabidopsis siliques, but not in wheat. In the third motif pattern mentioned above, with the exception of GKhib found in rhubarb , the motifs of Khib site flanked by A or G at different positions is only found in Arabidopsis siliques, while the motif of Khib site flanked by V is currently not found in other plants and is unique to soybean.
The several types of motif patterns extracted from plants suggest that there may be multiple enzymes responsible for catalyzing protein Khib modification in plants. Furthermore, the above comparative analysis of motifs in different plants showed both similarities and differences in the motifs extracted in different plants, which seems to suggest that there may be similar enzymes to catalyze Khib modifications of proteins in different plants, but there may also be specific enzymes for each plant.
Analysis of protein interaction networks of Khib-modified proteins in soybean
To better understand the cellular processes regulated by Khib modification in soybean, we performed a protein–protein interaction network for all of the Khib-modified proteins using the STRING database and Cytoscape software . The results showed that a total of 1043 Khib-modified proteins were matched as the network nodes and 5501 interactions were obtained from the STRING database with a combined score ≥ 0.90. The global PPI network graph of these interactions is shown in Fig. S4, and the detailed interactions are shown in Table S8. In addition, 16 clusters of Khib proteins were retrieved from the global PPI network by using a clustering algorithm performed with the MCODE plug-in tool kit (Table S8), and 5 representatives of highly-connected clusters, including ribosome, photosynthesis, oxidative phosphorylation, carbon fixation in photosynthetic organisms and TCA cycle are shown in Fig. 6. These data suggest that Khib modifications tend to occur in proteins associated with specific functional clusters. Notably, these extracted clusters are consistent with both KEGG pathway and GO annotation enrichment analysis.
Prior to the present study, the complete interaction network of Khib-modified proteins in plants was identified only in rice flowers  and seedlings . Ribosomes, photosynthesis, and carbon metabolism were identified as abundant clusters in rice seedlings , and these clusters differ from the highly connected clusters identified in rice flowers which include ribosome, proteases, ATPases, heat shock protein 70 (HSP70), histones, mitochondrial proteins, TCA cycle and glycolytic enzymes . The similar PPI clusters, such as ribosomes, photosynthesis, TCA cycle and oxidation phosphorylation were obtained in this study, suggesting that Khib modification has an important regulatory role in these physiological processes in plants.
In the present study, the first systematic analysis of Khib-modified proteins in soybean, an important oil crop and industrial crop, was performed by using a highly sensitive proteomic method. A total of 4251 overlapping Khib-modified peptides matched to 1532 proteins were identified in three replicates, which showed diverse subcellular compartments distribution. Functional and pathway analysis of the Khib-modified proteins suggested that these proteins are involved in multiple cellular functions and metabolic pathways. Moreover, Khib modifications were found to occur preferentially in proteins associated with ribosome/biosynthesis, central carbon metabolism, photosynthesis and fatty acid biosynthesis, suggesting that Khib modifications play a particularly important regulatory role in these pathways in soybean leaves. Furthermore, three types of motif patterns extracted from 12 Khib sequence motifs showed that A/V/G, E, and K were strongly biased around the Khib site. Our global analysis of Khib-modified proteins in soybean indicates that Khib is an abundant and conserved PTM in plants, and the obtained Khib proteome data will contribute to a better understanding of the potential physiological activity of Khib modifications in soybean. Although the above findings were obtained in this study, quantitative analysis of Khib-modified proteins and further biological experiments are needed to confirm the role of the identified Khib modifications in soybean.
Materials and methods
Plant material and cultivation
The soybean used in this study was cultivar Qihuang 34, and the cultivation of soybean was carried out at the experimental site of Shandong Agricultural University, as described previously . Briefly, soybean seeds were sowed in the field and the type soil is Typic-Hapli-Udic Argosols according to the Chinese Soil Taxonomy . Inverted trifoliate leaves were collected from three soybean plants of similar growth at early blooming flowering (R1), early podding (R3), early bulging (R5), and early maturity (R7), respectively, and kept in liquid nitrogen for 24 h and then stored at -80 °C. The collection of plant specimens has been approved by College of Agronomy, Shandong Agricultural University. Experimental research and field study on soybean in this study has complied with the IUCN Policy Statement on Research Involving Species at Risk of Extinction.
Total protein extraction and trypsin digestion
The soybean leaves collected at different stages were mixed in equal amounts and used as biological replicates. Whole proteins of leaves were extracted by using a phenol isolation method  with a slight modification. Briefly, the 400 mg of leaves were grinded into powder in liquid nitrogen and 10% polyvinylpolypyrrolidone (PVPP) and then suspended in 1.6 mL lysis buffer 1 (1 M sucrose, 0.5 M Tris–HCl pH 8.0, 0.1 M KCl, 50 mM ascorbic acid, 1% NP40, 1% sodium deoxycolate, 10 mM EDTA, 10 mM DTT, and 1% plant protease inhibitor cocktail). The obtained sample solutions were sonicated on ice water and then placed on ice for 30 min. The resulting solutions were mixed with equal volumes of tris-saturated phenol (pH 8.0) and then placed on ice for 10 min, followed by centrifugation at 16,000 g for 10 min at 4 °C and collection of the upper phenol phase. The resulting phenol solutions were mixed with four volumes of precipitation buffer (methanol with 0.1 M ammonium acetate, precooled at -20 ℃ before use), and then placed overnight at -20 °C to precipitate the proteins. The precipitates were collected by centrifugation at 16,000 g for 10 min at 4 °C, and the pellets were rinsed twice with cold methanol. The protein pellets were dissolved in 1.6 mL lysis buffer 2 (8 M urea, 50 mM Tris–HCl pH 8.0, 1% NP40, 1% sodium deoxycolate, 10 mM EDTA, 5 mM DTT, 1% plant protease inhibitor cocktail and 1% phosphatase inhibitor cocktail), followed by sonication on ice water and centrifugation at 20,000 g for 10 min at 4 °C. The final supernatants were collected and the protein concentration was determined using the 2D Quant kit (GE Healthcare).
For tryptic digestion, 2 mg of proteins were first reduced with 5 mM DTT at 56 °C for 30 min, followed by alkylation with 30 mM iodoacetamide (IAM) at room temperature in the dark for 45 min. Then, proteins were precipitated by adding 5 volumes of -20 °C precooled methanol and placed at -20 °C overnight, followed by centrifugation at 20,000 g for 10 min at 4 °C, and the obtained precipitates were rinsed twice with precooled methanol and then placed at -20 °C for 1 h. Finally, the resulting proteins were suspended in 0.1 M triethylammonium bicarbonate (TEAB) buffer and digested overnight with 50 µg of trypsin (Promega). After digestion, the tryptic peptides were desalted using a Strata X C18SPE column (Phenomenex) and vacuum dried before affinity enrichment.
Affinity enrichment of Khib-modified peptides
For affinity enrichment of Khib-modified peptides, the dried peptides were dissolved in NETN buffer (100 mM NaCl, 1 mM EDTA, 50 mM tris–HCl pH 8.0, 0.5% NP-40) and then incubated with anti-2-hydroxyisobutyryllysine pan antibody agarose beads (Micrometer Biotech) at 4 °C overnight with gentle shaking. Then, after washing the agarose beads four times with NETN buffer and once with double distilled water, the Khib-modified peptides were eluted with 0.1% trifluoroacetic acid (TFA) and desalted with a C18 ZipTips column, followed by vacuum drying for further use.
Liquid chromatography-mass spectrometry analysis
The dried Khib-modified peptides were firstly dissolved in buffer A (0.1% formic acid (FA)) and then separated by liquid chromatography (LC) with a reversed-phase C18analytical column (Thermo Acclaim PepMap RSLC C18column, 75 μm × 500 mm, 2 μm particles) at a flow rate of 250 nL/min on an UltiMate RSLCnano 3000 system (Thermo Scientific) as described previously [37, 38]. The gradient was set to 2–10% buffer B (0.1% FA and 80% acetonitrile in water) for 6 min, 10–20% buffer B for 45 min, 20–80% buffer B for 7 min and finally hold at 80% for 4 min.
The eluted peptides were ionized and electrosprayed (2.0 kV voltage) into the mass spectrometer (Thermo Scientific Q Exactive HFX) coupled online to the LC system. The mass spectrometric analysis was performed in data-dependent mode with an automatic switch between a full MS scan and a MS/MS scan in the Orbitrap. The peptides were detected in the MS at a resolution of 60,000 with a scan range of 350–1800 m/z, and with automatic gain control (AGC) of 1E6. The top 15 precursor ions were selected for MS/MS by higher-energy collision dissociation (HCD) using NCE setting as 26%, and the MS/MS spectra were detected at a resolution of 30,000. The dynamic exclusion for the data-dependent scan was 6 s, and AGC was set at 1E5. LC–MS/MS analysis was performed by Micrometer Biotech Company (Hangzhou, China).
The resulting MS/MS raw data were searched using MaxQuant search engine (v.188.8.131.52) against Glycine max (soybean) database from Uniprot (74,863 proteins, https://www.uniprot.org/proteomes/UP000008827). False discovery rate (FDR) thresholds for peptides, proteins and modifications adjusted to < 1%. Trypsin/P was specified as a cleavage enzyme with up to 4 missing cleavages and set the minimum number of amino acids for peptide as 7. Carbamidomethylation on Cys was specified as a fixed modification and oxidation on Met, 2-hydroxyisobutyrylation on Lys were specified as variable modifications. The mass tolerances for precursor ions and fragment ions were set at 10 ppm and 0.02 Da, respectively. Khib sites localization probability was set as > 0.75, minimum score for modified peptides was set > 40 and minimum delta score for modified peptides was set > 8. The identification of modified peptides requires at least 1 MS/MS. The MS data were deposited to ProteomeXchange Consortium via the PRIDE partner repository. The accession number is PXD03650.
Functional annotation of Khib-modified peptides and proteins
To explore the potential roles of the identified Khib-modified peptides and proteins, detailed bioinformatics analyses were performed. Functional annotation of Gene Ontology (GO) was performed for the functional classification and functional enrichment analysis . For functional annotation of Kyoto Encyclopedia of Genes and Genomes (KEGG), the KAAS tool (v.2.0) was used for pathway enrichment analysis . GO term or KEGG pathway enrichment analysis were carried out with the DAVID tool , and annotation terms with a corrected p-value < 0.05 by Fisher’s exact test were considered significantly enriched. The motif-X software  was used to analyse amino acid sequence motifs of Khib peptides. Motif-based clustering analysis was also performed and heatmap visualization of cluster members was performed using the “heatmap.2” function in the “gplots” R package. Subcellular location analysis of identified proteins was conducted using the WoLF PSORT platform (https://wolfpsort.hgc.jp). Secondary structure and surface accessibility analysis were performed using NetSurfP online software (v.2.0). STRING (v.11.0)  was used to evaluate potential protein–protein interaction relationships among those identified proteins, and only confidence score ≥ 0.7 were selected as significant. PPI networks were constructed and visualized using Cytoscape software (v.3.7) , and modules of PPI networks were screened using the Molecular Complex Detection (MCODE) plug-in tool in Cytoscape.
Availability of data and materials
The raw data of MS spectrometry were deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier of PXD036505 (http://www.ebi.ac.uk/pride/archive/projects/PXD036505).
Walsh CT, Garneau-Tsodikova S, Gatto GJ Jr. Protein posttranslational modifications: the chemistry of proteome diversifications. Angew Chem Int Ed Engl. 2005;44:7342–72.
Huang S, Tang D, Dai Y. Metabolic Functions of Lysine 2-Hydroxyisobutyrylation. Cureus. 2020;12:e9651.
Mann M, Jensen ON. Proteomic analysis of post-translational modifications. Nat Biotechnol. 2003;21:255–61.
Beltrao P, Bork P, Krogan NJ, van Noort V. Evolution and functional cross-talk of protein post-translational modifications. Mol Syst Biol. 2013;9:714.
Bo F, Shengdong L, Zongshuai W, Fang C, Zheng W, Chunhua G, Geng L, Ling’an K. Global analysis of lysine 2-hydroxyisobutyrylation in wheat root. Sci Rep. 2021;11:6327.
Pan J, Chen R, Li C, Li W, Ye Z. Global Analysis of Protein Lysine Succinylation Profiles and Their Overlap with Lysine Acetylation in the Marine Bacterium Vibrio parahemolyticus. J Proteome Res. 2015;14:4309–18.
Hong G, Su X, Xu K, Liu B, Wang G, Li J, Wang R, Zhu M, Li G. Salt stress downregulates 2-hydroxybutyrylation in Arabidopsis siliques. J Proteomics. 2022;250:104383.
Huang H, Luo Z, Qi S, Huang J, Xu P, Wang X, Gao L, Li F, Wang J, Zhao W, et al. Landscape of the regulatory elements for lysine 2-hydroxyisobutyrylation pathway. Cell Res. 2018;28:111–25.
Dai L, Peng C, Montellier E, Lu Z, Chen Y, Ishii H, Debernardi A, Buchou T, Rousseaux S, Jin F, et al. Lysine 2-hydroxyisobutyrylation is a widely distributed active histone mark. Nat Chem Biol. 2014;10:365–70.
Huang H, Tang S, Ji M, Tang Z, Shimada M, Liu X, Qi S, Locasale JW, Roeder RG, Zhao Y, et al. p300-Mediated Lysine 2-Hydroxyisobutyrylation Regulates Glycolysis. Mol Cell. 2018;70:663-78.e6.
Ge H, Li B, Chen W, Xu Q, Chen S, Zhang H, Wu J, Zhen Q, Li Y, Yong L, et al. Differential occurrence of lysine 2-hydroxyisobutyrylation in psoriasis skin lesions. J Proteomics. 2019;205:103420.
Xie T, Dong J, Zhou X, Tang D, Li D, Chen J, Chen Y, Xu H, Xue W, Liu D et al. Proteomics analysis of lysine crotonylation and 2-hydroxyisobutyrylation reveals significant features of systemic lupus erythematosus. Clin Rheumatol. 2022.
Du R, Liu G, Huang H. Deep 2-Hydroxyisobutyrylome in mouse liver expands the roles of lysine 2-hydroxyisobutyrylation pathway. Bioorg Med Chem. 2022;57:116634.
Lv Y, Wang J, Yang H, Li N, Farzaneh M, Wei S, Zhai H, Zhang S, Hu Y. Lysine 2-hydroxyisobutyrylation orchestrates cell development and aflatoxin biosynthesis in Aspergillus flavus. Environ Microbiol. 2022.
Zhao Y, Zhang L, Ju C, Zhang X, Huang J. Quantitative multiplexed proteomics analysis reveals reshaping of the lysine 2-hydroxyisobutyrylome in Fusarium graminearum by tebuconazole. BMC Genomics. 2022;23:145.
Zheng H, Song N, Zhou X, Mei H, Li D, Li X, Liu W. Proteome-Wide Analysis of Lysine 2-Hydroxyisobutyrylation in Candida albicans. mSystems. 2021; 6.
Dong H, Guo Z, Feng W, Zhang T, Zhai G, Palusiak A, Rozalski A, Tian S, Bai X, Shen L, et al. Systematic Identification of Lysine 2-hydroxyisobutyrylated Proteins in Proteus mirabilis. Mol Cell Proteomics. 2018;17:482–94.
Yin D, Jiang N, Zhang Y, Wang D, Sang X, Feng Y, Chen R, Wang X, Yang N, Chen Q. Global Lysine Crotonylation and 2-Hydroxyisobutyrylation in Phenotypically Different Toxoplasma gondii Parasites. Mol Cell Proteomics. 2019;18:2207–24.
Meng X, Xing S, Perez LM, Peng X, Zhao Q, Redoña ED, Wang C, Peng Z. Proteome-wide Analysis of Lysine 2-hydroxyisobutyrylation in Developing Rice (Oryza sativa) Seeds. Sci Rep. 2017;7:17486.
Xue C, Qiao Z, Chen X, Cao P, Liu K, Liu S, Ye L, Gong Z. Proteome-Wide Analyses Reveal the Diverse Functions of Lysine 2-Hydroxyisobutyrylation in Oryza sativa. Rice (N Y). 2020;13:34.
Zhang N, Zhang L, Li L, Geng J, Zhao L, Ren Y, Dong Z, Chen F. Global Profiling of 2-hydroxyisobutyrylome in Common Wheat. Genomics Proteomics Bioinformatics. 2021.
Feng B, Li S, Wang Z, Cao F, Wang Z, Li G, Liu K. Systematic analysis of lysine 2-hydroxyisobutyrylation posttranslational modification in wheat leaves. PLoS ONE. 2021;16:e0253325.
Qi T, Li J, Wang H, Han X, Li J, Du J. Global analysis of protein lysine 2-hydroxyisobutyrylation (K(hib)) profiles in Chinese herb rhubarb (Dahuang). BMC Genomics. 2021;22:542.
Yu Z, Ni J, Sheng W, Wang Z, Wu Y. Proteome-wide identification of lysine 2-hydroxyisobutyrylation reveals conserved and novel histone modifications in Physcomitrella patens. Sci Rep. 2017;7:15553.
Zhou X, Wang D, Mao Y, Zhou Y, Zhao L, Zhang C, Liu Y, Chen J. The Organ Size and Morphological Change During the Domestication Process of Soybean. Front Plant Sci. 2022;13:913238.
Wilson RFJSNY. Soybean: Market Driven Research Needs. 2008.
Gonzalez N, Vanhaeren H, Inzé D. Leaf size control: complex coordination of cell division and expansion. Trends Plant Sci. 2012;17:332–40.
Li G, Zheng B, Zhao W, Ren T, Zhang X, Ning T, Liu P. Global analysis of lysine acetylation in soybean leaves. Sci Rep. 2021;11:17858.
Chen X, Xu Q, Duan Y, Liu H, Chen X, Huang J, Luo C, Zhou DX, Zheng L. Ustilaginoidea virens modulates lysine 2-hydroxyisobutyrylation in rice flowers during infection. J Integr Plant Biol. 2021;63:1801–14.
Smith-Hammond CL, Swatek KN, Johnston ML, Thelen JJ, Miernyk JA. Initial description of the developing soybean seed protein Lys-N(ε)-acetylome. J Proteomics. 2014;96:56–66.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium Nat Genet. 2000;25:25–9.
Suga M, Shen J-R. Structural variations of photosystem I-antenna supercomplex in response to adaptations to different light environments. Curr Opin Struct Biol. 2020;63:10–7.
Chou MF, Schwartz D. Biological sequence motif discovery using motif-x. Curr Protoc Bioinformatics. 2011; Chapter 13:Unit 13.5–24.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27:431–2.
Li Z, Liu Z, Zhang M, Li C, Li YC, Wan Y, Martin CG. Long-term effects of controlled-release potassium chloride on soil available potassium, nutrient absorption and yield of maize plants. Soil and Tillage Research. 2020;196:104438.
Chitteti BR, Tan F, Mujahid H, Magee BG, Bridges SM, Peng Z. Comparative analysis of proteome differential regulation during cell dedifferentiation in Arabidopsis. Proteomics. 2008;8:4303–16.
Ye C, Ge Y, Zhang Y, Zhou L, Chen W, Zhu X, Pan J. Deletion of vp0057, a Gene Encoding a Ser/Thr Protein Kinase, Impacts the Proteome and Promotes Iron Uptake and Competitive Advantage in Vibrio parahaemolyticus. J Proteome Res. 2021;20:250–60.
Zhu X, Feng C, Zhou L, Li Z, Zhang Y, Pan J. Impacts of Ser/Thr Protein Kinase Stk1 on the Proteome, Twitching Motility, and Competitive Advantage in Pseudomonas aeruginosa. Front Microbiol. 2021;12:738690.
Dimmer EC, Huntley RP, Alam-Faruque Y, Sawford T, O’Donovan C, Martin MJ, Bely B, Browne P, Mun Chan W, Eberhardt R, et al. The UniProt-GO Annotation database in 2011. Nucleic Acids Res. 2012;40:D565–70.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4:P3.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
We thank Micrometer Biotech Company (Hangzhou, China) for performing the LC-MS/MS analysis.
This study was financially supported by the National Natural Science Foundation of China (31401339), the Shandong Agricultural Science and Technology Fund (2019YQ014), and funds from the Shandong “Double Tops” Program.
Ethics approval and consent to participate
All methods involving soybean samples were carried out in accordance with relevant guidelines of “List of National Key Protected Wild Plants of China”.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Fig. S1. The bubble chart shows the GO enrichment analysesof Khib-modifiedproteins in soybean leaves. Fig. S2. Thebubble chart shows the enrichment of KEGG pathway of Khib-modifiedproteins in soybean leaves. Fig. S3. Representative pathway mapsof Khib-modified proteins. (A) Khib-modified proteins in central carbon metabolism, includingpentose phosphate pathway and oxidative phosphorylation. (B) Khib-modifiedproteins in carbon fixation in photosynthetic organisms. (C) Khib-modifiedproteins in fatty acid biosynthesis. The identified Khib-modifiedproteins are indicated in boxes with a red background. Fig. S4. The global PPI network ofidentified Khib-modified proteins in soybean leaves. Table S1. All identified Khib sites and proteins that overlap in three replicates in soybean leaves. Table S2. GO functional classification of Khib-modified proteins in soybean leaves. Table S3. Subcellular localization distribution of Khib-modified proteins in soybean leaves. Table S4. Secondary structure and surface accessibility distribution of Khib-modified peptides in soybean leaves. Table S5. GO enrichment analysis of Khib-modified proteins in soybean leaves. Table S6. KEGG pathway enrichment of Khib-modified proteins in soybean leaves. Table S7. Motif analysis of Khib-modified peptides in soybean leaves. Table S8. The PPI network of Khib-modified proteins in soybean leaves.
About this article
Cite this article
Zhao, W., Ren, TH., Zhou, YZ. et al. Proteomic analysis of protein lysine 2-hydroxyisobutyrylation (Khib) in soybean leaves. BMC Plant Biol 23, 23 (2023). https://doi.org/10.1186/s12870-022-04033-6