- Research article
- Open Access
Label-free quantitative proteomics of Corynebacterium pseudotuberculosis isolates reveals differences between Biovars ovis and equi strains
BMC Genomics volume 18, Article number: 451 (2017)
Corynebacterium pseudotuberculosis is a pathogen classified into two biovars: C. pseudotuberculosis biovar ovis, the etiologic agent of caseous lymphadenitis and C. pseudotuberculosis biovar equi, which causes ulcerative lymphangitis. The available whole genome sequences of different C. pseudotuberculosis strains have enabled identify difference of genes related both virulence and physiology of each biovar. To evaluate be this difference could reflect at proteomic level and to better understand the shared factors and the exclusive ones of biovar ovis and biovar equi strains, we applied the label-free quantitative proteomic to characterize the proteome of the strains: 1002_ovis and 258_equi, isolated from goat (Brazil) and equine (Belgium), respectively.
From this analysis, we characterized a total of 1230 proteins in 1002_ovis and 1220 in 258_equi with high confidence. Moreover, the core-proteome between 1002_ovis and 258_equi obtained here is composed of 1122 proteins involved in different cellular processes, which could be necessary for the free living of C. pseudotuberculosis. In addition, 120 proteins from this core-proteome presented change in abundant with statistically significant differences. Considering the exclusive proteome, we detected strain-specific proteins to each strain. When correlated, the exclusive proteome of each strain and proteome with change in abundant, the proteomic differences, between the 1002_ovis and 258_equi, this related to proteins involved in cellular metabolism, information storage and processing, cellular processes and signaling.
This study reports the first comparative proteomic study of the biovars ovis and equi of C. pseudotuberculosis. The results generated in this study provide information about factors which can contribute to understanding both the physiology and the virulence of this pathogen.
Corynebacterium pseudotuberculosis is a Gram-positive facultative intracellular pathogen of the Corynebacterium, Mycobacterium, Nocardia, and Rhodococcus (CMNR) group. The CMNR group of pathogens has high G + C content in their genomes and shows a specific cell wall organization composed of peptidoglycan, arabinogalactan, and mycolic acids . C. pseudotuberculosis is subdivided into two biovars: (i) C. pseudotuberculosis biovar ovis (nitrate negative) which is the etiologic agent of caseous lymphadenitis in small ruminants  and mastitis in dairy cattle  and (ii) C. pseudotuberculosis biovar equi (nitrate positive) that causes ulcerative lymphangitis and abscesses in internal organs of equines  and oedematous skin disease in buffalos . C. pseudotuberculosis infection is reported worldwide and causes significant economic losses by affecting wool, meat, and milk production [6,7,8,9].
Various studies at genome level have been carried out by our research group in order to explore the molecular basis of specific and shared factors among different strains of C. pseudotuberculosis that could contribute to such biovar specific pathogenicity. Our studies on whole-genome sequencing and analysis of several C. pseudotuberculosis strains belonging to biovar ovis and equi, isolated from different hosts showed an average genome size of approximately 2,3 Mb, a core-genome having approximately 1504 genes across several C. pseudotuberculosis species, and accessory genomes of biovar equi and ovis composed of 95 and 314 genes, respectively [10,11,12]. According with pan-genome analysis, C. pseudotuberculosis biovar ovis presented a more clonal-like behavior, than the C. pseudotuberculosis biovar equi. In addition, in this in silico study was observed a variability most interesting related to pilus genes, where biovar ovis strain presented high similarity, while, biovar equi strains have a great variability, suggesting that this variability could influence in the adhesion and invasion cellular of each biovar .
Apart from the structural genome informatics studies of C. pseudotuberculosis, some proteomic studies were conducted to explore the functional genome of this pathogen [13,14,15,16,17,18,19]. However, all these proteomic studies were performed using only strains belonging to biovar ovis. Until the present time, no proteomic studies were performed between biovar equi strains or between biovar ovis and biovar equi strains. Therefore, to provide insights on shared and exclusive proteins among biovar ovis and biovar equi strains and to complement the previous studies on functional and structural genomics of C. pseudotuberculosis biovars, using LC-MSE approach [13, 18] this study reports for the first time a comparative proteomic analysis of two C. pseudotuberculosis strains, 1002_ovis and 258_equi, isolated from caprine (Brazil) and equine (Belgium), respectively. Our proteomic dataset promoted the validations of previous work in silico of C. pseudotuberculosis; in addition, the qualitative and quantitative differences in the proteins identified in this present work have potential to help understand the factors that might contribute for pathogenic process of biovar ovis and equi strains.
Bacterial strain and growth condition
C. pseudotuberculosis biovar ovis 1002, isolated from a goat in Brazil, and C. pseudotuberculosis biovar equi 258, isolated from a horse in Belgium, were maintained in brain–heart infusion broth or agar (1.5%) (BHI-HiMedia Laboratories Pvt. Ltd., India) at 37 °C. For proteomic analysis, overnight cultures (three biological replicate to each strain) in BHI were inoculated with a 1:100 dilution in fresh BHI at 37 °C and cells were harvested during the exponential growth at DO600 = 0.8 (Additional file 1: Figure S1).
Protein extraction and preparation of whole bacterial lysates for LC-MS/MS
After bacterial growth, the protein extraction was performed according to Silva et al. . The cultures were centrifuged at 4000 x g at 4 °C for 20 min. The cell pellets were washed in phosphate buffered saline (PBS) and then resuspended in 1 mL of lysis buffer (7 M Urea, 2 M Thiourea, CHAPS 4% and 1 M DTT) and 10 μL of Protease Inhibitor Mix (GE Healthcare, Piscataway, NJ, USA) was added. The cells were broken by sonication at 5 × 1 min cycles on ice and the lysates were centrifuged at 14,000 x g for 30 min at 4 °C. Subsequently, samples were concentrated and lysis buffer was replaced by 50 mM ammonium bicarbonate at pH 8.0 using a 10 kDa ultra-filtration device (Millipore, Ireland). All centrifugation steps were performed at room temperature. Finally the protein concentration was determined by Bradford method . A total of 50 μg proteins from each biological replicate of 1002_ovis and 258_equi were denatured by using RapiGEST SF [(0.1%) (Waters, Milford, CA, USA)] at 60 °C for 15 min, reduced with DTT [(10 mM) (GE Healthcare)], and alkylated with iodoacetamide [(10 mM) (GE Healthcare)]. For enzymatic digestion, trypsin [(0.5 μg/μL) (Promega, Sequencing Grade Modified Trypsin, Madison, WI, USA)] was added and placed in a thermomixer at 37 °C overnight. The digestion process was stopped by the addition of 10 μL of 5% TFA (Sigma-Aldrich, St. Louis, Missouri, USA) and glycogen phosphorilase (Sigma-Aldrich) was added to the digests to give 20 fmol.uL−1 as an internal standard for scouting normalization prior to each replicate injection into label-free quantitation .
LC-HDMSE analysis and data processing
Qualitative and quantitative analysis were performed using 2D RPxRP (two-dimensional reversed phase) nanoUPLC-MS (Nano Ultra Performance Liquid Chromatography Mass Spectrometry) approach with multiplexed Nano Electrospray High Definition Mass Spectrometry (nanoESI-HDMSE). To ensure that all samples were injected with the same amount into the columns and to ensure standardized molar values across all conditions, stoichiometric measurements based on scouting runs of the integrated total ion account (TIC) were performed prior to analysis. The experiments were conducted using both a 1 h reversed phase gradient from 7% to 40% (v/v) acetonitrile (0.1% v/v formic acid) and a 500 nL.min−1 on a 2D nanoACQUITY UPLC technology system . A nanoACQUITY UPLC HSS (High Strength Silica) T3 1.8 μm, 75 μm × 15 cm column (pH 3) was used in conjunction with a reverse phase (RP) XBridge BEH130 C18 5 μm 300 μm × 50 mm nanoflow column (pH 10). Typical on-column sample loads were 250 ng of total protein digests for each 5 fractions (250 ng/fraction/load). For all measurements, the mass spectrometer was operated in the resolution mode with a typical m/z resolving power of at least 35,000 FMHW and an ion mobility cell filled with nitrogen gas and a cross-section resolving power at least 40 Ω/ΔΩ. All analyses were performed using nano-electrospray ionization in the positive ion mode nanoESI (+) and a NanoLockSpray (Waters, Manchester, UK) ionization source.
The lock mass channel was sampled every 30 s. The mass spectrometer was calibrated with a MS/MS spectrum of [Glu1]-Fibrinopeptide B human (Glu-Fib) solution (100 fmol.uL−1) delivered through the reference sprayer of the NanoLockSpray source.The doubly- charged ion ([M + 2H]2+ = 785.8426) was used for initial single-point calibration and MS/MS fragment ions of Glu-Fib were used to obtain the final instrument calibration. Multiplexed data-independent (DIA) scanning with added specificity and selectivity of a non-linear ‘T-wave’ ion mobility (HDMSE) experiments were performed with a Synapt G2-S HDMS mass spectrometer (Waters), which was automatically planned to switch between standard MS (3 eV) and elevated collision energies HDMSE (19–45 eV) applied to the transfer ‘T-wave’ CID (collision-induced dissociation) cell with argon gas. The trap collision cell was adjusted for 1 eV, using a mili-seconds scan time previously adjusted based on the linear velocity of the chromatography peak delivered through nanoACQUITY UPLC to get a minimum of 20 scan points for each single peak, both in low energy and at high-energy transmission at an orthogonal acceleration time-of-flight (oa-TOF) from m/z 50 to 2000. The RF offset (MS profile) was adjusted is such a way that the nanoUPLC-HDMSE data are effectively acquired from m/z 400 to 2000, which ensured that any masses observed in the high energy spectra with less than m/z 400 arise from dissociations in the collision cell.
Database searching and quantification
Following the identification of proteins, the quantitative data were packaged using dedicated algorithms [23, 24] and searching against a database with default parameters to account for ions . The databases used were reversed “on-the fly” during the database queries and appended to the original database to assess the false positive rate (FDR) during identification. For proper spectra processing and database searching conditions, the Protein Lynx Global Server v.2.5.2 (PLGS) with IdentityE and ExpressionE informatics v.2.5.2 (Waters) were used. UniProtKB (release 2013_01) with manually reviewed annotations was used, and the search conditions were based on taxonomy (Corynebacterium pseudotuberculosis). We have utilized a database from genome annotation of 1002_ovis CP001809.2 version and 258_equi CP003540.2 version. These databases were randomized within PLGS v.2.5.2 for generate a concatenated database from both genomes. Thus, the measured MS/MS spectra from proteomic datasets of 1002_ovis and 258_equi were searched against this concatenated database. The maximum allowed missed cleavages by trypsin were up to one, and variable modifications by carbamidomethyl (C), acetyl N-terminal, phosphoryl (STY) and oxidation (M) were allowed and peptide mass tolerance value of 10 ppm was used . Peptides as source fragments, peptides with a charge state of at least [M + 2H]2+ and the absence of decoys were the factors we considered to increase the data quality. The collected proteins were organized by the PLGS ExpressionE tool algorithm into a statistically significant list that corresponded to higher or lower regulation ratios among the different groups. For protein quantitation, the PLGS v2.5.2 software was used with the IdentityE algorithm using the Hi3 methodology. The search threshold to accept each spectrum was the default value in the program with a false discovery rate value of 4%. The quantitative values were averaged over all samples, and the standard deviations at p < 0.05 were determined using the Expression software. Only proteins with a differential expression log2 ratio between the two conditions greater than or equal to 1.2 were considered .
The identified proteins in 1002_ovis and 258_equi were subjected to the bioinformatics analysis using the various prediction tools. SurfG+ v1.0  was used to predict sub-cellular localization, SignalP 4.1.0 server  to predict the presence of N-terminal signal peptides for secretory proteins, SecretomeP 2.0 server  to identify exported proteins from non-classical systems (positive prediction score greater than to 0.5), LipoP server  to determine lipoproteins, Blast2GO  and COG database  were used for functional annotations. The protein-protein interaction network was generated using Cytoscape version 2.8.3  with a spring-embedded layout.
Results and discussion
Characterization of the proteome of C. pseudotuberculosis biovar ovis and equi
In this study, we applied the 2D nanoUPLC-HDMSE approach to characterize the proteome of the strains 1002_ovis and 258_equi. Both strains were grown in BHI media, subsequently proteins were extracted and digested in solution, and then the peptides were analyzed by LC/MSE. Our proteomic analysis identified a total of 1227 non-redundant proteins in 1002_ovis (Additional file 2: Table S1 and Additional file 3: Table S2) and 1218 in 258_equi (Additional file 2: Table S1 and Additional file 4: Table S3) (Fig. 1a). The information about sequence coverage and a number of identified peptides for each protein sequence identified, as well as the information about the native peptide are available at Additional file 5: Table S4 and Additional file 6: Table S5. Altogether from the proteome of these two biovars, we identified a total of 1323 different proteins of C. pseudotuberculosis with high confidence (Fig. 1a) and characterized approximately 58% of the predicted proteome of 1002_ovis  (Fig. 1b). In the case of 258_equi, we characterized approximately 57% of the predicted proteome  (Fig. 1b). The proteins identified in both proteomes were analyzed by SurfG+ tool  to predict the subcellular localization into four categories: cytoplasmic (CYT), membrane (MEM), potentially surface-exposed (PSE) and secreted (SEC) (Fig. 1c). Further, we identified 83% (43 proteins) of the lipoproteins predicted in 1002_ovis and 79% (41 proteins) in 258_equi. Considering proteins with LPxTG motif which are involved in covalent linkage with peptidoglycan, we identified 6 proteins in 1002_ovis and 4 proteins in 258_equi that correspond to approximately 38% and 34% of the LPxTG proteins predicted in each strain, respectively.
The biovar equi and biovar ovis core proteome
The core-proteome, between 258_equi and 1002_ovis is composed of 1122 proteins (Fig. 1) (Additional file 2: Table S1). Interestingly, when correlated these 1122 proteins with in silico data of the C. pseudotuberculosis core-genome , we observed that 86% (960 proteins) of the Open Reading Frame (ORF) that encodes these proteins are part of the core-genome (Additional file 2: Table S1), what represents approximately 64% of the predicted core-genome of this pathogen. In addition, these data show a set of proteins involved in different cellular processes which could be necessary for the free living of C. pseudotuberculosis. The other 14% (262 proteins) of the proteins that constitute the core-proteome are shared by at least one of the 15 strains used in the core-genome study. According to Gene Ontology analysis [31, 32], the 1122 proteins were classified into four important functional groups: (i) metabolism, (ii) information storage and processing, (iii) cellular processes and signaling, and (iv) poorly characterized (Fig. 2a). As observed in the study of C. pseudotuberculosis  core genome in the categories “metabolism” and “information storage and processing” were detected a large number of proteins.
The label-free quantification was applied to evaluate the relative abundance of the core-proteome of 258_equi and 1002_ovis. The ProteinLynx Global Server (PLGS) v2.5.2 software with ExpressionE algorithm tool was used to identify proteins with p ≤ 0.05 (Additional file 2: Table S1). Among these proteins, 120 proteins between 258_equi and 1002_ovis showed difference in level of abundance (log2 ratios equal or greater than a factor of 1.2)  (Table 1). In this group of proteins that have presented different abundance level (258_equi:1002_ovis), 49 proteins were more abundant and 71 less abundant (Table 1). To visualize this differential distribution of the core-proteome a volcano plot of the log2 ratio of 258_equi/1002_ovis versus Log (e) Variance was generated (Fig. 2b). Interestingly, the Phospholipase D (Pld), the major virulence factor of C. pseudotuberculosis, was more abundant in 258_equi, than in 1002_ovis (Table 1). The Pld have an important play role in the pathogenic process of C. pseudotuberculosis, due to the sphingomyelinase activity of the Pld, this exotoxin increases vascular permeability through the exchange of polar groups attached to membrane-bound lipids and helps the bacteria in spread inside the host [34, 35]. In addition, this exotoxin is able to reduce the viability of both macrophages and neutrophils [34, 36]. In comparative proteomic studies between 1002_ovis and C231_ovis exoproteome, Pld was detected only in the C231_ovis supernatant [13, 15, 16]. A study performed with pld mutant strains presented decreased virulence . Thus, in relation to 258_equi, 1002_ovis could present a low potential of virulence.
The 120 differential proteins were organized by cluster of orthologous groups, and when evaluated the different biological processes that comprise each category listed above, we observed that 19 process were differentials between 258_equi and 1002_ovis (Fig. 2c, Additional file 7: Figure S2 and Additional file 8: Figure S3). The majority of the more abundant proteins (258_equi:1002_ovis) are related to cellular metabolism. On other hand, the majority of the less abundant proteins (258_equi:1002_ovis) are classified as poorly characterized or of unknown function. However, when proteins of known or predicted function are evaluated the majority of the less abundant proteins are related to cellular processes and signaling.
Difference among the major functional classes identified from the core-proteome analysis of 1002_ovis and 258_equi
During the infection process, pathogens need to adjust their metabolism in response to nutrient availability inside and outside the host. In our proteomic study, we identified several proteins related to different metabolic pathways. To determine the metabolic network of each strain, the proteins identified in this study were analyzed using Kyoto Encyclopedia of Genes pathways and Genomes (KEGG) . A total of 321 and 320 proteins, corresponding to 1002_ovis and 258_equi respectively, were mapped onto different metabolic pathways (Additional file 9: Figure S4 and Additional file 10: Figure S5). We observed differences in the metabolism of the biovars, related to Amino acid transport and metabolism, Carbohydrate transport and metabolism, Coenzyme metabolism, Energy metabolism, Lipid transport and metabolism, Nucleotide metabolism and Secondary metabolites biosynthesis, transport and catabolism. Difference in the metabolism cellular, also already were observed in others comparative proteomic study of C. pseudotuberculosis [13, 16, 17, 19], as well as in the Mycobacterium tuberculosis pathogen .
Interestingly, the PTS system fructose-specific EIIABC component (PstF) related to carbohydrate metabolism was more abundant in 258_equi, than in 1002_ovis (Table 1). This protein showed increased abundance in field isolates of C. pseudotuberculosis biovar ovis grown in BHI when compared to C231_ovis, a reference strain . This increased abundance of PstF in 258_equi, suggests that this protein could be important to the transport of carbon source both biovar ovis and biovar equi strains. On the other hand, the Precorrin 8X methyl mutase involved in cobalamin and vitamin B12 synthesis can be required only in biovar ovis strains, this protein beside being more abundant in 1002_ovis (Table 1), was also detected with greater abundance in the field isolates of C. pseudotuberculosis biovar ovis after having been grown in BHI . Glutamate dehydrogenase (GDH) was detected more abundant in 258_equi (Table 1). A study performed with the M. bovis pathogen showed that GDH contributes to the survival of this pathogen during macrophage infection .
In C. pseudotuberculosis, it was demonstrated that genes related the iron-acquisition are involved in the virulence of this pathogen . In the core-proteome of 1002_ovis and 258_equi, we detected proteins involved in this process, like CiuA, FagC and FagD; however, all these proteins were not differentially regulated between the two strains (Additional file 2: Table S1). On the other hand, HmuT protein, related to hemin uptake, was more abundant in 258_equi (Table 1). Additionally, we have also detected a cell surface hemin receptor in the exclusive proteome of this strain. Heme represents the major reservoir of iron source for many bacterial pathogens that rely on surface-associated heme-uptake receptors . The HmuT is a lipoprotein that acts as a hemin receptor. The hmuT gene is part of the operon hmuTUV, an ABC transport system (haemin transport system), which is normally present in pathogenic Corynebacterium [43, 44]. In addition, in the pathogen C. ulcerans, HmuT is required for normal hemin utilization .
Information storage and processing
Of the total protein of proteins identify in the category “information storage and processing” the majority of the differential proteins were less abundant in 258_equi (Table 1). Only, Metallophosphoesterase involved in DNA repair, SAM dependent methyltransferase related to transcriptional process and Ribosomal RNA small subunit methyltransferase I involved in translation process were more induced in 258_equi. In 1002_ovis the Exodeoxyribonuclease 7 important protein related to the DNA-damage pathway was more induced in this strain. In addition, we identified the TetR family regulatory protein as more abundant in 1002_ovis, this result was also observed in field isolates of C. pseudotuberculosis from sheep infected naturally . TerR proteins are related to regulation of multidrug efflux pumps, antibiotic biosynthesis, catabolic process and cellular differentiation process . Others important transcriptional regulators also were induced in 1002_ovis such as PvdS and GreA regulators.
Cellular processes and signaling
Our proteomic analyses detected differentially regulated proteins belonging to different antioxidant systems. These could contribute to the survival of C. pseudotuberculosis in various stress conditions, such as reactive oxygen species (ROS) and reactive nitrogen species (RNS), which are generally found in macrophage. The three major thiol-dependent antioxidant systems in prokaryotic pathogens are the thioredoxin system (Trx), the glutathione system (GSH-system) and the catalase system . Thioredoxin TrxA and Thiol-disulfide isomerase thioredoxin were more abundant in 258_equi (Table 1). These proteins are involved in the Trx-system, which has a major role against oxidative stress . However, proteins like catalase and glutaredoxin (nrdH) were less abundant in 258_equi (Table 1), being more active in 1002_ovis. Catalase plays an important role in resistance to ROS and RNS, as well as in the virulence of M. tuberculosis . The protein NrdH has a glutaredoxin amino acid sequence and thioredoxin activity. It is present in Escherichia coli  and C. ammoniagenes , as well as in bacteria where the GSH system is absent, such as M. tuberculosis . Thus, the presence of NrdH may represent one more factor that contributes to the resistance of C. pseudotuberculosis against ROS and RNS during the infection process, as well as to the maintenance of the balance of intracellular redox potential. Proteins like NorB and Glyoxalase/Bleomycin, which play roles in the nitrosative stress response of 1002_ovis, were identified in the exclusive proteome of this strain (Additional file 3: Table S2) [14, 18]. These results shown that beside of present proteins with difference in abundance both strains present a set of proteins that could contribute to adaptive process under stress conditions.
Difference proteomic observed in the exclusive proteome of 258_equi and 1002_ovis
We found respectively 105 and 96 proteins in the exclusive proteome of 1002_ovis and 258_equi (Fig. 1) (Additional file 3: Table S2 and Additional file 4: Table S3), related to different biological process (Additional file 7: Figure S2 and Additional file 8: Figure S3). Interestingly, in this exclusive proteome of 1002_ovis and 258_equi, we detected specific proteins in each strain (Table 2, Additional file 3: Table S2 and Additional file 4: Table S3). In the exclusive proteome of 258_equi, the ORFs that codify twenty proteins are annotated as pseudogene in 1002_ovis (Table 2, Additional file 3: Table S2 and Additional file 4: Table S3). On the other hand, the ORFs that encode six proteins were not detected in the genome of 1002_ovis. These proteins are two CRISPR, MoeB, and three unknown function proteins. CRISPR is an important bacterial defense system against infections by viruses or plasmids, this immunity is obtained from the integration of short sequences of invasive DNA ‘spacers’ into the CRISPR loci .
The distinction between the biovar ovis and biovar biovar equi strains is based on a biochemical assay, where biovar ovis strains are negative for nitrate reduction, whereas biovar equi strains are positive . However, to date, there is no available information regarding the molecular basis underlying nitrate reduction in C. pseudotuberculosis biovar equi. MoeB is involved in the molybdenum cofactor (Moco) biosynthesis, which plays an important role in anaerobic respiration in bacteria and also are required to activation of nitrate reductase (NAR) . In the closely related pathogen M. tuberculosis several studies have showed the great importance of molybdenum cofactor in its virulence and pathogenic process, mainly macrophage intracellular environmental . Therefore, more studies are necessary to explore the true role of Moco both physiology and virulence of biovar equi strains. Other protein that also could contribute to resistance of 258_equi macrophage is NADPH dependent nitro/flavin reductase (NfrA), a pseudogene in 1002_ovis. In addition, studies performed in Bacillus subtilis showed that NfrA is involved in both oxidative stress  and heat shock resistance .
In 1002_ovis, only the ORF that encodes a DNA methylase was not found in the 258_equi genome (Table 2, Additional file 3: Table S2 and Additional file 4: Table S3). In addition, the ORFs that codifies seven proteins identified in the exclusive proteome of the strain 1002_ovis are annotated like pseudogene in 258_equi (Table 2, Additional file 3: Table S2 and Additional file 4: Table S3). Inside this group, we have identified important proteins involved in the process of adhesion and invasion cellular, which might contribute in the pathogenesis of 1002_ovis. Adhesion to host cells is a crucial step that favors the bacterial colonization; this process is mediated by different adhesins . We identified proteins such as: collagen binding surface protein Cna-like and Sdr family related adhesin, which are members of the collagen-binding microbial surface components recognizing adhesive matrix molecules (MSCRAMMs) (Table 2). This class of proteins is present in several Gram positive pathogens and plays an important role in bacterial virulence by acting mainly in the cellular adhesion process [58,59,60,61].
Another detected protein that might contribute to the virulence of 1002_ovis is Neuraminidase (NanH) (Table 2). This protein belongs to a class of glycosyl hydrolases that contributes to the recognition of sialic acids exposed on host cell surfaces . In C. diphtheriae, it was demonstrated that a protein with trans-sialidase activity promotes cellular invasion [63, 64]. In addition, NanH was reported to be immunoreactive in the immunoproteome of 1002_ovis, showing the antigenicity of this protein . Interestingly, genomic difference in relation to gene involved in the adhesion and invasion process, also already were observed between biovar ovis strain and biovar equi strains, mainly in genes related to pilus [10, 12]. According to pathogenic process of each biovar, unlike biovar equi strains, which rarely causes visceral lesions , biovar ovis strains, are responsible mainly by visceral lesions [2, 35], what requires a high ability to adhere and invade the host cell, thus these protein could be responsible by this ability of biovar ovis strain in attacks visceral organs.
In our proteomic analysis, the measured MS/MS spectra from the proteomic datasets of 1002_ovis and 258_equi were searched against a concatenated database composed by genome annotation of 1002_ovis CP001809.2 version and 258_equi CP003540.2 version for identify possible errors or unannotated genes. Thus, by adopting more stringent criteria of considering only proteins with a minimum representative of two peptides and a FDR < 1%, we identified five proteins in 1002_ovis and seven proteins in 258_equi, which were not previously annotated. All parameters, as well as, the peptides sequence which were used for identification of these proteins are shown in Additional file 11: Table S6 and Additional file 12: Table S7. The proteins identified in this proteogenomic analysis are associated to different biological processes. For instance, the Aminopeptidase N involved in the amino acid metabolism was detected in 1002_ovis, whereas the Cobaltochelatase (cobN), associated to cobalt metabolism, glutamate dehydrogenase (gdh) involved in the L-glutamate metabolism, the PTS system fructose specific EIIABC related to fructose metabolism and the Phosphoribosylglycinamide formyltransferase involved in the purine biosynthesis were all detected in 258_equi. Proteins involved in DNA processes, such as Uracil DNA glycosylase in 258_equi; and Exodeoxyribonuclease 7 small subunit in 1002_ovis were also detected in both strains. Proteins with general function prediction only and unknown function were also identified in both strains.
In conclusion, we used a label-free quantitative approach to compare, for the first time, the proteome of C. pseudotuberculosis strains belonging to both ovis and equi biovars. Taken together, the findings reported here show a set of shared and exclusive factors of 1002_ovis and 258_equi at the protein level, which can contribute to understanding both the physiology and the virulence of these strains. In addition, the functional analysis of the genome of 1002_ovis and 258_equi allows the in silico validation of data of the genome of these strains. Thus, the proteins identified here may be used as potential new targets for the development of vaccines against ovis and equi C. pseudotuberculosis in future investigations.
Availability of supporting data
The datasets supporting the results of this article were then concatenated into a *xlsx file at peptide and protein level to fulfill the requirements and is available at supplemental material including sequence coverage and a number of identified peptides for each protein sequence identified. It also includes the native peptide information.
Two-dimensional reversed phase
False discovery rate
High strength silica
Liquid chromatograph mass spectrometry
Nano Electrospray High Definition Mass Spectrometry
Nano Ultra Performance Liquid Chromatography
Nano Ultra Performance Liquid Chromatography Mass Spectrometry
Phosphate buffered saline
Protein Lynx Global Server
Reactive nitrogen species
Reactive oxygen species
Dorella FA, Pacheco LG, Oliveira SC, Miyoshi A, Azevedo V. Corynebacterium pseudotuberculosis: microbiology, biochemical properties, pathogenesis and molecular studies of virulence. Vet Res. 2006;37:201–18.
Baird GJ, Fontaine MC. Corynebacterium pseudotuberculosis and its role in ovine caseous lymphadenitis. J Comp Pathol. 2007;137:179–210.
Shpigel NY, Elad D, Yeruham I, Winkler M, Saran A. An outbreak of Corynebacterium pseudotuberculosis infection in an Israeli dairy herd. Vet Rec. 1993;133:89–94.
Britz E, Spier SJ, Kass PH, Edman JM, Foley JE. The relationship between Corynebacterium pseudotuberculosis biovar equi phenotype with location and extent of lesions in horses. Vet J. 2014;200:282–6.
Selim SA. Oedematous skin disease of buffalo in Egypt. J Vet Med B Infect Dis Vet Public Health. 2001;48:241–58.
Paton MW, Walker SB, Rose IR, Watt GF. Prevalence of caseous lymphadenitis and usage of caseous lymphadenitis vaccines in sheep flocks. Aust Vet J. 2002;81:91–5.
Foley JE, Spier SJ, Mihalyi J, Drazenovich N, Leutenegger CM. Molecular epidemiologic features of Corynebacterium pseudotuberculosis isolated from horses. Am J Vet Res. 2004;65:1734–7.
Seyffert N, Guimarães AS, Pacheco LG, Portela RW, Bastos BL, Dorella FA, et al. High seroprevalence of caseous lymphadenitis in Brazilian goat herds revealed by Corynebacterium pseudotuberculosis secreted proteins-based ELISA. Res Vet Sci. 2010;88:50–5.
Kumar J, Singh F, Tripathi BN. Kumar 454 R, Dixit SK, Sonawane GG. Epidemiological, bacteriological and molecular studies on caseous lymphadenitis in Sirohi goats of Rajasthan, India. Trop Anim Health Prod. 2012;44:1319–22.
Soares SC, Silva A, Trost E, Blom J, Ramos R, Carneiro A, et al. The pan-genome of the animal pathogen Corynebacterium pseudotuberculosis reveals differences in genome plasticity between the biovar ovis and biovar equi strains. PLoS One. 2013;8:e53818.
Ruiz JC, D’Afonseca V, Silva A, Ali A, Pinto AC, Santos AR, et al. Evidence for reductive genome evolution and lateral acquisition of virulence functions in two Corynebacterium pseudotuberculosis strains. PLoS One. 2011;6:e18551.
Soares SC, Trost E, Ramos RTJ, Carneiro AR, Santos AR, Pinto AC, et al. Genome sequence of Corynebacterium pseudotuberculosis biovar equi strain 258 and prediction of antigenic targets to improve biotechnological vaccine production. J Biotechnol. 2012;20:135–41.
Pacheco LG, Slade SE, Seyffert N, Santos AR, Castro TL, Silva WM, et al. A combined approach for comparative exoproteome analysis of Corynebacterium pseudotuberculosis. BMC Microbiol. 2011;17:12.
Pacheco LG, Castro TL, Carvalho RD, Moraes PM, Dorella FA, Carvalho NB, et al. A Role for Sigma Factor σ(E) in Corynebacterium pseudotuberculosis Resistance to Nitric Oxide/Peroxide Stress. Front Microbiol. 2012;3:126.
Silva WM, Seyffert N, Santos AV, Castro TL, Pacheco LG, Santos AR, et al. Identification of 11 new exoproteins in Corynebacterium pseudotuberculosis by comparative analysis of the exoproteome. Microb Pathog. 2013a;16:37–42.
Silva WM, Seyffert N, Ciprandi A, Santos AV, Castro TL, Pacheco LG, et al. Differential Exoproteome analysis of two Corynebacterium pseudotuberculosis biovar ovis strains isolated from goat (1002) and sheep (C231). Curr Microbiol. 2013;67:460–5.
Rees MA, Kleifeld O, Crellin PK, Ho B, Stinear TP, Smith AI, et al. Proteomic characterization of a natural host-pathogen interaction: repertoire of in vivo expressed bacterial and host surface-associated proteins. J Proteome Res. 2015;14:120–32.
Silva WM, Carvalho RD, Soares SC, Bastos IF, Folador EL, Souza GH, et al. Label free proteomic analysis to confirm the predicted proteome of Corynebacterium pseudotuberculosis under nitrosative stress mediated by nitric oxide. BMC Genomics. 2014;15:1065.
Rees MA, Stinear TP, Goode RJA, Coppel RL, Smith AI, Kleifeld O. Changes in protein abundance are observed in bacterial isolates from a natural host. Front Cell Infect Microbiol. 2015;5:71.
Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72:248–54.
Silva JC, Gorenstein MV, Li GZ, Vissers JP, Geromanos SJ. Absolute quantification of proteins by LC-MSE: a virtue of parallel MS acquisition. Mol Cell Proteomics. 2006;5:144–56.
Gilar M, Olivova P, Daly AE, Gebler JC. Two-dimensional separation of peptide using RP-RP-HPLC system with different pH in first and second separation dimensions. J Sep Sci. 2005;28:1694–703.
Silva JC, Denny R, Dorschel CA, Gorenstein M, Kass IJ, Li GZ, et al. Quantitative proteomic analysis by accurate mass retention time pairs. Anal Chem. 2005;1:2187–200.
Geromanos SJ, Vissers JP, Silva JC, Dorschel CA, Li GZ, Gorenstein MV, et al. The detection, correlation, and comparison of peptide precursor and product ions from data independent LC-MS with data dependant LC-MS/MS. Proteomics. 2009;9:1683–95.
Li GZ, Vissers JP, Silva JC, Golick D, Gorenstein MV, Geromanos SJ. Database searching and accounting of multiplexed precursor and product ion spectra from the data independent analysis of simple and complex peptide mixtures. Proteomics. 2009;9:1696–719.
Levin Y, Hradetzky E, Bahn S. Quantification of proteins using data-independent analysis (MSE) in simple andcomplex samples: a systematic evaluation. Proteomics. 2011;11:3273–87.
Barinov A, Loux V, Hammani A, Nicolas P, Langella P, Ehrlich D, et al. Prediction of surface exposed proteins in Streptococcus pyogenes, with a potential application to other Gram-positive bacteria. Proteomics. 2009;9:61–73.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Bendtsen JD, Kiemer L, Fausboll A, Brunak S. Non-classical protein secretion in bacteria. BMC Microbiol. 2005;5:58.
Juncker AS, Willenbrock H, Von Heijne G, Brunak S, Nielsen H, Krogh A. Prediction of lipoprotein signal peptides in Gram-negative bacteria. Protein Sci. 2003;12:1652–62.
Conesa A, Gotz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;15:3674–6.
Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, et al. The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res. 2001;29:22–8.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
McKean SC, Davies JK, Moore RJ. Expression of phospholipase D, the major virulence factor of Corynebacterium pseudotuberculosis, is regulated by multiple environmental factors and plays a role in macrophage death. Microbiology. 2007;7:2203–11.
Batey RG. Pathogenesis of caseous lymphadenitis in sheep and goats. Aust Vet J. 1986;63:269–72.
Yozwiak ML, Songer JG. Effect of Corynebacterium pseudotuberculosis phospholipase D on viability and chemotactic responses of ovine neutrophils. Am J Vet Res. 1993;54:392–7.
McNamara PJ, Bradley GA, Songer JG. Targeted mutagenesis of the phospholipase D gene results in decreased virulence of Corynebacterium pseudotuberculosis. Mol Microbiol. 1994;12:921–30.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Gunawardena HP, Feltcher ME, Wrobel JA, Gu S, Braunstein M, Chen X. Comparison of the membrane proteome of virulent Mycobacterium tuberculosis and the attenuated Mycobacterium bovis BCG vaccine strain by label-free quantitative proteomics. J Proteome Res. 2013;6:5463–74.
Gallant JL, Viljoen AJ, van Helden PD, Wiid IJ. Glutamate Dehydrogenase Is Required by Mycobacterium bovis BCG for Resistance to Cellular Stress. PLoS One. 2016;29:e0147706.
Billington SJ, Esmay PA, Songer JG, Jost BH. Identification and role in virulence of putative iron acquisition genes from Corynebacterium pseudotuberculosis. FEMS Microbiol Lett. 2002;19:41–4.
Contreras H, Chim N, Credali A, Goulding CW. Heme uptake in bacterial pathogens. Curr Opin Chem Biol. 2014;19:34–41.
Drazek ES, Hammack AC, Schmitt PM. Corynebacterium diphtheriae genes required for acquisition of iron from hemin and hemoglobin are homologous to ABC hemin transporters. Mol Microbiol. 2000;36:68–84.
Schmitt MP, Drazek SE. Construction and consequences of directed mutations affecting the hemin receptor in pathogenic Corynebacterium species. J Bacteriol. 2001;183:1476–81.
Ramos JL, Martínez-Bueno M, Molina-Henares AJ, Terán W, Watanabe K, Zhang X, et al. The TetR family of transcriptional repressors. Microbiol Mol Biol Rev. 2005;69:326–56.
Lu J, Holmgren A. The thioredoxin antioxidant system. Free Radic Biol Med. 2014;8:75–87.
Ng VH, Cox JS, Sousa AO, MacMicking JD, McKinney JD. Role of KatG catalase peroxidase in mycobacterial pathogenesis: countering the phagocyte oxidative burst. Mol Microbiol. 2004;52:1291–302.
Jordan A, Aslund F, Pontis E, Reichard P, Holmgren A. Characterization of Escherichia coli NrdH. A glutaredoxin-like protein with a thioredoxin-like activity profile. J Biol Chem. 1997;18:18044–50.
Stehr M, Lindqvist Y. NrdH-redoxin of Corynebacterium ammoniagenes forms a domain-swapped dimer. Proteins Struct Funct Bioinform. 2004;55:613–9.
Leiting WU, Jianping XI. Comparative genomics analysis of Mycobacterium NrdH redoxins. Microb Pathog. 2010;48:97–102.
Jiang W, Marraffini LA. CRISPR-Cas: New Tools for Genetic Manipulations from Bacterial Immunity Systems. Annu Rev Microbiol. 2015;15:209–28.
Biberstein EL, Knight HD, Jang S. Two biotypes of Corynebacterium pseudotuberculosis. Vet Rec. 1971;25:691–2.
Leimkühler S, Iobbi-Nivol C. Bacterial molybdoenzymes: old enzymes for new purposes. FEMS Microbiol Rev. 2015;13:1–18.
Williams M, Mizrahi V, Kana BD. Molybdenum cofactor: a key component of Mycobacterium tuberculosis pathogenesis? Crit Rev Microbiol. 2014;40:18–29.
Mostertz J, Scharf C, Hecker M, Homuth G. Transcriptome and proteome analysis of Bacillus subtilis gene expression in response to superoxide and peroxide stress. Microbiology. 2004;150:497–512.
Moch C, Schrögel O, Allmansberger R. Transcription of the nfrA-ywcH operon from Bacillus subtilis is specifically induced in response to heat. J Bacteriol. 2000;182:4384–93.
Rogers EA, Das A, Ton-That H. Adhesion by pathogenic corynebacteria. Adv Exp Med Biol. 2011;715:91–103.
Patti JM, Allen BL, McGavin MJ, Hook M. MSCRAMM-mediated adherence of microorganisms to host tissues. Annu Rev Microbiol. 1994;48:585–617.
Lannergard J, Frykberg L, Guss B. CNE, a collagen-binding protein of Streptococcus equi. FEMS Microbiol Lett. 2003;222:69–74.
Nallapareddy SR, Weinstock GM, Murray BE. Clinical isolates of Enterococcus faecium exhibit strain-specific collagen binding mediated by Acm, a new member of the MSCRAMM family. Mol Microbiol. 2003;47:1733–47.
Kang M, Ko YP, Liang X, Ross CL, Liu Q, Murray BE, et al. Collagen-binding Microbial Surface Components Recognizing Adhesive Matrix Molecule (MSCRAMM) of Gram-positive Bacteria Inhibit Complement Activation via the Classical Pathway. J Biol Chem. 2013;12:20520–31.
Vimr ER, Kalivoda KA, Deszo EL, Steenbergen SM. Diversity of microbial sialic acid metabolism. Microbiol Mol Biol Rev. 2004;68:132–53.
Mattos-Guaraldi AL, Duarte Formiga LC, Pereira GA. Cell surface components and adhesion in Corynebacterium diphtheriae. Microbes Infect. 2000;2:1507–12.
Kim S, Oh DB, Kwon O, Kang HA. Identification and functional characterization of the NanH extracellular sialidase from Corynebacterium diphtheriae. J Biochem. 2010;147:523–33.
Seyffert N, Silva RF, Jardin J, Silva WM, Castro TL, Tartaglia NR, et al. Serological proteome analysis of Corynebacterium pseudotuberculosis isolated from different hosts reveals novel candidates for prophylactics to control caseous lymphadenitis. Vet Microbiol. 2014;7:255–60.
This work involved the collaboration of various institutions, including the Genomics and Proteomics Network of the State of Pará of the Federal University of Pará, the Amazon Research Foundation (FAPESPA), the National Council for Scientific and Technological Development (CNPq), the Brazilian Federal Agency for the Support and Evaluation of Graduate Education (CAPES), the Minas Gerais Research Foundation (FAPEMIG) and the Waters Corporation, Brazil.
WMS performed microbiological analyses and sample preparation for proteomic analysis. GHMFS and WMS conducted the proteomic analysis. SCS and ELF performed bioinformatics analysis of the data. CSS, AVS, AM and HF contributed substantially to data interpretation and revisions. VA, AS and YLL participated in all steps of the project as coordinators, and critically reviewed the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
No ethics approval was required for any aspect of this study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Growth rates in BHI media of 1002_ovis (blue circles) and 258_equi (red triangles). (JPEG 278 kb)
Total list of proteins identified in the core-proteome of 1002_ovis and 258_equi. (XLSX 215 kb)
Total list of proteins identified in the exclusive proteome of 1002_ovis. (XLSX 20 kb)
Total list of proteins identified in the exclusive proteome of 258_equi. (XLSX 21 kb)
Total list of peptide and proteins identified 1002_ovis. (XLSB 31769 kb)
Total list of peptide and proteins identified 258_equi. (XLSB 33204 kb)
The protein-protein interaction network of 1002_ovis. (A) General interactome of differentially regulated proteins, identified in the exclusive proteome of 1002_ovis. The proteins are marked with different shapes: exclusive proteome, circle; more abundant, square; less abundant, rhombus. The biological processes were marked with different colors: amino acid transport and metabolism, yellow; secondary metabolites biosynthesis, transport and catabolism, aquamarine; inorganic ion transport and metabolism, orange; coenzyme metabolism, brown; carbohydrate transport and metabolism, chartreuse green; nucleotide metabolism, cerulean; energy metabolism, olive; lipid transport and metabolism, viridian; adhesion and motility cell, crimson; iuntracellular trafficking secretion and vesicular transport, persian blue; signal transduction mechanisms, maroon; cell wall/membrane and envelope, gray; defense mechanism, red; post-translational modification, protein turnover, chaperones, electric blue; DNA metabolism, replication, recombination and repair, violet; translation, ribosomal structure and biogenesis, amber; transcription, regulation, degradation and RNA processing, salmon; poorly characterized, white. (JPEG 3310 kb)
The protein-protein interaction network of 258_equi. (A) General interactome of the differentially regulated proteins, identified in the exclusive proteome of 258_equi. The proteins are marked with different shapes: exclusive proteome, circle; more abundant, square; less abundant, rhombus. The biological processes are marked with different colors: amino acid transport and metabolism, yellow; secondary metabolites biosynthesis, transport and catabolism, aquamarine; inorganic ion transport and metabolism, orange; coenzyme metabolism, brown; carbohydrate transport and metabolism, chartreuse green; nucleotide metabolism, cerulean; energy metabolism, olive; lipid transport and metabolism, viridian; adhesion and motility cell, crimson; intracellular trafficking secretion and vesicular transport, persian blue; signal transduction mechanisms, maroon; cell wall/membrane and envelope, gray; defense mechanism, red; post-translational modification, protein turnover, chaperones, electric blue; DNA metabolism, replication, recombination and repair, violet; translation, ribosomal structure and biogenesis, amber; transcription, regulation, degradation and RNA processing, salmon; poorly characterized, white. (JPEG 4178 kb)
Metabolic network of 1002_ovis. Red line, proteins identified in the proteomic analysis, other colors represent proteins not identified in this study. (JPEG 8633 kb)
Metabolic network of 258_equi. Red line, proteins identified in the proteomic analysis, other colors represent proteins not identified in this study. (JPEG 1267 kb)
Proteins identified in 1002_ovis by Proteogenomics. (XLSX 216 kb)
Proteins identified in 258_equi by Proteogenomics. (XLSX 266 kb)
About this article
Cite this article
Silva, W.M., Folador, E.L., Soares, S.C. et al. Label-free quantitative proteomics of Corynebacterium pseudotuberculosis isolates reveals differences between Biovars ovis and equi strains. BMC Genomics 18, 451 (2017). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-017-3835-y
- Corynebacterium pseudotuberculosis
- Caseous lymphadenitis
- Ulcerative lymphangitis
- Proteomic bacterial
- Label-free proteomics, proteomic