Sexual transmission of urogenital bacteria: whole metagenome sequencing evidence from a sexual network study.

Sexual transmission of the urogenital microbiota may contribute to adverse sexual and reproductive health outcomes. The extent of sexual transmission of the urogenital microbiota is unclear as prior studies largely investigated specific pathogens. We used epidemiologic data and whole metagenome sequencing to characterize urogenital microbiota strain concordance between participants of a sexual network study. Individuals who screened positive for genital Chlamydia trachomatis were enrolled and referred their sexual contacts from the prior 60-180 days. Snowball recruitment of sexual contacts continued for up to four waves. Vaginal swabs and penile urethral swabs were collected for whole metagenome sequencing. We evaluated bacterial strain concordance using inStrain and network analysis. We defined concordance as ≥99.99% average nucleotide identity over ≥50% shared coverage; we defined putative sexual transmission as concordance between sexual contacts with <5 single-nucleotide polymorphisms per megabase. Of 138 participants, 74 (54%) were female; 120 (87%) had genital chlamydia; and 43 (31%) were recruited contacts. We identified 115 strain-concordance events among 54 participants representing 25 bacterial species. Seven events (6%) were between sexual contacts including putative heterosexual transmission of Fannyhessea vaginae, Gardnerella leopoldii, Prevotella amnii, Sneathia sanguinegens, and Sneathia vaginalis (one strain each), and putative sexual transmission of Lactobacillus iners between female contacts. Most concordance events (108, 94%) were between non-contacts, including eight female participants connected through 18 Lactobacillus crispatus and 3 Lactobacillus jensenii concordant strains, and 14 female and 2 male participants densely interconnected through 52 Gardnerella swidsinskii concordance events.IMPORTANCEEpidemiologic evidence consistently indicates bacterial vaginosis (BV) is sexually associated and may be sexually transmitted, though sexual transmission remains subject to debate. This study is not capable of demonstrating BV sexual transmission; however, we do provide strain-level metagenomic evidence that strongly supports heterosexual transmission of BV-associated species. These findings strengthen the evidence base that supports ongoing investigations of concurrent male partner treatment for reducing BV recurrence. Our data suggest that measuring the impact of male partner treatment on F. vaginae, G. leopoldii, P. amnii, S. sanguinegens, and S. vaginalis may provide insight into why a regimen does or does not perform well. We also observed a high degree of strain concordance between non-sexual-contact female participants. We posit that this may reflect limited dispersal capacity of vaginal bacteria coupled with individuals' comembership in regional transmission networks where transmission may occur between parent and child at birth, cohabiting individuals, or sexual contacts.

Integrating compositional and functional content to describe vaginal microbiomes in health and disease.

BACKGROUND: A Lactobacillus-dominated vaginal microbiome provides the first line of defense against adverse genital tract health outcomes. However, there is limited understanding of the mechanisms by which the vaginal microbiome modulates protection, as prior work mostly described its composition through morphologic assessment and marker gene sequencing methods that do not capture functional information. To address this gap, we developed metagenomic community state types (mgCSTs) which use metagenomic sequences to describe and define vaginal microbiomes based on both composition and functional potential. RESULTS: MgCSTs are categories of microbiomes classified using taxonomy and the functional potential encoded in their metagenomes. MgCSTs reflect unique combinations of metagenomic subspecies (mgSs), which are assemblages of bacterial strains of the same species, within a microbiome. We demonstrate that mgCSTs are associated with demographics such as age and race, as well as vaginal pH and Gram stain assessment of vaginal smears. Importantly, these associations varied between mgCSTs predominated by the same bacterial species. A subset of mgCSTs, including three of the six predominated by Gardnerella vaginalis mgSs, as well as mgSs of L. iners, were associated with a greater likelihood of bacterial vaginosis diagnosed by Amsel clinical criteria. This L. iners mgSs, among other functional features, encoded enhanced genetic capabilities for epithelial cell attachment that could facilitate cytotoxin-mediated cell lysis. Finally, we report a mgSs and mgCST classifier for which source code is provided and may be adapted for use by the microbiome research community. CONCLUSIONS: MgCSTs are a novel and easily implemented approach to reduce the dimension of complex metagenomic datasets while maintaining their functional uniqueness. MgCSTs enable the investigation of multiple strains of the same species and the functional diversity in that species. Future investigations of functional diversity may be key to unraveling the pathways by which the vaginal microbiome modulates the protection of the genital tract. Importantly, our findings support the hypothesis that functional differences between vaginal microbiomes, including those that may look compositionally similar, are critical considerations in vaginal health. Ultimately, mgCSTs may lead to novel hypotheses concerning the role of the vaginal microbiome in promoting health and disease, and identify targets for novel prognostic, diagnostic, and therapeutic strategies to improve women's genital health. Video Abstract.

Evaluation of vaginal microbiome equilibrium states identifies microbial parameters linked to resilience after menses and antibiotic therapy.

The vaginal microbiome (VMB) is a complex microbial community that is closely tied to reproductive health. Optimal VMB communities have compositions that are commonly defined by the dominance of certain Lactobacillus spp. and can remain stable over time or transition to non-optimal states dominated by anaerobic bacteria and associated with bacterial vaginosis (BV). The ability to remain stable or undergo transitions suggests a system with either single (mono-stable) or multiple (multi-stable) equilibrium states, though factors that contribute to stability have been difficult to determine due to heterogeneity in microbial growth characteristics and inter-species interactions. Here, we use a computational model to determine whether differences in microbial growth and interaction parameters could alter equilibrium state accessibility and account for variability in community composition after menses and antibiotic therapies. Using a global uncertainty and sensitivity analysis that captures parameter sets sampled from a physiologically relevant range, model simulations predicted that 79.7% of microbial communities were mono-stable (gravitate to one composition type) and 20.3% were predicted to be multi-stable (can gravitate to more than one composition type, given external perturbations), which was not significantly different from observations in two clinical cohorts (HMP cohort, 75.2% and 24.8%; Gajer cohort, 78.1% and 21.9%, respectively). The model identified key microbial parameters that governed equilibrium state accessibility, such as the importance of non-optimal anaerobic bacteria interactions with Lactobacillus spp., which is largely understudied. Model predictions for composition changes after menses and antibiotics were not significantly different from those observed in clinical cohorts. Lastly, simulations were performed to illustrate how this quantitative framework can be used to gain insight into the development of new combinatorial therapies involving altered prebiotic and antibiotic dosing strategies. Altogether, dynamical models could guide development of more precise therapeutic strategies to manage BV.

Bacterial amylases enable glycogen degradation by the vaginal microbiome.

The human vaginal microbiota is frequently dominated by lactobacilli and transition to a more diverse community of anaerobic microbes is associated with health risks. Glycogen released by lysed epithelial cells is believed to be an important nutrient source in the vagina. However, the mechanism by which vaginal bacteria metabolize glycogen is unclear, with evidence implicating both bacterial and human enzymes. Here we biochemically characterize six glycogen-degrading enzymes (GDEs), all of which are pullanases (PulA homologues), from vaginal bacteria that support the growth of amylase-deficient Lactobacillus crispatus on glycogen. We reveal variations in their pH tolerance, substrate preferences, breakdown products and susceptibility to inhibition. Analysis of vaginal microbiome datasets shows that these enzymes are expressed in all community state types. Finally, we confirm the presence and activity of bacterial and human GDEs in cervicovaginal fluid. This work establishes that bacterial GDEs can participate in the breakdown of glycogen, providing insight into metabolism that may shape the vaginal microbiota.

High-resolution functional description of vaginal microbiomes in health and disease.

Background: A Lactobacillus-dominated vaginal microbiome provides the first line of defense against numerous adverse genital tract health outcomes. However, there is limited understanding of the mechanisms by which the vaginal microbiome modulates protection, as prior work mostly described its composition through morphologic assessment and marker gene sequencing methods that do not capture functional information. To address this limitation, we developed metagenomic community state types (mgCSTs) which uses metagenomic sequences to describe and define vaginal microbiomes based on both composition and function. Results: MgCSTs are categories of microbiomes classified using taxonomy and the functional potential encoded in their metagenomes. MgCSTs reflect unique combinations of metagenomic subspecies (mgSs), which are assemblages of bacterial strains of the same species, within a microbiome. We demonstrate that mgCSTs are associated with demographics such as age and race, as well as vaginal pH and Gram stain assessment of vaginal smears. Importantly, these associations varied between mgCSTs predominated by the same bacterial species. A subset of mgCSTs, including three of the six predominated by Gardnerella mgSs, as well as a mgSs of L. iners, were associated with a greater likelihood of Amsel bacterial vaginosis diagnosis. This L. iners mgSs, among other functional features, encoded enhanced genetic capabilities for epithelial cell attachment that could facilitate cytotoxin-mediated cell lysis. Finally, we report a mgSs and mgCST classifier as an easily applied, standardized method for use by the microbiome research community. Conclusions: MgCSTs are a novel and easily implemented approach to reducing the dimension of complex metagenomic datasets, while maintaining their functional uniqueness. MgCSTs enable investigation of multiple strains of the same species and the functional diversity in that species. Future investigations of functional diversity may be key to unraveling the pathways by which the vaginal microbiome modulates protection to the genital tract. Importantly, our findings support the hypothesis that functional differences between vaginal microbiomes, including those that may look compositionally similar, are critical considerations in vaginal health. Ultimately, mgCSTs may lead to novel hypotheses concerning the role of the vaginal microbiome in promoting health and disease, and identify targets for novel prognostic, diagnostic, and therapeutic strategies to improve women's genital health.

Identification of shared bacterial strains in the vaginal microbiota of related and unrelated reproductive-age mothers and daughters using genome-resolved metagenomics.

It has been suggested that the human microbiome might be vertically transmitted from mother to offspring and that early colonizers may play a critical role in development of the immune system. Studies have shown limited support for the vertical transmission of the intestinal microbiota but the derivation of the vaginal microbiota remains largely unknown. Although the vaginal microbiota of children and reproductive age women differ in composition, the vaginal microbiota could be vertically transmitted. To determine whether there was any support for this hypothesis, we examined the vaginal microbiota of daughter-mother pairs from the Baltimore metropolitan area (ages 14-27, 32-51; n = 39). We assessed whether the daughter's microbiota was similar in composition to their mother's using metataxonomics. Permutation tests revealed that while some pairs did have similar vaginal microbiota, the degree of similarity did not exceed that expected by chance. Genome-resolved metagenomics was used to identify shared bacterial strains in a subset of the families (n = 22). We found a small number of bacterial strains that were shared between mother-daughter pairs but identified more shared strains between individuals from different families, indicating that vaginal bacteria may display biogeographic patterns. Earlier-in-life studies are needed to demonstrate vertical transmission of the vaginal microbiota.

Complete Genome Sequences of Ezakiella coagulans C0061C1 and Fenollaria massiliensis C0061C2.

Ezakiella coagulans and Fenollaria massiliensis are two obligate anaerobic bacteria in the family Peptoniphilaceae and are both uncommon members of the human vaginal microbiota. We isolated a strain of each bacterium from the same vaginal swab specimen and here report the first complete genome sequences of the two species.

Insight into the ecology of vaginal bacteria through integrative analyses of metagenomic and metatranscriptomic data.

BACKGROUND: Vaginal bacterial communities dominated by Lactobacillus species are associated with a reduced risk of various adverse health outcomes. However, somewhat unexpectedly, many healthy women have microbiota that are not dominated by lactobacilli. To determine the factors that drive vaginal community composition we characterized the genetic composition and transcriptional activities of vaginal microbiota in healthy women. RESULTS: We demonstrate that the abundance of a species is not always indicative of its transcriptional activity and that impending changes in community composition can be predicted from metatranscriptomic data. Functional comparisons highlight differences in the metabolic activities of these communities, notably in their degradation of host produced mucin but not glycogen. Degradation of mucin by communities not dominated by Lactobacillus may play a role in their association with adverse health outcomes. Finally, we show that the transcriptional activities of L. crispatus, L. iners, and Gardnerella vaginalis vary with the taxonomic composition of the communities in which they reside. Notably, L. iners and G. vaginalis both demonstrate lower expression of their cholesterol-dependent cytolysins when co-resident with Lactobacillus spp. and higher expression when co-resident with other facultative and obligate anaerobes. The pathogenic potential of these species may depend on the communities in which they reside and thus could be modulated by interventional strategies. CONCLUSIONS: Our results provide insight to the functional ecology of the vaginal microbiota, demonstrate the diagnostic potential of metatranscriptomic data, and reveal strategies for the management of these ecosystems.

Quantitative modeling predicts mechanistic links between pre-treatment microbiome composition and metronidazole efficacy in bacterial vaginosis.

Bacterial vaginosis is a condition associated with adverse reproductive outcomes and characterized by a shift from a Lactobacillus-dominant vaginal microbiota to a polymicrobial microbiota, consistently colonized by strains of Gardnerella vaginalis. Metronidazole is the first-line treatment; however, treatment failure and recurrence rates remain high. To understand complex interactions between Gardnerella vaginalis and Lactobacillus involved in efficacy, here we develop an ordinary differential equation model that predicts bacterial growth as a function of metronidazole uptake, sensitivity, and metabolism. The model shows that a critical factor in efficacy is Lactobacillus sequestration of metronidazole, and efficacy decreases when the relative abundance of Lactobacillus is higher pre-treatment. We validate results in Gardnerella and Lactobacillus co-cultures, and in two clinical cohorts, finding women with recurrence have significantly higher pre-treatment levels of Lactobacillus relative to bacterial vaginosis-associated bacteria. Overall results provide mechanistic insight into how personalized differences in microbial communities influence vaginal antibiotic efficacy.

VALENCIA: a nearest centroid classification method for vaginal microbial communities based on composition.

BACKGROUND: Taxonomic profiles of vaginal microbial communities can be sorted into a discrete number of categories termed community state types (CSTs). This approach is advantageous because collapsing a hyper-dimensional taxonomic profile into a single categorical variable enables efforts such as data exploration, epidemiological studies, and statistical modeling. Vaginal communities are typically assigned to CSTs based on the results of hierarchical clustering of the pairwise distances between samples. However, this approach is problematic because it complicates between-study comparisons and because the results are entirely dependent on the particular set of samples that were analyzed. We sought to standardize and advance the assignment of samples to CSTs. RESULTS: We developed VALENCIA (VAginaL community state typE Nearest CentroId clAssifier), a nearest centroid-based tool which classifies samples based on their similarity to a set of reference centroids. The references were defined using a comprehensive set of 13,160 taxonomic profiles from 1975 women in the USA. This large dataset allowed us to comprehensively identify, define, and characterize vaginal CSTs common to reproductive age women and expand upon the CSTs that had been defined in previous studies. We validated the broad applicability of VALENCIA for the classification of vaginal microbial communities by using it to classify three test datasets which included reproductive age eastern and southern African women, adolescent girls, and a racially/ethnically and geographically diverse sample of postmenopausal women. VALENCIA performed well on all three datasets despite the substantial variations in sequencing strategies and bioinformatics pipelines, indicating its broad application to vaginal microbiota. We further describe the relationships between community characteristics (vaginal pH, Nugent score) and participant demographics (race, age) and the CSTs defined by VALENCIA. CONCLUSION: VALENCIA provides a much-needed solution for the robust and reproducible assignment of vaginal community state types. This will allow unbiased analysis of both small and large vaginal microbiota datasets, comparisons between datasets and meta-analyses that combine multiple datasets. Video abstract.

Complete Genome Sequences of Six Lactobacillus iners Strains Isolated from the Human Vagina.

Lactobacillus iners is a common member of the human vaginal microbiota, with a genome size smaller than that of other lactobacilli. Here, we report the complete genome sequences of six L. iners strains isolated from different vaginal swab specimens. Three strains were found to harbor ∼100-kbp plasmids, which were not known previously.

Comparative Metagenome-Assembled Genome Analysis of "Candidatus Lachnocurva vaginae", Formerly Known as Bacterial Vaginosis-Associated Bacterium-1 (BVAB1).

Bacterial vaginosis-associated bacterium 1 (BVAB1) is an as-yet uncultured bacterial species found in the human vagina that belongs to the family Lachnospiraceae within the order Clostridiales. As its name suggests, this bacterium is often associated with bacterial vaginosis (BV), a common vaginal disorder that has been shown to increase a woman's risk for HIV, Chlamydia trachomatis, and Neisseria gonorrhoeae infections as well as preterm birth. BVAB1 has been further associated with the persistence of BV following metronidazole treatment, increased vaginal inflammation, and adverse obstetrics outcomes. There is no available complete genome sequence of BVAB1, which has made it difficult to mechanistically understand its role in disease. We present here a circularized metagenome-assembled genome (cMAG) of BVAB1 as well as a comparative analysis including an additional six metagenome-assembled genomes (MAGs) of this species. These sequences were derived from cervicovaginal samples of seven separate women. The cMAG was obtained from a metagenome sequenced with long-read technology on a PacBio Sequel II instrument while the others were derived from metagenomes sequenced on the Illumina HiSeq platform. The cMAG is 1.649 Mb in size and encodes 1,578 genes. We propose to rename BVAB1 to "Candidatus Lachnocurva vaginae" based on phylogenetic analyses, and provide genomic and metabolomic evidence that this candidate species may metabolize D-lactate, produce trimethylamine (one of the chemicals responsible for BV-associated odor), and be motile. The cMAG and the six MAGs are valuable resources that will further contribute to our understanding of the heterogeneous etiology of bacterial vaginosis.

A comprehensive non-redundant gene catalog reveals extensive within-community intraspecies diversity in the human vagina.

Analysis of metagenomic and metatranscriptomic data is complicated and typically requires extensive computational resources. Leveraging a curated reference database of genes encoded by members of the target microbiome can make these analyses more tractable. In this study, we assemble a comprehensive human vaginal non-redundant gene catalog (VIRGO) that includes 0.95 million non-redundant genes. The gene catalog is functionally and taxonomically annotated. We also construct a vaginal orthologous groups (VOG) from VIRGO. The gene-centric design of VIRGO and VOG provides an easily accessible tool to comprehensively characterize the structure and function of vaginal metagenome and metatranscriptome datasets. To highlight the utility of VIRGO, we analyze 1,507 additional vaginal metagenomes, and identify a high degree of intraspecies diversity within and across vaginal microbiota. VIRGO offers a convenient reference database and toolkit that will facilitate a more in-depth understanding of the role of vaginal microorganisms in women's health and reproductive outcomes.

Spatial structure facilitates the accumulation and persistence of antibiotic-resistant mutants in biofilms.

The emergence and spread of antibiotic resistance in bacterial pathogens are a global crisis. Because many bacterial infections are caused by pathogens that reside in biofilms, we sought to investigate how biofilms influence the evolution of antibiotic resistance. We hypothesize that the inherent spatial structure of biofilms facilitates the accumulation and persistence of spontaneously evolved antibiotic-resistant mutants. To test this, we tracked the frequency of mutants resistant to kanamycin and rifampicin in biofilm populations of Escherichia coli before, during, and after an antibiotic treatment regimen. Our results show that biofilms accumulate resistant mutants even in the absence of antibiotics. This resistance was found to be heritable and thus unlike the phenotypic plasticity of so-called "persister cells" that have been shown to occur in biofilms. Upon exposure to an antibiotic, resistant mutants swept to high frequency. Following the conclusion of treatment, these resistant mutants remained at unexpectedly high frequencies in the biofilms for over 45 days. In contrast, when samples from kanamycin-treated biofilms were used to found well-mixed liquid cultures and propagated by serial transfer, the frequency of resistant cells dramatically decreased as they were outcompeted by sensitive clones. These observations suggest that the emergence of antibiotic resistance through spontaneous mutations in spatially structured biofilms may significantly contribute to the emergence and persistence of mutants that are resistant to antibiotics used to treat bacterial infections.

The Relationship between Spatial Structure and the Maintenance of Diversity in Microbial Populations.

Spatial structure is pervasive in the microbial world, yet we know little about how it influences the evolution of microbial populations. It is thought that spatial structure limits the scale of competitive interactions and protracts selective sweeps. This may allow microbial populations to simultaneously explore multiple evolutionary paths. But how structured a microbial population must be before this effect is realized is not known. We used empirical and simulation studies to explore the relationship between spatial structure and the maintenance of diversity. The degree of spatial structure experienced by Escherichia coli metapopulations was manipulated by varying the migration rate between its component subpopulations. Each subpopulation was inoculated with an equal number of two equally fit genotypes, and their frequencies in 12 subpopulations were determined during 150 generations of evolution. We observed that the frequency of the "loser" genotypes decreased exponentially as the migration rate between the subpopulations was increased and that higher frequencies of the loser genotypes were maintained in structured metapopulations. These results demonstrate that structured microbial populations can evolve along multiple evolutionary trajectories even when migration rates between the subpopulations are relatively high.

The composition of the vaginal microbiome in first trimester pregnant women influences the level of autophagy and stress in vaginal epithelial cells.

Epithelial cells lining the vagina are major components of genital tract immunity. The influence of the vaginal microbiome on properties of host epithelial cells is largely unexplored. We evaluated whether differences in the most abundant lactobacilli species or bacterial genera in the vagina of first trimester pregnant women were associated with variations in the extent of stress and autophagy in vaginal epithelial cells. Vaginal swabs from 154 first trimester pregnant women were analyzed for bacterial composition by amplification and sequencing of the V1-V3 region of bacterial 16S rRNA genes. Vaginal epithelial cells were lysed and autophagy quantitated by measurement of p62. Intracellular levels of the inducible 70kDa heat shock protein (hsp70), an indicator of cell stress and an autophagy inhibitor, were determined. When Lactobacillus crispatus was the most abundant member of the vaginal microbiota, epithelial p62 and hsp70 levels were lowest as compared to when other bacterial taxa were most abundant. The highest concentrations of p62 and hsp70 were associated with Streptococcus and Bifidobacterium abundance. The p62 level associated with Gardnerella abundance was lower than that observed when lactobacilli other than L. crispatus were most abundant. In conclusion, in the first trimester of pregnancy the abundance of different bacterial taxa is associated with variations in autophagy and magnitude of the stress response in vaginal epithelial cells.

Genomic Comparisons of Lactobacillus crispatus and Lactobacillus iners Reveal Potential Ecological Drivers of Community Composition in the Vagina.

Lactobacillus crispatus and Lactobacillus iners are common inhabitants of the healthy human vagina. These two species are closely related and are thought to perform similar ecological functions in the vaginal environment. Temporal data on the vaginal microbiome have shown that nontransient instances of cooccurrence are uncommon, while transitions from an L. iners-dominated community to one dominated by L. crispatus, and vice versa, occur often. This suggests that there is substantial overlap in the fundamental niches of these species. Given this apparent niche overlap, it is unclear how they have been maintained as common inhabitants of the human vagina. In this study, we characterized and compared the genomes of L. iners and L. crispatus to gain insight into possible mechanisms driving the maintenance of this species diversity. Our results highlight differences in the genomes of these two species that may facilitate the partitioning of their shared niche space. Many of the identified differences may impact the protective benefits provided to the host by these two species. IMPORTANCE: The microbial communities that inhabit the human vagina play a critical role in the maintenance of vaginal health through the production of lactic acid and lowering the environmental pH. This precludes the growth of nonindigenous organisms and protects against infectious disease. The two most common types of vaginal communities are dominated by either Lactobacillus iners or Lactobacillus crispatus, while some communities alternate between the two over time. We combined ecological theory with state-of-the-art genome analyses to characterize how these two species might partition their shared niche space in the vagina. We show that the genomes of L. iners and L. crispatus differ in many respects, several of which may drive differences in their competitive abilities in the vagina. Our results provide insight into factors that drive the complicated temporal dynamics of the vaginal microbiome and demonstrate how closely related microbial species partition shared fundamental niche space.

Interference Competition Among Household Strains of Pseudomonas.

Bacterial species exhibit biogeographical patterns like those observed in larger organisms. The distribution of bacterial species is driven by environmental selection through abiotic and biotic factors as well dispersal limitations. We asked whether interference competition, a biotic factor, could explain variability in habitat use by Pseudomonas species in the human home. To answer this question, we screened almost 8000 directional, pairwise interactions between 89 Pseudomonas strains including members of the Pseudomonas aeruginosa (n = 29), Pseudomonas fluorescens (n = 21), and Pseudomonas putida (n = 39) species groups for the presence of killing. This diverse set of Pseudomonas strains includes those isolated from several different habitats within the home environment and includes combinations of strains that were isolated from different spatial scales. The use of this strain set not only allowed us to analyze the commonality and phylogenetic scale of interference competition within the genus Pseudomonas but also allowed us to investigate the influence of spatial scale on this trait. Overall, the probability of killing was found to decrease with increasing phylogenetic distance, making it unlikely that interference competition accounts for previously observed differential habitat use among Pseudomonas species and species groups. Strikingly, conspecific P. aeruginosa killing accounted for the vast majority of the observed killing, and this killing was found to differ across the habitat type and spatial scale of the strains' isolation. These data suggest that interference competition likely plays a large role in the within-species dynamics of P. aeruginosa but not other household Pseudomonas species.

Pseudomonas putida and Pseudomonas fluorescens Species Group Recovery from Human Homes Varies Seasonally and by Environment.

By shedding light on variation in time as well as in space, long-term biogeographic studies can help us define organisms' distribution patterns and understand their underlying drivers. Here we examine distributions of Pseudomonas in and around 15 human homes, focusing on the P. putida and P. fluorescens species groups. We describe recovery from 10,941 samples collected during up to 8 visits per home, occurring on average 2.6 times per year. We collected a mean of 141 samples per visit, from sites in most rooms of the house, from the surrounding yards, and from human and pet occupants. We recovered Pseudomonas in 9.7% of samples, with the majority of isolates being from the P. putida and P. fluorescens species groups (approximately 62% and 23% of Pseudomonas samples recovered respectively). Although representatives of both groups were recovered from every season, every house, and every type of environment sampled, recovery was highly variable across houses and samplings. Whereas recovery of P. putida group was higher in summer and fall than in winter and spring, P. fluorescens group isolates were most often recovered in spring. P. putida group recovery from soils was substantially higher than its recovery from all other environment types, while higher P. fluorescens group recovery from soils than from other sites was much less pronounced. Both species groups were recovered from skin and upper respiratory tract samples from healthy humans and pets, although this occurred infrequently. This study indicates that even species that are able to survive under a broad range of conditions can be rare and variable in their distributions in space and in time. For such groups, determining patterns and causes of stochastic and seasonal variability may be more important for understanding the processes driving their biogeography than the identity of the types of environments in which they can be found.