The relaxing ori-ter balance of Mycoplasma genomes

Mycoplasma are wall-less bacteria with small genomes, which are thought to have resulted from massive genome reductive processes, during which the ori-ter balance may be disrupted. For technical difficulties, ori and ter have been located only in a few Mycoplasma strains. Using the Z curve method, we were able to locate turning points on the Mycoplasma genomes, with the minimum and maximum points co-locating with ori or ter in the reference genomes. Assuming Z curve correctly located ori and ter, we calculated the distances from ori to ter in both directions on the circular genome and calculated the ori-ter balance status. The Mycoplasma genomes were not balanced, possibly as a result of close association of Mycoplasma with hosts, where there would be no other microbes for Mycoplasma to compete with for nutrients, so fastest possible growth related to balanced genomes might not be needed by Mycoplasma, leading to a relaxing ori-ter balance.     

Analysis of distribution indicates diverse functions of simple sequence repeats in Mycoplasma genomes

Simple sequence repeats (SSRs) composed of extensive tandem iterations of a single nucleotide or a short oligonucleotide are rare in most bacterial genomes, but they are common among Mycoplasma. Some of these repeats act as contingency loci in association with families of surface antigens. By contraction or expansion during replication, these SSRs increase genetic variance of the population and facilitate avoidance of the immune response of the host. Occurrence and distribution of SSRs are analyzed in complete genomes of 11 Mycoplasma and 3 related Mollicutes in order to gain insights into functional and evolutionary diversity of the SSRs in Mycoplasma. The results revealed an unexpected variety of SSRs with respect to their distribution and composition and suggest that it is unlikely that all SSRs function as contingency loci or recombination hot spots. Various types of SSRs are most abundant in Mycoplasma hyopneumoniae, whereas Mycoplasma penetrans, Mycoplasma mobile, and Mycoplasma synoviae do not contain unusually long SSRs. Mycoplasma hyopneumoniae and Mycoplasma pulmonis feature abundant short adenine and thymine runs periodically spaced at 11 and 12 bp, respectively, which likely affect the supercoiling propensities of the DNA molecule. Physiological roles of long adenine and thymine runs in M. hyopneumoniae appear independent of location upstream or downstream of genes, unlike contingency loci that are typically located in protein-coding regions or upstream regulatory regions. Comparisons among 3 M. hyopneumoniae strains suggest that the adenine and thymine runs are rarely involved in genome rearrangements. The results indicate that the SSRs in the Mycoplasma genomes play diverse roles, including modulating gene expression as contingency loci, facilitating genome rearrangements via recombination, affecting protein structure and possibly protein-protein interactions, and contributing to the organization of the DNA molecule in the cell.     

Causal analysis of CpG suppression in the Mycoplasma genome

Some bacterial genomes are known to have low CpG dinucleotide frequencies. While their causes are not clearly understood, the frequency of CpG is suppressed significantly in the genome of Mycoplasma genitalium, but not in that of Mycoplasma pneumoniae. We compared orthologous gene pairs of the two closely related species to analyze CpG substitution patterns between these two genomes. We also divided genome sequences into three regions: protein-coding, noncoding, and RNA-coding, and obtained the CpG frequencies for each region for each organism. It was found that the observed/expected ratio of CpG dinucleotides is low in both the protein-coding and noncoding regions; while that ratio is in the normal range in the RNA-coding region. Our results indicate that CpG suppression of the Mycoplasma genome is not caused by (1) biased usage amino acid; (2) biased usage of synonymous codon; or (3) methylation effects by the CpG methyltransferase in the genomes of their hosts. Instead, we consider it likely that a certain global pressure, such as genome-wide pressure for the advantages of DNA stability or replication, has the effect of decreasing CpG over the entire genome, which, in turn, resulted in the biased codon usage.     

Sequence analysis of 56 kb from the genome of the bacterium Mycoplasma pneumoniae comprising the dnaA region, the atp operon and a cluster of ribosomal protein genes.

To sequence the entire 800 kilobase pair genome of the bacterium Mycoplasma pneumoniae, a plasmid library was established with contained the majority of the EcoR1 fragments from M.pneumoniae. The EcoR1 fragments were subcloned from an ordered cosmid library comprising the complete M.pneumoniae genome. Individual plasmid clones were sequenced in an ordered fashion mainly by primer walking. We report here the initial results from the sequence analysis of -56 kb comprising the dnaA region as a potential origin of replication, the ATPase operon and a region coding for a cluster of ribosomal protein genes. The data were compared with the corresponding genes/operons from Bacillus subtilis, Escherichia coli, Mycoplasma capricolum and Mycoplasma gallisepticum.     

Genome size of Mycoplasma genitalium.

The genome size of Mycoplasma genitalium was determined by using restriction enzymes that infrequently cut its DNA. The calculated value of 577 to 590 kilobases is one-fourth smaller than the genome of Mycoplasma pneumoniae, which is considered among the smallest genomes of self-replicating organisms.     

Identification of a ribonuclease H gene in both Mycoplasma genitalium and Mycoplasma pneumoniae by a new method for exhaustive identification of ORFs in the complete genome sequences

Exhaustive identification of open reading frames in complete genome sequences is a difficult task. It is possible that important genes are missed. In our efforts to reanalyze the intergenic regions of Mycoplasma genitalium and Mycoplasma pneumoniae, we have newly identified a number of new open reading frames (ORFs) in both M. genitalium and M. pneumoniae. The most significant identification was that of a ribonuclease H enzyme in both species which until now has not been identified or assumed absent and interpreted as such. In this paper we discuss the biological importance of RNase H and its evolutionary implication. We also stress the usefulness of our method for identifying new ORFs by reanalyzing intergenic regions of existing ORFs in complete genome sequences.     

The influence of regulatory elements on Mycoplasma hyopneumoniae 7448 transcriptional response during oxidative stress and heat shock

Molecular Biology Reports - The comprehension of genome organization and gene modulation is essential for understanding pathogens’ infection mechanisms. Mycoplasma hyopneumoniae 7448 genome...     

Complete nucleotide sequence of the mycoplasma virus P1 genome

Mycoplasma virus P1 is one of only four viruses isolated from the genus Mycoplasma. The host for P1, Mycoplasma pulmonis, possesses complex, phase-variable restriction and modification enzymes and the Vsa family of phase-variable surface proteins. The ability of P1 virus to infect host cells is influenced by these phase-variable systems, rendering P1 a valuable tool for assessing host properties. The double-stranded P1 DNA genome was sequenced (11,660 bp) and 11 ORFs were identified. The predicted P1 DNA polymerase is similar to that of phages that are known to have terminal protein (TP) attached to the 5' end of their genome, consistent with previous studies indicating that P1 DNA has covalently attached TP. Most of the other predicted P1 proteins have little sequence similarity to known proteins, and P1 virus is unrelated to the other mycoplasma virus, MAV1, for which the genome sequence is known. One of the predicted P1 proteins, the ORF 8 gene product, contains a repetitive collagen-like motif characteristic of some bacteriophage tail fiber proteins and is a candidate for interacting with the Vsa proteins.     

Chaperone networks in bacteria: analysis of protein homeostasis in minimal cells

The prevention of aberrant behavior of proteins is fundamental to cellular life. Protein homeostatic processes are present in cells to stabilize protein conformations, refold misfolded proteins, and degrade proteins that might be detrimental to the cell. Molecular chaperones and proteases perform a major role in these processes. In bacteria, the main cytoplasmic components involved in protein homeostasis include the chaperones trigger factor, DnaK/DnaJ/GrpE, GroEL/GroES, HtpG, as well as ClpB and the proteases ClpXP, ClpAP, HslUV, Lon, and FtsH. Based on recent genome sequencing efforts, it was surprising to find that the Mycoplasma, a genus proposed to include a minimal form of cellular life, do not contain certain major members of the protein homeostatic network, including GroEL/GroES. We propose that, in mycoplasmas, there has been a fundamental shift towards favoring processes that promote protein degradation rather than protein folding. The arguments are based on two different premises: (1) the regulation of stress response in Mycoplasma and (2) the unique characteristics of the Mycoplasma proteome.     

Identification of a plant-derived mollicute as a strain of an avian pathogen, Mycoplasma iowae, and its implications for mollicute taxonomy.

Strain PPAV, a filamentous but nonhelical mollicute, was isolated from aborted apple seeds in France in late 1979. This organism grew well in SP-4 broth, fermented glucose, and required sterol for growth, and most of its properties suggested that it belonged to the genus Mycoplasma. However, it was serologically distinct; in addition, unlike other Mycoplasma species, genome measurements consistently yielded values of about 1,000 MDa (ca. 1,500 kbp), and the organism had a growth temperature optimum of 43 degrees C. A comparison of strain PPAV 16S rRNA sequences with those of other mollicutes revealed a high degree of sequence similarity to a strain of Mycoplasma iowae, which is commonly encountered in poultry. This relationship was confirmed by performing a restriction endonuclease pattern analysis and DNA-DNA hybridization tests. The genome size of type strain 695 of M. iowae was determined to be about 1,000 MDa (1,500 kbp) by renaturation kinetics, a value which is much higher than any other value known in the genus. Additional measurements by pulsed-field gel electrophoresis yielded values of 1,300 kbp for both strain PPAV and M. iowae. Subsequent phenotypic comparisons supported this relationship. Serologic tests with strain PPAV and other strains of M. iowae confirmed the findings of other investigators that this species is serologically heterogeneous. The high optimum temperature for growth of strain PPAV was also shared by a number of M. iowae isolates. Genome size is an inappropriate character for taxonomic assignment to the family Mycoplasmataceae because strain PPAV and other established species in this family are now known to have genomes ranging in size from 1,000 to 1,400 kbp.     

Evolutionary genetics of <Emphasis Type="Italic">Carpodacus mexicanus</Emphasis>, a recently colonized host of a bacterial pathogen, <Emphasis Type="Italic">Mycoplasma gallisepticum</Emphasis>

We present molecular data documenting how introduction to the eastern United States and an epizootic involving a bacterial pathogen has affected the genetic diversity of house finches, a cardueline songbird. Population bottlenecks during introduction can cause loss of genetic variation and may negatively affect a population's ability to adapt to novel stressors such as disease. Although a genome-wide survey using Amplified Fragment Length Polymorphism (AFLP) markers suggests little loss of genetic diversity in introduced populations, an epizootic of bacterial Mycoplasma has nonetheless caused dramatic declines in the eastern US population. Sequence analysis of a candidate gene for pathogen resistance in the Major Histocompatibity Complex (MHC) in pre- and post-epizootic population samples reveals allele frequency shifts since introduction of the pathogen, but similar shifts are also observed in control populations not exposed to the bacteria, and in a neutral non-coding locus. Expression studies using a novel subtractive hybridization approach indicate decreased expression of the class II MHC locus upon exposure to Mycoplasma, a pattern also seen in MHC class I loci in mice infected with cytomegalovirus and consistent with manipulation of the finch immune system by Mycoplasma. These results will be further expanded using experimental studies as well as examination of evolution of the pathogen genome itself.     

Complete Genome Sequence of Mycoplasma wenyonii Strain Massachusetts

Mycoplasma wenyonii is a hemotrophic mycoplasma that causes acute and chronic infections in cattle. Here, we announce the first complete genome sequence of this organism. The genome is a single circular chromosome with 650,228 bp and G+C% of 33.9. Analyses of M. wenyonii genome will provide insights into its biology.     

The complete genome sequence of Mycoplasma bovis strain Hubei-1

Infection by Mycoplasma bovis (M. bovis) can induce diseases, such as pneumonia and otitis media in young calves and mastitis and arthritis in older animals. Here, we report the finished and annotated genome sequence of M. bovis strain Hubei-1, a strain isolated in 2008 that caused calf pneumonia on a Chinese farm. The genome of M. bovis strain Hubei-1 contains a single circular chromosome of 953,114 bp with a 29.37% GC content. We identified 803 open reading frames (ORFs) that occupy 89.5% of the genome. While 34 ORFs were Hubei-1 specific, 662 ORFs had orthologs in the M. bovis type strain PG45 genome. Genome analysis validated lateral gene transfer between M. bovis and the Mycoplasma mycoides subspecies mycoides, while phylogenetic analysis found that the closest M. bovis neighbor is Mycoplasma agalactiae. Glycerol may be the main carbon and energy source of M. bovis, and most of the biosynthesis pathways were incomplete. We report that 47 lipoproteins, 12 extracellular proteins and 18 transmembrane proteins are phase-variable and may help M. bovis escape the immune response. Besides lipoproteins and phase-variable proteins, genomic analysis found two possible pathogenicity islands, which consist of four genes and 11 genes each, and several other virulence factors including hemolysin, lipoate protein ligase, dihydrolipoamide dehydrogenase, extracellular cysteine protease and 5'-nucleotidase.     

Genomic maps of some strains within the Mycoplasma mycoides cluster

Genomic restriction maps for the small colony (SC) strains (PG1, KH3J, Gladysdale, and V5) of Mycoplasma mycoides subsp. mycoides (the agent of contagious bovine pleuropneumonia) and for Mycoplasma strain PG50 (classified as bovine serogroup 7), with respective sizes of 1,280, 1,280, 1,260, 1,230, and 1,040 kbp, were compared with the map (1,200 kbp) for a large colony strain (Y goat) of M. mycoides subsp. mycoides. The number and order of all mapped restriction sites were fully conserved in the SC genomes, as were the approximate positions of mapped loci. A number of these restriction sites in the Y genome and some, but fewer, in the PG50 genome appeared to be conserved. The SC and large colony strains shared conservation in the relative positions of the mapped loci, except for rpoC.     

Prevalence of novel repeat sequences in and around the P1 operon in the genome of Mycoplasma pneumoniae

The presence of numerous different repetitive elements in the genome of Mycoplasma pneumoniae has been documented by several laboratories. One which we previously identified, denoted as SDC1, has now been further characterized, verified to be distinct from those discussed in previous publications and shown to lack homology to several other species of Mycoplasma when tested under our stringency conditions. As many as eight versions of the SDC1-type repeat, which is more than 400 bp long, are scattered throughout the genome of M. pneumoniae. The prototype for SDC1 is found within a gene encoding a putative 130-kDa membrane-binding protein lying just downstream from the gene encoding the cytadhesin protein P1. In fact, all of the reported M. pneumoniae repetitive elements have at least one representative either within or adjacent to the P1 operon; many if not all of these lie within open reading frames. The function of these repetitive elements is still unclear.     

Remarkable sequence signatures in archaeal genomes

Complete archaeal genomes were probed for the presence of long (> or = 25 bp) oligonucleotide repeats (words). We detected the presence of many words distributed in tandem with narrow ranges of periodicity (i.e., spacer length between repeats). Similar words were not identified in genomes of non-archaeal species, namely Escherichia coli, Bacillus subtilis, Haemophilus influenzae, Mycoplasma genitalium and Mycoplasma pneumoniae. BLAST similarity searches against the GenBank nucleotide sequence database revealed that these words were archaeal species-specific, indicating that they are of a signature character. Sequence analysis and genome viewing tools showed these repeats to be restricted to non-coding regions. Thus, archaea appear to possess a non-coding genomic signature that is absent in bacterial species. The identification of a species-specific genomic signature would be of great value to archaeal genome mapping, evolutionary studies and analyses of genome complexity.     

Genomic repeats, genome plasticity and the dynamics of Mycoplasma evolution

Mycoplasmas evolved by a drastic reduction in genome size, but their genomes contain numerous repeated sequences with important roles in their evolution. We have established a bioinformatic strategy to detect the major recombination hot-spots in the genomes of Mycoplasma pneumoniae, Mycoplasma genitalium, Ureaplasma urealyticum and Mycoplasma pulmonis. This allowed the identification of large numbers of potentially variable regions, as well as a comparison of the relative recombination potentials of different genomic regions. Different trends are perceptible among mycoplasmas, probably due to different functional and structural constraints. The largest potential for illegitimate recombination in M.pulmonis is found at the vsa locus and its comparison in two different strains reveals numerous changes since divergence. On the other hand, the main M.pneumoniae and M.genitalium adhesins rely on large distant repeats and, hence, homologous recombination for variation. However, the relation between the existence of repeats and antigenic variation is not necessarily straightforward, since repeats of P1 adhesin were found to be anti-correlated with epitopes recognized by patient antibodies. These different strategies have important consequences for the structures of genomes, since large distant repeats correlate well with the major chromosomal rearrangements. Probably to avoid such events, mycoplasmas strongly avoid inverse repeats, in comparison to co-oriented repeats.     

Comparative genomics analysis of Mycoplasma capricolum subsp. capripneumoniae 87001

Mycoplasma capricolum subsp. capripneumoniae (Mccp), belongs to Mycoplasma mycoides cluster and is a causal pathogen of contagious caprine pleuropneumonia (CCPP). This paper presents the complete annotated genome sequence of Mccp Strain 87001—a strain that was isolated from pneumonia affected goats on a farm in China, and comparative genomics analysis of five Mccp genomes in addition to comparative genomics within Mycoplasma mycoides cluster. The Mccp strain 87001 genome consists of a single circular chromosome 1017333 bp in length and encodes 898 open reading frames (orfs) averaging 944 bp in length. Fifty eight potential virulence genes were identified, including variable surface lipoproteins, hemolysin A, and P60 surface lipoprotein. Comparative genomic analysis revealed eight virulence genes and four extracellular genes which remained unchanged in five Mccp genomes for forty years, which can be used as potential target for drug development and vaccine design. We revealed 183 Mccp unique genes as markers to distinguish Mccp with other mycoplasma strains from goats, and different virulence factors contributing to host specificity and different syndrome of bovine pathogens and caprine pathogens.     

Complete genome sequence of Mycoplasma hyopneumoniae strain 168

Mycoplasma hyopneumoniae strain 168, a pathogenic strain prevalent in China, was isolated in 1974. Although this strain has been widespread for a long time, the genome sequence had not been determined. Here, we announce the complete genome sequence of M. hyopneumoniae strain 168.     

The genome of Erysipelothrix rhusiopathiae, the causative agent of swine erysipelas, reveals new insights into the evolution of firmicutes and the organism's intracellular adaptations

Erysipelothrix rhusiopathiae is a Gram-positive bacterium that represents a new class, Erysipelotrichia, in the phylum Firmicutes. The organism is a facultative intracellular pathogen that causes swine erysipelas, as well as a variety of diseases in many animals. Here, we report the first complete genome sequence analysis of a member of the class Erysipelotrichia. The E. rhusiopathiae genome (1,787,941 bp) is one of the smallest genomes in the phylum Firmicutes. Phylogenetic analyses based on the 16S rRNA gene and 31 universal protein families suggest that E. rhusiopathiae is phylogenetically close to Mollicutes, which comprises Mycoplasma species. Genome analyses show that the overall features of the E. rhusiopathiae genome are similar to those of other Gram-positive bacteria; it possesses a complete set of peptidoglycan biosynthesis genes, two-component regulatory systems, and various cell wall-associated virulence factors, including a capsule and adhesins. However, it lacks many orthologous genes for the biosynthesis of wall teichoic acids (WTA) and lipoteichoic acids (LTA) and the dltABCD operon, which is responsible for d-alanine incorporation into WTA and LTA, suggesting that the organism has an atypical cell wall. In addition, like Mollicutes, its genome shows a complete loss of fatty acid biosynthesis pathways and lacks the genes for the biosynthesis of many amino acids, cofactors, and vitamins, indicating reductive genome evolution. The genome encodes nine antioxidant factors and nine phospholipases, which facilitate intracellular survival in phagocytes. Thus, the E. rhusiopathiae genome represents evolutionary traits of both Firmicutes and Mollicutes and provides new insights into its evolutionary adaptations for intracellular survival.