BioMed Research International

BioMed Research International / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 4165735 | 18 pages |

Quantitative Proteomic Analyses of a Pathogenic Strain and Its Highly Passaged Attenuated Strain of Mycoplasma hyopneumoniae

Academic Editor: Clara G. de los Reyes-Gavilan
Received17 Dec 2018
Revised14 Apr 2019
Accepted27 May 2019
Published01 Jul 2019


Mycoplasma hyopneumoniae is the causative agent of porcine enzootic pneumonia, a chronic respiratory disease in swine resulting in enormous economic losses. To identify the components that contribute to virulence and unveil those biological processes potentially related to attenuation, we used isobaric tags for relative and absolute quantification technology (iTRAQ) to compare the protein profiles of the virulent M. hyopneumoniae strain 168 and its attenuated highly passaged strain 168L. We identified 489 proteins in total, 70 of which showing significant differences in level of expression between the two strains. Remarkably, proteins participating in inositol phosphate metabolism were significantly downregulated in the virulent strain, while some proteins involved in nucleoside metabolism were upregulated. We also mined a series of novel promising virulence-associated factors in our study compared with those in previous reports, such as some moonlighting adhesins, transporters, lipoate-protein ligase, and ribonuclease and several hypothetical proteins with conserved functional domains, deserving further research. Our survey constitutes an iTRAQ-based comparative proteomic analysis of a virulent M. hyopneumoniae strain and its attenuated strain originating from a single parent with a well-characterized genetic background and lays the groundwork for future work to mine for potential virulence factors and identify candidate vaccine proteins.

1. Introduction

M. hyopneumoniae is the causative agent of porcine enzootic pneumonia, which is a worldwide epidemic that can cause enormous economic losses as a result of retarded growth in pigs and the cost of disease control [1]. Despite its low direct mortality, M. hyopneumoniae increases the host’s susceptibility to secondary respiratory infections by damaging cilia and epithelial cell, resulting in aggravated lung lesions and fatal respiratory diseases [2]. To date, the determination of the potential molecular mechanisms of the pathogenicity of M. hyopneumoniae has been hampered by the fastidious growth requirements of this organism and the lack of tools for its genetic manipulation.

Comparative genomic analyses have been performed to reveal the genetic basis of virulence attenuation of M. hyopneumoniae and to predict potential virulence factors [35]. However, in the case of M. hyopneumoniae, the different strains share highly similar genome structures and gene orders, and few differences have been detected among them at the genomic level. Thus, we decided to turn our attention to the variations in expression levels of the predicted virulence factors.

Previous studies have described the transcriptome changes under different growth conditions [69] and with changes in gene regulation among different mycoplasma species [10]. No transcriptomic studies comparing virulent and avirulent M. hyopneumoniae strains have been published. The comparative proteomic reports of pathogenic and nonpathogenic strains based on two-dimensional gel electrophoresis have revealed a few differentially expressed proteins [1113]. However, gel-based proteomic methods are usually hindered by their low-throughput and the difficulty of direct quantitative comparison between samples. Recently, Paes et al. adopted cell fractioning technology coupled with mass spectrometry, providing much more potential virulence factor/vaccine candidates than conventional ones through comparative whole cell proteome profiles analysis of two M. hyopneumoniae strains and M. flocculare [14].

In this study, we undertook a large-scale proteomic comparison between the virulent M. hyopneumoniae strain 168 and its highly passaged attenuated vaccine strain 168L and identified differentially displayed proteins using the nongel-based isobaric tags for relative and absolute quantification (iTRAQ) approach. This survey supplied a comprehensive iTRAQ-based comparative proteomic analysis of a pathogenic and a nonpathogenic M. hyopneumoniae strain, both originating from one parent strain with a well-characterized genetic background. This approach eliminated the effects caused by differences at the genome level, and it demonstrated some novel promising virulence-related proteins in these strains compared with previous researches.

2. Materials and Methods

2.1. Ethics Statement

The animal experiments in this study were approved by The Scientific Ethic Committee of Huazhong Agricultural University, Wuhan, China (Approval Number: HZAURAB-2017-008) and conducted in accordance with the Hubei Regulations for the Administration of Affairs Concerning Experimental Animals.

2.2. Mycoplasma Strains and Their Cultivation

M. hyopneumoniae strains 168 and 168L were acquired from the Jiangsu Academy of Agricultural Sciences (Nanjing, China). M. hyopneumoniae strain 168 is a pathogenic strain isolated in China in 1974, leading typical mycoplasmal pneumonia of swine [5, 15]. The stable attenuated strain 168L was obtained after 380 continuous serial passages in KM2 cell-free medium (a modified Friis medium) and has been developed into a commercially available vaccine against M. hyopneumoniae in China [5, 16]. Cultures were maintained in KM2 cell-free medium at 37°C. For the proteomic study, 1.5 L cultures of strains 168 and 168L were grown to late-log phase as described by Calus [17]. Three replicate samples of each culture were harvested by centrifugation at 1,2000×g for 30 minutes at 4°C and washed three times with ice-cold PBS.

2.3. Protein Extraction, Digestion, and Labeling with iTRAQ Reagents

Bacterial cells were resuspended in lysis buffer (7 M urea, 2 M thiourea, 4% CHAPS, 40 mM Tris-HCl, pH 8.5) containing 1 mM PMSF, 2 mM EDTA, and 10 mM DTT. The samples were then sonicated for further protein solubilization. After centrifugation (25,000×g, 4°C, 20 min), the protein concentration of the supernatants was determined using the Bradford assay. Each sample (100 μg) was digested with trypsin (Promega, Madison, WI, USA) at 37°C for 12 h. The three lysate replicates of M. hyopneumoniae strain 168 were labeled with iTRAQ reagents 113, 114, and 115, and of strain 168L with iTRAQ reagents 116, 117, and 119.

2.4. SCX Chromatography

iTRAQ-labeled mixed peptides were fractionated by SCX chromatography using a 20AB HPLC system (Shimadzu, Kyoto, Japan). Peptide mixture was reconstituted in buffer A (25 mM NaH2PO4 in 25% acetonitrile, pH2.7) and loaded onto an UltremexSCX column (4.6 × 250 mm). The loaded peptides were eluted at a flow rate of 1 ml/min with 5% buffer B (25 mM NaH2PO4, 1 M KCl in 25% acetonitrile, pH 2.7) for 7 min, a linear gradient of 5–60% and 60-100% buffer B for 20 min and 2 min, respectively, and 100% buffer B for 1 min and then 5% buffer B for 10 min. The chromatograms were recorded at 214 nm. In total, 20 fractions were collected, desalted with a StrataX desalting column, and dried using vacuum centrifugation.

2.5. LC-MS/MS Analysis

The lyophilized fractions were resuspended in buffer A (5% acetonitrile, 0.1% formic acid) and loaded on a Symmetry C18 trapping column (5 μm particles, 180 μm inner diameter × 20 mm; Waters) for online trapping and desalting. Then the peptides were eluted from the trapping column over a BEH 130 C18 column (1.7 μm particles, 100 μm inner diameter × 100 mm) with a gradient buffer B (95% acetonitrile, 0.1% formic acid) at 300 nL/min (5-35% B, 35 min; 35-60% B, 5 min; 80% B, 2 min; 80-5% B, 1 min; 5% B, 10 min).

MS/MS was operated with a Triple TOF 5600 mass spectrometer (AB Sciex, Concord, ON). Data were acquired using a nanospray voltage of 2.5 kV, an interface heater temperature of 150°C, curtain gas at 30 psi, and nebulizer gas at 15 psi. Survey scans were acquired in 250 ms and as many as 30 product ion scans were collected if exceeding the threshold of 120 counts per second (cps) with +2 to +5 charge states. Total cycle time was fixed at 3.3 s. Dynamic exclusion was set for 1/2 of peak width (18 s) and the fragmentation energy was set to 35 ± 5 eV.

2.6. Protein Identification and Quantification

Tandem mass spectra were extracted by AB SCIEX MS Data Converter version 1.3. All MS/MS data were analyzed using the Mascot search algorithm (Matrix Science, London, UK; version b2.5.1) against the NCBI M. hyopneumoniae 168 & 168L database (201607, 3408 entries) with trypsin as the digestion enzyme. For protein identification, Mascot was searched with a parent ion tolerance of 15 ppm and a fragment ion mass tolerance of 0.050 Da. Carbamidomethyl of cysteine and iTRAQ8plex of lysine and the n-terminus were specified in Mascot as fixed modifications, while oxidation of methionine and iTRAQ8plex of tyrosine were allowed as variable modifications.

MS/MS-based peptide and protein identifications were validated using Scaffold (version Scaffold_4.6.1, Proteome Software Inc., Portland, OR). Peptide identifications were accepted if they could be established with less than 1% false discovery rate (FDR) by the Scaffold Local FDR algorithm, and each identified protein involved at least one unique peptide.

Scaffold Q+ (version Scaffold_4.6.1, Proteome Software Inc., Portland, OR) was used to quantitate the isobaric tag peptide and protein identifications. Normalization was performed iteratively (across samples and spectra) on intensities. Medians were used for averaging. Spectra data were log-transformed, pruned of those matched to multiple proteins, and weighted by an adaptive intensity weighting algorithm. A protein was considered to be differentially displayed if it contained at least two unique peptides with a minimum fold change of 0.585 on the scale, and Mann–Whitney Test with an unadjusted significance level p<0.05 corrected by the Benjamini-Hochberg method.

2.7. Bioinformatic Analysis of Proteins

Gene Ontology (GO) analysis was performed to map all of the identified proteins to GO terms in database ( The metabolic pathway of identified proteins was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database ( The significantly enriched GO terms or pathways of differentially expressed proteins in the background of the identified proteins were determined by a hypergeometric test (p<0.05).

The physical and chemical parameters of the given proteins were computed via ProtParam tool ( [18]. Lipoprotein signal peptides were predicted by LIPO ( [19] and LipoP 1.0 Server ( [20]. Subcellular location of proteins was predicted by using PSORTb version 3.0.2 ( [21] and CELLO v.2.5 ( [22]. TMHMM v2.0 ( and SignalP v.4.1 ( were used for predicting transmembrane helices and the presence and location of signal cleavage sites, respectively [23, 24]. SecretomeP 2.0 server ( was used to identify nonclassical protein secretion, i.e., not signal peptide triggered protein secretion [25]. Uncharacterized/hypothetical proteins were subjected to BlastP (, Pfam 31.0 web servers ( [26] and InterPro ( [27] to examine conserved functional domains for putative function assignment. For identifying proteins associated with virulence we used VirulentPred web server ( [28].

2.8. Gene Cloning, Expression, and Purification of Recombinant Proteins

Full-length genes of the four selected proteins were amplified directly from M. hyopneumoniae strain 168 genomic DNA using primer pairs in Table 1. Mycoplasmas use UGA as tryptophan codon while E. coli retain it as a stop codon. In order to mutate TGA to TGG and achieve mycoplasma proteins’ heterologous expression, we used specific primers (Table 1) to conduct overlapping PCR for site-directed mutagenesis. Then the products were cloned into pET-28a vector and transformed into E. coli DH5α. After checking the inserts by sequencing, the reconstructed plasmids were transformed into E. coli BL21 (DE3) to express the N-terminal 6×His-tagged recombinant proteins. Obtained proteins were purified by nickel affinity chromatograph (GE Healthcare, Piscataway, NJ, USA), dialyzed in PBS, and concentrated with Amicon Ultra centrifugal filter units (Millipore, Darmstadt, Germany). Protein concentration was determined by BCA Protein Assay Kit (Beyotime, Nanjing, China), and the purity was verified on 12% Coomassie-stained SDS-PAGE gels.

TargetPrimersSequence (5′ to 3′)Description

MHP168_174-M1RCGCATTTCCcCAAATTCCGCSite directed mutagenesis
MHP168_174-M1FGCGGAATTTGgGGAAATGCGSite directed mutagenesis

MHP168_271-M1RCTTTTTGACCATAcCAATTATCTGCSite directed mutagenesis
MHP168_271-M1FGCAGATAATTGgTATGGTCAAAAAGSite directed mutagenesis
MHP168_271-M2RGATAAAACCCTGcCAATCAGATTCSite directed mutagenesis
MHP168_271-M2FGAATCTGATTGgCAGGGTTTTATCSite directed mutagenesis

Translation elongation factor TuMHP168_533-FCTCGGATCCATGGCAGTTGTTAAAACGAmplification

Methylmalonate-semialdehyde dehydrogenaseMHP168_244-FCGGGGATCCATGAATTTTAAAAGAAATTACAmplification
MHP168_244-M1RCATTTGGAATTGAcCAATTTGGCACSite directed mutagenesis
MHP168_244-M1FGTGCCAAATTGgTCAATTCCAAATGSite directed mutagenesis
MHP168_244-M2RCGGCTTGTTTcCAAAGATCAGSite directed mutagenesis
MHP168_244-M2FCTGATCTTTGgAAACAAGCCGSite directed mutagenesis
MHP168_244-M3RGAATCTCTGAcCAGTCAATAATTGSite directed mutagenesis
MHP168_244-M3FCAATTATTGACTGgTCAGAGATTCSite directed mutagenesis
MHP168_244-M4RAATCAAATTAGTcCAGAAATTAATACTGGSite directed mutagenesis
MHP168_244-M4FCCAGTATTAATTTCTGgACTAATTTGATTSite directed mutagenesis
MHP168_244-M5RCTTCCGGGTTcCATTTTGTAGSite directed mutagenesis
MHP168_244-M5FCTACAAAATGgAACCCGGAAGSite directed mutagenesis
MHP168_244-RCTTCTCGAGTTATGACATcCAATCCTGAmplification and site directed mutagenesis

2.9. Production of Polyclonal Antibodies

Polyclonal antisera were produced in New Zealand White rabbits via four times subcutaneous immunization at 2-week intervals. Rabbits were first immunized with 1mg of purified proteins mixed with equal volumes of complete Freund’s adjuvant (Sigma, St. Louis, MI, USA). The rest three booster immunizations used 1mg of purified proteins suspended in equal volumes of incomplete Freund’s adjuvant (Sigma, St. Louis, MI, USA). Ten days after the last immunization, blood samples were collected via cardiac bleeding and centrifuged to obtain serum.

2.10. Western Blot Analysis

Western workflow utilizing stain-free technology was performed according to the supplier protocol. Samples of 10 ug were separated in 12% TGX Stain-Free Fastcast acrylamide gels (Bio-Rad, Hercules, CA, USA). After the electrophoresis, gels were placed in the ChemiDoc™ Touch Imaging System (Bio-Rad, Hercules, CA, USA) and activated for 5 min with UV treatment to visualize total proteins. Then the proteins were transferred to PVDF membranes. After blocking with 10% milk for 3 h at room temperature, membranes were incubated at room temperature for 4 h with rabbit polyclonal antibodies (1:500) against methylmalonate-semialdehyde dehydrogenase, endo-1,4-beta-glucanase, enolase, and translation elongation factor Tu. Goat anti-rabbit IgGs conjugated with horseradish peroxidase (1:4000; Beyotime, Nanjing, China) were used as secondary antibodies, and, after antibodies incubation, stain-free blot images were captured for total protein loading control and normalization. Then proteins were detected with ClarityTM ECL reagents (Bio-Rad, Hercules, CA, USA), and the intensity values of the target bands were normalized to that of the total proteins in the lane. The 168/168L intensity ratios were calculated for each protein, and three replicates were subjected to statistical analysis.

3. Results

3.1. Global Protein Profiles Overview

To establish more detailed comparable proteome profiles of the virulent strain 168 and its attenuated strain 168L, whole proteins were extracted from the two strains, digested and labeled with iTRAQ reagents, followed by downstream analysis LC-MS/MS.

Results of the iTRAQ assay showed that 489 proteins, 70% (489/695) of the coding sequences in M. hyopneumoniae, were identified from 30,214 spectra from the whole output of 322,266 spectra with an FDR of 1% (Table S1) [15]. The molecular weights of most of the identified proteins ranged between 10 and 60 kDa (68%, Figure 1(a)). 50.5% of the identified proteins possessed sequence coverage of the identified peptides of higher than 30% (Figure 1(b)).

Of the 489 proteins identified, 273 proteins were annotated into 9 groups based on the KEGG categories. The categories with the highest abundance were translation (70 proteins), carbohydrate metabolism (45 proteins), replication and repair (41 proteins), nucleotide metabolism (40 proteins), and membrane transport (40 proteins) (Figure 2(a)). Meanwhile, further classification indicated that the top ten pathways in which the identified proteins were involved were ribosome biosynthesis, microbial metabolism in diverse environments, biosynthesis of secondary metabolites, biosynthesis of antibiotics, carbon metabolism, ABC transporters, aminoacyl-tRNA biosynthesis, pyrimidine metabolism, purine metabolism, and biosynthesis of amino acids (Figure 2(b)).

3.2. Differentially Expressed Proteins

We detected and quantified 489 proteins, of which 35 were more abundant in the Mhp168 strain and 35 in the Mhp168L strain, based on having a ratio ≥ 0.585 or ≤-0.585 and a statistically significant differences at p<0.05 (Table 2).

Identified ProteinsAccession NumberGene_locusGene ratios (168/168L)Unique PeptidePercentage Coverage

proteins which were more abundant in 168L

Myo-inositol catabolism proteinADQ90457MHP168_246--3.041045%

Myo-inositol 2-dehydrogenaseADQ90458MHP168_247--2.891643%

Periplasmic sugar-binding proteinADQ90463MHP168_252rbsB-2.762161%

Methylmalonate-semialdehyde dehydrogenaseADQ90455MHP168_244--2.532148%

Serine hydroxymethyltransferaseADQ90454MHP168_243glyA-2.41436%


Myo-inositol 2-dehydrogenaseADQ90464MHP168_253--2.071139%

Myo-inositol catabolismADQ90459MHP168_248iolD-2.041727%

Myo-inositol catabolism proteinADQ90460MHP168_249iolE-2.04520%

Site-specific DNA-methyltransferaseADQ90525MHP168_318mod-1.65727%

Hypothetical proteinADQ90289MHP168_066--1.5455.80%

Ribose ABC transport ATP-binding proteinADQ90461MHP168_250mglA-1.391435%


Putative uncharacterized proteinADQ90727MHP168_531--1.216333%

Protein P97-copy 2ADQ90409MHP168_195--1.123438%

Glucose kinaseADQ90722MHP168_526glcK-1.04725%

Protein p97, cilium adhesinADQ90328MHP168_110--1.037064%

tRNA uridine 5-carboxymethylaminomethyl modification enzyme mnmGADQ90231MHP168_003gidA-1.011220%

Neutrophil activating factorADQ90913MHP168_736--0.98641%

Xylose ABC transporter ATP-binding proteinADQ90719MHP168_523xylG-0.931535%

Xylose ABC transporter permease proteinADQ90720MHP168_524xylH-0.9348.20%

ABC transporter ATP-binding proteinADQ90246MHP168_019--0.9168.50%


GTPase obgADQ90265MHP168_042obgE-0.88617%

Protein P102-copy 2ADQ90410MHP168_196--0.832125%

50S ribosomal protein L24ADQ90339MHP168_121rplX-0.81648%

Heat-inducible transcription repressor hrcAADQ90238MHP168_011hrcA-0.781336%

50S ribosomal protein L23ADQ90331MHP168_113rplW-0.771156%


Putative uncharacterized proteinADQ90524MHP168_317--0.761316%

DNA topoisomerase IADQ90512MHP168_305topA-0.711628%

50S ribosomal protein L14ADQ90338MHP168_120rplN-0.711064%

Holliday junction ATP-dependent DNA helicase ruvAADQ90632MHP168_430ruvA-0.67836%

Transcription termination-antitermination factor nusAADQ90786MHP168_593nusA-0.652134%

50S ribosomal protein L3ADQ90329MHP168_111rplC-0.61458%

proteins which were more abundant in 168

Triacylglycerol lipaseADQ90483MHP168_274lip0.591346%

Putative uncharacterized proteinADQ90691MHP168_495-0.59830%


Translation elongation factor TuADQ90729MHP168_533tuf0.612366%

Hypoxanthine phosphoribosyltransferaseADQ90492MHP168_283hpt0.61950%

Transport protein sgaTADQ90622MHP168_420ulaA0.62713%

Cation-transporting P-type ATPaseADQ90498MHP168_289pacL0.63813%

Hypothetical proteinADQ90530MHP168_323-0.63514%

Putative uncharacterized proteinADQ90824MHP168_634-0.662228%

Putative uncharacterized proteinADQ90745MHP168_549-0.681146%

Acetate kinase AckAADQ90711MHP168_515ackA0.691754%

Cysteinyl-tRNA synthetaseADQ90838MHP168_650cysS0.69715%

Putative uncharacterized proteinADQ90859MHP168_673-0.71128%

Pts system, lichenan-specific IIa componentADQ90264MHP168_041licA0.71933%

Probable O-sialoglycoprotein endopeptidaseADQ90833MHP168_644gcp0.731037%

CTP synthaseADQ90509MHP168_302pyrG0.741425%

Lipoate-protein ligase AADQ90531MHP168_324lplA-10.751433%

ABC transporter ATP-binding-Pr1ADQ90875MHP168_691pr10.75716%

Hypothetical proteinADQ90910MHP168_731-0.76435%

Thymidine kinaseADQ90808MHP168_618tdk0.78840%

Lipoate-protein ligase AADQ90484MHP168_275lplA0.78516%

Purine-nucleoside phosphorylaseADQ90309MHP168_087deoD0.8729%


NADH dependent flavin oxidoreductaseADQ90532MHP168_325baiH0.85411%

ABC transporter proteinADQ90806MHP168_616-0.87611%

Sugar ABC transporter ATP-binding proteinADQ90805MHP168_615-0.91026%

MG2+ ion transporterADQ90693MHP168_497mgtE0.9923%

Ribonuclease IIIADQ90602MHP168_400rnc0.94934%

Ribose-phosphate pyrophosphokinaseADQ90852MHP168_664prsA0.981239%

Glycerol transporter subunit AADQ90619MHP168_417gtsA1.13411%

Putative uncharacterized proteinADQ90827MHP168_638-1.1647.50%

hypothetical proteinADQ90324MHP168_104-1.591923%

Glycerol kinaseADQ90587MHP168_385glpK1.592249%


PTS system galactitol-specific enzyme IIB componentADQ90760MHP168_564-1.78457%

Gene Ontology (GO) enrichment analysis was then performed for the differentially expressed proteins for the categories “molecular function” (F), “cellular components” (C), and “biological process” (P) (Table S2). Figure 3(a) shows the result of the enrichment analysis for these functions (GO terms). A total of 12 GO terms were enriched for molecular function and 30 for biological processes, while there’s no statistically significant results for cellular components. For the GO category “molecular function”, significant synthesis of the proteins related to catalytic activity were detected. In addition, the GO terms “N-methyltransferase activity”, “kinase activity”, and “magnesium ion binding” were also highly enriched. The highest enrichment for biological processes was associated with small molecule catabolism. In addition, the GO terms “inositol catabolic process”, “organic hydroxy compound catabolic process”, “cellular carbohydrate catabolic process”, “alcohol catabolic process”, “polyol catabolic process”, “DNA methylation or demethylation”, “DNA alkylation”, “DNA methylation”, and “inositol metabolic process” were also highly enriched.

To examine which biological pathways were apparently altered, we carried out a pathway analysis. Enriched pathways were grouped into 2 categories with P value <0.05 (Table S3, Figure 3(b)). The two pathways were involved in microbial metabolism in diverse environments and in inositol phosphate metabolism.

3.3. Expression of Recombinant Proteins and Preparation of Polyclonal Antibodies

Four differently expressed proteins, endo-1,4-beta-glucanase, enolase, translation elongation factor Tu, and methylmalonate-semialdehyde dehydrogenase were chosen to verify the proteomic differences. Within the virulent strain 168 samples, endo-1,4-beta-glucanase, enolase, and translation elongation factor Tu were more abundant, which might conduce to the virulence of M. hyopneumoniae. All these three proteins are cytosolic enzymes, lacking transmembrane domains or traditional signal sequences (Table 3), yet intensive researches during the past years have suggested that they could be detected on microbe surface and moonlight as adhesins [2934]. On one hand we expect to obtain sera to validate the comparative results and on the other hand we are calculated to conduct further experiments around these proteins. We expressed these recombinant proteins in prokaryotic system (Figure 4(a)) then got the rabbit sera against the proteins, respectively (Figure 4(b)).

ProteinNumber of amino acidsMolecular weightTheoretical PINumber of predicted TMHsSignal peptideNumber of TGA codons

Endo-1,4-beta-glucanase356 aa39.2kDa6.210No1

Enolase452 aa49.5kDa5.830No2

Translation elongation factor Tu402 aa44.1kDa5.610No0

Methylmalonate-semialdehyde dehydrogenase489 aa53.9kDa8.460No6

3.4. Validation of Selected Proteins by Western Blot

We performed western blot analyses using a stain-free technology to assess the levels of four selected differentially expressed proteins. As shown in Figure 5, the alterations in expression levels of these proteins between strain 168 and 168L were consistent with the results from quantitative proteomic analysis, suggesting the credibility of our proteomic data.

4. Discussion

In this study, we compared the protein profiles of the pathogenic M. hyopneumoniae strain 168 and its highly passaged attenuated strain 168L using iTRAQ strategy for the first time. We identified 70 differentially expressed proteins to mine candidate virulence determinants and proteins or biological processes leading to attenuation.

4.1. Proteins Involved in Inositol Phosphate Metabolism

Among all the sequenced mycoplasma species, M. hyopneumoniae is the only one having a gene cluster for myo-inositol utilization [35]. A previous report, based on the reconstruction of a genome-scaled metabolic model for three mycoplasmas in silico, suggested that the myo-inositol metabolism may be one of the reasons for the high virulence of M. hyopneumoniae compared to that of M. flocculare and M. hyorhinis [35]. In our study, seven related enzymes of inositol phosphate metabolism, namely, myo-inositol 2-dehydrogenase (MHP168_253), myo-inositol catabolism protein (MHP168_246), myo-inositol 2-dehydrogenase (MHP168_247), methylmalonate-semialdehyde dehydrogenase (MHP168_244), 5-dehydro-2-deoxygluconokinase (IolC), myo-inositol catabolism (IolD), and myo-inositol catabolism protein (IolE), were more abundant in the vaccine strain. These proteins covered all the enzymes comprising the classical myo-inositol bacterial catabolic pathway in M. hyopneumoniae, with the exception of the enzyme acting as 5-dehydro-2-deoxyphosphogluconate aldolase (IolJ) whose gene was also absent in its genome. Nevertheless, Ferrarini et al. [36] hypothesized that one gene copy among those annotated as fructose-biphosphate aldolase (Fba) from M. hyopneumoniae might function as IolJ from other organisms and confirmed that M. hyopneumoniae had been able to utilize myo-inositol from the culture medium.

The increments of those protein levels might facilitate the catabolism of myo-inositol to produce dihydroxyacetone-phosphate (DHAP) and acetyl-coenzyme A (CoA). DHAP can enter glycolysis while acetyl CoA can be widely used in macromolecular biosynthesis and energy production to support cell growth and proliferation. The presence of inositol in mammalian hosts’ bloodstream [37] and the degradation of phosphatidylinositol from host’ pulmonary surfactant [38] make myo-inositol available for the strains in vivo. Considering the fact that myo-inositol has been one stable readily abundant component in the culture medium containing swine serum, we may assume that M. hyopneumoniae retains the ability of degrading myo-inositol in vivo, then the capability is increased with passage in vitro. The vaccine strain might have developed an enhanced ability to cope with differences in environmental carbon sources. However, the myo-inositol utilization gene cluster existed in both strains, and there were no genetic variations in the coding sequences or intergenic regions of those genes [5]. It may be suggested that strains which possess this gene cluster could express it constitutively, while the regulatory mechanism for the enhancements is unclear and still need future research.

4.2. Proteins Involved in Nucleotide Metabolism

M. hyopneumoniae cannot synthesize purines and pyrimidines de novo. However, this organism is able to take up exogenous nucleobases and nucleosides, as well as those produced internally by DNA and RNA degradation, and then synthesize nucleotides through salvage pathways and interconversions [39]. In the virulent strain, the upregulation of transport protein SgaT (UlaA) and PTS system galactitol-specific enzyme IIB component (MHP168_564) participated in ascorbate transport suggests the enhanced influx of D-Xylulose-5-phosphate (Xyl5P) through ascorbate metabolism into the pentose phosphate pathway to generate phosphoribosyl pyrophosphate (PRPP) (Figure 6). Meanwhile, the increased expression of ribose-phosphate pyrophosphokinase (PrsA) of the pentose phosphate pathway, catalyzing ribose-5-phosphate (R5P) to PRPP, implies the enriched production of PRPP which is a key intermediate to synthesize purine and pyrimidine nucleotides, as well as nicotinamide adenine dinucleotide (NAD) (Figure 6) [40]. We found that in our data four enzymes participated in nucleoside metabolism: purine-nucleoside phosphorylase (DeoD), hypoxanthine-guanine phosphoribosyl transferase (Hpt), thymidine kinase (Tdk), and CTP synthase (PyrG), all upregulated in the virulent strain. NAD is required for the generation of NADPH that responds to oxidative burst, while synthesis of nucleotides can also promote DNA damage repair caused by oxidative stress [41]. Thus, the increased synthesis of PRPP, NAD, and the upregulation of the nucleotide synthesis in virulent strain 168 could be identified as defense to reactive oxygen species (ROS) damage, fitting well with the idea of an elevated virulence.

4.3. Potential Virulence Factors

Mycoplasmas contain a range of virulence factors in their pathogenic machinery, including enzymes, transporters, transcriptional regulators, lipoproteins [42]. Among the overrepresented proteins in virulent strain 168, proteins associated with virulence usually play important roles in M. hyopneumoniae offence and defense, meriting more attention. To predict putative novel virulence factors and search for potential drug/vaccine targets, we used VirulentPred to perform the in silico analysis for the 35 upexpressed proteins in the virulent strain 168. The server categorized 18 proteins to be associated with virulence. Out of these, 10 proteins are ones with known functions, and the others are termed as “hypothetical” or “uncharacterized”. As a complement to the bioinformatics analyses, we also searched manually the 35 proteins and/or their paralogues for the references that support their virulence-related activities in M. hyopneumoniae, or in other pathogens. Then additional 8 proteins classified as “nonvirulent” were recruited into the group of putative virulence proteins. The results are summarized in Table 4 and the listed proteins are discussed in the main text below.

Identified ProteinsAccession NumberPrediction resultsReference

Putative uncharacterized proteinADQ90691VirulentOkan et al. 2013; Shahbaaz et al. 2015

Hypothetical proteinADQ90530VirulentCheng et al. 2017

Putative uncharacterized proteinADQ90824VirulentHeithoff et al. 1999; Shahbaaz et al. 2015

Putative uncharacterized proteinADQ90745VirulentGay et al. 1991; Kobe et al. 2001

Cysteinyl-tRNA synthetaseADQ90838Virulent

Putative uncharacterized proteinADQ90859VirulentCheng et al. 1996; Monnerat et al. 1999

Pts system, lichenan-specific IIa componentADQ90264Virulent

ABC transporter ATP-binding-Pr1ADQ90875VirulentGarmory et al. 2004

Hypothetical proteinADQ90910VirulentTwine et al. 2006

ABC transporter proteinADQ90806VirulentGarmory et al. 2004

Sugar ABC transporter ATP-binding proteinADQ90805VirulentGarmory et al. 2004

MG2+ ion transporterADQ90693VirulentGroisman et al. 2013; Merino et al. 2001; O'Connor et al. 2009

Ribonuclease IIIADQ90602VirulentHaddad et al. 2013; Mao et al. 2016

Ribose-phosphate pyrophosphokinaseADQ90852VirulentDonini et al. 2017

Putative uncharacterized proteinADQ90827Virulent

Hypothetical proteinADQ90324Virulent


PTS system galactitol-specific enzyme IIB componentADQ90760Virulent

EnolaseADQ90480NonvirulentBao et al. 2014; Pancholi 2001; Schreiner et al. 2012

Translation elongation factor TuADQ90729NonvirulentDallo et al. 2002; Jiang et al. 2016

Cation-transporting P-type ATPaseADQ90498NonvirulentNovoa-Aponte et al. 2014

Probable O-sialoglycoprotein endopeptidaseADQ90833NonvirulentAruni et al. 2011; Davis et al. 2003

Lipoate-protein ligase AADQ90531NonvirulentO'Riordan et al. 2003

Lipoate-protein ligase AADQ90484NonvirulentO'Riordan et al. 2003

Endo-1,4-beta-glucanaseADQ90389NonvirulentRobinson et al. 2013

Acetate kinase AckAADQ90711Nonvirulent

NADH dependent flavin oxidoreductaseADQ90532Nonvirulent

Glycerol transporter subunit AADQ90619Nonvirulent

Glycerol kinaseADQ90587Nonvirulent

Triacylglycerol lipaseADQ90483Nonvirulent

Hypoxanthine phosphoribosyltransferaseADQ90492Nonvirulent

Transport protein sgaTADQ90622Nonvirulent

CTP synthaseADQ90509Nonvirulent

Thymidine kinaseADQ90808Nonvirulent

Purine-nucleoside phosphorylaseADQ90309Nonvirulent

Lipoproteins. Three lipoproteins were identified in this study. MHP168_392 and MHP168_393 were downexpressed and MHP168_418 were upexpressed in the virulent strain 168. Most lipoproteins are surface exposed and important components of mycoplasma membranes, and some in M. hyopneumoniae have been identified to play important roles in pathogenesis [4349]. Mhp378, a homolog of MHP168_392 from M. hyopneumoniae 232, has been identified as a species-specific, highly immunogenic membrane-associated protein [50]. The homolog of MHP168_393 in M. hyopneumoniae 232 (mhp379) is a surface-exposed exonuclease with probable function in importing nucleic acid precursors [47]. Furthermore, surface located MHP168_418 is shown with the ability to induce apoptosis of porcine peripheral blood mononuclear cells in vitro [51]. M. hyopneumoniae usually cannot employ the similar mechanisms used by other mycoplasma species to generate antigenic diversity through genetic variation [52]. Nonetheless, the posttranslational proteolytic processing, which targeted mycoplasmal membrane lipoproteins and other surface-associated proteins, could create a dynamic surface topograph [53]. The identified lipoproteins with differential expression levels, as well as their probable posttranslational cleavage, could lead to antigenic variations, thereby affecting its immune evasion in host.

Adhesion-Related Proteins. The adherence of M. hyopneumoniae to ciliated respiratory epithelium is mainly mediated by the membrane protein P97 [54]. The finding of higher expression levels of P97, protein P97-copy 2, and protein P102-copy 2 in the vaccine strain compared with those of the virulent strain was unexpected, but similar results were also found in the transcriptome comparison between M. flocculare and M. hyopneumoniae [10]. A previous study indicated that the cilium binding domain of P97 is found exclusively in the R1 region, the functional site requiring a minimum of eight tandem repeating units (AAKPV/E) [54]. However, three transversion mutations (E863V) occurred in the tandem repeating units in P97 of 168L, which might partly affect the adhesion of vaccine strain [5]. P102 accompanies P97 to form a two-gene operon. Both the P97 and P102 genes have several paralogs within the M. hyopneumoniae genome [5]. The sequences of the paralogs are uncompleted. Furthermore, the p97-copy 2 protein of both strains lacks the R1 domain while protein P102-copy 2 was also truncated [5]. Thus, a higher abundance of these proteins might not mean a stronger adhesion ability of the vaccine strain. Adherence to ciliated respiratory epithelium is a multifactorial process that also involves other proteins.

In our proteomic data, we detected several upregulated moonlighting proteins in strain 168 that could be used to invade host cells: elongation factor Tu [30, 31], enolase [29, 32, 34], and endo-1,4-beta-glucanase [33]. Intensive reports have described some bacterial metabolic enzymes not only performing key metabolic functions in the cytosol of bacterial cell but also locating on the bacterial surface without a signal sequence, moonlighting as an adhesion contributor to host cells [55]. These multifunctional proteins located on the mycoplasma surface could adhere to the swine tracheal cilia and bind to host factors plasminogen and fibronectin. Mycoplasma surface-bound plasminogen is converted to plasmin by tissue plasminogen activator. Plasmin cleaves host extracellular matrix proteins and activates matrix metalloproteases, assisting the pathogen to invade host issues and providing amino acids for growth of M. hyopneumoniae [33]. Fibronectin is widespread in the ciliary borders of the bronchioles and could bind to glycosaminoglycans, collagens, DNA, fibrin, and cell surface integrins, etc. These properties of fibronectin make it a physical bridge between pathogens and host cells [56]. Those mycoplasma surface proteins binding to plasminogen and fibronectin have been suggested to be associated with virulence and warrant further investigation.

Significantly expressed protein O-sialoglycoprotein endopeptidase in strain 168 also regulates the invasion of pathogen via degrading sialoglycosylated host-cell proteins and mediating adhesion to host cells [57, 58]. The traits of these proteins, playing roles in the pathogen-host interactions, affect the pathogenesis of M. hyopneumoniae in part.

Transporters. Transporter proteins play important roles in bacteria, transporting various molecules to support survival and growth in different niches [59]. Cation-transporting P-type ATPase (PacL) is a member of transmembrane P-type ATPases, which are involved in transportation of ions and phospholipids, using the energy derived from ATP hydrolysis [6062]. Mg2+ ion transporter (MgtE) is a highly Mg2+-selective channel gated by Mg2+, transporting substrates across the cytoplasmic membrane by utilizing the electrochemical gradient [63]. Both of MgtE and PacL function in maintaining metal homeostasis in pathogen and survival in host, thus considered as virulence determinants in some bacteria [6265]. Inactivation of ABC transporters often has deleterious effects on the virulence in bacteria, resulting in attenuated phenotypes and decreased adherence to host cells [59]. ABC transporter ATP-binding-Pr1 (Pr1), ABC transporter protein (MHP168_616), and sugar ABC transporter ATP-binding protein (MHP168_615) may be associated with virulence via participating in the regulation of cation homeostasis or adhesive ability. Pts system lichenan-specific IIa component (LicA) and PTS system galactitol-specific enzyme IIB component (MHP168_564) are predicted to be virulent factors on VirulentPred webserver, and their expression levels in strain 168 were upregulated as compared to 168L, suggesting that these proteins may be suitable targets for antibacterial vaccine and therapies.

Enzymes. Two lipoate-protein ligase A (LplA-1 and LplA) were upregulated in the pathogenic strain 168. LplA ligates exogenous lipoic acid to lipoyl domains of certain metabolic enzymes complexes involved in oxidative metabolism [66]. Previous study showed that growth of LplA1-deficient L. monocytogenes was damaged specifically in the host and virulence was 1/300th as that of wild-type in animals [67]. Possibly, we speculate that LplA could also act as a virulence factor in M. hyopneumoniae.

Enzymes PrsA have been treated as potential targets for therapeutic and vaccine candidates according to previous reports [68]. Cysteinyl-tRNA synthetase is essential for bacteria growth and classified as a “virulent” protein in our predicted results. Increased level of these proteins in virulent strain 168 could support this idea; nonetheless, what kind of part do these proteins play in the pathogenesis and whether they could indeed cause the predicted effects in this strain still needs experimental verification.

Transcriptional Regulator. In strain 168, we found the expression of ribonuclease III (Rnc) was significantly elevated. This ribonuclease is known to be assisting in pathogenesis of other bacteria via regulating the synthesis of virulence factors through RNase III-dependent posttranscriptional manner [69, 70]; its upexpression in virulent strain made us guess that it might also modulate gene expression of virulent factors in M. hyopneumoniae, therefore contributing to high virulence.

Uncharacterized/Hypothetical Virulence-Related Proteins. According to the previously published genome data, 268 of 695 (~39%) M. hyopneumoniae proteins were not assigned functions and annotated as “uncharacterized” or “hypothetical” proteins [15]. In many sequenced bacterial genomes, the uncharacterized/hypothetical proteins account for around 20-40% of the total genome [42, 71], which are important for complementing the genomic and proteomic framework theory. Understanding the functional properties of these proteins will be crucial for a more profound understanding of the microbes at the molecular machinery level. For pathogens, those proteins with unknown functions probably involve in virulence which help elucidating pathogenesis and are virgin area for developing novel drug/vaccine targets.

In our study, 45 hypothetical proteins and 78 uncharacterized proteins were identified. Among them, 11 uncharacterized/hypothetical proteins displayed differential expression levels, 8 were overrepresented, and 3 were underrepresented in strain 168. It is reasonable for us to suppose that the differentially expressed proteins of unknown function might give a clue for pathogenic mechanisms of M. hyopneumoniae or novel drug/vaccine targets. The VirulentPred results showed that all the overrepresented proteins in strain 168 among the differentially expressed uncharacterized/hypothetical proteins are predicted to be associated with virulence. Thus the protein sequences were submitted to BlastP, Pfam, and InterPro web servers for putative function annotation. Table 5 shows the results of the BlastP, Pfam, and InterPro. We also performed the prediction of subcellular localization as a complement to facilitate our knowledge of these uncharacterized/hypothetical proteins (Table S4). However, experimental studies should be carried out to assess and confirm their functions in biological processes and pathogenesis.

Accession NumberBlastPpfamInterPro

proteins which were more abundant in 168

ADQ90691Cof-type HAD-IIB family hydrolaseHydrolase_3 DomainNone predicted.

ADQ90530deacetylase SIR2No resultSirtuin family, catalytic core domain

ADQ90824SAM-dependent DNA methyltransferaseN6_Mtase FamilyDNA methylase, adenine-specific

ADQ90745Leucine-rich repeat domain-containing proteinLRR_5 RepeatBspA type Leucine rich repeat region

ADQ90859LppA family lipoproteinNo resultNone predicted

ADQ90910Glycine cleavage system protein HGCV_H DomainNone predicted

ADQ90827DUF31 domain-containing proteinDUF31 Domain (Putative peptidase)Putative peptidase DUF31

ADQ90324No resultNo resultNone predicted

proteins which were more abundant in 168L

ADQ90289SGNH/GDSL hydrolase family proteinNo resultBromodomain

ADQ90727SGNH/GDSL hydrolase family proteinLipase_GDSL FamilyGDSL lipase/esterase

ADQ90524DEAD/DEAH box helicaseResIII FamilyHelicase/UvrB, N-terminal

Proteins ADQ90691, ADQ90727, and ADQ90524 were annotated as hydrolase enzymes. Members of this class are involved in various significant biological processes, including virulence mechanisms [42, 72]. ADQ90691 was significantly expressed in virulent strain and predicted as a virulence factor, probably playing an important role in the pathogenesis of M. hyopneumoniae. This protein was annotated to be a Cof-like hydrolase, while might function as the phosphatase Cof in E.coli, catalyzing the hydrolysis of 4-amino-2-methyl-5-hydroxymethylpyrimidine pyrophosphate to 4-amino-2-methyl-5-hydroxymethylpyrimidine phosphate [73].

The hypothetical protein ADQ90530 was categorized as a sirtuin (also known as Sir2), which is responsible for NAD+-dependent deacetylation. As the deacetylase of M. fermentans is expressed inside mammalian cells, it inhibits cell proliferation but promotes their antioxidation and antistarvation capacities, and alters gene expression, affecting physiological functions and the corresponding signal transduction pathways in host cells [74]. The long-held view is that M. hyopneumoniae is an extracellular pathogen, thus it remains unclear whether ADQ90530 in M. hyopneumoniae functions similarly with secreted deacetylase in intracellular pathogen M. fermentans. Elucidating this issue will require further investigation.

Protein ADQ90824 was predicted to be an N-6 adenine-specific DNA methylase, which controls methylation at adenine residues of important biological processes. The methylation process plays an important role in bacterial pathogenesis by regulating the synthesis of virulence factors, and DNA adenine methylases serve as promising antimicrobials and vaccines targets [42, 75].

Leucine-rich repeats (LRRs) are found to be present in a number of proteins with diverse functions, including cell-adhesion molecules, virulence factors, and extracellular matrix-binding glycoproteins, and function in signal transduction, cell-adhesion, and protein-protein interactions [76, 77]. Protein ADQ90745 containing a leucine-rich repeat domain is predicted to be related to virulence, yet more detailed information about its function is unavailable so far.

Uncharacterized protein ADQ90859 belongs to the mycoides cluster lipoprotein, LppA/p72 family; members of this protein family are predicted lipoproteins with a typical prokaryotic signal peptidase II processing and lipid attachment site [78]. Paralogues in other mycoplasmas have been identified as specific antigenic proteins with potential for use in development of diagnostic reagents [78, 79].

Hypothetical protein ADQ90910 was classified as one member (H-protein) of the glycine cleavage system composed of four proteins: the T-, P-, L-, and H-protein. This system catalyzes the reversible reaction: Glycine + H4folate + NAD+ <==> 5, 10-methylene-H4folate + CO2 + NH3 + NADH + H+, and H-protein shuttles some of the intermediate products [80]. H-protein was expressed exclusively in Francisella tularensis isolated from mouse spleens compared with in vitro grown controls, suggesting that H-protein may play an important role in the metabolic fitness of Francisella tularensis [81]. A similar result of predicted H-protein was detected in proteomic analysis of M. hyopneumoniae; we assumed hypothetical protein ADQ90910 also play a part in the adaptive response of pathogen to in vivo environment.

The “DUF31” domain that has no known function was found in uncharacterized protein ADQ90827, and none of the domains was predicted in hypothetical protein ADQ90324. The function of them cannot be identified but is predicted to be associated with virulence.

5. Conclusions

In conclusion, this survey identified mycoplasmal proteins and unveiled biological processes potentially related to the difference in virulence between the pathogenic and attenuated strains of M. hyopneumoniae. The components that play roles in these critical mechanisms are natural targets for specific drugs. These potential virulence factors are usually immunogenic and can be treated as vaccine candidates [82, 83]. Our future experimental work will focus on those proteins to delineate the pathogenic mechanism(s) and identify those which can be targeted for drug design and vaccine development.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


This work was supported by the Key Technology R&D Programme of China 2015BAD12B02 and the National Natural Science Foundation of China 31170160 and 31502069. We would like to thank Dr. Guoqing Shao at Jiangsu Academy of Agricultural Sciences for providing the strains.

Supplementary Materials

Supplementary 1 . Table S1: information of all proteins identified in this study.

Supplementary 2 . Table S2: Gene Ontology enrichment analysis for differentially expressed proteins.

Supplementary 3 . Table S3: pathway enrichment analysis of the annotated differentially expressed proteins.

Supplementary 4 . Table S4: bioinformatic analyses of the differentially expressed uncharacterized/hypothetical proteins.


  1. D. Maes, J. Segales, T. Meyns, M. Sibila, M. Pieters, and F. Haesebrouck, “Control of Mycoplasma hyopneumoniae infections in pigs,” Veterinary Microbiology, vol. 126, no. 4, pp. 297–309, 2008. View at: Publisher Site | Google Scholar
  2. S. Simionatto, S. B. Marchioro, D. Maes, and O. A. Dellagostin, “Mycoplasma hyopneumoniae: From disease to vaccine development,” Veterinary Microbiology, vol. 165, no. 3-4, pp. 234–242, 2013. View at: Publisher Site | Google Scholar
  3. F. M. Siqueira, C. E. Thompson, V. G. Virginio et al., “New insights on the biology of swine respiratory tract mycoplasmas from a comparative genome analysis,” BMC Genomics, vol. 14, no. 175, pp. 1471–2164, 2013. View at: Google Scholar
  4. A. T. Vasconcelos, H. B. Ferreira, C. V. Bizarro, and etal., “Swine and poultry pathogens: the complete genome sequences of two strains of mycoplasma hyopneumoniae and a strain of mycoplasma synoviae,” Journal of Bacteriology, vol. 187, no. 16, pp. 5568–5577, 2005. View at: Publisher Site | Google Scholar
  5. W. Liu, S. Xiao, M. Li et al., “Comparative genomic analyses of Mycoplasma hyopneumoniae pathogenic 168 strain and its high-passaged attenuated strain,” BMC Genomics, vol. 14, no. 80, pp. 1471–2164, 2013. View at: Publisher Site | Google Scholar
  6. M. L. Madsen, D. Nettleton, E. L. Thacker, R. Edwards, and F. C. Minion, “Transcriptional profiling of mycoplasma hyopneumoniae during heat shock using microarrays,” Infection and Immunity, vol. 74, no. 1, pp. 160–166, 2005. View at: Publisher Site | Google Scholar
  7. M. L. Madsen, D. Nettleton, E. L. Thacker, and F. C. Minion, “Transcriptional profiling of Mycoplasma hyopneumoniae during iron depletion using microarrays,” Microbiology, vol. 152, no. 4, pp. 937–944, 2006. View at: Publisher Site | Google Scholar
  8. M. L. Madsen, S. Puttamreddy, E. L. Thacker, M. D. Carruthers, and F. C. Minion, “Transcriptome changes in mycoplasma hyopneumoniae during infection,” Infection and Immunity, vol. 76, no. 2, pp. 658–663, 2008. View at: Publisher Site | Google Scholar
  9. M. J. Oneal, E. R. Schafer, M. L. Madsen, and F. C. Minion, “Global transcriptional analysis of Mycoplasma hyopneumoniae following exposure to norepinephrine,” Microbiology, vol. 154, no. 9, pp. 2581–2588, 2008. View at: Publisher Site | Google Scholar
  10. F. M. Siqueira, A. L. Gerber, R. L. M. Guedes et al., “Unravelling the transcriptome profile of the swine respiratory tract mycoplasmas,” PLoS ONE, vol. 9, no. 10, 2014. View at: Google Scholar
  11. Y. Li, Y. Ho, S. Chen, T. Chiou, Z. Li, and D. Shiuan, “Proteomic comparative analysis of pathogenic strain 232 and avirulent strain J of Mycoplasma hyopneumoniae,” Biochemistry (Moscow), vol. 74, no. 2, pp. 215–220, 2009. View at: Publisher Site | Google Scholar
  12. P. M. Pinto, C. S. Klein, A. Zaha, and H. B. Ferreira, “Comparative proteomic analysis of pathogenic and non-pathogenic strains from the swine pathogen Mycoplasma hyopneumoniae,” Proteome Science, vol. 7, no. 45, pp. 1477–5956, 2009. View at: Google Scholar
  13. Y. Yu, M. Liu, L. Hua et al., “Fructose-1,6-bisphosphate aldolase encoded by a core gene of Mycoplasma hyopneumoniae contributes to host cell adhesion,” Veterinary Research, vol. 49, no. 1, pp. 018–0610, 2018. View at: Publisher Site | Google Scholar
  14. J. A. Paes, L. D. Machado, F. M. dos Anjos Leal et al., “Comparative proteomics of two Mycoplasma hyopneumoniae strains and Mycoplasma flocculare identified potential porcine enzootic pneumonia determinants,” Virulence, vol. 9, no. 1, pp. 1230–1246, 2018. View at: Publisher Site | Google Scholar
  15. W. Liu, Z. Feng, L. Fang et al., “Complete genome sequence of mycoplasma hyopneumoniae strain 168,” Journal of Bacteriology, vol. 193, no. 4, pp. 1016-1017, 2011. View at: Publisher Site | Google Scholar
  16. Q. Xiong, Y. Wei, H. Xie et al., “Effect of different adjuvant formulations on the immunogenicity and protective effect of a live Mycoplasma hyopneumoniae vaccine after intramuscular inoculation,” Vaccine, vol. 32, no. 27, pp. 3445–3451, 2014. View at: Publisher Site | Google Scholar
  17. D. Calus, D. Maes, K. Vranckx, I. Villareal, F. Pasmans, and F. Haesebrouck, “Validation of ATP luminometry for rapid and accurate titration of Mycoplasma hyopneumoniae in Friis medium and a comparison with the color changing units assay,” Journal of Microbiological Methods, vol. 83, no. 3, pp. 335–340, 2010. View at: Publisher Site | Google Scholar
  18. M. R. Wilkins, E. Gasteiger, A. Bairoch et al., “Protein identification and analysis tools in the ExPASy server,” Methods in Molecular Biology, vol. 112, pp. 531–552, 1999. View at: Google Scholar
  19. F. S. Berven, O. A. Karlsen, A. H. Straume et al., “Analysing the outer membrane subproteome of Methylococcus capsulatus (Bath) using proteomics and novel biocomputing tools,” Archives of Microbiology, vol. 184, no. 6, pp. 362–377, 2006. View at: Publisher Site | Google Scholar
  20. A. S. Juncker, H. Willenbrock, G. von Heijne, S. Brunak, H. Nielsen, and A. Krogh, “Prediction of lipoprotein signal peptides in Gram-negative bacteria,” Protein Science, vol. 12, no. 8, pp. 1652–1662, 2003. View at: Publisher Site | Google Scholar
  21. N. Y. Yu, J. R. Wagner, M. R. Laird et al., “PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes,” Bioinformatics, vol. 26, no. 13, pp. 1608–1615, 2010. View at: Publisher Site | Google Scholar
  22. C. S. Yu, C. J. Lin, and J. K. Hwang, “Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions,” Protein Science, vol. 13, no. 5, pp. 1402–1406, 2004. View at: Publisher Site | Google Scholar
  23. A. Krogh, B. Larsson, G. Von Heijne, and E. L. L. Sonnhammer, “Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes,” Journal of Molecular Biology, vol. 305, no. 3, pp. 567–580, 2001. View at: Publisher Site | Google Scholar
  24. T. N. Petersen, S. Brunak, G. Von Heijne, and H. Nielsen, “SignalP 4.0: discriminating signal peptides from transmembrane regions,” Nature Methods, vol. 8, no. 10, pp. 785-786, 2011. View at: Publisher Site | Google Scholar
  25. J. D. Bendtsen, L. Kiemer, A. Fausbøll, and S. Brunak, “Non-classical protein secretion in bacteria,” BMC Microbiology, vol. 5, no. 58, pp. 1471–2180, 2005. View at: Publisher Site | Google Scholar
  26. S. El-Gebali, J. Mistry, A. Bateman et al., “The Pfam protein families database in 2019,” Nucleic Acids Research, vol. 47, no. D1, pp. D427–D432, 2019. View at: Publisher Site | Google Scholar
  27. A. L. Mitchell, T. K. Attwood, P. C. Babbitt et al., “InterPro in 2019: improving coverage, classification and access to protein sequence annotations,” Nucleic Acids Research, vol. 47, no. D1, pp. D351–D360, 2019. View at: Google Scholar
  28. S. Saha and G. P. S. Raghava, “VICMpred: an SVM-based method for the prediction of functional proteins of gram-negative bacteria using amino acid patterns and composition,” Genomics, Proteomics & Bioinformatics, vol. 4, no. 1, pp. 42–47, 2006. View at: Publisher Site | Google Scholar
  29. S. Bao, X. Guo, S. Yu et al., “Mycoplasma synoviaeenolase is a plasminogen/fibronectin binding protein,” BMC Veterinary Research, vol. 10, no. 223, pp. 014–0223, 2014. View at: Publisher Site | Google Scholar
  30. S. F. Dallo, T. R. Kannan, M. W. Blaylock, and J. B. Baseman, “Elongation factor Tu and E1 subunit of pyruvate dehydrogenase complex act as fibronectin binding proteins in Mycoplasma pneumoniae,” Molecular Microbiology, vol. 46, no. 4, pp. 1041–1051, 2002. View at: Publisher Site | Google Scholar
  31. F. Jiang, J. He, N. Navarro-Alvarez et al., “Elongation factor Tu and heat shock protein 70 are membrane-associated proteins from Mycoplasma ovipneumoniae capable of inducing strong immune response in mice,” Plos One, vol. 11, no. 8, 2016. View at: Google Scholar
  32. V. Pancholi, “Multifunctional alpha-enolase: its role in diseases,” Cellular and Molecular Life Sciences, vol. 58, no. 7, pp. 902–920, 2001. View at: Google Scholar
  33. M. W. Robinson, K. A. Buchtmann, C. Jenkins et al., “MHJ-0125 is an M42 glutamyl aminopeptidase that moonlights as a multifunctional adhesin on the surface of Mycoplasma hyopneumoniae,” Open Biology, vol. 3, no. 4, Article ID 130017, 2013. View at: Google Scholar
  34. S. A. Schreiner, A. Sokoli, K. M. Felder et al., “The surface-localised α-enolase of Mycoplasma suis is an adhesion protein,” Veterinary Microbiology, vol. 156, no. 1-2, pp. 88–95, 2012. View at: Publisher Site | Google Scholar
  35. M. G. Ferrarini, F. M. Siqueira, S. G. Mucha et al., “Insights on the virulence of swine respiratory tract mycoplasmas through genome-scale metabolic modeling,” BMC Genomics, vol. 17, no. 353, pp. 016–264, 2016. View at: Google Scholar
  36. M. G. Ferrarini, S. G. Mucha, D. Parrot et al., “Hydrogen peroxide production and myo-inositol metabolism as important traits for virulence of Mycoplasma hyopneumoniae,” Molecular Microbiology, vol. 108, no. 6, pp. 683–696, 2018. View at: Publisher Site | Google Scholar
  37. T. B. Reynolds, “Strategies for acquiring the phospholipid metabolite inositol in pathogenic bacteria, fungi and protozoa: making it and taking it,” Microbiology, vol. 155, no. 5, pp. 1386–1396, 2009. View at: Publisher Site | Google Scholar
  38. R. Veldhuizen, K. Nag, S. Orgeig, and F. Possmayer, “The role of lipids in pulmonary surfactant,” Biochimica et Biophysica Acta, vol. 19, pp. 2-3, 1998. View at: Google Scholar
  39. C. V. Bizarro and D. C. Schuck, “Purine and pyrimidine nucleotide metabolism in Mollicutes,” Genetics and Molecular Biology, vol. 30, no. 1, pp. 1271–1295, 2007. View at: Publisher Site | Google Scholar
  40. B. Hove-Jensen, “Mutation in the phosphoribosylpyrophosphate synthetase gene (prs) that results in simultaneous requirements for purine and pyrimidine nucleosides, nicotinamide nucleotide, histidine, and tryptophan in Escherichia coli.,” Journal of Bacteriology, vol. 170, no. 3, pp. 1148–1152, 1988. View at: Publisher Site | Google Scholar
  41. A. Kuehne, H. Emmert, J. Soehle et al., “Acute activation of oxidative pentose phosphate pathway as first-line response to oxidative stress in human skin cells,” Molecular Cell, vol. 59, no. 3, pp. 359–371, 2015. View at: Publisher Site | Google Scholar
  42. M. Shahbaaz, K. Bisetty, F. Ahmad, and M. I. Hassan, “In silico approaches for the identification of virulence candidates amongst hypothetical proteins of Mycoplasma pneumoniae 309,” Computational Biology and Chemistry, vol. 59, pp. 67–80, 2015. View at: Publisher Site | Google Scholar
  43. K. C. S. Bouh, F. Shareck, and S. Dea, “Monoclonal antibodies to Escherichia coli-expressed P46 and P65 membranous proteins for specific immunodetection of Mycoplasma hyopneumoniae in lungs of infected pigs,” Clinical and Diagnostic Laboratory Immunology, vol. 10, no. 3, pp. 459–468, 2003. View at: Google Scholar
  44. G. Browning, M. Marenda, A. Noormohammadi, and P. Markham, “The central role of lipoproteins in the pathogenesis of mycoplasmoses,” Veterinary Microbiology, vol. 153, no. 1-2, pp. 44–50, 2011. View at: Publisher Site | Google Scholar
  45. S. Futo, Y. Seto, M. Okada et al., “Recombinant 46-kilodalton surface antigen (P46) of Mycoplasma hyopneumoniae expressed in Escherichia coli can be used for early specific diagnosis of mycoplasmal pneumonia of swine by enzyme-linked immunosorbent assay,” Journal of Clinical Microbiology, vol. 33, no. 3, pp. 680–683, 1995. View at: Google Scholar
  46. J. A. Schmidt, G. F. Browning, and P. F. Markham, “Mycoplasma hyopneumoniae p65 surface lipoprotein is a lipolytic enzyme with a preference for shorter-chain fatty acids,” Journal of Bacteriology, vol. 186, no. 17, pp. 5790–5798, 2004. View at: Publisher Site | Google Scholar
  47. J. A. Schmidt, G. F. Browning, and P. F. Markham, “Mycoplasma hyopneumoniae mhp379 Is a Ca2+-dependent, sugar-nonspecific exonuclease exposed on the cell surface,” Journal of Bacteriology, vol. 189, no. 9, pp. 3414–3424, 2007. View at: Publisher Site | Google Scholar
  48. T. Stakenborg, J. Vicca, D. Maes et al., “Comparison of molecular techniques for the typing of Mycoplasma hyopneumoniae isolates,” Journal of Microbiological Methods, vol. 66, no. 2, pp. 263–275, 2006. View at: Publisher Site | Google Scholar
  49. V. G. Virginio, T. Gonchoroski, J. A. Paes, D. C. Schuck, A. Zaha, and H. B. Ferreira, “Immune responses elicited by Mycoplasma hyopneumoniae recombinant antigens and DNA constructs with potential for use in vaccination against porcine enzootic pneumonia,” Vaccine, vol. 32, no. 44, pp. 5832–5838, 2014. View at: Publisher Site | Google Scholar
  50. J. Meens, M. Selke, and G. Gerlach, “Identification and immunological characterization of conserved Mycoplasma hyopneumoniae lipoproteins Mhp378 and Mhp651,” Veterinary Microbiology, vol. 116, no. 1-3, pp. 85–95, 2006. View at: Publisher Site | Google Scholar
  51. W. Liu, D. Zhou, F. Yuan et al., “Surface proteins mhp390 (P68) contributes to cilium adherence and mediates inflammation and apoptosis in Mycoplasma hyopneumoniae,” Microbial Pathogenesis, vol. 126, pp. 92–100, 2019. View at: Publisher Site | Google Scholar
  52. F. C. Minion, E. J. Lefkowitz, M. L. Madsen, B. J. Cleary, S. M. Swartzell, and G. G. Mahairas, “The genome sequence of mycoplasma hyopneumoniae strain 232, the agent of swine mycoplasmosis,” Journal of Bacteriology, vol. 186, no. 21, pp. 7123–7133, 2004. View at: Publisher Site | Google Scholar
  53. J. L. Tacchi, B. B. A. Raymond, P. A. Haynes et al., “Post-translational processing targets functionally diverse proteins in Mycoplasma hyopneumoniae,” Open Biology, vol. 6, no. 2, Article ID 150210, 2016. View at: Google Scholar
  54. F. C. Minion, C. Adams, and T. Hsu, “R1 region of P97 mediates adherence of mycoplasma hyopneumoniae to swine cilia,” Infection and Immunity, vol. 68, no. 5, pp. 3056–3060, 2000. View at: Publisher Site | Google Scholar
  55. B. Henderson and A. Martin, “Bacterial virulence in the moonlight: multitasking bacterial moonlighting proteins are virulence determinants in infectious disease,” Infection and Immunity, vol. 79, no. 9, pp. 3476–3491, 2011. View at: Publisher Site | Google Scholar
  56. B. Henderson, S. Nair, J. Pallas, and M. A. Williams, “Fibronectin: a multidomain host adhesin targeted by bacterial fibronectin-binding proteins,” FEMS Microbiology Reviews, vol. 35, no. 1, pp. 147–200, 2011. View at: Publisher Site | Google Scholar
  57. W. Aruni, E. Vanterpool, D. Osbourne et al., “Sialidase and Sialoglycoproteases Can Modulate Virulence in Porphyromonas gingivalis,” Infection and Immunity, vol. 79, no. 7, pp. 2779–2791, 2011. View at: Publisher Site | Google Scholar
  58. R. E. Davis, R. Jomantiene, Y. Zhao, and E. L. Dally, “Folate biosynthesis pseudogenes, PsifolP and PsifolK, and an O-sialoglycoprotein endopeptidase gene homolog in the phytoplasma genome,” DNA and Cell Biology, vol. 22, no. 11, pp. 697–706, 2003. View at: Publisher Site | Google Scholar
  59. H. S. Garmory and R. W. Titball, “ATP-binding cassette transporters are targets for the development of antibacterial vaccines and therapies,” Infection and Immunity, vol. 72, no. 12, pp. 6757–6763, 2004. View at: Publisher Site | Google Scholar
  60. T. Berkelman, P. Garret-Engele, and N. E. Hoffman, “The pacL gene of Synechococcus sp. strain PCC 7942 encodes a Ca(2+)-transporting ATPase.,” Journal of Bacteriology, vol. 176, no. 14, pp. 4430–4436, 1994. View at: Publisher Site | Google Scholar
  61. K. L. Hein, P. Nissen, and J. P. Morth, “Purification, crystallization and preliminary crystallographic studies of a PacL homologue from Listeria monocytogenes,” Acta Crystallographica Section F: Structural Biology Communications, vol. 68, no. 4, pp. 424–427, 2012. View at: Publisher Site | Google Scholar
  62. L. Novoa-Aponte and C. Y. Soto Ospina, “Mycobacterium tuberculosis P-type ATPases: possible targets for drug or vaccine development,” BioMed Research International, vol. 2014, Article ID 296986, 9 pages, 2014. View at: Publisher Site | Google Scholar
  63. E. A. Groisman, K. Hollands, M. A. Kriner, E. Lee, S. Park, and M. H. Pontes, “Bacterial Mg2+ homeostasis, transport, and virulence,” Annual Review of Genetics, vol. 47, no. 1, pp. 625–646, 2013. View at: Publisher Site | Google Scholar
  64. K. O'Connor, S. A. Fletcher, and L. N. Csonka, “Increased expression of Mg2+ transport proteins enhances the survival of Salmonella enterica at high temperature,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 106, no. 41, pp. 17522–17527, 2009. View at: Publisher Site | Google Scholar
  65. S. Merino, R. Gavin, M. Altarriba et al., “The MgtE Mg2+ transport protein is involved in Aeromonas hydrophila adherence,” FEMS Microbiology Letters, vol. 198, no. 2, pp. 189–195. View at: Publisher Site | Google Scholar
  66. K. M. Keeney, J. A. Stuckey, and M. X. O'Riordan, “LplA1-dependent utilization of host lipoyl peptides enables Listeria cytosolic growth and virulence,” Molecular Microbiology, vol. 66, no. 3, pp. 758–770, 2007. View at: Publisher Site | Google Scholar
  67. M. O'Riordan, M. A. Moors, and D. A. Portnoy, “Listeria intracellular growth and virulence require host-derived lipoic acid,” Science, vol. 302, no. 5644, pp. 462–464, 2003. View at: Publisher Site | Google Scholar
  68. S. Donini, S. Garavaglia, D. M. Ferraris et al., “Biochemical and structural investigations on phosphoribosylpyrophosphate synthetase from mycobacterium smegmatis,” PLoS ONE, vol. 12, no. 4, 2017. View at: Google Scholar
  69. N. Haddad, M. Saramago, R. G. Matos, H. Prevost, and C. M. Arraiano, “Characterization of the biochemical properties of Campylobacter jejuni RNase III,” Bioscience Reports, vol. 33, no. 6, 2013. View at: Google Scholar
  70. M.-Y. Mao, Y.-M. Yang, K.-Z. Li et al., “The rnc gene promotes exopolysaccharide synthesis and represses the vicRKX gene expressions via microRNA-size small RNAs in Streptococcus mutans,” Frontiers in Microbiology, vol. 7, no. 687, 2016. View at: Google Scholar
  71. G. Singh, D. Sharma, V. Singh et al., “In silico functional elucidation of uncharacterized proteins of Chlamydia abortus strain LLG,” Future Science OA, vol. 3, no. 1, pp. 2016–2066, 2017. View at: Publisher Site | Google Scholar
  72. N. A. Okan, S. Chalabaev, T.-H. Kim, A. Fink, R. A. Ross, and D. L. Kasper, “Kdo hydrolase is required for francisella tularensis virulence and evasion of TLR2-mediated innate immunity,” mBio, vol. 4, no. 1, Article ID 00638-12, 2013. View at: Google Scholar
  73. B. G. Lawhorn, S. Y. Gerdes, and T. P. Begley, “A genetic screen for the identification of thiamin metabolic genes,” The Journal of Biological Chemistry, vol. 279, no. 42, pp. 43555–43559, 2004. View at: Publisher Site | Google Scholar
  74. Q. Cheng, L. Wu, R. Tu et al., “Mycoplasma fermentans deacetylase promotes mammalian cell stress tolerance,” Microbiological Research, vol. 201, pp. 1–11, 2017. View at: Publisher Site | Google Scholar
  75. D. M. Heithoff, R. L. Sinsheimer, D. A. Low, and M. J. Mahan, “An essential role for DNA adenine methylation in bacterial virulence,” Science, vol. 284, no. 5416, pp. 967–970. View at: Publisher Site | Google Scholar
  76. B. Kobe and A. V. Kajava, “The leucine-rich repeat as a protein recognition motif,” Current Opinion in Structural Biology, vol. 11, no. 6, pp. 725–732, 2001. View at: Publisher Site | Google Scholar
  77. N. J. Gay, L. C. Packman, M. A. Weldon, and J. C. J. Barna, “A leucine-rich repeat peptide derived from the Drosophila Toll receptor forms extended filaments with a β-sheet structure,” FEBS Letters, vol. 291, no. 1, pp. 87–91, 1991. View at: Publisher Site | Google Scholar
  78. X. Cheng, J. Nicolet, R. Miserez et al., “Characterization of the gene for an immunodominant 72 kDa lipoprotein of Mycoplasma mycoides subsp. mycoides small colony type,” Microbiology, vol. 142, no. 12, pp. 3515–3524, 1996. View at: Publisher Site | Google Scholar
  79. M. Monnerat, F. Thiaucourt, J. Nicolet, and J. Frey, “Comparative analysis of the lppA locus in Mycoplasma capricolum subsp. capricolum and Mycoplasma capricolum subsp. capripneumoniae,” Veterinary Microbiology, vol. 69, no. 3, pp. 157–172, 1999. View at: Publisher Site | Google Scholar
  80. G. Kikuchi, Y. Motokawa, T. Yoshida et al., “Glycine cleavage system, reaction mechanism, physiological significance, and hyperglycinemia,” Proceedings of the Japan Academy, Series B, vol. 84, no. 7, pp. 246–263, 2008. View at: Publisher Site | Google Scholar
  81. S. M. Twine, N. C. Mykytczuk, M. D. Petit et al., “In vivo proteomic analysis of the intracellular bacterial pathogen, Francisella tularensis, isolated from mouse spleen,” Biochemical and Biophysical Research Communications, vol. 345, no. 4, pp. 1621–1633, 2006. View at: Publisher Site | Google Scholar
  82. E. Allan and B. W. Wren, “Genes to genetic immunization: identification of bacterial vaccine candidates,” Methods, vol. 31, no. 3, pp. 193–198, 2003. View at: Publisher Site | Google Scholar
  83. H. B. Ferreira and L. A. Castro, “A preliminary survey of M. hyopneumoniae virulence factors based on comparative genomic analysis,” Genetics and Molecular Biology, vol. 30, no. 1, pp. 245–255, 2007. View at: Publisher Site | Google Scholar

Copyright © 2019 Sha Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

822 Views | 283 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.