Phylogeny and Strain Typing of Escherichia coli, Inferred from Variation at Mononucleotide Repeat Loci

Multilocus sequencing of housekeeping genes has been used previously for bacterial strain typing and for inferring evolutionary relationships among strains of Escherichia coli. In this study, we used shorter intergenic sequences that contained simple sequence repeats (SSRs) of repeating mononucleotide motifs (mononucleotide repeats [MNRs]) to infer the phylogeny of pathogenic and commensal E. coli strains. Seven noncoding loci (four MNRs and three non-SSRs) were sequenced in 27 strains, including enterohemorrhagic (six isolates of O157:H7), enteropathogenic, enterotoxigenic, B, and K-12 strains. The four MNRs were also sequenced in 20 representative strains of the E. coli reference (ECOR) collection. Sequence polymorphism was significantly higher at the MNR loci, including the flanking sequences, indicating a higher mutation rate in the sequences flanking the MNR tracts. The four MNR loci were amplifiable by PCR in the standard ECOR A, B1, and D groups, but only one (yaiN) in the B2 group was amplified, which is consistent with previous studies that suggested that B2 is the most ancient group. High sequence compatibility was found between the four MNR loci, indicating that they are in the same clonal frame. The phylogenetic trees that were constructed from the sequence data were in good agreement with those of previous studies that used multilocus enzyme electrophoresis. The results demonstrate that MNR loci are useful for inferring phylogenetic relationships and provide much higher sequence variation than housekeeping genes. Therefore, the use of MNR loci for multilocus sequence typing should prove efficient for clinical diagnostics, epidemiology, and evolutionary study of bacteria.

Related papers:

title: Standard reference strains of Escherichia coli from natural populations
abstract: A set of 72 reference strains of Escherichia coli isolated from a variety of hosts and geographical locations has been established for use in studies of variation and genetic structure in natural populations. The strains, which have been characterized by multilocus enzyme electrophoresis, are representative of the range of genotypic variation in the species as a whole.

title: Clonal relationships among Escherichia coli strains that cause hemorrhagic colitis and infantile diarrhea
abstract: The genetic relationships among 1,300 isolates of Escherichia coli representing 16 serotypes associated with enteric disease, including O157:H7 strains recovered from patients with hemorrhagic colitis and hemolytic uremic syndrome and O26:H11, H7, O111:H2, and O128:H2 strains, many of which were isolated originally from infants with diarrhea, were estimated from allelic variation among 20 enzyme-encoding genes detected by multilocus enzyme electrophoresis. Multiple electrophoretic types were observed among isolates of each serotype, with isolates of the same O serogroup differing on average at 28% of the enzyme loci. Comparisons of the multilocus enzyme profiles revealed that 72% of the isolates belong to 15 major electrophoretic types, each of which corresponds to a bacterial clone with a wide geographic distribution. Genetically, the O157:H7 clone is most closely related to a clone of O55:H7 strains that has long been associated with worldwide outbreaks of infantile diarrhea. We propose that the new pathogen emerged when an O55:H7-like progenitor, already possessing a mechanism for adherence to intestinal cells, acquired secondary virulence factors (Shiga-like cytotoxins and plasmid-encoded adhesins) via horizontal transfer and recombination.

title: Escherichia coli that cause diarrhea: enterotoxigenic, enteropathogenic, enteroinvasive, enterohemorrhagic, and enteroadherent.
abstract: There are four major categories of diarrheagenic Escherichia coli: enterotoxigenic (a major cause of travelers' diarrhea and infant diarrhea in less-developed countries), enteroinvasive (a cause of dysentery), enteropathogenic (an important cause of infant diarrhea), and enterohemorrhagic (a cause of hemorrhagic colitis and hemolytic uremic syndrome). Besides manifesting distinct clinical patterns, these categories of E. coli differ in their epidemiology and pathogenesis and in their O:H serotypes. Common features (albeit distinct for each category) include plasmid-encoded virulence properties, characteristic interactions with intestinal mucosa, and elaboration of various types of enterotoxins or cytotoxins. A less-well-defined fifth category of diarrheagenic E. coli is that of enteroadherent E. coli, so far identifiable only by their pattern of adherence to Hep-2 cells in tissue culture.

title: Enzyme polymorphism and genetic population structure in Escherichia coli and Shigella.
abstract: Electrophoretically demonstrable variation in 12 enzymes was studied in more than 1 600 isolates of Escherichia coli from human and animal sources and in 123 strains of the four species of Shigella. All 12 enzymes were polymorphic; and the number of allozymes (mobility variants), which were equated with alleles, averaged 9.3 per locus in E. coli. For Shigella species, the mean number of alleles was 2.9 per locus. Some 77% of the allozymes recorded in Shigella were shared with E. coli. A total of 302 unique genotypic combinations of alleles over the 12 loci (electrophoretic types, ETs) was distinguished, of which 279 represented E. coli and 23 were Shigella. Among electrophoretic types, mean allelic diversity per locus was 0.52 for E. coli and 0.29 for Shigella. It was estimated that there are, on the average, about 0.3 detectable codon differences per locus between pairs of strains of E. coli and Shigella, which is roughly equivalent to 1.2 amino acid differences per enzyme. Evidence that the enzyme loci studied are a random sample of the genome is provided by a significant positive correlation between estimates of genetic divergence between pairs of strains obtained by DNA reassociation tests and estimates of genetic distance between the same strains based on electrophoresis. A principal components analysis of allozyme profiles revealed that the 302 ETs fall into three overlapping clusters, reflecting strong non-random associations of alleles, largely at four loci. Each of the four ETs of E. coli that have been most frequently recovered from natural populations has an allozyme profile that is very similar to, or identical with, the hypothetical modal ET of one of the groups. ETs of Shigella fall into two of the groups. No biological significance can at present bbe attributed to the genetic structure revealed by Multilocus electrophoretic techniques. The electrophoretic data are fully compatible with other molecular and more conventional evidence of a close affinity between E. coli and Shigella, and they raise questions regarding the present assignments of certain strains to species. In support of evidence from DNA reassociation tests and serotyping, the present study suggests that S. sonnei is homogeneous in chromosomal genotype.

title: Evolutionary relationships among pathogenic and nonpathogenic Escherichia coli strains inferred from multilocus enzyme electrophoresis and mdh sequence studies
abstract: Within the species Escherichia coli, there are commensal strains and a variety of pathogenic strains, including enteropathogenic E. coli (EPEC), enterohemorrhagic E. coli (EHEC), enterotoxigenic E. coli (ETEC), enteroinvasive E. coli (EIEC), and urinary tract infection (UTI) strains. The pathogenic strains are identified by serotype and by possession of specific virulence determinants (toxins and adhesions, etc.) encoded by either monocistronic genes, plasmids, or pathogenicity islands. Although there are studies on the relationships between selected pathogenic strains, the relatedness among the majority of the pathogenic forms to each other, to commensal E. coli, and to the genus Shigella (which has often been suggested to be part of E. coli) has not been determined. We used multilocus enzyme electrophoresis (MLEE) at 10 enzyme loci and the sequence of the mdh housekeeping gene to study the genetic relationships of pathogenic E. coli strains (including Shigella clones), namely, 5 EPEC strains (serotypes O111 and O55), 3 EHEC strains (serotype O157), 6 ETEC strains (serotypes O78, O159, and O148), 5 EIEC strains (serotypes O124, O28, and O112), and 13 Shigella strains representing clones Flexneri, Dysenteriae, Boydii, and Sonnei, to commensal E. coli strains. Both the MLEE and mdh sequence trees reveal that EPEC, EHEC, ETEC, EIEC, and UTI strains are distributed among the ECOR set groups, with no overall clustering of EPEC, ETEC, EIEC, or UTI strains. The genus Shigella is shown to comprise a group of closely related pathogenic E. coli strains. Six pathogenic strains, i.e., M502 (EIEC; O112ac:NM), M503 (EPEC; O111:H12), M526 (ETEC; O159:H4), M522 (EPEC; O111ac:H12), M524 (ETEC; O78:H11), and M506 (ETEC; O78:H11), were found to have mdh sequences identical to those of five ECOR group A strains (ECOR5, ECOR10, ECOR14, ECOR6, and K-12). All 11 strains are closely related by MLEE. The results indicate that pathogenic strains of E. coli do not have a single evolutionary origin within E. coli but have arisen many times. The results also suggest the possibility that any E. coli strain acquiring the appropriate virulence factors may give rise to a pathogenic form.

title: Methods of multilocus enzyme electrophoresis for bacterial population genetics and systematics
abstract: Methodes d'extraction d'enzymes, d'electrophorese en gel et de coloration specifique des enzymes utilisees pour etudier la variation genetique chez Escherichia coli et d'autres bacteries. Les procedures decrites peuvent s'appliquer, avec quelques modifications mineures, a toute espece de bacterie

title: Escherichia coli molecular phylogeny using the incongruence length difference test.
abstract: Molecular phylogeny of the species Escherichia coli using the E. coli reference (ECOR) collection strains has been hampered by (1) the absence of rooting in the commonly used phenogram obtained from multilocus enzyme electrophoresis (MLEE) data and (2) the existence of recombination events between strains that scramble phylogenetic trees reconstructed from the nucleotide sequences of genes. We attempted to determine the phylogeny for E. coli based on the ECOR strain data by extracting from GenBank the nucleotide sequences of 11 chromosomal structural and 2 plasmid genes for which the Salmonella enterica homologous gene sequences were available. For each of the 13 DNA data sets studied, incongruence with a nonnucleotide whole-genome data set including MLEE, random amplified polymorphic DNA, and rrn restriction fragment length polymorphism data was measured using the incongruence length difference (ILD) test of Farris et al. As previously reported, the incongruence observed between the gnd and plasmid gene data and the whole-genome data was multiple, indicating numerous horizontal transfer and/or recombination events. In five cases, the incongruence detected by the ILD test was punctual, and the donor group was identified. Congruence was not rejected for the remaining data sets. The strains responsible for incongruences with the whole-genome data set were removed, leading to a "prior-agreement" approach, i.e., the determination of a phylogeny for E. coli based on several genes, excluding (1) the genes with multiple incongruences with the whole genome data, (2) the strains responsible for punctual incongruences, and (3) the genes incongruent with each other. The obtained phylogeny shows that the most basal group of E. coli strains is the B2 group rather than the A group, as generally thought. The D group then emerges as the sister group of the rest. Finally, the A and B1 groups are sister groups. Interestingly, the most primitive taxon within E. coli in terms of branching pattern, i.e., the B2 group, includes highly virulent extraintestinal strains with derived characters (extraintestinal virulence determinants) occurring on its own branch.

title: Potential genetic functions of tandem repeated DNA sequence blocks in the human genome are based on a highly conserved “chromatin folding code”
abstract: SummaryThis review is based on a thorough description of the structure and sequence organization of tandemly organized repetitive DNA sequence families in the human genome; it is aimed at revealing the locusspecific sequence organization of tandemly repetitive sequence structures as a highly conserved DNA sequence code. These repetitive so-called “super-structures” or “higher-order” structures are able to attract specific nuclear proteins. I shall define this code therefore as a “chromatin folding code”. Since locus-specific superstructures of tandemly repetitive sequence units are present not only in the chromosome centromere or telomere region but also on the arms of the chromosomes, I assume that their chromatin folding code may contribute to, or even organize, the folding pathway of the chromatin chain in the nucleus. The “chromatin folding code” is based on its specific “chromatin code”, which describes the sequence dependence of the helical pathway of the DNA primary sequence (i.e., secondary structure) entrapping the histone octamers in preferential positions. There is no periodicity in the distribution of the nucleosomes along the DNA chain. The folding pathway of the nucleosomal chromatin chain is however still flexible and determined by e.g., the length of the DNA chain between the nucleosomes. The fixation and stabilization of the chromatin chain in the space of the nucleus (i.e., its “functional state”) may be mediated by additionally unique DNA protein interactions that are dictated by the “chromatin folding code”. The unique DNA-protein interactions around the centromeres of human chromosomes are revealed for example by their “C-banding”. I wish to stress that it is not my aim to relate each block of repetitive DNA sequences to a specific “chromatin folding code”, but I shall demonstrate that there is an inherent potential for tandem repeated sequence units to develop a locus-specific repetitive higher order structure; this potential may create a specific chromatin folding code whenever a selection force exists at the position of this repetitive DNA structure in the genome.

title: Hypervariability of simple sequences as a general source for polymorphic DNA markers.
abstract: Short simple sequence stretches occur as highly repetitive elements in all eukaryotic genomes and partially also in prokaryotes and eubacteria. They are thought to arise by slippage like events working on randomly occurring internally repetitive sequence stretches. This predicts that they should be generally hypervariable in length. I have used the polymerase chain reaction (PCR) process to show that several randomly chosen simple sequence loci with different nucleotide composition and from different species show extensive length polymorphisms. These simple sequence length polymorphisms (SSLP) may be usefully exploited for identity testing, population studies, linkage analysis and genome mapping.

title: Mutation and migration in models of microsatellite evolution
abstract: Dynamic and statistical properties of generalized stepwise mutation models are described and used to compare data on human di-, tri-, and tetranucleotide polymorphisms. The time-dependent behavior of an island model with stepwise mutation is analyzed and its equilibrium properties used to estimate the product Nm of the population size and the migration rate. Population statistics are derived from the complete equilibrium analysis of the model, and these are combined to give an estimate of Nm that is related to Slatkin’s 1995 estimate which used RST . When this new statistic is applied to a set of 85 human microsatellite polymorphisms, the resulting population clusters match the tree of Bowcock et al. 1994 quite well.  Because of their high level of polymorphism and ubiquity in eukaryotic genomes, microsatellites are widely preferred markers in evolutionary and ecological studies. Databases of microsatellites isolated for population-level analyses are under development. Statistical evaluation and evolutionary interpretation of microsatellite polymorphisms demand models for the origin and maintenance of this variability, and the models used so far are variants of the stepwise mutation model (SMM), which was originally introduced by Ohta and Kimura (1973) to model electrophoretically detectable enzyme variation in finite populations. The first part of this paper describes these models and illustrates how stepwise mutation, in combination with genetic drift, affects important measures of withinand betweenpopulation variation. Using empirical estimates of overall mutation rates, these measures may then provide estimates of average population size or divergence times between populations. In the second part of the paper, the dynamics of a set of populations subject to migration according to Wright’s (1943) island model are developed. From the equilibrium structure of this model, estimators of the extent of migration based on population statistics are suggested. We shall assume throughout that the microsatellites are perfect; that is, each allele at a locus is completely specified by the number of times a motif is repeated. Generations are non-overlapping and each is produced from the previous by multinomial sampling from the parental generation’s array of alleles. Our analysis follows that of Moran (1975), who developed recursions for the population central moments under the simplest stepwise mutation model permitting alleles of arbitrary (positive and negative) repeat number. The modifications to this model that we discuss below include asymmetric mutations of arbitrary size, a simple model of linear bias in which the rate of mutation depends on the number of repeats in the allele, and range constraints on the permitted repeat score. Moran’s  A Generalized SMM. Consider a population of N diploid individuals and a locus at which each allele is characterized by its repeat score, which may take any positive or negative integer values. Let μc be the probability that mutation changes an allele by c repeat units irrespective of its original count. The total mutation rate is then μ = ∑ c 6=0 μc and the expected change in repeat score for any allele is c̄ = ∑ cμc. The variance in repeat score change is σ m = ∑ c(μc/μ)− c̄. We write w = μσ m. The standard one-step symmetric model has μ+1 = μ−1 = μ/2, say, and μ2 = 0 otherwise, in which case c̄ = 0 and w = μ. We shall refer to this special case as the one-step symmetric SMM. Let pi be the frequency of the allele carrying i repeats in the parental generation. Then, following mutation the frequency of this allele is

title: Abundant microsatellite polymorphism in Saccharomyces cerevisiae, and the different distributions of microsatellites in eight prokaryotes and S. cerevisiae, result from strong mutation pressures and a variety of selective forces.
abstract: We examined the distributions of short tandemly repeated DNAs (microsatellites) in nine complete microbial genomes (Saccharomyces cerevisiae, Archaeoglobus fulgidus, Escherichia coli, Haemophilus influenzae, Helicobacter pylori, Methanococcus jannaschii, Mycoplasma pneumoniae, M. genitalium, and Synechocystis PCC6803.) These repeats contribute differently to the global features of these genomes, and we explore the evolutionary implications of these differences by empirical examination of length polymorphisms at 20 long triplet-repeats repeats in S. cerevisiae, and by comparison of observed and expected repeat distributions. All of a sample of 20 microsatellites found in S. cerevisiae are highly polymorphic in length, suggesting that mutation pressure overcomes overall selection for small genome size that will tend to shorten or eliminate unnecessary DNA. By comparison, prokaryotes have fewer long repeats than expected, except for a few statistically improbable repeats that appear to function in gene regulation. Finally, we find that in all these genomes there is an excess of repeats shorter than those traditionally considered to be microsatellites. This finding suggests that even in prokaryotes these repeats are being generated by mutational pressures. These results have important potential implications for understanding genome stability and evolution in these microbial species.

title: Simple sequence repeats in Escherichia coli: abundance, distribution, composition, and polymorphism.
abstract: Computer-based genome-wide screening of the DNA sequence of Escherichia coli strain K12 revealed tens of thousands of tandem simple sequence repeat (SSR) tracts, with motifs ranging from 1 to 6 nucleotides. SSRs were well distributed throughout the genome. Mononucleotide SSRs were over-represented in noncoding regions and under-represented in open reading frames (ORFs). Nucleotide composition of mono- and dinucleotide SSRs, both in ORFs and in noncoding regions, differed from that of the genomic region in which they occurred, with 93% of all mononucleotide SSRs proving to be of A or T. Computer-based analysis of the fine position of every SSR locus in the noncoding portion of the genome relative to downstream ORFs showed SSRs located in areas that could affect gene regulation. DNA sequences at 14 arbitrarily chosen SSR tracts were compared among E. coli strains. Polymorphisms of SSR copy number were observed at four of seven mononucleotide SSR tracts screened, with all polymorphisms occurring in noncoding regions. SSR polymorphism could prove important as a genome-wide source of variation, both for practical applications (including rapid detection, strain identification, and detection of loci affecting key phenotypes) and for evolutionary adaptation of microbes.

title: The microsatellites of Escherichia coli: rapidly evolving repetitive DNAs in a non‐pathogenic prokaryote
abstract: We have demonstrated hypervariability of native short‐motif repeats (microsatellites) in Escherichia coli. Twenty‐five of the longest microsatellites in the E. coli genome were identified. These were analysed for length variability among 22 wild‐type (non‐mutator) isolates from the E. coli collection of reference (ECOR). Non‐coding mononucleotide repeats are consistently polymorphic among these genetically diverse E. coli. Length differences in variable microsatellites allowed all E. coli strains examined to be uniquely differentiated. Phylogenetic analysis of the variable repeats shows ubiquitous homoplasy at the level of divergence represented by the sample set, suggesting that these markers are hypermutable and should prove valuable for the discrimination of closely related strains that are not otherwise genetically differentiable. Genomic analyses suggest that similar markers are also likely to be found in all other prokaryotes.

title: Phylogenetic distribution of branched RNA-linked multicopy single-stranded DNA among natural isolates of Escherichia coli
abstract: Multicopy single-stranded DNA (msDNA), a branched DNA-RNA molecule, has been shown in Escherichia coli B and clinical strain Cl-1 to be synthesized by reverse transcriptase. We report that 13% of the strains of the ECOR collection, a sample of 72 E. coli isolates representing the breadth of genetic variation of the species, produce msDNA. Three of the four major subspecific groups include msDNA-producing strains. Screening of 25 isolates that are genetically related to msDNA-producing clinical strains uncovered 22 additional msDNA-producing strains. A phylogenetic tree based on allelic variation detected electrophoretically at 20 enzyme-encoding loci revealed two major clusters and several deep branches composed of strains that synthesize msDNA. Although E. coli K-12 does not harbor msDNA, other closely related strains of the K-12 family do. The results support the hypothesis that msDNA-synthesizing systems, including reverse transcriptase genes, were acquired recently and independently in different lineages of E. coli.

title: Parallel evolution of virulence in pathogenic Escherichia coli
abstract: 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 × 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.

title: A program for calculating and displaying compatibility matrices as an aid in determining reticulate evolution in molecular sequences
abstract: Reticulate evolution in molecular sequences is caused by recombination or gene conversion, and may interfere with the reconstruction of evolutionary history. This paper presents a program which calculates compatibility matrices for detecting reticulate evolution. In addition to visual inspection of matrices, they can be analysed statistically for clustering. The method is demonstrated using human and chimpanzee gamma-globin sequences.

title: How clonal are bacteria?
abstract: Data from multilocus enzyme electrophoresis of bacterial populations were analyzed using a statistical test designed to detect associations between genes at different loci. Some species (e.g., Salmonella) were found to be clonal at all levels of analysis. At the other extreme, Neisseria gonorrhoeae is panmictic, with random association between loci. Two intermediate types of population structure were also found. Neisseria meningitidis displays what we have called an "epidemic" structure. There is significant association between loci, but this arises only because of the recent, explosive, increase in particular electrophoretic types; when this effect is eliminated the population is found to be effectively panmictic. In contrast, linkage disequilibrium in a population of Rhizobium meliloti exists because the sample consisted of two genetically isolated divisions, often fixed for different alleles: within each division association between loci was almost random. The method of analysis is appropriate whenever there is doubt about the extent of genetic recombination between members of a population. To illustrate this we analyzed data on protozoan parasites and again found panmictic, epidemic, and clonal population structures.

title: Molecular Evolution and Phylogenetics
abstract: 1. Molecular basis of evolution 2. Evolutionary changes of amino acid sequences 3. Evolutionary changes of DNA sequences 4. Synonymous and nonsynonymous nucleotide substitutions 5. Phylogenetic trees 6. Phylogenetic inference: Distance methods 7. Phylogenetic inference: Maximum parsimony methods 8. Phylogenetic inference: Maximum likelihood methods 9. Accuracies and statistical tests of phylogenetic trees 10. Molecular clocks and linearized trees 11. Ancestral nucleotide and amino acid sequences 12. Genetic polymorphism and evolution 13. Population trees from genetic markers 14. Perspectives Appendices A. Mathematical sumbols and notations B. Geological timescale C. Geological events in the Cenozoic and Meszoic eras D. Evolution of organisms based on the fossil record

title: Clonal divergence in Escherichia coli as a result of recombination, not mutation.
abstract: Nucleotide sequence analysis was performed on 12 natural isolates of Escherichia coli in four loci located in close proximity on the chromosome. A comparison of gene genealogies indicated that three recombination events have occurred in a subset of the strains (ECOR group A) in the time since their divergence from a common ancestor, while during the same time, no mutational divergence has occurred. The common ancestor of this subset existed no more than 2400 years ago, and recombination was shown to occur at a rate of 5.0 x 10(-9) changes per nucleotide per generation--50-fold higher than the mutation rate. Thus, recombination has been the dominant force driving the clonal divergence of the ECOR group A strains and must be considered a significant factor in structuring E. coli populations.

title: Molecular archaeology of the Escherichia coli genome.
abstract: The availability of the complete sequence of Escherichia coli strain MG1655 provides the first opportunity to assess the overall impact of horizontal genetic transfer on the evolution of bacterial genomes. We found that 755 of 4,288 ORFs (547.8 kb) have been introduced into the E. coli genome in at least 234 lateral transfer events since this species diverged from the Salmonella lineage 100 million years (Myr) ago. The average age of introduced genes was 14.4 Myr, yielding a rate of transfer 16 kb/Myr/lineage since divergence. Although most of the acquired genes subsequently were deleted, the sequences that have persisted ( approximately 18% of the current chromosome) have conferred properties permitting E. coli to explore otherwise unreachable ecological niches.

title: Linkage disequilibrium due to random genetic drift in finite subdivided populations.
abstract: In order to clarify the mechanisms responsible for the observed linkage disequilibrium, such as that found between markers in the major histocompatibility complex of man and mouse, linkage disequilibrium between two linked loci was studied for a finite population with a subdivided structure. The infinite allele model was used. In analogy with the subdivision of the inbreeding coefficient, the linkage disequilibrium coefficient was subdivided, and various variance components of disequilibrium were defined. It was found that the disequilibrium components may get very large when migration is limited if the correlation of alleles at the two loci within a colony is taken relative to that of the entire population. In other words, with limited migration, random genetic drift of gamete types prevails in each colony. A possible test in which the variance components of disequilibrium are compared is suggested; the test discriminates between epistatic natural selection and limited migration, showing the former as the main cause of the observed linkage disequilibrium. It is pointed out that the major histocompatibility complex polymorphism is in accord with the population genetics model of multigene families that incorporate gene conversion or double unequal crossing-over between the loci in the supergene family.

title: Strains of Escherichia coli O157:H7 Differ Primarily by Insertions or Deletions, Not Single-Nucleotide Polymorphisms
abstract: Escherichia coli O157:H7 (O157) strains demonstrate varied pulsed-field gel electrophoresis patterns following XbaI digestion, which enable epidemiological surveillance of this important human pathogen. The genetic events underlying PFGE differences between strains, however, are not defined. We investigated the mechanisms for strain variation in O157 by recovering and examining nucleotide sequences flanking each of the XbaI restriction enzyme sites in the genome. Our analysis demonstrated that differences between O157 strains were due to discrete insertions or deletions that contained the XbaI sites polymorphic between strains rather than single-nucleotide polymorphisms in the XbaI sites themselves. These insertions and deletions were found to be uniquely localized within the regions of the genome that are specific to O157 compared to E. coli K-12 (O islands), suggesting that strain-to-strain variation occurs in these O islands. These results may be utilized to devise novel strain-typing tools for this pathogen.

title: Multilocus Sequence Typing Reveals a Lack of Diversity among Escherichia coli O157:H7 Isolates That Are Distinct by Pulsed-Field Gel Electrophoresis
abstract: Escherichia coli O157:H7 is a major cause of foodborne illness in the United States. Pulsed-field gel electrophoresis (PFGE) is the molecular epidemiologic method mostly commonly used to identify food-borne outbreaks. Although PFGE is a powerful epidemiologic tool, it has disadvantages that make a DNA sequence-based approach potentially attractive. Multilocus sequence typing (MLST) analyzes the internal fragments of housekeeping genes to establish genetic relatedness between isolates. We sequenced selected portions of seven housekeeping genes and two membrane protein genes (ompA and espA) of 77 isolates that were diverse by PFGE to determine whether there was sufficient sequence variation to be useful as an epidemiologic tool. There was no DNA sequence diversity in the sequenced portions of the seven housekeeping genes and espA. For ompA, all but five isolates had sequence identical to that of the reference strains. E. coli O157:H7 has a striking lack of genetic diversity in the genes we explored, even among isolates that are clearly distinct by PFGE. Other approaches to identify improved molecular subtyping methods for E. coli 0157:H7 are needed.

title: Recombination and clonality in natural populations of Escherichia coli.
abstract: Bacteria such as Escherichia coli have been commonly viewed as being primarily clonal organisms. As such, the genetic variation within clones was thought to be almost exclusively the result of the mutational process. This conclusion has recently been challenged by data from DNA sequencing studies of natural isolates that are incompatible with a primarily clonal structure. Molecular population genetic analyses of these data, including gene genealogical comparisons, have raised the possibility of a much more complex population structure that may encompass relatively frequent recombination, recurrent selective sweeps and extensive ecological population subdivision.

title: Genome-scale approaches to resolving incongruence in molecular phylogenies
abstract: One of the most pervasive challenges in molecular phylogenetics is the incongruence between phylogenies obtained using different data sets, such as individual genes. To systematically investigate the degree of incongruence, and potential methods for resolving it, we screened the genome sequences of eight yeast species and selected 106 widely distributed orthologous genes for phylogenetic analyses, singly and by concatenation. Our results suggest that data sets consisting of single or a small number of concatenated genes have a significant probability of supporting conflicting topologies. By contrast, analyses of the entire data set of concatenated genes yielded a single, fully resolved species tree with maximum support. Comparable results were obtained with a concatenation of a minimum of 20 genes; substantially more genes than commonly used but a small fraction of any genome. These results have important implications for resolving branches of the tree of life.

title: Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms.
abstract: Traditional and molecular typing schemes for the characterization of pathogenic microorganisms are poorly portable because they index variation that is difficult to compare among laboratories. To overcome these problems, we propose multilocus sequence typing (MLST), which exploits the unambiguous nature and electronic portability of nucleotide sequence data for the characterization of microorganisms. To evaluate MLST, we determined the sequences of approximately 470-bp fragments from 11 housekeeping genes in a reference set of 107 isolates of Neisseria meningitidis from invasive disease and healthy carriers. For each locus, alleles were assigned arbitrary numbers and dendrograms were constructed from the pairwise differences in multilocus allelic profiles by cluster analysis. The strain associations obtained were consistent with clonal groupings previously determined by multilocus enzyme electrophoresis. A subset of six gene fragments was chosen that retained the resolution and congruence achieved by using all 11 loci. Most isolates from hyper-virulent lineages of serogroups A, B, and C meningococci were identical for all loci or differed from the majority type at only a single locus. MLST using six loci therefore reliably identified the major meningococcal lineages associated with invasive disease. MLST can be applied to almost all bacterial species and other haploid organisms, including those that are difficult to cultivate. The overwhelming advantage of MLST over other molecular typing methods is that sequence data are truly portable between laboratories, permitting one expanding global database per species to be placed on a World-Wide Web site, thus enabling exchange of molecular typing data for global epidemiology via the Internet.

title: Simple sequence repeats as a source of quantitative genetic variation.
abstract: Most traits in biological populations appear to be under stabilizing selection, which acts to eliminate quantitative genetic variation. Yet, virtually all measured traits in biological populations continue to show significant quantitative genetic variation. The paradox can be resolved by postulating the existence of an abundant, though unspecified, source of mutations that has quantitative effects on phenotype, but does not reduce fitness. Does such a source actually exist? We propose that it does, in the form of repeat-number variation in SSRs (simple sequence repeats, of which the triplet repeats of human neurodegenerative diseases are a special case). Viewing SSRs as a major source of quantitative mutation has broad implications for understanding molecular processes of evolutionary adaptation, including the evolutionary control of the mutation process itself.

title: Pathological, physiological, and evolutionary aspects of short unstable DNA repeats in the human genome.
abstract: One of the salient features of the mammalian genome is the vast excess of DNA without obvious function, such as repetitive DNAs, spacers, and introns. In recent years, microsatellites, which include short triplet repeats (mostly CAGn and CGGn) and dinucleotide repeats (notably CAn) have gained widespread attention, along with minisatellites which consist of somewhat longer repeat units. Micro- and minisatellites, collectively called variable number tandem repeats (VNTRs), can be highly unstable and display an amazing degree of polymorphism. This property is exploited for gene mapping, for tumor diagnosis, and in forensic medicine. Undue expansion of gene-associated microsatellites is also responsible for some severe genetic diseases, such as fragile X syndrome. Most or all of these diseases are caused by expansion of CAG and CGG triplets. Within protein-coding regions these triplets usually code for polymers of glutamine, serine, alanine or proline. Physiologically, such amino acid repeats are often found in transcription factors and can increase or decrease their activity, depending on the repeat number. Alone or in conjunction with DNA methylation, such repeats may offer a unique opportunity for subtle, semi-stable modulation of gene activity. Also, at least in some plants and perhaps other organisms, a quasi-Lamarckian inheritance is mediated by repetitive DNA. Generally, repetitive DNA sequences, whether represented by short or by long DNA segments, may be beneficial for the evolution of a species.

title: Variable numbers of tandem repeat loci in genetically homogeneous Haemophilus influenzae strains alter during persistent colonisation of cystic fibrosis patients.
abstract: Serial sputum isolates of Haemophilus influenzae (n = 69) were obtained from eight patients suffering from cystic fibrosis. For two of these patients all strains were analysed for polymorphism in the major outer membrane protein profile. For all patients the strains were genetically characterised by random amplification of polymorphic DNA analysis. All strains were included in a survey for polymorphism in regions containing moieties of repetitive DNA as well. A single locus containing trinucleotide repeat units, three loci harbouring tetranucleotides, one region comprising pentanucleotide units and two hexanucleotide repeat unit-containing loci were analysed for repeat number variability. Most of the regions were previously shown to be directly adjacent to or even within virulence genes. All regions behaved as genuine variable number of tandem repeat loci in the sense that genetic polymorphism based on the presence of varying numbers of repeat units could be demonstrated among different strains. Interestingly, several of the repeats showed variation in the absence of the variability as assessed by major outer membrane protein or random amplification of polymorphic DNA analysis. These observations indicate that the repeat loci may vary independently from major chromosomal polymorphism. Consequently, H. influenzae appears to modify its virulence gene regions of the chromosome during persistent colonisation of the lung in cystic fibrosis patients.

title: Adaptive mutation by deletions in small mononucleotide repeats.
abstract: Adaptive reversion of a +1 frameshift mutation in Escherichia coli, which requires homologous recombination functions, is shown here to occur by -1 deletions in regions of small mononucleotide repeats. This pattern makes improbable recombinational mechanisms for adaptive mutation in which blocks of sequences are transferred into the mutating gene, and it supports mechanisms that use DNA polymerase errors. The pattern appears similar to that of mutations found in yeast cells and in hereditary colon cancer cells that are deficient in mismatch repair. These results suggest a recombinational mechanism for adaptive mutation that functions through polymerase errors that persist as a result of a deficiency in post-synthesis mismatch repair.

title: Short sequence repeats in microbial pathogenesis and evolution
abstract: DNA is ubiquitous in microbial genomes. Different classes of short sequence repeats (SSRs) have been identified and demonstrated to be generally heterogeneous in a locus-dependent manner, reflected in variation in the number of repeat units present at a given genomic site or by sequence heterogeneity among individual units. Both types of variability can be used to assess intra-species genetic diversity. Repeat variability often affects the coding potential of the region in which the repetitive element is located. This implies that determination of the primary structure of variable numbers of tandem repeats can be used for epidemiological identification purposes, and also for the analysis of gene function. Precise assessment of SSR structure can also generate insight into the regulation of gene expression. Together, DNA repeat analysis in microbial species provides information on both functional and evolutionary aspects of genetic diversity among microbial isolates.