The Mycoplasma conjunctivae genome sequencing, annotation and analysis

<p>Abstract</p> <p>Background</p> <p>The mollicute <it>Mycoplasma conjunctivae </it>is the etiological agent leading to infectious keratoconjunctivitis (IKC) in domestic sheep and wild caprinae. Although this pathogen is relatively benign for domestic animals treated by antibiotics, it can lead wild animals to blindness and death. This is a major cause of death in the protected species in the Alps (e.g., <it>Capra ibex</it>, <it>Rupicapra rupicapra</it>).</p> <p>Methods</p> <p>The genome was sequenced using a combined technique of GS-FLX (454) and Sanger sequencing, and annotated by an automatic pipeline that we designed using several tools interconnected via PERL scripts. The resulting annotations are stored in a MySQL database.</p> <p>Results</p> <p>The annotated sequence is deposited in the EMBL database (<ext-link ext-link-type="embl" ext-link-id="FM864216">FM864216</ext-link>) and uploaded into the mollicutes database MolliGen <url>http://cbi.labri.fr/outils/molligen/</url> allowing for comparative genomics.</p> <p>Conclusion</p> <p>We show that our automatic pipeline allows for annotating a complete mycoplasma genome and present several examples of analysis in search for biological targets (e.g., pathogenic proteins).</p

Mycoplasma agalactiae ST35: a new sequence type with a minimal accessory genome primarily affecting goats

Background: Mycoplasma agalactiae, causing agent of contagious agalactia, infects domestic small ruminants such as sheep and goats but also wild Caprinae. M. agalactiae is highly contagious and transmitted through oral, respiratory, and mammary routes spreading rapidly in an infected herd. Results: In an outbreak of contagious agalactia in a mixed herd of sheep and goats, 80% of the goats were affected displaying swollen udders and loss of milk production but no other symptom such as kerato-conjunctivitis, arthritis or pulmonary distress commonly associated to contagious agalactia. Surprisingly, none of the sheep grazing on a common pasture and belonging to the same farm as the goats were affected. Whole genome sequencing and analysis of M. agalactiae strain GrTh01 isolated from the outbreak, revealed a previously unknown sequence type, ST35, and a particularly small, genome size of 841′635 bp when compared to others available in public databases. Overall, GrTh01 displayed a reduced accessory genome, with repertoires of gene families encoding variable surface proteins involved in host-adhesion and variable antigenicity being scaled down. GrTh01 was also deprived of Integrative Conjugative Element or prophage, and had a single IS element, suggesting that GrTh01 has a limited capacity to adapt and evolve. Conclusions: The lack of most of the variable antigens and the Integrative Conjugative Element, both major virulence- and host specificity factors of a M. agalactiae strain isolated from an outbreak affecting particularly goats, indicates the implication of these factors in host specificity. Whole genome sequencing and full assembly of bacterial pathogens provides a most valuable tool for epidemiological and virulence studies of M. agalactiae without experimental infections. Keywords: Mycoplasma agalactiae, Sequence type 35, Goats, Full genome, Contagious agalacti

Comparative genomic analyses of Mycoplasma hyopneumoniae pathogenic 168 strain and its high-passaged attenuated strain

Background: Mycoplasma hyopneumoniae is the causative agent of porcine enzootic pneumonia (EP), a mild, chronic pneumonia of swine. Despite presenting with low direct mortality, EP is responsible for major economic losses in the pig industry. To identify the virulence-associated determinants of M. hyopneumoniae, we determined the whole genome sequence of M. hyopneumoniae strain 168 and its attenuated high-passage strain 168-L and carried out comparative genomic analyses. Results: We performed the first comprehensive analysis of M. hyopneumoniae strain 168 and its attenuated strain and made a preliminary survey of coding sequences (CDSs) that may be related to virulence. The 168-L genome has a highly similar gene content and order to that of 168, but is 4,483 bp smaller because there are 60 insertions and 43 deletions in 168-L. Besides these indels, 227 single nucleotide variations (SNVs) were identified. We further investigated the variants that affected CDSs, and compared them to reported virulence determinants. Notably, almost all of the reported virulence determinants are included in these variants affected CDSs. In addition to variations previously described in mycoplasma adhesins (P97, P102, P146, P159, P216, and LppT), cell envelope proteins (P95), cell surface antigens (P36), secreted proteins and chaperone protein (DnaK), mutations in genes related to metabolism and growth may also contribute to the attenuated virulence in 168-L. Furthermore, many mutations were located in the previously described repeat motif, which may be of primary importance for virulence. Conclusions: We studied the virulence attenuation mechanism of M. hyopneumoniae by comparative genomic analysis of virulent strain 168 and its attenuated high-passage strain 168-L. Our findings provide a preliminary survey of CDSs that may be related to virulence. While these include reported virulence-related genes, other novel virulence determinants were also detected. This new information will form the foundation of future investigations into the pathogenesis of M. hyopneumoniae and facilitate the design of new vaccines

The mysterious orphans of Mycoplasmataceae

Background: The length of a protein sequence is largely determined by its function, i.e. each functional group is associated with an optimal size. However, comparative genomics revealed that proteins length may be affected by additional factors. In 2002 it was shown that in bacterium Escherichia coli and the archaeon Archaeoglobus fulgidus, protein sequences with no homologs are, on average, shorter than those with homologs. Most experts now agree that the length distributions are distinctly different between protein sequences with and without homologs in bacterial and archaeal genomes. In this study, we examine this postulate by a comprehensive analysis of all annotated prokaryotic genomes and focusing on certain exceptions. Results: We compared lengths distributions of having homologs proteins (HHPs) and non-having homologs proteins (orphans or ORFans) in all currently annotated completely sequenced prokaryotic genomes. As expected, the HHPs and ORFans have strikingly different length distributions in almost all genomes. As previously established, the HHPs, indeed, are, on average, longer than the ORFans, and the length distributions for the ORFans have a relatively narrow peak, in contrast to the HHPs, whose lengths spread over a wider range of values. However, about thirty genomes do not obey these rules. Practically all genomes of Mycoplasma and Ureaplasma have atypical ORFans distributions, with the mean lengths of ORFan larger than the mean lengths of HHPs. These genera constitute over 80% of atypical genomes. Conclusions: We confirmed on a ubiquitous set of genomes the previous observation that HHPs and ORFans have different gene length distributions. We also showed that Mycoplasmataceae genomes have distinctive distributions of ORFans lengths. We offer several possible biological explanations of this phenomenon

Comparative Genomics of Mycoplasma: Analysis of Conserved Essential Genes and Diversity of the Pan-Genome

Mycoplasma, the smallest self-replicating organism with a minimal metabolism and little genomic redundancy, is expected to be a close approximation to the minimal set of genes needed to sustain bacterial life. This study employs comparative evolutionary analysis of twenty Mycoplasma genomes to gain an improved understanding of essential genes. By analyzing the core genome of mycoplasmas, we finally revealed the conserved essential genes set for mycoplasma survival. Further analysis showed that the core genome set has many characteristics in common with experimentally identified essential genes. Several key genes, which are related to DNA replication and repair and can be disrupted in transposon mutagenesis studies, may be critical for bacteria survival especially over long period natural selection. Phylogenomic reconstructions based on 3,355 homologous groups allowed robust estimation of phylogenetic relatedness among mycoplasma strains. To obtain deeper insight into the relative roles of molecular evolution in pathogen adaptation to their hosts, we also analyzed the positive selection pressures on particular sites and lineages. There appears to be an approximate correlation between the divergence of species and the level of positive selection detected in corresponding lineages

Complexity of the Mycoplasma fermentans M64 Genome and Metabolic Essentiality and Diversity among Mycoplasmas

Recently, the genomes of two Mycoplasma fermentans strains, namely M64 and JER, have been completely sequenced. Gross comparison indicated that the genome of M64 is significantly bigger than the other strain and the difference is mainly contributed by the repetitive sequences including seven families of simple and complex transposable elements ranging from 973 to 23,778 bps. Analysis of these repeats resulted in the identification of a new distinct family of Integrative Conjugal Elements of M. fermentans, designated as ICEF-III. Using the concept of “reaction connectivity”, the metabolic capabilities in M. fermentans manifested by the complete and partial connected biomodules were revealed. A comparison of the reported M. pulmonis, M. arthritidis, M. genitalium, B. subtilis, and E. coli essential genes and the genes predicted from the M64 genome indicated that more than 73% of the Mycoplasmas essential genes are preserved in M. fermentans. Further examination of the highly and partly connected reactions by a novel combinatorial phylogenetic tree, metabolic network, and essential gene analysis indicated that some of the pathways (e.g. purine and pyrimidine metabolisms) with partial connected reactions may be important for the conversions of intermediate metabolites. Taken together, in light of systems and network analyses, the diversity among the Mycoplasma species was manifested on the variations of their limited metabolic abilities during evolution

Microbial genomic taxonomy

A need for a genomic species definition is emerging from several independent studies worldwide. In this commentary paper, we discuss recent studies on the genomic taxonomy of diverse microbial groups and a unified species definition based on genomics. Accordingly, strains from the same microbial species share >95% Average Amino Acid Identity (AAI) and Average Nucleotide Identity (ANI), >95% identity based on multiple alignment genes, 70% in silico Genome-to-Genome Hybridization similarity (GGDH). Species of the same genus will form monophyletic groups on the basis of 16S rRNA gene sequences, Multilocus Sequence Analysis (MLSA) and supertree analysis. In addition to the established requirements for species descriptions, we propose that new taxa descriptions should also include at least a draft genome sequence of the type strain in order to obtain a clear outlook on the genomic landscape of the novel microbe. The application of the new genomic species definition put forward here will allow researchers to use genome sequences to define simultaneously coherent phenotypic and genomic groups