- Open Access
First comprehensive proteome analysis of lysine crotonylation in Streptococcus agalactiae, a pathogen causing meningoencephalitis in teleosts
Proteome Science volume 19, Article number: 14 (2021)
Streptococcus agalactiae is a common colonizer of the rectovaginal tract and lead to infectious diseases of neonatal and non-pregnant adults, which also causes infectious disease in fish and a zoonotic risk as well. Lysine crotonylation (Kcr) is a kind of histone post-translational modifications discovered in 2011. In yeast and mammals, Kcr function as potential enhancers and promote gene expression. However, lysine crotonylation in S. agalactiae has not been studied yet.
In this study, the crotonylation profiling of fish pathogen, S. agalactiae was investigated by combining affinity enrichment with LC MS/MS. The Kcr modification of several selected proteins were further validated by Western blotting.
In the present study, we conducted the proteome-wide profiling of Kcr in S. agalactiae and identified 241 Kcr sites from 675 screened proteins for the first time. Bioinformatics analysis showed that 164 sequences were matched to a total of six definitively conserved motifs, and many of them were significantly enriched in metabolic processes, cellular process, and single-organism processes. Moreover, four crotonylation modified proteins were predicted as virulence factors or to being part of the quorum sensing system PTMs on bacteria. The data are available via ProteomeXchange with identifier PXD026445.
These data provide a promising starting point for further functional research of crotonylation in bacterial virulence in S. agalactiae.
Lysine crotonylation is a newly discovered post-translational modification, which is structurally and functionally different from the widely studied lysine acetylation . Crotonylation occurs on the ε-amino group of lysine but differ from acetylation in its four-carbon length and planar orientation . In order to characterize the global crotonylation proteome, a proteomic method based on sensitive immune-affinity purification and high-resolution liquid chromatography-tandem mass spectrometry (LC-MS/MS) was applied to identify new crotonylated proteins and modification sites. With the utilization of high specificity antibodies, increasing crotonylated lysine residues and crotonylated proteins were identified so far. Also, the involvement of crotonylated proteins in multiple cellular processes were well-recorded .
In bacterial cells, lysine acylation (Kace) modification has been well studied in many species, including Vibrio alginolyticus, Vibrio parahaemolyticus, Escherichia coli and Mycobacterium tuberculosis [4,5,6,7]. The Kace modification was involved in many basal physiological functions such as chemotaxis, virulence, and antibiotic resistance [8,9,10]. However, the crotonylome profile of bacterial pathogen of fish is still unclear. Streptococcus agalactiae or group B streptococcus (GBS) is an important bacterial in fish that result in abnormal swimming behaviour, popeye, haemorrhages, and severe mortalities. GBS is also an important zoonotic pathogen due to the fact that it can infect human as well . Recently, several proteins of GBS such as LuxS  and RbsB  participated in AI-2-mediated quorum sensing system. However, the intrinsic quorum sensing mechanisms are still largely unknown, especially the roles of Kac modification.
In present study, we investigated the global lysine crotonylation proteome of S. agalactiae using high-resolution LC-MS/MS coupled with highly sensitive immune-affinity purification. We conducted proteome-wide profiling of Kcr in S. agalactiae and identified 241 Kcr sites from 675 proteins, indicating the Kcr exist in S. agalactiae. Bioinformatics analysis showed that 164 sequences were matched to a total of six definitively conserved motifs, and many of them were significantly enriched in metabolic processes, cellular process, and single-organism processes. Moreover, four contonylation modified proteins was associated with bacterial QS system and virulence. Our data expand the understanding of Kcr modification profiling in prokaryotes and provides new insights into the potential role of PTMs in fish pathogens.
Bacterial materials and growth conditions
S.agalactiae ZQ0910, a virulent strain isolated from intensive tilapia farm with a typical streptococcicosis outbreak in the Guangdong province of China in 2009 . It belongs to serotype Ia and used for challenge study. The strain was grown aerobically overnight at 28 °C in a shaker bath, and then overnight cultured cells were diluted into 1:100 in BHI medium. Bacteria were then harvested by centrifugation (5000 g, 5 min). Cell pellets were washed and resuspended in PBS to obtain the desired concentrations for the infection doses. The final bacterial concentrations were confirmed by plating ten-fold serial dilutions onto blood-agar plates.
Sample was sonicated three times on ice using a high intensity ultrasonic processor (Scientz) in lysis buffer (8 M urea, 3 μM TSA and 50 mM NAM and 1% Protease Inhibitor Cocktail, Millipore). The remaining debris was removed by centrifugation at 12,000 g at 4 °C for 10 min. Finally, the supernatant was collected and the protein concentration was determined with BCA kit according to the manufacturer’s instructions.
For digestion, the protein solution was reduced with 5 mM dithiothreitol (final concentration) for 30 min at 56 °C and alkylated with 11 mM iodoacetamide (final concentration) for 45 min at room temperature in darkness. The protein sample was then diluted by adding 100 mM TEAB to urea concentration less than 2 M. Finally, trypsin was added at 1:50 trypsin-to-protein mass ratio for the first digestion overnight and 1:100 trypsin-to-protein mass ratio for a second 4 h-digestion .
To enrich Kcro peptides, tryptic peptides dissolved in NETN buffer (100 mM NaCl, 1 mM EDTA, 50 mM Tris-HCl, 0.5% NP-40, pH 8.0) were incubated with pre-washed antibody beads (PTM Biolabs) at 4 °C overnight with gentle shaking. Then the beads were washed four times with NETN buffer and twice with H2O. The bound peptides were eluted from the beads with 0.1% trifluoroacetic acid. Finally, 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, followed by LC-MS/MS analysis.
Quantitative proteomic analysis by LC-MS/MS
The tryptic peptides were dissolved in 0.1% formic acid (solvent A), directly loaded onto a reversed-phase pre-column (Acclaim PepMap 100; Thermo Fisher Scientific, Inc.). Peptide separation was performed using a reversed-phase analytical column (Acclaim PepMap RSLC; Thermo Fisher Scientific, Inc.) with gradient that comprised of an increase from 6 to 23% solvent B (0.1% formic acid in 98% acetonitrile) over 26 min, 23 to 35% in 8 min and climbing to 80% in 3 min then holding at 80% for the last 3 min, all at a constant flow rate of 400 nL/min on an EASY-nLC 1000 UPLC system .
The peptides were subjected to nanospray ion source (NSI) source followed by tandem mass spectrometry (MS/MS) in Q ExactiveTM Plus (Thermo) coupled online to the 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 in the Orbitrap at a resolution of 70,000. Peptides were then selected for MS/MS using normalized collision energy (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 s dynamic exclusion. The electrospray voltage applied was 2.0 kV. Automatic gain control (AGC) was used to prevent overfilling of the Qrbitrap; 50,000 ions were accunulated to generate MS/MS spectra. For MS scans, the m/z scan range was 350 to 1800. The fixed first mass was set as 100 m/z.
The resulting MS/MS data were processed using Maxquant search engine (v.22.214.171.124). Tandem mass spectra were searched against 1123 database concatenated with reverse decoy database. Trypsin/P was specified as cleavage enzyme allowing up to 4 missing cleavages. The mass tolerance for precursor ions was set as 20 ppm in First search and 5 ppm in Main search, and the mass tolerance for fragment ions was set as 0.02 Da. Carbamidomethyl on Cys was specified as fixed modification and crotonylation on lysine were specified as variable modifications. FDR was adjusted to < 1% and minimum score for modified peptides was set > 40. For selected specific Kcr sites, site localization probability was set to > 0.75. All other parameters in MaxQuant were used as default.
Bioinformatics analysis for gene ontology annotation
The Gene Ontology, or GO, is a major bioinformatics initiative to unify the representation of gene and gene product attributes across all species. Gene Ontology (GO) annotation proteome was derived from the UniProt-GOA database (www. http://www.ebi.ac.uk/GOA/) . Firstly, Converting identified protein ID to UniProt ID and then mapping to GO IDs by protein ID. If some identified proteins were not annotated by UniProt-GOA database, the InterProScan soft would be used to annotated protein’s GO functional based on protein sequence alignment method. Then proteins were classified by Gene Ontology annotation based on three categories: biological process, cellular component and molecular function.
For subcelluar localization, we used wolfpsort a subcellular localization predication soft to predict subcellular localization. Wolfpsort is an updated version of PSORT/PSORT II for the prediction of eukaryotic sequences. Special for protokaryon species, Subcellular localization prediction soft CELLO was used.
Kyoto Encyclopedia of Genes and Genomes (KEGG ) database was used to annotate protein pathway. Firstly, using KEGG online service tools KAAS to annotated protein’s KEGG database description. Then mapping the annotation result on the KEGG pathway database using KEGG online service tools KEGG mapper.
Bioinformatics analysis for enrichment GO and KEGG pathway analysis
Proteins were classified by GO annotation into three categories: biological process, cellular compartment and molecular function. For each category, a two-tailed Fisher’s exact test was employed to test the enrichment of the identified modified protein against all proteins of the species database. The GO with a corrected p-value < 0.05 is considered significant.
Encyclopedia of Genes and Genomes (KEGG) database was used to identify enriched pathways by a two-tailed Fisher’s exact test to test the enrichment of the identified modified protein against all proteins of the species database. The pathway with a corrected p-value < 0.05 was considered significant. These pathways were classified into hierarchical categories according to the KEGG website.
Soft MoMo (motif-x algorithm) was used to analysis the model of sequences constituted with amino acids in specific positions of modify-21-mers (10 amino acids upstream and downstream of the site, but phosphorylation with modify-13-mers that 6 amino acids upstream and downstream of the site) in all protein sequences. And all the database protein sequences were used as background database parameter. Minimum number of occurrences was set to 20. Emulate original motif-x was ticked, and other parameters with default.
For further hierarchical clustering based on differentially modified protein functional classification (such as: GO, Domain, Pathway, Complex). We first collated all the categories obtained after enrichment along with their P values, and then filtered for those categories which were at least enriched in one of the clusters with P value < 0.05. This filtered P value matrix was transformed by the function x = −log10 (P value). Finally these x values were z-transformed for each functional category. These z scores were then clustered by one-way hierarchical clustering (Euclidean distance, average linkage clustering) in Genesis. Cluster membership were visualized by a heat map using the “heatmap.2” function from the “gplots” R-package.
Protein-protein interaction network
All differentially expressed modified protein database accession or sequence were searched against the STRING database version 10.5 for protein-protein interactions . Only interactions between the proteins belonging to the searched data set were selected, thereby excluding external candidates. STRING defines a metric called “confidence score” to define interaction confidence; we fetched all interactions that had a confidence score > 0.7 (high confidence). Interaction network form STRING was visualized in R package “networkD3”.
Proteins were run 12% SDS-PAGE gels, then adjust the voltage to the 100 V and electrophoresis for 15 min. After the bromophenol blue migrated out of the concentrated gel, the voltage was switched to 200 V for 40 min. And transferred to a polyvinylidene fluoride (PVDF) membrane by semi-dry technique and set the power to 100 V (constant voltage) for 1 h at 4 °C. The membranes were blocked in TBS (Tris buffered saline) containing 0.05% (v/v) Tween 20 with 5% (w/v) skim milk and incubated 1 h at room temperature. The primary antibodies used in the western blot were anti-Kcr (1:2000), anti-YebC (1:5000), anti-CsbD (1:4000), anti-Sal (1:4000), anti-Coldshock (1:4000) and incubated overnight at 4 °C. After washing in 10 mM PBS (pH 7.4) containing 0.1% Tween-20, the membranes were incubated in horseradish peroxidase (HRP) conjugated goat anti-rat IgG (CWBIO, China) diluted 1:5000 at room temperature for 3 h. Finally, the membrane was visualized using the ECL system (Bio-Rad, Hercules, CA, USA), and recorded by the ChemiDoc™ MP (Bio-Rad, Hercules, CA, USA) imaging system .
Results and discussion
Detection of lysine-crotonylated proteins in S. agalactiae
To obtain the global crotonylation proteome of S. agalactiae, a proteomic method based on sensitive immune-affinity purification and high-resolution LC-MS/MS was applied to identify crotonylated proteins and their modification sites in S. agalactiae. An overview of the experimental procedures is shown in Fig. 1. A total of 241 lysine crotonylation sites distributed in 675 proteins were identified. Detailed information for all identified crotonylated peptides and their corresponding proteins was shown in (Fig. 2a, Supplementary Table S1). Among the 675 crotonylated proteins, 145 (21%) proteins contained one or two crotonylation sites, and 13 (1%) proteins had 7 or more crotonylation sites (Fig. 2b). Most peptides ranged from 7 to 15 amino acids in length (Fig. 2c). For quality control validation of the MS data, we evaluated the mass error of all identified peptides. The mass error of most crotonylated peptides ranged from − 5 to 5 ppm, indicating an expected error control from the MS dataset (Fig. 2d).
Motifs and secondary structures of lysine crotonylated peptides
To evaluate the nature of the crotonylated lysines in S. agalactiae, the sequence motifis in all identified crotonylated peptides were investigated using the Motif-X programme . Of all surrounding sequences, 164 sequences were matched to a total of two definitively conserved motifs, there motifs were K***Akcr, AKkcr shown in Fig. 3a. These motifs are likely to represent a feature of crotonylation in S. agalactiae. Hierarchical cluster analysis was also performed to further analyze these motifs. As shown in the heat map (Fig. 3b), the enrichment of negatively charged A residues was observed in the − 10 to − 6, − 4 to + 4 and − 6 to + 10 position, while negatively charged residues E was markedly enriched in the − 6 to + 6 position, the residues K was markedly enriched in the − 10 to − 6 and + 3 to + 9 position.
To explore the relationship between lysine crotonylation and protein secondary structures, a structural analysis of all crotonylated proteins was performed using the algorithm NetSurfP. As shown in Fig. 3c, approximately 28% of the crotonylated sites were located in ɑ-helices, and 7% of the sites were located in β-strands. The remaining 65% of the crotonylated sites were located in disordered coils. The surface accessibility of the crotonylated lysine sites was also evaluated. The results showed that 41.12% of the crotonylated lysine sites were exposed to the protein surface, close to that of all lysine residues (Fig. 3d). Therefore, lysine crotonylation likely does not affect the surface properties of modified proteins.
Functional annotation and subcellular localization of crotonylated proteins
To obtain an overview of crotonylated proteins in S.agalactiae, the Gene Ontology (GO) functional classification of all crotonylated proteins based on their biological processes, cellular component, molecular function and subcellular locations was investigated (Supplementary Table S2). Within the biological processes category, the majority of crotonylated proteins were associated with metabolic processes, cellular process, and single-organism processes, respectively accounting for 36, 32 and 19% of all the crotonylated proteins (Fig. 4a). For the cellular component category, 40, 26 and 22% of the crotonylated proteins were associated with cell, macromolecular complex and organelle, respectively (Fig. 4b). In the molecular function category, the majority of catalytic activity, binding and structural molecule activity, respectively accounting for 40, 34 and 16% of all the crotonylated proteins (Fig. 4c). Subcellular localization analysis revealed that most of the crotonylated proteins were localized in cytoplasmic (79%), membrane (11%) and extracellular (10%) (Fig. 4d).
Functional enrichment analysis
To better comprehend the biological function of these crotonylated proteins, we performed an enrichment analysis of the GO (Supplementary Table S3), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway (Supplementary Table S4), and Pfam domain databases (Supplementary S5). The enrichment analysis of the biological process revealed that the crotonylated proteins were significantly enriched in the peptide biosynthetic process, peptide metabolic process, amide biosynthetic process, cellular protein metabolic process and cellular macromolecule biosynthetic process (Fig. 5a). Based on the enrichment results of the cellular component category, most crotonylated proteins were related to cell, cell part, intracellular and intracellular part (Fig. 5b). In the molecular function category, structural molecule activity, glycerone kinase activity and RNA binding (Fig. 5c). The KEGG pathway enrichment analysis showed that a majority of the crotonylated proteins were ribosome-related proteins, Glycolysis/Gluconeogenesis, methane metabolism, pyruvate metabolism, propanoate metabolism, citrate cycle (TCA cycle) and fructose and mannose metabolism (Fig. 5d, Supplement Figs. 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10). Consistent with these observations, Pfam domains, including the Aldehyde/histidinol dehydrogenase, Aldehyde dehudrogenase (C-terminal and N-terminal), Aldehyde dehydrogenase domain and translation protein (beta-barrel domain) (Fig. 5d), were significantly found in crotonylated proteins, implying an important role for lysine crotonylation in these processes.
Protein interaction network of the crotonylated proteins in S. agalactiae
To further identify the cellular processes regulated through crotonylation in S. agalactiae, the crotonylated protein interaction network was established using an algorithm in Cytoscape software. A total of 193 crotonylated proteins were mapped to the protein interaction database (Supplementary Table S4). As shown in Fig. 6, ribosomal, arginine biosynthesis, fructose and mannose metabolism, Glycolysis/Gluconeogenesis and RNA degradation proteins were highly interconnected, indicating that these translation process may be regulated by crotonylated in S. agalactiae. The physiological interactions among these crotonylated protein complexes may contribute to their cooperation and coordination in S. agalactiae. The results of PPI networks suggest that lysine crotonylation probably play a role in regulating post-translational modification of protein in S. agalactiae, which contributes to cooperation and coordination of metabolic pathways.
Validation of ten lysine-crotonylated proteins using Western blotting
To further validate the identified lysine-crotonylated results, four Kcr proteins (YebC/PmpR family DNA-binding transcriptional regulator , surface antigen-like protein , cold-shock protein , diacylglycerol kinase, CsbD ) were selected and analyzed by western blotting. The four proteins were identified with anti-crotonylation and anti-target protein antibody, respectively (Fig. 7). Consistent with lysine-crotonylated proteomic data, the all four detected proteins exhibited crotonylation modifications.
In this study, using high-resolution LC-MS/MS coupled with highly sensitive immune-affinity purificationwe, we investigated the global lysine crotonylation proteome of S. agalactiae. We reported the proteome-wide profiling of Kcr in S. agalactiae and identified 241 Kcr sites of 675 proteins for the first time. The bio-informatics analysis showed that 164 sequences were matched to a total of six definitively conserved motifs, and many of them were enriched in metabolic processes, cellular process, and single-organism processes. Moreover, we also found a considerable four crotonylation modified proteins involve in sensing system and bacterial virulence, potentially providing targets for vaccine development. The present data provides new insights into the pathogenesis of S. agalactiae from fish, and the virulence mechanisms should be further investigated.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Wan J, Liu H, Chu J, Zhang H. Functions and mechanisms of lysine crotonylation. J Cell Mol Med. 2019;23(11):7163–9.
Sun H, Liu X, Li F, Li W, Zhang J, Xiao Z, et al. First comprehensive proteome analysis of lysine crotonylation in seedling leaves of Nicotiana tabacum. Sci Rep. 2017;7(1):3013.
Wei W, Mao A, Tang B, Zeng Q, Gao S, Lu L, et al. Large-scale identification of protein Crotonylation reveals its role in multiple cellular functions. J Proteome Res. 2017;16(4):1743–52.
Pang H, Li W, Zhang W, et al. Acetylome profiling of Vibrio alginolyticus reveals its role in bacterial virulence. J Proteom. 2019;211:103543.
Pan J, Ye Z, Cheng Z, Peng X, Wen L, Zhao FJ. Systematic analysis of the lysine acetylome in Vibrio parahemolyticus. J Proteome Res. 2014;13(7):3294–302.
Zhang K, Zheng S, Yang JS, Chen Y, Cheng ZJ. Comprehensive Profiling of Protein Lysine Acetylation in Escherichia Coli. J Proteome Res. 2013;12(2):844.
Liu F, Yang M, Wang X, Yang S, Gu J, Zhou J, et al. Acetylome analysis reveals diverse functions of lysine acetylation in Mycobacterium tuberculosis; 2014.
Yan J, Barak R, Liarzi O, Shainskaya A, Eisenbach MJ. In Vivo Acetylation of CheY, a Response Regulator in Chemotaxis of Escherichia coli. J Mol Biol. 2008;376(5):1260–71.
Ren J, Sang Y, Lu J, Yao YF. Protein Acetylation and Its Role in Bacterial Virulence; 2017. p. 768–79.
Liu L, Wang G, Song L, Lv B, Liang W. Acetylome analysis reveals the involvement of lysine acetylation in biosynthesis of antibiotics in Bacillus amyloliquefaciens. Sci Rep. 2016;6:20108.
Delannoy CMJ, Crumlish M, Fontaine MC, Pollock J, Foster G, Dagleish MP, et al. Human Streptococcus agalactiae strains in aquatic mammals and fish. BMC Microbiol. 2013;13(1):41.
Balestrino D, Haagensen JAJ, Rich C, Forestier CJ. Characterization of Type 2 Quorum Sensing in Klebsiella pneumoniae and Relationship with Biofilm Formation. J Bacteriol. 2005;187(8):2870–80.
Demuth DR, Lamont RJ, James D, Shao HJI. Immunity: The Actinobacillus actinomycetemcomitans ribose binding protein RbsB interacts with cognate and heterologous autoinducer 2 signals. Infect Immun. 2006;74(7):4021–9.
Wang B, Jian J, Lu Y, Cai S, Huang Y, Tang J, et al. Complete genome sequence of Streptococcus agalactiae ZQ0910, a pathogen causing meningoencephalitis in the GIFT strain of Nile tilapia (Oreochromis niloticus). J Bacteriol. 2012;194(18):5132–3.
Chen Z, Luo L, Chen R, Hu H, Pan Y, Jiang H, et al. Acetylome profiling reveals extensive lysine acetylation of the fatty acid metabolism pathway in the DiatomPhaeodactylum tricornutum; 2017.
Pan Z, Wang M, Ye Z, Zhang S, Xu XJOL. Global analysis of histone lysine acetylation and proteomic changes in EC109 cells treated with the histone deacetylase inhibitor FK228; 2018.
Li W, Yao Z, Zhang X, Huang F, Lin W, Li XJ. Biotechnology: Global protein expression profile response of planktonic Aeromonas hydrophila exposed to chlortetracycline; 2017.
Kohl M, Wiese S, Warscheid BJ. Cytoscape: software for visualization and analysis of biological networks. Methods Mol Biol. 2011;696:291–303.
Wang Y, Wang X, Ali F, Li Z, Fu Y, Yang X, Lin W, Lin X. Comparative Extracellular Proteomics of Aeromonas hydrophila Reveals Iron-Regulated Secreted Proteins as Potential Vaccine Candidates. Front Immunol. 2019;10:256. https://doi.org/10.3389/fimmu.2019.00256. PMID: 30833947; PMCID: PMC6387970.
Chou MF, Schwartz D. Biological sequence motif discovery using motif-x.. Curr Protoc Bioinform. 2011;35(1).
Wei L, Wu Y, Qiao H, Xu W, Zhang Y, Liu X, Wang Q. YebC controls virulence by activating T3SS gene expression in the pathogen Edwardsiella piscicida. FEMS Microbiol Lett. 2018;365(14). https://doi.org/10.1093/femsle/fny137. PMID: 29901702.
Goel S, Muthusamy A, Miao J, Cui L, Salanti A, Winzeler EA, et al. Targeted disruption of a ring-infected erythrocyte surface antigen (resa)-like export protein gene in plasmodium falciparum confers stable chondroitin 4-sulfate cytoadherence capacity. J Biol Chem. 2014;289(49):34408–21.
Petrosian SA, Makhatadze GI. Contribution of proton linkage to the thermodynamic stability of the major coldshock protein of Escherichia coli CspA. Protein Sci. 2000;9(2):387-94. https://doi.org/10.1110/ps.9.2.387. PMID: 10716191; PMCID: PMC2144560.
Mérida I, ávila-Flores A, Merino E. Diacylglycerol kinases: at the hub of cell signalling. Biochem J. 2008;409(1):1–18.
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD026445.
This study was supported by the National Natural Science Foundation of China (Grant No. 32002426, U20A2065, 32073006). The Natural Science Foundation of Guangdong Province, China (2019A1515011875), International Science and Technology Cooperation Project of Science and Technology Planning Project of Guangdong Province, China (2019A050510044).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Guangdong Ocean University.
Consent for publication
The submission is not under review at any other publication, in whole or in part, and all the authors listed have approved the enclosed manuscript. All the authors declare no financial conflict of interest and have fulfilled the criteria of authorship for the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, X., Fan, B., Fan, C. et al. First comprehensive proteome analysis of lysine crotonylation in Streptococcus agalactiae, a pathogen causing meningoencephalitis in teleosts. Proteome Sci 19, 14 (2021). https://doi.org/10.1186/s12953-021-00182-y
- Streptococcus agalactiae
- Virulence factors