Secondary Logo

Journal Logo

Comprehensive analysis of lysine crotonylation in proteome of maintenance hemodialysis patients

Chen, Wenbiao PhDa; Tang, Donge PhDb; Xu, Yong PhDb; Zou, Yaoshuang PhDc; Sui, Weiguo PhDc; Dai, Yong PhDb,*; Diao, Hongyan PhDa,*

Section Editor(s): Bahous., Sola Aoun

doi: 10.1097/MD.0000000000012035
Research Article: Clinical Trial/Experimental Study
Open
SDC

Introduction: Histone post-translational modifications (PTMs) carry epigenetic information to regulate diverse cellular processes at the chromatin level. Crotonylation, one of the most important and common PTMs, plays a key role in the regulation of various biological processes. However, no study has evaluated the role of lysine crotonylation in maintenance hemodialysis patients (MHP).

Methods: Here, we comparatively evaluated the crotonylation proteome of normal controls (NC) and MHP using liquid chromatography tandem mass spectrometry (LC-MS/MS) coupled with highly sensitive immune-affinity purification.

Results: A total of 1109 lysine modification sites distributed on 347 proteins were identified, including 93 and 252 crotonylated upregulated and downregulated proteins, respectively. Thus, a decrease in crotonylation of histone proteins was observed in patients with kidney failure undergoing maintenance hemodialysis. Intensive bioinformatic analysis revealed that most of the crotonylated proteins were distributed in the cytoplasm, nucleus, mitochondria, and extracellular region. Gene ontology enrichment analysis showed that the crotonylated proteins were significantly enriched in the platelet alpha granule lumen, platelet degranulation, and cell adhesion molecule binding. In addition, protein domain, including fibrinogen alpha/beta/gamma chain, zinc finger, and WD40-repeat-containing domain, were significantly enriched in crotonylated proteins. Kyoto Encyclopedia of Genes and Genomes (KEGG)-based functional enrichment analysis revealed that crotonylated proteins were enriched in complement and coagulation cascades, cardiac muscle contraction, and hematopoietic cell lineage, all of which have important associations with hemodialysis complications.

Conclusions: This is the first report on the global crotonylation proteome of MHP. Lysine crotonylation was found to play important regulatory roles in pathophysiological processes in MHP.

aState Key Laboratory for Diagnosis and Treatment of Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou

bClinical Medical Research Center, the Second Clinical Medical College of Jinan University (Shenzhen People's Hospital), Shenzhen, Guangdong

cNephrology Department of Guilin No.181 Hospital, Guangxi Key Laboratory of Metabolic Diseases Research, Guilin Key laboratory of Kidney Diseases Research, Guilin, Guangxi, P.R. China.

Corresponding author: Yong Dai, Clinical Medical Research Center, The Second Clinical Medical College, Jinan University, Shenzhen People's Hospital, Shenzhen, 518020, P.R. China (e-mail: daiyong22@aliyun.com); Hongyan Diao, State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, 31003, P.R. China (e-mail: diaohy@zju.edu.cn).

Abbreviations: CKD = chronic kidney disease, DPF = double PHD finger, GO = gene ontology, H3K4Cr = lysine-4 crotonylated histone H3, HD = hemodialysis, HDACs = class I histone deacetylases, Kcr = Lysine crotonylation, KEGG = Kyoto Encyclopedia of Genes and Genomes, KF = kidney failure, LC-MS/MS = liquid chromatography tandem mass spectrometry, MHP = maintenance hemodialysis patients, MS = mass spectrometry, NC = normal controls, PBMCs = peripheral blood mononuclear cells, PTMs = histone post-translational modifications, SAHA = suberoylanilide hydroxamic acid, TEAB = triethylammonium bicarbonate.

This project was supported by the National Natural Science Foundation of Guangxi (No. 2017JJA130240), Guangxi Key laboratory of Metabolic Diseases Research (No.2017-181h-03), science and technology Planning Project of guilin (No. 20170117-1).

Availability of data: Out of the need of scientific research, reader can contact corresponding authors with email.

There is no competing financial interest.

Supplemental Digital Content is available for this article.

Supplemental digital content is available for this article. Direct URL citations appear in the printed text and are provided in the HTML and PDF versions of this article on the journal's Website (www.md-journal.com).

This is an open access article distributed under the Creative Commons Attribution License 4.0 (CCBY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. http://creativecommons.org/licenses/by/4.0.

Received March 22, 2018

Accepted August 1, 2018

Back to Top | Article Outline

1 Introduction

Protein post-translational modifications (PTMs) play essential roles in various cellular processes that modulates the physical and chemical properties, folding, conformation, stability, and activity of proteins, all of which may modify protein functions.[1,2] The improved throughput of mass spectrometry (MS) technology has not only brought about a surge in proteome-scale studies but also identified several PTMs.[3] Among all PTMs, lysine acetylation is one of the most studied PTMs. Protein lysine acetylation is a ubiquitous and reversible PTM first reported in histones that mainly regulated gene transcription and expression by the catalysis of histone acetyltransferases and histone deacetylases in the nucleus.[4] The discovery of lysine acetylation on tubulin and mitochondrial proteins suggested an important role for lysine acetylation in cellular biology in addition to chromatin biology.[5] Aside from lysine acetylation, some new types of PTMs, such as malonylation and lysine succinylation are identified as important regulatory mechanisms for a broad array of eukaryotic and prokaryotic physiological functions.[6,7]

Lysine crotonylation (Kcr) is a new protein PTM that was initially detected by Zhao's group in histones of human cell lines and mouse sperm.[8] The discovery of Kcr has attracted attention of the scientific community and several researches have been directed to study Kcr. Kcr was first demonstrated as a robust indicator of active promoters and thought to be an important signal in the control of men germ cell differentiation.[8] Sabari et al[9] showed that histone acetyltransferase p300 possesses histone crotonyltransferase activity and demonstrated that p300-catalyzed histone crotonylation directly stimulates transcription to a level greater than that mediated by p300-catalyzed histone acetylation. Bao et al[10] profiled “eraser” enzymes that recognize the lysine-4 crotonylated histone H3 (H3K4Cr) and found that Sirt1, Sirt2, and Sirt3 may catalyze the hydrolysis of lysine-crotonylated histone peptides and proteins. Moreover, Wei et al[11] demonstrated that class I histone deacetylases (HDACs) are the major histone decrotonylases in mammals. Suberoylanilide hydroxamic acid (SAHA), a well-known HDAC, may inhibit decrotonylation on both histones and nonhistone proteins through the inhibition of HDACs.[12] Li et al[13] defined the evolutionarily conserved YEATS domain as a family of Kcr-favorable readers. Furthermore, Xiong et al[14] demonstrated that histone acetylation-binding double PHD finger (DPF) domains of MOZ (KAT6A) and DPF2 (BAF45d) are engaged in a wide range of acylations at H3K14, with the highest affinity for crotonylation, which may act as a new regulatory mechanism of transcription control centered on the readout of histone Kcr. This discovery opened a new chapter in the regulation of Kcr. The study by Liu et al[15] uncovered a biochemical pathway in the regulation of histone Kcr and implicated chromodomain Y-like protein-regulated histone Kcr in spermatogenesis. This study improved the understanding of the physiology of male reproduction and the mechanism underlying spermatogenic failure in azoospermia Factor c-deleted infertile men. Wei et al[16] suggested that protein crotonylation may lead to the inhibition of DNA replication and, hence, influence cell cycle. Taken together, Kcr controls the interpretation of genetic information at the chromatin level and plays a key role in processes such as gene expression and cell fate.

Histone acetylation has been particularly well characterized and the purification and identification of histone acetyltransferases have improved our understanding of the role of targeted lysine acetylation in gene regulation.[17,18] The most important role of histone acetylation is its involvement in gene regulation, and the aberrant regulation of histone acetylation is correlated with major human diseases.[19] Until now, studies on histone acetylation have primarily focused on several diseases, including tumor,[20] neuropsychiatric disorders,[21] lupus,[22] cardiovascular diseases,[23] acute lymphoblastic leukemia,[24] diabetic nephropathy,[25] acute kidney injury, and chronic kidney disease (CKD).[26] “The “writer” and “eraser” for this new histone modification have been reported based on the structural similarity between histone acetylation and histone crotonylation.”[8,9] Furthermore, the same set of enzyme systems were shared by histone acetylation and histone crotonylation.[9,11] Histone crotonylation is physiologically significant and functionally distinct from or redundant to histone acetylation.[11] Hence, we speculate that Kcr exhibits a crucial role in a wide range of biological processes and may be critically implicated in the pathogenesis of diseases. However, very little is known on the involvement of Kcr in diseases and the underlying mechanism.

The prevalence of CKD and kidney failure (KF) is steadily increasing and becoming an important public health problem worldwide.[27] The World Health Organization (WHO) report indicates that kidney diseases were responsible for 2,993,000 years of life lost and 38,104,000 disability-adjusted life years lost globally. Kidney transplantation is the best form of treatment. However, the rate of kidney transplant is very low, owing to the shortage of organ donors. Hemodialysis (HD) and peritoneal dialysis are the most common therapeutic options chosen by the patients when kidney transplantation is limited or contraindicated.[28,29] It was verified that PTMs may form an epigenetic layer over the genetic layer that can respond to environmental cues and external stimuli to alter the expression of genes associated with CKD.[26] The importance of epigenetic alterations, including PTMs, in fibrosis, inflammation, and immunity associated with various renal disorders and in kidney development is well appreciated.[30] KF patients undergoing HD exhibit evident physiopathological changes in the internal environment, which could contributed to water and electrolyte balance disorder, resistance or hyperactivity of small molecules in vivo, organic lesions of various organs, and a decline in the immune system protection. A series of clinical symptoms could followed, such as hypertension or hypotension, hypokalemia or hypocalcemia, hyposecretion erythropoietin and hyperparathyroid parathormone, vascular fibrosis, and pigmentation of skin and so on. And then, HD complications, such as hemorrhage, uremic cardiomyopathy, stroke, anemia, poor nutrition, infection, and pneumonedema could happen and brought severe impairment in the quality of life and higher mortality to HD patients.[31] Therefore, we speculate that Kcr has close relationship with the pathogenesis of KF and pathophysiological changes in patient undergoing maintenance HD.

Here, we investigated the global Kcr of proteome from normal control (NC) and maintenance hemodialysis patients (MHP) using high-resolution liquid chromatography tandem mass spectrometry (LC-MS/MS) coupled with highly sensitive immune-affinity purification. In total, we identified 1109 Kcr sites in 347 proteins. The identified crotonylated proteins were primarily localized in the cytoplasm, nucleus, mitochondria, and extracellular region. Bioinformatic analysis was performed to reveal the biological functions of crotonylated proteins; these proteins were significantly enriched in the platelet alpha granule lumen, platelet degranulation, and cell adhesion molecule binding and participated in complement and coagulation cascades, cardiac muscle contraction, and hematopoietic cell lineage pathways. To our knowledge, this is the first study to describe Kcr in the global proteome of MHP and expand our current understanding of Kcr roles in pathophysiological processes of KF patient undergoing maintenance HD.

Back to Top | Article Outline

2 Materials and methods

2.1 Samples collection

Peripheral blood used to validate candidate samples was collected from 16 patients with KF on maintenance HD and 16 NC. Patients were definitely diagnosed with KF and underwent HD for >6 months in Shenzhen People's Hospital. All 16 patients (age between 35 and 50) had chronic glomerulonephritis as baseline kidney disease and were receiving hemodialysis treatment 3 times a week. Moreover, all 16 patients had no obvious complications of hemodialysis for the time being. The peripheral blood of 16 patients were collected before the first hemodialysis in a week. The NC group (age between 30 and 50 years) corresponded to healthy subjects who presented to the hospital for routine health check-up and who had no kidney disease. All participants were free of clinical signs of any diseases, including cancer, allergic diseases, immune-compromised conditions, diabetes, and other infectious diseases. All patients and healthy controls provided written informed consent. The study was approved by the Clinical Research Ethics Committee of the Shenzhen People's Hospital. After sample collection, peripheral blood mononuclear cells (PBMCs) were isolated by density gradient centrifugation using Ficoll-Hypaque (GE Healthcare Bio-sciences AB, Uppsala, Sweden) according to the manufacturer's protocol. PBMCs were lysed with TRIzol reagent (Invitrogen, Carlsbad, CA) and stored at −80 °C.

Back to Top | Article Outline

2.2 Protein extraction

Sample was sonicated thrice on ice using a high intensity ultrasonic processor (Scientz) in a lysis buffer (8 M urea, 1% protease inhibitor cocktail). Cell debris were removed by centrifugation at 12,000 × g and 4 °C for 10 minutes. The supernatant was collected and protein concentration was determined with the bicinchoninic acid kit (Sigma Chemical Co., St. Louis, MO) according to the manufacturer's instructions.

Back to Top | Article Outline

2.3 Trypsin digestion

For digestion, the protein solution was reduced with 5 mM dithiothreitol for 30 minutes at 56 °C and alkylated with 11 mM iodoacetamide for 15 minutes at room temperature in the dark. The protein sample was diluted by adding 100 mM triethylammonium bicarbonate (TEAB) to urea concentration <2 M. The sample was overnight treated with trypsin at 1:50 trypsin-to-protein mass ratio for the first digestion, followed by the second digestion at 1:100 trypsin-to-protein mass ratio for 4 hours.

Back to Top | Article Outline

2.4 Tandem mass tag/isobaric tags for relative and absolute quantitation (TMT/iTRAQ) labeling

After trypsin digestion, the peptide was desalted with Strata X C18 SPE column (Phenomenex, Torrance, CA) and subjected to vacuum drying. The peptide was reconstituted in 0.5 M TEAB and processed according to the manufacturer's protocol for TMT kit/iTRAQ kit (AB Sciex, Foster City, CA). Briefly, 1 unit of TMT/iTRAQ reagent was thawed and reconstituted in acetonitrile. The peptide mixtures were incubated for 2 hours at room temperature and pooled, desalted, and dried by vacuum centrifugation.

Back to Top | Article Outline

2.5 High-performance liquid chromatography (HPLC) fractionation

The tryptic peptides were fractionated into fractions with high pH reverse-phase HPLC using Thermo Betasil C18 column (5 μm particles, 10 mm ID, 250 mm length). Briefly, peptides were first separated with a gradient of 8% to 32% acetonitrile (pH 9.0) over 60 minutes into 60 fractions. The peptides were then combined into 6 fractions and dried by vacuum centrifugation.

Back to Top | Article Outline

2.6 Affinity enrichment

To enrich Kcr-modified peptides, tryptic peptides dissolved in NETN buffer (100 mM sodium chloride [NaCl], 1 mM ethylenediaminetetraacetic acid [EDTA], 50 mM Tris–HCl, and 0.5% NP-40, pH 8.0) were incubated with pre-washed antibody beads (PTM Bio) at 4 °C overnight with gentle shaking. The beads were washed 4 times with NETN buffer and twice with water. The bound peptides were eluted from the beads with 0.1% trifluoroacetic acid. The eluted fractions were combined and vacuum dried. For LC-MS/MS analysis, the resulting peptides were desalted with C18 ZipTips (Millipore) according to the manufacturer's instructions.

Back to Top | Article Outline

2.7 LC-MS/MS analysis

The tryptic peptides were dissolved in 0.1% formic acid (solvent A) and directly loaded onto a home-made reversed-phase analytical column (15 cm length, 75 μm i.d.). The gradient comprised an increase from 6% to 23% of solvent B (0.1% formic acid in 98% acetonitrile) over 26 minute, 23% to 35% in 8 minute, and increasing to 80% in 3 minute, followed by a hold at 80% for the last 3 minute. All steps were performed at a constant flow rate of 400 nL/min on an EASY-nLC 1000 ultra-performance liquid chromatography (UPLC) system. The peptides were subjected to NSI source, followed by MS/MS in Q ExactiveTM Plus (Thermo) coupled online to UPLC. The electrospray voltage applied was 2.0 kV. The m/z scan range was 350 to 1800 for full scan and intact peptides were detected with an Orbitrap at a resolution of 70,000. Peptides were selected for MS/MS using NCE setting as 28 and the fragments were detected in the Orbitrap at a resolution of 17,500. A data-dependent procedure that alternated between one MS scan followed by 20 MS/MS scans with 15.0 second dynamic exclusion duration was used. Automatic gain control (AGC) was set at 5E4 and the fixed first mass was set as 100 m/z.

Back to Top | Article Outline

2.8 Database search

The resulting MS/MS data were processed using MaxQuant search engine (v.1.5.2.8). Tandem mass spectra were searched against the human database concatenated with a reverse decoy database. Trypsin/P was specified as the cleavage enzyme, allowing up to 4 missing cleavages. The mass tolerance for precursor ions was set as 20 ppm in the first search and 5 ppm in the main search, while that for the fragment ions was set as 0.02 Da. Carbamidomethylation on Cys was specified as fixed modification and oxidation on Met was specified as variable modification. False discovery rate (FDR) was adjusted to <1% and the minimum score for modified peptides was set as >40.

Back to Top | Article Outline

2.9 Gene ontology (GO) annotation

Gene ontology annotation was derived from the UniProt-GOA database (http://www.ebi.ac.uk/GOA/). The identified protein ID was converted to UniProt ID, followed by its mapping to GO IDs using the protein ID. If some identified proteins were not annotated by UniProt-GOA database, InterProScan software (http://www.ebi.ac.uk/interpro/search/sequence-search) was used to annotate protein's GO functions, based on protein sequence alignment method. The proteins were GO annotated based on the following 3 categories: biological process, cellular component, and molecular function. For each category, a two-tailed Fisher exact test was employed to test the enrichment of the differentially expressed protein against all identified proteins. GO with a corrected value of P < .05 was considered significant.

Back to Top | Article Outline

2.10 Domain annotation

The domain functional description of the identified proteins was annotated by InterProScan (a sequence analysis application) based on the protein sequence alignment method using InterPro (http://www.ebi.ac.uk/interpro/) domain database. For each category of proteins, InterPro database was searched and a two-tailed Fisher exact test was used to test the enrichment of the differentially expressed protein against all identified proteins. Protein domains with a corrected value of P < .05 were considered significant.

Back to Top | Article Outline

2.11 Kyoto encyclopedia of genes and genomes (KEGG) pathway annotation

We used KEGG online service tool KAAS (http://www.genome.jp/tools/kaas/) to annotate protein's KEGG database description. The mapping of the annotation result on the KEGG pathway database was performed using KEGG online service tool KEGG mapper. The pathway with a corrected value of P < .05 was considered significant. These pathways were classified into hierarchical categories according to KEGG website.

Back to Top | Article Outline

2.12 Subcellular localization

We used WoLF PSORT (http://www.genscript.com/wolf-psort.html), a subcellular localization predication software to predict subcellular localization.

Back to Top | Article Outline

2.13 Motif analysis

Soft motif-x (http://scan-x.med.harvard.edu/scan-x.html) was used to analyze the model of sequences comprising amino acids at specific positions of modify-21-mers (10 amino acids upstream and downstream of the site) in all protein sequences. All database protein sequences were used as the background database parameter, while other parameters were set as default.

Back to Top | Article Outline

3 Result

3.1 Comparative analysis of whole proteome and Kcr in MHP and NC

The complete experimental procedure included 9 steps (Fig. 1A). Altogether, 1209 Kcr sites distributed in 377 proteins were identified and 1109 sites from 347 proteins were quantifiable (Table S1, http://links.lww.com/MD/C469). A fold change of >1.2 was considered as an upregulation, while that <1/1.2 was termed as downregulation. Differential expression in MHP and NC was observed at 860 sites for 345 proteins; these comprised 155 sites for 93 proteins that were upregulated and 705 sites for 252 proteins that were downregulated. All data are listed in Table S2, http://links.lww.com/MD/C470. To confirm the validation of MS data, the mass error of all identified peptides was assessed. The mass error was centered around 0 and below 10 ppm, suggesting that the mass accuracy of MS data met the requirement (Fig. 1B). Among 347 crotonylated proteins, most proteins contained 1 or 2 crotonylation sites and few proteins had ≥7 or more crotonylation sites (Fig. 1C). Most peptides ranged from 8 to 20 amino acids in length, consistent with the rule of trypsin digestion (Fig. 1D).

Figure 1

Figure 1

Back to Top | Article Outline

3.2 Analysis of crotonylation site motif

To investigate the sequence commonalities surrounding crotonylation sites and compare them with acetylation sites, sequence motifs in all identified crotonylated peptides were investigated using the motif-X program. A total of 7 conserved motifs (KK, KD, AK, EK, K.D, KE, and K.......K) were retrieved (Table S3, http://links.lww.com/MD/C471, Fig. 2A). In particular, motifs KE and KD were strikingly conserved. The significantly conserved amino acids in these motifs, namely E and D, were both negatively charged and rarely identified in other PTMs. These motifs are likely to represent a feature of crotonylation in MHP. In addition, hierarchical cluster analysis was performed to further analyze these motifs (Fig. 2B). The enrichment of positively charged K residues was observed at −10 to −5 and +10 to +5 positions, while negatively charged residues D and E were markedly enriched at −1 to +4 position. Short aliphatic A residues were frequently observed at −7 to +8 position, while the sulfur-containing C residue was undetected.

Figure 2

Figure 2

Back to Top | Article Outline

3.3 Functional classification of Kcr in GO and Kcr subcellular localization

The subcellular localization of the identified Kcr was characterized (Fig. 3A), and the majority of upregulated crotonylated proteins was distributed in the cytoplasm (46%), nucleus (14%), mitochondria (12%), and extracellular region (11%). In contrast, the majority of the downregulated proteins was distributed in the cytoplasm (57%), extracellular region (11%), mitochondria (9%), and nucleus (8%) (Fig. 3B). According to the data, a significant difference was not observed in the localization of the upregulated and downregulated crotonylated proteins.

Figure 3

Figure 3

To obtain an overview of crotonylated proteins in MHP, GO functional classification was evaluated for all crotonylated proteins based on their biological processes, molecular functions, and cellular components (Table S4, http://links.lww.com/MD/C472). Within the biological processes category, most of the upregulated crotonylated proteins were related to cellular process, single-organism process, biological regulation, and response to stimulus (Fig. 3C). In contrast, most of the downregulated proteins were related to cellular process, single-organism process, biological regulation, and response to stimulus (Fig. 3D). Within the cellular component category, most of the upregulated crotonylated proteins were related to cell, organelle, extracellular region, and membrane (Fig. 3E), while most of the downregulated proteins were related to cell, organelle, extracellular region, and membrane (Fig. 3F). Within the molecular function category, most of the upregulated crotonylated proteins were related to binding, catalytic activity, structural molecule activity, and molecular function regulator (Fig. 3G), while most of the downregulated proteins were related to binding, catalytic activity, molecular function regulator, and structural molecule activity (Fig. 3H). GO functional classification revealed no significant difference between the upregulated and downregulated crotonylated proteins and suggested that Kcr may have extensive biological functions.

Back to Top | Article Outline

3.4 Functional enrichment of Kcr in GO, KEGG, and protein domain

The result of GO-based functional enrichment analysis is shown in Table S5, http://links.lww.com/MD/C473. The upregulated crotonylated proteins were highly enriched in the cellular component of platelet alpha granule, secretory granule, and secretory vesicle and molecular functions, including cell adhesion molecule binding and actin filament binding, and biological processes, including platelet degranulation, regulated exocytosis, and exocytosis. However, no downregulated crotonylated proteins could be enriched in GO analysis (Fig. 4A). KEGG-based functional enrichment analysis (Table S6, http://links.lww.com/MD/C474) showed that the upregulated crotonylated proteins were only enriched in 3 important processes, including complement and coagulation cascades, cardiac muscle contraction, and hematopoietic cell lineage, all of which had important relationship with HD complications (Fig. 4B). Crotonylated protein involved in complement and coagulation cascades were retrieved, comprising a dense protein interaction network (Fig. 4C). However, no downregulated crotonylated proteins could be enriched in KEGG analysis. Furthermore, upregulated crotonylated protein domains such as fibrinogen, alpha/beta/gamma chain, C-terminal globular, subdomain 2; zinc finger, LIM-type; fibrinogen, alpha/beta/gamma chain, C-terminal globular domain; fibrinogen, alpha/beta/gamma chain, C-terminal globular, subdomain 1; fibrinogen, alpha/beta/gamma chain, coiled coil domain; and WD40-repeat-containing domain were enriched (Table S7, http://links.lww.com/MD/C475, Fig. 4D), suggestive of the important role of Kcr in these processes. No downregulated crotonylated proteins were observed in the protein domain.

Figure 4

Figure 4

Back to Top | Article Outline

3.5 Cluster analyses in GO, KEGG, and protein domain

To evaluate further details of the functions of Kcr, cluster analyses based on GO, KEGG, and protein domain enrichment were performed. All quantified crotonylated proteins were divided into the following 4 quantiles (Q1–Q4) according to fold changes of Kcr sites: Q1 (0 < ratio < 0.77), Q2 (0.77 < ratio < 0.77), Q3 (1.2 < ratio < 1.3), and Q4 (ratio > 1.3) and P value <.05. We identified 381, 80, 57, and 252 crotonylated proteins in Q1, Q2, Q3, and Q4, respectively. The quantifiable proteins from the 4 categories were plotted for the enrichment-based cluster analysis. Q1 and Q2 were regarded as downregulated proteins, while Q3 and Q4 were regarded as upregulated proteins (Fig. 5A).

Figure 5

Figure 5

For GO analysis, crotonylated proteins in Q2 were localized on the cell substrate, junction, neuron, and cytoskeleton, while the crotonylated proteins in Q4 were mainly enriched in cellular secretory, cellular membrane, and cellular microstructure in the cellular component (Fig. 5B). Biological process enrichment of crotonylation was performed and crotonylated proteins in Q2 were found to be localized in cellular functions such as migration and development. Crotonylated proteins in Q3 were mainly enriched in cellular metabolic process and cellular response to outside. Moreover, crotonylated proteins in Q4 were mainly enriched in cellular response to the body pathophysiology (Fig. 5C). Crotonylated proteins in molecular function were highly enriched in the cellular binding process and cellular activity of enzymes in Q1, Q2, Q3, and Q4 (Fig. 5D).

KEGG enrichment of the crotonylated proteins revealed some pathways that were related to diseases in Q2 and Q4, such as pancreatic cancer, type II diabetes mellitus, dilated cardiomyopathy, and Salmonella infection, while Q1 crotonylated proteins were enriched in cellular signaling pathway and cellular metabolism (Fig. 5E).

Crotonylated proteins of protein domain were analyzed. These were highly enriched in Q4 with fibrinogen, alpha/beta/gamma chain, C-terminal globular, subdomain 2; zinc finger, LIM-type; fibrinogen, alpha/beta/gamma chain, C-terminal globular domain; fibrinogen, alpha/beta/gamma chain, C-terminal globular, subdomain 1; fibrinogen, alpha/beta/gamma chain, coiled coil domain; and WD40-repeat-containing domain, similar to the functional classification of protein domain (Fig. 5F).

Back to Top | Article Outline

3.6 Protein–protein interaction network of the Kcr proteins

To further identify the cellular processes regulated through crotonylation in MHP, protein−protein interaction network of the Kcr proteins was established (Fig. 6). A total of 888 protein–protein pairs mapped to the protein interaction database, presenting a global view of the diverse cellular functions of crotonylated proteins in MHP. The physiological interactions among these crotonylated protein complexes likely contribute to their cooperation and coordination in MHP.

Figure 6

Figure 6

Back to Top | Article Outline

4 Discussion

To the best of our knowledge, this is the first study to describe Kcr in KF patients on maintenance HD. Several evidences support the role of variations in PTMs in KF. A recent study has associated key events related to KF such as dysregulated expression of genes related to fibrosis, cell cycle, and inflammation with changes in PTMs.[26] Sun et al[32] found that transferrin growth factor (TGF)-β upregulated the expression of H3K4-methyltransferase SET7, which enhanced the expression of TGF-β-induced fibrotic genes. These events may contribute to renal fibrosis changes in KF, implicating the key role of histone methylation in KF. During histone acetylation, TGF-β promotes p300 recruitment with a concomitant increase in H3K9/14Ac and chromatin relaxation at fibrotic gene promoters and Smad2/3 acetylation/activation, leading to an increase in fibrosis in diabetic renal failure.[33] We found an increase in the expression of HDAC-2, -4, and -5 in kidney biopsies from patients with type 1 diabetes. HDAC4-signal transducer and activator of transcription 1 (STAT1) signaling promotes podocyte injury, leading to CKD or KF.[34] Although PTMs have been studied in a variety of model systems, the biochemical function of Kcr in KF patients on maintenance HD is poorly understood.

In this study, we comparatively evaluated the crotonylation proteome of both NC and MHP and identified crotonylated peptides by antibody affinity enrichment, followed by high-resolution mass spectrometry. This approach facilitated the identification of 1109 sites from 347 proteins. Among the quantitative Kcr data set, 155 and 705 sites, corresponding to 93 and 252 proteins were significantly upregulated and downregulated, respectively. According to above data, Kcr levels in histone proteins decreased in MHP compared with NC. Kidney tissue histone crotonylation was shown to increase during acute kidney injury. The elevated histone crotonylation may have a beneficial effect on acute kidney injury, as evident from the increase in the level of peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1α) expression in cultured tubular cells and healthy kidneys.[35] PGC-1α, a major regulator of mitochondrial biogenesis and metabolism, not only suppresses the degree of renal impairment but also helps in the recovery from this disorder.[36] Acute kidney injury is a cause of chronic KF. However, acute kidney injury may be reversed by early treatment to prevent chronic KF. Histone crotonylation may be helpful for the recovery from acute kidney injury, as evident from the increase in the level of histone crotonylation during acute kidney injury. However, it was impossible to restore kidney function in patients on maintenance HD. This observation explains why Kcr levels were lower in histone proteins from MHP as compared with NC.

Functional enrichment of Kcr in GO analysis revealed the crotonylated proteins were associated with diverse biological processes, including cellular structure components, cellular molecular binding, and pathophysiological processes of cellular involvement. Thus, a wide range of interactions involved in these biological processes may be modulated through protein crotonylation. The functional enrichment of Kcr in KEGG analysis revealed only 3 pathways, including complement and coagulation cascades, cardiac muscle contraction, and hematopoietic cell lineage, all of which have important relationship with HD complications. The activation of complement and coagulation cascades system is known to be common during extracorporeal procedures such as HD, wherein the human blood is exposed to biomaterial surfaces.[37] HD is associated with complement and coagulation cascades reactions that cause whole body inflammation and may contribute to accelerated arteriosclerosis, ultimately leading to cardiovascular events.[38] Patients undergoing maintenance HD develop both structural and functional cardiovascular abnormalities.[39] Despite improvement in the dialysis technology, cardiovascular mortality in this population is high. A large proportion of patients on HD dying of cardiac issues fail to show cardiac muscle contraction.[40] Postmortem and biopsy studies have demonstrated that patients with CKD and end-stage renal disease on HD have high levels of interstitial myocardial fibrosis, which was the main reason for the absence of cardiac muscle contraction.[41] Furthermore, HD influences the cardiovascular system and results in hemodynamic disturbances as well as electrolyte shifts, which alter the myocardial electrophysiology. The loss of hematopoietic cell lineage was a commonly diagnosed complication among MHP. Several factors such as the deterioration of kidney function, iron deficiency, poor nutrition, and uremic toxins are known to impair the function of hematopoietic stem cells.[42] Therefore, additional treatment with erythropoietin that drives the bone marrow to produce blood is essential for MHP. At the same time, these patients need HD (for the removal of toxins) and iron supplements as well as highly nutritive substances as the raw material for blood cell.[43] Our study results indicate the possible involvement of Kcr in HD complications. Kcr may reduce HD complications; further research is needed to reveal the role of Kcr in the global proteome of patients on maintenance HD.

Back to Top | Article Outline

Acknowledgments

The authors thank maintenance hemodialysis patients and healthy persons for participating in the basic research of this project.

Back to Top | Article Outline

Author contributions

Conceptualization: Donge Tang.

Data curation: Donge Tang.

Funding acquisition: Dai Yong.

Investigation: Weiguo Sui.

Methodology: Yong Xu.

Project administration: Aoshuang Zou, Hongyan Diao.

Resources: Aoshuang Zou.

Supervision: Dai Yong, Weiguo Sui, Hongyan Diao.

Visualization: Yong Xu.

Writing – original draft: Wenbiao Chen.

Writing – review & editing: Wenbiao Chen.

Back to Top | Article Outline

References

[1]. Kouzarides T. Chromatin modifications and their function. Cell 2007;128:693–705.
[2]. Ruthenburg AJ, Li H, Patel DJ, et al. Multivalent engagement of chromatin modifications by linked binding modules. Nat Rev Mol Cell Biol 2007;8:983–94.
[3]. Bensimon A, Heck AJ, Aebersold R. Mass spectrometry-based proteomics and network biology. Annu Rev Biochem 2012;81:379–405.
[4]. Choudhary C, Kumar C, Gnad F, et al. Lysine acetylation targets protein complexes and co-regulates major cellular functions. Science 2009;325:834–40.
[5]. Kim SC, Sprung R, Chen Y, et al. Substrate and functional diversity of lysine acetylation revealed by a proteomics survey. Mol Cell 2006;23:607–18.
[6]. Zhang Z, Tan M, Xie Z, et al. Identification of lysine succinylation as a new post-translational modification. Nat Chem Biol 2011;7:58–63.
[7]. Peng C, Lu Z, Xie Z, et al. The first identification of lysine malonylation substrates and its regulatory enzyme. Mol Cell Proteomics 2011;10:M111.012658.
[8]. Tan M, Luo H, Lee S, et al. Identification of 67 histone marks and histone lysine crotonylation as a new type of histone modification. Cell 2011;146:1016–28.
[9]. Sabari BR, Tang Z, Huang H, et al. Intracellular crotonyl-CoA stimulates transcription through p300-catalyzed histone crotonylation. Mol Cell 2015;58:203–15.
[10]. Bao X, Wang Y, Li X, et al. Identification of ’erasers’ for lysine crotonylated histone marks using a chemical proteomics approach. Elife 2014;3:e02999.
[11]. Wei W, Liu X, Chen J, et al. Class I histone deacetylases are major histone decrotonylases: evidence for critical and broad function of histone crotonylation in transcription. Cell Res 2017;27:898–915.
[12]. Wu Q, Li W, Wang C, et al. Ultradeep lysine crotonylome reveals the crotonylation enhancement on both histones and nonhistone proteins by SAHA treatment. J Proteome Res 2017;16:3664–71.
[13]. Li Y, Sabari BR, Panchenko T, et al. Molecular coupling of histone crotonylation and active transcription by AF9 YEATS domain. Mol Cell 2016;62:181–93.
[14]. Xiong X, Panchenko T, Yang S, et al. Selective recognition of histone crotonylation by double PHD fingers of MOZ and DPF2. Nat Chem Biol 2016;12:1111–8.
[15]. Liu S, Yu H, Liu Y, et al. Chromodomain protein CDYL acts as a Crotonyl-CoA hydratase to regulate histone crotonylation and spermatogenesis. Mol Cell 2017;67:853.e5–66.e5.
[16]. Wei W, Mao A, Tang B, et al. Large-scale identification of protein crotonylation reveals its role in multiple cellular functions. J Proteome Res 2017;16:1743–52.
[17]. Roth SY, Denu JM, Allis CD. Histone acetyltransferases. Annu Rev Biochem 2001;70:81–120.
[18]. Sauve AA, Wolberger C, Schramm VL, et al. The biochemistry of sirtuins. Annu Rev Biochem 2006;75:435–65.
[19]. Portela A, Esteller M. Epigenetic modifications and human disease. Nat Biotechnol 2010;28:1057–68.
[20]. Gil J, Ramirez-Torres A, Encarnacion-Guevara S. Lysine acetylation and cancer: a proteomics perspective. J Proteomics 2017;150:297–309.
[21]. Kazantsev AG, Thompson LM. Therapeutic application of histone deacetylase inhibitors for central nervous system disorders. Nat Rev Drug Discov 2008;7:854–68.
[22]. Pieterse E, Hofstra J, Berden J, et al. Acetylated histones contribute to the immunostimulatory potential of neutrophil extracellular traps in systemic lupus erythematosus. Clin Exp Immunol 2015;179:68–74.
[23]. Wang Y, Miao X, Liu Y, et al. Dysregulation of histone acetyltransferases and deacetylases in cardiovascular diseases. Oxid Med Cell Longev 2014;2014:641979.
[24]. Zhang C, Zhong JF, Stucky A, et al. Histone acetylation: novel target for the treatment of acute lymphoblastic leukemia. Clin Epigenetics 2015;7:117.
[25]. Li X, Li C, Sun G. Histone acetylation and its modifiers in the pathogenesis of diabetic nephropathy. J Diabetes Res 2016;2016:4065382.
[26]. Reddy MA, Natarajan R. Recent developments in epigenetics of acute and chronic kidney diseases. Kidney Int 2015;88:250–61.
[27]. Levey AS, Atkins R, Coresh J, et al. Chronic kidney disease as a global public health problem: approaches and initiatives - a position statement from Kidney Disease Improving Global Outcomes. Kidney International 2007;72:247–59.
[28]. Terasaki PI, Cecka JM, Gjertson DW, et al. High survival rates of kidney transplants from spousal and living unrelated donors. N Engl J Med 1995;333:333–6.
[29]. Jha V, Garcia-Garcia G, Iseki K, et al. Chronic kidney disease: global dimension and perspectives. Lancet 2013;382:260–72.
[30]. Beckerman P, Ko YA, Susztak K. Epigenetics: a new way to look at kidney diseases. Nephrol Dial Transplant 2014;29:1821–7.
[31]. Himmelfarb J. Hemodialysis complications. Am J Kidney Dis 2005;45:1122–31.
[32]. Sun G, Reddy MA, Yuan H, et al. Epigenetic histone methylation modulates fibrotic gene expression. J Am Soc Nephrol 2010;21:2069–80.
[33]. Yuan H, Reddy MA, Sun G, et al. Involvement of p300/CBP and epigenetic histone acetylation in TGF-beta1-mediated gene transcription in mesangial cells. Am J Physiol Renal Physiol 2013;304:F601–13.
[34]. Wang X, Liu J, Zhen J, et al. Histone deacetylase 4 selectively contributes to podocyte injury in diabetic nephropathy. Kidney Int 2014;86:712–25.
[35]. Ruiz-Andres O, Sanchez-Nino MD, Cannata-Ortiz P, et al. Histone lysine crotonylation during acute kidney injury in mice. Dis Model Mech 2016;9:633–45.
[36]. Tran M, Tam D, Bardia A, et al. PGC-1alpha promotes recovery after acute kidney injury during systemic inflammation in mice. J Clin Invest 2011;121:4003–14.
[37]. Janatova J. Activation and control of complement, inflammation, and infection associated with the use of biomedical polymers. ASAIO J 2000;46:S53–62.
[38]. Nilsson B, Ekdahl KN, Mollnes TE, et al. The role of complement in biomaterial-induced inflammation. Mol Immunol 2007;44:82–94.
[39]. Chirakarnjanakorn S, Navaneethan SD, Francis GS, et al. Cardiovascular impact in patients undergoing maintenance hemodialysis: clinical management considerations. Int J Cardiol 2017;232:12–23.
[40]. de Jager DJ, Grootendorst DC, Jager KJ, et al. Cardiovascular and noncardiovascular mortality among patients starting dialysis. JAMA 2009;302:1782–9.
[41]. Aoki J, Ikari Y, Nakajima H, et al. Clinical and pathologic characteristics of dilated cardiomyopathy in hemodialysis patients. Kidney Int 2005;67:333–40.
[42]. Mikhail A, Brown C, Williams JA, et al. Renal association clinical practice guideline on Anaemia of Chronic Kidney Disease. BMC Nephrol 2017;18:345.
[43]. Bonomini M, Del Vecchio L, Sirolli V, et al. New treatment approaches for the anemia of CKD. Am J Kidney Dis 2016;67:133–42.
Keywords:

lysine crotonylation; maintenance hemodialysis; mass spectrometry

Supplemental Digital Content

Back to Top | Article Outline
Copyright © 2018 The Authors. Published by Wolters Kluwer Health, Inc. All rights reserved.