Pathogenomes of Shiga Toxin Positive and Negative Escherichia coli O157:H7 Strains TT12A and TT12B: Comprehensive Phylogenomic Analysis Using Closed Genomes.

Shiga toxin-producing Escherichia coli are zoonotic pathogens that cause food-borne human disease. Among these, the O157:H7 serotype has evolved from an enteropathogenic O55:H7 ancestor through the displacement of the somatic gene cluster and recurrent toxigenic conversion by Shiga toxin-converting bacteriophages. However, atypical strains that lack the Shiga toxin, the characteristic virulence hallmark, are circulating in this lineage. For this study, we analyzed the pathogenome and virulence inventories of the stx+ strain, TT12A, isolated from a patient with hemorrhagic colitis, and its respective co-isolated stx- strain, TT12B. Sequencing the genomes to closure proved critical to the cataloguing of subtle strain differentiating sequence and structural polymorphisms at a high-level of phylogenetic accuracy and resolution. Phylogenomic profiling revealed SNP and MLST profiles similar to the near clonal outbreak isolates. Their prophage inventories, however, were notably different. The attenuated atypical non-shigatoxigenic status of TT12B is explained by the absence of both the ΦStx1a- and ΦStx2a-prophages carried by TT12A, and we also recorded further alterations in the non-Stx prophage complement. Phenotypic characterization indicated that culture growth was directly impacted by the strains' distinct lytic phage complement. Altogether, our phylogenomic and phenotypic analyses show that these intimately related isogenic strains are on divergent Stx(+/stx-) evolutionary paths.

Pathogenomes and virulence profiles of representative big six non-O157 serogroup Shiga toxin-producing Escherichia coli.

Shiga toxin (Stx)-producing Escherichia coli (STEC) of non-O157:H7 serotypes are responsible for global and widespread human food-borne disease. Among these serogroups, O26, O45, O103, O111, O121, and O145 account for the majority of clinical infections and are colloquially referred to as the "Big Six." The "Big Six" strain panel we sequenced and analyzed in this study are reference type cultures comprised of six strains representing each of the non-O157 STEC serogroups curated and distributed by the American Type Culture Collection (ATCC) as a resource to the research community under panel number ATCC MP-9. The application of long- and short-read hybrid sequencing yielded closed chromosomes and a total of 14 plasmids of diverse functions. Through high-resolution comparative phylogenomics, we cataloged the shared and strain-specific virulence and resistance gene content and established the close relationship of serogroup O26 and O103 strains featuring flagellar H-type 11. Virulence phenotyping revealed statistically significant differences in the Stx-production capabilities that we found to be correlated to the strain's individual stx-status. Among the carried Stx1a, Stx2a, and Stx2d phages, the Stx2a phage is by far the most responsive upon RecA-mediated phage mobilization, and in consequence, stx2a + isolates produced the highest-level of toxin in this panel. The availability of high-quality closed genomes for this "Big Six" reference set, including carried plasmids, along with the recorded genomic virulence profiles and Stx-production phenotypes will provide a valuable foundation to further explore the plasticity in evolutionary trajectories in these emerging non-O157 STEC lineages, which are major culprits of human food-borne disease.

Pathogenome comparison and global phylogeny of Escherichia coli ST1485 strains.

Escherichia coli ST1485 strains belong to the clinically important phylogroup F and have disseminated worldwide in humans, animals, and the environment. Here, we elucidated the pathogenome of a global collection of E. coli ST1485 isolates from diverse sources retrieved from public databases and a high-quality sequenced complete genome of colistin-resistant E. coli strain CFSAN061771 isolated from raw milk cheese which designated as a reference strain. CFSAN061771 belongs to O83:H42-ST1485 pathotype and carries a conjugative ColV plasmid, pCFSAN061771_01, combining extraintestinal virulence genes (ompt, sitA, iroN, etsC, traT, cvaC, hylF, iss, tsh, mchf, iucC, iutA) with a multidrug resistance island (blaTEM-1, aph(6)-Id, aph(3″)-Ib, sul2, dfrA14). Comparative genomic analysis revealed a high frequency of pCFSAN061771_01-like plasmids in E. coli ST1485. A notable evolutionary genetic event in E. coli ST1485 strains is the acquisition of a pCFSAN061771_02-like plasmid, which confers resistance to several antimicrobials, tellurium, and quaternary ammonium compounds. The identical virulence and antibiotic resistance profiles identified in some human and animal strains are worrisome. This is the first study to emphasize the significance of E. coli ST1485 as a global high-risk virulent and multidrug-resistant clone with zoonotic potential.

Pathogenomes and variations in Shiga toxin production among geographically distinct clones of Escherichia coli O113:H21.

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains' featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.

Pathogenomes of Atypical Non-shigatoxigenic Escherichia coli NSF/SF O157:H7/NM: Comprehensive Phylogenomic Analysis Using Closed Genomes.

The toxigenic conversion of Escherichia coli strains by Shiga toxin-converting (Stx) bacteriophages were prominent and recurring events in the stepwise evolution of enterohemorrhagic E. coli (EHEC) O157:H7 from an enteropathogenic (EPEC) O55:H7 ancestor. Atypical, attenuated isolates have been described for both non-sorbitol fermenting (NSF) O157:H7 and SF O157:NM serotypes, which are distinguished by the absence of Stx, the characteristic virulence hallmark of Stx-producing E. coli (STEC). Such atypical isolates either never acquired Stx-phages or may have secondarily lost stx during the course of infection, isolation, or routine subculture; the latter are commonly referred to as LST (Lost Shiga Toxin)-isolates. In this study we analyzed the genomes of 15 NSF O157:H7 and SF O157:NM strains from North America, Europe, and Asia that are characterized by the absence of stx, the virulence hallmark of STEC. The individual genomic basis of the Stx (-) phenotype has remained largely undetermined as the majority of STEC genomes in public genome repositories were generated using short read technology and are in draft stage, posing a major obstacle for the high-resolution whole genome sequence typing (WGST). The application of LRT (long-read technology) sequencing provided us with closed genomes, which proved critical to put the atypical non-shigatoxigenic NSF O157:H7 and SF O157:NM strains into the phylogenomic context of the stepwise evolutionary model. Availability of closed chromosomes for representative Stx (-) NSF O157:H7 and SF O157:NM strains allowed to describe the genomic basis and individual evolutionary trajectories underlying the absence of Stx at high accuracy and resolution. The ability of LRT to recover and accurately assemble plasmids revealed a strong correlation between the strains' featured plasmid genotype and chromosomally inferred clade, which suggests the coevolution of the chromosome and accessory plasmids. The identified ancestral traits in the pSFO157 plasmid of NSF O157:H7 strain LSU-61 provided additional evidence for its intermediate status. Taken together, these observations highlight the utility of LRTs for advancing our understanding of EHEC O157:H7/NM pathogenome evolution. Insights into the genomic and phenotypic plasticity of STEC on a lineage- and genome-wide scale are foundational to improve and inform risk assessment, biosurveillance, and prevention strategies.

Closed Genome Sequence of Escherichia coli K-12 Group Strain C600.

Escherichia coli strain C600 is a prototypical K-12 derived laboratory strain which has been broadly used for molecular microbiology and bacterial physiology studies since its isolation in 1954. Here, we present the closed genome sequence of E. coli strain C600, retrieved from the American Type Culture Collection (ATCC 23724).

Closed Genome and Comparative Phylogenetic Analysis of the Clinical Multidrug Resistant Shigella sonnei Strain 866.

Shigella sonnei is responsible for the majority of shigellosis infections in the US with over 500,000 cases reported annually. Here, we present the complete genome of the clinical multidrug resistant (MDR) strain 866, which is highly susceptible to bacteriophage infections. The strain has a circular chromosome of 4.85 Mb and carries a 113 kb MDR plasmid. This IncB/O/K/Z-type plasmid, termed p866, confers resistance to five different classes of antibiotics including ß-lactamase, sulfonamide, tetracycline, aminoglycoside, and trimethoprim. Comparative analysis of the plasmid architecture and gene inventory revealed that p866 shares its plasmid backbone with previously described IncB/O/K/Z-type Shigella spp. and Escherichia coli plasmids, but is differentiated by the insertion of antibiotic resistance cassettes, which we found associated with mobile genetic elements such as Tn3, Tn7, and Tn10. A whole genome-derived phylogenetic reconstruction showed the evolutionary relationships of S. sonnei strain 866 and the four established Shigella species, highlighting the clonal nature of S. sonnei.

Closed Genome Sequence of Vibrio cholerae O1 El Tor Inaba Strain A1552.

Vibrio cholerae is a Gram-negative waterborne human pathogen and the causative agent of cholera. Here, we present the complete genome sequence of the seventh pandemic O1 biovar El Tor Inaba strain A1552 isolated in 1992. This clinical strain has served as an important model strain for studying cholera pathogenicity traits.

Whole Genome Sequencing for Genomics-Guided Investigations of Escherichia coli O157:H7 Outbreaks.

Multi isolate whole genome sequencing (WGS) and typing for outbreak investigations has become a reality in the post-genomics era. We applied this technology to strains from Escherichia coli O157:H7 outbreaks. These include isolates from seven North America outbreaks, as well as multiple isolates from the same patient and from different infected individuals in the same household. Customized high-resolution bioinformatics sequence typing strategies were developed to assess the core genome and mobilome plasticity. Sequence typing was performed using an in-house single nucleotide polymorphism (SNP) discovery and validation pipeline. Discriminatory power becomes of particular importance for the investigation of isolates from outbreaks in which macrogenomic techniques such as pulse-field gel electrophoresis or multiple locus variable number tandem repeat analysis do not differentiate closely related organisms. We also characterized differences in the phage inventory, allowing us to identify plasticity among outbreak strains that is not detectable at the core genome level. Our comprehensive analysis of the mobilome identified multiple plasmids that have not previously been associated with this lineage. Applied phylogenomics approaches provide strong molecular evidence for exceptionally little heterogeneity of strains within outbreaks and demonstrate the value of intra-cluster comparisons, rather than basing the analysis on archetypal reference strains. Next generation sequencing and whole genome typing strategies provide the technological foundation for genomic epidemiology outbreak investigation utilizing its significantly higher sample throughput, cost efficiency, and phylogenetic relatedness accuracy. These phylogenomics approaches have major public health relevance in translating information from the sequence-based survey to support timely and informed countermeasures. Polymorphisms identified in this work offer robust phylogenetic signals that index both short- and long-term evolution and can complement currently employed typing schemes for outbreak ex- and inclusion, diagnostics, surveillance, and forensic studies.

Genome Sequence of Bacillus thuringiensis Strain Btm27, an Egyptian Isolate Highly Toxic to Cotton Leafworm.

Bacillus thuringiensis is a potent microbial control agent against insect pests. Here, we present the draft genome of the Egyptian strain Btm27 that shows high toxicity toward the cotton leafworm. The genome contains three insecticidal genes cry1Ac9, cry2Ab1, and vip3V that have been implicated in conferring toxicity toward lepidoptera.

Characterization of the pathogenome and phylogenomic classification of enteropathogenic Escherichia coli of the O157:non-H7 serotypes.

Escherichia coli of the O157 serogroup are comprised of a diverse collection of more than 100 O157:non-H7 serotypes that are found in the environment, animal reservoir and infected patients and some have been linked to severe outbreaks of human disease. Among these, the enteropathogenic E. coli O157:non-H7 serotypes carry virulence factors that are hallmarks of enterohemorrhagic E. coli, such as causing attaching and effacing lesions during human gastrointestinal tract infections. Given the shared virulence gene pool between O157:H7 and O157:non-H7 serotypes, our objective was to examine the prevalence of virulence traits of O157:non-H7 serotypes within and across their H-serotype and when compared to other E. coli pathovars. We sequenced six O157:non-H7 genomes complemented by four genomes from public repositories in an effort to determine their virulence state and genetic relatedness to the highly pathogenic enterohemorrhagic O157:H7 lineage and its ancestral O55:H7 serotype. Whole-genome-based phylogenomic analysis and molecular typing is indicative of a non-monophyletic origin of the heterogeneous O157:non-H7 serotypes that are only distantly related to the O157:H7 serotype. The availability of multiple genomes enables robust phylogenomic placement of these strains into their evolutionary context, and the assessment of the pathogenic potential of the O157:non-H7 strains in causing human disease.

Genomic epidemiology of the Haitian cholera outbreak: a single introduction followed by rapid, extensive, and continued spread characterized the onset of the epidemic.

UNLABELLED: For centuries, cholera has been one of the most feared diseases. The causative agent Vibrio cholerae is a waterborne Gram-negative enteric pathogen eliciting a severe watery diarrheal disease. In October 2010, the seventh pandemic reached Haiti, a country that had not experienced cholera for more than a century. By using whole-genome sequence typing and mapping strategies of 116 serotype O1 strains from global sources, including 44 Haitian genomes, we present a detailed reconstructed evolutionary history of the seventh pandemic with a focus on the Haitian outbreak. We catalogued subtle genomic alterations at the nucleotide level in the genome core and architectural rearrangements from whole-genome map comparisons. Isolates closely related to the Haitian isolates caused several recent outbreaks in southern Asia. This study provides evidence for a single-source introduction of cholera from Nepal into Haiti followed by rapid, extensive, and continued clonal expansion. The phylogeographic patterns in both southern Asia and Haiti argue for the rapid dissemination of V. cholerae across the landscape necessitating real-time surveillance efforts to complement the whole-genome epidemiological analysis. As eradication efforts move forward, phylogeographic knowledge will be important for identifying persistent sources and monitoring success at regional levels. The results of molecular and epidemiological analyses of this outbreak suggest that an indigenous Haitian source of V. cholerae is unlikely and that an indigenous source has not contributed to the genomic evolution of this clade. IMPORTANCE: In this genomic epidemiology study, we have applied high-resolution whole-genome-based sequence typing methodologies on a comprehensive set of genome sequences that have become available in the aftermath of the Haitian cholera epidemic. These sequence resources enabled us to reassess the degree of genomic heterogeneity within the Vibrio cholerae O1 serotype and to refine boundaries and evolutionary relationships. The established phylogenomic framework showed how outbreak isolates fit into the global phylogeographic patterns compared to a comprehensive globally and temporally diverse strain collection and provides strong molecular evidence that points to a nonindigenous source of the 2010 Haitian cholera outbreak and refines epidemiological standards used in outbreak investigations for outbreak inclusion/exclusion following the concept of genomic epidemiology. The generated phylogenomic data have major public health relevance in translating sequence-based information to assist in future diagnostic, epidemiological, surveillance, and forensic studies of cholera.

Genome Sequence of Escherichia coli O157:H7 Strain 2886-75, Associated with the First Reported Case of Human Infection in the United States.

First identified in 1982 as a human pathogen, enterohemorrhagic Escherichia coli of the O157:H7 serotype is a major cause of food-borne acquired human infections. Here, we report the genome sequence of the first known strain of this serotype isolated in the United States.

Twice-daily application of HIV microbicides alter the vaginal microbiota.

UNLABELLED: Vaginal HIV microbicides offer great promise in preventing HIV transmission, but failures of phase 3 clinical trials, in which microbicide-treated subjects had an increased risk of HIV transmission, raised concerns about endpoints used to evaluate microbicide safety. A possible explanation for the increased transmission risk is that the agents shifted the vaginal bacterial community, resulting in loss of natural protection and enhanced HIV transmission susceptibility. We characterized vaginal microbiota, using pyrosequencing of bar-coded 16S rRNA gene fragments, in samples from 35 healthy, sexually abstinent female volunteer subjects (ages 18 to 50 years) with regular menses in a repeat phase 1 study of twice-daily application over 13.5 days of 1 of 3 gel products: a hydroxyethylcellulose (HEC)-based "universal" placebo (10 subjects), 6% cellulose sulfate (CS; 13 subjects), and 4% nonoxynol-9 (N-9; 12 subjects). We used mixed effects models inferred using Bayesian Markov chain Monte Carlo methods, which showed that treatment with active agents shifted the microbiota toward a community type lacking significant numbers of Lactobacillus spp. and dominated by strict anaerobes. This state of the vaginal microbiota was associated with a low or intermediate Nugent score and was not identical to bacterial vaginosis, an HIV transmission risk factor. The placebo arm contained a higher proportion of communities dominated by Lactobacillus spp., particularly L. crispatus, throughout treatment. The data suggest that molecular evaluation of microbicide effects on vaginal microbiota may be a critical endpoint that should be incorporated in early clinical assessment of microbicide candidates. IMPORTANCE: Despite large prevention efforts, HIV transmission and acquisition rates remain unacceptably high. In developing countries, transmission mainly occurs through heterosexual intercourse, where women are significantly more vulnerable to infection than men. Vaginal microbicides are considered to be one of the most promising female-controlled products, in that women themselves insert the microbicides into the vagina to prevent HIV transmission during sexual intercourse. The failure of several microbicides in clinical trials has raised questions concerning the low in vivo efficacy of such anti-HIV molecules. This study was designed to gain insights into the failures of two microbicides by testing the hypothesis that the microbicides negatively affect a critical line of defense against HIV, the vaginal microbiota. The results suggest that in the early assessment of candidate microbicides, culture-independent evaluation of their effect on the vaginal microbiota should be considered and may constitute a critical endpoint.

Temporal dynamics of the human vaginal microbiota.

Elucidating the factors that impinge on the stability of bacterial communities in the vagina may help in predicting the risk of diseases that affect women's health. Here, we describe the temporal dynamics of the composition of vaginal bacterial communities in 32 reproductive-age women over a 16-week period. The analysis revealed the dynamics of five major classes of bacterial communities and showed that some communities change markedly over short time periods, whereas others are relatively stable. Modeling community stability using new quantitative measures indicates that deviation from stability correlates with time in the menstrual cycle, bacterial community composition, and sexual activity. The women studied are healthy; thus, it appears that neither variation in community composition per se nor higher levels of observed diversity (co-dominance) are necessarily indicative of dysbiosis.

Proof of concept of microbiome-metabolome analysis and delayed gluten exposure on celiac disease autoimmunity in genetically at-risk infants.

Celiac disease (CD) is a unique autoimmune disorder in which the genetic factors (DQ2/DQ8) and the environmental trigger (gluten) are known and necessary but not sufficient for its development. Other environmental components contributing to CD are poorly understood. Studies suggest that aspects of gluten intake might influence the risk of CD occurrence and timing of its onset, i.e., the amount and quality of ingested gluten, together with the pattern of infant feeding and the age at which gluten is introduced in the diet. In this study, we hypothesize that the intestinal microbiota as a whole rather than specific infections dictates the switch from tolerance to immune response in genetically susceptible individuals. Using a sample of infants genetically at risk of CD, we characterized the longitudinal changes in the microbial communities that colonize infants from birth to 24 months and the impact of two patterns of gluten introduction (early vs. late) on the gut microbiota and metabolome, and the switch from gluten tolerance to immune response, including onset of CD autoimmunity. We show that infants genetically susceptible to CD who are exposed to gluten early mount an immune response against gluten and develop CD autoimmunity more frequently than at-risk infants in which gluten exposure is delayed until 12 months of age. The data, while derived from a relatively small number of subjects, suggest differences between the developing microbiota of infants with genetic predisposition for CD and the microbiota from infants with a non-selected genetic background, with an overall lack of bacteria of the phylum Bacteriodetes along with a high abundance of Firmicutes and microbiota that do not resemble that of adults even at 2 years of age. Furthermore, metabolomics analysis reveals potential biomarkers for the prediction of CD. This study constitutes a definite proof-of-principle that these combined genomic and metabolomic approaches will be key to deciphering the role of the gut microbiota on CD onset.

Genome sequences of the biotechnologically important Bacillus megaterium strains QM B1551 and DSM319.

Bacillus megaterium is deep-rooted in the Bacillus phylogeny, making it an evolutionarily key species and of particular importance in understanding genome evolution, dynamics, and plasticity in the bacilli. B. megaterium is a commercially available, nonpathogenic host for the biotechnological production of several substances, including vitamin B(12), penicillin acylase, and amylases. Here, we report the analysis of the first complete genome sequences of two important B. megaterium strains, the plasmidless strain DSM319 and QM B1551, which harbors seven indigenous plasmids. The 5.1-Mbp chromosome carries approximately 5,300 genes, while QM B1551 plasmids represent a combined 417 kb and 523 genes, one of the largest plasmid arrays sequenced in a single bacterial strain. We have documented extensive gene transfer between the plasmids and the chromosome. Each strain carries roughly 300 strain-specific chromosomal genes that account for differences in their experimentally confirmed phenotypes. B. megaterium is able to synthesize vitamin B(12) through an oxygen-independent adenosylcobalamin pathway, which together with other key energetic and metabolic pathways has now been fully reconstructed. Other novel genes include a second ftsZ gene, which may be responsible for the large cell size of members of this species, as well as genes for gas vesicles, a second ß-galactosidase gene, and most but not all of the genes needed for genetic competence. Comprehensive analyses of the global Bacillus gene pool showed that only an asymmetric region around the origin of replication was syntenic across the genus. This appears to be a characteristic feature of the Bacillus spp. genome architecture and may be key to their sporulating lifestyle.

Vaginal microbiome of reproductive-age women.

The means by which vaginal microbiomes help prevent urogenital diseases in women and maintain health are poorly understood. To gain insight into this, the vaginal bacterial communities of 396 asymptomatic North American women who represented four ethnic groups (white, black, Hispanic, and Asian) were sampled and the species composition characterized by pyrosequencing of barcoded 16S rRNA genes. The communities clustered into five groups: four were dominated by Lactobacillus iners, L. crispatus, L. gasseri, or L. jensenii, whereas the fifth had lower proportions of lactic acid bacteria and higher proportions of strictly anaerobic organisms, indicating that a potential key ecological function, the production of lactic acid, seems to be conserved in all communities. The proportions of each community group varied among the four ethnic groups, and these differences were statistically significant [χ(2)(10) = 36.8, P < 0.0001]. Moreover, the vaginal pH of women in different ethnic groups also differed and was higher in Hispanic (pH 5.0 ± 0.59) and black (pH 4.7 ± 1.04) women as compared with Asian (pH 4.4 ± 0.59) and white (pH 4.2 ± 0.3) women. Phylotypes with correlated relative abundances were found in all communities, and these patterns were associated with either high or low Nugent scores, which are used as a factor for the diagnosis of bacterial vaginosis. The inherent differences within and between women in different ethnic groups strongly argues for a more refined definition of the kinds of bacterial communities normally found in healthy women and the need to appreciate differences between individuals so they can be taken into account in risk assessment and disease diagnosis.

Comparison of self-collected and physician-collected vaginal swabs for microbiome analysis.

To our knowledge, no data are available on whether the microbial species composition and abundance sampled with self-collected vaginal swabs are comparable to those of swabs collected by clinicians. Twenty healthy women were recruited to the study during a routine gynecological visit. Eligible women were between 18 and 40 years old with regular menstrual cycles. Participants self-collected a vaginal swab using a standardized protocol and then were examined by a physician, who collected an additional five swabs from the lateral wall of the mid-vagina. In this study, the self-collected and three physician-obtained swabs were analyzed and compared using terminal restriction fragment length polymorphism and sequence analyses of the 16S rRNA genes. Vaginal microbial community comparative statistical analyses of both T-RFLP and 16S rRNA gene sequence datasets revealed that self-collected vaginal swabs sampled the same microbial diversity as physician collected swabs of the mid-vagina. These findings enable large-scale, field-based studies of the vaginal microbiome.