Disruption of the S41 Peptidase Gene in Mycoplasma mycoides capri Impacts Proteome Profile, H2O2 Production, and Sensitivity to Heat Shock

doi:10.1371/journal.pone.0051345

PLoS One. 2012; 7(12): e51345.

Published online 2012 December 31. doi: 10.1371/journal.pone.0051345

PMCID: PMC3534093

Disruption of the S41 Peptidase Gene in Mycoplasma mycoides capri Impacts Proteome Profile, H₂O₂ Production, and Sensitivity to Heat Shock

Ayman B. Allam, Mary B. Brown,^* and Leticia Reyes

Department of Infectious Diseases and Pathology, College of Veterinary Medicine, University of Florida, Gainesville, Florida, United States of America

Ilse D. Jacobsen, Editor

Leibniz Institute for Natural Products Research and Infection Biology- Hans Knoell Institute, Germany

* E-mail: mbbrown/at/ufl.edu

Competing Interests: The authors have declared that no competing interests exist.

Conceived and designed the experiments: ABA MB LR. Performed the experiments: ABA LR. Analyzed the data: ABA MB LR. Contributed reagents/materials/analysis tools: ABA MB LR. Wrote the paper: ABA MB LR.

Received July 9, 2012; Accepted November 7, 2012.

This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Members of the Mycoplasma mycoides cluster are among the most virulent of the mycoplasmas, causing worldwide economically significant diseases of cattle and goats. A distinguishing phenotype among the members of the cluster is the ability to degrade casein. The MMCAP2_0241 gene, an S41 peptidase, confers the proteolytic phenotype in Mycoplasma mycoides subsp. capri GM12. In order to determine the impact of disruption of the gene, we used differential proteome profiling to compare the M. mycoides subsp. capri wild type with a mutant lacking the proteolytic phenotype. Disruption of MMCAP2_0241 resulted in altered phenotypes reminiscent of M. mycoides subsp. mycoides SC and had significant impacts on the proteome profile of the microbe. The mutant exhibited increased production of hydrogen peroxide, decreased lactate dehydrogenase activity, and increased sensitivity to heat shock.

Introduction

Members of the Mycoplasma mycoides cluster are among the most virulent of the mycoplasmas, causing worldwide economically significant diseases of cattle and goats [1], [2]. Two members of the cluster, M. mycoides subsp. mycoides Small Colony type and M. capricolum subsp. capripneumoniae (formerly F38 group), cause contagious bovine and caprine pleuropneumonia, respectively and are listed (http://www.selectagents.gov/select agents and Toxins list.html) by USDA and APHIS (Federal Register 67, No. 155, 9 CFR 121.2b) and the World Organization for Animal Health (http://www.oie.int). Two closely related members of the M. mycoides cluster [3], [4], M. mycoides subsp. capri (formerly M. mycoides subsp. mycoides Large Colony) and M. capricolum, are associated with respiratory disease in small ruminants worldwide and may also cause extrapulmonary complications and sepsis, but these pathogens are not listed agents.

Prior to the advent of genomic analysis, members of the M. mycoides cluster were speciated based on key phenotypic differences [5], [6], [7], [8], including colony size on agar, biochemical activities, and thermal stability. A key phenotypic difference between the two species causing contagious pleuropneumonia and other members of the cluster, including M. mycoides subsp. capri GM12, is the ability to degrade casein [5], [6], [7], an activity that is quite uncommon in mycoplasmas. In addition to the differences in proteolytic activity, M. mycoides subsp. capri and M. mycoides subsp. mycoides SC also have different phenotypes with respect to carbohydrate utilization and thermal stability [8]. Further, production of H₂O₂ is a key virulence factor in the M. mycoides cluster [9], [10], [11], [12], [13], [14], [15], and quantitative differences in H₂O₂ production are directly linked to disease severity in mycoplasmas [9], [10], [13].

We identified the gene (MMCAP2_0241) responsible for the proteolytic phenotype in M. mycoides subsp. capri GM12 ATCC 35297 and created a mutant that has the coding sequence of the gene interrupted by insertion of tetM via homologous recombination [16]. In silico analysis of the predicted protein product of MMCAP2_0241 indicated a 651 amino acid protein with an unknown N-terminal domain (aa 25–339) and the S41 peptidase tail specific protease (TSP) domain (aa 340 to 544). Within the TSP domain, two carboxyl-tail processing (Ctp) motifs were found from aa 411–422 and aa 477–507. Also present in MMCAP2_0241 was a signal peptide domain (aa 1–24) and two transmembrane domains (aa 7–29 and 631–650). The alignment of MMCAP2_0241 showing conserved domains with other selected CtpA proteins is provided in Fig. S1, supplementary material. Because of the presence of the TSP domain, the presence of Ctp motifs, and homology with similar Ctp proteins in other bacterial species, we previously referred to this gene as ctpA [16]. Subsequently, we used the PHYRE (http://www.sbg.bio.ic.ac.uk/phyre2) automatic fold recognition server [17] to evaluate MMCAP2_0241. Based on this analysis, 266 residues (100% confidence, 41% coverage) modeled to the d1fc6a4 template for the superfamily ClpP, which also contains the TSP domain and is a carboxyl-tail processing S41 peptidase [18]. Because clpP is known to encode a caseinolytic protease [18], [19], we suggest that clpP-like rather than ctpA may be the more appropriate designation. However, for consistency with the previous published designation [16], we will refer to the mutant as the M. mycoides subsp. capri ctpA(clpP)::tetM mutant and to the protein as MMCAP2_0241 (ClpP-like).

Our M. mycoides subsp. capri ctpA(clpP)::tetM mutant, like the etiologic agents of contagious pleuropneumonia, lacks the proteolytic phenotype [8], [16]. Therefore, it was of interest to determine if disruption of MMCAP2_0241 had significant impacts on the proteome profile of the microbe as well as other phenotypes that differentiate among members of the M. mycoides cluster [8]. Here we report that disruption of MMCAP2_0241 resulted in altered phenotypes reminiscent of M. mycoides subsp. mycoides SC. Specifically, the mutant exhibited increased sensitivity to heat shock as well as increased production of H₂O₂. We also observed decreased lactate dehydrogenase (LDH) activity and significant changes in the proteome profile in the M. mycoides subsp. capri ctpA(clpP)::tetM mutant.

Results

Disruption of MMCAP2_0241 resulted in increased sensitivity to heat shock and increased H₂O₂ production

As previously reported [16], the M. mycoides subsp. capri ctpA(clpP)::tetM mutant lost the proteolytic phenotype, as evidenced by lack of enzymatic activity on casein agar; MMCAP2_0241 gene disruption also was demonstrated by both Northern and Southern blots. With respect to carbohydrate use, the M. mycoides subsp. capri ctpA(clpP)::tetM mutant did not differ from the wild type with respect to fermentation of sorbitol or oxidation of maltose, trehalose, mannose or glucosamine (data not shown).

The disruption of MMCAP2_0241 altered the ability of the mutant to respond to heat stress (Fig. 1), resulting in a phenotype that was more similar to that of M. mycoides subsp. mycoides SC. At 37°C, the log CFU and growth rates of the wild type and M. mycoides subsp. capri ctpA(clpP)::tetM mutant were comparable at all time points (Fig. 1). However, after 3 hr at 42°C, the mutant grew more slowly than did the wild type strain (P<0.01). At subsequent time points, the difference was even more dramatic (P<0.001). The CFU of the wild type continued to increase, but the mutant failed to grow, indicating an increased sensitivity to heat shock. The M. mycoides subsp. capri ctpA(clpP)::tetM mutant exhibited a significant (P<0.01) increase in the production of H₂O₂ (Fig. 2), with the mutant producing approximately 40% more H₂O₂ than the wild type.

Figure 1

Heat shock response in the Mycoplasma mycoides capri GM12 wild type and the ctpA(clpP)::tetM mutant.

Figure 2

H₂0₂ production in the Mycoplasma mycoides capri GM12 wild type and the ctpA(clpP)::tetM mutant.

Disruption of MMCAP2_0241 had a significant impact on the proteome profile

In other microbes, deletion of ctpA or clpP is associated with pleiotropic effects and proteome differences [20], [21], [22], [23]. Therefore, we next used two complementary approaches, 2D gel electrophoresis/differential gel electrophoresis (DIGE) and amine specific peptide-based labeling (iTRAQ™) followed by tandem mass spectrometry, to identify proteins that were perturbed in the M. mycoides subsp. capri ctpA(clpP)::tetM mutant. Both methods assess relative differences in the protein concentrations between the M. mycoides subsp. capri wild type and the M. mycoides subsp. capri ctpA(clpP)::tetM mutant. However, these two methodologies differ in that DIGE is protein-centric whereas iTRAQ is peptide-centric [24], [25], [26]. There are also discrepancies based on cellular compartmentalization (cytosol vs. membrane), abundance, and limitations inherent to each assay.

In our 2-D gels, eight proteins were found to differ significantly (P<0.006) between the M. mycoides subsp. capri wild type and M. mycoides subsp. capri ctpA(clpP)::tetM mutant (Table 1 and Fig. S2, supplementary material). Proteins that were significantly decreased in M. mycoides subsp. capri ctpA(clpP)::tetM were preprotein translocase SecA, adenylsuccinate synthase, phosphoglycerate kinase, and two conserved hypothetical proteins. The only proteins significantly increased in the mutant were the transcription anti-termination protein NusG and serine-tRNA ligase. TetM was present only in the M. mycoides subsp. capri ctpA(clpP)::tetM mutant, as expected.

Table 1

Proteins that significantly differed (P<0.006) in the M. mycoides subsp. capri ctpA::tetM mutant as determined by differential 2-dimensional gel electrophoresis.

The iTRAQ™ analysis was more revealing than 2D-DIGE analysis in that 221 proteins were identified with 95% confidence and an error factor <2 [27] (for complete dataset, see Tables S1 and S2 in supplementary material). Based on the Protein Pilot™ algorithm [27], 61 proteins were present in significantly changed concentrations (P<0.01) in the M. mycoides subsp. capri ctpA(clpP):tetM mutant when compared to the wild type (Table 2; Table S2). In order to derive biological meaning from these differences, proteins were grouped according to global biologic functions as assigned in the Molligen 3.0 database [28] (http://cbi.labri.fr/outils/molligen) or UNIPROT [29], [30] (www.uniprot.org/).

Table 2

Proteins that were significantly altered by the disruption of MMCAP2_0241 in Mycoplasma mycoides subsp. mycoides as determined by iTRAQ™.

An interesting feature in the M. mycoides subsp. capri ctpA(clpP)::tetM mutant is that most of the significantly altered proteins involved in metabolism were decreased, with the exception of putrescine carbamoyltransferase. Conversely, proteins that were significantly increased in the mutant are associated with genetic and environmental information processing. Proteins with an unassigned unknown function also tended to be decreased. Many of the proteins involved in glycolysis, including LDH, had reduced levels in the mutant. We chose to determine functional activity of LDH because this enzyme is required to convert pyruvate to lactate, thereby recycling NADH. Because mycoplasma lacks the electron respiratory chain, recycling NADH is a key component to maintaining redox potential in the microbe. The M. mycoides subsp. capri ctpA(clpP)::tetM mutant exhibited about a 50% reduction in LDH activity in comparison to the wild type strain, P<0.0001 (Fig. 3). Thus, the observed functional activity was in agreement with the iTRAQ data for LDH.

Figure 3

Lactate dehydrogenase production in the Mycoplasma mycoides subsp. capri GM12 wild type and the ctpA(clpP)::tetM mutant.

Discussion

Ctps are a group of heterologous serine proteases characterized by resistance to conventional protease inhibitors and a catalytic center that differs from classical serine proteases [31], [32], [33]. Initial studies in bacteria demonstrated that Ctp proteins appear to be multifunctional [34], may be important for degrading damaged or aberrant proteins [35], coping with environmental stress [36], impacting known virulence factors of bacteria [37], [38], [39], [40], [41], and modulating host response to infection [39], [41]. Because Ctp selectively targets and cleaves the non-polar C-terminal of many precursor proteins [40], these enzymes play a role in protein maturation, protein transport into other cellular organelles, or export to the periplasm of Gram negative bacteria [42], [43]. The ClpPs are a S41 peptidase superfamily, also containing the TSP domain. In Gram positive bacteria with low GC content [19], ClpPs are essential in removing heat-damaged proteins. In addition, clpP deletion mutants exhibit diverse and pleiotropic effects, ranging from impacts on general stress response, metabolism, sensitivity to high temperatures, metal ion transport, and virulence [19], [20], [22], [23], [44]. In our current study, the disruption of the MMCAP2_0241 gene clearly had diverse, pleiotropic effect, including changes in the proteome profile, decreased proteins involved in metabolism, decreased LDH activity, increased H₂O₂ production, and decreased thermal stability. These pleiotropic changes are, for the most part, consistent with the observed effects of other microbial Ctp and ClpP proteins. One notable exception was that in our mutant, production of H₂O₂ was increased, whereas in clpP mutants increased sensitivity to H₂O₂ was observed [22], [44].

In part because of their limited genome size, mycoplasmas have restricted catabolic pathways, and many of the pathways commonly present in other bacteria are missing or truncated in mycoplasmas [45]. Although glycolysis is the most common central pathway, glycerol also can be used by many mycoplasmas as a carbon and energy source. Glycerol is transported into mycoplasmas, phosphorylated, and then metabolized to dihyroxyacetone phosphate with the concomitant release of H₂O₂ [14], [45]. Glycerol oxidation and the subsequent release of H₂O₂ plays a central role in the cytotoxicity of M. mycoides subsp. mycoides SC and the M. mycoides cluster [10], [11], [12], [14], [15] as well as in M. pneumonia [9], [46]. Increased production of H₂O₂ also has been associated with increased virulence and toxicity in other ruminant mycoplasmas, including M. agalactiae [13], M. arginini [47], M. bovis [13], and M. ovipneumoniae [47]. European strains of M. mycoides subsp. mycoides SC are significantly less virulent than African strains and have reduced H₂O₂ production [12], [14], [48]. Genomic analysis [14], [48], [49] has shown that the European strains lack an 8.8 Kb DNA segment in the gtsABC operon and are less efficient at glycerol uptake, thus explaining the reduced H₂O₂ production. A similar gtsABC operon was described in M. leachii [15], which produces levels of H₂O₂ comparable to African strains of M. mycoides subsp. mycoides SC. However, the production of H₂O₂ by itself is not enough to elicit the cytotoxic effect. For example, vaccine strains of M. mycoides subsp. mycoides SC are able to produce H₂O₂ comparable to virulent strains, but have attenuated virulence [10], suggesting that additional virulence factors are required. Based on in vitro studies with adhesin-deficient mutants [10], [12], it is likely that strong contact between the mycoplasma and the host cell is required for H₂O₂ toxicity, and that this direct contact facilitates the translocation of H₂O₂ into the cytoplasm of the host cell. Thus, both adherence to the host cell surface as well as the release of H₂O₂ is required for toxicity.

The disruption of MMCAP2_0241 in M. mycoides subsp. capri impacted a number of proteins, including LDH. The M. mycoides subsp. capri ctpA(clpP)::tetM mutant had a significant reduction (about 50%) in LDH activity. Although the full biological impact of this reduced activity is not known, it could potentially contribute to oxidative stress. Under normal conditions, M. mycoides subsp. capri catabolizes sugars to pyruvate via glycolysis [11], [45], [50]. Pyruvate can be reduced to lactate by LDH, which results in reoxidation of NADH to yield NAD⁺ [45]. Because these microbes do not possess an electron transport chain that can be used for this purpose, the recycling of NADH is critical to maintain cellular redox balance [51] and likely to be a crucial point for the adjustment of mycoplasmal metabolism [51], [52].

Another potential source of oxidative stress is reactive oxygen species (ROS) like H₂O₂. We observed an approximately 40% relative increase in production of H₂O₂ in the mutant. Increased H₂O₂ may result in damage not only to the host but also the microbe via oxidation of macromolecules like proteins and lipids. Although most mycoplasmas are deficient in superoxide dismutase and catalase [53], both thioredoxin (TrxA) and methionine sulfoxide reductase (Msr) are present and can modulate oxidative stress and ameliorate the toxic effects of ROS [53], [54]. Interestingly, however, significant decreases as measured by iTRAQ™ were observed in the M. mycoides subsp. capri ctpA(clpP)::tetM mutant for both Msr (37% decrease) and Trx (65%) proteins, suggesting that the mutant might be under oxidative stress and less able to ameliorate the toxic effects of ROS. The observed increases in H₂O₂ levels could be the result of metabolic changes. Most notably, the twofold increase in dihydroxyacetone kinase in the ctpA(clpP)::tetM mutant and the decrease in several enzymes in the glycolytic pathway are suggestive of a shift to glycerol metabolism. Because catalase and peroxyredoxin are not present in the genome sequence of M. mycoides subsp. capri, it is highly unlikely that the increased levels of H₂O₂ we observed are due to a decrease in these activities.

The upregulation of ribosomal protein expression under stress conditions has previously been reported in mycoplasmas. The increased expression of ribosomal proteins may reflect the need to increase the translation and protein synthesis. Since the ctpA(clpP)::tetM mutant is more susceptible to heat shock and more exposed to oxidative stress, it may experience more protein damage and misfolding. This is consistent with the observed upregulation of ribosomal proteins in response to heat shock in both M. hyopneumoniae [55] and M. pneumoniae [56]. Additionally, loss of the ctpA/clpP protein could impact the ability of the mutant to degrade damaged or aberrant proteins [35]. Thus, the increased translation could be a compensatory mechanism to cope with these consequences.

Our M. mycoides subsp. capri ctpA(clpP)::tetM mutant, like the M. mycoides subsp. mycoides SC, was sensitive to heat shock [8]. Although the mechanism by which MMCAP2_0241 (ClpP-like) is involved in thermal stability in M. mycoides subsp. capri is unknown, it is interesting that disruption of an E. coli gene with a similar TSP-processing domain also resulted in sensitivity to high temperature [57], suggesting that the observed effect might be a result of loss of protein processing.

In our current study, the disruption of the MMCAP2_0241 gene clearly had a pleiotropic effect. Based on our results, it is reasonable to suggest that MMCAP2_0241 (ClpP-like) plays a role in stress response in M. mycoides subsp. capri. Although the specific mechanisms are not known and may be direct or indirect, it is important to note that by addressing changes in the proteome profile, we were able to identify proteins of interest that would not be predicted based on the loss of the proteolytic phenotype alone. In addition to the loss of the proteolytic phenotype, the mutant also exhibited reduced LDH activity, increased H₂O₂ production, and increased susceptibility to heat stress. The increased H₂O₂ production is particularly intriguing, as it may have implications for virulence in mycoplasmas [9], [10], [12], [13], [14], [46], [47] and also is a known source of oxidative stress in other bacteria [58].

Materials and Methods

Mycoplasma strains and cultivation

M. mycoides subsp. capri GM12 type ATCC 35297 wild type [1] and the M. mycoides subsp. capri ctpA(clpP)::tetM mutant that was generated by double cross-over homologous recombination [16] were used in this study. M. mycoides subsp. capri GM12 type ATCC 35297 has been fully sequenced by The J. Craig Venter Institute, and the full sequences (both gene and protein) are available at GenBank (http://www.ncbi.nlm.nih.gov/bioproject?term=PRJNA39245) and via Genomes on Line (http://www.genomesonline.org/cgi-bin/GOLD/index.cgi). For all experiments, both strains were cultivated in parallel at 37°C in the same batch of SP4 medium, with the exception that the medium was supplemented with 10 µg/ml of tetracycline for growth of the mutant. Microbial growth was monitored by optical density at 640 nm, and cultures were harvested at late log phase (OD₆₄₀

0.08). All cultures contained 10⁹ CFU per ml of medium, which was confirmed by direct colony counts. For proteomic and gel electrophoresis studies, cultures were concentrated by centrifugation and used at a final concentration of 10¹² CFU per ml.

Preparation of protein extracts for proteomics

Bacterial suspensions were divided into two aliquots: one aliquot was used for 2D gel electrophoresis and the second was used for iTRAQ™ followed by tandem mass spectrometry. Bacterial suspensions were pelleted by centrifugation and washed with wash buffer solution (Calbiochem ProteoExtract® kit, San Diego, CA). For 2D gel electrophoresis, proteins were extracted with Trizol (Invitrogen Corp., Carlsbad, CA) according to the manufacturer's protocol. Pelleted protein extracts were allowed to air dry and were stored at −20°C. For iTRAQ™ analysis, proteins were extracted with ProteoExtract® Complete Mammalian Proteome Extraction Kit (Calbiochem, San Diego, CA) according to the manufacturer's protocol. The total protein concentration of all samples was determined with the Non-Interfering Protein Assay™ Kit (Calbiochem, San Diego, CA).

2-Dimensional differential gel electrophoresis (2D-DIGE) and protein identification

Protein pellets were prepared for 2D electrophoresis as previously described [59]. Three 2-D electrophoresis experiments were performed with samples obtained from three independent experiments. The Cy2 internal standard was used to co-detect, match and normalize protein spots in all three gels. Gel images were obtained with Typhoon 9600 Variable Mode Imager (GE Healthcare, Piscataway, NJ) and images were analyzed with DeCyder 2D version 7.0 software (GE Healthcare, Piscataway, NJ). Protein ratios for each gel spot were generated by dividing the total area of M. mycoides subsp. capri ctpA(clpP)::tetM spot by the total area of the corresponding wild type M. mycoides subsp. capri spot. Only protein ratios that were 2-fold or greater in difference were considered for further analysis. An automated spot picker (ProPic Workstation, Digilab Genomic Solutions Inc., Ann Arbor, MI) selected protein targets for identification. The same protein spot from each biological replicate was pooled for processing and identification by tandem mass spectrometry as previously described [59]. Tandem mass spectrometric data was searched against NCBI nr bacterial database using Mascot (Matrix Science, Boston, MA) database search engine. Protein identification was performed with Scaffold version 2_01_02 (Proteome Software Inc., Portland, OR). Protein identifications were accepted if they could be established at greater than 80.0% probability and contained at least one identified peptide. Protein probabilities were assigned by the Protein Prophet algorithm [60]. Proteins that contained similar peptides and could not be differentiated based on MS/MS analysis alone were grouped to satisfy the principles of parsimony.

Quantitative proteomics using peptide labeling and 2D-LC-MS/MS

In order to minimize variability, the protein extracts from all three biological replicates of wild type M. mycoides subsp. capri (control) were combined and the total protein concentration of the pooled sample was adjusted to match the total protein concentration of each biological replicate of M. mycoides subsp. capri ctpA(clpP)::tetM mutant. Bacterial protein extracts were prepared as previously described [59] and labeled with an amine specific peptide-based labeling system, iTRAQ™, according to the manufacturer's instructions (Applied Biosystems, Foster City, CA). Labeled samples were analyzed inline with a hybrid quadrupole-TOF mass spectrometer QSTAR (Applied Biosystems Inc) as previously described [59]. Tandem mass spectra were extracted by Analyst (v 1.1.; Applied Biosystem Inc) and the NCBI bacterial protein database (concatenation of the forward and random sequences) was used for protein identification. Searches were performed using MS/MS data interpretation algorithms from Protein Pilot™ (Paragon™ algorithm, v 3.0, Applied Biosystem Inc) and Mascot (v 2.2, Matrix Science, London, UK) as already described [59]. Quantification of protein ratios was based on a minimum of three spectra. Only protein ratios with an error factor (EF) <2 were retained for further analysis. EF is a measure of the variation among the different iTRAQ™ ratios (the greater the variation, the greater the uncertainty) and represents the 95% uncertainty range for a reported ratio. Proteins were grouped according to global biologic function as assigned in the Molligen 2.0 database [28] and/or UniProt [29], [30]. Protein ratios were considered significantly different if they had P values<0.01 as determined by the Pro Group™ algorithm [27].

Carbohydrate use

We compared the carbohydrate use of M. mycoides subsp. capri and the ctpA(clpP)::tetM mutant using the BBL Crystal™ system (Becton Dickinson Microbiology Systems, Cockeysville, Md.) according to the manufacturer's instructions. Both strains were grown to midlogarithmic growth in 2 ml of SP4 broth. Cells were pelleted, suspended in the inoculation medium, the kit test plate inoculated and incubated at 37°C for 24 hr before reading results. Tests were performed in triplicate.

Response to heat shock

Both wild type and mutant strains were grown to midlogarithmic phase at 37°C as described above. Each culture was divided into two aliquots. One aliquot was kept at 37°C and the other was immediately transferred to 42°C. Each treatment was done in triplicate. At selected time points, an aliquot was removed, serially diluted, and the CFU determined by direct colony count.

Measurement of H₂O₂ production

H₂O₂ production was determined as previously described for M. mycoides subsp. mycoides SC [10]. Six replicates were performed for each strain. Briefly, wild type and mutant M. mycoides subsp. capri were grown to midlogarithmic phase and harvested by centrifugation at 10,000 rpm for 15 min at 4°C. Bacterial pellets were washed 3 times in incubation buffer (67.7 mM HEPES, pH 7.3; 140 mM NaCl, and 7 mM MgCl₂). After the final wash, bacteria were resuspended in incubation buffer to a cell density of 10⁸ cells/ml and incubated for 20 min at 37°C. Immediately after the incubation period, glycerol was added to give a final concentration of 100 µM. H₂O₂ production was measured at 30 min after glycerol addition using an H₂O₂ assay kit (Cayman Chemical, Ann Arbor, MI). The assay, standard curve determination, and specificity quality controls were performed according to the instructional manual.

Measurement of lactate dehydrogenase activity

Lactate dehydrogenase (LDH) activity and concentration were measured using the Lactate Dehydrogenase Colorimetric Assay Kit (Abcam Inc., Cambridge, MA). Five replicates were performed for each strain. Briefly, 4 ml of both wild type and mutant M. mycoides subsp. capri were grown to midlogarithmic phase (10⁶ cells/ml) and harvested by centrifugation as described above. The cell pellets were homogenized in 0.5 ml of cold assay buffer, centrifuged at 10,000 rpm for 15 minutes at 4°C, and the supernatant collected. Twenty µl of the supernatant was used in a total 200 µl reaction. A standard curve was constructed, and the LDH activity and concentration were measured and calculated according to the manufacturer instructions.

Statistical analysis

The 2D gel data, H₂O₂ production data, LDH activity, and heat shock growth data were analyzed by unpaired students T test or ANOVA. CFU data was log transformed prior to analysis by ANOVA. For statistical analysis of iTRAQ™ data, protein ratios were generated with Pro Group™ algorithm and automatically corrected for bias. The calculated P-value obtained with the ProGroup™ algorithm is based on 95% confidence interval. A P value<0.05 was accepted as significant.

Supporting Information

Figure S1

Alignment and conserved signature sequences of M. mycoides subsp. capri (MMCAP2_0241) and M. capricolum (MCAP_0240) compared with representative bacterial CTP proteins. Bacillus subtilis (ZP_03591706.1), Bartonella bacilliformis (YP_988644.1), Burkholderia pseudomallei (YP_107067.1), Escherichia coli (Prc; M75634.1), Legionella pneumophila (YP_122899.1), Neisseria gonorrheae (ZP_04720812.1), Synechococcus sp. (NP_898059.1), and Staphylococcus aureus (BAB42513.1) were used for comparative purposes. Identical and similar amino acid (aa) residues are marked red and blue, respectively. For M. mycoides capri, note the signal peptide domain (aa 1–24) followed by the N-terminal domain (aa 25–339), two transmembrane domains (aa 7–29 and aa 631–650) and the S41 peptidase tail specific protease (TSP) domain (aa 340–544). Within the TSP domain and denoted by yellow boxes are the two carboxyl-tail processing (Ctp) motifs (aa 411–422 and aa 477–507). Protein sequences were aligned using Clustal Omega (http://www.clustal.org/omega) followed by analysis using the Sequence Manipulation suite: multiple align show (www.bioinformatics.org/sms). Highly conserved, identical amino acids within the TSP domain are highlighted in red; highly similar amino acids are highlighted in blue. Note that both the amino- and carboxyl- terminal portions of the proteins are quite diverse.

(TIFF)

Click here for additional data file.^{(12M, tiff)}

Figure S2

Differential 2-dimensional electrophoresis of M. mycoides subsp. capri GM12 wild type and M. mycoides capri ctpA::tetM mutant. A. The distribution of proteins in the wild type (Cy3 label, red) and mutant (Cy5 label, green). B. The distribution of protein spots that were significantly increased in the mutant. C. The distribution of protein spots that were significantly increased in M. mycoides capri.

(TIFF)

Click here for additional data file.^{(375K, tiff)}

Table S1

Proteins that were not significantly altered by the disruption of ctpA in Mycoplasma mycoides subsp. capri GM12. The gene locus and gene name, if annotated, were identified using the Molligen database and the M. mycoides subsp. capri GM12, taxon:436113 genome sequence. GI accession numbers refer to GenBank (http://www.ncbi.nlm.nih.gov/genbank). The protein sequences are available at GenBank (http://www.ncbi.nlm.nih.gov/bioproject?term=PRJNA39245) and via Genomes on Line (http://www.genomesonline.org/cgi-bin/GOLD/index.cgi). Data is shown for three independent biological replicates. Protein ratios were generated by dividing the spectral intensity of a specific protein in the M. mycoides subsp. capri ctpA (clpP)::tetM mutant by the spectral intensity of the specific protein in M. mycoides subsp. capri GM12. Protein ratios were calculated with the Pro Group™ algorithm (Applied Biosystems/MDS SCIEX). Only ratios from the spectra that are distinct to each protein (or protein form) were used for the calculation. The total protein score is the measurement of all the peptide evidence for a protein and is analogous to protein scores reported by other protein identification software. Each identified peptide within a protein was assigned a score based on confidence (95% confidence

2, 95% confidence

1.3). Raw peptide identification was performed using the Paragon™ database searching algorithm (Applied Biosystems/MDS SCIEX). Raw peptide identification was further processed with the Pro Group Algorithm™ (Applied Biosystems/MDS SCIEX). Only proteins with 95% or greater confidence and an EF factor <2 were considered. Proteins were decreased in the mutant if the ratio was <1.0; proteins were increased in the mutant if the ratio was >1.0. Proteins were considered to be unchanged in the mutant if the P value was >0.001. Further, significance levels had to be obtained in all 3 biological replicates.

(XLSX)

Click here for additional data file.^{(92K, xlsx)}

Table S2

Proteins that were significantly altered by the disruption of ctpA in M. mycoides subsp. capri GM12. The gene locus and gene name, if annotated, were identified using the Molligen database and the M. mycoides subsp. capri GM12, taxon:436113 genome sequence. GI accession numbers refer to GenBank (http://www.ncbi.nlm.nih.gov/genbank). The protein sequences are available at GenBank (http://www.ncbi.nlm.nih.gov/bioproject?term=PRJNA39245) and via Genomes on Line (http://www.genomesonline.org/cgi-bin/GOLD/index.cgi). Data is shown for three independent biological replicates. Protein ratios were generated by dividing the spectral intensity of a specific protein in the M. mycoides subsp. capri ctpA (clpP)::tetM mutant by the spectral intensity of the specific protein in M. mycoides subsp. capri GM12. Protein ratios were calculated with the Pro Group™ algorithm (Applied Biosystems/MDS SCIEX). Only ratios from the spectra that are distinct to each protein (or protein form) were used for the calculation. The total protein score is the measurement of all the peptide evidence for a protein and is analogous to protein scores reported by other protein identification software. Each identified peptide within a protein was assigned a score based on confidence (95% confidence

2, 95% confidence

1.3). Raw peptide identification was performed using the Paragon™ database searching algorithm (Applied Biosystems/MDS SCIEX). Raw peptide identification was further processed with the Pro Group Algorithm™ (Applied Biosystems/MDS SCIEX). Only proteins with 95% or greater confidence and an EF factor <2 were considered. Proteins were decreased in the mutant if the ratio was <1.0; proteins were increased in the mutant if the ratio was >1.0. Proteins were significantly different in the mutant if the P value was <0.001. Further, significance levels had to be obtained in all 3 biological replicates.

(XLSX)

Click here for additional data file.^{(73K, xlsx)}

Acknowledgments

Portions of the work were funded by USDA CSREES grant 99-35204-8367 to M.B.B. and subcontract 6415-100607 to M.B.B. from the University of South Florida Center for Biological Defense (contract W911SR-06-C-0023). L.R. was supported in part by Public Health Service grant 5K08DK07651 from the National Institute of Diabetes, Digestive and Kidney Diseases.

Funding Statement

The majority of the work was funded by USDA CSREES grant 99-35204-8367 to M.B.B. and subcontract 6415-100607 to M.B.B. from the University of South Florida Center for Biological Defense (contract W911SR-06-C-0023). L.R. was supported by Public Health Service grant 5K08DK07651 from the National Institute of Diabetes and Digestive and Kidney Diseases. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. DaMassa AJ, Brooks DL, Adler HE (1983) Caprine mycoplasmosis: widespread infection in goats with Mycoplasma mycoides subsp mycoides (large-colony type). Am J Vet Res 44: 322–325. [PubMed]

2. Rodriguez JL, Poveda JB, Oros J, Herraez P, Sierra MA, et al. (1995) High mortality in goats associated with the isolation of a strain of Mycoplasma mycoides subsp. mycoides (Large Colony type). Zentralbl Veterinarmed B 42: 587–593. [PubMed]

3. Manso-Silvan L, Perrier X, Thiaucourt F (2007) Phylogeny of the Mycoplasma mycoides cluster based on analysis of five conserved protein-coding sequences and possible implications for the taxonomy of the group. Int J Syst Evol Microbiol 57: 2247–2258. [PubMed]

4. Manso-Silvan L, Vilei EM, Sachse K, Djordjevic SP, Thiaucourt F, et al. (2009) Mycoplasma leachii sp. nov. as a new species designation for Mycoplasma sp. bovine group 7 of Leach, and reclassification of Mycoplasma mycoides subsp. mycoides LC as a serovar of Mycoplasma mycoides subsp. capri. Int J Syst Evol Microbiol 59: 1353–1358. [PubMed]

5. Abu-Groun EA, Taylor RR, Varsani H, Wadher BJ, Leach RH, et al. (1994) Biochemical diversity within the “Mycoplasma mycoides” cluster. Microbiology 140 (Pt 8) 2033–2042. [PubMed]

6. Cottew GS, Breard A, DaMassa AJ, Erno H, Leach RH, et al. (1987) Taxonomy of the Mycoplasma mycoides cluster. Isr J Med Sci 23: 632–635. [PubMed]

7. Cottew GS, Yeats FR (1978) Subdivision of Mycoplasma mycoides subsp. mycoides from cattle and goats into two types. Aust Vet J 54: 293–296. [PubMed]

8. Thiaucourt F, Manso-Silvan L, Salah W, Barbe V, Vacherie B, et al. (2011) Mycoplasma mycoides, from “mycoides Small Colony” to “capri”. A microevolutionary perspective. BMC Genomics 12: 114. [PMC free article] [PubMed]

9. Hames C, Halbedel S, Hoppert M, Frey J, Stulke J (2009) Glycerol metabolism is important for cytotoxicity of Mycoplasma pneumoniae. J Bacteriol 191: 747–753. [PMC free article] [PubMed]

10. Bischof DF, Janis C, Vilei EM, Bertoni G, Frey J (2008) Cytotoxicity of Mycoplasma mycoides subsp. mycoides small colony type to bovine epithelial cells. Infect Immun 76: 263–269. [PMC free article] [PubMed]

11. Vilei EM, Correia I, Ferronha MH, Bischof DF, Frey J (2007) Beta-D-glucoside utilization by Mycoplasma mycoides subsp. mycoides SC: possible involvement in the control of cytotoxicity towards bovine lung cells. BMC Microbiol 7: 31. [PMC free article] [PubMed]

12. Pilo P, Vilei EM, Peterhans E, Bonvin-Klotz L, Stoffel MH, et al. (2005) A metabolic enzyme as a primary virulence factor of Mycoplasma mycoides subsp. mycoides small colony. J Bacteriol 187: 6824–6831. [PMC free article] [PubMed]

13. Khan LA, Miles RJ, Nicholas RA (2005) Hydrogen peroxide production by Mycoplasma bovis and Mycoplasma agalactiae and effect of in vitro passage on a Mycoplasma bovis strain producing high levels of H2O2. Vet Res Commun 29: 181–188. [PubMed]

14. Vilei EM, Frey J (2001) Genetic and biochemical characterization of glycerol uptake in mycoplasma mycoides subsp. mycoides SC: its impact on H(2)O(2) production and virulence. Clin Diagn Lab Immunol 8: 85–92. [PMC free article] [PubMed]

15. Djordjevic SP, Vilei EM, Frey J (2003) Characterization of a chromosomal region of Mycoplasma sp. bovine group 7 strain PG50 encoding a glycerol transport locus (gtsABC). Microbiology 149: 195–204. [PubMed]

16. Allam AB, Reyes L, Assad-Garcia N, Glass JI, Brown MB (2010) Enhancement of targeted homologous recombination in Mycoplasma mycoides subsp. capri by inclusion of heterologous recA. Appl Environ Microbiol 76: 6951–6954. [PMC free article] [PubMed]

17. Kelley LA, Sternberg MJ (2009) Protein structure prediction on the Web: a case study using the Phyre server. Nature Protocols 4: 363–371. [PubMed]

18. Maurizi MR, Clark WP, Kim SH, Gottesman S (1990) Clp P represents a unique family of serine proteases. J Biol Chem 265: 12546–12552. [PubMed]

19. Frees D, Savijoki K, Varmanen P, Ingmer H (2007) Clp ATPases and ClpP proteolytic complexes regulate vital biological processes in low GC, Gram-positive bacteria. Mol Microbiol 63: 1285–1295. [PubMed]

20. Chattoraj P, Banerjee A, Biswas S, Biswas I (2010) ClpP of Streptococcus mutans differentially regulates expression of genomic islands, mutacin production, and antibiotic tolerance. J Bacteriol 192: 1312–1323. [PMC free article] [PubMed]

21. Ostberg Y, Carroll JA, Pinne M, Krum JG, Rosa P, et al. (2004) Pleiotropic effects of inactivating a carboxyl-terminal protease, CtpA, in Borrelia burgdorferi. J Bacteriol 186: 2074–2084. [PMC free article] [PubMed]

22. Robertson GT, Ng WL, Foley J, Gilmour R, Winkler ME (2002) Global transcriptional analysis of clpP mutations of type 2 Streptococcus pneumoniae and their effects on physiology and virulence. J Bacteriol 184: 3508–3520. [PMC free article] [PubMed]

23. Zhang Q, Huang Y, Wang H, Xu W, Liu L, et al. (2012) The role of ClpP in protein expression of Streptococcus pneumoniae. Current microbiology 64: 294–299. [PubMed]

24. Gorg A, Weiss W, Dunn MJ (2004) Current two-dimensional electrophoresis technology for proteomics. Proteomics 4: 3665–3685. [PubMed]

25. Lilley KS, Friedman DB (2004) All about DIGE: quantification technology for differential-display 2D-gel proteomics. Expert Rev proteomics 1: 401–409. [PubMed]

26. Thon JN, Schubert P, Duguay M, Serrano K, Lin S, et al. (2008) Comprehensive proteomic analysis of protein changes during platelet storage requires complementary proteomic approaches. Transfusion 48: 425–435. [PubMed]

27. Shilov IV, Seymour SL, Patel AA, Loboda A, Tang WH, et al. (2007) The Paragon Algorithm, a next generation search engine that uses sequence temperature values and feature probabilities to identify peptides from tandem mass spectra. Mol Cell Proteomics 6: 1638–1655. [PubMed]

28. Barre A, de Daruvar A, Blanchard A (2004) MolliGen, a database dedicated to the comparative genomics of Mollicutes. Nucleic Acids Res 32: D307–310. [PMC free article] [PubMed]

29. Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, et al. (2004) UniProt: the Universal Protein knowledgebase. Nucleic acids research 32: D115–119. [PMC free article] [PubMed]

30. Dimmer EC, Huntley RP, Alam-Faruque Y, Sawford T, O'Donovan C, et al. (2012) The UniProt-GO Annotation database in 2011. Nucleic acids research 40: D565–570. [PMC free article] [PubMed]

31. Ekici OD, Paetzel M, Dalbey RE (2008) Unconventional serine proteases: variations on the catalytic Ser/His/Asp triad configuration. Protein Sci 17: 2023–2037. [PubMed]

32. Liao DI, Qian J, Chisholm DA, Jordan DB, Diner BA (2000) Crystal structures of the photosystem II D1 C-terminal processing protease. Nat Struct Biol 7: 749–753. [PubMed]

33. Paetzel M, Dalbey RE (1997) Catalytic hydroxyl/amine dyads within serine proteases. Trends Biochem Sci 22: 28–31. [PubMed]

34. Baumler AJ, Kusters JG, Stojiljkovic I, Heffron F (1994) Salmonella typhimurium loci involved in survival within macrophages. Infect Immun 62: 1623–1630. [PMC free article] [PubMed]

35. Keiler KC, Waller PR, Sauer RT (1996) Role of a peptide tagging system in degradation of proteins synthesized from damaged messenger RNA. Science 271: 990–993. [PubMed]

36. Gilbert KB, Vanderlinde EM, Yost CK (2007) Mutagenesis of the carboxy terminal protease CtpA decreases desiccation tolerance in Rhizobium leguminosarum. FEMS Microbiol Lett 272: 65–74. [PubMed]

37. Bandara AB, DeShazer D, Inzana TJ, Sriranganathan N, Schurig GG, et al. (2008) A disruption of ctpA encoding carboxy-terminal protease attenuates Burkholderia mallei and induces partial protection in CD1 mice. Microb Pathog 45: 207–216. [PubMed]

38. Bandara AB, Sriranganathan N, Schurig GG, Boyle SM (2005) Carboxyl-terminal protease regulates Brucella suis morphology in culture and persistence in macrophages and mice. J Bacteriol 187: 5767–5775. [PMC free article] [PubMed]

39. Lad SP, Yang G, Scott DA, Wang G, Nair P, et al. (2007) Chlamydial CT441 is a PDZ domain-containing tail-specific protease that interferes with the NF-kappaB pathway of immune response. J Bacteriol 189: 6619–6625. [PMC free article] [PubMed]

40. Silber KR, Sauer RT (1994) Deletion of the prc (tsp) gene provides evidence for additional tail-specific proteolytic activity in Escherichia coli K-12. Mol Gen Genet 242: 237–240. [PubMed]

41. Zhong G, Fan P, Ji H, Dong F, Huang Y (2001) Identification of a chlamydial protease-like activity factor responsible for the degradation of host transcription factors. J Exp Med 193: 935–942. [PMC free article] [PubMed]

42. Diner BA, Ries DF, Cohen BN, Metz JG (1988) COOH-terminal processing of polypeptide D1 of the photosystem II reaction center of Scenedesmus obliquus is necessary for the assembly of the oxygen-evolving complex. J Biol Chem 263: 8972–8980. [PubMed]

43. Nagasawa H, Sakagami Y, Suzuki A, Suzuki H, Hara H, et al. (1989) Determination of the cleavage site involved in C-terminal processing of penicillin-binding protein 3 of Escherichia coli. J Bacteriol 171: 5890–5893. [PMC free article] [PubMed]

44. Kwon HY, Kim SW, Choi MH, Ogunniyi AD, Paton JC, et al. (2003) Effect of heat shock and mutations in ClpL and ClpP on virulence gene expression in Streptococcus pneumoniae. Infect Immun 71: 3757–3765. [PMC free article] [PubMed]

45. Halbedel S, Hames C, Stulke J (2007) Regulation of carbon metabolism in the mollicutes and its relation to virulence. J Mol Microbiol Biotechnol 12: 147–154. [PubMed]

46. Schmidl SR, Otto A, Lluch-Senar M, Pinol J, Busse J, et al. (2011) A trigger enzyme in Mycoplasma pneumoniae: impact of the glycerophosphodiesterase GlpQ on virulence and gene expression. PLoS pathogens 7: e1002263. [PMC free article] [PubMed]

47. Niang M, Rosenbusch RF, DeBey MC, Niyo Y, Andrews JJ, et al. (1998) Field isolates of Mycoplasma ovipneumoniae exhibit distinct cytopathic effects in ovine tracheal organ cultures. Zentralbl Veterinarmed A 45: 29–40. [PubMed]

48. Vilei EM, Frey J (2001) Genetic and biochemical characterization of glycerol uptake in mycoplasma mycoides subsp. mycoides SC: its impact on H(2)O(2) production and virulence. Clin Diagn Lab Immunol 8: 85–92. [PMC free article] [PubMed]

49. Vilei EM, Abdo EM, Nicolet J, Botelho A, Goncalves R, et al. (2000) Genomic and antigenic differences between the European and African/Australian clusters of Mycoplasma mycoides subsp. mycoides SC. Microbiology 146 (Pt 2) 477–486. [PubMed]

50. Miles RJ (1992) Catabolism in mollicutes. J Gen Microbiol 138: 1773–1783. [PubMed]

51. Pollack JD, Merola AJ, Platz M, Booth RL Jr (1981) Respiration-associated components of Mollicutes. J Bacteriol 146: 907–913. [PMC free article] [PubMed]

52. Kannan TR, Baseman JB (2000) Hemolytic and hemoxidative activities in Mycoplasma penetrans. Infect Immun 68: 6419–6422. [PMC free article] [PubMed]

53. Ben-Menachem G, Himmelreich R, Herrmann R, Aharonowitz Y, Rottem S (1997) The thioredoxin reductase system of mycoplasmas. Microbiology 143 (Pt 6) 1933–1940. [PubMed]

54. Pollack JD, Williams MV, McElhaney RN (1997) The comparative metabolism of the mollicutes (Mycoplasmas): the utility for taxonomic classification and the relationship of putative gene annotation and phylogeny to enzymatic function in the smallest free-living cells. Crit Rev Microbiol 23: 269–354. [PubMed]

55. Madsen ML, Nettleton D, Thacker EL, Edwards R, Minion FC (2006) Transcriptional profiling of Mycoplasma hyopneumoniae during heat shock using microarrays. Infect Immun 74: 160–166. [PMC free article] [PubMed]

56. Weiner J 3rd, Zimmerman CU, Gohlmann HW, Herrmann R (2003) Transcription profiles of the bacterium Mycoplasma pneumoniae grown at different temperatures. Nucleic Acids Res 31: 6306–6320. [PMC free article] [PubMed]

57. Hara H, Yamamoto Y, Higashitani A, Suzuki H, Nishimura Y (1991) Cloning, mapping, and characterization of the Escherichia coli prc gene, which is involved in C-terminal processing of penicillin-binding protein 3. J Bacteriol 173: 4799–4813. [PMC free article] [PubMed]

58. Mishra S, Imlay J (2012) Why do bacteria use so many enzymes to scavenge hydrogen peroxide? Arch Biochem Biophys 525: 145–160. [PubMed]

59. Reyes L, Alvarez S, Allam A, Reinhard M, Brown MB (2009) Complicated urinary tract infection is associated with uroepithelial expression of proinflammatory protein S100A8. Infect Immun 77: 4265–4274. [PMC free article] [PubMed]

60. Nesvizhskii AI, Keller A, Kolker E, Aebersold R (2003) A statistical model for identifying proteins by tandem mass spectrometry. Anal Chem 75: 4646–4658. [PubMed]

Articles from PLoS ONE are provided here courtesy of
Public Library of Science