Amelioration of Bacterial Genomes: Rates of Change and Exchange

Although bacterial species display wide variation in their overall GC contents, the genes within a particular species' genome are relatively similar in base composition. As a result, sequences that are novel to a bacterial genome—i.e., DNA introduced through recent horizontal transfer—often bear unusual sequence characteristics and can be distinguished from ancestral DNA. At the time of introgression, horizontally transferred genes reflect the base composition of the donor genome; but, over time, these sequences will ameliorate to reflect the DNA composition of the new genome because the introgressed genes are subject to the same mutational processes affecting all genes in the recipient genome. This process of amelioration is evident in a large group of genes involved in host-cell invasion by enteric bacteria and can be modeled to predict the amount of time required after transfer for foreign DNA to resemble native DNA. Furthermore, models of amelioration can be used to estimate the time of introgression of foreign genes in a chromosome. Applying this approach to a 1.43-megabase continuous sequence, we have calculated that the entire Escherichia coli chromosome contains more than 600 kb of horizontally transferred, protein-coding DNA. Estimates of amelioration times indicate that this DNA has accumulated at a rate of 31 kb per million years, which is on the order of the amount of variant DNA introduced by point mutations. This rate predicts that the E. coli and Salmonella enterica lineages have each gained and lost more than 3 megabases of novel DNA since their divergence. Received: 7 July 1996 / Accepted: 27 September 1996     

A survey of the Mycoplasma genitalium genome by using random sequencing.

A total of 508 random clones from five Mycoplasma genitalium genomic libraries were partially sequenced and analyzed. This resulted in the identification of 291 unique contigs. Sequence information from these clones (100,993 nucleotides), representing approximately 17% of this pathogen's genome, was analyzed by comparison to the DNA and protein sequence data bases. The frequency with which clones could be identified, by virtue of possessing homology to another data base entry, was 46%. Sequence analysis indicated the following. (i) The M. genitalium genome contains many genes involved in various metabolic processes. (ii) Repetitive DNA may comprise as much as 4% of this genome. (iii) The MgPa adhesin gene may be the result of horizontal transfer from an unknown origin. (iv) Not all dinucleotide pairs are present in this genome at the expected frequency. (v) This genome potentially encodes approximately 390 proteins and makes very efficient use of its limited amount of DNA. In addition, this study allowed us to estimate the number of genes involved with various cellular functions.     

Analysis of genomic characters reveals that four distinct gene clusters are correlated with different functions in Burkholderia cenocepacia AU 1054

Possessing three circular chromosomes is a distinct genomic characteristic of Burkholderia cenocepacia AU 1054, a clinically important pathogen in cystic fibrosis. In this study, base composition, codon usage and functional role category were analyzed in the B. cenocepacia AU 1054 genome. Although no bias in the base and codon usage was detected between any two chromosomes, function differences did exist in the genes of each chromosome. Similar base composition and differential functional role categories indicated that genes on these three chromosomes were relatively stable and that a proper division of labor was established. Based on variations in the base or codon usage, four small gene clusters were observed in all of the genes. Multivariate analysis revealed that protein hydrophobicity played a predominant role in shaping base usage bias, while horizontal gene transfer and the gene expression level were the two most important factors that affected the codon usage bias. Interestingly, we also found that these gene clusters were correlated with different biological functions: (i) 45 pyrimidine-leading-codon preferred genes were predominantly involved in regulatory function; (ii) most drug resistance-related genes involved in 826 genes that coding for hydrophobic proteins; (iii) most of the 111 horizontal transfer genes were responsible for genomic plasticity; and (iv) 73 highly expressed genes (predicted by their codon adaptation index values) showed environmental adaptation to cystic fibrosis. Our results showed that genes with base or codon usage bias were affected by mutational pressure and natural selection, and their functions could contribute to drug assistance and transmissible activity in B. cenocepacia.     

Strain-specific genes of Helicobacter pylori: genome evolution driven by a novel type IV secretion system and genomic island transfer

The availability of multiple bacterial genome sequences has revealed a surprising extent of variability among strains of the same species. The human gastric pathogen Helicobacter pylori is known as one of the most genetically diverse species. We have compared the genome sequence of the duodenal ulcer strain P12 and six other H. pylori genomes to elucidate the genetic repertoire and genome evolution mechanisms of this species. In agreement with previous findings, we estimate that the core genome comprises about 1200 genes and that H. pylori possesses an open pan-genome. Strain-specific genes are preferentially located at potential genome rearrangement sites or in distinct plasticity zones, suggesting two different mechanisms of genome evolution. The P12 genome contains three plasticity zones, two of which encode type IV secretion systems and have typical features of genomic islands. We demonstrate for the first time that one of these islands is capable of self-excision and horizontal transfer by a conjugative process. We also show that excision is mediated by a protein of the XerD family of tyrosine recombinases. Thus, in addition to its natural transformation competence, conjugative transfer of genomic islands has to be considered as an important source of genetic diversity in H. pylori.     

The Mitochondrial Genome of Paramphistomum cervi (Digenea), the First Representative for the Family Paramphistomidae

We determined the complete mitochondrial DNA (mtDNA) sequence of a fluke, Paramphistomum cervi (Digenea: Paramphistomidae). This genome (14,014 bp) is slightly larger than that of Clonorchis sinensis (13,875 bp), but smaller than those of other digenean species. The mt genome of P. cervi contains 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and 2 non-coding regions (NCRs), a complement consistent with those of other digeneans. The arrangement of protein-coding and ribosomal RNA genes in the P. cervi mitochondrial genome is identical to that of other digeneans except for a group of Schistosoma species that exhibit a derived arrangement. The positions of some transfer RNA genes differ. Bayesian phylogenetic analyses, based on concatenated nucleotide sequences and amino-acid sequences of the 12 protein-coding genes, placed P. cervi within the Order Plagiorchiida, but relationships depicted within that order were not quite as expected from previous studies. The complete mtDNA sequence of P. cervi provides important genetic markers for diagnostics, ecological and evolutionary studies of digeneans.     

The 203 kbp Mitochondrial Genome of the Phytopathogenic Fungus Sclerotinia borealis Reveals Multiple Invasions of Introns and Genomic Duplications

Here we report the complete sequence of the mitochondrial (mt) genome of the necrotrophic phytopathogenic fungus Sclerotinia borealis, a member of the order Helotiales of Ascomycetes. The 203,051 bp long mtDNA of S. borealis represents one of the largest sequenced fungal mt genomes. The large size is mostly determined by the presence of mobile genetic elements, which include 61 introns. Introns contain a total of 125,394 bp, are scattered throughout the genome, and are found in 12 protein-coding genes and in the ribosomal RNA genes. Most introns contain complete or truncated ORFs that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Integrations of mobile elements are also evidenced by the presence of two regions similar to fragments of inverton-like plasmids. Although duplications of some short genome regions, resulting in the appearance of truncated extra copies of genes, did occur, we found no evidences of extensive accumulation of repeat sequences accounting for mitochondrial genome size expansion in some other fungi. Comparisons of mtDNA of S. borealis with other members of the order Helotiales reveal considerable gene order conservation and a dynamic pattern of intron acquisition and loss during evolution. Our data are consistent with the hypothesis that horizontal DNA transfer has played a significant role in the evolution and size expansion of the S. borealis mt genome.     

The genome of ‘Candidatus Phytoplasma solani’ strain SA-1 is highly dynamic and prone to adopting foreign sequences

Bacteria of the genus ‘Candidatus Phytoplasma’ are uncultivated intracellular plant pathogens transmitted by phloem-feeding insects. They have small genomes lacking genes for essential metabolites, which they acquire from either plant or insect hosts. Nonetheless, some phytoplasmas, such as ‘Ca. P. solani’, have broad plant host range and are transmitted by several polyphagous insect species. To understand better how these obligate symbionts can colonize such a wide range of hosts, the genome of ‘Ca. P. solani’ strain SA-1 was sequenced from infected periwinkle via a metagenomics approach. The de novo assembly generated a draft genome with 19 contigs totalling 821,322 bp, which corresponded to more than 80% of the estimated genome size. Further completion of the genome was challenging due to the high occurrence of repetitive sequences. The majority of repeats consisted of gene arrangements characteristic of phytoplasma potential mobile units (PMUs). These regions showed variation in gene orders intermixed with genes of unknown functions and lack of similarity to other phytoplasma genes, suggesting that they were prone to rearrangements and acquisition of new sequences via recombination. The availability of this high-quality draft genome also provided a foundation for genome-scale genotypic analysis (e.g., average nucleotide identity and average amino acid identity) and molecular phylogenetic analysis. Phylogenetic analyses provided evidence of horizontal transfer for PMU-like elements from various phytoplasmas, including distantly related ones. The ‘Ca. P. solani’ SA-1 genome also contained putative secreted protein/effector genes, including a homologue of SAP11, found in many other phytoplasma species.     

Positive Selection and Horizontal Gene Transfer in the Genome of a Male-Killing Wolbachia

Wolbachia are a genus of widespread bacterial endosymbionts in which some strains can hijack or manipulate arthropod host reproduction. Male killing is one such manipulation in which these maternally transmitted bacteria benefit surviving daughters in part by removing competition with the sons for scarce resources. Despite previous findings of interesting genome features of microbial sex ratio distorters, the population genomics of male-killers remain largely uncharacterized. Here, we uncover several unique features of the genome and population genomics of four Arizonan populations of a male-killing Wolbachia strain, wInn, that infects mushroom-feeding Drosophila innubila. We first compared the wInn genome with other closely related Wolbachia genomes of Drosophila hosts in terms of genome content and confirm that the wInn genome is largely similar in overall gene content to the wMel strain infecting D. melanogaster. However, it also contains many unique genes and repetitive genetic elements that indicate lateral gene transfers between wInn and non-Drosophila eukaryotes. We also find that, in line with literature precedent, genes in the Wolbachia prophage and Octomom regions are under positive selection. Of all the genes under positive selection, many also show evidence of recent horizontal transfer among Wolbachia symbiont genomes. These dynamics of selection and horizontal gene transfer across the genomes of several Wolbachia strains and diverse host species may be important underlying factors in Wolbachia’s success as a male-killer of divergent host species.     

An Uncultivated Nitrate-Reducing Member of the Genus Herminiimonas Degrades Toluene

Stable isotope probing (SIP) is a cultivation-free methodology that provides information about the identity of microorganisms participating in assimilatory processes in complex communities. In this study, a Herminiimonas-related bacterium was identified as the dominant member of a denitrifying microcosm fed [¹³C]toluene. The genome of the uncultivated toluene-degrading bacterium was obtained by applying pyrosequencing to the heavy DNA fraction. The draft genome comprised ∼3.8 Mb, in 131 assembled contigs. Metabolic reconstruction of aromatic hydrocarbon (toluene, benzoate, p-cresol, 4-hydroxybenzoate, phenylacetate, and cyclohexane carboxylate) degradation indicated that the bacterium might specialize in anaerobic hydrocarbon degradation. This characteristic is novel for the order Burkholderiales within the class Betaproteobacteria. Under aerobic conditions, the benzoate oxidation gene cluster (BOX) system is likely involved in the degradation of benzoate via benzoyl coenzyme A. Many putative genes for aromatic hydrocarbon degradation were closely related to those in the Rhodocyclaceae (particularly Aromatoleum aromaticum EbN1) with respect to organization and sequence similarity. Putative mobile genetic elements associated with these catabolic genes were highly abundant, suggesting gene acquisition by Herminiimonas via horizontal gene transfer.     

<Emphasis Type="Italic">In silico</Emphasis> prediction of horizontal gene transfer in <Emphasis Type="Italic">Streptococcus thermophilus</Emphasis>

A combination of gene loss and acquisition through horizontal gene transfer (HGT) is thought to drive Streptococcus thermophilus adaptation to its niche, i.e. milk. In this study, we describe an in silico analysis combining a stochastic data mining method, analysis of homologous gene distribution and the identification of features frequently associated with horizontally transferred genes to assess the proportion of the S. thermophilus genome that could originate from HGT. Our mining approach pointed out that about 17.7% of S. thermophilus genes (362 CDSs of 1,915) showed a composition bias; these genes were called ‘atypical’. For 22% of them, their functional annotation strongly support their acquisition through HGT and consisted mainly in genes encoding mobile genetic recombinases, exopolysaccharide (EPS) biosynthesis enzymes or resistance mechanisms to bacteriophages. The distribution of the atypical genes in the Firmicutes phylum as well as in S. thermophilus species was sporadic and supported the HGT prediction for more than a half (52%, 189). Among them, 46 were found specific to S. thermophilus. Finally, by combining our method, gene annotation and sequence specific features, new genome islands were suggested in the S. thermophilus genome.     

Genomic Features of Lactococcus lactis IO-1, a Lactic Acid Bacterium That Utilizes Xylose and Produces High Levels of L-Lactic Acid

Lactococcus lactis IO-1 (JCM7638) produces L-lactic acid predominantly when grown at high xylose concentrations, and its utilization is highly desired in the green plastics industry. Therefore it is worthwhile studying its genomic traits. In this study, we focused on (i) genes of possible horizontal transfer derivation (prophages, the nisin-sucrose transposon, and several restriction-modification systems), and (ii) genes for the synthetic pathways of amino acids and vitamins in the IO-1 genome. In view of the results of this analysis, we consider their meanings in strain IO-1.     

The Aryl Hydrocarbon Receptor-interacting Protein (AIP) Is Required for Dioxin-induced Hepatotoxicity but Not for the Induction of the Cyp1a1 and Cyp1a2 Genes

The aryl hydrocarbon receptor (AHR) plays an essential role in the toxic response to environmental pollutants such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (dioxin), in the adaptive up-regulation of xenobiotic metabolizing enzymes, and in hepatic vascular development. In our model of AHR signaling, the receptor is found in a cytosolic complex with a number of molecular chaperones, including Hsp90, p23, and the aryl hydrocarbon receptor-interacting protein (AIP), also known as ARA9 and XAP2. To understand the role of AIP in adaptive and toxic aspects of AHR signaling, we generated a conditional mouse model where the Aip locus can be deleted in hepatocytes. Using this model, we demonstrate two important roles for the AIP protein in AHR biology. (i) The expression of AIP in hepatocytes is essential to maintain high levels of functional cytosolic AHR protein in the mammalian liver. (ii) Expression of the AIP protein is essential for dioxin-induced hepatotoxicity. Interestingly, classical AHR-driven genes show differential dependence on AIP expression. The Cyp1b1 and Ahrr genes require AIP expression for normal up-regulation by dioxin, whereas Cyp1a1 and Cyp1a2 do not. This differential dependence on AIP provides evidence that the mammalian genome contains more than one class of AHR-responsive genes and suggests that a search for AIP-dependent, AHR-responsive genes may guide us to the targets of the dioxin-induced hepatotoxicity.     

Staphylococcus aureus mobile genetic elements

Among the bacteria groups, most of them are known to be beneficial to human being whereas only a minority is being recognized as harmful. The pathogenicity of bacteria is due, in part, to their rapid adaptation in the presence of selective pressures exerted by the human host. In addition, through their genomes, bacteria are subject to mutations, various rearrangements or horizontal gene transfer among and/or within bacterial species. Bacteria’s essential metabolic functions are generally encoding by the core genes. Apart of the core genes, there are several number of mobile genetic elements (MGE) acquired by horizontal gene transfer that might be beneficial under certain environmental conditions. These MGE namely bacteriophages, transposons, plasmids, and pathogenicity islands represent about 15 % Staphylococcus aureus genomes. The acquisition of most of the MGE is made by horizontal genomic islands (GEI), recognized as discrete DNA segments between closely related strains, transfer. The GEI contributes to the wide spread of microorganisms with an important effect on their genome plasticity and evolution. The GEI are also involve in the antibiotics resistance and virulence genes dissemination. In this review, we summarize the mobile genetic elements of S. aureus.     

Complete Mitochondrial Genome of Echinostoma hortense (Digenea: Echinostomatidae)

Echinostoma hortense (Digenea: Echinostomatidae) is one of the intestinal flukes with medical importance in humans. However, the mitochondrial (mt) genome of this fluke has not been known yet. The present study has determined the complete mt genome sequences of E. hortense and assessed the phylogenetic relationships with other digenean species for which the complete mt genome sequences are available in GenBank using concatenated amino acid sequences inferred from 12 protein-coding genes. The mt genome of E. hortense contained 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 non-coding region. The length of the mt genome of E. hortense was 14,994 bp, which was somewhat smaller than those of other trematode species. Phylogenetic analyses based on concatenated nucleotide sequence datasets for all 12 protein-coding genes using maximum parsimony (MP) method showed that E. hortense and Hypoderaeum conoideum gathered together, and they were closer to each other than to Fasciolidae and other echinostomatid trematodes. The availability of the complete mt genome sequences of E. hortense provides important genetic markers for diagnostics, population genetics, and evolutionary studies of digeneans.     

Genome-wide survey and genetic characteristics of Ophichthus evermanni based on Illumina sequencing platform

Ophichthidae fishes limit to continental shelf of all tropical and subtropical oceans and contain more than 350 species, representing the greatest specialization diversity in the order Anguiliformes. In the present study, we conducted a genome survey sequencing (GSS) analysis of Ophichthus evermanni by Illumina sequencing platform to briefly reveal its genomic characteristics and phylogenetic relationship. The first de novo assembled 1.97 Gb draft genome of O. evermanni was predicted based on K-mer analysis without obvious nucleotide bias. The heterozygosity ratio was 0.70%, and the sequence repeat ratio was calculated to be 43.30%. A total of 9016 putative coding genes were successfully predicted, in which 3587 unigenes were identified by gene ontology (GO) analysis and 4375 unigenes were classified into cluster of orthologous groups for enkaryotic complete genomes (KOG) functional categories. About 2,812,813 microsatellite motifs including mono-, di-, tri-, tetra-, penta- and hexanucleotide motifs were identified, with an occurrence frequency of 23.32%. The most abundant type was dinucleotide repeat motifs, accounting for 49.19% of the total repeat types. The mitochondrial genome, as a byproduct of GSS, was assembled to investigate the evolutionary relationships between O. evermanni and its relatives. Bayesian inference (BI) phylogenetic tree inferring from concatenated 12 protein-coding genes (PCGs) showed complicated relationships among Ophichthidae species, indicating a polyphyletic origin of the family. The results would achieve more thorough genetic information of snake eels and provide a theoretical basis and reference for further genome-wide analysis of O. evermanni.     

Transformation competence and type-4 pilus biogenesis in Neisseria gonorrhoeae – a review

In Neisseria gonorrhoea (Ngo), the processes of type-4 pilus biogenesis and DNA transformation are functionally linked and play a pivotal role in the life style of this strictly human pathogen. The assembly of pili from its main subunit pilin (PilE) is a prerequisite for gonococcal infection since it allows the first contact to epithelial cells in conjunction with the pilus tip-associated PilC protein. While the components of the pilus and its assembly machinery are either directly or indirectly involved in the transport of DNA across the outer membrane, other factors unrelated to pilus biogenesis appear to facilitate further DNA transfer across the murein layer (ComL, Tpc) and the inner membrane (ComA) before the transforming DNA is rescued in the recipient bacterial chromosome in a RecA-dependent manner. Interestingly, PilE is essential for the first step of transformation, i.e., DNA uptake, and is itself also subject to transformation-mediated phase and antigenic variation. This short-term adaptive mechanism allows Ngo to cope with changing micro-environments in the host as well as to escape the immune response during the course of infection. Given the fact that Ngo has no ecological niche other than man, horizontal genetic exchange is essential for a successful co-evolution with the host. Horizontal exchange gives rise to heterogeneous populations harboring clones which better withstand selective forces within the host. Such extended horizontal exchange is reflected by a high genome plasticity, the existence of mosaic genes and a low linkage disequilibrium of genetic loci within the neisserial population. This led to the concept that rather than regarding individual Neisseria species as independent traits, they comprise a collective of species interconnected via horizontal exchange and relying on a common gene pool.     

Disruption of Mycobacterium avium subsp. paratuberculosis-specific genes impairs in vivo fitness

Background

Mycobacterium avium subsp. paratuberculosis (MAP) is an obligate intracellular pathogen that infects many ruminant species. The acquisition of foreign genes via horizontal gene transfer has been postulated to contribute to its pathogenesis, as these genetic elements are absent from its putative ancestor, M. avium subsp. hominissuis (MAH), an environmental organism with lesser pathogenicity. In this study, high-throughput sequencing of MAP transposon libraries were analyzed to qualitatively and quantitatively determine the contribution of individual genes to bacterial survival during infection.

Results

Out of 52384 TA dinucleotides present in the MAP K-10 genome, 12607 had a MycoMarT7 transposon in the input pool, interrupting 2443 of the 4350 genes in the MAP genome (56%). Of 96 genes situated in MAP-specific genomic islands, 82 were disrupted in the input pool, indicating that MAP-specific genomic regions are dispensable for in vitro growth (odds ratio = 0.21). Following 5 independent in vivo infections with this pool of mutants, the correlation between output pools was high for 4 of 5 (R = 0.49 to 0.61) enabling us to define genes whose disruption reproducibly reduced bacterial fitness in vivo. At three different thresholds for reduced fitness in vivo, MAP-specific genes were over-represented in the list of predicted essential genes. We also identified additional genes that were severely depleted after infection, and several of them have orthologues that are essential genes in M. tuberculosis.

Conclusions

This work indicates that the genetic elements required for the in vivo survival of MAP represent a combination of conserved mycobacterial virulence genes and MAP-specific genes acquired via horizontal gene transfer. In addition, the in vitro and in vivo essential genes identified in this study may be further characterized to offer a better understanding of MAP pathogenesis, and potentially contribute to the discovery of novel therapeutic and vaccine targets.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-415) contains supplementary material, which is available to authorized users.     

A genomic perspective on the relationship between the Aquificales and the epsilon-Proteobacteria

The goal of this study was to determine to what extent the Aquificales are related to the epsilon-Proteobacteria. The genome sequence of several members of this group as well as the genome sequence of Aquifex aeolicus are available. In this study we used information extracted from those whole-genome sequences to gain further insights into the relationships between these organisms, including the fraction of shared putative orthologous protein-encoding genes, dinucleotide relative abundance values and the sequences of the 16S rRNA gene and 20 housekeeping genes. The results of our analyses show that it is not straightforward to come to a consistent picture of the phylogenetic position of the order Aquificales but our data clearly show that there is no particularly close relationship between A. aeolicus and the epsilon-Proteobacteria as (i) they do not share more genes with each other than do other distantly related organisms and (ii) they do not share significant sequence similarity in many macromolecules. In addition, there is considerable evidence that confirms the placement of the Aquificales near the root of the bacterial tree.     

The first complete mitochondrial genome from Bostrychus genus (Bostrychus sinensis) and partitioned Bayesian analysis of Eleotridae fish phylogeny

To understand the phylogenetic position of Bostrychus sinensis in Eleotridae and the phylogenetic relationships of the family, we determined the nucleotide sequence of the mitochondrial (mt) genome of Bostrychus sinensis. It is the first complete mitochondrial genome sequence of Bostrychus genus. The entire mtDNA sequence was 16508 bp in length with a standard set of 13 protein-coding genes, 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a noncoding control region. The mitochondrial genome of B. sinensis had common features with those of other bony fishes with respect to gene arrangement, base composition, and tRNA structures. Phylogenetic hypotheses within Eleotridae fish have been controversial at the genus level. We used the mitochondrial cytochrome b (cytb) gene sequence to examine phylogenetic relationships of Eleotridae by using partitioned Bayesian method. When the specific models and parameter estimates were presumed for partitioning the total data, the harmonic mean –lnL was improved. The phylogenetic analysis supported the monophyly of Hypseleotris and Gobiomorphs. In addition, the Bostrychus were most closely related to Ophiocara, and the Philypnodon is also the sister to Microphlypnus, based on the current datasets. Further, extensive taxonomic sampling and more molecular information are needed to confirm the phylogenetic relationships in Eleotridae.     

Life on Arginine for Mycoplasma hominis: Clues from Its Minimal Genome and Comparison with Other Human Urogenital Mycoplasmas

Mycoplasma hominis is an opportunistic human mycoplasma. Two other pathogenic human species, M. genitalium and Ureaplasma parvum, reside within the same natural niche as M. hominis: the urogenital tract. These three species have overlapping, but distinct, pathogenic roles. They have minimal genomes and, thus, reduced metabolic capabilities characterized by distinct energy-generating pathways. Analysis of the M. hominis PG21 genome sequence revealed that it is the second smallest genome among self-replicating free living organisms (665,445 bp, 537 coding sequences (CDSs)). Five clusters of genes were predicted to have undergone horizontal gene transfer (HGT) between M. hominis and the phylogenetically distant U. parvum species. We reconstructed M. hominis metabolic pathways from the predicted genes, with particular emphasis on energy-generating pathways. The Embden–Meyerhoff–Parnas pathway was incomplete, with a single enzyme absent. We identified the three proteins constituting the arginine dihydrolase pathway. This pathway was found essential to promote growth in vivo. The predicted presence of dimethylarginine dimethylaminohydrolase suggested that arginine catabolism is more complex than initially described. This enzyme may have been acquired by HGT from non-mollicute bacteria. Comparison of the three minimal mollicute genomes showed that 247 CDSs were common to all three genomes, whereas 220 CDSs were specific to M. hominis, 172 CDSs were specific to M. genitalium, and 280 CDSs were specific to U. parvum. Within these species-specific genes, two major sets of genes could be identified: one including genes involved in various energy-generating pathways, depending on the energy source used (glucose, urea, or arginine) and another involved in cytadherence and virulence. Therefore, a minimal mycoplasma cell, not including cytadherence and virulence-related genes, could be envisaged containing a core genome (247 genes), plus a set of genes required for providing energy. For M. hominis, this set would include 247+9 genes, resulting in a theoretical minimal genome of 256 genes.