Bacterial community structure alterations within the colorectal cancer gut microbiome.

BACKGROUND: Colorectal cancer is a leading cause of cancer-related deaths worldwide. The human gut microbiome has become an active area of research for understanding the initiation, progression, and treatment of colorectal cancer. Despite multiple studies having found significant alterations in the carriage of specific bacteria within the gut microbiome of colorectal cancer patients, no single bacterium has been unequivocally connected to all cases. Whether alterations in species carriages are the cause or outcome of cancer formation is still unclear, but what is clear is that focus should be placed on understanding changes to the bacterial community structure within the cancer-associated gut microbiome. RESULTS: By applying a novel set of analyses on 252 previously published whole-genome shotgun sequenced fecal samples from healthy and late-stage colorectal cancer subjects, we identify taxonomic, functional, and structural changes within the cancer-associated human gut microbiome. Bacterial association networks constructed from these data exhibited widespread differences in the underlying bacterial community structure between healthy and colorectal cancer associated gut microbiomes. Within the cancer-associated ecosystem, bacterial species were found to form associations with other species that are taxonomically and functionally dissimilar to themselves, as well as form modules functionally geared towards potential changes in the tumor-associated ecosystem. Bacterial community profiling of these samples revealed a significant increase in species diversity within the cancer-associated gut microbiome, and an elevated relative abundance of species classified as originating from the oral microbiome including, but not limited to, Fusobacterium nucleatum, Peptostreptococcus stomatis, Gemella morbillorum, and Parvimonas micra. Differential abundance analyses of community functional capabilities revealed an elevation in functions linked to virulence factors and peptide degradation, and a reduction in functions involved in amino-acid biosynthesis within the colorectal cancer gut microbiome. CONCLUSIONS: We utilize whole-genome shotgun sequenced fecal samples provided from a large cohort of late-stage colorectal cancer and healthy subjects to identify a number of potentially important taxonomic, functional, and structural alterations occurring within the colorectal cancer associated gut microbiome. Our analyses indicate that the cancer-associated ecosystem influences bacterial partner selection in the native microbiota, and we highlight specific oral bacteria and their associations as potentially relevant towards aiding tumor progression.

Peroxidized Linoleic Acid, 13-HPODE, Alters Gene Expression Profile in Intestinal Epithelial Cells.

Lipid peroxides (LOOHs) abound in processed food and have been implicated in the pathology of diverse diseases including gut, cardiovascular, and cancer diseases. Recently, RNA Sequencing (RNA-seq) has been widely used to profile gene expression. To characterize gene expression and pathway dysregulation upon exposure to peroxidized linoleic acid, we incubated intestinal epithelial cells (Caco-2) with 100 µM of 13-hydroperoxyoctadecadienoic acid (13-HPODE) or linoleic acid (LA) for 24 h. Total RNA was extracted for library preparation and Illumina HiSeq sequencing. We identified 3094 differentially expressed genes (DEGs) in 13-HPODE-treated cells and 2862 DEGs in LA-treated cells relative to untreated cells. We show that 13-HPODE enhanced lipid metabolic pathways, including steroid hormone biosynthesis, PPAR signaling, and bile secretion, which alter lipid uptake and transport. 13-HPODE and LA treatments promoted detoxification mechanisms including cytochrome-P450. Conversely, both treatments suppressed oxidative phosphorylation. We also show that both treatments may promote absorptive cell differentiation and reduce proliferation by suppressing pathways involved in the cell cycle, DNA synthesis/repair and ribosomes, and enhancing focal adhesion. A qRT-PCR analysis of representative DEGs validated the RNA-seq analysis. This study provides insights into mechanisms by which 13-HPODE alters cellular processes and its possible involvement in mitochondrial dysfunction-related disorders and proposes potential therapeutic strategies to treat LOOH-related pathologies.

Identification of Nasal Gammaproteobacteria with Potent Activity against Staphylococcus aureus: Novel Insights into the "Noncarrier" State.

Staphylococcus aureus nasal carriage provides the bacterial reservoir for opportunistic infection. In comparing the nasal microbiomes of culture-defined persistent S. aureus carriers versus noncarriers, we detected S. aureus DNA in all noses, including those with an established history of S. aureus negativity based on culture. Colonization with Gammaproteobacteria, including Klebsiella aerogenes, Citrobacter koseri, Moraxella lincolnii, and select Acinetobacter spp., was associated with S. aureus noncarriage. We next developed physiological competition assays for testing anti-S. aureus activity of isolated nasal species, utilizing medium modeling the nutrient-limited fluid of the nasal mucosa, polarized primary nasal epithelia, and nasal secretions. K. aerogenes from the nose of an S. aureus noncarrier demonstrated >99% inhibition of S. aureus recovery in all assays, even when S. aureus was coincubated in 9-fold excess. Secreted S. aureus inhibitory proteins from K. aerogenes and M. lincolnii were heat-stable and <30 kDa, fitting the profile of antimicrobial peptides. C. koseri, Acinetobacter haemolyticus, Acinetobacter junii, and Acinetobacter schindleri inhibited S. aureus recovery on nasal epithelia in a contact-dependent manner, while several other species either had no effect or promoted S. aureus growth. Collectively, this project is one of the first to identify resident nasal microbial species that impede S. aureus survival, and it implies that detectable nasal S. aureus results from shifts in microbial community composition.IMPORTANCE Nasal carriage of Staphylococcus aureus is a risk factor for infection, but it is not yet understood why some individuals carry nasal S. aureus persistently, intermittently, or seemingly not at all when tested via culture methods. This study compared the nasal microbiomes of established S. aureus carriers and noncarriers, identified species associated with noncarriage, and tested them for anti-S. aureus activity using assays developed to model the nutrient-limited nasal mucosa. We determined that all nostril swabs contain S. aureus DNA, even swabs from hosts considered to be long-term noncarriers. Select members of the Gammaproteobacteria class were more prevalent in noncarrier than carrier nostrils and demonstrated potent activity against multiple strains of S. aureus The results described here provide a better understanding of how the nasal microbiome controls S. aureus growth and viability and may be useful in the design of improved S. aureus decolonization strategies.

Gut microbiota dependent anti-tumor immunity restricts melanoma growth in Rnf5-/- mice.

Accumulating evidence points to an important role for the gut microbiome in anti-tumor immunity. Here, we show that altered intestinal microbiota contributes to anti-tumor immunity, limiting tumor expansion. Mice lacking the ubiquitin ligase RNF5 exhibit attenuated activation of the unfolded protein response (UPR) components, which coincides with increased expression of inflammasome components, recruitment and activation of dendritic cells and reduced expression of antimicrobial peptides in intestinal epithelial cells. Reduced UPR expression is also seen in murine and human melanoma tumor specimens that responded to immune checkpoint therapy. Co-housing of Rnf5-/- and WT mice abolishes the anti-tumor immunity and tumor inhibition phenotype, whereas transfer of 11 bacterial strains, including B. rodentium, enriched in Rnf5-/- mice, establishes anti-tumor immunity and restricts melanoma growth in germ-free WT mice. Altered UPR signaling, exemplified in Rnf5-/- mice, coincides with altered gut microbiota composition and anti-tumor immunity to control melanoma growth.

A Microbiomic Analysis in African Americans with Colonic Lesions Reveals Streptococcus sp.VT162 as a Marker of Neoplastic Transformation.

Increasing evidence suggests a role of the gut microbiota in colorectal carcinogenesis (CRC). To detect bacterial markers of colorectal cancer in African Americans a metabolomic analysis was performed on fecal water extracts. DNA from stool samples of adenoma and healthy subjects and from colon cancer and matched normal tissues was analyzed to determine the microbiota composition (using 16S rDNA) and genomic content (metagenomics). Metagenomic functions with discriminative power between healthy and neoplastic specimens were established. Quantitative Polymerase Chain Reaction (q-PCR) using primers and probes specific to Streptococcus sp. VT_162 were used to validate this bacterium association with neoplastic transformation in stool samples from two independent cohorts of African Americans and Chinese patients with colorectal lesions. The metabolomic analysis of adenomas revealed low amino acids content. The microbiota in both cancer vs. normal tissues and adenoma vs. normal stool samples were different at the 16S rRNA gene level. Cross-mapping of metagenomic data led to 9 markers with significant discriminative power between normal and diseased specimens. These markers identified with Streptococcus sp. VT_162. Q-PCR data showed a statistically significant presence of this bacterium in advanced adenoma and cancer samples in an independent cohort of CRC patients. We defined metagenomic functions from Streptococcus sp. VT_162 with discriminative power among cancers vs. matched normal and adenomas vs. healthy subjects' stools. Streptococcus sp. VT_162 specific 16S rDNA was validated in an independent cohort. These findings might facilitate non-invasive screening for colorectal cancer.

Gut Microbiome-Based Metagenomic Signature for Non-invasive Detection of Advanced Fibrosis in Human Nonalcoholic Fatty Liver Disease.

The presence of advanced fibrosis in nonalcoholic fatty liver disease (NAFLD) is the most important predictor of liver mortality. There are limited data on the diagnostic accuracy of gut microbiota-derived signature for predicting the presence of advanced fibrosis. In this prospective study, we characterized the gut microbiome compositions using whole-genome shotgun sequencing of DNA extracted from stool samples. This study included 86 uniquely well-characterized patients with biopsy-proven NAFLD, of which 72 had mild/moderate (stage 0-2 fibrosis) NAFLD, and 14 had advanced fibrosis (stage 3 or 4 fibrosis). We identified a set of 40 features (p < 0.006), which included 37 bacterial species that were used to construct a Random Forest classifier model to distinguish mild/moderate NAFLD from advanced fibrosis. The model had a robust diagnostic accuracy (AUC 0.936) for detecting advanced fibrosis. This study provides preliminary evidence for a fecal-microbiome-derived metagenomic signature to detect advanced fibrosis in NAFLD.

Effect of Macondo Prospect 252 Oil on Microbiota Associated with Pelagic Sargassum in the Northern Gulf of Mexico.

The environmental impact of major oil spills on marine microorganisms has yet to be thoroughly investigated using molecular biology techniques. The Deepwater Horizon (DWH) drilling rig explosion of 2010 affected an approximately 176,000 km2 surface area of the Gulf of Mexico (GOM) when an estimated 210 million gallons of oil from the Macondo Prospect spilled into the environment. Pelagic Sargassum, a complex of two surface drifting species (Sargassum natans and Sargassum fluitans) of marine brown macroalgae and a critically important habitat in the GOM ecosystem, was suffused by Macondo Prospect 252 oil released during the DWH event. Using 16S rRNA PCR and Roche 454 pyrosequencing, the effect of the oil on the bacterial population associated with pelagic Sargassum and contiguous waters was examined by comparing sequence data generated from samples collected from oiled and non-oiled locations in the northern GOM. Sequence data showed similar microbial composition in Sargassum regardless of exposure to oil primarily dominated by five phyla; Proteobacteria, Bacteroidetes, Actinobacteria, Verrucomicrobia, and unclassified bacteria. The microbial composition in water samples was significantly less diverse than for Sargassum and consisted primarily of Proteobacteria, Firmicutes, and Bacteroidetes. Due to the evenly distributed abundance of microbial species on oiled and non-oiled pelagic Sargassum, study findings indicate that DWH spilled oil had minimal effect on the composition and diversity of the microbial community associated with Sargassum and contiguous waters. However, higher abundances of Sulfitobacter and one species of Psychrobacter were found in oiled water samples when compared to non-oiled water samples indicating some effect of DHW oil in the microbial composition of seawater. Though there are a number of marine studies using molecular biology approaches, this is the first molecular examination of the impact of the DWH oil spill on bacterial communities associated with pelagic Sargassum and contiguous waters from the GOM.

GRASPx: efficient homolog-search of short peptide metagenome database through simultaneous alignment and assembly.

BACKGROUND: Metagenomics is a cultivation-independent approach that enables the study of the genomic composition of microbes present in an environment. Metagenomic samples are routinely sequenced using next-generation sequencing technologies that generate short nucleotide reads. Proteins identified from these reads are mostly of partial length. On the other hand, de novo assembly of a large metagenomic dataset is computationally demanding and the assembled contigs are often fragmented, resulting in the identification of protein sequences that are also of partial length and incomplete. Annotation of an incomplete protein sequence often proceeds by identifying its homologs in a database of reference sequences. Identifying the homologs of incomplete sequences is a challenge and can result in substandard annotation of proteins from metagenomic datasets. To address this problem, we recently developed a homology detection algorithm named GRASP (Guided Reference-based Assembly of Short Peptides) that identifies the homologs of a given reference protein sequence in a database of short peptide metagenomic sequences. GRASP was developed to implement a simultaneous alignment and assembly algorithm for annotation of short peptides identified on metagenomic reads. The program achieves significantly improved recall rate at the cost of computational efficiency. In this article, we adopted three techniques to speed up the original version of GRASP, including the pre-construction of extension links, local assembly of individual seeds, and the implementation of query-level parallelism. RESULTS: The resulting new program, GRASPx, achieves >30X speedup compared to its predecessor GRASP. At the same time, we show that the performance of GRASPx is consistent with that of GRASP, and that both of them significantly outperform other popular homology-search tools including the BLAST and FASTA suites. GRASPx was also applied to a human saliva metagenome dataset and shows superior performance for both recall and precision rates. CONCLUSIONS: In this article we present GRASPx, a fast and accurate homology-search program implementing a simultaneous alignment and assembly framework. GRASPx can be used for more comprehensive and accurate annotation of short peptides. GRASPx is freely available at http://graspx.sourceforge.net/ .

Gastrointestinal microbial populations can distinguish pediatric and adolescent Acute Lymphoblastic Leukemia (ALL) at the time of disease diagnosis.

BACKGROUND: An estimated 15,000 children and adolescents under the age of 19 years are diagnosed with leukemia, lymphoma and other tumors in the USA every year. All children and adolescent acute leukemia patients will undergo chemotherapy as part of their treatment regimen. Fortunately, survival rates for most pediatric cancers have improved at a remarkable pace over the past three decades, and the overall survival rate is greater than 90 % today. However, significant differences in survival rate have been found in different age groups (94 % in 1-9.99 years, 82 % in ≥10 years and 76 % in ≥15 years). ALL accounts for about three out of four cases of childhood leukemia. Intensive chemotherapy treatment coupled with prophylactic or therapeutic antibiotic use could potentially have a long-term effect on the resident gastrointestinal (GI) microbiome. The composition of GI microbiome and its changes upon chemotherapy in pediatric and adolescent leukemia patients is poorly understood. In this study, using 16S rRNA marker gene sequences we profile the GI microbial communities of pediatric and adolescent acute leukemia patients before and after chemotherapy treatment and compare with the microbiota of their healthy siblings. RESULTS: Our study cohort consisted of 51 participants, made up of matched pediatric and adolescent patients with ALL and a healthy sibling. We elucidated and compared the GI microbiota profiles of patients and their healthy sibling controls via analysis of 16S rRNA gene sequencing data. We assessed the GI microbiota composition in pediatric and adolescent patients with ALL during the course of chemotherapy by comparing stool samples taken before chemotherapy with stool samples collected at varying time points during the chemotherapeutic treatment. The microbiota profiles of both patients and control sibling groups are dominated by members of Bacteroides, Prevotella, and Faecalibacterium. At the genus level, both groups share many taxa in common, but the microbiota diversity of the patient group is significantly lower than that of the control group. It was possible to distinguish between the patient and control groups based on their microbiota profiles. The top taxa include Anaerostipes, Coprococcus, Roseburia, and Ruminococcus2 with relatively higher abundance in the control group. The observed microbiota changes are likely the result of several factors including a direct influence of therapeutic compounds on the gut flora and an indirect effect of chemotherapy on the immune system, which, in turn, affects the microbiota. CONCLUSIONS: This study provides significant information on GI microbiota populations in immunocompromised children and opens up the potential for developing novel diagnostics based on stool tests and therapies to improve the dysbiotic condition of the microbiota at the time of diagnosis and in the earliest stages of chemotherapy.

SFA-SPA: a suffix array based short peptide assembler for metagenomic data.

UNLABELLED: The determination of protein sequences from a metagenomic dataset enables the study of metabolism and functional roles of the organisms that are present in the sampled microbial community. We had previously introduced algorithm and software for the accurate reconstruction of protein sequences from short peptides identified on nucleotide reads in a metagenomic dataset. Here, we present significant computational improvements to the short peptide assembly algorithm that make it practical to reconstruct proteins from large metagenomic datasets containing several hundred million reads, while maintaining accuracy. The improved computational efficiency is achieved using a suffix array data structure that allows for fast querying during the assembly process, and a significant redesign of assembly steps that enables multi-threaded execution. AVAILABILITY AND IMPLEMENTATION: The program is available under the GPLv3 license from sourceforge.net/projects/spa-assembler.

GRASP: guided reference-based assembly of short peptides.

Protein sequences predicted from metagenomic datasets are annotated by identifying their homologs via sequence comparisons with reference or curated proteins. However, a majority of metagenomic protein sequences are partial-length, arising as a result of identifying genes on sequencing reads or on assembled nucleotide contigs, which themselves are often very fragmented. The fragmented nature of metagenomic protein predictions adversely impacts homology detection and, therefore, the quality of the overall annotation of the dataset. Here we present a novel algorithm called GRASP that accurately identifies the homologs of a given reference protein sequence from a database consisting of partial-length metagenomic proteins. Our homology detection strategy is guided by the reference sequence, and involves the simultaneous search and assembly of overlapping database sequences. GRASP was compared to three commonly used protein sequence search programs (BLASTP, PSI-BLAST and FASTM). Our evaluations using several simulated and real datasets show that GRASP has a significantly higher sensitivity than these programs while maintaining a very high specificity. GRASP can be a very useful program for detecting and quantifying taxonomic and protein family abundances in metagenomic datasets. GRASP is implemented in GNU C++, and is freely available at http://sourceforge.net/projects/grasp-release.

Cultivation of a human-associated TM7 phylotype reveals a reduced genome and epibiotic parasitic lifestyle.

The candidate phylum TM7 is globally distributed and often associated with human inflammatory mucosal diseases. Despite its prevalence, the TM7 phylum remains recalcitrant to cultivation, making it one of the most enigmatic phyla known. In this study, we cultivated a TM7 phylotype (TM7x) from the human oral cavity. This extremely small coccus (200-300 nm) has a distinctive lifestyle not previously observed in human-associated microbes. It is an obligate epibiont of an Actinomyces odontolyticus strain (XH001) yet also has a parasitic phase, thereby killing its host. This first completed genome (705 kb) for a human-associated TM7 phylotype revealed a complete lack of amino acid biosynthetic capacity. Comparative genomics analyses with uncultivated environmental TM7 assemblies show remarkable conserved gene synteny and only minimal gene loss/gain that may have occurred as TM7x adapted to conditions within the human host. Transcriptomic and metabolomic profiles provided the first indications, to our knowledge, that there is signaling interaction between TM7x and XH001. Furthermore, the induction of TNF-α production in macrophages by XH001 was repressed in the presence of TM7x, suggesting its potential immune suppression ability. Overall, our data provide intriguing insights into the uncultivability, pathogenicity, and unique lifestyle of this previously uncharacterized oral TM7 phylotype.

Human papillomavirus community in healthy persons, defined by metagenomics analysis of human microbiome project shotgun sequencing data sets.

UNLABELLED: Human papillomavirus (HPV) causes a number of neoplastic diseases in humans. Here, we show a complex normal HPV community in a cohort of 103 healthy human subjects, by metagenomics analysis of the shotgun sequencing data generated from the NIH Human Microbiome Project. The overall HPV prevalence was 68.9% and was highest in the skin (61.3%), followed by the vagina (41.5%), mouth (30%), and gut (17.3%). Of the 109 HPV types as well as additional unclassified types detected, most were undetectable by the widely used commercial kits targeting the vaginal/cervical HPV types. These HPVs likely represent true HPV infections rather than transitory exposure because of strong organ tropism and persistence of the same HPV types in repeat samples. Coexistence of multiple HPV types was found in 48.1% of the HPV-positive samples. Networking between HPV types, cooccurrence or exclusion, was detected in vaginal and skin samples. Large contigs assembled from short HPV reads were obtained from several samples, confirming their genuine HPV origin. This first large-scale survey of HPV using a shotgun sequencing approach yielded a comprehensive map of HPV infections among different body sites of healthy human subjects. IMPORTANCE: This nonbiased survey indicates that the HPV community in healthy humans is much more complex than previously defined by widely used kits that are target selective for only a few high- and low-risk HPV types for cervical cancer. The importance of nononcogenic viruses in a mixed HPV infection could be for stimulating or inhibiting a coexisting oncogenic virus via viral interference or immune cross-reaction. Knowledge gained from this study will be helpful to guide the designing of epidemiological and clinical studies in the future to determine the impact of nononcogenic HPV types on the outcome of HPV infections.

Microbiome analysis of stool samples from African Americans with colon polyps.

BACKGROUND: Colonic polyps are common tumors occurring in ~50% of Western populations with ~10% risk of malignant progression. Dietary agents have been considered the primary environmental exposure to promote colorectal cancer (CRC) development. However, the colonic mucosa is permanently in contact with the microbiota and its metabolic products including toxins that also have the potential to trigger oncogenic transformation. AIM: To analyze fecal DNA for microbiota composition and functional potential in African Americans with pre-neoplastic lesions. MATERIALS & METHODS: We analyzed the bacterial composition of stool samples from 6 healthy individuals and 6 patients with colon polyps using 16S ribosomal RNA-based phylogenetic microarray; the Human intestinal Tract Chip (HITChip) and 16S rRNA gene barcoded 454 pyrosequencing. The functional potential was determined by sequence-based metagenomics using 454 pyrosequencing. RESULTS: Fecal microbiota profiling of samples from the healthy and polyp patients using both a phylogenetic microarraying (HITChip) and barcoded 454 pyrosequencing generated similar results. A distinction between both sets of samples was only obtained when the analysis was performed at the sub-genus level. Most of the species leading to the dissociation were from the Bacteroides group. The metagenomic analysis did not reveal major differences in bacterial gene prevalence/abundances between the two groups even when the analysis and comparisons were restricted to available Bacteroides genomes. CONCLUSION: This study reveals that at the pre-neoplastic stages, there is a trend showing microbiota changes between healthy and colon polyp patients at the sub-genus level. These differences were not reflected at the genome/functions levels. Bacteria and associated functions within the Bacteroides group need to be further analyzed and dissected to pinpoint potential actors in the early colon oncogenic transformation in a large sample size.

SPA: a short peptide assembler for metagenomic data.

The metagenomic paradigm allows for an understanding of the metabolic and functional potential of microbes in a community via a study of their proteins. The substrate for protein identification is either the set of individual nucleotide reads generated from metagenomic samples or the set of contig sequences produced by assembling these reads. However, a read-based strategy using reads generated by next-generation sequencing (NGS) technologies, results in an overwhelming majority of partial-length protein predictions. A nucleotide assembly-based strategy does not fare much better, as metagenomic assemblies are typically fragmented and also leave a large fraction of reads unassembled. Here, we present a method for reconstructing complete protein sequences directly from NGS metagenomic data. Our framework is based on a novel short peptide assembler (SPA) that assembles protein sequences from their constituent peptide fragments identified on short reads. The SPA algorithm is based on informed traversals of a de Bruijn graph, defined on an amino acid alphabet, to identify probable paths that correspond to proteins. Using large simulated and real metagenomic data sets, we show that our method outperforms the alternate approach of identifying genes on nucleotide sequence assemblies and generates longer protein sequences that can be more effectively analysed.

Metagenomic investigation of microbes and viruses in patients with jaw osteonecrosis associated with bisphosphonate therapy.

OBJECTIVE: The goal of this preliminary study was to use metagenomic approaches to investigate the taxonomic diversity of microorganisms in patients with bisphosphonate-related osteonecrosis of the jaw (BRONJ). STUDY DESIGN: Samples of saliva for planktonic microbial analysis and biofilm cultivation were collected from 10 patients (5 with BRONJ and 5 non-BRONJ control subjects) who met all ascertainment criteria. Prophage induction experiments-16S rRNA polymerase chain reaction and 454 pyrosequencing-and epifluorescent microscopy were performed for characterization and enumeration of microbes and viruses. RESULTS: Three phyla of microbes-Proteobacteria (70%), Firmicutes (26.9%), and Actinobacteria (1.95%)-dominated all BRONJ samples and accounted for almost 99% of the total data. Viral abundance was ∼1 order of magnitude greater than microbial cell abundance and comprised mainly phage viruses. CONCLUSIONS: Individuals with jaw osteonecrosis harbored different microbial assemblages than nonaffected patients, and in general viral abundance and prophage induction increased with biofilm formation, suggesting that biofilm formation encouraged lysogenic interactions between viruses and microbial hosts and may contribute to pathogenicity.

Comparative genome analysis of 19 Ureaplasma urealyticum and Ureaplasma parvum strains.

BACKGROUND: Ureaplasma urealyticum (UUR) and Ureaplasma parvum (UPA) are sexually transmitted bacteria among humans implicated in a variety of disease states including but not limited to: nongonococcal urethritis, infertility, adverse pregnancy outcomes, chorioamnionitis, and bronchopulmonary dysplasia in neonates. There are 10 distinct serotypes of UUR and 4 of UPA. Efforts to determine whether difference in pathogenic potential exists at the ureaplasma serovar level have been hampered by limitations of antibody-based typing methods, multiple cross-reactions and poor discriminating capacity in clinical samples containing two or more serovars. RESULTS: We determined the genome sequences of the American Type Culture Collection (ATCC) type strains of all UUR and UPA serovars as well as four clinical isolates of UUR for which we were not able to determine serovar designation. UPA serovars had 0.75-0.78 Mbp genomes and UUR serovars were 0.84-0.95 Mbp. The original classification of ureaplasma isolates into distinct serovars was largely based on differences in the major ureaplasma surface antigen called the multiple banded antigen (MBA) and reactions of human and animal sera to the organisms. Whole genome analysis of the 14 serovars and the 4 clinical isolates showed the mba gene was part of a large superfamily, which is a phase variable gene system, and that some serovars have identical sets of mba genes. Most of the differences among serovars are hypothetical genes, and in general the two species and 14 serovars are extremely similar at the genome level. CONCLUSIONS: Comparative genome analysis suggests UUR is more capable of acquiring genes horizontally, which may contribute to its greater virulence for some conditions. The overwhelming evidence of extensive horizontal gene transfer among these organisms from our previous studies combined with our comparative analysis indicates that ureaplasmas exist as quasi-species rather than as stable serovars in their native environment. Therefore, differential pathogenicity and clinical outcome of a ureaplasmal infection is most likely not on the serovar level, but rather may be due to the presence or absence of potential pathogenicity factors in an individual ureaplasma clinical isolate and/or patient to patient differences in terms of autoimmunity and microbiome.