Large-scale phylogenetic analyses reveal the causes of high tropical amphibian diversity

Many groups show higher species richness in tropical regions but the underlying causes remain unclear. Despite many competing hypotheses to explain latitudinal diversity gradients, only three processes can directly change species richness across regions: speciation, extinction and dispersal. These processes can be addressed most powerfully using large-scale phylogenetic approaches, but most previous studies have focused on small groups and recent time scales, or did not separate speciation and extinction rates. We investigate the origins of high tropical diversity in amphibians, applying new phylogenetic comparative methods to a tree of 2871 species. Our results show that high tropical diversity is explained by higher speciation in the tropics, higher extinction in temperate regions and limited dispersal out of the tropics compared with colonization of the tropics from temperate regions. These patterns are strongly associated with climate-related variables such as temperature, precipitation and ecosystem energy. Results from models of diversity dependence in speciation rate suggest that temperate clades may have lower carrying capacities and may be more saturated (closer to carrying capacity) than tropical clades. Furthermore, we estimate strikingly low tropical extinction rates over geological time scales, in stark contrast to the dramatic losses of diversity occurring in tropical regions presently.     

Vaginal microbiome variances in sample groups categorized by clinical criteria of bacterial vaginosis

Background

One of the most common and recurrent vaginal infections is bacterial vaginosis (BV). The diagnosis is based on changes to the “normal” vaginal microbiome; however, the normal microbiome appears to differ according to reproductive status and ethnicity, and even among individuals within these groups. The Amsel criteria and Nugent score test are widely used for diagnosing BV; however, these tests are based on different criteria, and so may indicate distinct changes in the vaginal microbial community. Nevertheless, few studies have compared the results of these test against metagenomics analysis.

Methods

Vaginal flora samples from 77 participants were classified according to the Amsel criteria and Nugent score test. The microbiota composition was analyzed using 16S ribosome RNA gene amplicon sequencing. Bioinformatics analysis and multivariate statistical analysis were used to evaluate the microbial diversity and function.

Results

Only 3 % of the participants diagnosed BV negative using the Amsel criteria (A−) were BV-positive according to the Nugent score test (N+), while over half of the BV-positive patients using the Amsel criteria (A+) were BV-negative according to the Nugent score test (N−). Thirteen genera showed significant differences in distribution among BV status defined by BV tests (e.g., A − N−, A + N− and A + N+). Variations in the four most abundant taxa, Lactobacillus, Gardnerella, Prevotella, and Escherichia, were responsible for most of this dissimilarity. Furthermore, vaginal microbial diversity differed significantly among the three groups classified by the Nugent score test (N−, N+, and intermediate flora), but not between the Amsel criteria groups. Numerous predictive microbial functions, such as bacterial chemotaxis and bacterial invasion of epithelial cells, differed significantly among multiple BV test, but not between the A− and A+ groups.

Conclusions

Metagenomics analysis can greatly expand our current understanding of vaginal microbial diversity in health and disease. Metagenomics profiling may also provide more reliable diagnostic criteria for BV testing.

     

Microbial communities of urban stormwater sediments: the phylogenetic structure of bacterial communities varies with porosity

This study focuses on the distribution of bacterial and fungal communities within the microstructure of a multi-contaminated sedimentary layer resulting from urban stormwater infiltration. Fractionation was performed on the basis of differential porosity and aggregate grain size, resulting in five fractions: leachable fitting macroporosity, < 10, 10-160, 160-1000 μm fitting aggregates, > 1000 μm. Amounts of both bacterial and fungal biomasses are greater in the < 10 μm and leachable fractions. The aggregates contain numerous bacteria but very low amounts of fungal biomass. Single-strand conformational polymorphism molecular profiles highlighted the differences between bacterial and fungal communities of the leachable fraction and those of the aggregates. Random Sanger sequencing of ssu clones revealed that these differences were mainly because of the presence of Epsilonproteobacteria and Firmicutes in the leachable fractions, while the aggregates contained more Cyanobacteria. The Cyanobacteria phylotypes in the aggregates were dominated by the sequences related to Microcoleus vaginatus while the leachable fractions presented the sequences of chloroplastic origin. Therefore, more than 50% of the phylotypes observed were related to Proteobacteria while 40% were related to Cyanobacteria and Bacteroidetes. Preferential distribution of clades in almost all the phyla or classes detected was observed. This study provides insight into the identities of dominant members of the bacterial communities of urban sediments. Microcoleus vaginatus appeared to predominate in pioneer soils.     

The complex vaginal flora of West African women with bacterial vaginosis

Background

The spectrum of bacteria associated with bacterial vaginosis (BV) has recently expanded through taxonomic changes and the use of molecular methods. These methods have yet to be used in large-scale epidemiological studies in Africa where BV is highly prevalent.

Methods

An analysis of samples obtained during a clinical trial of the management of vaginal discharge in four West African countries. Samples were available from 1555 participants; 843 (54%) had BV. Nucleic acids of 13 bacterial genera or species potentially associated with BV were detected through the polymerase chain reaction.

Results

The associations between various components of the vaginal flora were complex. Excluding Lactobacillus, the other 12 micro-organisms were all associated with each other at the p≤0.001 level. The prevalence of various bacterial genera or species varied according to age, sexual activity and HIV status. In multivariate analysis, the presence of Gardnerella vaginalis, Bifidobacterium, Megasphaera elsdenii, Dialister, Mycoplasma hominis, Leptotrichia, and Prevotella were independently associated with BV as was the absence of Lactobacillus and Peptoniphilus. However, Mobiluncus, Atopobium vaginae, Anaerococcus, and Eggerthella were not independently associated with BV. Unexpectedly, after treatment with a regimen that included either metronidazole or tinidazole, the proportion of patients with a complete resolution of symptoms by day 14 increased with the number of bacterial genera or species present at enrolment.

Conclusions

Numerous bacterial genera or species were strongly associated with each other in a pattern that suggested a symbiotic relationship. BV cases with a simpler flora were less likely to respond to treatment. Overall, the vaginal flora of West African women with BV was reminiscent of that of their counterparts in industrialized countries.     

A survey for biochemical synapomorphies to reveal phylogenetic relationships of halichondrid demosponges (Metazoa: Porifera)

Biochemical compounds are a suggested alternative for morphological characters in sponge systematics. The monophyly of the crucial demosponge order Halichondrida is still largely unresolved. Here, we performed a search for synapomorphies, which might aid to define the order and its five families Axinellidae, Bubaridae, Dictyonellidae, Desmoxyidae and Halichondriidae. We also investigated possible links of the order Halichondrida with other demosponge taxa such as suberitids (Hadromerida). Our survey of published compounds revealed hydroxy-nor-sterols as a potential synapomorphy for axinellid taxa, 2,3,6 sulfated sterol nuclei as possible connection between Agelasida and Halichondrida, besides the known pyrrol-2-carboxylic acids and galactosylceramides, furthermore carbonimidic dichlorides, pupukeananes, diterpene isocyanides for the order Halichondrida and linear and cyclic diterpenes for the family Desmoxyidae. The suitability of biochemical characters in systematics is discussed.     

Increased prostaglandin concentrations in the cervical mucus of pregnant women with bacterial vaginosis.

Microorganisms associated with bacterial vaginosis are commonly recovered from the amniotic fluid and chorion-amnion of patients who deliver prematurely. Bacteria closely related to those causing bacterial vaginosis may play a role in the initiation of uterine contractions, ripening of the cervix and weakening of the fetal membranes by stimulating prostaglandin synthesis. In the present investigation, cervical mucus was collected by brush from early pregnant women with and without bacterial vaginosis. The concentrations of PGE2, PGF2 alpha and 6-keto-PGF1 alpha were determined in the mucus samples by methyl oximation and then radioimmunoassay, utilizing antibodies raised against oximated prostaglandins. It was found that the concentration of PGE2 and PGF2 alpha was significantly higher in the mucus of women with bacterial vaginosis compared with healthy women. The concentration of 6-keto-PGF1 alpha was similar in both study groups. All patients had been instructed to abstain from sexual intercourse for 24 hours before sampling. However, it may be that women with high concentrations in their mucus may have had intercourse anyway. However, it is fairly well possible that the significant differences in the PGE2 and PGF2 alpha values are causally related to the higher rate of preterm labor in women with the commonplace infection of bacterial vaginosis.     

Resolving deep phylogenetic relationships in salamanders: analyses of mitochondrial and nuclear genomic data

Phylogenetic relationships among salamander families illustrate analytical challenges inherent to inferring phylogenies in which terminal branches are temporally very long relative to internal branches. We present new mitochondrial DNA sequences, approximately 2,100 base pairs from the genes encoding ND1, ND2, COI, and the intervening tRNA genes for 34 species representing all 10 salamander families, to examine these relationships. Parsimony analysis of these mtDNA sequences supports monophyly of all families except Proteidae, but yields a tree largely unresolved with respect to interfamilial relationships and the phylogenetic positions of the proteid genera Necturus and Proteus. In contrast, Bayesian and maximum-likelihood analyses of the mtDNA data produce a topology concordant with phylogenetic results from nuclear-encoded rRNA sequences, and they statistically reject monophyly of the internally fertilizing salamanders, suborder Salamandroidea. Phylogenetic simulations based on our mitochondrial DNA sequences reveal that Bayesian analyses outperform parsimony in reconstructing short branches located deep in the phylogenetic history of a taxon. However, phylogenetic conflicts between our results and a recent analysis of nuclear RAG-1 gene sequences suggest that statistical rejection of a monophyletic Salamandroidea by Bayesian analyses of our mitochondrial genomic data is probably erroneous. Bayesian and likelihood-based analyses may overestimate phylogenetic precision when estimating short branches located deep in a phylogeny from data showing substitutional saturation; an analysis of nucleotide substitutions indicates that these methods may be overly sensitive to a relatively small number of sites that show substitutions judged uncommon by the favored evolutionary model.     

Microbial fingerprinting detects unique bacterial communities in the faecal microbiota of rats with experimentally-induced colitis

An abnormal composition of the gut microbiota is believed to be associated with the pathogenesis of inflammatory bowel disease (IBD). We utilized terminal restriction fragment length polymorphism (T-RFLP) analysis to quantify faecal bacterial communities from rats with experimental colitis. Male Sprague Dawley rats (n=10/group) ingested 2% dextran sulfate sodium (DSS) or water for up to 7 days. Rats were killed and colonic tissues collected for histological analysis. Damage severity score in the distal colon was significantly greater (P<0.001) following DSS consumption compared to controls. T-RFLP faecal bacterial profiles generated with either MspI or CfoI revealed a significant difference (P<0.001) in community composition between healthy and colitic rats, with bacterial composition in healthy rats more variable than in rats with colitis. Operational taxonomic units (OTU: taxonomically related groups of bacteria) associated with either the healthy or colitic state were identified. OTU (116, 226, 360, and 948; CfoI) and (118 and 188; MspI) were strongly associated with untreated healthy rats, while OTU (94, 98, 174, and 384; CfoI) and (94 and 914; MspI) were predominantly associated with DSS-treated colitic rats. Phylogenetic OTU assignment suggested that Bacteroidales and Lactobacillus sp. were predominantly associated with the colitic and healthy rats, respectively. These results show that faecal bacterial profiling is a rapid, sensitive and non-invasive tool for detecting and identifying changes in gut microbiota associated with colitis. Restoring microbial homeostasis by targeting colitis-associated OTU through specific microbiological interventions could form the basis of novel therapeutic strategies for IBD.     

Bacterial vaginosis among a group of married Jordanian women: occurrence and laboratory diagnosis

A total of 310 vaginal swabs collected from a group of married Jordanian women complaining of vaginal discharge were examined for bacterial vaginosis. The scoring system of Nugent for the interpretation of Gram staining was employed. This system revealed the presence of the condition in 29.7% of patients. Results obtained using the scoring system correlated significantly with the detection of clue cells and the scarcity of white blood cells in the vaginal discharge. An inverse relationship was found between bacterial vaginosis and Lactobacillus morphotypes determined by Gram staining. No definite relationship was detected between bacterial vaginosis and the recovery of Gardnerella vaginalis by culture as this organism was isolated from swabs which according to the Nugent criterion were negative for bacterial vaginosis. Bacterial vaginosis among the women investigated was more prevalent than vaginitis caused by Trichomonas vaginalis or yeasts.     

Cultivation-independent methods reveal differences among bacterial gut microbiota in triatomine vectors of Chagas disease

Background

Chagas disease is a trypanosomiasis whose agent is the protozoan parasite Trypanosoma cruzi, which is transmitted to humans by hematophagous bugs known as triatomines. Even though insecticide treatments allow effective control of these bugs in most Latin American countries where Chagas disease is endemic, the disease still affects a large proportion of the population of South America. The features of the disease in humans have been extensively studied, and the genome of the parasite has been sequenced, but no effective drug is yet available to treat Chagas disease. The digestive tract of the insect vectors in which T. cruzi develops has been much less well investigated than blood from its human hosts and constitutes a dynamic environment with very different conditions. Thus, we investigated the composition of the predominant bacterial species of the microbiota in insect vectors from Rhodnius, Triatoma, Panstrongylus and Dipetalogaster genera.

Methodology/Principal Findings

Microbiota of triatomine guts were investigated using cultivation-independent methods, i.e., phylogenetic analysis of 16s rDNA using denaturing gradient gel electrophoresis (DGGE) and cloned-based sequencing. The Chao index showed that the diversity of bacterial species in triatomine guts is low, comprising fewer than 20 predominant species, and that these species vary between insect species. The analyses showed that Serratia predominates in Rhodnius, Arsenophonus predominates in Triatoma and Panstrongylus, while Candidatus Rohrkolberia predominates in Dipetalogaster.

Conclusions/Significance

The microbiota of triatomine guts represents one of the factors that may interfere with T. cruzi transmission and virulence in humans. The knowledge of its composition according to insect species is important for designing measures of biological control for T. cruzi. We found that the predominant species of the bacterial microbiota in triatomines form a group of low complexity whose structure differs according to the vector genus.     

Dynamics of vaginal bacterial communities in women developing bacterial vaginosis, candidiasis, or no infection, analyzed by PCR-denaturing gradient gel electrophoresis and real-time PCR

The microbial flora of the vagina plays a major role in preventing genital infections, including bacterial vaginosis (BV) and candidiasis (CA). An integrated approach based on PCR-denaturing gradient gel electrophoresis (PCR-DGGE) and real-time PCR was used to study the structure and dynamics of bacterial communities in vaginal fluids of healthy women and patients developing BV and CA. Universal eubacterial primers and Lactobacillus genus-specific primers, both targeted at 16S rRNA genes, were used in DGGE and real-time PCR analysis, respectively. The DGGE profiles revealed that the vaginal flora was dominated by Lactobacillus species under healthy conditions, whereas several potentially pathogenic bacteria were present in the flora of women with BV. Lactobacilli were the predominant bacterial population in the vagina for patients affected by CA, but changes in the composition of Lactobacillus species were observed. Real-time PCR analysis allowed the quantitative estimation of variations in lactobacilli associated with BV and CA diseases. A statistically significant decrease in the relative abundance of lactobacilli was found in vaginal fluids of patients with BV compared to the relative abundance of lactobacilli in the vaginal fluids of healthy women and patients with CA.     

Genetic and biochemical diversity of Gardnerella vaginalis strains isolated from women with bacterial vaginosis

Gardnerella vaginalis is considered a substantial player in the progression of bacterial vaginosis (BV). We analysed 17 G. vaginalis strains isolated from the genital tract of women diagnosed with BV to establish a potential link between genotypes/biotypes and the expression of virulence factors, vaginolysin (VLY) and sialidase, which are assumed to play a substantial role in the pathogenesis of BV. Amplified ribosomal DNA restriction analysis revealed two G. vaginalis genotypes. Gardnerella vaginalis isolates of genotype 2 appeared more complex than genotype 1 and were subdivided into three subtypes. Biochemical typing allowed us to distinguish four different biotypes. A great diversity of the level of VLY production among the isolates of G. vaginalis may be related to a different cytotoxicity level of the strains. We did not find any correlation between VLY production level and G. vaginalis genotype/biotype. In contrast, a link between G. vaginalis genotype and sialidase production was established. Our findings on the diversity of VLY expression level in different clinical isolates and linking sialidase activity with the genotype of G. vaginalis could help to evaluate the pathogenic potential of different G. vaginalis strains.     

Bacterial vaginosis (BV) candidate bacteria: associations with BV and behavioural practices in sexually-experienced and inexperienced women

Background

In recent years several new fastidious bacteria have been identified that display a high specificity for BV; however no previous studies have comprehensively assessed the behavioural risk associations of these bacterial vaginosis-candidate organisms (BV-COs).

Methods

We examined the associations between 8 key previously described BV-COs and BV status established by Nugent''s score (NS). We also examined the sexual practices associated with each BV-CO. We incorporated 2 study populations: 193 from a sexually-inexperienced university population and 146 from a highly sexually-active clinic population. Detailed behavioural data was collected by questionnaire and vaginal smears were scored by the Nugent method. Stored samples were tested by quantitative PCR assays for the 8 BV-COs: Atopobium vaginae, Gardnerella vaginalis, Leptotrichia spp., Megasphaera type I, Sneathia spp., and the Clostridia-like bacteria BVAB1, BVAB2 and BVAB3. Associations between BV-COs and BV and behaviours were examined by univariate and multivariable analyses.

Results

On univariate analysis, all BV-COs were more common in BV compared to normal flora. However, only Megasphaera type I, BVAB2, A. vaginae and G. vaginalis were significantly independently associated with BV by multivariable analysis. Six of the eight BV-COs (Megasphaera type I, BVAB2, BVAB3, Sneathia, Leptotrichia and G. vaginalis) were rare or absent in sexually-unexposed women, and demonstrated increasing odds of detection with increasing levels of sexual activity and/or numbers of lifetime sexual partners. Only G. vaginalis and A. vaginae were commonly detected in sexually-unexposed women. Megasphaera type I was independently associated with women-who-have-sex-with women (WSW) and lifetime sexual partner numbers, while unprotected penile-vaginal-sex was associated with BVAB2 detection by multivariate analysis.

Conclusions

Four of eight key BV-COs were significantly associated with BV after adjusting for the presence of other BV-COs. The majority of BV-COs were absent or rare in sexually-unexposed women, and associated with increasing sexual exposure, suggesting potential sexual transmission of BV-COs.     

Phytoplankton communities and antecedent conditions: high resolution sampling in Esthwaite Water

1. The succession of a phytoplankton community was investigated through an intensive period of sampling and related to physical, chemical and biological conditions sampled at an equal, or higher, temporal resolution. 2. Phytoplankton samples were taken on a weekly basis from June to September 2004 and analysed for diversity, species composition, and contribution of different functional groups to total biomass. Physical and chemical data were collected on the sampling days, and physical environmental factors were also logged continuously throughout the period by automatic measuring stations. This continuous logging allowed community structure to be compared with physical data averaged over periods from a day to a week before each sampling date. 3. The Schmidt stability of the lake, a measure of the strength of stratification calculated from thermal data, showed a negative correlation with phytoplankton species diversity. This is consistent with the hypothesis that mixing was preventing exclusion by species that would otherwise dominate in stratified conditions. 4. At a functional level, stress tolerant (S‐type) species dominated during the stratified summer conditions, with small, colonising species (C types) and ruderal, disturbance tolerant species (R types) contributing little to the overall biomass. Of the stress tolerant species, the faster growing (S_C) phytoplankters were significantly favoured by more stable, stratified conditions and higher solar radiation. Increased abundance of this group resulted in decreased species diversity. Correlations were generally strongest when using the 6‐ to 7‐day averaged physical data, stressing the importance of continuous measurements of these drivers in phytoplankton studies.     

Assessing the complex sponge microbiota: core,variable and species-specific bacterial communities in marine sponges

Marine sponges are well known for their associations with highly diverse, yet very specific and often highly similar microbiota. The aim of this study was to identify potential bacterial sub-populations in relation to sponge phylogeny and sampling sites and to define the core bacterial community. 16S ribosomal RNA gene amplicon pyrosequencing was applied to 32 sponge species from eight locations around the world''s oceans, thereby generating 2567 operational taxonomic units (OTUs at the 97% sequence similarity level) in total and up to 364 different OTUs per sponge species. The taxonomic richness detected in this study comprised 25 bacterial phyla with Proteobacteria, Chloroflexi and Poribacteria being most diverse in sponges. Among these phyla were nine candidate phyla, six of them found for the first time in sponges. Similarity comparison of bacterial communities revealed no correlation with host phylogeny but a tropical sub-population in that tropical sponges have more similar bacterial communities to each other than to subtropical sponges. A minimal core bacterial community consisting of very few OTUs (97%, 95% and 90%) was found. These microbes have a global distribution and are probably acquired via environmental transmission. In contrast, a large species-specific bacterial community was detected, which is represented by OTUs present in only a single sponge species. The species-specific bacterial community is probably mainly vertically transmitted. It is proposed that different sponges contain different bacterial species, however, these bacteria are still closely related to each other explaining the observed similarity of bacterial communities in sponges in this and previous studies. This global analysis represents the most comprehensive study of bacterial symbionts in sponges to date and provides novel insights into the complex structure of these unique associations.     

Lipid biomarker and phylogenetic analyses to reveal archaeal biodiversity and distribution in hypersaline microbial mat and underlying sediment

This study has utilized the tools of lipid biomarker chemistry and molecular phylogenetic analyses to assess the archaeal contribution to diversity and abundance within a microbial mat and underlying sediment from a hypersaline lagoon in Baja California. Based on abundance of ether-linked isoprenoids, archaea made up from 1 to 4% of the cell numbers throughout the upper 100 mm of mat and sediment core. Below this depth archaeal lipid was two times more abundant than bacterial. Archaeol was the primary archaeal lipid in all layers. Relatively small amounts of caldarchaeol (dibiphytanyl glyceroltetraether) were present at most depths with phytanyl to biphytanyl molar ratios lowest (~10 : 1) in the 4–17 mm and 100–130 mm horizons, and highest (132 : 1) in the surface 0–2 mm. Lipids with cyclic biphytanyl cores were only detected below 100 mm. A novel polar lipid containing a C₃₀ isoprenoid (squalane) moiety was isolated from the upper anoxic portion of the core and partially characterized. Hydrocarbon biomarker lipids included pentamethylicosane (2–10 mm) and crocetane (primarily below 10 mm). Archaeal molecular diversity varied somewhat with depth. With the exception of samples at 0–2 mm and 35–65 mm, Thermoplasmatales of marine benthic group D dominated clone libraries. A significant number of phylotypes representing the Crenarchaeota from marine benthic group B were generally present below 17 mm and dominated the 35–65 mm sample. Halobacteriaceae family made up 80% of the clone library of the surface 2 mm, and consisted primarily of sequences affiliated with the haloalkaliphilic Natronomonas pharaonis .     

Deep sequencing analyses of low density microbial communities: working at the boundary of accurate microbiota detection

Introduction

Accurate analyses of microbiota composition of low-density communities (10³–10⁴ bacteria/sample) can be challenging. Background DNA from chemicals and consumables, extraction biases as well as differences in PCR efficiency can significantly interfere with microbiota assessment. This study was aiming to establish protocols for accurate microbiota analysis at low microbial density.

Methods

To examine possible effects of bacterial density on microbiota analyses we compared microbiota profiles of serial diluted saliva and low (nares, nasopharynx) and high-density (oropharynx) upper airway communities in four healthy individuals. DNA was extracted with four different extraction methods (Epicentre Masterpure, Qiagen DNeasy, Mobio Powersoil and a phenol bead-beating protocol combined with Agowa-Mag-mini). Bacterial DNA recovery was analysed by 16S qPCR and microbiota profiles through GS-FLX-Titanium-Sequencing of 16S rRNA gene amplicons spanning the V5–V7 regions.

Results

Lower template concentrations significantly impacted microbiota profiling results. With higher dilutions, low abundant species were overrepresented. In samples of <10⁵ bacteria per ml, e.g. DNA <1 pg/µl, microbiota profiling deviated from the original sample and other dilutions showing a significant increase in the taxa Proteobacteria and decrease in Bacteroidetes. In similar low density samples, DNA extraction method determined if DNA levels were below or above 1 pg/µl and, together with lysis preferences per method, had profound impact on microbiota analyses in both relative abundance as well as representation of species.

Conclusion

This study aimed to interpret microbiota analyses of low-density communities. Bacterial density seemed to interfere with microbiota analyses at < than 10⁶ bacteria per ml or DNA <1 pg/µl. We therefore recommend this threshold for working with low density materials. This study underlines that bias reduction is crucial for adequate profiling of especially low-density bacterial communities.     

Use of cluster and discriminant analyses to compare rhizosphere bacterial communities following biological perturbation

We present an approach to comparing the diversity and composition of bacterial communities from different habitats and for identifying which members of a community are most affected by an introduced bacterium. We use this method to explore both previously published and new data from field and growth chamber experiments in which we isolated heterotrophic bacteria from samples of root-free soil, roots of nontreated soybean seedlings, and from the roots of soybean seedlings grown from Bacillus cereus UW85nl—treated seeds. We characterize bacterial isolates for 40 physiological attributes, and grouped the isolates hierarchically using two-stage density-linkage cluster analysis. Multivariate analysis of variance and discriminant analysis of the relative frequencies of the clusters in the soil and rhizosphere habitats were then used to determine whether there were differences among the bacterial communities from the various habitats, and which of the clusters were most useful in discriminating among the communities. We used rarefied estimates of richness as a measure of community diversity in the various habitats. Introduction of UW85n 1 affected the composition and/or diversity of rhizosphere communities in three of four experiments. Present address: Department of Environmental Science, Policy and Management, 108 Hilgard Hall, University of California, Berkeley, CA 94720 Correspondence to: G.S. Gilbert.