Impact Factor 4.076

The 3rd most cited journal in Microbiology

This article is part of the Research Topic

Foodborne Pathogens: Hygiene and Safety

Original Research ARTICLE

Front. Microbiol., 11 January 2018 | https://doi.org/10.3389/fmicb.2017.02674

Comprehensive Proteomic Analysis of Lysine Acetylation in the Foodborne Pathogen Trichinella spiralis

Yong Yang1,2†, Mingwei Tong3†, Xue Bai1†, Xiaolei Liu1, Xuepeng Cai4,5, Xuenong Luo5, Peihao Zhang2, Wei Cai2, Isabelle Vallée6, Yonghua Zhou7* and Mingyuan Liu1*
  • 1Key Laboratory of Zoonosis Research, Ministry of Education, Institute of Zoonosis/College of Veterinary Medicine, Jilin University, Changchun, China
  • 2Wu Xi Medical School, Jiangnan University, Affiliated Hospital of Jiangnan University, Wuxi, China
  • 3State Key Laboratory for Molecular Biology of Special Economic Animals, Institute of Special Economic Animal and Plant Sciences, Chinese Academy of Agricultural Sciences, Changchun, China
  • 4China Institute of Veterinary Drug Control, Beijing, China
  • 5State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou, China
  • 6JRU BIPAR, ANSES, École Nationale Vétérinaire d’Alfort, INRA, Université Paris-Est, Animal Health Laboratory, Maisons-Alfort, France
  • 7Jiangsu Institute of Parasitic Disease, Wuxi, China

Lysine acetylation is a dynamic and highly conserved post-translational modification that plays a critical role in regulating diverse cellular processes. Trichinella spiralis is a foodborne parasite with a considerable socio-economic impact. However, to date, little is known regarding the role of lysine acetylation in this parasitic nematode. In this study, we utilized a proteomic approach involving anti-acetyl lysine-based enrichment and highly sensitive mass spectrometry to identify the global acetylated proteome and investigate lysine acetylation in T. spiralis. In total, 3872 lysine modification sites were identified in 1592 proteins that are involved in a wide variety of biological processes. Consistent with the results of previous studies, a large number of the acetylated proteins appear to be involved in metabolic and biosynthetic processes. Interestingly, according to the functional enrichment analysis, 29 acetylated proteins were associated with phagocytosis, suggesting an important role of lysine acetylation in this process. Among the identified proteins, 15 putative acetylation motifs were detected. The presence of serine downstream of the lysine acetylation site was commonly observed in the regions surrounding the sites. Moreover, protein interaction network analysis revealed that various interactions are regulated by protein acetylation. These data represent the first report of the acetylome of T. spiralis and provide an important resource for further explorations of the role of lysine acetylation in this foodborne pathogen.

Introduction

Trichinellosis is a common food-borne parasitic zoonosis worldwide. Infection occurs by consuming raw or inadequately cooked meat containing muscle larvae of the Trichinella parasite. This parasite can infect a wide variety of hosts, including humans. Trichinella infection has been reported in 66 countries around the world. Trichinella spiralis is the major causative agent of trichinellosis. Outbreaks of trichinellosis in humans have been regularly reported worldwide, and serological analysis has indicated that more than 11 million people may be infected with Trichinella (Murrell and Pozio, 2011). Trichinellosis is regarded as an emerging or re-emerging infectious disease, particularly in developing countries, where cooking habits, poor sanitary conditions, and a lack of veterinary inspection facilitate infection (Gottstein et al., 2009). This zoonosis is both a public health hazard and an economic issue in terms of the porcine breeding industry and food safety (Bai et al., 2017). Therefore, in order to control and eradicate this disease, it is necessary to develop novel diagnostic and therapeutic methods for trichinellosis, which will require a comprehensive understanding of the parasite’s biology.

In recent years, the genome, transcriptome, epigenome, and proteome of T. spiralis have been reported and released in public databases, providing a greater understanding of the molecular basis of T. spiralis biology, host–parasite interactions, and pathogenesis (Mitreva et al., 2011; Gao et al., 2012; Liu et al., 2012, 2016). However, to the best of our knowledge, no post-translational modifications (PTMs), which occur during or after protein biosynthesis, have been described in T. spiralis, limiting further functional studies. PTMs regulate diverse protein properties, including folding, stability, localization, binding affinities, and activity, and they play important roles in a variety of cellular pathways and processes. Because of their substantial effects on protein functions and cellular processes, studies of protein PTMs have become more important in the post-genomic era (Walsh et al., 2005; Witze et al., 2007). To date, more than 461 different PTMs, such as phosphorylation, methylation, ubiquitination, malonylation, succinylation, glycosylation, and acetylation, have been identified to play a role in the functional regulation of many prokaryotic and eukaryotic proteins. These PTMs form a complex regulatory system for the dynamic control of cellular processes under various conditions (Yang and Seto, 2008).

Among these PTMs, lysine acetylation (Kac) is a highly dynamic and regulated modification that is defined as the reversible transfer of an acetyl group from the acetyl donor acetyl-CoA to the 𝜀-amino group of a lysine residue of a protein molecule. The Kac process is regulated by lysine acetyltransferases and deacetylases (KDACs), resulting in a regulatory range that is extensive and comparable with those of other major PTMs (Choudhary et al., 2009; Drazic et al., 2016). KDACs have been identified as drug target to develop novel anti-parasite treatments, and some selective KDAC ligands showed the potential for the treatment of parasitic diseases (Wang et al., 2015).

Recent studies have indicated that non-enzymatic chemical acetylation occurs broadly in eukaryotic cells, suggesting that this PTM may have various roles in protein function and cellular processes (Wagner and Hirschey, 2014; Olia et al., 2015). Kac was first identified in eukaryotic histones, where it is involved in regulating chromatin structure and transcription. In addition to its role in histones, Kac were found in a number of non-nuclear proteins which are located in almost every compartment of the cell (Friedmann and Marmorstein, 2013; Drazic et al., 2016). Protein acetylation outside the nucleus participates in diverse cellular physiological processes, including cell-cycle regulation, cytoskeletal dynamics, apoptosis, autophagy, and metabolism (Kim et al., 2006; Close et al., 2010; Wang et al., 2010; Zhao et al., 2010; Sun et al., 2015). In addition, acetylation has also been implicated to influence the persistence, virulence and antibiotic resistance of pathogenic bacteria (Xie et al., 2015).

Recent advances in antibody-based affinity purification of acetylated peptides and highly sensitive mass spectrometry (MS) have made significant contributions to the study of Kac at the whole-proteome level. Global Kac has been characterized in a number of eukaryotic and prokaryotic organisms (Jeffers and Sullivan, 2012; Zhang et al., 2013; Pan et al., 2014; Nie et al., 2015; Xie et al., 2015). These proteome-wide analyses of Kac have mapped large datasets of lysine-acetylated proteins and have confirmed that Kac is involved in various cellular events, especially in modulating metabolic processes. Based on these studies, it is proposed that the functional roles of acetylation in these processes are evolutionarily conserved in both eukaryotes and prokaryotes, and acetylation may thus rival phosphorylation.

Although Kac has been characterized in many organisms, the study of the lysine acetylome in helminths has been relatively limited, with only one nematode species (Schistosoma japonicum) being previously characterized (Hong et al., 2016). Given the general use of acetylation as a regulatory mechanism in different organisms, it can be speculated that T. spiralis also employs acetylation as a modification to modulate various aspects of its biology. Lending evidence to this hypothesis, T. spiralis appears to contain a great number of acetyltransferases and deacetylases in its genome (Mitreva et al., 2011).

To bridge this knowledge gap, in this study, proteomics was used in combination with immunoprecipitation in the first large-scale analysis of lysine-acetylated proteins in T. spiralis, identifying 3872 Kac sites in 1592 proteins. These lysine-acetylated proteins are involved in a variety of biological functions and are localized to multiple cellular compartments. Bioinformatics analysis showed that acetylation in T. spiralis modulates a wide range of cell processes. This study provides a global picture of the T. spiralis acetylome, which not only greatly expands our knowledge of helminth protein acetylation but also provides important information for the design of new effective drugs or vaccines for controlling trichinellosis.

Materials and Methods

Parasite Collection and Protein Extraction

Animal experiments were performed according to the guidelines of the National Institute of Health (publication No. 85–23, revised 1996) and Animal Experimentation Guidelines of Jilin University and were approved by the Animal Experimental Ethics Committee of Jilin University, China. T. spiralis (ISS534) parasites were maintained by serial passage in female BALB/c mice. T. spiralis muscle larvae (ML) were recovered from infected BALB/c mice at 35 days post-infection (dpi) using a standard HCl-pepsin digestion method. ML were washed manually in phosphate-buffered saline (PBS) at 37°C to remove any residual host proteins and were then collected for further lysine acetylome analysis.

Total proteins were extracted as previously described (Zhou et al., 2016) with minor modifications. In brief, parasites were lysed with lysis buffer containing 8 M urea, 10 mM dithiothreitol (DTT) and 0.1% protease inhibitor cocktail, followed by sonication. Cellular debris was then removed by centrifugation at 20,000 × g and 4°C for 10 min. Proteins contained in the supernatant were precipitated with cold 15% trichloroacetic acid (TCA) for 2 h at -20°C. After centrifugation, the pellet was washed three times with cold acetone, followed by dissolution in buffer containing 8 M urea and 100 mM NH4CO3 (pH 8.0). Protein concentration was measured using a 2-D Quant kit (GE Healthcare) according to the manufacturer’s instructions.

Trypsin Digestion and High-Performance Liquid Chromatography (HPLC) Fractionation

Proteins were reduced with 10 mM DTT, alkylated with 20 mM iodoacetamide (IAA), and digested with trypsin (Promega) overnight at a 1:50 trypsin:protein mass ratio. To ensure complete digestion, additional trypsin was added at a 1:100 trypsin:protein mass ratio for 4 h for a second digestion. Tryptic peptides were separated into six fractions using high-pH reverse-phase HPLC, as previously described (Zhou et al., 2016). Separated peptides were then dried completely in a SpeedVac (Thermo Scientific) and stored at -80°C for further enrichment of acetylated peptides.

Affinity Enrichment and Liquid Chromatography Electrospray Ionization Tandem Mass Spectrometry (LC-ESI–MS/MS) Analysis

Fractionated peptides were dissolved in NETN buffer containing 100 mM NaCl, 1 mM EDTA, 50 mM Tris–HCl, and 0.5% NP-40 (pH 8.0) and incubated with anti-acetyl lysine antibody beads (PTM Biolabs) at 4°C overnight with gentle shaking. After washing with NETN buffer and deionized H2O, the bound peptides were eluted with 0.1% trifluoroacetic acid, dried in a SpeedVac, and cleaned using C18 ZipTips (Millipore) according to the manufacturer’s instructions.

The obtained peptides were then analyzed by LC-ESI–MS/MS on a Q ExactiveTM Plus hybrid quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific) coupled online to the EASY-nLC 1000 ultra-performance liquid chromatography (UPLC) system. The resulting MS/MS data were processed using MaxQuant with an integrated Andromeda search engine (v.1.5). Tandem mass spectra were queried against the UniProt T. spiralis database using the following parameters (Zhou et al., 2016): trypsin/P was specified as the cleavage enzyme and allowed up to four missed cleavages, five modifications per peptide, and five charges. Mass error was set to 10 ppm for precursor ions and 0.02 Da for fragment ions. Carbamidomethylation on Cys was specified as a fixed modification, and oxidation on Met, acetylation on Lys, and acetylation on the protein N-terminus were specified as variable modifications. False discovery rate (FDR) thresholds for the protein, peptide, and modification sites were set at 1%. The minimum peptide length was set at 7. The site localization probability was set as > 0.75.

Bioinformatics Analysis

Gene Ontology (GO) annotation of the proteome was derived from the UniProt-GOA database1 according to Zhou et al. (2016). All identified acetylated proteins were classified using three categories of GO annotations: biological process, cellular component, and molecular function. The subcellular localization predication program, Wolfpsort (Horton et al., 2007), was used to analyze cellular localization data, and the InterPro database2 and InterProScan were used to compile functional descriptions of protein domains. Protein pathways were annotated using the Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2004). Amino acid sequence motifs in identified protein were analyzed using Soft motif-x, and a position-specific heat map was generated using WebLogo 3.43. Secondary structures of all acetylated proteins were predicted using NetSurfP software. Fisher’s exact test was used to analyze GO and KEGG pathway and domain enrichment. Terms with adjusted p-values < 0.05 were considered significantly enriched. The Search Tool for Retrieval of Interacting Genes/Proteins (STRING) database4 was used for analysis of the protein–protein interaction (PPI) network. The interaction network form STRING was visualized using Cytoscape (version 3.0) software (Bauer-Mehren, 2013).

Results and Discussion

Identification and Analysis of Lysine-Acetylated Proteins in T. spiralis

Sequencing of the T. spiralis genome revealed the existence of a number of acetyltransferase and deacetylase orthologs, suggesting that Kac may play a vital role in T. spiralis development and metabolism (Mitreva et al., 2011). To comprehensively characterize protein acetylation in T. spiralis, analysis of Kac was carried out at the proteome level using tryptic digestion, anti-acetyl lysine antibody enrichment, and high-resolution LC-MS/MS. Generally, the distribution of mass errors was near zero, and most were ≤ 10 ppm (Figure 1A). Moreover, all lysine-acetylated peptides were between 7 and 37 amino acids in length, which is consistent with tryptic peptides (Figure 1B). These results indicated that the sample preparation method was adequate and that the modified peptide data obtained from MS was highly accurate.

FIGURE 1
www.frontiersin.org

FIGURE 1. Proteome-wide identification of lysine acetylation sites in T. spiralis. (A) Mass error distribution of all identified peptides. (B) Distribution of acetylated peptides based on length. (C) Distribution of acetylated proteins based on number of acetylation sites.

In total, 3872 Kac sites distributed across 1592 acetylated proteins were identified (Supplementary Table S1). This exceeds the level of Kac that has so far been reported in S. japonicum (Hong et al., 2016), indicating the essential role of Kac in T. spiralis. The acetylated proteins of T. spiralis contained various numbers of acetylation sites, ranging from 1 to 47 in the putative myosin head protein (E5SJK9) (Figure 1C), with an average frequency of 2.4 sites per protein. Among all acetylated proteins, approximately half (795/1592) contained only one acetylated site, and the percentages of proteins with two, three, four, or > 4 modification sites were 20.4, 11.4, 6.5, and 11.7%, respectively. Acetylation of non-histone proteins at multiple sites have shown to modulate every step of cellular processes from signaling to transcription to protein degradation through various molecular mechanisms (Huang and Berger, 2008; Spange et al., 2009). The role of multiple acetylation sites in regulating protein function remains to be elucidated.

Analysis of Kac Motifs

Previous studies in both eukaryotic and prokaryotic cells have shown preferences for amino acid residues at specific positions surrounding acetylated lysines (Nie et al., 2015; Zhou et al., 2016). In our study, a total of 15 conserved amino acid sequence motifs from -10 to +10 surrounding the acetylated lysine were defined from among 3440 peptides (accounting for 88.8% of the total identified peptides). These motifs included KacS, KacR, KacH, KacN, E∗∗KacK, KacK, KacT, KKac, KacR, TKacV, KacD, Kac∗∗R, KacE, KacV, and KacD, an asterisk indicates a random amino acid residue. These motifs exhibited different abundances (Figure 2A and Supplementary Table S2, with the KacS, KacR, KacH, and KacN motifs being the most frequent and accounting for 18.7, 16.1, 11.4, and 11.3% of all identified peptides, respectively (Figure 2B). In addition, most of the conserved residues were located at the ± 1 or ± 2 positions of the Kac sites. The motif analysis results suggest that two particular types of residues are found near Kac sites. Enrichment of residues with major hydrophilic side chain groups, such as serine (S), arginine (R), and aspartic acid (D), was observed at the +1 position, and enrichment of threonine (T) with hydrophobic side chain groups was observed at the -1 position. These results suggest that amino acid residues with hydrophilic and hydrophobic side chains may be more easily attacked by Kac.

FIGURE 2
www.frontiersin.org

FIGURE 2. Bioinformatics study of lysine-acetylated sites. (A) Probable sequence motifs of T. spiralis acetylation sites consisting of 20 residues surrounding the targeted lysine residue, identified using Motif-X. (B) Number of identified peptides containing acetylated lysines and each motif. (C) Heat map showing the relative frequencies of amino acids in specific positions, including enrichment (red) or depletion (green) of amino acids flanking the acetylated lysine in T. spiralis proteins. (D) Probabilities of lysine acetylation in different protein secondary structures (alpha-helix, beta-strand, and coil). (E) Predicted surface accessibility of acetylation sites. All lysine sites are shown in green, and acetylated lysine sites are in red.

Studies of acetylation motifs have shown that aromatic amino acids frequently surround modification sites (Hong et al., 2016; Wang et al., 2017). Unexpectedly, we noticed only one aromatic amino acid (histidine, H) surrounding acetylated lysines, suggesting that the lysine acetyltransferases of T. spiralis might not prefer polypeptides with aromatic amino acid residues as substrates. Furthermore, some acetylated lysine motifs found in T. spiralis have also been identified in other organisms, such as KacH in humans (Shao et al., 2012) and Mycobacterium tuberculosis (Xie et al., 2015); Kac∗∗R in humans (Shao et al., 2012); KacE in rats (Lundby et al., 2012); and KacR, KacT, and KacV in Toxoplasma gondii (Jeffers and Sullivan, 2012), confirming that Kac is a highly conserved modification among various species. Moreover, LceLogo heat maps of the amino acid compositions surrounding the acetylation sites showed that the frequencies of serine (S), threonine (T), and histidine (H) in positions -2 to + 2 were the highest, while the frequency of tryptophan (W) was the lowest (Figure 2C) in regions surrounding these sites. Based on these findings, we infer that proteins with S, T, and H but not W around lysine residues are the favored targets of lysine acetyltransferases in T. spiralis. A structural analysis of the acetylated proteins was carried out using NetSurfP software in order to understand the relationship between acetylation and the presence of specific protein secondary structures in T. spiralis. As shown in Figure 2D, acetylation sites (56.8%) were more common in regions with ordered secondary structures.

Among the acetylated sites, 43.8% of the acetylated sites were located in alpha-helices, and 13.2% were located in beta-strands. The remaining 43.2% of the acetylated sites were distributed in unstructured protein regions. These results suggest that acetylases may prefer ordered structural regions to disordered regions in T. spiralis proteins, which differs from findings of the acetyl proteomes of humans and S. japonicum (Hong et al., 2016). The acetylation sites in these species have been found to be located primarily in disordered protein regions, with less acetylation in ordered protein structures, suggesting that secondary structure preferences at lysine-acetylated sites may vary among species. However, a similar distribution pattern of secondary structures was observed for non-acetylated lysines in the proteome of T. spiralis, indicating that there is no tendency toward acetylation in T. spiralis. Finally, the surface accessibility of the acetylated lysine sites was also studied. The results showed that the exposure of acetylated residues on the protein surface is very similar to that of non-acetylated lysine residues (Figure 2E). Therefore, Kac may not be affected by the surface properties of proteins in T. spiralis.

Evolutionary Conservation of Kac in T. spiralis

Increasing evidence has indicated that Kac maybe evolutionarily conserved among various species (Zhou et al., 2016). To further improve our understanding of the evolution of acetylation in T. spiralis, we compared the acetylated proteins in T. spiralis with sacetylated proteins that have been reported in other organisms based on the quantity of orthologs. The published acetylated proteins of Escherichia coli (Zhang et al., 2013), budding yeast (Saccharomyces cerevisiae) (Henriksen et al., 2012), Candida albicans (Zhou et al., 2016), T. gondii (Jeffers and Sullivan, 2012), S. japonicum (Hong et al., 2016), mouse (liver tissue), and human [cervical cancer (HeLa) cells] (Weinert et al., 2013) were selected. These acetylome data are comparable, as similar methods were used for the systematic identification of Kac in each of these species. Among the 1592 acetylated proteins in T. spiralis, 408 (35.5%) proteins had orthologs in the S. japonicum acetylome, 479 (35.2%) had orthologs in the mouse acetylome, 471 (36.3%) had orthologs in the human acetylome, 159 (31.6%) had orthologs in the S. cerevisiae acetylome, 190 (39.8%) had orthologs in the C. albicans acetylome, 87 (30%) had orthologs in the T. gondii acetylome, and 43 (12.3%) had orthologs in the E. coli acetylome (Figure 3 and Supplementary Table S3). The lowest ortholog rate was between T. spiralis and E. coli, and this could reflect differences in cellular structures between prokaryotes and eukaryotes. The highest overlap was between T. spiralis and C. albicans, which can be explained by the genetic relatedness between them and is likely due to a similar subcellular distribution of acetylated proteins. These results indicate that acetylation occurs frequently in both eukaryotes and prokaryotes with certain species.

FIGURE 3
www.frontiersin.org

FIGURE 3. Analysis comparing the T. spiralis acetyl proteome with those of E. coli, budding yeast (S. cerevisiae), C. albicans, T. gondii, S. japonicum, mouse, and human. Venn diagram describing the relationships among lysine-acetylated proteins from these species.

Functional Annotation and Cellular Localization of Acetylated Proteins in T. spiralis

To functionally characterize the identified lysine-acetylated proteins, we first investigated the GO functional classifications of all acetylated proteins (Figure 4 and Supplementary Table S4). As shown in Figure 4A, acetylated proteins were distributed among a variety of cellular compartments, primarily in cells (37%), organelles (24%), macromolecular complexes (22%), and membranes (13%). According to the molecular function classification, the largest group of acetylated proteins was composed of those involved in the binding of various targets, accounting for 40% of all identified acetylated proteins (Figure 4B). This result suggests that acetylation may play essential roles in regulating PPIs in T. spiralis. The second largest group in the molecular function category included proteins related to catalytic activity (37%). Enzymes play important roles in regulating biochemical reactions in organisms because of their catalytic activity. Studies have indicated that reversible Kac can regulate metabolic enzymes by increasing or decreasing enzyme quantity, catalytic activity, or substrate accessibility (Xiong and Guan, 2012). Our results are consistent with those of previous studies, in which a considerable number of lysine-acetylated proteins with catalytic activity were enzymes involved in metabolic processes (Jeffers and Sullivan, 2012; Hong et al., 2016).

FIGURE 4
www.frontiersin.org

FIGURE 4. Gene Ontology functional classification and subcellular location information of identified acetylated proteins. (A) Acetylated proteins classified according to cellular components. (B) Acetylated proteins classified according to biological process. (C) Acetylated proteins classified according to molecular function. (D) Subcellular localizations of identified acetylated proteins.

In the biological process classification of the GO annotation of acetylated proteins, two main classes of acetylated proteins were involved in cellular and metabolic processes, accounting for 29 and 28% of identified proteins, respectively (Figure 4C). Previous studies have demonstrated that reversible Kac is an emerging regulatory mechanism for proteins involved in metabolism and cellular processes in bacteria (Pan et al., 2014), parasites (Jeffers and Sullivan, 2012), animals, and humans (Weinert et al., 2013). Here, our results indicate that Kac may also play a critical role in T. spiralis metabolism and cellular regulation.

The localization of proteins to various subcellular compartments is important for allowing proteins to perform their necessary functions. To obtain a profile of acetylated proteins throughout the cell, we analyzed the subcellular localization of the 1592 acetylated proteins in T. spiralis (Figure 4D). Analysis of the acetylated proteins showed that the majority of identified Kac proteins were localized to the cytoplasm (38%), followed by the nucleus (17%), and mitochondria (12%). This result is consistent with those of the biological process analysis, which showed that most acetylated proteins were involved in protein cellular metabolism and protein synthesis. Importantly, we also found that a number of acetylated proteins were cytoplasmic membrane (12%) or extracellular (12%) proteins, which could function in regulating cuticular layer-related proteins, impacting the interactions between T. spiralis and host cells. In addition, 1% of the acetylated proteins in this study were located in the peroxisome, where fatty acids are activated and converted into acyl-CoA before being transformed into sucrose.

Functional Enrichment Analysis

In order to reveal the preferred targets of Kac, we performed a GO enrichment analysis of the acetylation data. The GO enrichment analysis of the molecular function category further demonstrated that acetylated proteins were significantly enriched in those related to binding activity and the structural constituents of ribosomes (Figure 5A). GO enrichment analysis based on the biological process category showed that acetylated proteins were enriched in those involved in biosynthetic and metabolic processes (Figure 5A), suggesting that these processes may be precisely regulated by protein acetylation. As acetylated proteins may be involved in both of these fundamental processes, Kac appears to be a pivotal regulatory mechanism in T. spiralis.

FIGURE 5
www.frontiersin.org

FIGURE 5. Enrichment analysis of acetylated proteins in T. spiralis. (A) GO-based enrichment analysis of acetylated proteins in terms of molecular function, cellular component, and biological process. (B) KEGG pathway enrichment analysis. (C) Protein domain enrichment analysis.

To gain in-depth insights into the functions of Kac proteins in T. spiralis, KEGG enrichment analysis was also carried out (Figure 5B). Our data showed that the acetyl proteome of T. spiralis was highly enriched in proteins involved in protein biosynthesis and the metabolism of three types of biological macromolecules (Figure 5B). The significant enrichment of proteins involved in various pathways, including translation, transcription, and metabolism, among acetylated proteins has also been reported in other eukaryotes and prokaryotes (Nie et al., 2015; Hong et al., 2016; Zhou et al., 2016), confirming the essential role of Kac in various organisms.

Small protein motifs involved in protein interactions are also regulated by PTMs. Thus, we performed a domain enrichment analysis. In total, acetylated proteins were significantly enriched for 25 protein domains, including the NAD(P)-binding domain, AAA+ ATPase domain, and nucleic acid-binding (OB-fold) domain (Figure 5C). Proteins with these functional domains are involved in a variety of metabolic pathways. A previous study indicated that enzymes with an NAD(P)-binding domain catalyze redox reactions, suggesting that Kac plays an important role in the stress response (Hua et al., 2014). The OB-fold is a compact structural motif that recognizes single-stranded and unusually structured nucleic acids, and it plays a critical role in DNA replication, DNA repair, transcription, translation, and telomere maintenance (Theobald et al., 2003), making these proteins attractive anti-parasitic drug targets. In our study, a large number of proteins containing nucleic acid-binding domains, such as helicase, tyrosine–tRNA ligase, 40S ribosomal protein S11, GTP-binding protein, asparaginyl–tRNA synthetase, and nuclease domain-containing protein 1, were found to have different numbers of acetylation sites. Prevalent acetylation of proteins containing nucleic acid-binding domains may play an important role in directing the recognition of single-stranded and unusually structured nucleic acids. Overall, our results indicate that acetylated proteins are widely distributed throughout the cell, with important effects on a variety of processes in T. spiralis.

PPI Analysis

Protein–protein interactions occur widely in all cell types and are important in regulating cellular processes. To better understand how interactions among acetylated proteins may affect the regulation of T. spiralis physiology and development, we assembled a PPI network containing all of the identified acetylated proteins using Cytoscape. The results showed that 666 acetylated proteins represented network nodes that were connected by 9446 direct physical interactions, with a combined score > 0.90 obtained from the STRING database. The majority of proteins in this interaction network contain multiple Kac sites, and the network suggests that the physiological interactions among these acetylated proteins are involved in multiple protein interaction networks and control various signaling pathways in T. spiralis. A complete T. spiralis network of acetylated proteins is shown in Supplementary Table S5 and represents the first high-quantity interaction network of acetylated proteins in a foodborne parasite.

On the basis of the interaction network, we retrieved four highly interconnected interaction clusters of acetylated proteins using the MCODE algorithm in Cytoscape.

The top cluster, with 73 acetylated proteins (cluster I), consisted of ribosome-associated proteins that were connected in a PPI network with a relatively high density (Figure 6A). Clusters II and III consisted of proteins associated with the proteasome and RNA transport, whereas the fourth highest-ranking complex consisted of proteins involved in the tricarboxylic acid (TCA) cycle (Figures 6B–D). These subnetworks also have relatively high densities, and many of the acetylated proteins in the network contain multiple Kac sites. These findings suggest that the protein metabolic metabolite machinery is probably regulated by Kac in T. spiralis; the identified proteins were also involved in central metabolic pathways.

FIGURE 6
www.frontiersin.org

FIGURE 6. Interaction networks of acetylated proteins in T. spiralis. (A) Interaction network of acetylated proteins associated with the ribosome. (B) Interaction network of acetylated proteins associated with the proteasome. (C) Interaction network of acetylated proteins associated with RNA transport. (D) Interaction network of acetylated proteins associated with the TCA cycle.

Analysis of Histone Kac

Histone acetylation occurs in diverse organisms and is now known to play a critical role in chromatin remodeling, DNA repair, and the epigenetic regulation of gene expression (Campos and Reinberg, 2009). The acetylation of histones 3 and 4 is commonly associated with active transcription, and histone acetylation at H3K27 play critical roles in regulating enhancer activities in cooperation with histone methylation (Hyun et al., 2017; Shen et al., 2017; Youn, 2017) Histone acetylation is also essential for parasite survival. In S. mansoni, histone deacetylase inhibitors induce mortality and apoptosis, which are associated with the increased expression of the caspase 3 and 7 transcripts due to histone H4 hyperacetylation. Thus, the histone deacetylase of S. mansoni may be an important target for schistosomiasis chemotherapy (Dubois et al., 2009). Many specific histone acetylation sites have been characterized in various parasites. The acetylation of H4K8 and H4K16 in Giardia lamblia plays a crucial role in regulating encystation, which protects the parasite from harsh environmental conditions (Carranza et al., 2016). Similarly, variations in the levels of H4K8 acetylation are important in regulating stage-specific gene expression throughout the Plasmodium falciparum life cycle (Carranza et al., 2016). In the present study, we identified six histone proteins with multiple Kac sites (Supplementary Table S6), with H2B (E5RY73) being heavily acetylated at eight lysine residues. Some acetylated lysine sites identified in our study have also been found in other parasites, such as T. gondii (Jeffers and Sullivan, 2012) and P. falciparum (Miao et al., 2013). However, further research is required to determine the significance of these modifications in T. spiralis and their roles in development and infection.

Protein Acetylation Regulates Diverse Metabolic Pathways in T. spiralis

Recently, proteomic studies into acetylation have indicated that abundant proteins are acetylated in various organisms, suggesting that acetylation is a common mechanism of metabolic regulation (Hong et al., 2016; Zhou et al., 2016). Consistent with previous results, our study showed that 138 acetylated proteins were involved in multiple metabolic pathways. Here, we focused on two of these pathways: glycometabolism and protein metabolism.

Glycometabolism

The acetylation of proteins involved in glycometabolism (mainly enzymes) has been showed to play a crucial regulatory role in processes such as parasite development and growth (Miao et al., 2013; Hong et al., 2016). In our study, a large number of enzymes in these pathways were observed to undergo acetylation in T. spiralis (Figure 7), indicating that acetylation also participates in regulating glycometabolism in T. spiralis.

FIGURE 7
www.frontiersin.org

FIGURE 7. Lysine-acetylated proteins involved in the central carbon metabolic network in T. spiralis. Identified acetylated enzymes are shown in red.

Trichinella spiralis larvae ferment glycogen under aerobic or anaerobic conditions via phosphorylative glycolysis, converting it to a mixture of volatile fatty acids (Agosin and Aravena, 1959). We found that almost all glycolytic enzymes involved in the conversion of glucose to pyruvate were substrates of Kac (15 proteins). Previous studies have found that seven out of ten glycolytic enzymes are acetylated in the E. coli acetylome (Zhang et al., 2013), while only glucose-6-phosphate isomerase is not acetylated in P. falciparum (Miao et al., 2013). Similarly, all enzymes involved in glycolysis are acetylated in the S. japonicum acetylome (Hong et al., 2016). In our study, only glyceraldehyde-phosphate dehydrogenase (GAPDH) was not found to be acetylated. In glycolysis, GAPDH acetylation increases glycolytic activity. Thus, the acetylation of GAPDH is important for controlling carbon flux, which is beneficial for glycolysis in cells supplied with glucose (Li et al., 2014). As the GAPDH protein of T. spiralis was not acetylated in our study, other enzymes, such as fructose-bisphosphate aldolase or phosphoglycerate kinase (PGK), with six and eight Kac sites, respectively, may be the main targets for the regulation of glycolysis in T. spiralis. Previous study have indicated that acetylation of fructose 1,6-bisphosphate aldolase at K14 increased enzymatic activity (Kim et al., 2006). PGK is another central enzyme in the glycolysis pathway, catalyzing the conversion of 1,3-diphosphoglycerate to 3-phosphoglycerate and generating a molecule of ATP. The proliferation of cancer cells demands a lot of energy and is mainly dependent on anaerobic glycolysis for the conversion of glucose, a phenomenon known as the “Warburg effect.” Glycolysis is much less potent than mitochondrial oxidative phosphorylation in terms of ATP production, as only two enzymes, PGK and pyruvate kinase, control the production of ATP during anaerobic glycolysis in cancer cells. A recent study found that the acetylation of PGK at K323 enhances PGK activity and promotes cancer cell proliferation and tumorigenesis (Hu et al., 2017). In our study, the lysine sites in PGK were heavily acetylated, which may facilitate nutrient absorption and T. spiralis growth.

In addition, in our study, enolase exhibited the greatest number of acetylated sites (12 sites). In addition to their glycolytic function, enolases were recently shown to interact with plasminogen to induce fibrinolytic activity and to promote the migration of larval stages through tissues by plasmin-mediated proteolysis, such as the degradation of the host’s extracellular matrix. Enolase can induce some protection against schistosomal infection and also inhibits egg production in schistosomes (Yang et al., 2010). Therefore, enolase is considered a potential target for the development of new anti-parasitic vaccines and drugs. However, more work is required to fully elucidate the role of acetylation in controlling the function of enolase.

Another interesting pathway, the pentose phosphate pathway contained 10 acetylated proteins according to our results. This pathway provides reductive potential by producing NADPH, a cellular reductant that protects against antioxidant enzyme activity. It has also been shown to play an essential protective role against metalloid-mediated oxidative stress by stringently regulating the accumulation of reactive oxygen species (ROS) (Ghosh et al., 2017). Thus, the pentose phosphate pathway is considered to play an essential role in protecting cells against oxidative stress (Ghosh et al., 2015). Based on our results, we speculate that Kac plays a vital role in protecting the biomolecules of T. spiralis from oxidative damage.

Carbon is metabolized continuously in parasites and is essential for energy circulation and parasite survival. The TCA cycle is essential for supplying the energy needed for numerous cellular functions. Unlike Schistosoma, which completely relies upon glucose fermentation and the formation of lactate to meet its energy requirements, T. spiralis utilizes the oxidation of carbohydrates to CO2 and water via the TCA cycle as its major energy source (Goldberg, 1957).

A key entry point for carbon into the TCA cycle for energy production is the pyruvate dehydrogenase complex (PDC), which catalyzes the overall conversion of the glycolytic degradation products of carbohydrate metabolism and pyruvate to acetyl-CoA and CO2. The PDC is composed of three different subunits: pyruvate decarboxylase (E1, PDH), dihydrolipoamide transacetylase (E2, DLAT), and dihydrolipoamide dehydrogenase (E3, DLD). All of these subunits are required to form NADH from the oxidation of pyruvate. In our study, we identified five acetylated proteins in this complex. Of note, the PDC in glycolysis was acetylated at 22 sites among its three subunits, with six and nine lysine sites on DLAT and DLD, respectively, modified by Kac. In human cancer cells, K321 acetylation of PDHA1 and K202 acetylation of PDH phosphatase (PDP1) inhibit PDC activity (Fan et al., 2014). Thus, the intensive acetylation of acetyl-CoA metabolism-related enzymes in T. spiralis may regulate the catalytic activity of the complex, especially PDH, and control its functional balance in T. spiralis.

Many mitochondrial proteins that play an essential role in energy metabolism exhibit a high level of acetylation in many organisms (Hong et al., 2016; Zhou et al., 2016). Previous studies have shown that the acetylation of malate dehydrogenase 1 (MDH1), a key enzyme in the TCA cycle, significantly enhanced its enzymatic activity. This led to a subsequent increase in intracellular NADPH levels, promoting adipogenic differentiation in preadipocyte cells (Kim et al., 2012). We identified eight acetylation sites in malate dehydrogenase (MDH) in this study, including the acetylation of lysines 108, 119, 145, 168, 178, 228, 312, and 314. This protein was also found to be acetylated in S. japonicum (Hong et al., 2016). In our study, other enzymes involved in the TCA cycle were also identified as Kac targets, confirming that acetylation has a potential conserved function in the regulation of the TCA cycle in both eukaryotic and prokaryotic cells.

Protein Metabolism

In addition its regulation of glycometabolism, Kac was also observed in proteins regulating protein expression and proteolysis. In this study, we identified 120 lysine-acetylated proteins related to protein synthesis, accounting for 7.5% of the total lysine-acetylated proteins (Supplementary Table S1). Our results are consistent with those of previous studies showing that ribosomal proteins and aminoacyl–tRNA synthetases are the main targets of acetylation in both prokaryotes and eukaryotes. In E. coli, 22% of the acetylated proteins are involved in the translational machinery or processes (Zhang et al., 2013), while this figure is 15% in Toxoplasma (Jeffers and Sullivan, 2012). Among these proteins, elongation factor 1-α of T. spiralis was heavily acetylated at 15 lysine residues, similar to the acetylation of its orthologs in other eukaryotes (humans, rats, and Toxoplasma) (Jeffers and Sullivan, 2012; Lundby et al., 2012; Shao et al., 2012). tRNA synthetases represent potential drug targets for the prevention and control of parasitic infections. For example, inhibitors of isoleucine–tRNA synthetase and seryl–tRNA synthetase can suppress P. falciparum growth (Miao et al., 2013). In this study, 20 aminoacyl–tRNA synthetases were found to have acetylated lysine residues, reflecting a higher frequency of acetylation than among some other eukaryotes.

The ubiquitin–proteasome system provides a major mechanism for regulating the degradation of cell-cycle regulators and potentially toxic misfolded proteins, and it is required for maintaining cellular protein homeostasis in eukaryotic cells. The core of this large protease is the 20S proteasome, a barrel-shaped structure made of a stack of four heptameric rings. Kac has dual roles in regulating protein degradation in the proteasome. A direct role of protein acetylation involves competing with protein ubiquitination for lysine sites, preventing the proteasomal degradation of target proteins. In addition, in some cases, the acetylation of a specific lysine can unexpectedly create a high-affinity binding site for other proteins, leading to its subsequent ubiquitination and degradation (Caron et al., 2005). In our study, the 20S core particle of the proteasome was frequently modified by Kac (eight proteins), particularly the α subunits (Supplementary Figure S1). Accordingly, three types of amino acid degradation-related enzymes associated with the degradation of valine, leucine, and isoleucine were acetylated (Supplementary Table S1), suggesting that acetylation plays a direct regulatory role in the ubiquitin–proteasome system and influences proteolytic processes in T. spiralis.

Phagocytosis, a multistep process that is completely dependent on cytoskeletal complexes (Radulovic and Godovac-Zimmermann, 2011; Drazic et al., 2016), has been shown to play crucial roles in the acquisition of nutrients, cell proliferation, and virulence in parasites (Mansuri et al., 2014). After binding to a cell surface receptor, the actin cytoskeleton is locally remodeled to provide the necessary force for the formation of phagocytic cups and phagosomes. The spatial and temporal regulation of actin dynamics is required for controlling the various stages of phagocytosis, and interference with actin dynamics reduces phagocytosis (Mansuri et al., 2014). Kac have been suggested to play a crucial role in maintaining actin structure and dynamics and facilitating stress fiber formation and actomyosin interactions (Choudhary et al., 2009). Another group of important cytoskeleton-associated proteins are microtubules, which are formed from α- and β-tubulin heterodimers. The acetylation of α-tubulin has been shown to affect the kinetics of tubulin polymerization and the localization of the motor protein kinesin. The hyperacetylation of microtubules causes the misregulation of transport by kinesin motors (Yang and Seto, 2008; Friedmann and Marmorstein, 2013). In addition, investigations into the antitumor efficacies of histone acetyltransferase inhibitors have shown that the acetylation level of tubulin may be a predictive biomarker of its sensitivity to Histone acetyltransferase inhibition (Di Martile et al., 2016). In our study, the cytoskeleton-associated proteins actin and tubulin, which are also involved in phagocytosis, were found to be acetylated (Supplementary Figure S2), indicating that the acetylation of cytoskeleton-associated proteins may participate in the regulation of phagocytosis in T. spiralis.

Recently, it has become evident that two major endoplasmic reticulum (ER) proteins, calnexin and calreticulin, are crucial regulators of particle uptake into phagosomes, controlling the opening and closing of the cups in Dictyostelium discoideum (Muller-Taubenberger et al., 2001). During phagocytosis, an acidic phagosomal environment generated by V-ATPases is of pivotal importance for the degradation of macromolecules. Unexpectedly, we found that these proteins important in phagocytosis were heavily acetylated in our study. Among these, eight ATPases were found to be acetylated at multiple sites (Supplementary Figure S2). To the best of our knowledge, this is the first time that abundant proteins involved in phagocytosis have been found to be acetylated in an organism, confirming that Kac is of great significance in the regulation of phagocytosis.

The complete life-cycle of T. spiralis occurs in a single host and can be divided into three main stages: ML, adult worms, and newborn larvae (NBL). These occupy two distinct intracellular niches: the intestinal epithelium and skeletal muscle cells. Among these developmental stages, T. spiralis has developed fascinating strategies for adapting to different environments and successfully establishing infection, such as the switching of parasite surface antigens. Studies have indicated that acetylation can regulate antigenic variation and developmental transitions in parasites, and these are principally determined by strong selective pressures (Coleman et al., 2014; Carranza et al., 2016). Thus, it would be of great interest to investigate the global acetyl-proteome of T. spiralis in other developmental stages and compare the levels of acetylation of the proteomes. This will facilitate the dissection of the role of the lysine acetylome in regulating the developmental transitions of T. spiralis and its establishment of infection under the selective pressures of different living environments.

Conclusion

We have produced the first large-scale acetyl-proteome of T. spiralis, an important foodborne parasite, showing abundant Kac among T. spiralis proteins. The lysine acetylome in the present study provides a good foundation for in-depth studies of the functions of reversible Kac in the growth and development of T. spiralis and other nematodes. Identifying the functions of the acetylation of target proteins may also help to design specific and effective drugs or vaccines to prevent and control trichinellosis.

Author Contributions

ML and YZ designed the research. YY, MT, XB, XiL, and IV carried out the data acquisition, analysis, and interpretation. PZ and WC performed the T. spiralis infection experiments. YY, MT, and XB contributed to the writing of the manuscript. XC and XuL contributed to improvement of the manuscript. All authors contributed to the study concept and critical revision of the manuscript for important intellectual content.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This study was supported by the National Key Research and Development Program of China (2017YFD0501302 and 2016YFD0500707), Guangdong Innovative and Entrepreneurial Research Team Program (No. 2014ZT05S123), the scientific fund from Health and Family Planning Commission of Wuxi (No. Q201640), and the Fundamental Research Funds for the Central Universities (1282050205162430).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2017.02674/full#supplementary-material

Footnotes

  1. ^ http://www.ebi.ac.uk/GOA/
  2. ^ http://www.ebi.ac.uk/interpro/
  3. ^ http://weblogo.threeplusone.com/create.cgi
  4. ^ http://string-db.org/

References

Agosin, M., and Aravena, L. C. (1959). Anaerobic glycolysis in homogenates of Trichinella spiralis larvae. Exp. Parasitol. 8, 10–30. doi: 10.1016/0014-4894(59)90003-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bai, X., Hu, X. X., Liu, X. L., Tang, B., and Liu, M. Y. (2017). Current research of trichinellosis in China. Front. Microbiol. 8:1472. doi: 10.3389/fmicb.2017.01472

CrossRef Full Text | Google Scholar

Bauer-Mehren, A. (2013). Integration of genomic information with biological networks using Cytoscape. Methods Mol. Biol. 1021, 37–61. doi: 10.1007/978-1-62703-450-0_3

PubMed Abstract | CrossRef Full Text | Google Scholar

Campos, E. I., and Reinberg, D. (2009). Histones: annotating chromatin. Annu. Rev. Genet. 43, 559–599. doi: 10.1146/annurev.genet.032608.103928

CrossRef Full Text | Google Scholar

Caron, C., Boyault, C., and Khochbin, S. (2005). Regulatory cross-talk between lysine acetylation and ubiquitination: role in the control of protein stability. Bioessays 27, 408–415. doi: 10.1002/bies.20210

PubMed Abstract | CrossRef Full Text | Google Scholar

Carranza, P. G., Gargantini, P. R., Prucca, C. G., Torri, A., Saura, A., Svard, S., et al. (2016). Specific histone modifications play critical roles in the control of encystation and antigenic variation in the early-branching eukaryote Giardia lamblia. Int. J. Biochem. Cell Biol. 81(Pt A), 32–43. doi: 10.1016/j.biocel.2016.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Choudhary, C., Kumar, C., Gnad, F., Nielsen, M. L., Rehman, M., Walther, T. C., et al. (2009). Lysine acetylation targets protein complexes and co-regulates major cellular functions. Science 325, 834–840. doi: 10.1126/science.1175371

PubMed Abstract | CrossRef Full Text | Google Scholar

Close, P., Creppe, C., Gillard, M., Ladang, A., Chapelle, J. P., Nguyen, L., et al. (2010). The emerging role of lysine acetylation of non-nuclear proteins. Cell. Mol. Life. Sci. 67, 1255–1264. doi: 10.1007/s00018-009-0252-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman, B. I., Skillman, K. M., Jiang, R. H. Y., Childs, L. M., Altenhofen, L. M., Ganter, M., et al. (2014). A Plasmodium falciparum histone deacetylase regulates antigenic variation and gametocyte conversion. Cell Host Microbe 16, 177–186. doi: 10.1016/j.chom.2014.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Martile, M., Desideri, M., De Luca, T., Gabellini, C., Buglioni, S., Eramo, A., et al. (2016). Histone acetyltransferase inhibitor CPTH6 preferentially targets lung cancer stem-like cells. Oncotarget 7, 11332–11348. doi: 10.18632/oncotarget.7238

PubMed Abstract | CrossRef Full Text | Google Scholar

Drazic, A., Myklebust, L. M., Ree, R., and Arnesen, T. (2016). The world of protein acetylation. Biochim. Biophys. Acta 1864, 1372–1401. doi: 10.1016/j.bbapap.2016.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubois, F., Caby, S., Oger, F., Cosseau, C., Capron, M., Grunau, C., et al. (2009). Histone deacetylase inhibitors induce apoptosis, histone hyperacetylation and up-regulation of gene transcription in Schistosoma mansoni. Mol. Biochem. Parasitol. 168, 7–15. doi: 10.1016/j.molbiopara.2009.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, J., Shan, C. L., Kang, H. B., Elf, S., Xie, J. X., Tucker, M., et al. (2014). Tyr phosphorylation of PDP1 toggles recruitment between ACAT1 and SIRT3 to regulate the pyruvate dehydrogenase complex. Mol. Cell 53, 534–548. doi: 10.1016/j.molcel.2013.12.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedmann, D. R., and Marmorstein, R. (2013). Structure and mechanism of non-histone protein acetyltransferase enzymes. FEBS J. 280, 5570–5581. doi: 10.1111/febs.12373

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, F., Liu, X., Wu, X. P., Wang, X. L., Gong, D., Lu, H., et al. (2012). Differential DNA methylation in discrete developmental stages of the parasitic nematode Trichinella spiralis. Genome Biol. 13:R100. doi: 10.1186/gb-2012-13-10-r100

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosh, A. K., Saini, S., Das, S., Mandal, A., Sardar, A. H., Ansari, M. Y., et al. (2017). Glucose-6-phosphate dehydrogenase and Trypanothione reductase interaction protects Leishmania donovani from metalloid mediated oxidative stress. Free Radic. Biol. Med. 106, 10–23. doi: 10.1016/j.freeradbiomed.2017.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosh, A. K., Sardar, A. H., Mandal, A., Saini, S., Abhishek, K., Kumar, A., et al. (2015). Metabolic reconfiguration of the central glucose metabolism: a crucial strategy of Leishmania donovani for its survival during oxidative stress. FASEB J. 29, 2081–2098. doi: 10.1096/fj.14-258624

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldberg, E. (1957). Studies on the intermediary metabolism of Trichinella spiralis. Exp. Parasitol. 6, 367–382. doi: 10.1016/0014-4894(57)90028-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gottstein, B., Pozio, E., and Nockler, K. (2009). Epidemiology, diagnosis, treatment, and control of trichinellosis. Clin. Microbiol. Rev. 22, 127–145. doi: 10.1128/CMR.00026-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Henriksen, P., Wagner, S. A., Weinert, B. T., Sharma, S., Bacinskaja, G., Rehman, M., et al. (2012). Proteome-wide analysis of lysine acetylation suggests its broad regulatory scope in Saccharomyces cerevisiae. Mol. Cell. Proteomics 11, 1510–1522. doi: 10.1074/mcp.M112.017251

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, Y., Cao, X., Han, Q., Yuan, C., Zhang, M., Han, Y., et al. (2016). Proteome-wide analysis of lysine acetylation in adult Schistosoma japonicum worm. J. Proteomics 148, 202–212. doi: 10.1016/j.jprot.2016.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Horton, P., Park, K. J., Obayashi, T., Fujita, N., Harada, H., Adams-Collier, C. J., et al. (2007). WoLF PSORT: protein localization predictor. Nucleic Acids Res. 35, W585–W587. doi: 10.1093/nar/gkm259

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, H. G., Zhu, W. W., Qin, J., Chen, M., Gong, L. Y., Li, L., et al. (2017). Acetylation of PGK1 promotes liver cancer cell proliferation and tumorigenesis. Hepatology 65, 515–528. doi: 10.1002/hep.28887

PubMed Abstract | CrossRef Full Text | Google Scholar

Hua, Y. H., Wu, C. Y., Sargsyan, K., and Lim, C. (2014). Sequence-motif detection of NAD(P)-binding proteins: discovery of a unique antibacterial drug target. Sci. Rep. 4:6471. doi: 10.1038/srep06471

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., and Berger, S. L. (2008). The emerging field of dynamic lysine methylation of non-histone proteins. Curr. Opin. Genet. Dev. 18, 152–158. doi: 10.1016/j.gde.2008.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Hyun, K., Jeon, J., Park, K., and Kim, J. (2017). Writing, erasing and reading histone lysine methylations. Exp. Mol. Med. 49:e324. doi: 10.1038/emm.2017.11

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeffers, V., and Sullivan, W. J. (2012). Lysine acetylation is widespread on proteins of diverse function and localization in the protozoan parasite Toxoplasma gondii. Eukaryot. Cell 11, 735–742. doi: 10.1128/Ec.00088-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Goto, S., Kawashima, S., Okuno, Y., and Hattori, M. (2004). The KEGG resource for deciphering the genome. Nucleic Acids Res. 32, D277–D280. doi: 10.1093/nar/gkh063

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, E. Y., Kim, W. K., Kang, H. J., Kim, J. H., Chung, S. J., Seo, Y. S., et al. (2012). Acetylation of malate dehydrogenase 1 promotes adipogenic differentiation via activating its enzymatic activity. J. Lipid Res. 53, 1864–1876. doi: 10.1194/jlr.M026567

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. C., Sprung, R., Chen, Y., Xu, Y., Ball, H., Pei, J., et al. (2006). Substrate and functional diversity of lysine acetylation revealed by a proteomics survey. Mol. Cell. 23, 607–618. doi: 10.1016/j.molcel.2006.06.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T. T., Liu, M. X., Feng, X., Wang, Z., Das, I., Xu, Y. P., et al. (2014). Glyceraldehyde-3-phosphate dehydrogenase is activated by lysine 254 acetylation in response to glucose signal. J. Biol. Chem. 289, 3775–3785. doi: 10.1074/jbc.M113.531640

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J. Y., Zhang, N. Z., Li, W. H., Li, L., Yan, H. B., Qu, Z. G., et al. (2016). Proteomic analysis of differentially expressed proteins in the three developmental stages of Trichinella spiralis. Vet. Parasitol. 231, 32–38. doi: 10.1016/j.vetpar.2016.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X. L., Song, Y. X., Jiang, N., Wang, J. L., Tang, B., Lu, H. J., et al. (2012). Global gene expression analysis of the zoonotic parasite Trichinella spiralis revealed novel genes in host parasite interaction. PLOS Negl. Trop. Dis. 6:e1794. doi: 10.1371/journal.pntd.0001794

PubMed Abstract | CrossRef Full Text | Google Scholar

Lundby, A., Lage, K., Weinert, B. T., Bekker-Jensen, D. B., Secher, A., Skovgaard, T., et al. (2012). Proteomic analysis of lysine acetylation sites in rat tissues reveals organ specificity and subcellular patterns. Cell Rep. 2, 419–431. doi: 10.1016/j.celrep.2012.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Mansuri, M. S., Bhattacharya, S., and Bhattacharya, A. (2014). A novel alpha kinase EhAK1 phosphorylates actin and regulates phagocytosis in Entamoeba histolytica. PLOS Pathog. 10:e1004411. doi: 10.1371/journal.ppat.1004411

PubMed Abstract | CrossRef Full Text | Google Scholar

Miao, J., Lawrence, M., Jeffers, V., Zhao, F., Parker, D., Ge, Y., et al. (2013). Extensive lysine acetylation occurs in evolutionarily conserved metabolic pathways and parasite-specific functions during Plasmodium falciparum intraerythrocytic development. Mol. Microbiol. 89, 660–675. doi: 10.1111/mmi.12303

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitreva, M., Jasmer, D. P., Zarlenga, D. S., Wang, Z., Abubucker, S., Martin, J., et al. (2011). The draft genome of the parasitic nematode Trichinella spiralis. Nat. Genet. 43, 228–235. doi: 10.1038/ng.769

PubMed Abstract | CrossRef Full Text | Google Scholar

Muller-Taubenberger, A., Lupas, A. N., Li, H. W., Ecke, M., Simmeth, E., and Gerisch, G. (2001). Calreticulin and calnexin in the endoplasmic reticulum are important for phagocytosis. EMBO J. 20, 6772–6782. doi: 10.1093/emboj/20.23.6772

CrossRef Full Text | Google Scholar

Murrell, K. D., and Pozio, E. (2011). Worldwide occurrence and impact of human trichinellosis, 1986-2009. Emerg. Infect. Dis. 17, 2194–2202. doi: 10.3201/eid1712.110896

PubMed Abstract | CrossRef Full Text | Google Scholar

Nie, Z., Zhu, H., Zhou, Y., Wu, C., Liu, Y., Sheng, Q., et al. (2015). Comprehensive profiling of lysine acetylation suggests the widespread function is regulated by protein acetylation in the silkworm, Bombyx mori. Proteomics 15, 3253–3266. doi: 10.1002/pmic.201500001

PubMed Abstract | CrossRef Full Text | Google Scholar

Olia, A. S., Barker, K., McCullough, C. E., Tang, H. Y., Speicher, D. W., Qiu, J., et al. (2015). Non-enzymatic protein acetylation detected by NAPPA protein arrays. Acs Chem. Biol. 10, 2034–2047. doi: 10.1021/acschembio.5b00342

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, J. Y., Ye, Z. C., Cheng, Z. Y., Peng, X. J., Wen, L. Y., and Zhao, F. K. (2014). Systematic analysis of the lysine acetylome in Vibrio parahemolyticus. J. Proteome Res. 13, 3294–3302. doi: 10.1021/pr500133t

PubMed Abstract | CrossRef Full Text | Google Scholar

Radulovic, M., and Godovac-Zimmermann, J. (2011). Proteomic approaches to understanding the role of the cytoskeleton in host-defense mechanisms. Expert Rev. Proteomics 8, 117–126. doi: 10.1586/epr.10.91

PubMed Abstract | CrossRef Full Text | Google Scholar

Shao, J. L., Xu, D., Hu, L. D., Kwan, Y. W., Wang, Y. F., Kong, X. Y., et al. (2012). Systematic analysis of human lysine acetylation proteins and accurate prediction of human lysine acetylation through bi-relative adapted binomial score Bayes feature representation. Mol. Biosyst. 8, 2964–2973. doi: 10.1039/c2mb25251a

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, H., Xu, W., and Lan, F. (2017). Histone lysine demethylases in mammalian embryonic development. Exp. Mol. Med. 49:e325. doi: 10.1038/emm.2017.57

PubMed Abstract | CrossRef Full Text | Google Scholar

Spange, S., Wagner, T., Heinzel, T., and Kramer, O. H. (2009). Acetylation of non-histone proteins modulates cellular signalling at multiple levels. Int. J. Biochem. Cell Biol. 41, 185–198. doi: 10.1016/j.biocel.2008.08.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, T., Li, X., Zhang, P., Chen, W. D., Zhang, H. L., Li, D. D., et al. (2015). Acetylation of Beclin 1 inhibits autophagosome maturation and promotes tumour growth. Nat. Commun. 6:7215. doi: 10.1038/ncomms8215

PubMed Abstract | CrossRef Full Text | Google Scholar

Theobald, D. L., Mitton-Fry, R. M., and Wuttke, D. S. (2003). Nucleic acid recognition by OB-fold proteins. Annu. Rev. Biophys. Biomol. Struct. 32, 115–133. doi: 10.1146/annurev.biophys.32.110601.142506

CrossRef Full Text | Google Scholar

Wagner, G. R., and Hirschey, M. D. (2014). Non-enzymatic protein acylation as a carbon stress regulated by sirtuin deacylases. Mol. Cell. 54, 5–16. doi: 10.1016/j.molcel.2014.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Walsh, C. T., Garneau-Tsodikova, S., and Gatto, G. J. (2005). Protein posttranslational modifications: the chemistry of proteome diversifications. Angew. Chem. Int. Ed. 44, 7342–7372. doi: 10.1002/anie.200501023

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Rosa, B. A., Nare, B., Powell, K., Valente, S., Rotili, D., et al. (2015). Targeting lysine deacetylases (KDACs) in parasites. PLOS Negl. Trop. Dis. 9:e0004026. doi: 10.1371/journal.pntd.0004026

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Zhang, Y., Yang, C., Xiong, H., Lin, Y., Yao, J., et al. (2010). Acetylation of metabolic enzymes coordinates carbon source utilization and metabolic flux. Science 327, 1004–1007. doi: 10.1126/science.1179687

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z. K., Cai, Q., Liu, J., Ying, S. H., and Feng, M. G. (2017). Global insight into lysine acetylation events and their links to biological aspects in Beauveria bassiana, a fungal insect pathogen. Sci. Rep. 7:44360. doi: 10.1038/srep44360

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinert, B. T., Scholz, C., Wagner, S. A., Iesmantavicius, V., Su, D., Daniel, J. A., et al. (2013). Lysine succinylation is a frequently occurring modification in prokaryotes and eukaryotes and extensively overlaps with acetylation. Cell Rep. 4, 842–851. doi: 10.1016/j.celrep.2013.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Witze, E. S., Old, W. M., Resing, K. A., and Ahn, N. G. (2007). Mapping protein post-translational modifications with mass spectrometry. Nat. Methods 4, 798–806. doi: 10.1038/nmeth1100

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, L., Wang, X., Zeng, J., Zhou, M., Duan, X., Li, Q., et al. (2015). Proteome-wide lysine acetylation profiling of the human pathogen Mycobacterium tuberculosis. Int. J. Biochem. Cell Biol. 59, 193–202. doi: 10.1016/j.biocel.2014.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong, Y., and Guan, K. L. (2012). Mechanistic insights into the regulation of metabolic enzymes by acetylation. J. Cell Biol. 198, 155–164. doi: 10.1083/jcb.201202056

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J. M., Qiu, C. H., Xia, Y. X., Yao, L. X., Fu, Z. Q., Yuan, C. X., et al. (2010). Molecular cloning and functional characterization of Schistosoma japonicum enolase which is highly expressed at the schistosomulum stage. Parasitol. Res. 107, 667–677. doi: 10.1007/s00436-010-1913-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X. J., and Seto, E. (2008). Lysine acetylation: codified crosstalk with other posttranslational modifications. Mol. Cell. 31, 449–461. doi: 10.1016/j.molcel.2008.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Youn, H. D. (2017). Methylation and demethylation of DNA and histones in chromatin: the most complicated epigenetic marker. Exp. Mol. Med. 49:e321. doi: 10.1038/emm.2017.38

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, K., Zheng, S. Z., Yang, J. S., Chen, Y., and Cheng, Z. Y. (2013). Comprehensive profiling of protein lysine acetylation in Escherichia coli. J. Proteome Res. 12, 844–851. doi: 10.1021/pr300912q

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, S., Xu, W., Jiang, W., Yu, W., Lin, Y., Zhang, T., et al. (2010). Regulation of cellular metabolism by protein lysine acetylation. Science 327, 1000–1004. doi: 10.1126/science.1179689

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Qian, G., Yi, X., Li, X., and Liu, W. (2016). Systematic analysis of the lysine acetylome in Candida albicans. J. Proteome Res. 15, 2525–2536. doi: 10.1021/acs.jproteome.6b00052

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: lysine acetylation, post-translational modification, lysine acetylation motif, interaction network, Trichinella spiralis

Citation: Yang Y, Tong M, Bai X, Liu X, Cai X, Luo X, Zhang P, Cai W, Vallée I, Zhou Y and Liu M (2018) Comprehensive Proteomic Analysis of Lysine Acetylation in the Foodborne Pathogen Trichinella spiralis. Front. Microbiol. 8:2674. doi: 10.3389/fmicb.2017.02674

Received: 01 August 2017; Accepted: 21 December 2017;
Published: 11 January 2018.

Edited by:

Maria Schirone, Università di Teramo, Italy

Reviewed by:

Zhao Chen, University of California, Davis, United States
M. Guadalupe Ortega-Pierres, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional (CINVESTAV-IPN), Mexico

Copyright © 2018 Yang, Tong, Bai, Liu, Cai, Luo, Zhang, Cai, Vallée, Zhou and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mingyuan Liu, liumy36@163.com Yonghua Zhou, toxo2001@163.com

These authors have contributed equally to this work.