BACKGROUND
mycoplasmas  are among the smallest microorganisms capable of self-replication and autonomous life  <cit> .
the genus mycoplasma includes a large number of highly genomically-reduced species which in nature are associated with hosts either commensally or pathogenically  <cit> .
general features of the class mollicutes are small genome, lack of cell wall and low gc content.
indeed, the mycoplasma species have genomes of  <dig>  to  <dig>  mbp.
weisburg et al.
<cit>  and woese et al.
<cit>  revealed that mycoplasma have evolved from more classical bacteria of the firmicutes taxon by a so-called regressive evolution that resulted in massive genome reduction  <cit>  and minimal metabolic activities.
consequently, they adopted a strict parasitic life style, mainly occurring as extracellular parasites often restricted to a living host, with some species having the ability to invade host cells as described by sirand-pugnet et al.
<cit> , rosengarten et al.
<cit>  and citti et al.
<cit> .
they have a predilection for the mucosal surfaces, where they successfully compete for nutrients with many other organisms, establishing chronic infections  <cit> .
they do not show specific virulence factor as known in other bacteria, instead they seem to use toxic metabolic intermediates that they secrete and translocate to the host cells as virulence factors  <cit> .
additionally, due to the lack of cell wall, they are not affected by some antibiotics which target synthesis of cell wall such penicillin or other beta-lactam antibiotics making these organisms particularly interesting in medicine.
infectious keratoconjunctivitis
mycoplasma conjunctivae is considered as the major etiological agent of infectious keratoconjunctivitis  for both domestic and wild caprinae species.
in the european alps it affects several species such as alpine ibex , alpine chamois , and mouflon , as well as in domestic sheep and goat  <cit> .
in switzerland, m. conjunctivae is known to be the primary cause of this disease  <cit> .
the implied role of m. conjunctivae is based on the frequent isolation of this organism from inflamed eyes and on limited attempts to induce ocular disease experimentally showing that m. conjunctivae is one agent responsible for epidemic keratoconjunctivitis  <cit> .
nonetheless, even if the molecular epidemiology has been well described by belloy et al.
<cit> , the molecular infection mechanism is still not established and remains a mystery.
methods
bacterial strain
m. conjunctivae type strain hrc/581t   <cit>  was grown on standard mycoplasma broth medium enriched with 20% horse serum,  <dig> % yeast extract and 1% glucose .
the cells were harvested by centrifugation at  <dig>  × g for  <dig> min, washed three times in tes buffer , and then re-suspended in tes buffer to a concentration of approximately  <dig> bacteria/ml.
dna was extracted by the guanidium thiocyanate method  <cit> , extracted  <dig> times with pcia  and  <dig> times with cia , precipitated with 50% isoproanol, washed  <dig> times with 70% ethanol to remove salt, dried in the air for  <dig> min and re-suspended in double distilled h2o at a concentration of  <dig> μg/ml.
sequencing
sequencing and assembly of the genome was carried out by microsynth ag.
the quality of the isolated genomic dna was verified by gel electrophoresis and displayed a pure high molecular weight dna.
the dna was sheared by passing it several times through a needle, in order to construct two different libraries: a plasmid library and a fosmid library.
for the plasmid library , the genomic dna was passed  <dig> times through a 30-gauge needle and sonicated for  <dig> seconds .
for the fosmid library , the genomic dna was passed  <dig> times through a 23-gauge needle without sonication.
small fragments were ligated with a linker, fractionated twice through  <dig> % agarose gels.
fractions of  <dig> different sizes  were cut out from the gel and cloned into vector potw <dig> .
moreover, the large fragments were fractionated using a chef-dr ii system .
fragments of  <dig> kbp were cut out from the gel and ligated into pcc1fos .
from the plasmid library 11' <dig> clones and from the fosmid library  <dig> clones were end-sequenced on an abi  <dig> capillary sequencer.
a second part of the small fragments were sequenced using  <dig> life science flx technology leading to 263' <dig> reads that were reduced to 78' <dig> reads covering 20'569' <dig> bp after applying a quality cut-off filter .
assembly
the assembly was carried out using the seqman module of the dnastar lasergene version  <dig> combining both classical sanger sequences  and  <dig> flx reads.
a check was conducted with "amosvalidator" of the amos package  <cit> , allowing identifying suspicious regions in the assembly.
to help in the assembly process the  <dig> fosmids paired-end reads were aligned to the final sequence.
the reads display a nice spreading at regular intervals except for  <dig> clones that were absent from the results.
the  <dig> regions were analyzed for the presence of potentially lethal genes for e. coli.
the first region contains a homologue of the gene lepa that is known to be lethal when overexpressed in e. coli  <cit> .
this region also contains a restriction enzyme that might cut e. coli genome.
the second region contains a transposase and some phage genes.
this might explain the toxicity of these two fosmids in e. coli.
automatic annotation
the automatic annotation pipeline was entirely built locally using available software and linking them with perl scripts.
gene prediction was carried out using glimmer  <dig>   <cit>  and the genetic code specific for mycoplasma .
the interpolated context models  were calculated by self-training on the long orfs of the contigs.
the rnas were predicted using infernal with models obtained from rfam  <cit> , trnascan-se  <cit> , and blastn for 16s and 23s  <cit> .
predicted coding sequences  were translated using the emboss package   <cit>  and a similarity search was run by blastp against the uniprot/swiss-prot knowledgebase   <cit> .
the cds were also scanned against the hamap families  <cit>  to identify orthologous protein families.
in addition the cds were searched for potential known domains using interproscan  <cit> , and for biased compositional regions with seg  <cit>  and marcoil  <cit> .
the biological interest of an annotation project is to identify the gene products by designating a descriptive common name for the protein and its function with as much specificity as the evidence supports.
we use homology-based annotation transfer to assign the name and associated information of gene product: gene symbol, ec number if protein is identified as an enzyme and other features.
homology search is performed by blastp that allows finding the best matches with the highest significant sequence similarity appearing between the putative proteins sequences compared first to a database of known mycoplasma proteins and secondly to proteins from the uniprot/swiss-prot knowledgebase.
additionally, matches with hamap protein family permits to support homology annotation and raise the confidence level of annotation transfer.
other characterization features like functional domains constitute an additional support evidence of function assignment .
the results obtained from the various programs are parsed and stored in gff <dig> format in a local mysql database.
the embl format is produced from the data stored in this database and deposited at the ebi embl database, the accession number is fm <dig>
description of the criteria used to assign the genes products into the  <dig> following categories: known protein , putative protein  and hypothetical protein .
in order to assess the confidence of the results provided by our annotation pipeline, we used the known genome of mycoplasma hyopneumoniae  already annotated  for comparing pipeline results with those provided at ncbi.
the annotation pipeline identifies  <dig> cds and  <dig> non-coding rnas  whereas only  <dig> cds and  <dig> ncrnas are known.
the  <dig> % of total genes were correctly predicted leaving only  <dig> genes not found by the predictor programs.
regarding the functional annotation, our pipeline provides  <dig> % of correct gene annotations and even in some cases complements the existing one.
<dig> % of cds are wrongly annotated or in a different way.
the  <dig> genes predicted in addition to known genes could be considered as false positives even if they also could represent new potential genes of m. hyopneumoniae genome .
using results of m. hyopneumoniae  <dig> annotation, we evaluated the specificity and sensitivity of the pipeline and obtained a sensitivity of 92% and a specificity of 90% .
* cds were predicted using glimmer <dig> and ncrnas were predicted using blastn for 23s and 16s rrnas, trnascan-se for trnas and infernal for other ncrnas
fp = false positive, tp = true positive, fn = false negative
RESULTS
general genome features
composition and functional gene assignment, origin of replication
the genome consists of a single chromosome with a size estimated of about  <dig>  mbp.
currently, more than 95% of the genome is available as a single contig, but due to the presence of repeated sequences, we experienced difficulties in assembling the sequences to close the gap.
the contig has a size of 846' <dig> bp with a g+c content of 29% .
a total of  <dig> genes have been computationally predicted.
we found both 23s and 16s ribosomal rnas in unique copies located next to each other.
the 5s ribosomal rna is located remotely of 23s and 16s genes.
we identified  <dig> transfer rnas covering all  <dig> amino acids.
other non-coding rnas were found: bacterial rnase p class b, tpp riboswitch , tmrna  and the bacterial signal recognition particle rna.
additionally,  <dig> genes were predicted as coding sequences for proteins.
49% of those genes have a clear homologue with a known function.
<dig> % of those genes were annotated as "putative proteins" because the closest known protein was aligned with a marginal e-value.
while the remaining 45% have unknown function, and were named as "hypothetical" proteins.
from those hypothetical proteins, 75% matched non-significatively to other proteins and 25% are considered unique for m. conjunctivae because no match was obtained by blastp against both databases.
a summary is shown in table  <dig> and table  <dig>
it is important to note that our method  does not allow to distinguish between close homologues having different functions.
the origin of replication  was searched by comparing the sequences of  <dig> strains of m. hyopneumoniae with the sequence of m. conjunctivae in the region of dnaa gene .
we attempted to identify several features that have been associated with replication origins in other bacterial species, including mollicutes.
bacterial origins of replication are typically located in the vicinity of the dnaa and dnan genes.
usually several dnaa-box motifs are found within the intergenic regions around dnaa gene  <cit> .
we searched unsuccessfully for the presence of consensus dnaa-box motifs with the pattern ttatc a  using fuzznuc of the emboss package  <cit> .
when we used a slightly different, more relaxed dnaa-box consensus motifs tt   a a, two sequences matching each of these patterns were found between the dnaa and rpmh genes .
however, well over  <dig>  hits located throughout the rest of the genome were also seen.
therefore, the specificity of the pattern used to try to detect dnaa-box motifs was very low, decreasing our confidence in the significance of the sequences identified.
these findings are in contrast to the multiple dnaa-boxes found in the intergenic regions surrounding dnaa in other mollicutes  <cit> .
in addition to the presence of dnaa-box motifs, replication origins can also frequently be identified by looking for biases in strand composition through measures such as the cumulative gc skew  <cit> .
for m. conjunctivae, we found no significant asymmetries that can be readily detected with gc skew.
the lack of a clear bias in m. conjunctivae is similar to that observed for the m. hyopneumoniae  <cit> .
therefore, the only significant feature of the m. conjunctivae genome that provides any possible indication of the location of the origin of replication is the presence of the dnaa gene.
otherwise, there are no features that allow definitive mapping of the origin to the intergenic region upstream of the dnaa gene, as seen in other bacteria.
potential pathogenic features
bacteria have many ways to produce virulence that reside in the ability to adhere, invade and cause damage to host cells.
various strategies of pathogenicity such as cytolysins, toxins and invasins enable other bacteria to produce infection.
in mycoplasma species no such typical primary virulence genes have been found.
mycoplasmas seem rather to use intrinsic metabolic and catabolic functions to cause disease in the affected host and to ensure the microbe's survival.
our efforts to identify genes involved in the pathogenicity of mycoplasma conjunctivae were concentrated on the one hand, try to find those primary virulence genes, toxins principally, rare in other mycoplasmas.
on the other hand, on metabolic pathways that has been proposed by studies carried out in other mycoplasmas  <cit> .
glycerol pathway
we found using manual blastp queries by an expert, the genes for a glycerol-3-phosphate dehydrogenase , a glycerol kinase , a glycerol uptake facilitator protein  and an abc transporter system  that are implicated in the glycerol metabolism producing cell damage, inflammation and disease in mycoplasma mycoides subsp.
mycoides small colony   <cit> .
the pathway starts with the assimilation of glycerol by the abc glycerol transporter .
afterwards, the glycerol is phosphorylated into glycerol-3-phosphate, then oxidized by glpo in presence of o <dig> into dihydroxyactone-phosphate  producing one molecule of h2o <dig>  h2o <dig> is released directly inside the host cells by the transmembrane glpo protein leading to cell death  <cit> .
the absence of any gene having a catalase or dismutase activity favors this hypothesis.
the identification of those genes in mycoplasma conjunctivae constitutes an important discovery given that a relationship between the glycerol metabolism and cytotoxicity is established in the laboratory <cit> .
further work to validate this hypothesis in m. conjunctivae is required and has been started in collaboration with a laboratory of the institute for veterinary bacteriology .
toxins
toxins constitute an important type of virulence factors in several bacteria.
thereby, we searched for toxins in m. conjunctivae and we found  <dig> proteins highly similar with toxins of treponema hyodysenteriae .
those proteins are hemolysin a , hemolysin b  and hemolysin c .
the  <dig> genes are scattered on the genome.
those proteins are present in other mycoplasmas, particularly m. hyopneumoniae and m. capricolum, and even if in those species, these toxins are not essential for pathogenicity mechanisms, it can not be excluded that these toxins contribute to the pathogenicity of m. conjunctivae.
is elements
insertion sequences  are short dna elements that function as simple transposable elements by coding for proteins implicated in the transposition activity.
transposase and other regulatory protein are the proteins generally coded by is elements: the transposase catalyses the enzymatic reaction allowing the is to move.
regulatory proteins act by enhancing or inhibiting the transposition activity.
the coding region in an insertion sequence is usually flanked by inverted repeats  <cit> .
we found several genes coding for complete or partial transposases .
an is <dig> insertion element has particularly brought our attention.
is <dig> elements belong to is <dig> family are prevalent in other mycoplasmas and are the only  that have been demonstrated directly to undergo autonomous transposition  <cit> .
interestingly a transposase for one is <dig> insertion elements is followed by homologues of a methylase hpai and a type ii restriction enzyme hpai from haemophilus parainfluenzae forming a restriction-methylation cassette.
the hypothesis of a horizontal transfer from h. parainfluenza to m. conjunctivae was formulated.
we evaluated the g+c content of this cassette, but we did not observe a higher g+c content inside the cassette compared to the surrounding area.
if the g+c content inside the cassette would be different from that of m. conjunctivae  and similar to that of h. parainfluenzae  it could constitute an evidence of the transfer.
comparative genomics
the list of proteins was classified and compared to  <dig> other mycoplasma genomes as shown in table  <dig>  the main difference with other mycoplasma is an apparently low carbohydrate and transport metabolism that could explain the need for a strong glycerol pathway, as well as the large number of hypothetical proteins probably due to the fully automatic annotation process.
mycoplasma conjunctivae
discussion
mycoplasma conjunctivae is the fourteenth genome of a mycoplasma species that has been fully sequenced.
phylogenetically, the closest relative among the sequenced mycoplasmas is m. hyopneumoniae reflected by the high similarity of most of the proteins identified in m. conjunctivae.
the analysis of m. conjunctivae genome features, describes this organism as a typical mycoplasma, with a genome size and a g+c content within the range of other mycoplasma genomes.
the comparison of mycoplasma genome sizes demonstrates that the sequence length is variable not only within the same genus but even among strains of the same species as shown in table  <dig>  even if we do not know the final size of the genome, we expect a chromosome length of about 900' <dig> bp, size almost similar with the genome size of m. hyopneumoniae.
mycoplasma conjunctivae
globally the mycoplasma genomes have a characteristically low g+c content within the range of  <dig>  to  <dig> mol% .
the highest g+c content found in m. pneumoniae and the lowest in m. capricolum.
regarding mycoplasma conjunctivae, g+c content has a typical value of about 29%.
the codon usage is similar to that of m. hyopneumoniae and opposite to that of m. capricolum and m. mycoides  <cit> .
the presence of repeats across the genome was the principal difficulty for finishing the genome assembly.
insertion sequence  elements are reported in the majority of mycoplasmas and in m. conjunctivae we found transposases for is-elements in the genome.
some of those transposases genes are complete sequences and some other are fragmented showing a predicted length of less than  <dig> bp.
since those insertion elements are nearly identical they created difficulties for assembling the genome.
the findings highlighted by this project, principally the glycerol pathway, require further experimental confirmation.
in particular, the hypothesis for damaging the host cells by the glycerol metabolism need to be confirmed by demonstrating the localization of glpo in the membrane and the release of h2o <dig> outside the cell.
if this hypothesis can be verified, the possibility to block at any stage the glycerol pathway could constitute a candidate target for controlling the disease.
CONCLUSIONS
in conclusion, we created an automatic pipeline to annotate a prokaryotic genome sequence using various tools for the prediction and the identification of the genes.
this pipeline is customized for handling sequences of mycoplasma species.
we deposited the mycoplasma conjunctivae genome fully annotated in the embl database .
data stored into our local database can be searched and genome can be visualized through our website .
analysis of annotated genes gives new insights about potential mechanisms of pathogenicity as well as the possibility to go deeper into the knowledge of mycoplasma conjunctivae and the ikc disease and opens the way to finding methods to prevent m. conjunctivae infections of domestic animals as reservoir for this pathogen and hence prevent ikc in wild animals.
competing interests
the authors declare that they have no competing interests.
authors' contributions
jf proposed and supported the project.
spcc created the pipeline and did the analysis.
maq did the genomic library.
ts and maq sequenced the clones.
ts, gw and cw did the assembly.
spcc and lf wrote the manuscript.

 mycoplasmas are among the smallest microorganisms capable of self-replication and autonomous life. the genus mycoplasma includes a large number of highly genomically-reduced species which in nature are associated with hosts either commensally or pathogenically. general features of the class mollicutes are small genome, lack of cell wall and low gc content. indeed, the mycoplasma species have genomes of to mbp. weisburg et al. and woese et al. revealed that mycoplasma have evolved from more classical bacteria of the firmicutes taxon by a so-called regressive evolution that resulted in massive genome reduction and minimal metabolic activities. consequently, they adopted a strict parasitic life style, mainly occurring as extracellular parasites often restricted to a living host, with some species having the ability to invade host cells as described by sirand-pugnet et al., rosengarten et al. and citti et al.. they have a predilection for the mucosal surfaces, where they successfully compete for nutrients with many other organisms, establishing chronic infections. they do not show specific virulence factor as known in other bacteria, instead they seem to use toxic metabolic intermediates that they secrete and translocate to the host cells as virulence factors. additionally, due to the lack of cell wall, they are not affected by some antibiotics which target synthesis of cell wall such penicillin or other beta-lactam antibiotics making these organisms particularly interesting in medicine. infectious keratoconjunctivitis mycoplasma conjunctivae is considered as the major etiological agent of infectious keratoconjunctivitis for both domestic and wild caprinae species. in the european alps it affects several species such as alpine ibex, alpine chamois, and mouflon, as well as in domestic sheep and goat. in switzerland, m. conjunctivae is known to be the primary cause of this disease. the implied role of m. conjunctivae is based on the frequent isolation of this organism from inflamed eyes and on limited attempts to induce ocular disease experimentally showing that m. conjunctivae is one agent responsible for epidemic keratoconjunctivitis. nonetheless, even if the molecular epidemiology has been well described by belloy et al., the molecular infection mechanism is still not established and remains a mystery. methods bacterial strain m. conjunctivae type strain hrc/581t was grown on standard mycoplasma broth medium enriched with 20% horse serum, % yeast extract and 1% glucose. the cells were harvested by centrifugation at × g for min, washed three times in tes buffer, and then re-suspended in tes buffer to a concentration of approximately bacteria/ml. dna was extracted by the guanidium thiocyanate method, extracted times with pcia and times with cia, precipitated with 50% isoproanol, washed times with 70% ethanol to remove salt, dried in the air for min and re-suspended in double distilled h2o at a concentration of μg/ml. sequencing sequencing and assembly of the genome was carried out by microsynth ag. the quality of the isolated genomic dna was verified by gel electrophoresis and displayed a pure high molecular weight dna. the dna was sheared by passing it several times through a needle, in order to construct two different libraries: a plasmid library and a fosmid library. for the plasmid library, the genomic dna was passed times through a 30-gauge needle and sonicated for seconds. for the fosmid library, the genomic dna was passed times through a 23-gauge needle without sonication. small fragments were ligated with a linker, fractionated twice through % agarose gels. fractions of different sizes were cut out from the gel and cloned into vector potw. moreover, the large fragments were fractionated using a chef-dr ii system. fragments of kbp were cut out from the gel and ligated into pcc1fos. from the plasmid library 11' clones and from the fosmid library clones were end-sequenced on an abi capillary sequencer. a second part of the small fragments were sequenced using life science flx technology leading to 263' reads that were reduced to 78' reads covering 20'569' bp after applying a quality cut-off filter. assembly the assembly was carried out using the seqman module of the dnastar lasergene version combining both classical sanger sequences and flx reads. a check was conducted with "amosvalidator" of the amos package, allowing identifying suspicious regions in the assembly. to help in the assembly process the fosmids paired-end reads were aligned to the final sequence. the reads display a nice spreading at regular intervals except for clones that were absent from the results. the regions were analyzed for the presence of potentially lethal genes for e. coli. the first region contains a homologue of the gene lepa that is known to be lethal when overexpressed in e. coli. this region also contains a restriction enzyme that might cut e. coli genome. the second region contains a transposase and some phage genes. this might explain the toxicity of these two fosmids in e. coli. automatic annotation the automatic annotation pipeline was entirely built locally using available software and linking them with perl scripts. gene prediction was carried out using glimmer and the genetic code specific for mycoplasma. the interpolated context models were calculated by self-training on the long orfs of the contigs. the rnas were predicted using infernal with models obtained from rfam, trnascan-se, and blastn for 16s and 23s. predicted coding sequences were translated using the emboss package and a similarity search was run by blastp against the uniprot/swiss-prot knowledgebase. the cds were also scanned against the hamap families to identify orthologous protein families. in addition the cds were searched for potential known domains using interproscan, and for biased compositional regions with seg and marcoil. the biological interest of an annotation project is to identify the gene products by designating a descriptive common name for the protein and its function with as much specificity as the evidence supports. we use homology-based annotation transfer to assign the name and associated information of gene product: gene symbol, ec number if protein is identified as an enzyme and other features. homology search is performed by blastp that allows finding the best matches with the highest significant sequence similarity appearing between the putative proteins sequences compared first to a database of known mycoplasma proteins and secondly to proteins from the uniprot/swiss-prot knowledgebase. additionally, matches with hamap protein family permits to support homology annotation and raise the confidence level of annotation transfer. other characterization features like functional domains constitute an additional support evidence of function assignment. the results obtained from the various programs are parsed and stored in gff format in a local mysql database. the embl format is produced from the data stored in this database and deposited at the ebi embl database, the accession number is fm description of the criteria used to assign the genes products into the following categories: known protein, putative protein and hypothetical protein. in order to assess the confidence of the results provided by our annotation pipeline, we used the known genome of mycoplasma hyopneumoniae already annotated for comparing pipeline results with those provided at ncbi. the annotation pipeline identifies cds and non-coding rnas whereas only cds and ncrnas are known. the % of total genes were correctly predicted leaving only genes not found by the predictor programs. regarding the functional annotation, our pipeline provides % of correct gene annotations and even in some cases complements the existing one. % of cds are wrongly annotated or in a different way. the genes predicted in addition to known genes could be considered as false positives even if they also could represent new potential genes of m. hyopneumoniae genome. using results of m. hyopneumoniae annotation, we evaluated the specificity and sensitivity of the pipeline and obtained a sensitivity of 92% and a specificity of 90%. * cds were predicted using glimmer and ncrnas were predicted using blastn for 23s and 16s rrnas, trnascan-se for trnas and infernal for other ncrnas fp = false positive, tp = true positive, fn = false negative  general genome features composition and functional gene assignment, origin of replication the genome consists of a single chromosome with a size estimated of about mbp. currently, more than 95% of the genome is available as a single contig, but due to the presence of repeated sequences, we experienced difficulties in assembling the sequences to close the gap. the contig has a size of 846' bp with a g+c content of 29%. a total of genes have been computationally predicted. we found both 23s and 16s ribosomal rnas in unique copies located next to each other. the 5s ribosomal rna is located remotely of 23s and 16s genes. we identified transfer rnas covering all amino acids. other non-coding rnas were found: bacterial rnase p class b, tpp riboswitch, tmrna and the bacterial signal recognition particle rna. additionally, genes were predicted as coding sequences for proteins. 49% of those genes have a clear homologue with a known function. % of those genes were annotated as "putative proteins" because the closest known protein was aligned with a marginal e-value. while the remaining 45% have unknown function, and were named as "hypothetical" proteins. from those hypothetical proteins, 75% matched non-significatively to other proteins and 25% are considered unique for m. conjunctivae because no match was obtained by blastp against both databases. a summary is shown in table and table it is important to note that our method does not allow to distinguish between close homologues having different functions. the origin of replication was searched by comparing the sequences of strains of m. hyopneumoniae with the sequence of m. conjunctivae in the region of dnaa gene. we attempted to identify several features that have been associated with replication origins in other bacterial species, including mollicutes. bacterial origins of replication are typically located in the vicinity of the dnaa and dnan genes. usually several dnaa-box motifs are found within the intergenic regions around dnaa gene. we searched unsuccessfully for the presence of consensus dnaa-box motifs with the pattern ttatc a using fuzznuc of the emboss package. when we used a slightly different, more relaxed dnaa-box consensus motifs tt a a, two sequences matching each of these patterns were found between the dnaa and rpmh genes. however, well over hits located throughout the rest of the genome were also seen. therefore, the specificity of the pattern used to try to detect dnaa-box motifs was very low, decreasing our confidence in the significance of the sequences identified. these findings are in contrast to the multiple dnaa-boxes found in the intergenic regions surrounding dnaa in other mollicutes. in addition to the presence of dnaa-box motifs, replication origins can also frequently be identified by looking for biases in strand composition through measures such as the cumulative gc skew. for m. conjunctivae, we found no significant asymmetries that can be readily detected with gc skew. the lack of a clear bias in m. conjunctivae is similar to that observed for the m. hyopneumoniae. therefore, the only significant feature of the m. conjunctivae genome that provides any possible indication of the location of the origin of replication is the presence of the dnaa gene. otherwise, there are no features that allow definitive mapping of the origin to the intergenic region upstream of the dnaa gene, as seen in other bacteria. potential pathogenic features bacteria have many ways to produce virulence that reside in the ability to adhere, invade and cause damage to host cells. various strategies of pathogenicity such as cytolysins, toxins and invasins enable other bacteria to produce infection. in mycoplasma species no such typical primary virulence genes have been found. mycoplasmas seem rather to use intrinsic metabolic and catabolic functions to cause disease in the affected host and to ensure the microbe's survival. our efforts to identify genes involved in the pathogenicity of mycoplasma conjunctivae were concentrated on the one hand, try to find those primary virulence genes, toxins principally, rare in other mycoplasmas. on the other hand, on metabolic pathways that has been proposed by studies carried out in other mycoplasmas. glycerol pathway we found using manual blastp queries by an expert, the genes for a glycerol-3-phosphate dehydrogenase, a glycerol kinase, a glycerol uptake facilitator protein and an abc transporter system that are implicated in the glycerol metabolism producing cell damage, inflammation and disease in mycoplasma mycoides subsp. mycoides small colony. the pathway starts with the assimilation of glycerol by the abc glycerol transporter. afterwards, the glycerol is phosphorylated into glycerol-3-phosphate, then oxidized by glpo in presence of o into dihydroxyactone-phosphate producing one molecule of h2o h2o is released directly inside the host cells by the transmembrane glpo protein leading to cell death. the absence of any gene having a catalase or dismutase activity favors this hypothesis. the identification of those genes in mycoplasma conjunctivae constitutes an important discovery given that a relationship between the glycerol metabolism and cytotoxicity is established in the laboratory. further work to validate this hypothesis in m. conjunctivae is required and has been started in collaboration with a laboratory of the institute for veterinary bacteriology. toxins toxins constitute an important type of virulence factors in several bacteria. thereby, we searched for toxins in m. conjunctivae and we found proteins highly similar with toxins of treponema hyodysenteriae. those proteins are hemolysin a, hemolysin b and hemolysin c. the genes are scattered on the genome. those proteins are present in other mycoplasmas, particularly m. hyopneumoniae and m. capricolum, and even if in those species, these toxins are not essential for pathogenicity mechanisms, it can not be excluded that these toxins contribute to the pathogenicity of m. conjunctivae. is elements insertion sequences are short dna elements that function as simple transposable elements by coding for proteins implicated in the transposition activity. transposase and other regulatory protein are the proteins generally coded by is elements: the transposase catalyses the enzymatic reaction allowing the is to move. regulatory proteins act by enhancing or inhibiting the transposition activity. the coding region in an insertion sequence is usually flanked by inverted repeats. we found several genes coding for complete or partial transposases. an is insertion element has particularly brought our attention. is elements belong to is family are prevalent in other mycoplasmas and are the only that have been demonstrated directly to undergo autonomous transposition. interestingly a transposase for one is insertion elements is followed by homologues of a methylase hpai and a type ii restriction enzyme hpai from haemophilus parainfluenzae forming a restriction-methylation cassette. the hypothesis of a horizontal transfer from h. parainfluenza to m. conjunctivae was formulated. we evaluated the g+c content of this cassette, but we did not observe a higher g+c content inside the cassette compared to the surrounding area. if the g+c content inside the cassette would be different from that of m. conjunctivae and similar to that of h. parainfluenzae it could constitute an evidence of the transfer. comparative genomics the list of proteins was classified and compared to other mycoplasma genomes as shown in table the main difference with other mycoplasma is an apparently low carbohydrate and transport metabolism that could explain the need for a strong glycerol pathway, as well as the large number of hypothetical proteins probably due to the fully automatic annotation process. mycoplasma conjunctivae discussion mycoplasma conjunctivae is the fourteenth genome of a mycoplasma species that has been fully sequenced. phylogenetically, the closest relative among the sequenced mycoplasmas is m. hyopneumoniae reflected by the high similarity of most of the proteins identified in m. conjunctivae. the analysis of m. conjunctivae genome features, describes this organism as a typical mycoplasma, with a genome size and a g+c content within the range of other mycoplasma genomes. the comparison of mycoplasma genome sizes demonstrates that the sequence length is variable not only within the same genus but even among strains of the same species as shown in table even if we do not know the final size of the genome, we expect a chromosome length of about 900' bp, size almost similar with the genome size of m. hyopneumoniae. mycoplasma conjunctivae globally the mycoplasma genomes have a characteristically low g+c content within the range of to mol%. the highest g+c content found in m. pneumoniae and the lowest in m. capricolum. regarding mycoplasma conjunctivae, g+c content has a typical value of about 29%. the codon usage is similar to that of m. hyopneumoniae and opposite to that of m. capricolum and m. mycoides. the presence of repeats across the genome was the principal difficulty for finishing the genome assembly. insertion sequence elements are reported in the majority of mycoplasmas and in m. conjunctivae we found transposases for is-elements in the genome. some of those transposases genes are complete sequences and some other are fragmented showing a predicted length of less than bp. since those insertion elements are nearly identical they created difficulties for assembling the genome. the findings highlighted by this project, principally the glycerol pathway, require further experimental confirmation. in particular, the hypothesis for damaging the host cells by the glycerol metabolism need to be confirmed by demonstrating the localization of glpo in the membrane and the release of h2o outside the cell. if this hypothesis can be verified, the possibility to block at any stage the glycerol pathway could constitute a candidate target for controlling the disease.  in conclusion, we created an automatic pipeline to annotate a prokaryotic genome sequence using various tools for the prediction and the identification of the genes. this pipeline is customized for handling sequences of mycoplasma species. we deposited the mycoplasma conjunctivae genome fully annotated in the embl database. data stored into our local database can be searched and genome can be visualized through our website. analysis of annotated genes gives new insights about potential mechanisms of pathogenicity as well as the possibility to go deeper into the knowledge of mycoplasma conjunctivae and the ikc disease and opens the way to finding methods to prevent m. conjunctivae infections of domestic animals as reservoir for this pathogen and hence prevent ikc in wild animals. competing interests the authors declare that they have no competing interests. authors' contributions jf proposed and supported the project. spcc created the pipeline and did the analysis. maq did the genomic library. ts and maq sequenced the clones. ts, gw and cw did the assembly. spcc and lf wrote the manuscript. 
the annotated sequence is deposited in the embl database and uploaded into the mollicutes database molligen allowing for comparative genomics. the genome was sequenced using a combined technique of gs-flx and sanger sequencing, and annotated by an automatic pipeline that we designed using several tools interconnected via perl scripts. the mollicute mycoplasma conjunctivae is the etiological agent leading to infectious keratoconjunctivitis 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. we show that our automatic pipeline allows for annotating a complete mycoplasma genome and present several examples of analysis in search for biological targets.
 the mollicute mycoplasma conjunctivae is the etiological agent leading to infectious keratoconjunctivitis 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. methods the genome was sequenced using a combined technique of gs-flx 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.  the annotated sequence is deposited in the embl database and uploaded into the mollicutes database molligen allowing for comparative genomics.  we show that our automatic pipeline allows for annotating a complete mycoplasma genome and present several examples of analysis in search for biological targets. 18– september european molecular biology network conference 2008: 20th anniversary celebration martina franca, italy 