|Home | About | Journals | Submit | Contact Us | Français|
This work reports the results of analyses of three complete mycoplasma genomes, a pathogenic (7448) and a nonpathogenic (J) strain of the swine pathogen Mycoplasma hyopneumoniae and a strain of the avian pathogen Mycoplasma synoviae; the genome sizes of the three strains were 920,079 bp, 897,405 bp, and 799,476 bp, respectively. These genomes were compared with other sequenced mycoplasma genomes reported in the literature to examine several aspects of mycoplasma evolution. Strain-specific regions, including integrative and conjugal elements, and genome rearrangements and alterations in adhesin sequences were observed in the M. hyopneumoniae strains, and all of these were potentially related to pathogenicity. Genomic comparisons revealed that reduction in genome size implied loss of redundant metabolic pathways, with maintenance of alternative routes in different species. Horizontal gene transfer was consistently observed between M. synoviae and Mycoplasma gallisepticum. Our analyses indicated a likely transfer event of hemagglutinin-coding DNA sequences from M. gallisepticum to M. synoviae.
Mycoplasmas comprise a group of more than 180 species of wall-less bacteria that are obligate parasites of a wide range of organisms including humans, plants, and animals (46). Mycoplasmas typically exhibit strict host tissue specificities, probably due to their nutritional requirements (45), a direct consequence of the genome reduction that likely occurred as a consequence of the metabolic complementarity of their hosts (3). The evolutionary dynamics of these organisms involved population bottlenecks and asexual reproduction leading to accumulation of deleterious mutations, which resulted in further genome contraction (58). A predictable consequence of this process is preservation of a minimal genome comprising essential genes to maintain basic core functions and adaptation to specific environments.
Two species, Mycoplasma hyopneumoniae and Mycoplasma synoviae, have a significant adverse economic impact on animal production. The former is the infective agent of enzootic pneumonia in pigs, which results in deactivation of mucociliary functions (15) and increased susceptibility to secondary infections (12). The latter is responsible for respiratory tract disease and synovitis in chickens and turkeys. It can be transmitted vertically through contaminated eggs (28), resulting in considerable losses due to reduced egg production and meat quality as well as a lowered rate of viable hatchings. Thus, knowledge of their respective biological characteristics seems of paramount importance.
The genomes of several mycoplasmas have been sequenced and analyzed in recent years (11, 22, 25-27, 35, 44, 50, 61), but comparative analyses of species belonging to the Pneumoniae and Hominis clades have not been undertaken. In addition, interstrain, whole-genome comparisons have not yet been carried out, although the genes involved in DNA repair, including those of the organisms herein studied, have recently been analyzed (10).
Here we report the complete genome sequences of a pathogenic (7448) and a nonpathogenic (J [ATCC 25934]) strain of M. hyopneumoniae and the complete genome of M. synoviae strain 53. Comparative analyses of the M. hyopneumoniae strains allowed the identification of strain-specific regions that might be related to their variable pathogenicity. A detailed phylogenetic analysis of several mycoplasma species belonging to the Pneumoniae and Hominis clades was also carried out, comparing metabolic pathways and genes involved in the adhesion process. Comparisons of M. gallisepticum and M. synoviae genomes pointed to the evolutionary origin of the hemagglutinin gene family and showed evidence of horizontal transfer of other gene clusters.
M. hyopneumoniae strain J (ATCC 25934) was acquired from American Type Culture Collection by CNPSA, EMBRAPA (Concórdia, Santa Catarina, Brazil). This is a nonpathogenic strain with a reduced adhesion capacity to porcine cilia (62-64). M. hyopneumoniae strain 7448 was isolated from an infected swine in Lindóia do Sul, Santa Catarina, Brazil. Specific-pathogen-free pigs inoculated with strain 7448 consistently produced the characteristic symptoms of enzootic pneumonia. M. synoviae was isolated from a broiler breeder in the state of Paraná in Brazil (19).
Genomes were sequenced using the shotgun sequencing strategy (20). Sequencing, assembling, annotation, and comparative in silico analyses were carried out by the Brazilian National Genome Sequencing Consortium and the Southern Network for Genome Analysis (PIGS), involving a total of 28 sequencing laboratories, one bioinformatics center, and three coordinating laboratories. Template preparation was performed using standard protocols. DNA sequencing reactions were performed using the DYEnamic ET dye terminator cycle sequencing (MegaBACE) kit and run on MegaBACE 1000 capillary sequencers (Amersham Biosciences). Approximately 10,000 reads per genome with phred scores of >20 were generated from both ends of plasmid clones ranging from 2.0 to 4.0 kb, providing an approximately 13-fold genome coverage. Sequences were assembled using phred/phrap/consed (http://www.phrap.org). Sequencing gaps were closed using the information generated by autofinisher, while our recently developed strategy of PCR-assisted contig extension (PACE) (8) was used for physical gap closure. Annotation was carried out using the System for Automated Bacterial Integrated Annotation (SABIÁ) (2), developed to integrate public-domain and purpose-built software for the automated identification of genome landmarks, including tRNA and rRNA sequences, repetitive elements, and coding DNA sequences (CDSs) (which indicate regions likely to encode proteins). Paralogous gene families were defined using a cutoff E value of 10−5 with at least 60% query coverage and 50% identity.
Maximum likelihood (ML) phylogenies, based on individual orthologous proteins, were generated using ProtML (Molphy package [http://www.ism.ac.jp/ismlib/softother.e.html]) and TREE-PUZZLE 5.1 (51). A data set concatenating all proteins in a single sequence unit was analyzed using neighbor-joining (NJ) distance trees with MEGA 2.1 (30), maximum parsimony using PAUP* 4.0b10 (55), and by ML using ProtML, all with confidence estimates based on 100 bootstrap replicates. A bootstrap gene tree was calculated following 500 random, protein resamplings and concatenation, with subsequent analysis by NJ based on ML distances using the Molphy package. Divergence times were estimated by the linearized tree method using MEGA 2.1 and r8s 1.6 (http://ginger.ucdavis.edu/r8s), assuming 450 million years before the present (MYBP) as the time of divergence of the phytoplasmas from mycoplasmas (32). Orthologous clusters were identified using the bidirectional best hit method (43). Clustering of hemagglutinin CDSs was performed with Tribe-MCL (18), based on data from allXall National Center for Biotechnology Information (NCBI) BLASTp searches. Global genome alignments were carried out using Mauve (13). Genome rearrangements between M. hyopneumoniae strains J, 7448, and 232 were identified by combined analyses with GRIL (14) and Artemis (49). Genome duplications were inferred from self-BLAST searches. Visualization of local similarities between CDSs of complete mollicute genomes was carried out with PhyloGrapher (www.atgc.org/PhyloGrapher). Horizontal gene transfer (HGT) was initially detected with allXall BLAST searches in available mollicute genomes, supplemented by a compositional and codon bias scan. Further alignment was carried out with Mauve for genomes with best hits in the initial search. Phylogenetic distances between regions sharing at least 300 nucleotides, with Mauve alignment, were identified in each genome and were compared against distance estimates based on 16S rRNA sequence data. An HGT event was considered plausible when the estimated distance of the aligned region was lower than the estimated 16S rRNA distance.
All intergenic regions (IGRs) were extracted from the genome data of Mycoplasma hyopneumoniae strains J and 7448 and compared using MUMmer. The search for putative regulatory signals upstream of M. hyopneumoniae CDSs employed two different methodologies. Initially, a clusterization of the upstream region of M. hyopneumoniae genes was carried out (first 50 bases upstream of the translation starting point), on the basis of similarity, and employing the BLASTCLUST software (ftp://ftp.ncbi.nlm.nih.gov/blast/). Sequences were subsequently submitted to analysis with GLAM software (http://zlab.bu.edu/glam/), aiming to find conserved patterns among the initially clusterized sequences. A second strategy involved search of fuzzy motifs with the Self-Organizing Map, a neural-network algorithm generated by SOMBRERO software (http://bioinf.nuigalway.ie/sombrero/index.html).
Sequence data reported in this paper were deposited in GenBank (accession nos. AE017243, AE017244, and AE017245). Sequence and annotation data are available at http://www.brgene.lncc.br/finalMS (M. synoviae), http://www.genesul.lncc.br/finalMH (M. hyopneumoniae strain J), and http://www.genesul.lncc.br/finalMP (M. hyopneumoniae strain 7448).
The main features of the three newly sequenced genomes are shown in Table Table1.1. Clustering analysis of M. synoviae, M. hyopneumoniae strains J, 7448, and 232 (35), and eight other mycoplasma genomes (Mycoplasma pneumoniae M129, Mycoplasma pulmonis UAB CTIP, Mycoplasma penetrans HF-2, Mycoplasma genitalium G37, Mycoplasma gallisepticum R, Mycoplasma mycoides subsp. mycoides SC PG1, Mycoplasma mobile 163K, and Ureaplasma urealyticum serovar 3) revealed 235 orthologous clusters. A comparison of the number of CDSs in all sequenced mycoplasma genomes is shown in Table Table22.
Comparison of the three M. hyopneumoniae strains provided evidence of intraspecific rearrangements, resulting in strain-specific gene clusters (Fig. (Fig.1).1). This was the case for a 16-kb region of M. hyopneumoniae strain J, containing 15 CDSs, most of which encoded type III restriction-modification (R-M) system components and putative transposases. M. hyopneumoniae strain 7448 contained a specific 22.3-kb region similar to the integrative conjugal element (ICEF) of Mycoplasma fermentans (7), which was designated ICEH (for integrative conjugal element of M. hyopneumoniae). ICEH contained 14 CDSs, four of which similar to tra genes, usually associated with bacterial conjugative plasmids, and another encoding a single-strand binding protein (SSB), an essential protein for the transfer process. Direct repeat sequences (TAGATTTTT), generated by target site duplications, flanked ICEH. This target site was localized in the homologous region of M. hyopneumoniae strain J, pointing to the mobility of this element. Evidences for the presence of circular extrachromosomal forms of ICEH in M. hyopneumoniae strain 7448 and another unrelated, pathogenic Brazilian field isolate of M. hyopneumoniae were obtained by an inverted PCR assay (results not shown), indicating that this element might be functionally active in these isolates. Moreover, we also observed a region similar to ICEF in the genome of M. hyopneumoniae strain 232. It has recently been demonstrated that some pathogenic bacteria use the type IV secretion system, composed of subunits related to the conjugation machinery, for the delivery of effector molecules to host cells (16), and that this system may be involved in pathogenesis (52). However, the involvement of ICEH in pathogenesis through delivery of effector molecules into cells remains to be explored. Other strain-specific differences included an inverted region of 243,104 bp in M. hyopneumoniae strain 232 (Fig. (Fig.1)1) and less drastic rearrangements between M. hyopneumoniae strains at five genomic regions (Fig. (Fig.1;1; see Tables S1 to S3 in the supplemental material). The positions of regions 1, 2, 4, and 5 involved in rearrangements are conserved in all three strains. Region 3 is syntenic in M. hyopneumoniae strains J and 232. In M. hyopneumoniae strain J, ABC transporter-encoding region 1 lacks two CDSs that are present in M. hyopneumoniae strains 7448 (MHP0023 and MHP0024) and 232 (mhp025 and mhp026), which probably originated by duplication followed by divergence. Region 2 is characterized, in the three strains, by the presence of several unique CDSs. This region also contains short translocations and duplications. Region 3 shows short unique insertions in M. hyopneumoniae strain 7448 in comparison to strain J or 232 and a duplication of approximately 2.2 kb containing a hypothetical CDS and a serine protease-encoding CDS. Probably, this duplication was followed by deletion of part of the serine protease sequence. A third rearrangement involved a translocation that displaced the DNA segment containing the complete or partial serine protease-encoding CDS by approximately 300 kb in the two strains. This region was syntenic in both M. hyopneumoniae strains J and 232. Region 4, in the three genomes, shows short unique insertions and translocations of CDS-containing segments, probably mediated by a transposable element, as indicated by the presence of an ISMhp1 transposon-like element in this region. Region 5 presents a series of imperfect duplications of hypothetical CDSs, generating seven copies in M. hyopneumoniae strain J, five in strain 7448, and up to four copies in strain 232.
A comparison of 362 orthologous IGRs of M. hyopneumoniae strains showed that 21 to 29% of them shared identical sequences, 43 to 45% differed by 1 to 5 bp, 13 to 15% by 6 to 10 bp, and 13 to 20% by more than 10 bp. Moreover, in 34% of them, repeated sequences varied in length (see Table S4 in the supplemental material), as was the case for CDSs encoding DnaJ, thymidine phosphorylase, RpoB, a p97-like protein, methyltransferase, serine protease, and p65. This variation has been reported in both intergenic and coding regions in several genomes (31, 33, 47) and is associated with polymerase slippage during replication (57). The intergenic regions identified in the M. hyopneumoniae genomes present a higher A+T content (about 80%) than coding sequences (about 70%). It has been recently shown, in 152 genomes, that AT content is higher in the 200 bp upstream of translation start codons than in the 200 bp downstream (56). Only a few promoters have been identified and analyzed in mycoplasmas (36, 59, 60), and these studies identified strong consensus sequences at −10 regions but only weak consensus sequences at −35 regions. We searched for putative regulatory signals upstream of M. hyopneumoniae CDSs and found 17 significant clusters, representing only a small fraction of the CDSs of the M. hyopneumoniae genome (70 CDSs). The physiological relevance of genomic rearrangements, nucleotide substitutions, and changes in the lengths of IGRs in M. hyopneumoniae remains to be elucidated.
Genomic comparisons between pathogenic (M. hyopneumoniae strains 7448 and 232) and nonpathogenic (M. hyopneumoniae strain J) strains revealed important aspects related to pathogenicity. Among previously described adhesion-related proteins (29, 48) (Table (Table3),3), p97 is regarded as the major cell adhesion determinant, although other proteins, derived from either a mycoplasma or host, might also participate in membrane anchorage (17). In the three strains, the p97 CDS is linked to p102 CDS, comprising a two-CDS operon (p97-p102 operon I), as well as two other operons that are clearly related (p97-p102 operons II and III) that exhibit >80% identity between orthologous deduced amino acid sequences. The p97 orthologous CDSs in the p97-p102 operon I of M. hyopneumoniae strains 7448, 232, and J code for proteins with 10, 15, and 9 of the previously described R1 tandem repeats, respectively; all three strains had more than the minimum number of CDSs (8 CDSs) required for cilium binding (34). This indicates that other adhesion determinants, with different characteristics in M. hyopneumoniae strains 7448, 232, and J, must be responsible for their different adhesion properties. Differences between the encoded adhesins and other putative virulence factors in M. hyopneumoniae strains 7448, 232, and J include variations in the number of amino acid repeats between orthologous proteins (Table (Table4).4). These insertions/deletions result from variations in the number of tandem nucleotide repeats within coding regions, which is indicative of a molecular mechanism generating functional and/or antigenic variants. This variation in surface proteins is likely to be a key determinant of different pathogenic properties of each M. hyopneumoniae strain.
Phylogenetic relationships were established on the basis of a maximum of 206 single-copy CDSs (comprising ~86,000 aligned deduced amino acid positions). The results of phylogenetic analyses based on individual orthologous proteins were subsequently compared to the ML, maximum parsimony, and distance topologies based on the concatenated amino acid sequence alignment of a single sequence unit. In contrast to the highly incongruent topologies resulting from separate analyses of individual orthologous CDSs, the concatenated alignment generated a single, highly supported tree (Fig. (Fig.2A).2A). Phylogenies were also estimated with nucleotide data using a wide variety of methods, based on either gene content or gene order. Only the tree generated with the concatenated protein sequence unit was presented, as all methods gave essentially the same results. A tentative time frame of genome evolution was estimated (Fig. 2B and C), indicating that the Mycoplasmatales and Entomoplasmatales orders split into three clades between 600 to 400 MYBP, while most of the species diverged some 400 to 300 MYBP, approximately at the time of the emergence and diversification of tetrapods (4). M. pneumoniae and M. genitalium diverged more recently, between 75 to 35 MYBP, and M. hyopneumoniae strains J and 7448 diverged about 2 MYBP.
We analyzed the possibility of HGT in view of its putative role in determining the characteristics of prokaryotic genomes (41) and the transfer of genes related to pathogenesis (21). A comparison of mollicute genomes using both parametric and phylogenetic strategies showed that HGT was most likely to have occurred between M. synoviae and M. gallisepticum (Fig. (Fig.3).3). Fourteen putative transferred regions were identified (Table (Table5),5), the largest comprising 5.9 kb and encompassing not only hypothetical CDSs but also CDSs coding for an ABC transporter, a signal peptidase I, and a putative EF-G elongation factor. This region was almost identical in both genomes, indicating a recent transfer event. Additionally, another region containing a relevant, pathogenicity-related CDS (coding for a putative sialidase) might also have been horizontally transferred (Table (Table5).5). The presence of a putative sialidase in both genomes is noteworthy, since this enzyme had been identified only in Mycoplasma alligatoris elsewhere among the mycoplasmas (6). Sialidase cleaves terminal sialic acid residues from sialoglycoconjugates, generating free sialic acid (a likely nutrient). The M. synoviae genome region (216417 to 224099) contains five CDSs encoding enzymes involved in sialic acid catabolism (sialic acid lyase, N-acetylmannosamide kinase, N-acetylmannosamine-6-phosphate epimerase, glucosamine-6-phosphate isomerase, and N-acetylglucosamine-6-phosphate deacetylase). Some of these enzymes are found in selected mycoplasmas (e.g., M. hyopneumoniae, M. mycoides, and M. penetrans), but none of them contains the CDS cluster found in M. synoviae. This suggests that sialic acid could be a substrate for M. synoviae growth, a hypothesis that can be addressed experimentally.
Hemagglutinins play a fundamental role in the pathogenesis of M. synoviae and M. gallisepticum, and their genes could have been transferred between these species (5, 40). However, the organization of hemagglutinin genes differs sharply in M. synoviae and M. gallisepticum, the former with a single locus (1) comprising 70 CDSs (see Fig. Fig.11 in http://www.brgene.lncc.br/finalMS/fig1) and the latter containing 43 genes organized in five loci (44). Whole-genome alignments of M. synoviae and M. gallisepticum allowed the identification of three M. gallisepticum hemagglutinin genes showing a strong similarity to the M. synoviae hemagglutinin CDS cluster. Using the Tribe-MCL algorithm, based on similarity data of the whole set of hemagglutinin CDSs of M. gallisepticum and M. synoviae, three hemagglutinin CDS groups could be identified (see Table Table11 at http://www.brgene.lncc.br/finalMS/table1). The first one contained 41 M. synoviae CDSs, the second contained 35 M. gallisepticum CDSs, and the third contained 29 M. synoviae CDSs and the 3 M. gallisepticum highly homologous CDSs found in whole-genome alignments. These results, and additional phylogenetic and codon usage analyses, support the postulation of a likely transfer event of hemagglutinin genes from M. gallisepticum to M. synoviae.
A deeper understanding of the genomic diversity and evolution of Mycoplasmatales is now possible by analyzing the available sequenced genomes of species belonging to the Hominis clade (M. pulmonis, M. hyopneumoniae, M. synoviae, and M. mobile) and the Pneumoniae clade (M. genitalium, M. pneumoniae, M. gallisepticum, U. urealyticum, and M. penetrans). Apparently, the evolutionary reduction of mycoplasma genomes resulted in preservation of a minimum set of essential metabolic capabilities, rather than a minimum set of specific genes or pathways. The one-carbon pool tetrahydrofolate (C1-THF) metabolism clearly illustrates the alternative retention of redundant metabolic pathways (Fig. (Fig.4).4). The C1-THF metabolism, conserved in both Hominis and Pneumoniae clades, results in formylation of the methionyl-tRNA initiator (Met-tRNAi), as well as de novo dTMP synthesis. M. hyopneumoniae is the only exception, because it is apparently unable to formylate Met-tRNAi due to the absence of both methionyl-tRNA formyltransferase (Fmt) and peptide deformylase (Def).
M. hyopneumoniae is also devoid of enzymes involved in the C1-THF metabolism, except for glycine hydroxymethyltransferase (GlyA). Since all mycoplasmas, except M. hyopneumoniae, possess methylenetetrahydrofolate dehydrogenase (NADP+) (FolD), the activities of GlyA and formate-dihydrofolate ligase (FthS) may be considered metabolically redundant (Fig. (Fig.4).4). M. synoviae, M. pulmonis, M. pneumoniae, M. genitalium, and M. gallisepticum have retained GlyA but lost FthS, while the opposite has occurred in U. urealyticum. In both cases, a minimum set of metabolic interconversions essential for Met-tRNAi formylation and dTMP synthesis has been conserved, although the pathway itself has not. M. penetrans, which has the largest known mycoplasma genome, is the only species that has retained a redundant C1-THF pathway.
Formylation of initiator Met-tRNA provides selectivity for the initiation factor IF2 and also blocks the binding of the elongation factor EF-Tu to the initiator tRNA (39, 54). We noticed the absence of FolD, Fmt, and Def in a recently sequenced plant-pathogenic phytoplasma (42). As far as we know, M. hyopneumoniae and Phytoplasma sp. are the first naturally occurring eubacteria potentially unable to formylate . The physiological consequences, if any, to bacterial growth behavior remain to be studied. Interestingly, disruption of the fmt gene severely impairs growth of Escherichia coli (23). This could be related to the higher binding affinity of IF-2 to than to (24). However, a fmt-deficient strain of Pseudomonas aeruginosa can carry out formylation-independent initiation of protein synthesis, conceivably because IF-2 has a dual substrate specificity (38, 53).
Genome reduction resulted in a complex pattern of losses and retentions in the purine and pyrimidine metabolic pathways. Except for M. pneumoniae and M. genitalium, no other pair of mycoplasma species shares an identical set of pyrimidine metabolism pathways. Uridine kinase was found to be specific for the Pneumoniae clade, while 5′-nucleotidase is restricted to the Hominis clade. The lack of both enzymes, as previously reported for M. penetrans (50), was also observed in M. synoviae. Although UMP may be derived from either uracil or carbamoyl phosphate in M. penetrans (50), M. synoviae relies exclusively on uracil for UMP synthesis. The lack of both enzymes also implies that M. penetrans and M. synoviae cannot metabolize cytidine or uridine, only their deoxy derivatives. Thymidylate synthase (ThyA) and dihydrofolate reductase (FolA) are functionally coupled, since ThyA-mediated, de novo dTMP synthesis requires a continuous reduction of dihydrofolate reductase to THF (9). Interestingly, a novel class of flavin-dependent thymidylate synthases (ThyX) that preserve THF in its reduced form is widely distributed in bacterial genomes (37). Mycoplasmas lacking ThyA (M. hyopneumoniae, U. urealyticum, and M. mycoides) might have an alternative mechanism for de novo dTMP synthesis involving an unknown ThyX activity; alternatively, the dTMP pool could be exclusively dependent on thymidylate kinase, which is present in all mycoplasmas sequenced so far.
Different mechanisms of cell adhesion have evolved among mycoplasmas. Recent studies (29, 44) have shown that, in one branch of the Pneumoniae clade (leading to M. pneumoniae-M. genitalium-M. gallisepticum), mechanisms promoting attachment to host cells are mediated by proteins forming part of a tip organelle. The organization of operons encoding the major cytadhesins in these species (P1, MgPa, and GapA) and their cytadherence-related molecules have been evolutionarily conserved. This major adhesion mechanism is also conserved in M. penetrans, despite the facts that it diverged from M. pneumoniae and M. genitalium some 200 MYBP and that it switched from a mammal host to a bird host. Conversely, in all genomes of the more distantly related M. penetrans-U. urealyticum branch of the Pneumoniae clade and in the Hominis clade, MgPa-like protein CDSs were identified, including two in M. hyopneumoniae and four in M. synoviae. However, CDSs encoding most of the other components of the tip organelle were not found. These species also showed other exclusive or partially shared cell adhesion-related CDSs (Table (Table3),3), suggesting that cell adhesion mechanisms have followed different evolutionary pathways. The diversity of adhesion determinants is indicative of the plasticity of these small genomes (our data) (48), occasionally shuffled by HGT or intraspecific recombination, to generate smaller and adaptive arrangements to specific hosts or microhabitats.
Comparative analyses of mycoplasma genomes allowed the identification of adaptive mechanisms accounting for perpetuation across a wide host range. Despite their small genomes, mycoplasmas showed a heterogeneous gene composition, with several species-specific genes whose identification and functional characteristics might be helpful for the prevention and treatment of diseases caused by these bacteria. Our analyses confirmed the occurrence of high rates of genomic rearrangements in mycoplasmas, which was demonstrated by the presence of strain-specific regions in the three M. hyopneumoniae strains; some of these regions were probably involved in rearrangements or pathogenesis. In fact, a variety of CDSs coding for putative, outer membrane proteins, or adhesins, were found to contain motifs of repetitive sequences which might have a role in their biological function or in antigenic variation. The presence of an ICEH element restricted to pathogenic strains suggests its possible role in pathogenicity. For the first time, phylogenetic relationships of all sequenced mycoplasma genomes were established on the basis of a concatenated data set, which resulted in a single, highly supported tree. These data were used for estimating a time frame of genome evolution, which added new insights to the evolution of the Mycoplasmatales-Entomoplasmatales group. Our studies provide evidence pointing to HGT as the process that provided M. synoviae and M. gallisepticum with the capacity of infecting the same host.
The present and former staffs of the Ministério da Ciência e Tecnologia (MCT)/Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) are gratefully acknowledged for their strategic vision and enthusiastic support. We are also indebted to Juçara Parra (Ludwig Institute for Cancer Research) for administrative coordination. We thank the following individuals for technical and logistical expert assistance: Antonio Carneiro Kindermann (Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul); Marni Ramenzoni (CNPSA, EMBRAPA); José Fernando L. Machado, Jr., and João Francisco Valiati (Instituto de Informática, Universidade Federal do Rio Grande do Sul); Artur Luiz da Costa da Silva, Maria Silvanira Ribeiro Barbosa, and Juliana Simão Nina de Azevedo (Universidade Federal do Para, Belém, PA, Brazil); Jacqueline da Silva Batista, Jorge Ivan Rebelo Porto, José Antonio Alves Gomes, Alexandra Regina Bentes de Sousa, Naiara Alessandra Bertucchi Vogt, Tatiana Leite Marão, Audrey Alencar Arruda d′Assunção, and Kyara de Aquino Formiga (Instituto Nacional de Pesquisas da Amazônia); Mário Steindel (Universidade Federal de Santa Catarina); Janice Silva Sales (Universidade Federal de Alagoas); Emanuel Maltempi de Souza (Universidade Federal do Paraná); Raquel Liboredo Santos (Universidade Federal de Minas Gerais); Carlos Alfredo Galindo Blaha and Marbela Maria Fonseca (Universidade Federal do Rio Grande do Norte); Carla Cristina Jaremtchuk, Luiz Renato Kotlevski, Jr., and Rafael Andrzejewski (Pontifícia Universidade Católica do Paraná); Ebert Seixas Hanna (Universidade de São Paulo); Michelle Maia Jardim and Raquel Liboredo Santos (Universidade Federal de Minas Gerais); and Julia Araripe and Edson Rondinelli (Universidade Federal do Rio de Janeiro).
This work was undertaken by the Brazilian National Genome Program (Southern Network for Genome Analysis and Brazilian National Genome Project Consortium) with funding provided by MCT/CNPq and SCT/FAPERGS (RS).
†Supplemental material for this article may be found at http://jb.asm.org/.