BACKGROUND
the respiratory tract of swine is colonized by several bacteria, among which are three common mycoplasma species: mycoplasma flocculare, mycoplasma hyopneumoniae, and mycoplasma hyorhinis .
even though little information is available concerning the prevalence of bacteria in healthy lungs, these three mycoplasma species have been isolated from the respiratory tract of both healthy and diseased pigs .
while m. flocculare is usually described as a commensal bacterium  <cit> , m. hyopneumoniae and m. hyorhinis are considered pathogenic.
enzootic pneumonia, caused by m. hyopneumoniae, is widespread in pig populations and is considered a major cause of economic loss in the pig industry  <cit> ; m. hyorhinis is frequently present in cases of polyserositis and arthritis and m. flocculare has high prevalence in swine herds worldwide, but up to date, no disease has been per se associated with this species  <cit> .
in addition to these mycoplasmas, mycoplasma hyosynoviae, the primary agent of non-purulent arthritis, can occasionally colonize in the lower respiratory tract when pneumonic lesions pre-exist  <cit> .
besides mycoplasmal pneumonia, the porcine respiratory disease complex  has emerged as an economically signigicant respiratory disorder characterized by the slow growth, fever, cough, loss of appetite, lethargy and dyspnea in pigs  <cit> .
even though many species are related to prdc, it is essential to note that enzootic pneumonia caused by m. hyopneumoniae is by far the most costly disease in pig industry, and this bacteria is usually seen as an essential component to the successful establishment of a pathogenic community in the host  <cit> .
also, m. hyopneumoniae infections take longer to cause lesions and to be successfully eliminated than infections from other pathogens  <cit> .
while mycoplasmal diseases in swine have been extensively studied, their causative agents have not been explored from a mathematical and computational point of view, mostly because their genome sequences were not available until recently .
although recent studies have placed m. hyopneumoniae, m. flocculare and m. hyorhinis in close proximity within the hyopneumoniae clade by phylogenomic analysis  <cit> , which corroborates with their high 16s rrna sequence similarity  <cit> , it is not yet clear what causes the specific pathogenicity or lack thereof in each of them.
this elevated genomic resemblance combined with their different levels of pathogenicity is an indication that these species, as for most mycoplasmas, have unknown mechanisms of virulence and differential expression.
pathogenic determinants such as adhesion to the host cell and evasion from the immune response have already been well-described in the literature for both m. hyopneumoniae and m. hyorhinis .
the presence of a capsule in m. hyopneumoniae has been reported to be important for the interaction with the host cells in one single study  <cit> .
several studies also show that immunosuppressed animals experimentally infected with mycoplasma species develop less severe microscopic lesions of pneumonia if compared to normal animals .
this means that it is possible in some cases that a strong response from the host immune system might be the primary cause of pathogenesis.
however, up to date, it is not possible to draw any further conclusions due to lack of experimental data.
even if these topics are of utter importance for understanding swine respiratory tract mycoplasmal diseases, what has yet to be better understood is the direct participation of metabolism in the development of them.
for instance, although adhesion factors are related to pathogenicity, m. hyopneumoniae and m. flocculare harbor similar sets of adhesion proteins  <cit> , and have been shown to adhere to cilia in a similar way  <cit> .
thus, the ability of m. hyopneumoniae to cause disease if compared to m. flocculare might not be directly related to adhesion.
furthermore, the genome sizes of mycoplasma spp.
range from  <dig> kb  to more than  <dig>  kb , representing an important example of genome reduction  during the evolutionary process.
it is possible that in an initial symbiotic phase, the host provided a broad range of metabolites for these bacteria.
this, together with the ability of the bacteria to uptake such compounds, made several activities dispensable for the bacterial life.
over the course of evolution, these bacteria would have lost some of the genes that became unnecessary for life in an environment conditioned by another genome  <cit> .
in this way, we aimed at studying the reduced metabolism through genome-scale metabolic model reconstructions of m. hyopneumoniae, m. hyorhinis and m. flocculare to better understand their different life-styles.
based on the reconstructed networks, we propose that one of the mechanisms that may explain why m. hyorhinis and m. hyopneumoniae are pathogenic while m. flocculare is not, is their ability to use glycerol as a carbon source, thus enabling the production of the highly toxic hydrogen peroxide.
this trait may be directly involved in cytotoxicity, as already described for two lung pathogenic mycoplasmas, namely mycoplasma pneumoniae in human  <cit>  and mycoplasma mycoides subsp.
mycoides in ruminants  <cit> .
additionally, growth rates and fitness have been correlated in the past to virulence in several other organisms  and may be a key factor for the difference in pathogenicity among these three mycoplasmas.
indeed, the in silico models show that m. hyorhinis and m. hyopneumoniae have extra sets of enzymes that may help them overgrow the non-pathogenic species m. flocculare.
m. hyorhinis seems to have a wider carbohydrate uptake range, and probably for this reason is seen as a well-known contaminant of cell cultures  <cit> .
m. hyopneumoniae, on the other hand, might use myo-inositol as a carbon source producing acetyl-coenzyme-a as a byproduct.
this feature is of particular interest for the production of cofactor coenzyme-a in this species, since most enzymes from the biosynthetic pathway are missing in all species.
nevertheless, there is little metabolic experimental data available for the three species, which makes the reconstruction of a reliable metabolic model an extremely time-consuming work.
together with the fact that these mycoplasmas are generally grown in complex media, with high serum concentrations, we needed additional experimental data in order to compare the in silico reconstructed metabolic networks with the in vivo metabolic characteristics of the three species.
for this reason, we also performed nuclear magnetic resonance spectroscopy  analyses to detect metabolites consumed and produced in both complex and defined media.
these experiments corroborated with the reconstructed models and suggested two new features in particular:  the uptake of myo-inositol in m. hyopneumoniae might be related to a higher acetate production, and  m. hyorhinis showed a surprisingly reduced ability to convert pyruvate to acetate in the growth conditions used in this study.
all these in silico and in vivo metabolic differences might influence the different levels of pathogenicity in each of the species studied here.
we present in this work the metabolic networks reconstructed from the annotated genomes of the three species and the comparison done with growth rates and metabolomic experiments performed in vitro in order to better understand the basis of the pathologies caused by these bacteria, which might help prevent their development in the future.
methods
modeling methods
network reconstruction and refinement
we reconstructed the metabolic networks for all available strains isolated from swine so far of m. hyopneumoniae , m. hyorhinis  and m. flocculare .
the semi-automated reconstructions were generated by the pathologic tool from the pathway tools software  <cit>  using the complete genomes available online.
from now on, the species strains will be abbreviated according to table  <dig>
m. hyopneumoniae
m. hyorhinis
m. flocculare
notes:
na: not available
⋆a: strains available only in brazil are protected strains which are property of the government and cannot be used outside authorised laboratories.
strains available in france were purchased from the atcc repository; these strains were also available for testing in the laboratory in brazil
⋆b: m. hyorhinis strain gdl- <dig> was retrieved from a contaminated cell line and, to our knowledge it has never been reported neither as pathogenic nor as an attenuated strain
⋆c: m. hyorhinis strain atcc  <dig> genome is available online but not published in any paper up to now.
pathway tools automatically associated genes with reactions, based on the annotation names, gene ontology  terms and enzyme code  numbers contained in the genbank files.
the software assembled reactions into pathways by comparing them with the reference database, metacyc.
the software automatically added missing reactions from the reference database to the model, creating many orphan reactions .
in a gapfilling step, we used the pathway hole filler tool, included in the pathway tools software, to identify possible candidate genes associated with these orphan reactions.
orphan reactions in the metabolic network were allowed only to allow known functionalities of the organisms or biomass production in the reconstructed models.
we incorporated the minimal number of orphan reactions into the models to avoid adding a behavior that was never described.
in our work, the refinement of the networks was made in a subsystems approach  <cit>  simultaneously on all organisms.
as a result of different genomic annotations, slightly different reactions arose throughout the automatic models.
these reactions were checked for consistency, and after validation of only one of them, the duplicates were deleted.
gene-protein-reaction  associations were systematically validated or included based on experimental data, information provided by the pfam functional domain database  <cit> , synteny analysis  and reciprocal sequence homology  searches, using at first high e-value cutoffs of at least 10− <dig>  <cit> .
whenever no hits were achieved with these parameters, more relaxed parameters were assumed .
homolog proteins with equivalent functional domains were assigned as isozymes to a particular reaction, while proteins with distinct functional domains were assigned as subunits of a multi-protein complex.
non-metabolic reactions such as dna polymerization, protein synthesis and rna synthesis were explicitly deleted from the network but implicitly included in the biomass assembly.
generic reactions were either specified or excluded from the models.
the resulting models were further refined in accordance with a detailed protocol from palsson and thiele  <cit> .
reactions were computationally balanced for mass and charge, while cofactor usage was determined based on literature data for closely related species .
reaction directionality was thermodynamically checked and validated.
reaction directionality was determined based on the component contribution method  <cit> , which extends the group contribution method  <cit>  and achieves a significant improvement in the accuracy of the estimations of standard gibbs energies.
an online search and calculation interface called equilibrator along with metabolite and reaction thermodynamic databases are available at  <cit> .
heuristic rules were used to improve the directionality assignment.
we analyzed the network topology to identify compounds that were only produced or consumed in the network, the so-called dead-end metabolites .
whenever a dem was found, we either  added or validated an orphan reaction to reconnect it to the rest of the network, or  removed the reaction from the model when both substrate and product were disconnected from the network and/or did not affect the overall metabolism.
transporters were predicted at first from genome annotation using the transport identification parser, from pathway tools  <cit> .
transporters were predicted mainly based on sequence similarity to other known transporters, since experimental data for the three mycoplasmas were not available.
specific transport reactions and exchange reactions required for production of biomass components were manually added to the final version of the networks.
biomass composition and biomass equation
in order to simulate growth, we had first to estimate the average cell composition of these mycoplasmas .
to this end, we used a general macromolecule mycoplasma cell composition from razin and collaborators  <cit> .
depending on the species, the authors showed that the solid residue contained 54– <dig> % of proteins, 12– <dig> % of lipids, 3– <dig> % of carbohydrates, 8– <dig> % of rna, and 4– <dig> % of dna.
the membranes comprised around  <dig> % of dry weight of the organisms and contained 47– <dig> % of proteins, 35– <dig> % lipids, 4– <dig> % carbohydrates and small amounts of dna and rna.
since no information on metal ions and cofactors was available, we included them quantitatively based on the metabolic networks reconstructed for related mycoplasmas  <cit> .
membrane and lipid components were added based on literature composition of the selected species.
it is important to note that the only difference in composition between the species reported so far was related to the presence of glycolipids in m. hyopneumoniae and m. flocculare, which were reported to be absent in m. hyorhinis .
based on the previous information, we assumed the following fractions of macromolecules:  <dig> % proteins,  <dig> % lipids,  <dig>  % carbohydrates,  <dig> % rna,  <dig> % dna, and  <dig>  % of ions and cofactors.
to create a biomass elemental formula, we took into account the percentage contribution of each of the components to the overall cell.
the biomass reaction represented the drain of these components into biomass production.
it implicitly assembled dna replication, rna transcription, and protein synthesis into one single reaction.
amino acids were indirectly included in the biomass reaction: charged trnas were accounted as substrates and uncharged trnas as products.
growth and non-growth associated maintenance  were estimated based on the literature and were manually added to the models  <cit> .
the details for the complete assembly of the biomass reaction can be found in additional file  <dig>
model validation and fba analysis
the metabolic networks were exported as a mathematical systems biology markup language  models  <cit> .
they were uploaded to the cobra-toolbox v <dig>  matlab extension  <cit>  for flux balance analysis  testing.
minimum-maximum flux constraints were imposed based on literature information  <cit> .
growth simulations were achieved using biomass production as the fba objective function.
experimental methods
swine mycoplasmas cultivation
complex media comprised  friis media  <cit>  , and  a commercial mycoplasma broth .
defined media cultivation tests were performed in france and comprised  a medium described for m. pneumoniae strain  <dig> by yus and collaborators , and  commercial media cmrl with no glutathione .
since we had no information on metabolism of swine mycoplasmas, we decided to supplement the defined yus medium with all amino acids .
we also supplemented the cmrl- <dig> medium with other peptone and/or other cofactors .
information on the composition of all defined media is available in additional file 1: table s <dig>  cells were cultivated at  <dig> °c for different time periods, under gentle agitation .
cell concentration estimation and viability by color changing units  measurement
cell growth and viability was measured with triplicate time-matched samples of cells and culture media for ccu as described by stemke and robertson  <cit> .
viability of cells is visible by a change in medium color from red to yellow.
for cell concentration measurements, the cultures were subjected to a series of 10-fold dilutions in complex media and  <dig> ccu/ml was defined as the highest dilution of cells able to change the medium color  <cit> .
samples for nmr spectroscopy
nmr was performed with complex friis medium  and defined yus+ medium .
the medium was collected at the following time intervals:  <dig> h,  <dig> h,  <dig> h,  <dig> h,  <dig> h, and  <dig> h for friis medium and  <dig> h,  <dig> h,  <dig> h,  <dig> h,  <dig> h,  <dig> h and  <dig> h for yus+ medium.
cells were separated from growth media through sedimentation at  <dig> g for nmr analysis.
samples consisted of biological triplicates in complex medium and biological duplicates in defined medium.
since we had slight different time intervals, care was taken by not directly comparing the results from different experiments.
nmr analysis
sample preparation was as follows:  <dig> μl of a mixture containing  <dig>  mkh2po <dig> phosphate buffer  in d2o with  <dig> mmnan <dig> and  <dig>  % trimethylsilyl propionate  was added to  <dig> μl supernatant samples.
both solutions were mixed thoroughly and  <dig> μl were then transferred to  <dig> mm nmr tubes and sorted in 96-tubes racks.
all nmr experiments were carried out on a bruker  <dig> mhz nmr spectrometer equipped with a  <dig> mm txi probe and a samplejet autosampler, enabling high-throughput data acquisition for large collections of samples.
the temperature was controlled at  <dig> °c throughout the experiments, and the samples were kept refrigerated at  <dig> °c during a waiting time of less than  <dig> h in the autosampler, before the nmr analysis.
standard 1h 1d nmr pulse sequence nuclear overhauser effect spectroscopy  with z-gradient and carr-purcell-meiboom-gill  with water presaturation  were applied on each sample to obtain the corresponding metabolic profiles.
a total of  <dig> transient free induction decays  were collected for each experiment with a spectral width of  <dig> ppm.
the relaxation delay was set to  <dig> s. the noesy mixing time was set to  <dig> ms. the total acquisition time of each sample was  <dig> min  <dig> s. the cpmg spin-echo delay was set to  <dig> ms, for a total filter of  <dig> ms, allowing an efficient attenuation of the lipid and protein nmr signals.
the 90° pulse length was automatically calibrated for each sample at around  <dig> μs.
data processing: all fids were multiplied by an exponential function corresponding to a  <dig>  hz line-broadening factor, prior fourier transformation.
1h-nmr spectra were manually phased and referenced to the tsp signal  using topspin  <dig>  .
extraction of a data matrix for multivariate statistical analysis from the 1h nmr profiles was done using the statistics toolbox of amix .
spectra were integrated from  <dig>  to  <dig> ppm at a step of  <dig>  ppm but excluding the regions of residual water at  <dig> - <dig>  ppm.
no normalization of the intensity was performed.
the resulting data matrix contains  <dig> nmr variables.
multivariate data analysis: principal component analysis  and hierarchical clustering analysis  were performed using simca-p  <dig>  with scaling based on the pareto method.
metabolites identification and quantification: metabolite identification was achieved by comparing spectra with databases such as hmdb  <cit> .
identification of the metabolites was further verified with homonuclear and heteronuclear 2d nmr experiments such as 1h−13c hsqc, 1h−1h tocsy and j-resolved experiments.
absolute quantification of the metabolites was performed using chenomx nmr suite .
RESULTS
model reconstruction and refinement
based on the published genomes of  <dig> strains of m. hyopneumoniae,  <dig> strains of m. hyorhinis and  <dig> strains of m. flocculare, we reconstructed  <dig> genome-scale metabolic models: one model for each strain separately, one for each species and a pan-reconstruction for all three species .
the semi-automated reconstructions were manually refined according to the description given in the “methods” section.
the pan-swine mycoplasma network  was initially composed by  <dig> reactions.
duplicate reactions arose from the fact that pathway tools is based only on the annotations contained in the genbank files.
we removed generic and duplicate reactions from the models and replaced them with the specific and validated ones.
non-metabolic reactions were also excluded at this point, along with absent capabilities of these mycoplasmas, such as heme, quinone or cytochrome dependent reactions  <cit> .
the networks were also tested for the presence of dem and the ability to produce all biomass precursors.
dems were analyzed on a case-by-case manner.
from initially  <dig> dems,  <dig> disconnected metabolites  were excluded from the models:  <dig> dems came from  <dig> spontaneous reactions;  <dig> were carbohydrate substrates derived from wide range transport reactions and were not used by any other reaction in the network;  <dig> were excluded as they did not interfere with the overall metabolism and/or the enzyme had already been assigned to one or many other reactions.
orphan reactions were excluded unless they were essential for biomass growth.
all excluded reactions  and all dems  are listed in additional file 1: tables s2c and s2d.
it is interesting to point out that other unknown enzymes  from these organisms may indeed use some of the excluded dems; however, since we have no experimental evidence at present, we could not assess their interference in the metabolic models.
the remaining dems consisted mainly of cofactors and biomass precursors  disconnected from the rest of the network.
they were solved along with the biomass precursor check;  <dig> reactions were added to allow growth:  <dig> biomass reaction,  <dig> enzymatic reactions,  <dig> drain-synthetic reaction,  <dig> transport reactions and  <dig> spontaneous reactions.
transport reactions were considered as such even if an enzymatic activity was also present .
in order to correctly assign possible and specific transporters, we performed an extensive literature search and reciprocal blast alignments to characterized transporters in other species.
these results have to be experimentally confirmed, since assignment of the correct substrate based only on sequence homology remains an open problem in genomic annotation.
we also changed inconsistencies of reversibility and cofactor usage.
even after all efforts, four dems still remained: ttp, hexulose-6-phosphate, deoxyinosine and xanthosine-5-phosphate.
after the final addition of  <dig> exchange reactions, we ran fba tests to check the consistency of all models.
if a reaction was essential for biomass growth and no homolog gene was found in the genome, an orphan reaction was validated or added to the reconstruction for modeling reasons only.
the resulting refined pan-network had a total of  <dig> reactions and  <dig> gpr associations.
a comparison between the species models may be found in fig.
<dig>  these results indicate that all strains from all species are indeed metabolically similar.
fig.
<dig> venn diagrams representing the comparison of refined networks.
numbers represent the exclusive and common reactions present in the refined networks  between species, and  between selected strains of m. hyopneumoniae.
this analysis shows that most of the metabolism is common to all organisms.
mhr: m. hyorhinis; mhp: m. hyopneumoniae; mfl: m. flocculare
the overall characteristics of each reconstruction can be seen in table  <dig>  a list of all reactions and metabolites added to the models along with the corresponding genes can be seen in additional file 1: table s <dig>
notes:
∗a: gprs are gene-proteins-reaction associations present in each genome.
for imfmhp, imfmhr, imfmfl and imfpan that account for more than one species, we added a gpr to the species model when any of the strains harbored a gene responsible for a specific activity
∗b: reactions were considered as transport reactions even if the transporter was capable of performing a concomitant enzymatic activity
∗c: spontaneous reactions included diffusion of small molecules and spontaneous conversions
∗d: total number of reactions included also all exchange reactions needed for the mathematical modeling
∗e: genes in the pan-network and species models were added synthetically , to enable a reaction essentiality analysis
biomass composition and biomass equation
the biomass equation drained all precursors  into biomass.
biomass composition according to the “methods” section along with a detailed description of the assembly of the biomass reaction can be found in additional file 1: table s <dig> and additional file  <dig>  considering the percentage contribution of each of the components to the overall cell, an approximate biomass elementary composition for m. hyorhinis was computed as follows: ch <dig> o <dig> n <dig> p <dig> s <dig> , with traces of calcium, chlorum, cobalt, copper, iron, potassium, magnesium, manganese, molibdenium and zync and molecular weight  on a c-mole basis of  <dig>  g/ c−mol.
we had a slightly different composition for m. hyopneumoniae and m. flocculare: ch <dig> o <dig> n <dig> p <dig> s <dig> , with the previously mentioned trace compounds and mw on a c-mole basis of  <dig>  g/ c−mol.
metabolism of m. hyorhinis, m. hyopneumoniae and m. flocculare
similar to most mycoplasma species studied so far , all reconstructed networks exhibit low connectivity due to the simplicity of the biological model.
out of the  <dig> reactions in the final pan-network,  <dig> had in at least one species a gpr association.
from these,  <dig> were common to all species .
the overall metabolism from the models reconstructed consisted of  <dig> distinct subsystems: amino sugar metabolism, amino acid metabolism, carbohydrate metabolism , cofactor metabolism, lipid metabolism, and nucleotide metabolism .
while all enzymes were present in glycolysis, most metabolic pathways had major enzyme gaps.
in the product and cofactor metabolism, for instance, gaps accounted for up to  <dig> % of the reactions.
m. flocculare was the only species that did not show any exclusive metabolic activities in the models, although this could not be verified in vivo due to lack of biochemical studies for this species.
m. hyopneumoniae had  <dig> exclusive reactions, linked to the myo-inositol metabolism and alcohol dehydrogenase activity.
m. hyorhinis had  <dig> exclusive reactions; most of them corresponded to carbohydrate metabolism.
a global model comparing each species enzymatic capabilities can be seen in figs.
<dig> and  <dig>  the complete list of metabolite abbreviations and ec numbers for these two figures can be found in additional file  <dig>
fig.
<dig> distribution of the model reactions in the subsystems.
the  <dig> reactions present in the model imfpan were separated into  biological subsystems and  further into reaction types, with the exclusion of exchange reactions in this analysis
fig.
<dig> global energy and carbohydrate complete models.
metabolites are depicted in dark green and separate enzymatic activities for m. hyorhinis, m. hyopneumoniae and m. flocculare can be seen in yellow, pink and blue, respectively.
whenever an enzyme is missing from the three species, the enzyme rectangle is depicted in grey.
complete list of metabolite abbreviations and ec numbers can be found in additional file 4
fig.
<dig> lipid, amino acid, nucleotide and cofactor complete models.
metabolites are depicted in dark green and separate enzymatic activities for m. hyorhinis, m. hyopneumoniae and m. flocculare can be seen in yellow, pink and blue, respectively.
whenever an enzyme is missing from the three species, the enzyme rectangle is depicted in grey.
complete list of metabolite abbreviations and ec numbers can be found in additional file 4
transporters
transporter assignment was made in silico based on sequence similarity to known transporters.
we were able to identify in all species complete putative abc transporters for multidrug/toxin efflux and for import of sugars, oligopeptides, peptides, spermidine/putrescine, phosphonate/phosphate/thiamine, cobalt, manganese/zinc and glycerol.
in the search of a possible nucleotide transporter, we came across with an unpublished work from lee  <cit>  in which a system previously annotated as sugar abc transporter was experimentally validated as involved in nucleoside uptake in mycoplasma bovis.
these results may be highly speculative, and need further verification, but we were able to find homolog genes in the genomes of all species.
m. hyopneumoniae seems to have an extra abc transporter proposed for myo-inositol and m. hyorhinis has a unique predicted abc transporter for maltose/maltodextrins.
three complete phospho-transferase transport systems  were common to all species.
although specificity of pts is difficult to ascertain, we have tried to compare conserved domains from validated pts substrate-specific enzymes ii  and candidates from the three species.
genes coding enzyme i and phosphocarrier were found in all strains, except in mhr_hub <dig>  from sequence similarity and gene context we were able to annotate the correct gene.
based on these in silico predictions from domain retrieval, we propose the following: one pts non-specific for sugar, one with specificity for fructose and another for mannitol.
m. hyopneumoniae and m. flocculare also share two extra complete systems: one possibly for ascorbate and one for n-acetylglucosamine .
m. hyorhinis and m. hyopneumoniae seem to have an extra iib component specific for glucose, and along with the other components of the general sugar pts were proposed to form a complete glucose pts.
the protein glpu, coded by gene mpn <dig> in m. pneumoniae, was recently described to act on the uptake of glycerophosphodiesters  <cit> .
we found homologs of this gene in all strains of m. hyopneumoniae and m. hyorhinis adjacent to the gene coding for the enzyme responsible for the metabolization of these substrates   <cit> .
in mhr_hub <dig>  this gene was annotated as a pseudo gene, but from gene context, we were able to validate it as glpu.
all species also possess several unspecific amino acid permeases, a glycerol facilitator protein , which is thought to be less efficient than the glycerol abc transporter   <cit> ) and a major facilitator protein  with unknown specificity.
we could also find common transporters for cobalt and magnesium, chromate, magnesium, potassium, zinc and a cation atpase.
m. hyorhinis has two extra transport systems: one for sialic acid  and one sodium/phosphate co-transporter.
details of the transport search throughout all species is summarized in additional file  <dig>
carbohydrate and energy metabolism
glycolysis
as previously described in  <cit> , all genes for the glycolysis exist in all strains of all species.
m. hyorhinis may be able to convert dihydroxyacetone  directly to dihydroxyacetone-phosphate  with the action of dha kinase , but the mechanism of dha uptake is not clear.
pentose phosphate
in agreement with previously described literature  <cit> , we found no enzymes for the oxidative branch of the pentose phosphate pathway in the studied species, turning ribulose-5-phosphate into a dem.
this problem was overcome by the assumption that the reaction ribulose-phosphate-3-epimerase  can be reversible, as reported in lactobacillus casei  <cit> .
out of the  <dig> reactions present in the final models, we found homologs for  <dig> in all species.
although a gene coding for transaldolase  is missing in all species, we included this reaction in the models to prevent both sedoheptulose-7-phosphate and erithrose-4-phosphate from becoming dem.
ascorbate metabolism
only one enzyme is missing from the typical pathway: l-ascorbate-6-phosphate-lactonase .
this enzyme is responsible for the turnover of l-ascorbate-6-phosphate to 3-keto-gulonate-6-phosphate .
from the gene context, mhp_7448_ <dig> from mhp_ <dig> and its homologs from all other organisms may encode the missing enzyme.
the gene was annotated as a conserved hypothetical protein with a proposed phosphotriesterase activity.
uptake of other carbohydrates
while the uptake of fructose, mannitol and mannose is similar in all species and they all can be fed into glycolysis as fructose-6-phosphate  and dhap, m. hyorhinis has an extra set of enzymes enabling the conversion of fructose-1-phosphate into fructose- <dig> -biphosphate by the action of 1-phosphofructokinase .
only in m. hyopneumoniae the conversion of glyceraldehyde into glycerol might be possible through the moonlight activity of alcohol dehydrogenase , although up to date no experimental evidence can sustain that assumption.
the presence of a unique transcriptional unit  with carbohydrate metabolism-related genes in m. hyorhinis probably enables the uptake and metabolization of isomaltose , maltose , trehalose  and sucrose  into either glucose-6-phosphate or f6p.
this may be related to the fact that m. hyorhinis can overgrow the other two species in several growth media.
myo-inositol metabolism
a tu for the myo-inositol catabolism is present in all m. hyopneumoniae strains, with the exception of the enzyme 5-dehydro-2-deoxy-phosphogluconate aldolase .
the gene encoding this enzyme in other organisms is similar to the fructose-biphosphate aldolase  from glycolysis .
since there are two copies annotated for fba in m. hyopneumoniae, we proposed one of them as candidate for this activity based on sequence alignments with both genes from bacillus subtilis.
inositol can be used as a carbon source and also produces acetyl coenzyme-a , which can be a source of cofactor coenzyme-a .
amino sugar metabolism
m. hyopneumoniae and m. flocculare models could uptake and convert glcnac to n-acetyl-glucosamine-6-phosphate ; m. hyorhinis models imported sialic acid and converted it to glcnac6p in three steps that are unique to m. hyorhinis .
at first, the genes coding for these three enzymes were absent from mhr_hub <dig> and mhr_gdl <dig>  however, based on chromosome alignment and synteny with other genes, we were able to find the regions and correctly annotate all genes.
pyruvate metabolism
out of fifteen enzymatic reactions in this pathway, only three were not found in any organism.
conversion of pyruvate to lactate and acetate were possible in all species.
from literature data, we knew that m. hyopneumoniae was able to produce oxaloacetate, malate and 2-oxoglutarate   <cit> .
decarboxylation of 2kg to succinyl coenzyme-a  has been previously reported in m. hyorhinis extracts  <cit> .
although we did not find any candidate for the activity of aspartate transaminase , gap-filling proposed that this activity could be performed by the pyruvate dehydrogenase complex.
moreover, cordwell and collaborators  <cit>  suggested that in mollicutes lactate dehydrogenase  could also function as malate dehydrogenase .
m. hyopneumoniae is the only organism among the three to have a gene coding for an alcohol dehydrogenase, probably enabling the turnover of accoa to acetaldehyde  and ethanol .
conversion of pyruvate to formate by pyruvate formate-lyase , malate to fumarate by fumarate hydratase  and succoa to succinate by succinyl-coa synthetase  were added to the models in a final step to accommodate the results from the metabolomics experiments.
lipid metabolism
the lipid metabolism in the models of the three species included the uptake of glycerol, choline, glycerophosphodiesters and fatty acids to the production of the biomass precursors cardiolipin,  <dig> -diacyl-sn-glycerol, phosphatidylglycerol, phosphatidylcholine and galactosyl-diacylglycerols .
extracellular cholesterol and sphyngomyelin were incorporated unmodified directly into biomass.
while m. hyopneumoniae seems to have four different ways to uptake and metabolize glycerol  m. hyorhinis lacks one and m. flocculare lacks two of them .
the turnover of glycerol-3-phosphate into dhap by the action of glycerol-3-phosphate oxidase  allows the usage of glycerol as the sole carbon source, with the production of highly toxic hydrogen peroxide.
this was only possible in the m. hyorhinis and m. hyopneumoniae models.
the enzyme responsible for the turnover of glycerol-3-phosphate into acyl-sn-glycerol-3-phosphate  was found only in the genomes of m. hyorhinis.
although this can be seen as a possible difference between the species, we had to add the reaction to the models from the other species to enable growth.
it is possible that the other two species import this metabolite directly from the media, however further experiments are required to confirm this hypothesis.
no homologs for the enzyme that produces phosphatidylglycerol  were found; however, all mollicutes definitely synthesize this metabolite, and hence must use a hitherto undetected enzyme for this step  <cit> .
the acyl-carrier protein  was only found in the genomes of m. hyorhinis and may act as a fatty acid and coa donor in these organisms.
orphan reactions for the production of glycolipids  were added to models for m. hyopneumoniae andm.
flocculare based on the presence of these metabolites in vivo  <cit> .
amino acid metabolism
the import of amino acids was added in two forms: oligopeptide import and cleavage or single amino acid import.
trna charging accounts for most of the reactions in this pathway .
production of the biomass precursor s-adenosyl-methionine from methionine is possible in all strains .
moreover, since the enzyme adenosylhomocysteine nucleosidase  is present in all species, we added two reactions to restore connectivity between these two subpathways .
this addition resulted in the production of two metabolites essential for quorum sensing and cell communication in other species: s-rybosyl-homocysteine and autoinducer- <dig>
nucleotide metabolism
the mycoplasmas in this study cannot synthesize de novo purines and pyrimidines; therefore, they have only salvage pathways and interconversions to supply the cell with nucleic acid precursors.
the three species have the same enzymatic capabilities, except for the presence of thymidylate synthetase  in m. hyorhinis, which is also important for cofactor metabolism and is responsible for the conversion of dump and  <dig> -methylene-tetrahydrofolate  to dtmp and dihydrofolate.
overall, the nucleotide metabolism consists in the uptake of guanine, adenine, uracil, thymine and cytidine and produces all deoxy-ribonucleotides  and ribonucleotides .
cofactor metabolism
around  <dig> % of the enzymes that comprise the cofactor metabolism of the reconstructed models were not found in any species.
for instance, thiamine-pyrophosphate is imported unchanged directly into biomass, while pyridoxal is imported and converted to the cofactors pyridoxal-phosphate by the action of pyridoxal kinase .
from the pantothenate/coenzyme-a metabolism, we only found homologs for the conversion of 4’-phosphopantetheine into dephospho-coenzyme-a by pantetheine-phosphate adenylyltransferase .
three forms of folate are incorporated to the biomass: methf,  <dig> , <dig> -tetra-hydrofolate and 10-formyl-tetra-hydrofolate.
we assumed that folate was imported and converted to these three biomass cofactors by: dihydrofolate reductase , serine hydroxymethyltransferase  and formate-tetrahydrofolate ligase .
this may also account for the better fitness of m. hyorhinis when compared to the other two species.
although no homologs for nad kinase were found, two reactions were added  to allow the presence of nadp and nadph, which are cofactors for several essential reactions.
the nadph and nadh produced from the degradation of glucose and other carbohydrates are recycled by the action of two enzymes: nadh oxidase  and thioredoxin disulfide-reductase .
mycoplasmas cultivation
we cultivated the three species of mycoplasmas in complex and defined media to gather experimental information and compare to in silico growth .
in the cmrl+ medium, m. hyopneumoniae and m. flocculare remained viable only when peptone was present .
the presence of peptone appeared to have a negative effect on m. hyorhinis growth, but more tests should be performed to verify this hypothesis.
yus+ medium did not allow proliferation in any species, but maintained cell concentration and viability even after  <dig> days of culture .
these results are in agreement with a previous work performed by bertin and colleagues  <cit>  in m. mycoides subsp.
mycoides that shows that cmrl- <dig> contains all components to support cellular metabolism but not growth.
here, the supplemented versions cmrl+ or cmrl+/pep allowed proliferation only for m. hyorhinis during the first  <dig> h of growth.
it seems that after this period of time, one or more essential metabolites initially present were no longer available.
by comparing the models reconstructed in this work with the media composition tested for the three species, it might seem plausible that key cofactors might not be delivered to mycoplasmas in the correct form.
for example, thiamine pyrophosphate should be directly delivered to all species instead of its precursor thiamine, and pyridoxal-5-phosphate might be the cofactor of choice for the media composition in place of pyridoxal.
a list of predicted biomass precursors from the models versus the actual precursors found in the yus+ medium can be seen in additional file 1: table s <dig>  this new proposed defined medium will be submitted in the near future to cultivation tests with m. hyorhinis, m. hyopneumoniae and m. flocculare.
fig.
<dig> cultivation curves in defined and complex media by species.
cell concentrations were estimated by the ccu method and error bars were calculated as the standard deviation between triplicate time-matched samples.
as expected, the three species had better growth rates in the complex medium than in defined media.
m. hyopneumoniae strain j ; m. flocculare strain  <dig> ; m. hyorhinis strain atcc <dig>
metabolomics of swine mycoplasma
nmr analysis of the culture media was performed in order to detect possible differences in the metabolism of the three species.
all data measured in defined and complex media can be seen in additional file 1: tables s6a and s6b.
the major differences in both complex and defined media were related to the metabolism of pyruvate : m. hyorhinis produced higher quantities of pyruvate at the end of  <dig> h in complex medium.
as a result, pyruvate conversion to acetate was detected in low quantities for this species in both media.
m. hyopneumoniae and m. flocculare, on the other hand, produced higher amounts of acetate .
low amounts of formate, fumarate and succinate were also detected, and this indicates the presence of genes encoding the enzymes pyruvate formate-lyase , fumarate hydratase  and succinyl-coa synthetase  in these genomes.
as previously mentioned, all models were thus modified to accommodate these activities.
formate production seems to be independent from acetate production in m. hyopneumoniae and m. flocculare.
however, the similar low concentrations of acetate and formate in both growth media from m. hyorhinis might be an indication that both are produced concomitantly by the action of pyruvate formate-lyase , phosphate acetyl-trasnferase  and acetate kinase .
fig.
<dig> distinct products of the metabolism of pyruvate from growth in complex friis medium and yus+ medium of m. hyopneumoniae strains  <dig>  and j , m. flocculare strain  <dig>  and m. hyorhinis strain atcc <dig> .
in complex medium, we calculated the ratio between the peak signal in cultivated versus control medium and error bars were calculated as the standard deviation between triplicate time-matched samples.
for defined medium, we detected the actual concentration for the metabolites and error bars were calculated as the standard deviation between duplicate time-matched samples.
a in complex medium, m. hyopneumoniae  and m. flocculare can produce high amounts of acetate; the yields are even higher from m. hyopneumoniae.
m. hyorhinis, on the other hand, produces low concentrations of acetate in this medium.
the final glycolysis product for m. hyorhinis is thus pyruvate.
b in defined medium, m. hyopneumoniae  and m. flocculare produce similar amounts of acetate while m. hyorhinis contains only residual levels of this metabolite.
the three species can produce low amounts of formate in both media
to check if the lower concentrations of acetate were a result of an impaired activity of the pyruvate dehydrogenase complex in m. hyorhinis, we compared the predicted proteins pdha, pdhb, pdhc and pdhd from all species.
the known active sites of odp <dig>  and oldh <dig>  from bacillus subtilis were analyzed in search for a possible mutation.
no active sites have been described for pdha and pdhb so far, but the enzymes for all three swine mycoplasmas do not present significant differences.
although not all active sites from pdhd seem to exist in mycoplasmas, all predicted proteins from the species analyzed  are similar in terms of active sites, when these are present.
the only difference between m. hyorhinis and the other two species was found in pdhc, in the protein portion related to the binding to pdha and pdhd.
this portion is the most variable between the species.
an enzyme characterization should be made in order to verify if this complex is not active in this species, or if it is responsible for another unforeseen activity.
the presence of myo-inositol in swine serum  <cit>  might be directly related to the higher production of acetate in the growth of m. hyopneumoniae in complex medium, if compared to m. flocculare.
since the myo-inositol pathway produces accoa, the recovery of a molecule of coa from myo-inositol is possible and may be directly linked to the production of acetate.
since the myo-inositol catabolic pathway is one of the few distinctions between the metabolic models of m. hyopneumoniae and m. flocculare, it is possible to assume that the differences in acetate concentration in vitro may indeed arise from the ability of m. hyopneumoniae to uptake myo-inositol.
indeed, when no source of myo-inositol is present , no difference in the concentration of acetate is observed for these two species.
the analysis of the amino acid uptake was not trivial since both media contained peptone, which at first was not measured through the cpmg experiment but, in the course of time, was degraded into single amino acids, and thus changed the overall signal of the nmr spectra.
to address this issue, we only took into account the yus+ medium at a single time point after  <dig> h of growth.
m. hyorhinis seemed to have lower concentrations of most amino acids at the end of the growth curve if compared either to the control medium or to the other species.
this may be related to several factors, such as lower rates  of peptone degradation by membrane proteases or higher uptake rates of amino acids by m. hyorhinis.
at this point, we could no longer verify which hypothesis is more suitable for this particular distinction, however, since peptone was not necessary for the maintenance of the viability of m. hyorhinis in the defined medium, it is possible that this species may harbor transporters with higher specificity or effectiveness for single amino acids.
the results as concerns the other metabolites indicate for the three species the uptake of glucose and nucleotides.
flux balance analysis
fba was used to check the properties and capabilities to produce detected metabolites by nmr spectra of all models .
fba simulations were performed with an in silico medium that contained all biomass precursors.
similarly to the in vitro growth results, in the presence of oxygen, acetate was the main product and nadh was recycled by the conversion of molecular oxygen into water.
in the absence of oxygen , all models were able to grow but the main product became lactate.
as expected, the availability of myo-inositol in the in silico medium allowed an alternative source of coa and higher concentrations of acetate as a final product in the models for m. hyopneumoniae.
hydrogen peroxide production was linked to the metabolism of glycerol in both m. hyopneumoniae and m. hyorhinis, while the models for m. flocculare did not produce toxic levels of this metabolite.
since we have no information on gene essentiality for any species, we checked reaction essentiality in the pan-network.
this means that we deleted individual reactions  even if they did not have a gpr association.
a total of  <dig> reactions  were essential for growth in the pan-network .
discussion
in this work, we created a metabolism as realistic as possible for the three known mycoplasmas present in the respiratory tract of swine.
it is essential to point out that, although we had to include  <dig> % of orphan reactions in order to allow growth, we only used about half of the genomes of these organisms.
the other half consists either of hypothetical or of conserved hypothetical proteins  <cit> , which may in part fill the missing gaps of the models.
the number of essential reactions in the pan-network in this work is not directly comparable to the number of indispensable genes predicted for m. genitalium  or m. pneumoniae   <cit> .
gene essentiality accounts for more than metabolic enzymes which means that proteins not included in our models, related to protein synthesis, dna polymerization or rna turnover, were not accounted for in the essential reactions of our models.
if we added these proteins, we would possibly arrive at numbers closer to those of m. genitalium and m. pneumoniae.
the lack of experimental information on gene essentiality is also a setback for the validation of the models created here.
further experiments on this matter should help us better refine and complete the networks.
as for most of the mycoplasma species studied , all reconstructed networks exhibit low connectivity due to the simplicity of the biological model.
we were able to show in this work that the three swine mycoplasma species have similar metabolic capabilities, except for the metabolism of myo-inositol, amino sugar, carbohydrates and glycerol .
overall, the methods used here enabled us to address some of the main problems caused by most automatic reconstruction methods which are the permissive inclusion of pathways and over prediction of capabilities.
annotation errors arise with the attribution of ambiguous or partial ec numbers  <cit>  and the propagation of these errors may then lead to many other ones.
thus, it is known that the prediction of gpr associations based only on name-matching and ec-codes is not sufficient to add confidence to a model  <cit> .
when we decided to simultaneously refine the  <dig> models, we enhanced the confidence of each gpr association, by adding information on the synteny between the genomes, protein sequence alignments and phylogenetic distance between orthologs.
fig.
<dig> differential metabolism of m. hyorhinis, m. hyopneumoniae and m. flocculare.
metabolites are depicted in dark green and separate enzymatic activities for m. hyorhinis, m. hyopneumoniae and m. flocculare can be seen in yellow, pink and blue, respectively.
whenever an enzyme is missing from the three species, the enzyme rectangle is depicted in grey.
complete list of metabolite abbreviations and ec numbers can be found in additional file 4
our main objective was not only to reconstruct the metabolic models for m. hyorhinis, m. hyopneumoniae and m. flocculare.
we also wanted to compare the metabolism of the three species and find possible links to virulence, host-colonization capacity and life-style.
m. hyorhinis, for instance, is not only a pathogen, but also a common contaminant in many mammalian cell cultures  <cit> .
this may be explained by the wide range of carbohydrates which this organism can uptake and metabolize.
this may be a crucial factor to the ability of m. hyorhinis to grow in diverse sites inside the host  <cit>  and even to invade other hosts and potentially develop cancer in humans  <cit> .
m. hyopneumoniae has been detected in other sites such as brain, liver and spleen  <cit> , but m. flocculare has never been detected outside of the respiratory tract of swine.
one may wonder what could explain this if both organisms are so similar.
m. hyopneumoniae is up to date the only mycoplasma species with sequenced genome that has the genes for the catabolism of myo-inositol.
myo-inositol is readily abundant in the bloodstream of mammalian hosts, and can be used as a secondary carbon source for energy production  <cit> .
mycoplasma iguanae has been described to produce acid from inositol  <cit> .
however, there is no available complete genome sequence for this organism, making impossible any comparison of the genes involved in this pathway.
although at this point we cannot confirm that this pathway is functional in m. hyopneumoniae, the nmr results we obtained suggest that myo-inositol might be directly related to the higher production of acetate in complex medium of m. hyopneumoniae if compared to m. flocculare.
this assumption is based on two factors: first, a previous work has detected myo-inositol in swine serum, which is a component of the complex medium  <cit> ; and second, the only difference in the networks that could influence the acetate concentration is the myo-inositol catabolism.
this pathway may also be involved with an alternative production of accoa, precursor to coa, an essential cofactor for growth in all species.
these metabolic distinctions may help m. hyopneumoniae to grow in diverse sites if compared to m. flocculare.
this is particularly interesting when we take a close look into the composition of lung surfactant from several mammalian models  <cit> , which likely contains myo-inositol as a degradation product from phosphatidylinositol.
thus, it is possible that the myo-inositol catabolism in m. hyopneumoniae is one of the reasons for the high virulence of this species when compared to m. flocculare and m. hyorhinis.
a striking result from the nmr results was the reduced capability of m. hyorhinis to produce acetate.
this was not foreseen in the model reconstruction, but predicted protein sequences from the pyruvate dehydrogenase complex of the three species showed a particular distinction in pdhc, more specifically in the region responsible for binding to other complex components.
we could not verify at this point if this complex is not active in the particular growth conditions of our work in m. hyorhinis or if it is responsible for another unforeseen activity.
however, a transcriptome profiling of m. hyorhinis has detected both pdha and pdhb in the pool of genes with the highest number of transcript reads  <cit> , indicating that the complex might be translated in vivo.
moreover, pyruvate dehydrogenase activity has been previously detected in m. hyorhinis extracts  <cit> , and it seems to be correlated to oxygen availability.
in our cultivations, the cells are not grown in a complete aerobic system, which may explain the differences between our findings and those previously published.
cultivation tests in vitro also showed that none of the species was able to grow in none of the defined media.
this is probably due to a lack of the correct cofactors as previously mentioned.
for instance, thiamine, pyridoxal/pyridoxine, pantothenate, spermine and folate were the actual media components tested, but the reconstructed models were not able to convert them to biomass precursors .
only m. hyorhinis seems to be able to uptake folate directly while the other two species might need intermediate metabolites.
this might also explain why m. hyorhinis is extensively found as a cell culture contaminant.
even though we supplemented cmrl- <dig>  with all missing components present in the yus+ medium, cells maintained viability for longer periods in the latter.
cofactors present originally in both media existed in lower concentrations in the cmrl- <dig>  this means that possibly one or more cofactors with higher concentrations in the yus+ medium are essential for the viability of the three species.
common cofactors with higher concentration in yus+ medium included: choline, folate, pantothenate, pyridoxal, thiamine and spermine.
other differences among the species seem to indicate that m. hyopneumoniae and m. flocculare lack one or more amino acid transporters.
this may be due to two things:  both remained viable for longer periods whenever cmrl+ was supplemented with peptone, and  we detected a higher amino acid accumulation in the defined media if compared to m. hyorhinis.
m. hyorhinis on the other hand may have all the transporters for single amino acids, but this hypothesis needs to be further verified.
from sequence analyses, we did not find significant differences among the species that would explain this distinct behavior.
moreover, since m. hyorhinis seemed to overgrow in the first  <dig> h of culture in both cmrl media , it is possible that one or more compounds were missing in the yus+ medium in order to allow growth of this species.
possible candidates for such  would be: coa, nicotinamide, ascorbate and/or pyridoxine.
it is still not clear if any of the swine mycoplasmas studied here indeed have a polysaccharide capsule, given that only one study has assesed its presence in m. hyopneumoniae.
however, this is particularly important because in this study tajima and yagihashi  <cit> , as discussed before, reported that capsular polysaccharides from m. hyopneumoniae play a key role in the interaction between pathogen and host.
moreover, it has been reported that some strains of m. hyopneumoniae become less pathogenic in broth culture and, after serial passages, lose their ability to produce gross pneumonia in pigs  <cit> .
in several bacterial species, it has been established that the amount of capsular polysaccharide is a major factor in virulence  <cit>  and it decreases significantly with in vitro passages  <cit> .
furthermore, despite increasing evidence supporting the existence of a polysaccharide capsule in several mycoplasma species, it is not yet known if the material is synthesized by the bacteria or imported from the medium for most species  <cit> .
recent reports have even shown the capacity of species from the mycoides cluster to assemble a capsule from exogenous phosphorylated glucose  <cit> .
however, no homolog proteins related to these activities were found in none of the species studied in this work.
for the reasons discussed above, we were not able to introduce the capsule production itself in the models; however, we tried to relate its possible existence with the substrates that might be used to produce it.
the composition of capsular polysaccharides has been related to the level of pathogenicity in other bacteria.
for instance, highly virulent escherichia coli strains seem to have practically nonimmunogenic capsular material, due to the fact that these antigens are similar or identical to the ones found in the host  <cit> .
the structure of the k <dig> antigen from e. coli capsules is basically glcnac and glucuronic acid  in a molar ratio of 1: <dig>  <cit> .
erlinger and collaborators  <cit>  identified heparan sulfate  as the predominant glycosaminoglycan in the porcine respiratory tract and the most common disaccharide unit within hs is glcnac linked to glca.
since m. hyopneumoniae and m. flocculare might be able to import glcnac unchanged, they may use it for the composition of its capsule; m. hyorhinis, on the other hand, imports sialic acid, and might use it directly to produce its own.
of course this is highly speculative, and we are not sure if all these hypotheses indeed happen in vivo.
the metabolism of glycerol is the final major difference between the pathogenic m. hyopneumoniae and m. hyorhinis and the non-pathogenic m. flocculare.
such metabolism and the production of hydrogen peroxide by the action of glpo are essential for the cytotoxicity of m. pneumoniae  <cit>  and m. mycoides subsp.
mycoides  <cit> .
particularly, the cytotoxicity of m. mycoides subsp.
mycoides is due to the translocation of the hydrogen peroxide into the host cells.
this is only possible because of the close proximity to the host cells along with the membrane-bound enzyme.
highly conserved homolog genes to this enzyme were only found in the genomes of the pathogenic species studied here.
whether or not these enzymes from m. hyopneumoniae and m. hyorhinis are indeed capable of the same activities as glpo is not yet confirmed; however, the high similarity may be an indication that this trait is essential for the pathogenicity of these two species.
this may also explain why both m. hyopneumoniae and m. flocculare can adhere to the cilia of tracheal epithelial cells, but only the adhesion of m. hyopneumoniae results in tissue damage  <cit> .
moreover, while m. hyorhinis has in the reconstructed models three ways of uptaking glycerol, m. hyopneumoniae seems to have five, and this may reflect in its enhanced pathogenicity.
CONCLUSIONS
we presented in this work an overview of the differential metabolism of m. hyopneumoniae, m. hyorhinis and m. flocculare using different approaches.
the reconstructed models showed some distinctions among the species, namely the myo-inositol metabolism for m. hyopneumoniae, the wider uptake of carbohydrates for m. hyorhinis and the usage of glycerol as a carbon source for the two pathogenic species.
the models also served as a basis for all the assumptions made for the experimental data.
metabolic profiling of both complex and defined media pointed to new differences that we were not able to identify based solely on the sequenced genomes.
the major ones were related to the pyruvate conversion to acetate, which appeared to be higher in m. hyopneumoniae and m. flocculare than in m. hyorhinis.
after we performed growth tests using defined media that hypothetically would lead to bacterial growth in silico, we assessed that some cofactors or metabolites were probably being delivered in the wrong form to the mycoplasmas.
the environmental context may explain the differences between the in silico models and the in vivo behavior.
whether the main differences among the species we reported here  are related to virulence or pathogenicity have not yet been addressed experimentally, but it is tempting to speculate.
the same factors that may enhance virulence of m. hyopneumoniae may help the commensal species m. flocculare to better survive inside the host.
all these in silico and in vivo metabolic differences among m. hyopneumoniae, m. hyorhinis and m. flocculare might influence the different levels of pathogenicity in each of them.
however, it is highly likely that gene regulation may directly interfere in the metabolism and pathogenicity and may be related to many other aspects still unaccounted for.
one of our future goals is therefore to understand and integrate gene regulation into the metabolic models.
upcoming experiments will aim at testing the hypotheses formulated here, particularly those related to the metabolisms of glycerol and myo-inositol.
we also intend to better understand the habitat of these species, and the possible metabolic and genetic dialogues with the host and other bacteria present in this environment.
either way, this work serves as a basis for the study of the differential metabolism and pathologies caused by the swine lung mycoplasmas and may help to propose ways to prevent disease development in the future.
availability of data and material
all models from this work are available in additional file  <dig> and online at the metexplore database .
additional files
additional file  <dig> supplementary tables.
table s <dig>  defined media composition used for the cultivation of m. hyopneumoniae, m. hyorhinis and m. flocculare.
table s <dig>  complete list of reactions and metabolites from all models.
table s <dig>  list of all metabolites included in the final biomass composition.
table s <dig>  distribution of gpr associations and exclusive gpr associations by metabolic pathways.
table s <dig>  proposition of correct precursors and cofactors for a new defined medium.
table s <dig>  nmr metabolomics: metabolites detected in the complex friis and yus+ media for the three swine mycoplasmas.
table s <dig>  prediction of essential reactions in the pan-network.
table s <dig>  differential metabolism overview.
additional file  <dig> biomass reaction assembly from reconstructed models.
pdf file containing all steps necessary to the assembly of the biomass reaction for all species.
additional file  <dig> abbreviation of metabolites and ec numbers from figs.
<dig>   <dig> and  <dig>  spreadsheet containing the lists of all abbreviations and ec numbers used in figs.
<dig> and  <dig> separately.
additional file  <dig> prediction of transporters.
pdf file containing all assumptions for the in silico prediction of transporters.
additional file  <dig> models in sbml format.
compressed file  containing all sbml models from this work .
abbreviations
atp-binding cassette; acpacyl-carrier protein
blastpblast local alignment search tool
demdead-end metabolite
ccucolor changing unit
cpmgcarr-purcell-meiboom-gill
ec numberenzyme-code number
fbaflux balance analysis
fidfree induction decay
gamgrowth associated maintenance
gogene ontology
gprgene-protein-reaction association
hcahierarchical clustering analysis
hsqc spectroscopyheteronuclear single-quantum correlation spectroscopy
mifmajor facilitator protein
mwmolecular weight
ngamnon-growth associated maintenace
nmrnuclear magnetic resonance
noesynuclear overhauser effect spectroscopy
pcaprincipal component analysis
pfam databaseprotein-families database
ptsphospho-transferase transport system
sbmlsystems biology markup language
tmsptrimethilsilylpropionate
tocsytotal-quantum correlation spectroscopy
tutranscriptional unit.
enzyme abbreviations: abc-gly abc-glycerol transport system
glyfglycerol facilitator protein
glpoglycerol-3-phosphate
glpuglycerophosphodiester transporter
fbafructose-biphosphate-aldolase
iolj5-dehydro-2-deoxy-phosphogluconate aldolase
ldhlactate dehydrogenase
mdhmalate dehydrogenase
pdh a/b/c/dpyruvate dehydrogenase complex subunit a/b/c/d
ulagl-ascorbate-6-phosphate lactonase
metabolite abbreviations used in the text : 2kg2-oxoglutarate
accoaacetyl-coenzyme-a
coacoenzyme-a
dhadihydroxyacetone
dhapdihydroxyacetone-phosphate
f6pfructose-6-phosphate
glcnacn-acetyl-glucosamine
nad/nadhnicotinamide adenine dinucleotide
nadp/nadphnicotinamide adenine dinucleotide phosphate
methf <dig> -methylene-tetrahydrofolate
succoasuccinyl-coenzyme-a
competing interests
the authors declare that they have no competing interests.
authors’ contributions
mfs and az conceived and designed the work in close collaboration with all the other authors.
mgf performed all computational metabolic reconstructions and most of the experimental work.
fms, sgm, tlp and ej collaborated in the experimental work.
tlp performed the metabolomics statistic analysis.
beh, ft, atrv and iss collaborated analyzing all data.
mgf, mfs and az wrote the manuscript with inputs from all the other authors.
all authors read and approved the final manuscript.

 the respiratory tract of swine is colonized by several bacteria, among which are three common mycoplasma species: mycoplasma flocculare, mycoplasma hyopneumoniae, and mycoplasma hyorhinis. even though little information is available concerning the prevalence of bacteria in healthy lungs, these three mycoplasma species have been isolated from the respiratory tract of both healthy and diseased pigs. while m. flocculare is usually described as a commensal bacterium, m. hyopneumoniae and m. hyorhinis are considered pathogenic. enzootic pneumonia, caused by m. hyopneumoniae, is widespread in pig populations and is considered a major cause of economic loss in the pig industry; m. hyorhinis is frequently present in cases of polyserositis and arthritis and m. flocculare has high prevalence in swine herds worldwide, but up to date, no disease has been per se associated with this species. in addition to these mycoplasmas, mycoplasma hyosynoviae, the primary agent of non-purulent arthritis, can occasionally colonize in the lower respiratory tract when pneumonic lesions pre-exist. besides mycoplasmal pneumonia, the porcine respiratory disease complex has emerged as an economically signigicant respiratory disorder characterized by the slow growth, fever, cough, loss of appetite, lethargy and dyspnea in pigs. even though many species are related to prdc, it is essential to note that enzootic pneumonia caused by m. hyopneumoniae is by far the most costly disease in pig industry, and this bacteria is usually seen as an essential component to the successful establishment of a pathogenic community in the host. also, m. hyopneumoniae infections take longer to cause lesions and to be successfully eliminated than infections from other pathogens. while mycoplasmal diseases in swine have been extensively studied, their causative agents have not been explored from a mathematical and computational point of view, mostly because their genome sequences were not available until recently. although recent studies have placed m. hyopneumoniae, m. flocculare and m. hyorhinis in close proximity within the hyopneumoniae clade by phylogenomic analysis, which corroborates with their high 16s rrna sequence similarity, it is not yet clear what causes the specific pathogenicity or lack thereof in each of them. this elevated genomic resemblance combined with their different levels of pathogenicity is an indication that these species, as for most mycoplasmas, have unknown mechanisms of virulence and differential expression. pathogenic determinants such as adhesion to the host cell and evasion from the immune response have already been well-described in the literature for both m. hyopneumoniae and m. hyorhinis. the presence of a capsule in m. hyopneumoniae has been reported to be important for the interaction with the host cells in one single study. several studies also show that immunosuppressed animals experimentally infected with mycoplasma species develop less severe microscopic lesions of pneumonia if compared to normal animals. this means that it is possible in some cases that a strong response from the host immune system might be the primary cause of pathogenesis. however, up to date, it is not possible to draw any further conclusions due to lack of experimental data. even if these topics are of utter importance for understanding swine respiratory tract mycoplasmal diseases, what has yet to be better understood is the direct participation of metabolism in the development of them. for instance, although adhesion factors are related to pathogenicity, m. hyopneumoniae and m. flocculare harbor similar sets of adhesion proteins, and have been shown to adhere to cilia in a similar way. thus, the ability of m. hyopneumoniae to cause disease if compared to m. flocculare might not be directly related to adhesion. furthermore, the genome sizes of mycoplasma spp. range from kb to more than kb, representing an important example of genome reduction during the evolutionary process. it is possible that in an initial symbiotic phase, the host provided a broad range of metabolites for these bacteria. this, together with the ability of the bacteria to uptake such compounds, made several activities dispensable for the bacterial life. over the course of evolution, these bacteria would have lost some of the genes that became unnecessary for life in an environment conditioned by another genome. in this way, we aimed at studying the reduced metabolism through genome-scale metabolic model reconstructions of m. hyopneumoniae, m. hyorhinis and m. flocculare to better understand their different life-styles. based on the reconstructed networks, we propose that one of the mechanisms that may explain why m. hyorhinis and m. hyopneumoniae are pathogenic while m. flocculare is not, is their ability to use glycerol as a carbon source, thus enabling the production of the highly toxic hydrogen peroxide. this trait may be directly involved in cytotoxicity, as already described for two lung pathogenic mycoplasmas, namely mycoplasma pneumoniae in human and mycoplasma mycoides subsp. mycoides in ruminants. additionally, growth rates and fitness have been correlated in the past to virulence in several other organisms and may be a key factor for the difference in pathogenicity among these three mycoplasmas. indeed, the in silico models show that m. hyorhinis and m. hyopneumoniae have extra sets of enzymes that may help them overgrow the non-pathogenic species m. flocculare. m. hyorhinis seems to have a wider carbohydrate uptake range, and probably for this reason is seen as a well-known contaminant of cell cultures. m. hyopneumoniae, on the other hand, might use myo-inositol as a carbon source producing acetyl-coenzyme-a as a byproduct. this feature is of particular interest for the production of cofactor coenzyme-a in this species, since most enzymes from the biosynthetic pathway are missing in all species. nevertheless, there is little metabolic experimental data available for the three species, which makes the reconstruction of a reliable metabolic model an extremely time-consuming work. together with the fact that these mycoplasmas are generally grown in complex media, with high serum concentrations, we needed additional experimental data in order to compare the in silico reconstructed metabolic networks with the in vivo metabolic characteristics of the three species. for this reason, we also performed nuclear magnetic resonance spectroscopy analyses to detect metabolites consumed and produced in both complex and defined media. these experiments corroborated with the reconstructed models and suggested two new features in particular: the uptake of myo-inositol in m. hyopneumoniae might be related to a higher acetate production, and m. hyorhinis showed a surprisingly reduced ability to convert pyruvate to acetate in the growth conditions used in this study. all these in silico and in vivo metabolic differences might influence the different levels of pathogenicity in each of the species studied here. we present in this work the metabolic networks reconstructed from the annotated genomes of the three species and the comparison done with growth rates and metabolomic experiments performed in vitro in order to better understand the basis of the pathologies caused by these bacteria, which might help prevent their development in the future. methods modeling methods network reconstruction and refinement we reconstructed the metabolic networks for all available strains isolated from swine so far of m. hyopneumoniae, m. hyorhinis and m. flocculare. the semi-automated reconstructions were generated by the pathologic tool from the pathway tools software using the complete genomes available online. from now on, the species strains will be abbreviated according to table m. hyopneumoniae m. hyorhinis m. flocculare notes: na: not available ⋆a: strains available only in brazil are protected strains which are property of the government and cannot be used outside authorised laboratories. strains available in france were purchased from the atcc repository; these strains were also available for testing in the laboratory in brazil ⋆b: m. hyorhinis strain gdl- was retrieved from a contaminated cell line and, to our knowledge it has never been reported neither as pathogenic nor as an attenuated strain ⋆c: m. hyorhinis strain atcc genome is available online but not published in any paper up to now. pathway tools automatically associated genes with reactions, based on the annotation names, gene ontology terms and enzyme code numbers contained in the genbank files. the software assembled reactions into pathways by comparing them with the reference database, metacyc. the software automatically added missing reactions from the reference database to the model, creating many orphan reactions. in a gapfilling step, we used the pathway hole filler tool, included in the pathway tools software, to identify possible candidate genes associated with these orphan reactions. orphan reactions in the metabolic network were allowed only to allow known functionalities of the organisms or biomass production in the reconstructed models. we incorporated the minimal number of orphan reactions into the models to avoid adding a behavior that was never described. in our work, the refinement of the networks was made in a subsystems approach simultaneously on all organisms. as a result of different genomic annotations, slightly different reactions arose throughout the automatic models. these reactions were checked for consistency, and after validation of only one of them, the duplicates were deleted. gene-protein-reaction associations were systematically validated or included based on experimental data, information provided by the pfam functional domain database, synteny analysis and reciprocal sequence homology searches, using at first high e-value cutoffs of at least 10−. whenever no hits were achieved with these parameters, more relaxed parameters were assumed. homolog proteins with equivalent functional domains were assigned as isozymes to a particular reaction, while proteins with distinct functional domains were assigned as subunits of a multi-protein complex. non-metabolic reactions such as dna polymerization, protein synthesis and rna synthesis were explicitly deleted from the network but implicitly included in the biomass assembly. generic reactions were either specified or excluded from the models. the resulting models were further refined in accordance with a detailed protocol from palsson and thiele. reactions were computationally balanced for mass and charge, while cofactor usage was determined based on literature data for closely related species. reaction directionality was thermodynamically checked and validated. reaction directionality was determined based on the component contribution method, which extends the group contribution method and achieves a significant improvement in the accuracy of the estimations of standard gibbs energies. an online search and calculation interface called equilibrator along with metabolite and reaction thermodynamic databases are available at. heuristic rules were used to improve the directionality assignment. we analyzed the network topology to identify compounds that were only produced or consumed in the network, the so-called dead-end metabolites. whenever a dem was found, we either added or validated an orphan reaction to reconnect it to the rest of the network, or removed the reaction from the model when both substrate and product were disconnected from the network and/or did not affect the overall metabolism. transporters were predicted at first from genome annotation using the transport identification parser, from pathway tools. transporters were predicted mainly based on sequence similarity to other known transporters, since experimental data for the three mycoplasmas were not available. specific transport reactions and exchange reactions required for production of biomass components were manually added to the final version of the networks. biomass composition and biomass equation in order to simulate growth, we had first to estimate the average cell composition of these mycoplasmas. to this end, we used a general macromolecule mycoplasma cell composition from razin and collaborators. depending on the species, the authors showed that the solid residue contained 54– % of proteins, 12– % of lipids, 3– % of carbohydrates, 8– % of rna, and 4– % of dna. the membranes comprised around % of dry weight of the organisms and contained 47– % of proteins, 35– % lipids, 4– % carbohydrates and small amounts of dna and rna. since no information on metal ions and cofactors was available, we included them quantitatively based on the metabolic networks reconstructed for related mycoplasmas. membrane and lipid components were added based on literature composition of the selected species. it is important to note that the only difference in composition between the species reported so far was related to the presence of glycolipids in m. hyopneumoniae and m. flocculare, which were reported to be absent in m. hyorhinis. based on the previous information, we assumed the following fractions of macromolecules: % proteins, % lipids, % carbohydrates, % rna, % dna, and % of ions and cofactors. to create a biomass elemental formula, we took into account the percentage contribution of each of the components to the overall cell. the biomass reaction represented the drain of these components into biomass production. it implicitly assembled dna replication, rna transcription, and protein synthesis into one single reaction. amino acids were indirectly included in the biomass reaction: charged trnas were accounted as substrates and uncharged trnas as products. growth and non-growth associated maintenance were estimated based on the literature and were manually added to the models. the details for the complete assembly of the biomass reaction can be found in additional file model validation and fba analysis the metabolic networks were exported as a mathematical systems biology markup language models. they were uploaded to the cobra-toolbox v matlab extension for flux balance analysis testing. minimum-maximum flux constraints were imposed based on literature information. growth simulations were achieved using biomass production as the fba objective function. experimental methods swine mycoplasmas cultivation complex media comprised friis media, and a commercial mycoplasma broth. defined media cultivation tests were performed in france and comprised a medium described for m. pneumoniae strain by yus and collaborators, and commercial media cmrl with no glutathione. since we had no information on metabolism of swine mycoplasmas, we decided to supplement the defined yus medium with all amino acids. we also supplemented the cmrl- medium with other peptone and/or other cofactors. information on the composition of all defined media is available in additional file 1: table s cells were cultivated at °c for different time periods, under gentle agitation. cell concentration estimation and viability by color changing units measurement cell growth and viability was measured with triplicate time-matched samples of cells and culture media for ccu as described by stemke and robertson. viability of cells is visible by a change in medium color from red to yellow. for cell concentration measurements, the cultures were subjected to a series of 10-fold dilutions in complex media and ccu/ml was defined as the highest dilution of cells able to change the medium color. samples for nmr spectroscopy nmr was performed with complex friis medium and defined yus+ medium. the medium was collected at the following time intervals: h, h, h, h, h, and h for friis medium and h, h, h, h, h, h and h for yus+ medium. cells were separated from growth media through sedimentation at g for nmr analysis. samples consisted of biological triplicates in complex medium and biological duplicates in defined medium. since we had slight different time intervals, care was taken by not directly comparing the results from different experiments. nmr analysis sample preparation was as follows: μl of a mixture containing mkh2po phosphate buffer in d2o with mmnan and % trimethylsilyl propionate was added to μl supernatant samples. both solutions were mixed thoroughly and μl were then transferred to mm nmr tubes and sorted in 96-tubes racks. all nmr experiments were carried out on a bruker mhz nmr spectrometer equipped with a mm txi probe and a samplejet autosampler, enabling high-throughput data acquisition for large collections of samples. the temperature was controlled at °c throughout the experiments, and the samples were kept refrigerated at °c during a waiting time of less than h in the autosampler, before the nmr analysis. standard 1h 1d nmr pulse sequence nuclear overhauser effect spectroscopy with z-gradient and carr-purcell-meiboom-gill with water presaturation were applied on each sample to obtain the corresponding metabolic profiles. a total of transient free induction decays were collected for each experiment with a spectral width of ppm. the relaxation delay was set to s. the noesy mixing time was set to ms. the total acquisition time of each sample was min s. the cpmg spin-echo delay was set to ms, for a total filter of ms, allowing an efficient attenuation of the lipid and protein nmr signals. the 90° pulse length was automatically calibrated for each sample at around μs. data processing: all fids were multiplied by an exponential function corresponding to a hz line-broadening factor, prior fourier transformation. 1h-nmr spectra were manually phased and referenced to the tsp signal using topspin. extraction of a data matrix for multivariate statistical analysis from the 1h nmr profiles was done using the statistics toolbox of amix. spectra were integrated from to ppm at a step of ppm but excluding the regions of residual water at - ppm. no normalization of the intensity was performed. the resulting data matrix contains nmr variables. multivariate data analysis: principal component analysis and hierarchical clustering analysis were performed using simca-p with scaling based on the pareto method. metabolites identification and quantification: metabolite identification was achieved by comparing spectra with databases such as hmdb. identification of the metabolites was further verified with homonuclear and heteronuclear 2d nmr experiments such as 1h−13c hsqc, 1h−1h tocsy and j-resolved experiments. absolute quantification of the metabolites was performed using chenomx nmr suite.  model reconstruction and refinement based on the published genomes of strains of m. hyopneumoniae, strains of m. hyorhinis and strains of m. flocculare, we reconstructed genome-scale metabolic models: one model for each strain separately, one for each species and a pan-reconstruction for all three species. the semi-automated reconstructions were manually refined according to the description given in the “methods” section. the pan-swine mycoplasma network was initially composed by reactions. duplicate reactions arose from the fact that pathway tools is based only on the annotations contained in the genbank files. we removed generic and duplicate reactions from the models and replaced them with the specific and validated ones. non-metabolic reactions were also excluded at this point, along with absent capabilities of these mycoplasmas, such as heme, quinone or cytochrome dependent reactions. the networks were also tested for the presence of dem and the ability to produce all biomass precursors. dems were analyzed on a case-by-case manner. from initially dems, disconnected metabolites were excluded from the models: dems came from spontaneous reactions; were carbohydrate substrates derived from wide range transport reactions and were not used by any other reaction in the network; were excluded as they did not interfere with the overall metabolism and/or the enzyme had already been assigned to one or many other reactions. orphan reactions were excluded unless they were essential for biomass growth. all excluded reactions and all dems are listed in additional file 1: tables s2c and s2d. it is interesting to point out that other unknown enzymes from these organisms may indeed use some of the excluded dems; however, since we have no experimental evidence at present, we could not assess their interference in the metabolic models. the remaining dems consisted mainly of cofactors and biomass precursors disconnected from the rest of the network. they were solved along with the biomass precursor check; reactions were added to allow growth: biomass reaction, enzymatic reactions, drain-synthetic reaction, transport reactions and spontaneous reactions. transport reactions were considered as such even if an enzymatic activity was also present. in order to correctly assign possible and specific transporters, we performed an extensive literature search and reciprocal blast alignments to characterized transporters in other species. these results have to be experimentally confirmed, since assignment of the correct substrate based only on sequence homology remains an open problem in genomic annotation. we also changed inconsistencies of reversibility and cofactor usage. even after all efforts, four dems still remained: ttp, hexulose-6-phosphate, deoxyinosine and xanthosine-5-phosphate. after the final addition of exchange reactions, we ran fba tests to check the consistency of all models. if a reaction was essential for biomass growth and no homolog gene was found in the genome, an orphan reaction was validated or added to the reconstruction for modeling reasons only. the resulting refined pan-network had a total of reactions and gpr associations. a comparison between the species models may be found in fig. these results indicate that all strains from all species are indeed metabolically similar. fig. venn diagrams representing the comparison of refined networks. numbers represent the exclusive and common reactions present in the refined networks between species, and between selected strains of m. hyopneumoniae. this analysis shows that most of the metabolism is common to all organisms. mhr: m. hyorhinis; mhp: m. hyopneumoniae; mfl: m. flocculare the overall characteristics of each reconstruction can be seen in table a list of all reactions and metabolites added to the models along with the corresponding genes can be seen in additional file 1: table s notes: ∗a: gprs are gene-proteins-reaction associations present in each genome. for imfmhp, imfmhr, imfmfl and imfpan that account for more than one species, we added a gpr to the species model when any of the strains harbored a gene responsible for a specific activity ∗b: reactions were considered as transport reactions even if the transporter was capable of performing a concomitant enzymatic activity ∗c: spontaneous reactions included diffusion of small molecules and spontaneous conversions ∗d: total number of reactions included also all exchange reactions needed for the mathematical modeling ∗e: genes in the pan-network and species models were added synthetically, to enable a reaction essentiality analysis biomass composition and biomass equation the biomass equation drained all precursors into biomass. biomass composition according to the “methods” section along with a detailed description of the assembly of the biomass reaction can be found in additional file 1: table s and additional file considering the percentage contribution of each of the components to the overall cell, an approximate biomass elementary composition for m. hyorhinis was computed as follows: ch o n p s, with traces of calcium, chlorum, cobalt, copper, iron, potassium, magnesium, manganese, molibdenium and zync and molecular weight on a c-mole basis of g/ c−mol. we had a slightly different composition for m. hyopneumoniae and m. flocculare: ch o n p s, with the previously mentioned trace compounds and mw on a c-mole basis of g/ c−mol. metabolism of m. hyorhinis, m. hyopneumoniae and m. flocculare similar to most mycoplasma species studied so far, all reconstructed networks exhibit low connectivity due to the simplicity of the biological model. out of the reactions in the final pan-network, had in at least one species a gpr association. from these, were common to all species. the overall metabolism from the models reconstructed consisted of distinct subsystems: amino sugar metabolism, amino acid metabolism, carbohydrate metabolism, cofactor metabolism, lipid metabolism, and nucleotide metabolism. while all enzymes were present in glycolysis, most metabolic pathways had major enzyme gaps. in the product and cofactor metabolism, for instance, gaps accounted for up to % of the reactions. m. flocculare was the only species that did not show any exclusive metabolic activities in the models, although this could not be verified in vivo due to lack of biochemical studies for this species. m. hyopneumoniae had exclusive reactions, linked to the myo-inositol metabolism and alcohol dehydrogenase activity. m. hyorhinis had exclusive reactions; most of them corresponded to carbohydrate metabolism. a global model comparing each species enzymatic capabilities can be seen in figs. and the complete list of metabolite abbreviations and ec numbers for these two figures can be found in additional file fig. distribution of the model reactions in the subsystems. the reactions present in the model imfpan were separated into biological subsystems and further into reaction types, with the exclusion of exchange reactions in this analysis fig. global energy and carbohydrate complete models. metabolites are depicted in dark green and separate enzymatic activities for m. hyorhinis, m. hyopneumoniae and m. flocculare can be seen in yellow, pink and blue, respectively. whenever an enzyme is missing from the three species, the enzyme rectangle is depicted in grey. complete list of metabolite abbreviations and ec numbers can be found in additional file 4 fig. lipid, amino acid, nucleotide and cofactor complete models. metabolites are depicted in dark green and separate enzymatic activities for m. hyorhinis, m. hyopneumoniae and m. flocculare can be seen in yellow, pink and blue, respectively. whenever an enzyme is missing from the three species, the enzyme rectangle is depicted in grey. complete list of metabolite abbreviations and ec numbers can be found in additional file 4 transporters transporter assignment was made in silico based on sequence similarity to known transporters. we were able to identify in all species complete putative abc transporters for multidrug/toxin efflux and for import of sugars, oligopeptides, peptides, spermidine/putrescine, phosphonate/phosphate/thiamine, cobalt, manganese/zinc and glycerol. in the search of a possible nucleotide transporter, we came across with an unpublished work from lee in which a system previously annotated as sugar abc transporter was experimentally validated as involved in nucleoside uptake in mycoplasma bovis. these results may be highly speculative, and need further verification, but we were able to find homolog genes in the genomes of all species. m. hyopneumoniae seems to have an extra abc transporter proposed for myo-inositol and m. hyorhinis has a unique predicted abc transporter for maltose/maltodextrins. three complete phospho-transferase transport systems were common to all species. although specificity of pts is difficult to ascertain, we have tried to compare conserved domains from validated pts substrate-specific enzymes ii and candidates from the three species. genes coding enzyme i and phosphocarrier were found in all strains, except in mhr_hub from sequence similarity and gene context we were able to annotate the correct gene. based on these in silico predictions from domain retrieval, we propose the following: one pts non-specific for sugar, one with specificity for fructose and another for mannitol. m. hyopneumoniae and m. flocculare also share two extra complete systems: one possibly for ascorbate and one for n-acetylglucosamine. m. hyorhinis and m. hyopneumoniae seem to have an extra iib component specific for glucose, and along with the other components of the general sugar pts were proposed to form a complete glucose pts. the protein glpu, coded by gene mpn in m. pneumoniae, was recently described to act on the uptake of glycerophosphodiesters. we found homologs of this gene in all strains of m. hyopneumoniae and m. hyorhinis adjacent to the gene coding for the enzyme responsible for the metabolization of these substrates. in mhr_hub this gene was annotated as a pseudo gene, but from gene context, we were able to validate it as glpu. all species also possess several unspecific amino acid permeases, a glycerol facilitator protein, which is thought to be less efficient than the glycerol abc transporter ) and a major facilitator protein with unknown specificity. we could also find common transporters for cobalt and magnesium, chromate, magnesium, potassium, zinc and a cation atpase. m. hyorhinis has two extra transport systems: one for sialic acid and one sodium/phosphate co-transporter. details of the transport search throughout all species is summarized in additional file carbohydrate and energy metabolism glycolysis as previously described in, all genes for the glycolysis exist in all strains of all species. m. hyorhinis may be able to convert dihydroxyacetone directly to dihydroxyacetone-phosphate with the action of dha kinase, but the mechanism of dha uptake is not clear. pentose phosphate in agreement with previously described literature, we found no enzymes for the oxidative branch of the pentose phosphate pathway in the studied species, turning ribulose-5-phosphate into a dem. this problem was overcome by the assumption that the reaction ribulose-phosphate-3-epimerase can be reversible, as reported in lactobacillus casei. out of the reactions present in the final models, we found homologs for in all species. although a gene coding for transaldolase is missing in all species, we included this reaction in the models to prevent both sedoheptulose-7-phosphate and erithrose-4-phosphate from becoming dem. ascorbate metabolism only one enzyme is missing from the typical pathway: l-ascorbate-6-phosphate-lactonase. this enzyme is responsible for the turnover of l-ascorbate-6-phosphate to 3-keto-gulonate-6-phosphate. from the gene context, mhp_7448_ from mhp_ and its homologs from all other organisms may encode the missing enzyme. the gene was annotated as a conserved hypothetical protein with a proposed phosphotriesterase activity. uptake of other carbohydrates while the uptake of fructose, mannitol and mannose is similar in all species and they all can be fed into glycolysis as fructose-6-phosphate and dhap, m. hyorhinis has an extra set of enzymes enabling the conversion of fructose-1-phosphate into fructose- -biphosphate by the action of 1-phosphofructokinase. only in m. hyopneumoniae the conversion of glyceraldehyde into glycerol might be possible through the moonlight activity of alcohol dehydrogenase, although up to date no experimental evidence can sustain that assumption. the presence of a unique transcriptional unit with carbohydrate metabolism-related genes in m. hyorhinis probably enables the uptake and metabolization of isomaltose, maltose, trehalose and sucrose into either glucose-6-phosphate or f6p. this may be related to the fact that m. hyorhinis can overgrow the other two species in several growth media. myo-inositol metabolism a tu for the myo-inositol catabolism is present in all m. hyopneumoniae strains, with the exception of the enzyme 5-dehydro-2-deoxy-phosphogluconate aldolase. the gene encoding this enzyme in other organisms is similar to the fructose-biphosphate aldolase from glycolysis. since there are two copies annotated for fba in m. hyopneumoniae, we proposed one of them as candidate for this activity based on sequence alignments with both genes from bacillus subtilis. inositol can be used as a carbon source and also produces acetyl coenzyme-a, which can be a source of cofactor coenzyme-a. amino sugar metabolism m. hyopneumoniae and m. flocculare models could uptake and convert glcnac to n-acetyl-glucosamine-6-phosphate; m. hyorhinis models imported sialic acid and converted it to glcnac6p in three steps that are unique to m. hyorhinis. at first, the genes coding for these three enzymes were absent from mhr_hub and mhr_gdl however, based on chromosome alignment and synteny with other genes, we were able to find the regions and correctly annotate all genes. pyruvate metabolism out of fifteen enzymatic reactions in this pathway, only three were not found in any organism. conversion of pyruvate to lactate and acetate were possible in all species. from literature data, we knew that m. hyopneumoniae was able to produce oxaloacetate, malate and 2-oxoglutarate. decarboxylation of 2kg to succinyl coenzyme-a has been previously reported in m. hyorhinis extracts. although we did not find any candidate for the activity of aspartate transaminase, gap-filling proposed that this activity could be performed by the pyruvate dehydrogenase complex. moreover, cordwell and collaborators suggested that in mollicutes lactate dehydrogenase could also function as malate dehydrogenase. m. hyopneumoniae is the only organism among the three to have a gene coding for an alcohol dehydrogenase, probably enabling the turnover of accoa to acetaldehyde and ethanol. conversion of pyruvate to formate by pyruvate formate-lyase, malate to fumarate by fumarate hydratase and succoa to succinate by succinyl-coa synthetase were added to the models in a final step to accommodate the results from the metabolomics experiments. lipid metabolism the lipid metabolism in the models of the three species included the uptake of glycerol, choline, glycerophosphodiesters and fatty acids to the production of the biomass precursors cardiolipin, -diacyl-sn-glycerol, phosphatidylglycerol, phosphatidylcholine and galactosyl-diacylglycerols. extracellular cholesterol and sphyngomyelin were incorporated unmodified directly into biomass. while m. hyopneumoniae seems to have four different ways to uptake and metabolize glycerol m. hyorhinis lacks one and m. flocculare lacks two of them. the turnover of glycerol-3-phosphate into dhap by the action of glycerol-3-phosphate oxidase allows the usage of glycerol as the sole carbon source, with the production of highly toxic hydrogen peroxide. this was only possible in the m. hyorhinis and m. hyopneumoniae models. the enzyme responsible for the turnover of glycerol-3-phosphate into acyl-sn-glycerol-3-phosphate was found only in the genomes of m. hyorhinis. although this can be seen as a possible difference between the species, we had to add the reaction to the models from the other species to enable growth. it is possible that the other two species import this metabolite directly from the media, however further experiments are required to confirm this hypothesis. no homologs for the enzyme that produces phosphatidylglycerol were found; however, all mollicutes definitely synthesize this metabolite, and hence must use a hitherto undetected enzyme for this step. the acyl-carrier protein was only found in the genomes of m. hyorhinis and may act as a fatty acid and coa donor in these organisms. orphan reactions for the production of glycolipids were added to models for m. hyopneumoniae andm. flocculare based on the presence of these metabolites in vivo. amino acid metabolism the import of amino acids was added in two forms: oligopeptide import and cleavage or single amino acid import. trna charging accounts for most of the reactions in this pathway. production of the biomass precursor s-adenosyl-methionine from methionine is possible in all strains. moreover, since the enzyme adenosylhomocysteine nucleosidase is present in all species, we added two reactions to restore connectivity between these two subpathways. this addition resulted in the production of two metabolites essential for quorum sensing and cell communication in other species: s-rybosyl-homocysteine and autoinducer- nucleotide metabolism the mycoplasmas in this study cannot synthesize de novo purines and pyrimidines; therefore, they have only salvage pathways and interconversions to supply the cell with nucleic acid precursors. the three species have the same enzymatic capabilities, except for the presence of thymidylate synthetase in m. hyorhinis, which is also important for cofactor metabolism and is responsible for the conversion of dump and -methylene-tetrahydrofolate to dtmp and dihydrofolate. overall, the nucleotide metabolism consists in the uptake of guanine, adenine, uracil, thymine and cytidine and produces all deoxy-ribonucleotides and ribonucleotides. cofactor metabolism around % of the enzymes that comprise the cofactor metabolism of the reconstructed models were not found in any species. for instance, thiamine-pyrophosphate is imported unchanged directly into biomass, while pyridoxal is imported and converted to the cofactors pyridoxal-phosphate by the action of pyridoxal kinase. from the pantothenate/coenzyme-a metabolism, we only found homologs for the conversion of 4’-phosphopantetheine into dephospho-coenzyme-a by pantetheine-phosphate adenylyltransferase. three forms of folate are incorporated to the biomass: methf,, -tetra-hydrofolate and 10-formyl-tetra-hydrofolate. we assumed that folate was imported and converted to these three biomass cofactors by: dihydrofolate reductase, serine hydroxymethyltransferase and formate-tetrahydrofolate ligase. this may also account for the better fitness of m. hyorhinis when compared to the other two species. although no homologs for nad kinase were found, two reactions were added to allow the presence of nadp and nadph, which are cofactors for several essential reactions. the nadph and nadh produced from the degradation of glucose and other carbohydrates are recycled by the action of two enzymes: nadh oxidase and thioredoxin disulfide-reductase. mycoplasmas cultivation we cultivated the three species of mycoplasmas in complex and defined media to gather experimental information and compare to in silico growth. in the cmrl+ medium, m. hyopneumoniae and m. flocculare remained viable only when peptone was present. the presence of peptone appeared to have a negative effect on m. hyorhinis growth, but more tests should be performed to verify this hypothesis. yus+ medium did not allow proliferation in any species, but maintained cell concentration and viability even after days of culture. these results are in agreement with a previous work performed by bertin and colleagues in m. mycoides subsp. mycoides that shows that cmrl- contains all components to support cellular metabolism but not growth. here, the supplemented versions cmrl+ or cmrl+/pep allowed proliferation only for m. hyorhinis during the first h of growth. it seems that after this period of time, one or more essential metabolites initially present were no longer available. by comparing the models reconstructed in this work with the media composition tested for the three species, it might seem plausible that key cofactors might not be delivered to mycoplasmas in the correct form. for example, thiamine pyrophosphate should be directly delivered to all species instead of its precursor thiamine, and pyridoxal-5-phosphate might be the cofactor of choice for the media composition in place of pyridoxal. a list of predicted biomass precursors from the models versus the actual precursors found in the yus+ medium can be seen in additional file 1: table s this new proposed defined medium will be submitted in the near future to cultivation tests with m. hyorhinis, m. hyopneumoniae and m. flocculare. fig. cultivation curves in defined and complex media by species. cell concentrations were estimated by the ccu method and error bars were calculated as the standard deviation between triplicate time-matched samples. as expected, the three species had better growth rates in the complex medium than in defined media. m. hyopneumoniae strain j; m. flocculare strain; m. hyorhinis strain atcc metabolomics of swine mycoplasma nmr analysis of the culture media was performed in order to detect possible differences in the metabolism of the three species. all data measured in defined and complex media can be seen in additional file 1: tables s6a and s6b. the major differences in both complex and defined media were related to the metabolism of pyruvate: m. hyorhinis produced higher quantities of pyruvate at the end of h in complex medium. as a result, pyruvate conversion to acetate was detected in low quantities for this species in both media. m. hyopneumoniae and m. flocculare, on the other hand, produced higher amounts of acetate. low amounts of formate, fumarate and succinate were also detected, and this indicates the presence of genes encoding the enzymes pyruvate formate-lyase, fumarate hydratase and succinyl-coa synthetase in these genomes. as previously mentioned, all models were thus modified to accommodate these activities. formate production seems to be independent from acetate production in m. hyopneumoniae and m. flocculare. however, the similar low concentrations of acetate and formate in both growth media from m. hyorhinis might be an indication that both are produced concomitantly by the action of pyruvate formate-lyase, phosphate acetyl-trasnferase and acetate kinase. fig. distinct products of the metabolism of pyruvate from growth in complex friis medium and yus+ medium of m. hyopneumoniae strains and j, m. flocculare strain and m. hyorhinis strain atcc. in complex medium, we calculated the ratio between the peak signal in cultivated versus control medium and error bars were calculated as the standard deviation between triplicate time-matched samples. for defined medium, we detected the actual concentration for the metabolites and error bars were calculated as the standard deviation between duplicate time-matched samples. a in complex medium, m. hyopneumoniae and m. flocculare can produce high amounts of acetate; the yields are even higher from m. hyopneumoniae. m. hyorhinis, on the other hand, produces low concentrations of acetate in this medium. the final glycolysis product for m. hyorhinis is thus pyruvate. b in defined medium, m. hyopneumoniae and m. flocculare produce similar amounts of acetate while m. hyorhinis contains only residual levels of this metabolite. the three species can produce low amounts of formate in both media to check if the lower concentrations of acetate were a result of an impaired activity of the pyruvate dehydrogenase complex in m. hyorhinis, we compared the predicted proteins pdha, pdhb, pdhc and pdhd from all species. the known active sites of odp and oldh from bacillus subtilis were analyzed in search for a possible mutation. no active sites have been described for pdha and pdhb so far, but the enzymes for all three swine mycoplasmas do not present significant differences. although not all active sites from pdhd seem to exist in mycoplasmas, all predicted proteins from the species analyzed are similar in terms of active sites, when these are present. the only difference between m. hyorhinis and the other two species was found in pdhc, in the protein portion related to the binding to pdha and pdhd. this portion is the most variable between the species. an enzyme characterization should be made in order to verify if this complex is not active in this species, or if it is responsible for another unforeseen activity. the presence of myo-inositol in swine serum might be directly related to the higher production of acetate in the growth of m. hyopneumoniae in complex medium, if compared to m. flocculare. since the myo-inositol pathway produces accoa, the recovery of a molecule of coa from myo-inositol is possible and may be directly linked to the production of acetate. since the myo-inositol catabolic pathway is one of the few distinctions between the metabolic models of m. hyopneumoniae and m. flocculare, it is possible to assume that the differences in acetate concentration in vitro may indeed arise from the ability of m. hyopneumoniae to uptake myo-inositol. indeed, when no source of myo-inositol is present, no difference in the concentration of acetate is observed for these two species. the analysis of the amino acid uptake was not trivial since both media contained peptone, which at first was not measured through the cpmg experiment but, in the course of time, was degraded into single amino acids, and thus changed the overall signal of the nmr spectra. to address this issue, we only took into account the yus+ medium at a single time point after h of growth. m. hyorhinis seemed to have lower concentrations of most amino acids at the end of the growth curve if compared either to the control medium or to the other species. this may be related to several factors, such as lower rates of peptone degradation by membrane proteases or higher uptake rates of amino acids by m. hyorhinis. at this point, we could no longer verify which hypothesis is more suitable for this particular distinction, however, since peptone was not necessary for the maintenance of the viability of m. hyorhinis in the defined medium, it is possible that this species may harbor transporters with higher specificity or effectiveness for single amino acids. the results as concerns the other metabolites indicate for the three species the uptake of glucose and nucleotides. flux balance analysis fba was used to check the properties and capabilities to produce detected metabolites by nmr spectra of all models. fba simulations were performed with an in silico medium that contained all biomass precursors. similarly to the in vitro growth results, in the presence of oxygen, acetate was the main product and nadh was recycled by the conversion of molecular oxygen into water. in the absence of oxygen, all models were able to grow but the main product became lactate. as expected, the availability of myo-inositol in the in silico medium allowed an alternative source of coa and higher concentrations of acetate as a final product in the models for m. hyopneumoniae. hydrogen peroxide production was linked to the metabolism of glycerol in both m. hyopneumoniae and m. hyorhinis, while the models for m. flocculare did not produce toxic levels of this metabolite. since we have no information on gene essentiality for any species, we checked reaction essentiality in the pan-network. this means that we deleted individual reactions even if they did not have a gpr association. a total of reactions were essential for growth in the pan-network. discussion in this work, we created a metabolism as realistic as possible for the three known mycoplasmas present in the respiratory tract of swine. it is essential to point out that, although we had to include % of orphan reactions in order to allow growth, we only used about half of the genomes of these organisms. the other half consists either of hypothetical or of conserved hypothetical proteins, which may in part fill the missing gaps of the models. the number of essential reactions in the pan-network in this work is not directly comparable to the number of indispensable genes predicted for m. genitalium or m. pneumoniae. gene essentiality accounts for more than metabolic enzymes which means that proteins not included in our models, related to protein synthesis, dna polymerization or rna turnover, were not accounted for in the essential reactions of our models. if we added these proteins, we would possibly arrive at numbers closer to those of m. genitalium and m. pneumoniae. the lack of experimental information on gene essentiality is also a setback for the validation of the models created here. further experiments on this matter should help us better refine and complete the networks. as for most of the mycoplasma species studied, all reconstructed networks exhibit low connectivity due to the simplicity of the biological model. we were able to show in this work that the three swine mycoplasma species have similar metabolic capabilities, except for the metabolism of myo-inositol, amino sugar, carbohydrates and glycerol. overall, the methods used here enabled us to address some of the main problems caused by most automatic reconstruction methods which are the permissive inclusion of pathways and over prediction of capabilities. annotation errors arise with the attribution of ambiguous or partial ec numbers and the propagation of these errors may then lead to many other ones. thus, it is known that the prediction of gpr associations based only on name-matching and ec-codes is not sufficient to add confidence to a model. when we decided to simultaneously refine the models, we enhanced the confidence of each gpr association, by adding information on the synteny between the genomes, protein sequence alignments and phylogenetic distance between orthologs. fig. differential metabolism of m. hyorhinis, m. hyopneumoniae and m. flocculare. metabolites are depicted in dark green and separate enzymatic activities for m. hyorhinis, m. hyopneumoniae and m. flocculare can be seen in yellow, pink and blue, respectively. whenever an enzyme is missing from the three species, the enzyme rectangle is depicted in grey. complete list of metabolite abbreviations and ec numbers can be found in additional file 4 our main objective was not only to reconstruct the metabolic models for m. hyorhinis, m. hyopneumoniae and m. flocculare. we also wanted to compare the metabolism of the three species and find possible links to virulence, host-colonization capacity and life-style. m. hyorhinis, for instance, is not only a pathogen, but also a common contaminant in many mammalian cell cultures. this may be explained by the wide range of carbohydrates which this organism can uptake and metabolize. this may be a crucial factor to the ability of m. hyorhinis to grow in diverse sites inside the host and even to invade other hosts and potentially develop cancer in humans. m. hyopneumoniae has been detected in other sites such as brain, liver and spleen, but m. flocculare has never been detected outside of the respiratory tract of swine. one may wonder what could explain this if both organisms are so similar. m. hyopneumoniae is up to date the only mycoplasma species with sequenced genome that has the genes for the catabolism of myo-inositol. myo-inositol is readily abundant in the bloodstream of mammalian hosts, and can be used as a secondary carbon source for energy production. mycoplasma iguanae has been described to produce acid from inositol. however, there is no available complete genome sequence for this organism, making impossible any comparison of the genes involved in this pathway. although at this point we cannot confirm that this pathway is functional in m. hyopneumoniae, the nmr results we obtained suggest that myo-inositol might be directly related to the higher production of acetate in complex medium of m. hyopneumoniae if compared to m. flocculare. this assumption is based on two factors: first, a previous work has detected myo-inositol in swine serum, which is a component of the complex medium; and second, the only difference in the networks that could influence the acetate concentration is the myo-inositol catabolism. this pathway may also be involved with an alternative production of accoa, precursor to coa, an essential cofactor for growth in all species. these metabolic distinctions may help m. hyopneumoniae to grow in diverse sites if compared to m. flocculare. this is particularly interesting when we take a close look into the composition of lung surfactant from several mammalian models, which likely contains myo-inositol as a degradation product from phosphatidylinositol. thus, it is possible that the myo-inositol catabolism in m. hyopneumoniae is one of the reasons for the high virulence of this species when compared to m. flocculare and m. hyorhinis. a striking result from the nmr results was the reduced capability of m. hyorhinis to produce acetate. this was not foreseen in the model reconstruction, but predicted protein sequences from the pyruvate dehydrogenase complex of the three species showed a particular distinction in pdhc, more specifically in the region responsible for binding to other complex components. we could not verify at this point if this complex is not active in the particular growth conditions of our work in m. hyorhinis or if it is responsible for another unforeseen activity. however, a transcriptome profiling of m. hyorhinis has detected both pdha and pdhb in the pool of genes with the highest number of transcript reads, indicating that the complex might be translated in vivo. moreover, pyruvate dehydrogenase activity has been previously detected in m. hyorhinis extracts, and it seems to be correlated to oxygen availability. in our cultivations, the cells are not grown in a complete aerobic system, which may explain the differences between our findings and those previously published. cultivation tests in vitro also showed that none of the species was able to grow in none of the defined media. this is probably due to a lack of the correct cofactors as previously mentioned. for instance, thiamine, pyridoxal/pyridoxine, pantothenate, spermine and folate were the actual media components tested, but the reconstructed models were not able to convert them to biomass precursors. only m. hyorhinis seems to be able to uptake folate directly while the other two species might need intermediate metabolites. this might also explain why m. hyorhinis is extensively found as a cell culture contaminant. even though we supplemented cmrl- with all missing components present in the yus+ medium, cells maintained viability for longer periods in the latter. cofactors present originally in both media existed in lower concentrations in the cmrl- this means that possibly one or more cofactors with higher concentrations in the yus+ medium are essential for the viability of the three species. common cofactors with higher concentration in yus+ medium included: choline, folate, pantothenate, pyridoxal, thiamine and spermine. other differences among the species seem to indicate that m. hyopneumoniae and m. flocculare lack one or more amino acid transporters. this may be due to two things: both remained viable for longer periods whenever cmrl+ was supplemented with peptone, and we detected a higher amino acid accumulation in the defined media if compared to m. hyorhinis. m. hyorhinis on the other hand may have all the transporters for single amino acids, but this hypothesis needs to be further verified. from sequence analyses, we did not find significant differences among the species that would explain this distinct behavior. moreover, since m. hyorhinis seemed to overgrow in the first h of culture in both cmrl media, it is possible that one or more compounds were missing in the yus+ medium in order to allow growth of this species. possible candidates for such would be: coa, nicotinamide, ascorbate and/or pyridoxine. it is still not clear if any of the swine mycoplasmas studied here indeed have a polysaccharide capsule, given that only one study has assesed its presence in m. hyopneumoniae. however, this is particularly important because in this study tajima and yagihashi, as discussed before, reported that capsular polysaccharides from m. hyopneumoniae play a key role in the interaction between pathogen and host. moreover, it has been reported that some strains of m. hyopneumoniae become less pathogenic in broth culture and, after serial passages, lose their ability to produce gross pneumonia in pigs. in several bacterial species, it has been established that the amount of capsular polysaccharide is a major factor in virulence and it decreases significantly with in vitro passages. furthermore, despite increasing evidence supporting the existence of a polysaccharide capsule in several mycoplasma species, it is not yet known if the material is synthesized by the bacteria or imported from the medium for most species. recent reports have even shown the capacity of species from the mycoides cluster to assemble a capsule from exogenous phosphorylated glucose. however, no homolog proteins related to these activities were found in none of the species studied in this work. for the reasons discussed above, we were not able to introduce the capsule production itself in the models; however, we tried to relate its possible existence with the substrates that might be used to produce it. the composition of capsular polysaccharides has been related to the level of pathogenicity in other bacteria. for instance, highly virulent escherichia coli strains seem to have practically nonimmunogenic capsular material, due to the fact that these antigens are similar or identical to the ones found in the host. the structure of the k antigen from e. coli capsules is basically glcnac and glucuronic acid in a molar ratio of 1:. erlinger and collaborators identified heparan sulfate as the predominant glycosaminoglycan in the porcine respiratory tract and the most common disaccharide unit within hs is glcnac linked to glca. since m. hyopneumoniae and m. flocculare might be able to import glcnac unchanged, they may use it for the composition of its capsule; m. hyorhinis, on the other hand, imports sialic acid, and might use it directly to produce its own. of course this is highly speculative, and we are not sure if all these hypotheses indeed happen in vivo. the metabolism of glycerol is the final major difference between the pathogenic m. hyopneumoniae and m. hyorhinis and the non-pathogenic m. flocculare. such metabolism and the production of hydrogen peroxide by the action of glpo are essential for the cytotoxicity of m. pneumoniae and m. mycoides subsp. mycoides. particularly, the cytotoxicity of m. mycoides subsp. mycoides is due to the translocation of the hydrogen peroxide into the host cells. this is only possible because of the close proximity to the host cells along with the membrane-bound enzyme. highly conserved homolog genes to this enzyme were only found in the genomes of the pathogenic species studied here. whether or not these enzymes from m. hyopneumoniae and m. hyorhinis are indeed capable of the same activities as glpo is not yet confirmed; however, the high similarity may be an indication that this trait is essential for the pathogenicity of these two species. this may also explain why both m. hyopneumoniae and m. flocculare can adhere to the cilia of tracheal epithelial cells, but only the adhesion of m. hyopneumoniae results in tissue damage. moreover, while m. hyorhinis has in the reconstructed models three ways of uptaking glycerol, m. hyopneumoniae seems to have five, and this may reflect in its enhanced pathogenicity.  we presented in this work an overview of the differential metabolism of m. hyopneumoniae, m. hyorhinis and m. flocculare using different approaches. the reconstructed models showed some distinctions among the species, namely the myo-inositol metabolism for m. hyopneumoniae, the wider uptake of carbohydrates for m. hyorhinis and the usage of glycerol as a carbon source for the two pathogenic species. the models also served as a basis for all the assumptions made for the experimental data. metabolic profiling of both complex and defined media pointed to new differences that we were not able to identify based solely on the sequenced genomes. the major ones were related to the pyruvate conversion to acetate, which appeared to be higher in m. hyopneumoniae and m. flocculare than in m. hyorhinis. after we performed growth tests using defined media that hypothetically would lead to bacterial growth in silico, we assessed that some cofactors or metabolites were probably being delivered in the wrong form to the mycoplasmas. the environmental context may explain the differences between the in silico models and the in vivo behavior. whether the main differences among the species we reported here are related to virulence or pathogenicity have not yet been addressed experimentally, but it is tempting to speculate. the same factors that may enhance virulence of m. hyopneumoniae may help the commensal species m. flocculare to better survive inside the host. all these in silico and in vivo metabolic differences among m. hyopneumoniae, m. hyorhinis and m. flocculare might influence the different levels of pathogenicity in each of them. however, it is highly likely that gene regulation may directly interfere in the metabolism and pathogenicity and may be related to many other aspects still unaccounted for. one of our future goals is therefore to understand and integrate gene regulation into the metabolic models. upcoming experiments will aim at testing the hypotheses formulated here, particularly those related to the metabolisms of glycerol and myo-inositol. we also intend to better understand the habitat of these species, and the possible metabolic and genetic dialogues with the host and other bacteria present in this environment. either way, this work serves as a basis for the study of the differential metabolism and pathologies caused by the swine lung mycoplasmas and may help to propose ways to prevent disease development in the future. availability of data and material all models from this work are available in additional file and online at the metexplore database. additional files additional file supplementary tables. table s defined media composition used for the cultivation of m. hyopneumoniae, m. hyorhinis and m. flocculare. table s complete list of reactions and metabolites from all models. table s list of all metabolites included in the final biomass composition. table s distribution of gpr associations and exclusive gpr associations by metabolic pathways. table s proposition of correct precursors and cofactors for a new defined medium. table s nmr metabolomics: metabolites detected in the complex friis and yus+ media for the three swine mycoplasmas. table s prediction of essential reactions in the pan-network. table s differential metabolism overview. additional file biomass reaction assembly from reconstructed models. pdf file containing all steps necessary to the assembly of the biomass reaction for all species. additional file abbreviation of metabolites and ec numbers from figs. and spreadsheet containing the lists of all abbreviations and ec numbers used in figs. and separately. additional file prediction of transporters. pdf file containing all assumptions for the in silico prediction of transporters. additional file models in sbml format. compressed file containing all sbml models from this work. abbreviations atp-binding cassette; acpacyl-carrier protein blastpblast local alignment search tool demdead-end metabolite ccucolor changing unit cpmgcarr-purcell-meiboom-gill ec numberenzyme-code number fbaflux balance analysis fidfree induction decay gamgrowth associated maintenance gogene ontology gprgene-protein-reaction association hcahierarchical clustering analysis hsqc spectroscopyheteronuclear single-quantum correlation spectroscopy mifmajor facilitator protein mwmolecular weight ngamnon-growth associated maintenace nmrnuclear magnetic resonance noesynuclear overhauser effect spectroscopy pcaprincipal component analysis pfam databaseprotein-families database ptsphospho-transferase transport system sbmlsystems biology markup language tmsptrimethilsilylpropionate tocsytotal-quantum correlation spectroscopy tutranscriptional unit. enzyme abbreviations: abc-gly abc-glycerol transport system glyfglycerol facilitator protein glpoglycerol-3-phosphate glpuglycerophosphodiester transporter fbafructose-biphosphate-aldolase iolj5-dehydro-2-deoxy-phosphogluconate aldolase ldhlactate dehydrogenase mdhmalate dehydrogenase pdh a/b/c/dpyruvate dehydrogenase complex subunit a/b/c/d ulagl-ascorbate-6-phosphate lactonase metabolite abbreviations used in the text: 2kg2-oxoglutarate accoaacetyl-coenzyme-a coacoenzyme-a dhadihydroxyacetone dhapdihydroxyacetone-phosphate f6pfructose-6-phosphate glcnacn-acetyl-glucosamine nad/nadhnicotinamide adenine dinucleotide nadp/nadphnicotinamide adenine dinucleotide phosphate methf -methylene-tetrahydrofolate succoasuccinyl-coenzyme-a competing interests the authors declare that they have no competing interests. authors’ contributions mfs and az conceived and designed the work in close collaboration with all the other authors. mgf performed all computational metabolic reconstructions and most of the experimental work. fms, sgm, tlp and ej collaborated in the experimental work. tlp performed the metabolomics statistic analysis. beh, ft, atrv and iss collaborated analyzing all data. mgf, mfs and az wrote the manuscript with inputs from all the other authors. all authors read and approved the final manuscript. 
this activity, missing only in m. flocculare, may be directly involved in cytotoxicity, as already described for two lung pathogenic mycoplasmas, namely mycoplasma pneumoniae in human and mycoplasma mycoides subsp. the respiratory tract of swine is colonized by several bacteria among which are three mycoplasma species: mycoplasma flocculare, mycoplasma hyopneumoniae and mycoplasma hyorhinis. we were able to infer from the reconstructed networks that the lack of pathogenicity of m. flocculare if compared to the highly pathogenic m. hyopneumoniae may be related to its incapacity to produce cytotoxic hydrogen peroxide. the genomic resemblance among these three mycoplasma species combined with their different levels of pathogenicity is an indication that they have unknown mechanisms of virulence and differential expression, as for most mycoplasmas. while colonization by m. flocculare is virtually asymptomatic, m. hyopneumoniae is the causative agent of enzootic pneumonia and m. hyorhinis is present in cases of pneumonia, polyserositis and arthritis. in addition, the myo-inositol catabolism is exclusive to m. hyopneumoniae and may be an important trait for virulence. even though the refined models have similar metabolic capabilities, interesting differences include a wider range of carbohydrate uptake in m. hyorhinis, which in turn may also explain why this species is a widely contaminant in cell cultures. moreover, the ability of m. hyorhinis to grow in diverse sites and even in different hosts may be a reflection of its enhanced and wider carbohydrate uptake. altogether, the metabolic differences highlighted in silico and in vitro provide important insights to the different levels of pathogenicity observed in each of the studied species. cultivation tests and metabolomic experiments through nuclear magnetic resonance spectroscopy were also performed to acquire experimental data and further refine the models reconstructed in silico. however, the most important difference seems to be related to glycerol conversion to dihydroxyacetone-phosphate, which produces toxic hydrogen peroxide. metabolomic data suggest that even though these mycoplasmas are extremely similar in terms of genome and metabolism, distinct products and reaction rates may be the result of differential expression throughout the species.
 the respiratory tract of swine is colonized by several bacteria among which are three mycoplasma species: mycoplasma flocculare, mycoplasma hyopneumoniae and mycoplasma hyorhinis. while colonization by m. flocculare is virtually asymptomatic, m. hyopneumoniae is the causative agent of enzootic pneumonia and m. hyorhinis is present in cases of pneumonia, polyserositis and arthritis. the genomic resemblance among these three mycoplasma species combined with their different levels of pathogenicity is an indication that they have unknown mechanisms of virulence and differential expression, as for most mycoplasmas. methods in this work, we performed whole-genome metabolic network reconstructions for these three mycoplasmas. cultivation tests and metabolomic experiments through nuclear magnetic resonance spectroscopy were also performed to acquire experimental data and further refine the models reconstructed in silico.  even though the refined models have similar metabolic capabilities, interesting differences include a wider range of carbohydrate uptake in m. hyorhinis, which in turn may also explain why this species is a widely contaminant in cell cultures. in addition, the myo-inositol catabolism is exclusive to m. hyopneumoniae and may be an important trait for virulence. however, the most important difference seems to be related to glycerol conversion to dihydroxyacetone-phosphate, which produces toxic hydrogen peroxide. this activity, missing only in m. flocculare, may be directly involved in cytotoxicity, as already described for two lung pathogenic mycoplasmas, namely mycoplasma pneumoniae in human and mycoplasma mycoides subsp. mycoides in ruminants. metabolomic data suggest that even though these mycoplasmas are extremely similar in terms of genome and metabolism, distinct products and reaction rates may be the result of differential expression throughout the species.  we were able to infer from the reconstructed networks that the lack of pathogenicity of m. flocculare if compared to the highly pathogenic m. hyopneumoniae may be related to its incapacity to produce cytotoxic hydrogen peroxide. moreover, the ability of m. hyorhinis to grow in diverse sites and even in different hosts may be a reflection of its enhanced and wider carbohydrate uptake. altogether, the metabolic differences highlighted in silico and in vitro provide important insights to the different levels of pathogenicity observed in each of the studied species. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords mycoplasmamollicutesmetabolic networkmetabolismwhole-genome metabolic reconstructionhydrogen peroxidehttp://dx.doi.org/ /501100000781european research councilfp /2007-2013/ 10schrank irene s. capes-cofecub782/13schrank irene s. capes-cofecub782/13tardy florence tgir-rmn-thc fr cnrscapes-cofecub782/13elena-herrmann bénédicte issue-copyright-statement© the author 2016 