Skip to main content

Microscopy and genomic analysis of Mycoplasma parvum strain Indiana

Abstract

Mycoplasma parvum [Eperythrozoon parvum] is the second hemotrophic mycoplasma (hemoplasma) described in pigs. Unlike M. suis, its closest phylogenetic relative, M. parvum, is considered a non-pathogenic bacterium in this host species. Natural infection of a domestic, 6-month-old splenectomized pig with M. parvum strain Indiana is described herein. Light and scanning electron microscopy of the bacteria were performed in addition to whole genome sequencing, analysis, and comparison to the genome of M. suis strain Illinois. Neither clinical signs nor anemia were observed during the infection. Microscopy analyses revealed coccoid to rod- shaped organisms varying from 0.2 to 0.5Ā Ī¼m; they were observed individually or in short chains by both light and electron microscopy, however less than 30% of the red blood cells were infected at peak bacteremia. The single circular chromosome of M. parvum was only 564 395Ā bp, smaller than M. genitalium, previously considered the tiniest member of the Mollicutes. Its general genomic features were similar to others in this class and species circumscription was verified by phylogenomic analysis. A gene-by-gene comparison between M. suis and M. parvum revealed all protein coding sequences (CDS) with assigned functions were shared, including metabolic functions, transporters and putative virulence factors. However, the number of CDS in paralogous gene families was remarkably different with about half as many paralogs in M. parvum. The differences in paralogous genes may be implicated in the different pathogenic potential of these two species, however variable gene expression may also play a role. Both are areas of ongoing investigation.

Introduction

Hemotrophic mycoplasmas (class Mollicutes), trivially known as hemoplasmas, are small cell-wall-less bacteria with a unique tropism for the red blood cells (RBC) of a variety of vertebrate hosts. The genome of ten hemoplasmas, eight different species and two strains each of Mycoplasma suis and M. haemofelis, had their genomes sequenced to date; seven of these were completed by our research group [1]ā€“[10]. Their genomes have a small core of essential genes; however, they rely on the acquisition of most nutrients from their environment to survive. This feature likely plays a key role in the lack of an established in vitro cultivation of these organisms.

Mycoplasma parvum [Eperythrozoon parvum] was first reported in 1950 as a nonpathogenic bacterium of the pig [11]. Unlike M. suis, its closest relative, M. parvum often accumulates on a single RBC but only a few RBC are infected. Rod and coccoid forms are common in M. parvum infection, ranging in size from 0.2 to 0.5Ā Ī¼m, whereas ring forms are frequent in M. suis. The latter form is larger and may be up to 2.5Ā Ī¼m in diameter [11]. Infection by M. parvum is usually characterized by the absence of clinical signs even in splenectomized pigs; on the contrary, M. suis bacteremia, characterized by fever and icteroanemia, may result in death. In addition, cross-inoculation studies differentiate these two hemoplasmas as separate species as well as from other species of red cell bacteria that infect cattle, sheep, goats, and mice [11],[12].

More than 50Ā years have passed since its first description, and yet less than 10 publications can be found when searching for ā€œMycoplasma parvumā€ or ā€œEperythrozoon parvumā€ keywords in scientific databases. The two previous reports documenting M. parvum infection in pigs were based solely on the amplification and sequencing of its 16S rRNA or ribonuclease P RNA genes [13],[14]. The sequences of these genes were divergent from those of M. suis, sharing only 96% and 88% identity, respectively. These molecular studies corroborated previous morphology, pathogenicity, and cross-inoculation studies, suggesting that M. parvum and M. suis were distinct hemoplasma species. Nevertheless, new molecular tools for species circumscription are likely to provide a better resolution for this complex definition, especially when dealing with these uncultivable bacteria. The average nucleotide identity (ANI) analysis of conserved and shared genes between bacterial strains and supporting calculation of tetranucleotide signature correlation index (TETRA), has recently been proposed as a reliable substitute for DNA-DNA hybridization (DDH) [15]. A value for species delineation based on ANI is currently set at 95%-96%. The TETRA signature is an alignment-free parameter and generally correlates with ANI.

We report herein the course of infection with M. parvum in a naturally infected, splenectomized 6-month-old pig, and the complete genome of M. parvum strain Indiana. Furthermore, comparative genomic analyses of M. parvum with other hemotrophic mycoplasmas were performed, as well as a gene-by-gene comparison between M. parvum strain Indiana and M. suis strain Illinois to better understand their contrasting pathogenicity.

Material and methods

Animal and Mycoplasma parvum DNA isolation

A male domestic 6-month-old pig (Sus scrofa domesticus, mixed breed) was purchased from the Animal Sciences Research and Education Center-ASREC (West Lafayette, IN, USA). While at the farm, this animal was kept in close contact with other pigs fed a commercial feed containing chlortetracycline (100 grams/ton) during growing and finishing stages (from 7ā€“8 weeks to 6Ā months old), and vaccinated against Mycoplasma hyopneumoniae (RespiSureĀ®/ER Bac PlusĀ®, Pfizer), Erysipelothrix rhusiopathiae, and circovirus. The animal was tested for hemotrophic mycoplasma infection [16] on two separate occasions prior to its arrival at the Purdue Animal Housing Facilities (Purdue University, West Lafayette, IN, USA). A conventional PCR (cPCR) developed in our laboratory, which amplifies a larger product (~840Ā bp) of the 16S rRNA gene of swine hemoplasmas, was also performed [17]. The identity of both qPCR and cPCR products was confirmed by Sanger sequencing. The animal was treated according to the Purdue Animal Care and Use Committee (PACUC) protocol number 1111000223. The pig was fed antibiotic-free commercial feed, as well as water (via a self-controlled nipple waterer) ad libitum during the entire study. Eleven days after its arrival, the animal was splenectomized according to the PACUC protocol. The pig was monitored daily (minimum twice a day) for clinical signs (e.g. elevated body temperature, and direct observation of behavior- BAR status), and blood was collected into EDTA tubes before and following splenectomy (every 2ā€“8 days depending on clinical signs and blood smear evaluation) for monitoring M. parvum infection by qPCR [16].

M. parvum strain Indiana was harvested from the blood of the pig at the peak of bacteremia; EDTA-whole blood was centrifuged at 4000ā€‰g for 10Ā min and the buffy coat was removed. The remaining red blood cells with M. parvum organisms attached were used for the subsequent DNA extraction. Genomic DNA was extracted using Quick-gDNAā„¢ MidiPrep kit according to the manufacturerā€™s instructions (Zymo Research, Irvine, CA, USA).

Light and scanning electron microscopy (SEM) of M. parvum

Blood smears were prepared from fresh EDTA-blood every 2Ā days after splenectomy to follow M. parvum infection and stained with Giemsa. Photomicrographs were taken using a total magnification of 100 X.

For SEM, aliquots of 1.0Ā mL of blood infected with M. parvum were centrifuged and pellet was resuspended in 1.0Ā mL of glutaraldehyde 2.5% in 0.1Ā M sodium cacodylate buffer, pHĀ 7.4, for fixation. Samples were prepared for scanning electron microscopy using standard procedures at Purdue Universityā€™s Biological Electron Microscopy Facility, and visualized using FEI Titan Krios microscope.

Sequencing and assembly of M. parvum strain Indiana genome

The whole genome was sequenced from paired-end libraries (TruSeq DNA sample preparation kit, Illumina, San Diego, CA, USA) using Illumina HiSeq 2000 (Illumina, Inc., San Diego, CA, USA) at Purdue University Genomics Core Facility. Average reads of about 100 bases were assembled using ABySS 1.2.7 [18]. After assembly resulting from 1000ā€‰Ć—ā€‰genome coverage of the Illumina reads, a single remaining gap was closed using conventional PCR followed by Sanger sequencing in both directions.

Annotation and analyses of the complete genome of M. parvum strain Indiana

NCBIā€™s prokaryotic Genomes Annotation Pipeline 2.0 provided the first pass annotation of M. parvum genome. The annotation tool Manatee (Institute for Genome Sciences (IGS), School of Medicine, University of Maryland) was used to perform manual annotation of the genome and comparative analysis between the genomes of M. parvum and M. suis strain Illinois (described below). Genomic data of other mycoplasmas available in the NCBI database (NCBI, Bethesda, MD, USA) were also used for comparative analyses.

Metabolic pathways were predicted based on the KEGG pathway database [19] and the study reported by Yus et al. [20]. BLASTclust tool, by Max-Planck Institute for Developmental Biology [21], was used to assign the paralogous gene families with thresholds of 70% covered length and 30% sequence identity.

Analyses for species differentiation

JSpecies software was used to calculate the average nucleotide identity (ANI; MUMmer algorithm) and tetranucleotide signature correlation indexes between selected genomes as previously described [15]. The following genome sequences were used in the analyses: M. parvum strain Indiana [CP006771.1], M. suis strain Illinois [CP002525.1], M. suis strain KI3806 [FQ790233.1], M. haemofelis strain Ohio2 [CP002808.1], M. haemofelis strain Langford 1 [FR773153.2], M. haemocanis strain Illinois [NC_016638.1], M. wenyonii strain Massachusetts [NC_018149.1], ā€œCandidatus M. haemominutumā€ strain Birmingham 1 [HE613254.1], and ā€œCandidatus M. haemolamaeā€ strain Purdue [NC_018219]. The thresholds for species circumscription are 94% and 0.99 for ANIm and tetranucleotide indexes, respectively [15].

Comparative genomics of M. parvum and M. suis strain Illinois

Comparative analyses of the whole genome of M. parvum and M. suis strain Illinois were performed using BLASTp and/or BLASTn of each CDS or gene of one genome against the other and vice versa. BLASTp was applied to identify unique CDS of M. parvum or M. suis; a CDS was considered unique when there were no matching sequences withā€‰ā‰„ā€‰80% coverage andā€‰ā‰„ā€‰40% orā€‰ā‰„ā€‰90% coverage andā€‰ā‰„ā€‰30% identity to the query sequence.

Each unique CDS of M. parvum and M. suis was analyzed for the identification of the following parameters: subcellular localization and protein sorting signals using PSORTb v3.0.2 [22],[23]; signal peptide cleavage sites using SignalP 4.1 [24]; and the presence of lipoproteins using LipoP version 1.0 software [25].

Whole genome synteny (gene order) was compared between M. parvum strain Indiana and M. suis strain Illinois using SynMap from CoGe [26]. SynMap generates two-dimensional dot-plot synteny maps using a DAGchainer algorithm coupled with BLAST to identify syntenic homologous genes; each dot represents putative homologous genes between any two genomes [27].

Results

Mycoplasma parvum infection in a splenectomized pig: bacterial loads and clinical signs

A 6-month-old barrow, domestic pig naturally infected with M. parvum was used in this study. The dynamics of its M. parvum blood load is shown in FigureĀ 1. Briefly, this animal was initially identified as positive using a specific TaqMan qPCR for swine hemoplasmas [16] while still at the farm (ASREC) (days 1 and 3, FigureĀ 1). In addition, a cPCR employed for amplification of ~840Ā bp fragment of the 16S rRNA gene of swine hemoplasmas was also positive [17]. Sequencing of the qPCR (157Ā bp) and cPCR (840Ā bp) products showed 100% and 99% identity, respectively, to the 16S rRNA sequence described for M. parvum (GenBank: JX489599.1), in contrast to 97% identity when compared to the 16S rRNA sequence of M. suis strain Illinois (NCBI: NR_103930.1). Moreover, the sequencing chromatograms showed single peaks along the entire length of the amplicons indicating the presence of a single sequence. These results confirmed that the animal was infected with M. parvum and free of M. suis. At day 22, this pig was transferred to the Purdue Animal Housing Facilities, where it remained for 48Ā days, from April 11th (day 22) to May 29th (day 70) 2013 (FigureĀ 1). The animal was then splenectomized 11Ā days after its arrival (day 33, FigureĀ 1) according to the Purdue Animal Care and Use Committee protocol (1111000223). Bacterial loads varied from 105 to 1010 organisms/mL of blood throughout the course of infection with a peak of bacteremia (1010 organisms/mL of blood) occurring 12ā€“14 days after its splenectomy (days 45ā€“47, FigureĀ 1). Organisms were only detected on the peripheral blood smears at peak and for 24Ā h thereafter. While the pigā€™s body temperature minimally elevated to 103.7Ā Ā°F (reference range 101.6-103.6Ā Ā°F) [28] at the peak, no overt clinical signs of infection were observed. The pigā€™s hematocrit varied from 30% (only on one occasion, not at the peak of bacteremia; FigureĀ 1) to 43%, but otherwise remained within the reference interval (32-50%) throughout the course of infection [29].

Figure 1
figure 1

M. parvum loads and PCV of a naturally infected pig before and after splenectomy. The loads of M. parvum strain Indiana organisms were detected by TaqMan qPCR on blood samples and the (packed cell volume) PCV was determined before and after splenectomy in the naturally infected pig.

Light and Scanning Electron Microscopy (SEM) of M. parvum organisms

A photomicrograph of the peripheral blood smear obtained at the peak of bacteremia with M. parvum (day 45, FigureĀ 1) shows the shape and distribution of these bacteria (FigureĀ 2A). As described previously [11], the organisms are small (approximately 0.2-0.5Ā Ī¼m in diameter), and often infect a few RBC with large numbers of organisms on a single RBC (FigureĀ 2). The SEM shows three M. parvum organisms on the surface of a red blood cell; their size varied from 0.2-0.5Ā Ī¼m in diameter (FigureĀ 2B).

Figure 2
figure 2

Light microscopy and SEM of M. parvum strain Indiana organisms. A) Photomicrograph of a blood smear stained with Giemsa obtained from a pig infected with Mycoplasma parvum (black arrows) at the peak of bacteremia (100X magnification). B) Scanning electron micrograph (SEM) showing red blood cells (RBC) of the pig with M. parvum organisms (black arrow) on the RBC surface. Scale bar, 1.0Ā Ī¼m.

General features of the genome of M. parvum strain Indiana

The general features of the genome of M. parvum strain Indiana regarding size, G + C content, percentage of coding sequences and paralogous genes were within the range of those described for all hemoplasmas sequenced to date [1],[4]ā€“[8],[10] (TableĀ 1). M. parvumā€™s single circular chromosome of 564.395 Kb, the second smallest genome amongst all hemoplasmas, and even smaller than that of M. genitalium, which was previously considered the tiniest member of the Mollicutes (TableĀ 1). It has a G + C content of 27%, and like other mycoplasmas, it appears to use the opal stop codon (UGA) for tryptophan. The 16S, 23S, and 5S rRNA genes are represented as single copies; however, like that of Mycoplasma suis, the 16S rRNA gene is separated from the 5S-23S rRNA in a different operon [1],[2]. Thirty-two tRNA were identified covering all amino acids.

Table 1 General features of the genome of Mycoplasma parvum strain Indiana compared to other Mycoplasma species of the pneumoniae group.

A total of 581 protein-coding sequences (CDS) were predicted and putative functions were assigned and manually verified using the Manatee annotation pipeline: 287 CDS have putative functional identities, representing almost 50% of the CDS, while the other 50% were represented by hypothetical proteins. Further, 24.3% of the genome is dedicated to duplicated genes organized in paralogous families, mostly composed of hypothetical CDS (TableĀ 1).

Hemoplasmas and their phylogenomic relationship based on ANI and Tetranucleotide signature indexes

A phylogenomic comparison based on average nucleotide identity (ANI) and tetranucleotide signature correlation indexes (tetra) was performed amongst selected hemoplasmas with sequenced genomes (TableĀ 2). M. parvum had an ANI varying from 83.6% to 89.5% when compared to other hemoplasma genomes, including M. suis. The tetra correlation indexes of M. parvum with other genomes had a range of 0.47 to 0.88. Both ANI and tetra correlation indexes were below the cutoff values of 94% and 0.99, respectively, for species circumscription.

Table 2 Average nucleotide identity* (ANI) and tetranucleotide signature (Tetra) correlation indexes of selected hemotrophic mycoplasmas.

Comparing all these species of hemoplasmas, ANI and tetra correlation indexes were between 83.6%-90.25% and 0.36-0.96, respectively, correctly separating these organisms as different species of mycoplasmas. On the contrary, strains of the same species (M. suis Illinois and KI3806; M. haemofelis Ohio2 and Langford1) showed ANI and tetranucleotide correlation indexes above the proposed thresholds for species definition, as expected.

M. parvum versus M. suis: similarities and differences at the genomic level

A complete comparison at the genomic level was performed between the genomes of M. parvum strain Indiana and M. suis strain Illinois (FigureĀ 3). The genome of M. suis is 178 Kb bigger than that of M. parvum. SynMap analysis (Additional file 1) indicated a lack of overall gene synteny between M. suis and M. parvum. Although some conserved blocks of genes, mostly corresponding to operons (e.g. PTS system), were observed, they do not necessarily occur at the same relative position of the genomes. A more detailed, manual analysis shows that the 16S rRNA gene and 5S-23S rRNA operon of M. parvum are located in the first half of its genome (dnaA at position zero), while M. suis copies are in the second half of its genome. Moreover, most of the CDS classified by TIGR role category are organized in a different fashion when comparing both genomes (FigureĀ 3). Despite these rearrangements, both genomes share all CDS with known assigned metabolic or other functions, such as transporters and putative virulence factors; this represents 49.4% and 34.7% of the M. parvum and M. suis genomes, respectively. The remainder of both genomes consisted in hypothetical CDS, most of which are dedicated to paralogous gene families: M. parvum has 141 CDS (24.3%) distributed into 38 paralogous gene families, while M. suis has 361 CDS (42.8%) in 68 paralogous gene families (TableĀ 1). Twenty-three out of the 38 paralog families of M. parvum are represented in the genome of M. suis (with a variable number of members per family between the genomes). In contrast, M. suis has 40 out of 68 families that are found exclusively in its genome. This difference in paralogous genes is represented by 188 exclusive CDS in the genome of M. suis and 64 CDS unique of M. parvum. Considering the hypothetical CDS that are not in paralogous families, both genomes share 74 of these, whereas 63 CDS are found only in M. parvum and 153 CDS only in M. suis. In total, M. parvum has a set of 127 unique CDS, while M. suis has 341 CDS found exclusively in its genome.

Figure 3
figure 3

Circular representation of the genomes of M. parvum strain Indiana and M. suis strain Illinois showing the organization of the coding sequences (CDS). The dnaA gene is at position zero in both genome plots. Outer to inner circles: circle 1: predicted CDS on the positive strand; circle 2: predicted CDS on the negative strand. Each CDS is classified by TIGR role category according to the color designation in the legend below the plots; circle 3: CDS in the largest paralogous gene families with each family represented by a different color in each genome, non-paralogous CDS are light grey. Paralogous families of M. parvum with less than 5 CDS are represented in orange. Black and red marks represent the 16S rRNA gene and the 23S/5S rRNA gene operon, respectively. Circle 4: GC skew. The diagrams were generated using Artemis 12.0 - DNAPlotter version 1.4, Sanger Institute. (M. suis plot was extracted from Guimaraes AM, Santos AP, SanMiguel P, Walter T, Timenetsky J, Messick JB: Complete genome sequence of Mycoplasma suis and insights into its biology and adaption to an erythrocyte niche. PLoS One 2011, 6:e19574 [1], with permission from the copyright holder).

The set of unique CDS of M. parvum and M. suis were analyzed for subcellular localization and protein sorting signals [22],[23], and the presence of signal peptide cleavage sites [24]. Moreover, the presence of lipoproteins was predicted [25] (TableĀ 3). The majority (54.0-66.5%) of the unique CDS either in paralogous families or not from both M. parvum and M. suis have an unknown subcellular localization. The remaining CDS are distributed into cytoplasmic membrane, cytoplasmic, and extracellular localization. The unique CDS not in paralogous families from both bacteria are divided mainly into cytoplasmic membrane and cytoplasmic localization, with a smaller number of CDS having a predicted extracellular localization. In contrast, the unique CDS in paralogous families are mostly located in the cytosol, and just a few CDS have cytoplasmic membrane or extracellular localization. The most frequent protein sorting signal amongst all unique CDS is the presence of 1 internal helix, followed by a signal peptide. Analyses of the CDS in paralogous families (details in the Additional file 2) show that a higher number of CDS of M. suis (16.5% by PSORTb, and 4.25% by SignalP) have signal peptides compared to M. parvum (1.6% bt PSORTb, and 0% by SignalP). Only three lipoproteins (SpII) were identified using LipoP: one unique CDS not paralog of M. parvum, one unique CDS not paralog of M. suis, and one unique CDS paralog of M. suis.

Table 3 Unique CDS of M. parvum strain Indiana and M. suis strain Illinois: subcellular localization and protein sorting signals (PSORTb v3.0.2), presence of signal peptide cleavage sites (SignalP), and presence of lipoproteins (LipoP).

Discussion

The course of infection with M. parvum strain Indiana in the single splenectomized pig evaluated was distinguished from what is commonly seen in M. suis infection by the absence of clinical signs, even at the peak of bacteremia. While the bacteremia was apparently fleeting with organisms detected on peripheral blood smears for only a few days after the peak, its persistence at low levels were shown by qPCR. However, these observations are limited to the strain of M. parvum described herein and to the animal selected for this study. Further studies should be conducted to evaluate the actual virulence of this isolate, including possible variation among strains and individual pigs. Since the prevalence of M. parvum by age is unknown, it is not possible to speculate when this animal got infected with this organism. It is likely, however, that an immune response had developed, which effectively controlled the infection, reducing the bacterial loads. Another possibility is that the use of subtherapeutic doses of antibiotics in the feed at the farm may have controlled the infection. However, the bacteria were not completely eliminated and the animal developed a chronic infection, as is often observed with other hemoplasma species [30]. In contrast, the fever and bacteremia in a splenectomized pig that is infected with M. suis is unrelenting and without antibiotics the animal may die [31]ā€“[33], which was not observed after splenectomy in the animal herein. The impact of acute and chronic hemoplasma infection on the immune system of the host has been poorly explored [34], and is an ongoing area of investigation in our laboratory.

In this report, M. parvum strain Indiana was described morphologically by light and scanning electron microscopy and its genome was completely sequenced, analyzed and compared to the genome of M. suis strain Illinois. A transient bacteremia was demonstrated by light microscopy, whereas SEM confirmed the rod and coccoid morphology, epicellular location, and small size of M. parvum (FigureĀ 2B) [11]. Its size, 0.2 to 0.5Ā Ī¼m, was similar to that reported for other hemoplasmas [35], but remarkably smaller than that reported for M. suis with ring forms approaching 1.0 to 2.5Ā Ī¼m in diameter [11],[36]ā€“[38]. The level of bacteremia for M. parvum at its peak as shown by qPCR was considerably less (one log) than that previously reported for M. suis[16].

M. parvum has the smallest single, circular chromosome of all the Mycoplasma genomes sequenced to date. The characteristics of the M. parvum genome, including its small size, low G + C content, use of UGA codon to encode tryptophan, and number of tRNA and rRNA were in agreement with those reported for the genomes of other hemoplasmas and are typical of mycoplasmas [1],[4],[6]. In addition, the percentage of CDS dedicated to paralogous genes (24.3%) was similar to that reported for other hemoplasmas [1]ā€“[8],[10]. Our group and others have hypothesized that archived sequences of these genes distributed throughout the chromosome play a role in the ability of mycoplasma organisms to persist despite an active host immune response [39].

Phylogenetic studies of M. parvum and M. suis, based on sequence analyses of the 16S rRNA and RNase P genes, suggest that these bacteria are closely related [13],[14]. Some authors have reported that M. parvum is an immature developmental stage of M. suis that is present concurrently with the mature pathogen [33],[37]. A phylogenomic approach based on ANI and tetranucleotide signatures that uses whole genome sequence information to compare different organisms provides a better resolution than single gene sequence approaches. It is comparable to the method considered as the gold standard for prokaryotic species definition, the DDH [15]. ANI and tetranucleotide results in this study were below the cut-off values for species definition separating all the organisms included in the analyses as distinct species of hemotrophic mycoplasmas. M. parvum is, indeed, a distinct species infecting the pig. Further, strains of the same species were clearly distinguished by ANI and tetra scores above the proposed threshold.

M. parvum and M. suis, share the same host, the pig; the same environment, the blood; and are phylogenetically, the closest related relatives of one another [13],[14]. However, the interaction of these bacteria with their host is completely distinct. M. suis may cause life threatening hemolytic anemia during acute infection [30],[40], while the strain of M. parvum shown in this study, and in agreement with previous reports about this organism [11], caused no clinical signs even at the peak of bacteremia in a young, splenectomized pig. At the genomic level, these bacteria are remarkably similar: they share all CDS with known functions. Their genomic synteny, however, is not well conserved. This loss of conservation in the gene order is expected, as different species of hemoplasmas rarely share large syntenic regions, with the exception of M. haemofelis and M. haemocanis[41]. Nevertheless, the genome of M. parvum has orthologous for all the CDS with metabolic functions identified in the genome of M. suis, implying that their metabolic pathways work very similarly and they have the same requirements for many nutrients to be acquired from the blood environment [1]. The putative virulence factors described for the M. suis genome are also present in M. parvum, which raises the question about what is different between these hemoplasmas that reflects in their contrasting pathogenicity. A gene-by-gene comparison between these organisms revealed that their only difference at the genomic level relies in the hypothetical CDS, mostly dedicated to paralogous gene families. M. suis has 40 paralogous families that are exclusive to its genome, while M. parvum has 15 families that are not found in M. suis (Additional file 2). In addition, the largest paralogous family (37 CDS) of M. suis is unique. Interestingly, the set of unique paralogous CDS of M. suis have a higher percentage of signal peptides detected compared to M. parvum CDS; this might indicate that these CDS are secreted or inserted into the membrane and therefore they could play an important role in the pathogenicity of M. suis as reported for other bacteria [42],[43]. The differences regarding the percentages of signal peptides detected by all three software (PSORTb, SignalP, and LipoP) are due to their distinct prediction algorithms. PSORTb is the only tool amongst the three applied in this study that provides an output for Gram-negative bacteria with no outer membrane or cell wall. It is presumed more accurate for Mycoplasma species, which although more phylogenetically similar to Gram-positive organisms, lack a peptidoglycan cell wall.

The maintenance of paralogous gene families in their genomes appears to be a high priority among the hemoplasmas. A prominent role for paralogous genes in the development of antigenic variation has been reported for bacterial species of the genus Anaplasma, Ehrlichia, Borrelia and Mycoplasma and linked to virulence as well as persistence of these infections [44]. It is tempting to speculate that differences in paralogous families between M. parvum and M. suis may be implicated in their differing pathogenic potential. However, this might also be due to variable expression of core genes, especially genes encoding putative virulence factors or to the presence and variable expression of regulatory RNA. The latter has been previously reported for pathogenic versus non-pathogenic bacteria species of Listeria[45], while variable expression of virulence genes was observed for different strains of Escherichia coli[46], and Listeria monocytogenes[47]. An investigation on the differences at the RNA level between M. parvum and M. suis is warranted, and it is an ongoing project in our laboratory.

Nucleotide sequence accession number

The genome of M. parvum strain Indiana was deposited in GenBank under the accession number [CP006771.1].

Additional files

References

  1. Guimaraes AM, Santos AP, SanMiguel P, Walter T, Timenetsky J, Messick JB: Complete genome sequence of Mycoplasma suis and insights into its biology and adaption to an erythrocyte niche. PLoS One. 2011, 6: e19574-10.1371/journal.pone.0019574.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  2. Oehlerking J, Kube M, Felder KM, Matter D, Wittenbrink MM, Schwarzenbach S, Kramer MM, Hoelzle K, Hoelzle LE: Complete genome sequence of the hemotrophic Mycoplasma suis strain KI3806. J Bacteriol. 2011, 193: 2369-2370. 10.1128/JB.00187-11.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  3. Barker EN, Helps CR, Peters IR, Darby AC, Radford AD, Tasker S: Complete genome sequence of Mycoplasma haemofelis, a hemotropic mycoplasma. J Bacteriol. 2011, 193: 2060-2061. 10.1128/JB.00076-11.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  4. Santos AP, Guimaraes AMS, Do Nascimento NC, Sanmiguel PJ, Martin SW, Messick JB: Genome of Mycoplasma haemofelis, unraveling its strategies for survival and persistence. Vet Res. 2011, 42: 102-10.1186/1297-9716-42-102.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  5. Barker EN, Darby AC, Helps CR, Peters IR, Hughes MA, Radford AD, Novacco M, Boretti FS, Hofmann-Lehmann R, Tasker S: Genome sequence for ā€œCandidatus Mycoplasma haemominutum,ā€ a low-pathogenicity hemoplasma species. J Bacteriol. 2012, 194: 905-906. 10.1128/JB.06560-11.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  6. Do Nascimento NC, Santos AP, Guimaraes AM, Sanmiguel PJ, Messick JB:Mycoplasma haemocanis - the canine hemoplasma and its feline counterpart in the genomic era. Vet Res. 2012, 43: 66-10.1186/1297-9716-43-66.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  7. Dos Santos AP, Guimaraes AM, Do Nascimento NC, SanMiguel PJ, Messick JB: Complete genome sequence of Mycoplasma wenyonii strain Massachusetts. J Bacteriol. 2012, 194: 5458-5459. 10.1128/JB.01240-12.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  8. Guimaraes AM, Toth B, Santos AP, Do Nascimento NC, Kritchevsky JE, Messick JB: Genome sequence of ā€œCandidatus Mycoplasma haemolamaeā€ strain purdue, a red blood cell pathogen of alpacas (Vicugna pacos) and llamas (Lama glama). J Bacteriol. 2012, 194: 6312-6313. 10.1128/JB.01557-12.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  9. Do Nascimento NC, Dos Santos AP, Chu Y, Guimaraes AM, Pagliaro A, Messick JB: Genome sequence of Mycoplasma parvum (formerly Eperythrozoon parvum), a diminutive hemoplasma of the pig. Genome Announc. 2013, 1: pii:e00986-13. 10.1128/genomeA.00986-13.

    ArticleĀ  Google ScholarĀ 

  10. Deshuillers PL, Santos AP, Do Nascimento NC, Hampel JA, Bergin IL, Dyson MC, Messick JB: Complete genome sequence of Mycoplasma ovis strain Michigan, a hemoplasma of sheep with two distinct 16S rRNA genes. Genome Announc. 2014, 2: e01235-13. 10.1128/genomeA.01235-13.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  11. Splitter EJ:Eperythrozoon suis n. sp. and Eperythrozoon parvum n. sp., 2 new blood parasites of swine. Science. 1950, 111: 513-514. 10.1126/science.111.2889.513.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  12. Seamer J: Studies with Eperythrozoon parvum Splitter, 1950. Parasitology. 1960, 50: 67-80. 10.1017/S0031182000025208.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  13. Zhou RQ, Nie K, Huang HC, Hu SJ, Zhou ZY, Luo HL: Phylogenetic analysis of Mycoplasma suis isolates based on 16S rRNA gene sequence in China. Vet Res Commun. 2009, 33: 855-863. 10.1007/s11259-009-9234-3.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  14. Watanabe Y, Fujihara M, Obara H, Nagai K, Harasawa R: Two genetic clusters in swine hemoplasmas revealed by analyses of the 16S rRNA and RNase P RNA genes. J Vet Med Sci. 2011, 73: 1657-1661. 10.1292/jvms.11-0293.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  15. Richter M, RossellĆ³-MĆ³ra R: Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A. 2009, 106: 19126-19131. 10.1073/pnas.0906412106.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  16. Guimaraes AMS, Vieira RFC, Polleto R, Vemulapalli R, Santos AP, De Moraes W, Cubas ZS, Santos LC, Marchant-Forde JN, Timenetsky J, Biondo AW, Messick JB: A quantitative TaqMan PCR assay for the detection of Mycoplasma suis. J Appl Microbiol. 2011, 111: 417-425. 10.1111/j.1365-2672.2011.05053.x.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  17. Messick JB, Cooper SK, Huntley M: Development and evaluation of a polymerase chain reaction assay using the 16S rRNA gene for the detection or Eperythrozoon suis infection. J Vet Diagn Invest. 1999, 11: 229-236. 10.1177/104063879901100304.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  18. Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJ, Birol I: ABySS: a parallel assembler for short read sequence data. Genome Res. 2009, 19: 1117-1123. 10.1101/gr.089532.108.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  19. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999, 27: 29-34. 10.1093/nar/27.1.29.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  20. Yus E, Maier T, Michalodimitrakis K, Van Noort V, Yamada T, Chen WH, Wodke JA, GĆ¼ell M, MartĆ­nez S, Bourgeois R, KĆ¼hner S, Raineri E, Letunic I, Kalinina OV, Rode M, Herrmann R, GutiĆ©rrez-Gallego R, Russell RB, Gavin AC, Bork P, Serrano L: Impact of genome reduction on bacterial metabolism and its regulation. Science. 2009, 326: 1263-1268. 10.1126/science.1177263.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  21. Biegert A, Mayer C, Remmert M, Sƶding J, Lupas AN: The MPI bioinformatics toolkit for protein sequence analysis. Nucleic Acids Res. 2006, 34: W335-W339. 10.1093/nar/gkl217.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  22. Nakai K, Horton P: PSORT: a program for detecting sorting signals in proteins and predicting their subcellular localization. Trends Biochem Sci. 1999, 24: 34-36. 10.1016/S0968-0004(98)01336-X.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  23. Yu NY, Wagner JR, Laird MR, Melli G, Rey S, Lo R, Dao P, Sahinalp SC, Ester M, Foster LJ, Brinkman FS: PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics. 2010, 26: 1608-1615. 10.1093/bioinformatics/btq249.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  24. Petersen TN, Brunak S, Von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8: 785-786. 10.1038/nmeth.1701.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  25. Juncker AS, Willenbrock H, Von Heijne G, Brunak S, Nielsen H, Krogh A: Prediction of lipoprotein signal peptides in Gram-negative bacteria. Protein Sci. 2003, 12: 1652-1662. 10.1110/ps.0303703.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  26. CoGePedia [https://www.genomevolution.org/wiki/index.php/Main_Page]

  27. SynMap whole genome synteny [https://genomevolution.org/CoGe/SynMap.pl]

  28. Robertshaw D: Temperature regulation and thermal environment. Dukesā€™ Physiology of Domestic Animals. Edited by: Reece WO. 2004, Cornell University Press, New York, 12

    Google ScholarĀ 

  29. Thorn CE: Normal hematology of the pig. Schalmā€™s Veterinary Hematology. Edited by: Feldman BF, Zinkl JG, Jain NC. 2000, Lippincott Williams and Wilkins, Philadelphia, 1089-1095. 5

    Google ScholarĀ 

  30. Messick JB: Hemotrophic mycoplasmas (hemoplasmas): a review and new insights into pathogenic potential. Vet Clin Pathol. 2004, 33: 2-13. 10.1111/j.1939-165X.2004.tb00342.x.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  31. Heinritzi K, Peteranderl W, Plank G:Eperythrozoon infection in swine: effect on the acidā€“base balance and the glucose, lactate and pyruvate content of venous blood. Dtsch Tierarztl Wochenschr. 1990, 97: 31-34. (in German),

    CASĀ  PubMedĀ  Google ScholarĀ 

  32. Heinritzi K, Plank G: The effect of Eperythrozoon suis infection on the osmotic fragility of erythrocytes. Berl Munch Tierarztl Wochenschr. 1992, 105: 380-383.

    CASĀ  PubMedĀ  Google ScholarĀ 

  33. Dent BT, Stevens KA, Korvick DL, Clymer JW:Mycoplasma suis infection in pigs after splenectomy. Lab Anim (NY). 2013, 42: 125-128. 10.1038/laban.171. (in German),

    ArticleĀ  Google ScholarĀ 

  34. Novacco M, Boretti FS, Franchini M, Riond B, Meli ML, Hofmann-Lehmann R: Protection from reinfection in ā€œCandidatus Mycoplasma turicensisā€-infected cats and characterization of the immune response. Vet Res. 2012, 43: 82-10.1186/1297-9716-43-82.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  35. Neimark H, Hoff B, Ganter M:Mycoplasma ovis comb. nov. (formerly Eperythrozoon ovis), an epierythrocytic agent of haemolytic anaemia in sheep and goats. Int J Syst Evol Microbiol. 2004, 54: 365-371. 10.1099/ijs.0.02858-0.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  36. Pospischil A, Hoffmann R:Eperythrozoon suis in naturally infected pigs: a light and electron microscopic study. Vet Pathol. 1982, 19: 651-657. 10.1177/030098588201900609.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  37. Liebich HG, Heinritzi K: Light and electron microscopic studies of Eperythrozoon suis. Tierarztl Prax. 1992, 20: 270-274. (in German),

    CASĀ  PubMedĀ  Google ScholarĀ 

  38. Groebel K, Hoelzle K, Wittenbrink MM, Ziegler U, Hoelzle LE:Mycoplasma suis invades porcine erythrocytes. Infect Immun. 2009, 77: 576-584. 10.1128/IAI.00773-08.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  39. Citti C, Nouvel LX, Baranowski E: Phase and antigenic variation in mycoplasmas. Future Microbiol. 2010, 5: 1073-1085. 10.2217/fmb.10.71.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  40. Strait EL, Hawkins PA, Wilson WD: Dysgalactia associated with Mycoplasma suis infection in a sow herd. J Am Vet Med Assoc. 2012, 241: 1666-1667. 10.2460/javma.241.12.1666.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  41. Guimaraes AM, Santos AP, Do Nascimento NC, Timenetsky J, Messick JB: Comparative genomics and phylogenomics of hemotrophic mycoplasmas. PLoS ONE. 2014, 9: e91445-10.1371/journal.pone.0091445.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  42. Saleh MT, Fillon M, Brennan PJ, Belisle JT: Identification of putative exported/secreted proteins in prokaryotic proteomes. Gene. 2001, 269: 195-204. 10.1016/S0378-1119(01)00436-X.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  43. Leversen NA, De Souza GA, MĆ„len H, Prasad S, Jonassen I, Wiker HG: Evaluation of signal peptide prediction algorithms for identification of mycobacterial signal peptides using sequence data from proteomic methods. Microbiology. 2009, 155: 2375-2383. 10.1099/mic.0.025270-0.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  44. Palmer GH: The highest priority: what microbial genomes are telling us about immunity. Vet Immunol Immunopathol. 2002, 85: 1-8. 10.1016/S0165-2427(01)00415-9.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  45. Wurtzel O, Sesto N, Mellin JR, Karunker I, Edelheit S, BĆ©cavin C, Archambaud C, Cossart P, Sorek R: Comparative transcriptomics of pathogenic and non-pathogenic Listeria species. Mol Syst Biol. 2012, 8: 583-10.1038/msb.2012.11.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  46. Dowd SE, Ishizaki H: Microarray based comparison of two Escherichia coli O157:H7 lineages. BMC Microbiol. 2006, 6: 30-10.1186/1471-2180-6-30.

    ArticleĀ  PubMed CentralĀ  PubMedĀ  Google ScholarĀ 

  47. Severino P, Dussurget O, VĆŖncio RZ, Dumas E, Garrido P, Padilla G, Piveteau P, LemaĆ®tre JP, Kunst F, Glaser P, Buchrieser C: Comparative transcriptome analysis of Listeria monocytogenes strains of the two major lineages reveals differences in virulence, cell wall, and stress response. Appl Environ Microbiol. 2007, 73: 6078-6088. 10.1128/AEM.02730-06.

    ArticleĀ  PubMed CentralĀ  CASĀ  PubMedĀ  Google ScholarĀ 

Download references

Acknowledgements

We are extremely grateful to The Genomic Core Facility at Purdue University for constructing the library, running the HiSeq 2000 and performing the ABySS assembly. We also thank the University of Marylandā€™s Institute of Genomic Sciences, which provided us with the manual annotation pipeline Manatee. This work was kindly supported by the Morris Animal Foundation (grant D10FE-0042), and the Brazilian government through CNPq Fellowship Program (Conselho Nacional de Desenvolvimento CientĆ­fico e TecnolĆ³gico) for Postdoctoral Researcher to Dr do Nascimento.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to NaĆ­la C do Nascimento or Joanne B Messick.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authorsā€™ contributions

Conceived and designed experiments: NCN, APS, YC, AMSG, JBM. Performed experiments: NCN, APS, YC, AMSG, JBM. Analyzed the data: NCN, APS, YC, AMSG, JBM. Anesthesia and splenectomy: ANB, ABW. Wrote and reviewed the paper: NCN, APS, YC, AMSG, ANB, ABW, JBM. All authors read and approved the final manuscript.

Electronic supplementary material

13567_2014_86_MOESM1_ESM.pdf

Additional file 1: Syntenic map between Mycoplasma suis strain Illinois and M. parvum strain Indiana. Plots were generated using comparative genomics suite CoGe SynMap (ā€œlastā€ analysis). Each dot represents a matching gene pair. Figure shows the syntenic map (organization of the genes) between Mycoplasma suis strain Illinois and M. parvum strain Indiana. (PDF 38 KB)

13567_2014_86_MOESM2_ESM.pdf

Additional file 2: Unique CDS in paralogous families of M. parvum strain Indiana and M. suis strain Illinois: subcellular localization and protein sorting signals (PSORTb v3.0.2), presence of signal peptide cleavage sites (SignalP), and presence of lipoproteins (LipoP). Table shows the numbers of unique CDS of Mycoplasma parvum strain Indiana and M. suis strain Illinois organized in paralogous families. Unique CDS were analyzed for: subcellular localization and protein sorting signals according to PSORTb v3.0.2, presence of signal peptide cleavage sites using SignalP 4.1, and presence of lipoproteins using LipoP version 1.0 software. (PDF 137 KB)

Authorsā€™ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

do Nascimento, N.C., dos Santos, A.P., Chu, Y. et al. Microscopy and genomic analysis of Mycoplasma parvum strain Indiana. Vet Res 45, 86 (2014). https://doi.org/10.1186/s13567-014-0086-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13567-014-0086-7

Keywords