Rhizobial 16S rRNA and dnaK genes: mosaicism and the uncertain phylogenetic placement of Rhizobium galegae.

The phylogenetic relatedness among 12 agriculturally important species in the order Rhizobiales was estimated by comparative 16S rRNA and dnaK sequence analyses. Two groups of related species were identified by neighbor-joining and maximum-parsimony analysis. One group consisted of Mesorhizobium loti and Mesorhizobium ciceri, and the other group consisted of Agrobacterium rhizogenes, Rhizobium tropici, Rhizobium etli, and Rhizobium leguminosarum. Although bootstrap support for the placement of the remaining six species varied, A. tumefaciens, Agrobacterium rubi, and Agrobacterium vitis were consistently associated in the same subcluster. The three other species included Rhizobium galegae, Sinorhizobium meliloti, and Brucella ovis. Among these, the placement of R. galegae was the least consistent, in that it was placed flanking the A. rhizogenes-Rhizobium cluster in the dnaK nucleotide sequence trees, while it was placed with the other three Agrobacterium species in the 16S rRNA and the DnaK amino acid trees. In an effort to explain the inconsistent placement of R. galegae, we examined polymorphic site distribution patterns among the various species. Localized runs of nucleotide sequence similarity were evident between R. galegae and certain other species, suggesting that the R. galegae genes are chimeric. These results provide a tenable explanation for the weak statistical support often associated with the phylogenetic placement of R. galegae, and they also illustrate a potential pitfall in the use of partial sequences for species identification.

Sequence diversity and molecular evolution of the leukotoxin (lktA) gene in bovine and ovine strains of Mannheimia (Pasteurella) haemolytica.

The molecular evolution of the leukotoxin structural gene (lktA) of Mannheimia (Pasteurella) haemolytica was investigated by nucleotide sequence comparison of lktA in 31 bovine and ovine strains representing the various evolutionary lineages and serotypes of the species. Eight major allelic variants (1.4 to 15.7% nucleotide divergence) were identified; these have mosaic structures of varying degrees of complexity reflecting a history of horizontal gene transfer and extensive intragenic recombination. The presence of identical alleles in strains of different genetic backgrounds suggests that assortative (entire gene) recombination has also contributed to strain diversification in M. haemolytica. Five allelic variants occur only in ovine strains and consist of recombinant segments derived from as many as four different sources. Four of these alleles consist of DNA (52.8 to 96.7%) derived from the lktA gene of the two related species Mannheimia glucosida and Pasteurella trehalosi, and four contain recombinant segments derived from an allele that is associated exclusively with bovine or bovine-like serotype A2 strains. The two major lineages of ovine serotype A2 strains possess lktA alleles that have very different evolutionary histories and encode divergent leukotoxins (5.3% amino acid divergence), but both contain segments derived from the bovine allele. Homologous segments of donor and recipient alleles are identical or nearly identical, indicating that the recombination events are relatively recent and probably postdate the domestication of cattle and sheep. Our findings suggest that host switching of bovine strains from cattle to sheep, together with inter- and intraspecies recombinational exchanges, has played an important role in generating leukotoxin diversity in ovine strains. In contrast, there is limited allelic diversity of lktA in bovine strains, suggesting that transmission of strains from sheep to cattle has been less important in leukotoxin evolution.

Parallel evolution of virulence in pathogenic Escherichia coli.

The mechanisms underlying the evolution and emergence of new bacterial pathogens are not well understood. To elucidate the evolution of pathogenic Escherichia coli strains, here we sequenced seven housekeeping genes to build a phylogenetic tree and trace the history of the acquisition of virulence genes. Compatibility analysis indicates that more than 70% of the informative sites agree with a single phylogeny, suggesting that recombination has not completely obscured the remnants of ancestral chromosomes. On the basis of the rate of synonymous substitution for E. coli and Salmonella enterica (4.7 x 10(-9) per site per year), the radiation of clones began about 9 million years ago and the highly virulent pathogen responsible for epidemics of food poisoning, E. coli O157:H7, separated from a common ancestor of E. coli K-12 as long as 4.5 million years ago. Phylogenetic analysis reveals that old lineages of E. coli have acquired the same virulence factors in parallel, including a pathogenicity island involved in intestinal adhesion, a plasmid-borne haemolysin, and phage-encoded Shiga toxins. Such parallel evolution indicates that natural selection has favoured an ordered acquisition of genes and the progressive build-up of molecular mechanisms that increase virulence.

Molecular evolution and mosaic structure of alpha, beta, and gamma intimins of pathogenic Escherichia coli.

Two types of pathogenic Escherichia coli, enteropathogenic E. coli (EPEC) and enterohemorrhagic E. coli (EHEC), cause diarrheal disease by disrupting the intestinal environment through the intimate attachment of the bacteria to the intestinal epithelium. This process is mediated by intimin, an outer membrane protein that is homologous to the invasins of pathogenic Yersinia. The intimin (eae) gene is part of a pathogenicity island, a 35-kb segment of DNA that has been acquired independently in different groups of pathogens. Nucleotide sequences of eae of three EPEC and four EHEC strains representing distinct clonal lineages revealed an exceptionally high level of divergence (15%) in the amino acid sequences of alpha, beta, and gamma intimin molecules, most of which is concentrated in the C-terminal region. The gamma intimin sequences from E. coli strains with serotypes O157:H7, O55:H7, and O157:H- are virtually identical, supporting the hypothesis that these bacteria belong to a single clonal lineage. Sequences of beta intimin of EPEC strains of serotypes O111:H2 and O128:H2 show substantial differences from alpha and gamma intimins, indicating that these strains have evolved independently. Strong nonrandom clustering of polymorphic sites indicates that the intimin genes are mosaics, suggesting that protein divergence has been accelerated by recombination and diversifying selection.

Genetic diversity and population structure of Vibrio cholerae.

Multilocus enzyme electrophoresis (MLEE) of 397 Vibrio cholerae isolates, including 143 serogroup reference strains and 244 strains from Mexico and Guatemala, identified 279 electrophoretic types (ETs) distributed in two major divisions (I and II). Linkage disequilibrium was demonstrated in both divisions and in subdivision Ic of division I but not in subdivision Ia, which includes 76% of the ETs. Despite this evidence of relatively frequent recombination, clonal lineages may persist for periods of time measured in at least decades. In addition to the pandemic clones of serogroups O1 and O139, which form a tight cluster of four ETs in subdivision Ia, MLEE analysis identified numerous apparent clonal lineages of non-O1 strains with intercontinental distributions. A clone of serogroup O37 that demonstrated epidemic potential in the 1960s is closely related to the pandemic O1/O139 clones, but the nontoxigenic O1 Inaba El Tor reference strain is not. A strain of serogroup O22, which has been identified as the most likely donor of exogenous rfb region DNA to the O1 progenitor of the O139 clone, is distantly related to the O1/O139 clones. The close evolutionary relationships of the O1, O139, and O37 epidemic clones indicates that new cholera clones are likely to arise by the modification of a lineage that is already epidemic or is closely related to such a clone.

Sequence diversity of flagellin (fliC) alleles in pathogenic Escherichia coli.

To study the molecular evolution of flagellin, the protein subunit specifying flagellar (H) antigens, the fliC genes from 15 pathogenic strains of Escherichia coli were amplified by PCR and sequenced. Comparison of fliC sequences of H6 and H7 strains revealed that alleles have a mosaic structure indicating the occurrence of past horizontal transfer of DNA segments between strains. The close similarity of H7 sequences also indicates the exchange of an entire fliC H7 allele between distant clonal lineages. In addition, the ratio of silent substitutions to amino acid replacements suggests that a short segment in the central region of fliC has been under positive selection in the divergence of H6 and H7 alleles. Phylogenetic analysis demonstrates that the fliC sequences of O157:H7 and O55:H7 serotypes are nearly identical and highly divergent from those of E. coli strains expressing H6 and H2 flagellar antigens. A nonmotile clone of sorbitol-fermenting O157 has rapidly accumulated multiple mutations in fliC, presumably as a result of the silencing of flagellin expression.

Mutators and long-term molecular evolution of pathogenic Escherichia coli O157:H7.

It has been proposed that an increased mutation rate (indicated by the frequency of hypermutable isolates) has facilitated the emergence of Escherichia coli O157:H7. Analysis of the divergence of 12 genes shows no evidence that the pathogen has undergone an unusually high rate of mutation and molecular evolution.

Size and sequence polymorphism in the isocitrate dehydrogenase kinase/phosphatase gene (aceK) and flanking regions in Salmonella enterica and Escherichia coli.

The sequence of aceK, which codes for the regulatory catalytic enzyme isocitrate dehydrogenase kinase/phosphatase (IDH K/P), and sequences of the 5' flanking region and part or all of the 3' flanking region were determined for 32 strains of Salmonella enterica and Escherichia coli. In E. coli, the aceK gene was 1734 bp long in 13 strains, but in three strains it was 12 bp shorter and the stop codon was TAA rather than TGA. Strains with the shorter aceK lacked an open reading frame (f728) downstream between aceK and iclR that was present, in variable length, in the other strains. Among the 72 ECOR strains, the truncated aceK gene was present in all isolates of the B2 group and half of those of the D group. Other variant conditions included the presence of IS1 elements in two strains and large deletions in two strains. The aceK-aceA intergenic region varied in length from 48 to 280 bp in E. coli, depending largely on the number of repetitive extragenic palindromic (REP) sequences present. Among the ECOR strains, the number of REP elements showed a high degree of phylogenetic association, and sequencing of the region in the ECOR strains permitted partial reconstruction of its evolutionary history. In S. entica, the normal length of aceK was 1752 bp, but three other length variants, ranging from 1746 to 1785 bp, were represented in five of the 16 strains examined. The flanking intergenic regions showed relatively minor variation in length and sequence. The occurrence of several nonrandom patterns of distribution of polymorphic synonymous nucleotide sites indicated that intragenic recombination of horizontally exchanged DNA has contributed to the generation of allelic diversity at the aceK locus in both species.

Evolutionary genetics of the isocitrate dehydrogenase gene (icd) in Escherichia coli and Salmonella enterica.

Sequences of the icd gene, encoding isocitrate dehydrogenase (IDH), were obtained for 33 strains representing the major phylogenetic lineages of Escherichia coli and Salmonella enterica. Evolutionary relationships of the strains based on variation in icd are generally similar to those previously obtained for several other housekeeping and for invasion genes, but the sequences of S. enterica subspecies V strains are unusual in being almost intermediate between those of the other S. enterica subspecies and E. coli. For S. enterica, the ratio of synonymous (silent) to nonsynonymous (replacement) nucleotide substitutions between pairs of strains was larger than comparable values for 12 other housekeeping and invasion genes, reflecting unusually strong purifying selection against amino acid replacement in the IDH enzyme. All amino acids involved in the catalytic activity and conformational changes of IDH are strictly conserved within and between species. In E. coli, the level of variation at the 3' end of the gene is elevated by the presence in some strains of a 165-bp replacement sequence supplied by the integration of either lambdoid phage 21 or defective prophage element e14. The 72 members of the E. coli Reference Collection (ECOR) and five additional E. coli strains were surveyed for the presence of phage 21 (as prophage) by PCR amplification of a phage 21-specific fragment in and adjacent to the host icd, and the sequence of the phage 21 segment extending from the 3' end of icd through the integrase gene (int) was determined in nine strains of E. coli. Phage 21 was found in 39% of E. coli strains, and its distribution among the ECOR strains is nonrandom. In two ECOR strains, the phage 21 int gene is interrupted by a 1,313-bp insertion element that has 99.3% nucleotide sequence identity with IS3411 of E. coli. The phylogenetic relationships of phage 21 strains derived from sequences of two different genomic regions were strongly incongruent, providing evidence of frequent recombination.

Evolutionary genetics of Pasteurella haemolytica isolates recovered from cattle and sheep.

Genetic diversity and relationships among 194 Pasteurella haemolytica isolates, which were recovered predominantly from cattle (39%) and sheep (58%) suffering from pneumonic pasteurellosis in the United Kingdom, Germany, and the United States, were estimated by examination of allelic variation at 18 enzyme-encoding loci detected by multilocus enzyme electrophoresis. The isolates formed two major divisions. One included 178 Pasteurella haemolytica sensu stricto strains representing serotypes A1, A2, A5 to A9, A12 to A14, and A16; the other was composed of 16 isolates belonging to the A11 taxon. P. haemolytica isolates were classified into 22 electrophoretic types (ETs) that formed three primary phylogenetic lineages. One lineage was represented by ovine serotype A2 isolates, a second lineage consisted of bovine serotype A2, together with serotype A7 and A13 isolates, and the third lineage included isolates representing all of the other serotypes, as well as a second group of serotype A7 strains. Electrophoretic types were nonrandomly associated with specific capsular serotypes, lipopolysaccharide (LPS) types, outer membrane protein (OMP) types, and host species. Bovine isolates were represented by only three serotypes (A1, A2, and A6) in 5 ETs, whereas ovine isolates were represented by all of the serotypes in 19 ETs. The majority (76%) of bovine isolates were of serotypes A1 or A6 and belonged to a single ET that marked a virulent, cattle-specific clonal group. Among the ovine isolates, 40% were of serotype A2 and belonged to two ETs that represented two virulent, sheep-specific clonal groups. Bovine A1 and A6 isolates and bovine A2 isolates were phylogenetically distinct from ovine isolates of the same serotypes, indicating that different subpopulations of these serotypes are associated with disease in cattle and sheep. Consistent differences in the OMP profiles of strains of the bovine and ovine lineages of these three serotypes suggest that certain OMPs are involved in host specificity and virulence. Evolutionary relationships among P. haemolytica isolates indicate that the ancestral host is the sheep and that several distinct clonal lineages have crossed the species barrier into cattle. The A11 taxon is a heterogeneous group of opportunistic pathogens of sheep that represents a separate species.

Sex-related coping responses in mice selectively bred for aggression.

Sex differences in strategies of coping with novel situations were studied in three strains of mice with regard to metabolism and open-field and maze activity as well as learning-induced adjustment. The 140 mice were selectively bred for high (Turku Aggressive [TA]) and low (Turku Nonaggressive [TNA]) levels of aggressiveness and originated from a Swiss albino stock normally distributed [N] for aggressiveness. The results indicated that TNA sex differences are more similar to those of the control N mice as compared to those of TA mice. In maze learning, however, the sex differences of TA mice are more in agreement with those of the N strain. Recordings of metabolism and open-field as well as maze activity were correlates of both gender and strain. Sex differences in learning-induced open-field coping behavior were unrelated to strain.

Comparative genetics of the inv-spa invasion gene complex of Salmonella enterica.

The chromosomal region containing the Salmonella enterica pathogenic island inv-spa was present in the last common ancestor of all the contemporary lineages of salmonellae. For multiple strains of S. enterica, representing all eight subspecies, nucleotide sequences were obtained for five genes of the inv-spa invasion complex, invH, invE, invA, spaM, and spaN, al of which encode proteins that are required for entry of the bacteria into cultured epithelial cells. The invE, invA, spaM, and spaN genes were present in all eight subspecies of S. enterica, and for invE and invA and their products, levels of sequence variation among strains were within the ranges reported for housekeeping genes. In contrast, the InvH, SpaM, and SpaN proteins were unusually variable in amino acid sequence. Furthermore, invH was absent from the subspecies V isolates examined. The SpaM and SpaN proteins provide further evidence of a relationship (first detected by Li et al. [J. Li, H. Ochman, E. A. Groisman, E. F. Boyd, F. Solomon, K. Nelson, and R. K. Selander, Proc. Natl. Acad. Sci. USA 92:7252-7256, 1995]) between the cellular location of the products of the inv-spa genes and evolutionary rate, as reflected in the level of polymorphism within S. enterica. Invasion proteins that are membrane bound or membrane associated are relatively conserved in amino acid sequence, whereas those that are exported to the extracellular environment are hypervariable, possibly reflecting the action of diversifying selection.

Molecular population genetics of the alcohol dehydrogenase locus in the Hawaiian drosophilid D. mimica.

Sequence variation among 10 alleles of the alcohol dehydrogenase (Adh) gene of the Hawaiian drosophilid D. mimica was analyzed with reference to the evolutionary history of the Hawaiian subgroup as well as to levels and patterns of polymorphism of the Adh gene in continental drosophilid species. The Adh gene of D. mimica is less polymorphic than that of other drosophilid species, and no replacement substitutions were found. Statistical analyses of the Adh alleles suggested the action of balancing selection and revealed significant linkage disequilibrium among three of the variable sites. The effective population size was estimated to be only slightly smaller than that of continental species and, surprisingly, on the same order of magnitude as the actual size.

Molecular genetic relationships of the salmonellae.

A multilocus enzyme electrophoresis analysis of 96 strains of the salmonellae distinguished 80 electrophoretic types (ETs) and placed them in eight groups, seven of which correspond precisely to the seven taxonomic groups (I, II, IIIa, IIIb, IV, V, and VI) previously defined on the basis of biotype and genomic DNA hybridization. In addition, multilocus enzyme electrophoresis identified an eighth distinctive group (designated VII) composed of five strains that had been assigned to group IV on the basis of biotype. An analysis of variation in the combined nucleotide sequences of five housekeeping genes among 16 strains representing all eight groups yielded estimates of overall genetic relationships that are fully consistent with those indicated by DNA hybridization. However, the nucleotide sequences of seven invasion genes (inv/spa) in the strains of group VII were closely similar to those of strains of group IV. These findings are interpreted as evidence that group VII represents an old, differentiated lineage to which one or more large parts of the chromosomal genome of the group IV lineage, including the 40-kb segment on which the invasion genes are located, have been horizontally transferred. All lines of molecular genetic evidence indicate that group V is very strongly differentiated from all other groups, thus supporting its current taxonomic treatment as a species, Salmonella bongori, separate from S. enterica. The Salmonella Reference Collection C, composed of the 16 strains used in DNA sequence studies, has been established for research on variation in natural populations.

Characterization of Streptococcus agalactiae strains by multilocus enzyme genotype and serotype: identification of multiple virulent clone families that cause invasive neonatal disease.

The chromosomal genotypes of 277 isolates of 16 serotypes of Streptococcus agalactiae were characterized by analysis of electrophoretically demonstrable allele profiles at 12 metabolic enzyme loci. The collection comprised the type strain and 276 strains recovered from French symptomatic and asymptomatic subjects. Sixty-one distinctive electrophoretic types (ETs), representing multilocus clonal genotypes, were identified. Cluster analysis of the ETs revealed two primary phylogenetic divisions separated by a genetic distance of 0.62, Division I contained 67 isolates which could be assigned to 13 ETs. Twenty-seven of these isolates were from samples of cerebrospinal fluid (CSF) from neonatal meningitis patients. Two ETs, separated by a genetic distance of 0.217, contained 26 of these 27 isolates. Division II contained 210 isolates, of which 27 were isolated from CSF. This division was more polymorphic and included 48 ETs. Spanning a genetic distance of 0.3, three clusters and one ET were identified within this group. Twenty-four of 27 strains isolated from CSF belonged to one cluster, and 19 of them belonged to two adjacent ETs with a genetic distance of 0.083. Fifty-five of the 68 serotype Ia strains and 24 of the 26 serotype Ib strains were each confined to one of the evolutionary lineages, and 85 of the 86 strains which carried protein antigen c belonged to phylogenetic division II. Most of the type III organisms were assigned to two clone families. The characteristics of this French population argue for the existence of particular groups of strains responsible for neonatal meningitis and demonstrate that serotyping can supply information about the genetic distribution of strains.

Relationship between evolutionary rate and cellular location among the Inv/Spa invasion proteins of Salmonella enterica.

For 21 strains of Salmonella enterica, nucleotide sequences were obtained for three invasion genes, spaO, spaP, and spaQ, of the chromosomal inv/spa complex, the products of which form a protein export system required for entry of the bacteria into nonphagocytic host cells. These genes are present in all eight subspecies of the salmonellae, and homologues occur in a variety of other bacteria, including the enteric pathogens Shigella and Yersinia, in which they are plasmid borne. Evolutionary diversification of the invasion genes among the subspecies of S. enterica has been generally similar in pattern and average rate to that of housekeeping genes. However, the range of variation in evolutionary rate among the invasion genes is unusually large, and there is a relationship between the evolutionary rate and cellular location of the invasion proteins, possibly reflecting diversifying selection on exported proteins in adaptation to variable host factors in extracellular environments. The SpaO protein, which is hypervariable in S. enterica and exhibits only 24% sequence identity with its homologues in Shigella and Yersinia, is secreted. In contrast, the membrane-associated proteins SpaP, SpaQ, and InvA are weakly polymorphic and have > 60% sequence identity with the corresponding proteins of other enteric bacteria. Acquisition of the inv/spa genes may have been a key event in the evolution of the salmonellae as pathogens, following which the invention of flagellar phase shifting facilitated niche expansion to include warm-blooded vertebrates.

Species limits in Rhizobium populations that nodulate the common bean (Phaseolus vulgaris).

Evolutionary genetic relationships among 146 bean-nodulating Rhizobium strains, including 94 field isolates from three localities in Colombia and 36 strains from Mexico, were examined by multilocus enzyme electrophoresis and restriction fragment length polymorphism analysis of a PCR-amplified 260-bp segment of the 16S rRNA gene. Seventy-five electrophoretic types (ETs), corresponding to multilocus enzyme genotypes, were identified, including a genotypically diverse group of 18 ETs in Colombia that is strongly differentiated from the ETs of R. etli, which occur in Mexico, Colombia, and Brazil. Most strains of the distinctive Colombian ETs carried the same 16S rRNA allele as did strains of R. etli, but, surprisingly, 17 isolates of two of these ETs had the allele that is characteristic of R. leguminosarum, and strains of two other divergent groups of ETs were also polymorphic for the two alleles. No fully satisfactory explanation for the occurrence of the R. leguminosarum 16S rRNA allele in three distantly related groups of strains is available, but horizontal transfer and recombination of the gene, in whole or in part, would seem to be more plausible than convergence in nucleotide sequence.

Intergeneric transfer and recombination of the 6-phosphogluconate dehydrogenase gene (gnd) in enteric bacteria.

The gnd gene, encoding 6-phosphogluconate dehydrogenase (EC 1.1.1.44), was sequenced in 87 strains of 15 species assigned to five nominal genera of the Enterobacteriaceae, including 36 isolates of Salmonella enterica and 32 strains of Escherichia coli. In S. enterica, the effective (realized) rate of recombination of horizontally transferred gnd sequences is only moderately higher than the rates for other chromosomal housekeeping genes. In contrast, recombination at gnd has occurred with such high frequency in Escherichia coli that the indicated evolutionary relationships among strains are not congruent with those estimated by sequence analysis of other genes and by multilocus enzyme electrophoresis. E. coli and S. enterica apparently have not exchanged gnd sequences, but those of several strains of E. coli have been imported from species of Citrobacter and Klebsiella. The relatively frequent exchange of gnd within and among taxonomic groups of the Enterobacteriaceae, compared with other housekeeping genes, apparently results from its close linkage with genes that are subject to diversifying selection, including those of the rfb region determining the structure of the O antigen polysaccharide.

Open-field thigmotaxis during various phases of the reproductive cycle.

Three experiments were conducted on female Swiss albino mice in order to determine whether open-field (OF) ambulation and thigmotaxis (wall-seeking) alter both as functions of reproductive phases such as estrus, mating, gestation, partus, lactation, weaning, as well as of presence, absence and OF location of a newborn pup. The results revealed increased thigmotaxis as an effect of estrus (p < 0.001). From mating to one month postpartum a growing tendency towards maternal occupation of centre OF partitions was observed. Maternal ambulation reached its peak two days before partus but dropped to its lowest level during the lactation period at which time the mice were the least wall-seeking. Ambulation was enhanced (p < 0.001) in the presence of a pup regardless of OF starting point and the pup's placement. However, when the mouse dam was initially placed together with the pup by the OF wall, she hardly penetrated centre units at all.

Recombinational basis of serovar diversity in Salmonella enterica.

The fliC gene, which encodes phase 1 flagellin, was sequenced in strains of 15 Salmonella enterica serovars expressing flagellar antigenic factors of the g series. The occurrence of each of the flagellin serotypes g,m, m,t, and g,z51 in distantly related strains is the result of horizontal exchange of DNA, as indicated by identity or close similarity in nucleotide sequence of all or parts of the antigenic factor-determining central region of fliC. The flagellin genes of some serovars are complex mosaic structures composed of diverse segments derived through multiple recombination events. Thus, recombination of horizontally transferred segments (intragenic) or entire genes (assortative) within and among subspecies is identified as a major evolutionary mechanism generating both allelic variation at the fliC locus and serovar diversity in natural populations. Evidence that flagellar serological diversity is promoted by diversifying selection in adaptation to host immune defense system or flagellotropic phage is discussed.