Sexuality in a natural population of bacteria– Bacillus subtilis challenges the clonal paradigm

Forty Years without Family: Three Novel Bacteriophages with High Similarity to SPP1 Reveal Decades of Evolutionary Stasis since the Isolation of Their Famous Relative

SPP1, an extensively studied bacteriophage of the Gram-positive Bacillus subtilis, is a model system for the study of phage-host interactions. Despite progress in the isolation and characterization of Bacillus phages, no previously fully sequenced phages have shared more than passing genetic similarity to SPP1. Here, we describe three virulent phages very similar to SPP1; SPP1 has greater than 80% nucleotide sequence identity and shares more that 85% of its protein coding genes with these phages. This is remarkable, given more than 40 years between the isolation of SPP1 and these phages. All three phages have somewhat larger genomes and more genes than SPP1. We identified a new putative gene in SPP1 based on a conserved sequence found in all phages. Gene conservation connotes purifying selection and is observed in structural genes and genes involved in DNA metabolism, but also in genes of unknown function, suggesting an important role in phage survival independent of the environment. Patterns of divergence point to genes or gene domains likely involved in adaptation to diverse hosts or different environments. Ultimately, comparative genomics of related phages provides insight into the long-term selective pressures that affect phage-bacteria interactions and alter phage genome content.

FINE-SCALE GENETIC AND PHENOTYPIC STRUCTURE IN NATURAL POPULATIONS OF BACILLUS SUBTILIS AND BACILLUS LICHENIFORMIS: IMPLICATIONS FOR BACTERIAL EVOLUTION AND SPECIATION

The genetic and phenotypic structure of sympatric populations of wild bacteria traditionally identified as Bacillus subtilis and B. licheniformis was analyzed. Small soil samples were taken from a single, tiny site in the Sonoran Desert of Arizona, USA, to provide a true population analysis, in contrast to many analyses of genetic structure using bacterial strain collections of widely heterogeneous origin. Genetic analyses of isolates used multilocus enzyme electrophoresis, mismatches in restriction fragment length polymorphism, and variants from Southern hybridization with B. subtilis DNA probes. Phenotypic analyses of isolates used the API test system for detection of growth and acid production on specific carbon sources. The two species were distinct both phenotypically and genetically, despite their known potential for genetic exchange in laboratory experiments. Genic and genotypic diversity were high in both species, and only 16% of observed allozyme variants might possibly be common to both species. Hence, there is probably modest genetic exchange, if any, between the species in nature. Clear hierarchies of population-genetic structure were found for both species. Different types of genetic data yield concordant population structures for B. subtilis. For both species, two-locus and multilocus statistical analyses of linkage demonstrated modest to strong disequilibrium at the species level but truly panmictic subunits within each species. The evidence for extensive genetic recombination within these fine-scale subdivisions is unequivocal, indicating that the sexuality of these bacteria can be well expressed in nature. The relation of these results to processes of bacterial evolution and speciation is discussed.

LONG-TERM EXPERIMENTAL EVOLUTION IN ESCHERICHIA COLI. VII. MECHANISMS MAINTAINING GENETIC VARIABILITY WITHIN POPULATIONS

Six replicate populations of the bacterium Escherichia coli were propagated for more than 10,000 generations in a defined environment. We sought to quantify the variation among clones within these populations with respect to their relative fitness, and to evaluate the roles of three distinct population genetic processes in maintaining this variation. On average, a pair of clones from the same population differed from one another in their relative fitness by approximately 4%. This within-population variation was small compared with the average fitness gain relative to the common ancestor, but it was statistically significant. According to one hypothesis, the variation in fitness is transient and reflects the ongoing substitution of beneficial alleles. We used Fisher's fundamental theorem to compare the observed rate of each population's change in mean fitness with the extent of variation for fitness within that population, but we failed to discern any correspondence between these quantities. A second hypothesis supposes that the variation in fitness is maintained by recurrent deleterious mutations that give rise to a mutation-selection balance. To test this hypothesis, we made use of the fact that two of the six replicate populations had evolved mutator phenotypes, which gave them a genomic mutation rate approximately 100-fold higher than that of the other populations. There was a marginally significant correlation between a population's mutation rate and the extent of its within-population variance for fitness, but this correlation was driven by only one population (whereas two of the populations had elevated mutation rates). Under a third hypothesis, this variation is maintained by frequency-dependent selection, whereby genotypes have an advantage when they are rare relative to when they are common. In all six populations, clones were more fit, on average, when they were rare than when they were common, although the magnitude of the advantage when rare was usually small (~1% in five populations and ~5% in the other). These three hypotheses are not mutually exclusive, but frequency-dependent selection appears to be the primary force maintaining the fitness variation within these experimental populations.

Evolution in the Bacillaceae

The family Bacillaceae constitutes a phenotypically diverse and globally ubiquitous assemblage of bacteria. Investigation into how evolution has shaped, and continues to shape, this family has relied on several widely ranging approaches from classical taxonomy, ecological field studies, and evolution in soil microcosms to genomic-scale phylogenetics, laboratory, and directed evolution experiments. One unifying characteristic of the Bacillaceae , the endospore, poses unique challenges to answering questions regarding both the calculation of evolutionary rates and claims of extreme longevity in ancient environmental samples.

Vectored dispersal of Symbiodinium by larvae of a Caribbean gorgonian octocoral

The ability of coral reefs to recover from natural and anthropogenic disturbance is difficult to predict, in part due to uncertainty regarding the dispersal capabilities and connectivity of their reef inhabitants. We developed microsatellite markers for the broadcast spawning gorgonian octocoral Eunicea (Plexaura) flexuosa (four markers) and its dinoflagellate symbiont, Symbiodinium B1 (five markers), and used them to assess genetic connectivity, specificity and directionality of gene flow among sites in Florida, Panama, Saba and the Dominican Republic. Bayesian analyses found that most E. flexuosa from the Florida reef tract, Saba and the Dominican Republic were strongly differentiated from many E. flexuosa in Panama, with the exception of five colonies from Key West that clustered with colonies from Panama. In contrast, Symbiodinium B1 was more highly structured. At least seven populations were detected that showed patterns of isolation by distance. The symbionts in the five unusual Key West colonies also clustered with symbionts from Panama, suggesting these colonies are the result of long-distance dispersal. Migration rate tests indicated a weak signal of northward immigration from the Panama population into the lower Florida Keys. As E. flexuosa clonemates only rarely associated with the same Symbiodinium B1 genotype (and vice versa), these data suggest a dynamic host-symbiont relationship in which E. flexuosa is relatively well dispersed but likely acquires Symbiodinium B1 from highly structured natal areas prior to dispersal. Once vectored by host larvae, these symbionts may then spread through the local population, and/or host colonies may acquire different local symbiont genotypes over time.

Diversity of culturable thermo-resistant aquatic bacteria along an environmental gradient in Cuatro Ciénegas, Coahuila, México

At the desert oasis of Cuatro Ciénegas in Coahuila, México, more than 300 oligotrophic pools can be found and a large number of endemic species of plants and animals. The most divergent taxa of diatoms, snail and fishes are located in the Churince hydrological system, where we analyzed the local diversification of cultivable Firmicutes and Actinobacteria. The Churince hydrological system is surrounded by gypsum dunes and has a strong gradient for salinity, temperature, pH and dissolved oxygen. In August 2003, surface water samples were taken in 10 sites along the Churince system together with the respective environmental measurements. 417 thermo-resistant bacteria were isolated and DNA was extracted to obtain their BOX-PCR fingerprints, revealing 55 different patterns. In order to identify similarities and differences in the diversity of the various sampling sites, an Ordination Analysis was applied using Principal Component Analysis. This analysis showed that conductivity is the environmental factor that explains the distribution of most of the microbial diversity. Phylogenetic reconstruction from their 16S rRNA sequences was performed for a sample of 150 isolates. Only 17 sequences had a 100% match in the Gene Bank (NCBI), representing 10 well known cosmopolitan taxa. The rest of the sequences cluster in 22 clades for Firmicutes and another 22 clades for Actinobacteria, supporting the idea of high diversity and differentiation for this site.

Genetic Variation within a Lotic Population of Janthinobacterium lividum

An understanding of the genetic variation within and between populations should allow scientists to address many problems, including those associated with endangered species and the release of genetically modified organisms into the environment. With respect to microorganisms, the release of genetically engineered microorganisms is likely to increase dramatically given the current growth in the bioremediation industry. In this study, genetic variation within a lotic, bacterial population of Janthinobacterium lividum was measured with restriction fragment length polymorphism analysis. Chromosomal DNA from 10 Kettle Creek (Hawk Mountain Sanctuary, Kempton, Pa.) J. lividum isolates was digested with six restriction endonucleases and probed with a 7.5-kb pKK3535 fragment containing the E. coli rrnB rRNA operon. Genetic variation, as measured in terms of nucleotide diversity, was high within the population. The 0.0781 value for genetic variation was especially high given the conservative nature of the genetic probe. The average percent similarity among isolates within the population was 67.25%. Pairwise comparisons of nucleotide diversity values (pi) and similarity coefficients (F) yielded values ranging from 0.0032 to 0.1816 and 0.3363 to 0.9808, respectively. Putative clonemates were not present within the group of isolates; however, all isolates shared 14 fragments across a spectrum of six restriction enzymes. The presence of these common fragments indicates that restriction fragment length polymorphism analysis may provide population- or species-specific diagnostic markers for J. lividum. Data that suggest a plume effect with respect to the downstream movement of J. lividum are also presented. An increase in genetic variation within groups of isolates along the longitudinal gradient of Kettle Creek is also suggested.

RAPD analysis of genetic diversity and qualitative assessment of hydrolytic activities in a collection of Bacillus sp. isolate

Genetic diversity and production of hydrolytic enzymes of 205 Bacillus isolates from different geographical and ecological niches in Serbia were studied. Combining RAPD analysis and 16S DNA sequencing, we determined 13 different groups of RAPD profiles within four (five) species: B. subtilis, B. cereus/B. thuringiensis, B. pumilus, and B. firmus. Screening for production of hydrolytic enzymes showed that there was no correlation of enzyme production with species. Most of the isolates from all habitats produced amylase, protease, lipase, mannanase, and xylanase to some extent at 25?C and 37?C. The number of isolates that retained enzyme production ability at 55?C is considerably lower and they predominantly came from manure.

Social interactions and distribution of Bacillus subtilis pherotypes at microscale

Bacillus subtilis strains communicate through the comQXPA quorum sensing (QS) system, which regulates genes expressed during early stationary phase. A high polymorphism of comQXP' loci was found in closely related strains isolated from desert soil samples separated by distances ranging from meters to kilometers. The observed polymorphism comprised four communication groups (pherotypes), such that strains belonging to the same pherotype exchanged information efficiently but strains from different pherotypes failed to communicate. To determine whether the same level of polymorphism in the comQXP' QS system could be detected at microscale, B. subtilis isolates were obtained from two separate 1-cm(3) soil samples, which were progressively divided into smaller sections. Cross-activation studies using pherotype-responsive reporter strains indicated the same number of communication pherotypes at microscale as previously determined at macroscale. Sequencing of the housekeeping gene gyrA and the QS comQ gene confirmed different evolutionary rates of these genes. Furthermore, an asymmetric communication response was detected inside the two pherotype clusters, suggesting continuous evolution of the QS system and possible development of new languages. To our knowledge, this is the first microscale study demonstrating the presence of different QS languages among isolates of one species, and the implications of this microscale diversity for microbial interactions are discussed.

Ecology and genomics of Bacillus subtilis

Bacillus subtilis is a remarkably diverse bacterial species that is capable of growth within many environments. Recent microarray-based comparative genomic analyses have revealed that members of this species also exhibit considerable genomic diversity. The identification of strain-specific genes might explain how B. subtilis has become so broadly adapted. The goal of identifying ecologically adaptive genes could soon be realized with the imminent release of several new B. subtilis genome sequences. As we embark upon this exciting new era of B. subtilis comparative genomics we review what is currently known about the ecology and evolution of this species.

Multilocus sequence typing (MLST) reveals high genetic diversity and clonal population structure of the toxic cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is one of the most prevalent bloom-forming cyanobacteria and has been the cause of increasing public health concern due to the production of hepatotoxins (microcystins). To investigate the genetic diversity, clonality and evolutionary genetic background with regard to the toxicity of M. aeruginosa, a multilocus sequence typing (MLST) scheme was developed, based on seven selected housekeeping loci (ftsZ, glnA, gltX, gyrB, pgi, recA and tpi). Analysis of a collection of 164 isolates from Japan and other countries identified 79 unique sequence types (STs), revealing a high level of genetic diversity (H=0.951). Although recombination between loci was indicated to be substantial by Shimodaira-Hasegawa (SH) tests, multilocus linkage disequilibrium analyses indicated that recombination between strains probably occurs at some frequency but not to the extent at which alleles are associated randomly, suggesting that the population structure of M. aeruginosa is clonal. Analysis of subsets of strains also indicated that the clonal population structure is maintained even in a local population. Phylogenetic analysis based on the concatenated sequences of seven MLST loci demonstrated that microcystin-producing genotypes are not monophyletic, providing further evidence for the gain and loss of toxicity during the intraspecific diversification of M. aeruginosa. However, toxic strains are genetically distinct from non-toxic strains in MLST allelic profiles, and it was also shown that non-toxic strains harbouring toxin genes fall into a single monophyletic clade, except for one case. These results suggest that the toxicity of M. aeruginosa is relatively stable in the short term, and therefore can be unequivocally characterized by MLST. The MLST scheme established here will be of great help for future detailed population genetic studies of M. aeruginosa.

Causes of insertion sequences abundance in prokaryotic genomes

Insertion sequences (ISs) are the smallest and most frequent transposable elements in prokaryotes where they play an important evolutionary role by promoting gene inactivation and genome plasticity. Their genomic abundance varies by several orders of magnitude for reasons largely unknown and widely speculated. The current availability of hundreds of genomes renders testable many of these hypotheses, notably that IS abundance correlates positively with the frequency of horizontal gene transfer (HGT), genome size, pathogenicity, nonobligatory ecological associations, and human association. We thus reannotated ISs in 262 prokaryotic genomes and tested these hypotheses showing that when using appropriate controls, there is no empirical basis for IS family specificity, pathogenicity, or human association to influence IS abundance or density. HGT seems necessary for the presence of ISs, but cannot alone explain the absence of ISs in more than 20% of the organisms, some of which showing high rates of HGT. Gene transfer is also not a significant determinant of the abundance of IS elements in genomes, suggesting that IS abundance is controlled at the level of transposition and ensuing natural selection and not at the level of infection. Prokaryotes engaging in obligatory associations have fewer ISs when controlled for genome size, but this may be caused by some being sexually isolated. Surprisingly, genome size is the only significant predictor of IS numbers and density. Alone, it explains over 40% of the variance of IS abundance. Because we find that genome size and IS abundance correlate negatively with minimal doubling times, we conclude that selection for rapid replication cannot account for the few ISs found in small genomes. Instead, we show evidence that IS numbers are controlled by the frequency of highly deleterious insertion targets. Indeed, IS abundance increases quickly with genome size, which is the exact inverse trend found for the density of genes under strong selection such as essential genes. Hence, for ISs, the bigger the genome the better.

Bacillus subtilis genome diversity

Microarray-based comparative genomic hybridization (M-CGH) is a powerful method for rapidly identifying regions of genome diversity among closely related organisms. We used M-CGH to examine the genome diversity of 17 strains belonging to the nonpathogenic species Bacillus subtilis. Our M-CGH results indicate that there is considerable genetic heterogeneity among members of this species; nearly one-third of Bsu168-specific genes exhibited variability, as measured by the microarray hybridization intensities. The variable loci include those encoding proteins involved in antibiotic production, cell wall synthesis, sporulation, and germination. The diversity in these genes may reflect this organism's ability to survive in diverse natural settings.

Genetic diversity and biogeography of haloalkaliphilic sulphur-oxidizing bacteria belonging to the genus Thioalkalivibrio

A group of 85 isolates of haloalkaliphilic obligately chemolithoautotrophic sulphur-oxidizing bacteria belonging to the genus Thioalkalivibrio were recently obtained from soda lakes in Mongolia, Kenya, California, Egypt and Siberia. They have been analyzed by repetitive extragenic palindromic (rep)-PCR genomic fingerprinting technique with BOX- and (GTG)5-primer set. Cluster analysis was performed using combined fingerprint profiles and a dendrogram similarity value (r) of 0.8 was used to define the same genotype. Fifty-six genotypes were found among the isolates, revealing a high genetic diversity. The strains can be divided into two major clusters, including isolates from the Asiatic (Siberia and Mongolia) and the African (Kenya and Egypt) continents, respectively. The majority (85.9%) of the genotypes were found in only one area, suggesting an endemic character of the Thioalkalivibrio strains. Furthermore, a correlation between fingerprint clustering, geographic origin and the characteristics of the lake of origin was found.

Ecological separation and genetic isolation of Neisseria gonorrhoeae and Neisseria meningitidis

BACKGROUND

Classifying bacteria into species is problematic. Most microbiologists consider species to be groups of isolates that share some arbitrary degree of relatedness of biochemical or molecular (such as DNA sequence) features and that, ideally, are clearly delineated from all other groups of isolates. The main problem in applying to bacteria a biological concept of species based on the ability or inability of their genes to recombine, is that recombination appears to be rare in bacteria in nature, as indicated by the strong linkage disequilibrium between alleles found in most bacterial populations. However, there are some naturally transformable bacteria in which assortative recombination appears to be so frequent that alleles are in, or close to, linkage equilibrium. For these recombining populations a biological concept of species might be applicable.

RESULTS

Populations of Neisseria gonorrhoeae and Neisseria meningitidis from Spain were analysed by multilocus enzyme electrophoresis. The data indicate that assortative recombination occurs frequently within populations, but not between populations. Similarly, the sequences of two house-keeping genes show no evidence of intragenic recombination between N. gonorrhoeae and N. meningitidis.

CONCLUSIONS

N. gonorrhoeae and N. meningitidis represent extremely closely related 'sexual' populations that appear to be genetically isolated in nature, and thus conform to the biological concept of species. The extreme uniformity of N. gonorrhoeae house-keeping genes suggests that this species may have arisen recently as a clone of N. meningitidis that could colonize the genital tract. Ecological isolation - of populations that can colonize the genital tract from those that can colonize the nasopharynx - may have been an important component in speciation, leading to a lower frequency of recombination between species than within species.

Genetic diversity, recombination and cryptic clades in Pseudomonas viridiflava infecting natural populations of Arabidopsis thaliana

Species-level genetic diversity and recombination in bacterial pathogens of wild plant populations have been nearly unexplored. Pseudomonas viridiflava is a common natural bacterial pathogen of Arabidopsis thaliana, for which pathogen defense genes and mechanisms are becoming increasing well known. The genetic variation contained within a worldwide sample of P. viridiflava collected from wild populations of A. thaliana was investigated using five genomic sequence fragments totaling 2.3 kb. Two distinct and deeply diverged clades were found within the P. viridiflava sample and in close proximity in multiple populations, each genetically diverse with synonymous variation as high as 9.3% in one of these clades. Within clades, there is evidence of frequent recombination within and between each sequenced locus and little geographic differentiation. Isolates from both clades were also found in a small sample of other herbaceous species in Midwest populations, indicating a possibly broad host range for P. viridiflava. The high levels of genetic variation and recombination together with a lack of geographic differentiation in this pathogen distinguish it from other bacterial plant pathogens for which intraspecific variation has been examined.

Relationship between spatial and genetic distance in Agrobacterium spp. in 1 cubic centimeter of soil

The spatial and genetic unit of bacterial population structure is the clone. Surprisingly, very little is known about the spread of a clone (spatial distance between clonally related bacteria) and the relationship between spatial distance and genetic distance, especially at very short scale (microhabitat scale), where cell division takes place. Agrobacterium spp. Biovar 1 was chosen because it is a soil bacterial taxon easy to isolate. A total of 865 microsamples 500 microm in diameter were sampled with spatial coordinates in 1 cm(3) of undisturbed soil. The 55 isolates obtained yielded 42 ribotypes, covering three genomic species based on amplified ribosomal DNA restriction analysis (ARDRA) of the intergenic spacer 16S-23S, seven of which contained two to six isolates. These clonemates (identical ARDRA patterns) could be found in the same microsample or 1 cm apart. The genetic diversity did not change with distance, indicating the same habitat variability across the cube. The mixing of ribotypes, as assessed by the spatial position of clonemates, corresponded to an overlapping of clones. Although the population probably was in a recession stage in the cube (10(3) agrobacteria g(-1)), a high genetic diversity was maintained. In two independent microsamples (500 microm in diameter) at the invasion stage, the average genetic diversity was at the same level as in the cube. Quantification of the microdiversity landscape will help to estimate the probability of encounter between bacteria under realistic natural conditions and to set appropriate sampling strategies for population genetic analysis.

Identification of complex composition, strong strain diversity and directional selection in local Pseudomonas stutzeri populations from marine sediment and soils

Members of Pseudomonas stutzeri have been isolated world-wide from various habitats including aquatic and terrestrial ecosystems. The global population has a clonal structure, is of exceptionally high genetic diversity and has been grouped into eight genomovars. We have analysed four local populations (n = 89-125) from three geographically separated habitats (two from a marine sediment and two from different soils) by random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR), restriction fragment length polymorphism (RFLP) of the rpoB gene and 16S rDNA sequences in order to quantify the influence of evolutionary forces on closely related groups of proliferating cells in situ. All populations consisted of a complex structure of genomic subgroups with variable numbers of members. The analyses revealed that the two populations from marine sediment were rather similar. At least three of the populations were influenced by migrational input as concluded from the presence of members from different genomovars. All populations showed very high strain diversity suggesting strong influence of mutability. Neutrality tests indicated that two or possibly three of the populations were shaped by directional selection. Thus, the local populations of P. stutzeri reflect already the high genetic diversity of the global population and are influenced, to different extents, by migration, mutation and directional selection.

Genotypic and phenotypic diversity within species of purple nonsulfur bacteria isolated from aquatic sediments

To assess the extent of genotypic and phenotypic diversity within species of purple nonsulfur bacteria found in aquatic sediments, a total of 128 strains were directly isolated from agar plates that had been inoculated with sediment samples from Haren and De Biesbosch in The Netherlands. All isolates were initially characterized by BOX-PCR genomic DNA fingerprinting, and 60 distinct genotypes were identified. Analyses of 16S rRNA gene sequences of representatives of each genotype showed that five and eight different phylotypes of purple nonsulfur bacteria were obtained from the Haren and De Biesbosch sites, respectively. At the Haren site, 80.5% of the clones were Rhodopseudomonas palustris, whereas Rhodoferax fermentans and Rhodopseudomonas palustris were numerically dominant at the De Biesbosch site and constituted 45.9 and 34.4% of the isolates obtained, respectively. BOX-PCR genomic fingerprints showed that there was a high level of genotypic diversity within each of these species. The genomic fingerprints of Rhodopseudomonas palustris isolates were significantly different for isolates from the two sampling sites, suggesting that certain strains may be endemic to each sampling site. Not all Rhodopseudomonas palustris isolates could degrade benzoate, a feature that has previously been thought to be characteristic of the species. There were differences in the BOX-PCR genomic fingerprints and restriction fragment length polymorphisms of benzoate-coenzyme A ligase genes and form I and form II ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) genes between benzoate-degrading and non-benzoate-degrading genotypes. The ability to distinguish these two Rhodopseudomonas palustris groups based on multiple genetic differences may reflect an incipient speciation event resulting from adaptive evolution to local environmental conditions.

Genetic differentiation between sympatric populations of Bacillus cereus and Bacillus thuringiensis

Little is known about genetic exchanges in natural populations of bacteria of the spore-forming Bacillus cereus group, because no population genetics studies have been performed with local sympatric populations. We isolated strains of Bacillus thuringiensis and B. cereus from small samples of soil collected at the same time from two separate geographical sites, one within the forest and the other at the edge of the forest. A total of 100 B. cereus and 98 B. thuringiensis strains were isolated and characterized by electrophoresis to determine allelic composition at nine enzymatic loci. We observed genetic differentiation between populations of B. cereus and B. thuringiensis. Populations of a given Bacillus species--B. thuringiensis or B. cereus--were genetically more similar to each other than to populations of the other Bacillus species. Hemolytic activity provided further evidence of this genetic divergence, which remained evident even if putative clones were removed from the data set. Our results suggest that the rate of gene flow was higher between strains of the same species, but that exchanges between B. cereus and B. thuringiensis were nonetheless possible. Linkage disequilibrium analysis revealed sufficient recombination for B. cereus populations to be considered panmictic units. In B. thuringiensis, the balance between clonal proliferation and recombination seemed to depend on location. Overall, our data indicate that it is not important for risk assessment purposes to determine whether B. cereus and B. thuringiensis belong to a single or two species. Assessment of the biosafety of pest control based on B. thuringiensis requires evaluation of the extent of genetic exchange between strains in realistic natural conditions.

The structure of a local population of phytopathogenic Pseudomonas brassicacearum from agricultural soil indicates development under purifying selection pressure

Among the isolates of a bacterial community from a soil sample taken from an agricultural plot in northern Germany, a population consisting of 119 strains was obtained that was identified by 16S rDNA sequencing and genomic fingerprinting as belonging to the recently described species Pseudomonas brassicacearum. Analysis of the population structure by allozyme electrophoresis (11 loci) and random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR; four primers) showed higher resolution with the latter method. Both methods indicated the presence of three lineages, one of which dominated strongly. Stochastic tests derived from the neutral theory of evolution (including Slatkin's exact test, Watterson's homozygosity test and the Tajima test) indicated that the population had developed under strong purifying selection pressure. The presence of strains clearly divergent from the majority of the population can be explained by in situ evolution or by influx of strains as a result of migration or both. Phytopathogenicity of a P. brassicacearum strain determined with tomato plants reached the level obtained with the type strain of the known pathogen Pseudomonas corrugata. The results show that a selective sweep was identified in a local population. Previously, a local selective sweep had not been seen in several populations of different bacterial species from a variety of environmental habitats.

Biogeography and degree of endemicity of fluorescent Pseudomonas strains in soil

Fluorescent Pseudomonas strains were isolated from 38 undisturbed pristine soil samples from 10 sites on four continents. A total of 248 isolates were confirmed as Pseudomonas sensu stricto by fluorescent pigment production and group-specific 16S ribosomal DNA (rDNA) primers. These isolates were analyzed by three molecular typing methods with different levels of resolution: 16S rDNA restriction analysis (ARDRA), 16S-23S rDNA intergenic spacer-restriction fragment length polymorphism (ITS-RFLP) analysis, and repetitive extragenic palindromic PCR genomic fingerprinting with a BOX primer set (BOX-PCR). All isolates showed very similar ARDRA patterns, as expected. Some ITS-RFLP types were also found at every geographic scale, although some ITS-RFLP types were unique to the site of origin, indicating weak endemicity at this level of resolution. Using a similarity value of 0.8 or more after cluster analysis of BOX-PCR fingerprinting patterns to define the same genotypes, we identified 85 unique fluorescent Pseudomonas genotypes in our collection. There were no overlapping genotypes between sites as well as continental regions, indicating strict site endemism. The genetic distance between isolates as determined by degree of dissimilarity in BOX-PCR patterns was meaningfully correlated to the geographic distance between the isolates' sites of origin. Also, a significant positive spatial autocorrelation of the distribution of the genotypes was observed among distances of <197 km, and significant negative autocorrelation was observed between regions. Hence, strong endemicity of fluorescent Pseudomonas genotypes was observed, suggesting that these heterotrophic soil bacteria are not globally mixed.

The potential for intraspecific horizontal gene exchange by natural genetic transformation: sexual isolation among genomovars of Pseudomonas stutzeri

The potential for natural genetic transformation among the seven genomovars (gvs) of Pseudomonas stutzeri was investigated. Of the 12 strains originating from a variety of environments, six strains (50%) from five gvs were competent for DNA uptake (Rif(R) marker). The transformation frequencies varied over more than three orders of magnitude. With three highly transformable strains (ATCC 17587, ATCC 17641, JM300) from two gvs and all other strains as DNA donors, sexual isolation from other pseudomonad species (Pseudomonas alcaligenes, Pseudomonas mendocina) and also from other P. stutzeri gvs was observed (i.e. heterogamic transformation was reduced). For ATCC 17587 (gv 2) and ATCC 17641 (gv 8), heterogamic transformation was up to two and three orders of magnitude lower with other P. stutzeri gv and the other species employed, respectively, than in homogamic transformations. Interestingly, whereas with ATCC 17587 and ATCC 17641 heterogamic transformation with donors of the same gv was as high as homogamic transformation, JM300 (gv 8) was sexually isolated from its nearest relative (ATCC 17641). Also, sexual isolation of JM300 from other P. stutzeri gvs was most pronounced among the recipients tested, in some cases reaching the highest levels found with the other species as DNA donors (reduction of heterogamic transformation by 4000-fold). Results obtained here from nucleotide sequence analysis of part (422 nt) of the gene for the RNA polymerase ss subunit (rpoB) from various strains indicated that sexual isolation of ATCC 17641 increased with nucleotide sequence divergence. Implications of the observed great heterogeneity in transformability, competence levels and sexual isolation among strains are discussed with regard to the evolution of P. stutzeri.

Allele-specific PCR shows that genetic exchange occurs among genetically diverse Nodularia (cyanobacteria) filaments in the Baltic Sea

Some cyanobacteria have been shown to exchange genetic information under laboratory conditions, but it has not been clear whether such genetic exchange occurs in the natural environment. To address this, a population genetic study was carried out on the filamentous diazotrophic cyanobacterium Nodularia in the Baltic Sea. Nodularia filaments were collected from 20 widely distributed sampling stations in the Baltic Sea during June and July 1998. Allele-specific PCR (AS-PCR) was used to characterize over 2000 filaments at three loci: a non-coding spacer between adjacent copies of the main structural gas vesicle gene gvpA (gvpA-IGS), the phycocyanin intergenic spacer (PC-IGS) and the rDNA internal transcribed spacer (rDNA-ITS). The three loci were all found to be polymorphic in the 1998 population: two alternative alleles were distinguished at the gvpA-IGS and PC-IGS loci, and three at the rDNA-ITS locus. All 12 possible combinations of alleles were found in the filaments studied, but some were much more common than others. The index of association (I:(A)) for all possible pairwise combinations of isolates was found to differ significantly from zero, which implies that there is some linkage disequilibrium between loci. The I:(A) values for 16 out of 20 individual sampling stations also differed significantly from zero: this shows that the observed linkage disequilibrium is not due to pooling data from genetically distinct subpopulations. Monte-Carlo simulations with random subsets of the data confirmed that some combinations showed significantly more linkage disequilibrium than expected by chance alone. It is concluded that genetic exchange occurs in the natural Nodularia population, but the frequency is not high enough for the loci to be in linkage equilibrium. The distribution of the 12 genotypes across the Baltic Sea was found to be non-random, but did not correlate with temperature, salinity or major nutrient concentrations. A significant relationship was found between the gene diversity among filaments at each station and the distance of the station from the centre of the sampling area: possible reasons for this trend are discussed.

Estimating recombinational parameters in Streptococcus pneumoniae from multilocus sequence typing data

Multilocus sequence typing (MLST) is a highly discriminatory molecular typing method that defines isolates of bacterial pathogens using the sequences of approximately 450-bp internal fragments of seven housekeeping genes. This technique has been applied to 575 isolates of Streptococcus pneumoniae and identifies a number of discrete clonal complexes. These clonal complexes are typically represented by a single group of isolates sharing identical alleles at all seven loci, plus single-locus variants that differ from this group at only one out of the seven loci. As MLST is highly discriminatory, the members of each clonal complex can be assumed to have a recent common ancestor, and the molecular events that give rise to the single-locus variants can be used to estimate the relative contributions of recombination and mutation to clonal divergence. By comparing the sequences of the variant alleles within each clonal complex with the allele typically found within that clonal complex, we estimate that recombination has generated new alleles at a frequency approximately 10-fold higher than mutation, and that a single nucleotide site is approximately 50 times more likely to change through recombination than mutation. We also demonstrate how to estimate the average length of recombinational replacements from MLST data.

Genetic structure of natural populations of Escherichia coli in wild hosts on different continents

Current knowledge of genotypic and phenotypic diversity in the species Escherichia coli is based almost entirely on strains recovered from humans or zoo animals. In this study, we analyzed a collection of 202 strains obtained from 81 mammalian species representing 39 families and 14 orders in Australia and the Americas, as well as several reference strains; we also included a strain from a reptile and 10 from different families of birds collected in Mexico. The strains were characterized genotypically by multilocus enzyme electrophoresis (MLEE) and phenotypically by patterns of sugar utilization, antibiotic resistance, and plasmid profile. MLEE analysis yielded an estimated genetic diversity (H) of 0.682 for 11 loci. The observed genetic diversity in this sample is the greatest yet reported for E. coli. However, this genetic diversity is not randomly distributed; geographic effects and host taxonomic group accounted for most of the genetic differentiation. The genetic relationship among the strains showed that they are more associated by origin and host order than is expected by chance. In a dendrogram, the ancestral cluster includes primarily strains from Australia and ECOR strains from groups B and C. The most differentiated E. coli in our analysis are strains from Mexican carnivores and strains from humans, including those in the ECOR group A. The kinds and numbers of sugars utilized by the strains varied by host taxonomic group and country of origin. Strains isolated from bats were found to exploit the greatest range of sugars, while those from primates utilized the fewest. Toxins are more frequent in strains from rodents from both continents than in any other taxonomic group. Strains from Mexican wild mammals were, on average, as resistant to antibiotics as strains from humans in cities. On average, the Australian strains presented a lower antibiotic resistance than the Mexican strains. However, strains recovered from hosts in cities carried significantly more plasmids than did strains isolated from wild mammals. Previous studies have shown that natural populations of E. coli harbor an extensive genetic diversity that is organized in a limited number of clones. However, knowledge of this worldwide bacterium has been limited. Here, we suggest that the strains from a wide range of wild hosts from different regions of the world are organized in an ecotypic structure where adaptation to the host plays an important role in the population structure.

Recombination and clonal groupings within Helicobacter pylori from different geographical regions

A collection of 20 strains of Helicobacter pylori from several regions of the world was studied to better understand the population genetic structure and diversity of this species. Sequences of fragments from seven housekeeping genes (atpA, efp, mutY, ppa, trpC, ureI, yphC ) and two virulence-associated genes (cagA, vacA) showed high levels of synonymous sequence variation (mean percentage Ks of 10-27%) and lower levels of non-synonymous variation (mean percentage Ka of 0.2-5.6%). Cluster analysis of pairwise differences between alleles revealed the existence of two weakly clonal groupings, which included half of the strains investigated. All six strains isolated from Japanese and coastal Chinese were assigned to the 'Asian' clonal grouping, probably reflecting descent from a distinct common ancestor. The clonal groupings were not totally uniform; recombination, as measured by the homoplasy test and compatibility matrices, was extremely common within all genes tested, except cagA. The fact that clonal descent could still be discerned despite such frequent recombination possibly reflects founder effects and geographical separation and/or selection for particular alleles of these genes.

Detecting linkage disequilibrium in bacterial populations

The distribution of the number of pairwise differences calculated from comparisons between n haploid genomes has frequently been used as a starting point for testing the hypothesis of linkage equilibrium. For this purpose the variance of the pairwise differences, VD, is used as a test statistic to evaluate the null hypothesis that all loci are in linkage equilibrium. The problem is to determine the critical value of the distribution of VD. This critical value can be estimated either by Monte Carlo simulation or by assuming that VD is distributed normally and calculating a one-tailed 95% critical value for VD, L, L = EVD + 1.645 sqrt(VarVD), where E(VD) is the expectation of VD, and Var(VD) is the variance of VD. If VD (observed) > L, the null hypothesis of linkage equilibrium is rejected. Using Monte Carlo simulation we show that the formula currently available for Var(VD) is incorrect, especially for genetically highly diverse data. This has implications for hypothesis testing in bacterial populations, which are often genetically highly diverse. For this reason we derive a new, exact formula for Var(VD). The distribution of VD is examined and shown to approach normality as the sample size increases. This makes the new formula a useful tool in the investigation of large data sets, where testing for linkage using Monte Carlo simulation can be very time consuming. Application of the new formula, in conjunction with Monte Carlo simulation, to populations of Bradyrhizobium japonicum, Rhizobium leguminosarum, and Bacillus subtilis reveals linkage disequilibrium where linkage equilibrium has previously been reported.

Free recombination within Helicobacter pylori

Sequences of three gene fragments (flaA, flaB, and vacA) from Helicobacter pylori strains isolated from patients in Germany, Canada, and South Africa were analyzed for diversity and for linkage equilibrium by using the Homoplasy Test and compatibility matrices. Horizontal genetic exchange in H. pylori is so frequent that different loci and polymorphisms within each locus are all at linkage equilibrium. These results indicate that H. pylori is panmictic. Comparisons with sequences from Escherichia coli, Neisseria meningitidis, and Drosophila melanogaster showed that recombination in H. pylori was much more frequent than in other species. In contrast, when multiple family members infected with H. pylori were investigated, some strains were indistinguishable at all three loci. Thus, H. pylori is clonal over short time periods after natural transmission.

Random association among alleles in clonal populations of Sclerotinia sclerotiorum

Multiple loci identified in DNA fingerprints were used to test for random association in two agricultural populations of S. sclerotiorum. In linkage disequilibrium tests among pairs of loci with frequencies between 0.1 and 0.9, 44.5 and 80.5% of pairs of loci were consistent with random association in the clone-corrected samples of the Canadian canola and the North Carolina cabbage populations, respectively. In estimates of corrected (Bonferroni) P value, 70.66 and 98.89% of pairs of loci were in random association. All four possible genotypes for each pair of loci were observed in the Canadian canola sample, consistent with random association among loci. In multilocus association tests across all loci, however, significant association was observed in both populations. In the Canadian canola population, 40 possible heterokaryons were identified. Our data suggest that populations of S. sclerotiorum are predominantly clonal and that occasional genetic exchange and recombination, and not mutation alone, may be a source of new genotypes.

Biodiversity of a Burkholderia cepacia population isolated from the maize rhizosphere at different plant growth stages

A Burkholderia cepacia population naturally occurring in the rhizosphere of Zea mays was investigated in order to assess the degree of root association and microbial biodiversity at five stages of plant growth. The bacterial strains isolated on semiselective PCAT medium were mostly assigned to the species B. cepacia by an analysis of the restriction patterns produced by amplified DNA coding for 16S rRNA (16S rDNA) (ARDRA) with the enzyme AluI. Partial 16S rDNA nucleotide sequences of some randomly chosen isolates confirmed the ARDRA results. Throughout the study, B. cepacia was strictly associated with maize roots, ranging from 0.6 to 3.6% of the total cultivable microflora. Biodiversity among 83 B. cepacia isolates was analyzed by the random amplified polymorphic DNA (RAPD) technique with two 10-mer primers. An analysis of RAPD patterns by the analysis of molecular variance method revealed a high level of intraspecific genetic diversity in this B. cepacia population. Moreover, the genetic diversity was related to divergences among maize root samplings, with microbial genetic variability markedly higher in the first stages of plant growth; in other words, the biodiversity of this rhizosphere bacterial population decreased over time.

Differences in genetic diversity of nonecapsulated Haemophilus influenzae from various diseases

Genetic relationships among 80 isolates of nonencapsulated Haemophilus influenzae recovered from different disease types were determined by multilocus enzyme electrophoresis (MEE) at 13 enzyme loci in an attempt to assess the association between multilocus genotype and disease. The isolates were obtained from 15 patients with meningitis, 10 with otitis media, 19 with chronic bronchitis, 20 with cystic fibrosis, and 16 were obtained from healthy carriers. The 80 isolates were assigned to 69 electrophoretic types (ETs) falling into 5 groups. Isolates from each disease entity were represented by a variety of genotypes; however, cluster analysis from a matrix of genetic distances between ETs revealed that the ETs of the otitis media and meningitis isolates were all clustered within a genetic distance of 0.55 (group I). In addition, no genotypes were shared between H. influenzae carrier isolates and isolates from cases of disease, H. influenzae isolates from healthy individuals were distributed significantly differently from those from chronic bronchitis meningitis and otitis media patients. The genetic diversity (H) of carrier strains was greatest, although not statistically different from that of isolates from patients with disease. It was concluded that the genetic distribution of acute disease isolates is not random over the five ET groups, although the genetic diversity within the groups is not different. The effect of bacterial persistence in the host on the genetic diversity of H. influenzae is discussed.

Evidence of parasexual exchange of DNA in the rice blast fungus challenges its exclusive clonality

ABSTRACT We applied DNA markers to determine whether parasexual recombination may contribute to the extreme genetic diversity and variability observed in Magnaporthe grisea, the causal agent of rice blast disease. Dispersed repetitive elements and mapped, low-copy restriction fragment length polymorphism (RFLP) probes were used to detect transfers of DNA between cultured isolates of M. grisea. Low-copy RFLP probes also were used to detect putative recombinants among isolates from well-characterized field populations of the pathogen. Microscopic examination of tufted mycelium between cocultured isolates revealed frequent hyphal fusions. Hyphal tips and conidia were recovered without selection from tufted zones in two separate vegetative pairings involving isolates with dissimilar haplotypes, based on the repetitive element MGR586. Haplotypic changes were observed at a higher frequency in tuft derivatives than in subcultures of each isolate alone. From 136 tuft derivatives analyzed, 5 putative recombinant haplotypes were identified. Introgression was demonstrated with two independent repetitive elements, fosbury and MGR586, as probes on DNA digested with several restriction enzymes. Introgressions were characterized by addition of 1 to 10 MGR586 bands, and 1 to 3 fosbury bands from one parent into the background of the other. Polymorphic single-copy probes were used to analyze putative recombinants. One probe detected an introgression event as predicted by analysis with MGR586. To assess the possible role of parasexual recombination in field populations of the pathogen, isolates in the Philippines previously grouped based on DNA fingerprinting were analyzed with low-copy RFLP markers. Polymorphism in single-copy loci typically was seen between, but not within, putative pathogen lineages. One lineage (designated lineage 4), however, was polymorphic for several probes. For some isolates, alleles at these loci comigrated with alleles characteristic of other lineages, suggesting the transfer of DNA fragments between lineages. One isolate was apparently a merodiploid, carrying an allele typical of lineage 4 plus another allele characteristic of a different lineage. In a survey of isolates from the Indian Himalayas, a merodiploid also was found with single- or low-copy probes. Examination of MGR586 profiles of the putative recombinant and its putative donor strains showed the expected introgression of MGR586 bands. The detection of parasexual DNA exchanges in wild-type strains under unselected conditions and the existence of merodiploids in nature suggest that parasexual recombination occurs in field populations of M. grisea. This raises questions concerning exclusive clonality in the blast fungus.

A comparison of the nucleotide sequences of the adk and recA genes of pathogenic and commensal Neisseria species: evidence for extensive interspecies recombination within adk

The sequences of the adenylate kinase gene (adk) and the RecA gene (recA) were determined from the same isolates of Neisseria gonorrhoeae, N. meningitidis, N. lactamica, N. polysaccharea, N. cinerea, N. mucosa, N. pharyngis var. flava, N. flavescens, and N. animalis. The patterns of sequence divergence observed at adk and recA were very different. Dendrograms constructed from the recA data using two different algorithms were statistically robust and were congruent with each other and with the relationships between the species previously proposed using other data. In contrast, the dendrograms derived from the adk data were noncogruent with each other, and with those from the recA data, and were statistically poorly supported. These results, along with the uniform distribution of pairwise sequence divergences between the species at adk, suggest there has been a history of interspecies recombination within the adk gene of the human Neisseria species which has obscured the phylogenetic relationships between the species. This view was supported by Sawyer's runs test, and the Index of Association (IA) between codons, which provided significant evidence for interspecies recombination between the adk genes from the human Neisseria species, but no evidence of interspecies recombination between the recA sequences.

Towards a unified evolutionary genetics of microorganisms

I propose here that evolutionary genetics, apart from improving our basic knowledge of the taxonomy and evolution of microbes (either eukaryotes or prokaryotes), can also greatly contribute to applied research in microbiology. Evolutionary genetics provides convenient guidelines for better interpreting genetic and molecular data dealing with microorganisms. The three main potential applications of evolutionary genetics in microbiology are (a) epidemiological follow-up (with the necessity of evaluating the stability of microbial genotypes over space and time); (b) taxonomy in the broad sense (better definition and sharper delimitation of presently described taxa, research of hidden genetic subdivisions); and (c) evaluation of the impact of the genetic diversity of microbes on their relevant properties (pathogenicity, resistance to drugs, etc). At present, two main kinds of population structure can be distinguished in natural microbial populations: (a) species that are not subdivided into discrete phylogenetic lineages (panmictic species or basically sexual species with occasional bouts of short-term clonality fall into this category); (b) species that are strongly subdivided by either cryptic speciation or clonal evolution. Improvements in available statistical methods are required to refine these distinctions and to better quantify the actual impact of gene exchange in natural microbial populations. Moreover, a codified selection of markers with appropriate molecular clocks (in other words: adapted levels of resolution) is sorely needed to answer distinct questions that address different scales of time and space: experimental, epidemic, and evolutionary. The problems raised by natural genetic diversity are very similar for all microbial species, in terms of both basic and applied science. Despite this fact, a regrettable compartmentalization among specialists has hampered progress in this field. I propose a synthetic approach, relying on the statistical improvements and technical standardizations called for above, to settle a unified evolutionary genetics of microorganisms, valid whatever the species studied, whether eukaryotic (parasitic protozoa and fungi) or prokaryotic (bacteria). Apart from benefits for basic evolutionary research, the anticipated payoff from this synthetic approach is to render routine and common-place the use of microbial evolutionary genetics in the fields of epidemiology, medicine, and agronomy.

Temporal Variation in Genetic Diversity and Structure of a Lotic Population of Burkholderia (Pseudomonas) cepacia

The genetic structure and temporal patterns of genetic diversity in a population of Burkholderia (Pseudomonas) cepacia, isolated from a southeastern blackwater stream, were investigated by multilocus enzyme electrophoresis. Allelic variation in seven structural gene loci was monitored at a single stream location at 0, 6, 12, and 24 h and at 2, 4, 8, 16, and 32 days. Over the length of the study, 217 isolates were collected, from which 65 unique electrophoretic types (ETs) were identified. Most of these ETs were present at only one or two time periods and were considered transients; however, one resident ET was particularly abundant (64 of the 217 isolates [29.4%]) and was found at all time points except day 32. The mean genetic diversity of the entire population was 0.520, and the index of association (a measure of multilocus linkage disequilibrium) was 1.33. These results, taken in conjunction with a previous study focusing on spatial patterns of genetic diversity in lotic B. cepacia, show that these bacterial populations exhibit greater variability among sites than within a site over time, suggesting relative stability over short time periods.

A hierarchical analysis of population genetic structure in Rhizobium leguminosarum bv. trifolii

Little is known about the population processes that shape the genetic diversity in natural populations of rhizobia. A sample of 912 Rhizobium leguminosarum biovar trifolii isolates were collected from naturalized red clover populations (Trifolium pratense) and analysed for 15 allozyme loci to determine the levels and distribution of genetic diversity. Hierarchical analyses compared different sampling levels, geographical separation, and temporal separation. Total genetic diversity across all isolates was H = 0.426, with 57.6% of the total diversity found among isolates obtained from individual red clover plants. Relatively low genetic differentiation among populations and high differentiation among plants within populations was observed; this suggests that gene flow and founder effect act differently at geographical and local scales. Significant differences were observed in (i) allele frequencies among populations and among plants within populations, and (ii) the frequency distribution of the most widespread and the most abundant strains. When multilocus linkage disequilibrium was calculated, significant levels of disequilibrium were observed in the total sample and in three of the eight populations.