Biogeographic patterns and current distribution of molecular-genetic variation among populations of speckled dace, Rhinichthys osculus (Girard)

We examined genetic diversity within- and among-populations of speckled dace (Rhinichthys osculus) in five major drainage systems in the state of Oregon in western North America. Analysis of sequence variation in a 670-bp segment of the mitochondrial cytochrome b gene revealed deep divergence among basins and high genetic diversity within basins. Application of a molecular clock indicated that the divergence time among basins reflects vicariant events during the late Miocene to early Pliocene. The high levels of genetic diversity observed within basins is likely due to large historic population sizes, in particular, within the Klamath Basin. Two highly divergent mtDNA lineages were found to co-occur in populations in the Klamath Basin. This result may be indicative of a complex history of isolation and reconnection in this basin and/or multiple colonization events. Based on the observed level of mtDNA divergence these lineages may represent two reproductively isolated sympatric taxa. We recommend that major basins be regarded as distinct ESUs based on high levels of subdivision, deep divergences, and reciprocal monophyly among basins.     

Molecular signatures of divergence and selection in closely related pine taxa

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.     

Rapid postglacial diversification and long‐term stasis within the songbird genus Junco: phylogeographic and phylogenomic evidence

Natural systems composed of closely related taxa that vary in the degree of phenotypic divergence and geographic isolation provide an opportunity to investigate the rate of phenotypic diversification and the relative roles of selection and drift in driving lineage formation. The genus Junco (Aves: Emberizidae) of North America includes parapatric northern forms that are markedly divergent in plumage pattern and colour, in contrast to geographically isolated southern populations in remote areas that show moderate phenotypic divergence. Here, we quantify patterns of phenotypic divergence in morphology and plumage colour and use mitochondrial DNA genes, a nuclear intron, and genomewide SNPs to reconstruct the demographic and evolutionary history of the genus to infer relative rates of evolutionary divergence among lineages. We found that geographically isolated populations have evolved independently for hundreds of thousands of years despite little differentiation in phenotype, in sharp contrast to phenotypically diverse northern forms, which have diversified within the last few thousand years as a result of the rapid postglacial recolonization of North America. SNP data resolved young northern lineages into reciprocally monophyletic lineages, indicating low rates of gene flow even among closely related parapatric forms, and suggesting a role for strong genetic drift or multifarious selection acting on multiple loci in driving lineage divergence. Juncos represent a compelling example of speciation in action, where the combined effects of historical and selective factors have produced one of the fastest cases of speciation known in vertebrates.     

Fossil calibration of molecular clocks and the divergence times of geminate species pairs separated by the isthmus of panama

Calibration of nucleotide sequence divergence rates provides an important method by which to test many hypotheses of evolution. In the absence of an adequate fossil record, geological events, rather than the first appearances of sister taxa in the geological record, are often used to calibrate molecular clocks. The formation of the Isthmus of Panama, which isolated the tropical western Atlantic and eastern Pacific oceans, is one such event that is frequently used to infer rates of nucleotide sequence divergence. Isthmian calibrations assume that morphologically similar "geminate" species living now on either side of the isthmus were isolated geographically by the latest stages of seaway closure 3.1-3.5 MYA. Here, I have applied calibration dates from the fossil record to cytochrome c oxidase-1 (CO1) and nuclear histone-3 (H3) divergences among six pairs of geminates in the Arcidae to test this hypothesis. Analysis of CO1 first and third positions yield geminate divergences that predate final seaway closure, and on the basis of CO1 first positions, times for all six geminates are significantly greater than 3.5 Myr. H3 sequences produce much more recent geminate divergences, some that are younger than 3.1 Myr. But H3-derived estimates for all arcid geminates are not significantly different from both 0 and 15 Myr. According to CO1, one of the two most divergent pairs, Arca mutabilis and A. imbricata, split more than 30 MYA. This date is compatible with the fossil record, which indicates that these species were morphologically distinct at least 16-21 MYA. Across all CO1 nucleotide sites, divergence rates for arcids are slower than the rates reported for other taxa on the basis of isthmian calibrations, with the exception of rates determined from the least divergent species pair in larger surveys of multiple transisthmian pairs. Rate differences between arcids and some taxa may be real, but these data suggest that divergence rates can be greatly overestimated when dates corresponding to final closure of the Central American Seaway are used to calibrate the molecular clocks of marine organisms.     

Allozyme Variation and Genetic Relationships among Species of Cimicifuga (Ranunculaceae) from Korea

Allozyme investigation of the five Cimicifuga taxa in Korea was conducted to assess genetic and clonal diversity within populations and genetic divergence among populations and taxa. Levels of allozyme variation maintained in Korean Cimicifuga taxa were comparable to those for most herbaceous perennials. In general, samples excluding copies of the same multilocus genotype maintained higher levels of genetic diversity than the total samples within populations. Copies of homozygous genotypes at several loci resulting from clonal spread lead to decreased levels of genetic diversity within populations, indicating that clonal reproduction found in Cimicifuga affects population genetic structure. In general, more widely distributed species such as C. dahurica and C. japonica harbored higher levels of allozyme diversity than the other taxa examined. Although two varieties of C. heracleifolia are geographically and reproductively isolated, the genetic and clonal structure of var. bifida seems to resemble var. heracleifolia, indicating that the two varieties may have had a similar evolutionary history. However, the allozyme data strongly indicate that the two morphological types (Groups I and II) of C. simplex should be treated as separate species.     

The tempo of trait divergence in geographic isolation: Avian speciation across the Marañon Valley of Peru

Geographic isolation is considered essential to most speciation events, but our understanding of what controls the pace and degree of phenotypic divergence among allopatric populations remains poor. Why do some taxa exhibit phenotypic differentiation across barriers to dispersal, whereas others do not? To test factors controlling phenotypic divergence in allopatry, we employed a comparative phylogeographic approach consisting of replicates of ecologically similar Andean bird species isolated across a major biogeographic barrier, the Marañon Valley of Peru. Our study design leverages variation among codistributed taxa in their degree of plumage, morphometric, and vocal differentiation across the Marañon to examine the tempo of phenotypic evolution. We found that substantial plumage differences between populations required roughly two million years to evolve. In contrast, morphometric trait evolution showed greater idiosyncrasy and stasis. Our results demonstrate that despite a large degree of idiosyncrasy in the relationship between genetic and phenotypic divergence across taxa and environments, comparative studies within regions may reveal predictability in the pace of phenotypic divergence. Our results also suggest that social selection is important for driving differentiation of populations found in similar environments.     

Patterns and relative rates of nucleotide and insertion/deletion evolution at six chloroplast intergenic regions in new world species of the Lecythidaceae

Insertions and deletions (indels) in chloroplast noncoding regions are common genetic markers to estimate population structure and gene flow, although relatively little is known about indel evolution among recently diverged lineages such as within plant families. Because indel events tend to occur nonrandomly along DNA sequences, recurrent mutations may generate homoplasy for indel haplotypes. This is a potential problem for population studies, because indel haplotypes may be shared among populations after recurrent mutation as well as gene flow. Furthermore, indel haplotypes may differ in fitness and therefore be subject to natural selection detectable as rate heterogeneity among lineages. Such selection could contribute to the spatial patterning of cpDNA haplotypes, greatly complicating the interpretation of cpDNA population structure. This study examined both nucleotide and indel cpDNA variation and divergence at six noncoding regions (psbB-psbH, atpB-rbcL, trnL-trnH, rpl20-5'rps12, trnS-trnG, and trnH-psbA) in 16 individuals from eight species in the Lecythidaceae and a Sapotaceae outgroup. We described patterns of cpDNA changes, assessed the level of indel homoplasy, and tested for rate heterogeneity among lineages and regions. Although regression analysis of branch lengths suggested some degree of indel homoplasy among the most divergent lineages, there was little evidence for indel homoplasy within the Lecythidaceae. Likelihood ratio tests applied to the entire phylogenetic tree revealed a consistent pattern rejecting a molecular clock. Tajima's 1D and 2D tests revealed two taxa with consistent rate heterogeneity, one showing relatively more and one relatively fewer changes than other taxa. In general, nucleotide changes showed more evidence of rate heterogeneity than did indel changes. The rate of evolution was highly variable among the six cpDNA regions examined, with the trnS-trnG and trnH-psbA regions showing as much as 10% and 15% divergence within the Lecythidaceae. Deviations from rate homogeneity in the two taxa were constant across cpDNA regions, consistent with lineage-specific rates of evolution rather than cpDNA region-specific natural selection. There is no evidence that indels are more likely than nucleotide changes to experience homoplasy within the Lecythidaceae. These results support a neutral interpretation of cpDNA indel and nucleotide variation in population studies within species such as Corythophora alta.     

A comparative description of mitochondrial DNA differentiation in selected avian and other vertebrate genera

Levels of mitochondrial DNA (mtDNA) sequence divergence between species within each of several avian (Anas, Aythya, Dendroica, Melospiza, and Zonotrichia) and nonavian (Lepomis and Hyla) vertebrate genera were compared. An analysis of digestion profiles generated by 13-18 restriction endonucleases indicates little overlap in magnitude of mtDNA divergence for the avian versus nonavian taxa examined. In 55 interspecific comparisons among the avian congeners, the fraction of identical fragment lengths (F) ranged from 0.26 to 0.96 (F = 0.46), and, given certain assumptions, these translate into estimates of nucleotide sequence divergence (p) ranging from 0.007 to 0.088; in 46 comparisons among the fish and amphibian congeners, F values ranged from 0.00 to 0.36 (F = 0.09), yielding estimates of P greater than 0.070. The small mtDNA distances among avian congeners are associated with protein-electrophoretic distances (D values) less than approximately 0.2, while the mtDNA distances among assayed fish and amphibian congeners are associated with D values usually greater than 0.4. Since the conservative pattern of protein differentiation previously reported for many avian versus nonavian taxa now appears to be paralleled by a conservative pattern of mtDNA divergence, it seems increasingly likely that many avian species have shared more recent common ancestors than have their nonavian taxonomic counterparts. However, estimates of avian divergence times derived from mtDNA- and protein-calibrated clocks cannot readily be reconciled with some published dates based on limited fossil remains. If the earlier paleontological interpretations are valid, then protein and mtDNA evolution must be somewhat decelerated in birds. The empirical and conceptual issues raised by these findings are highly analogous to those in the long-standing debate about rates of molecular evolution and times of separation of ancestral hominids from African apes.      

Large-subunit ribosomal RNA systematics of symbiotic dinoflagellates: morphology does not recapitulate phylogeny

Biochemical, histological, physiological, and genetic evidence indicates that dinoflagellates symbiotic with marine invertebrates are a heterogeneous complex of taxa, representing at least five genera in three orders. Despite a wealth of data regarding morphological, biochemical, and behavioral differences among symbiotic dinoflagellates, knowledge concerning patterns of diversification is limited. I analyzed approximately 900 bp of the 5' end of the large-subunit ribosomal RNA gene from 14 dinoflagellate isolates: six cultured Symbiodinium specimens, two cultured symbiotic Gymnodinium, two algal samples isolated from reef-building corals, an algal sample obtained from cultures of the jellyfish Cassiopea xamachana, and three free-living Gymnodinium isolates. Results show that morphological similarities among the examined symbiotic taxa do not necessarily correspond with molecular phylogeny. The included Symbiodinium taxa represent a paraphyletic assemblage while Gymnodinium is reconstructed as a polyphyletic assemblage. Analysis indicates that all the included symbiotic dinoflagellates descended from a common, symbiotic ancestor (though within the dinoflagellates, symbiosis is a polyphyletic trait). Additionally, two free-living dinoflagellates emerge within the symbiotic clade, suggesting that the symbiotic lifestyle has been lost at least once in this group. It has been hypothesized that rates of evolution within mutualistic endosymbioses should be reduced relative to free-living taxa. However, results indicate that rates of molecular, morphological, biochemical and behavioral change are similar among branches leading to symbiotic and free-living dinoflagellates.     

Evolution of the genomic recombination rate in murid rodents

Although very closely related species can differ in their fine-scale patterns of recombination hotspots, variation in the average genomic recombination rate among recently diverged taxa has rarely been surveyed. We measured recombination rates in eight species that collectively represent several temporal scales of divergence within a single rodent family, Muridae. We used a cytological approach that enables in situ visualization of crossovers at meiosis to quantify recombination rates in multiple males from each rodent group. We uncovered large differences in genomic recombination rate between rodent species, which were independent of karyotypic variation. The divergence in genomic recombination rate that we document is not proportional to DNA sequence divergence, suggesting that recombination has evolved at variable rates along the murid phylogeny. Additionally, we document significant variation in genomic recombination rate both within and between subspecies of house mice. Recombination rates estimated in F(1) hybrids reveal evidence for sex-linked loci contributing to the evolution of recombination in house mice. Our results provide one of the first detailed portraits of genomic-scale recombination rate variation within a single mammalian family and demonstrate that the low recombination rates in laboratory mice and rats reflect a more general reduction in recombination rate across murid rodents.     

Population genetics and geometric morphometrics of the <Emphasis Type="Italic">Bombus ephippiatus</Emphasis> species complex with implications for its use as a commercial pollinator

Mexico and Central America are among the most biodiverse regions on Earth, harboring many species with high levels of interpopulation morphological and genetic diversity. The mountainous topography of this region contains isolated sky island habitats that have the potential to promote speciation. This has been studied in vertebrates, yet few studies have examined the phylogeographic and genetic structure of insect species encompassing this region. Here we investigate geographic patterns of genetic and morphological divergence and speciation among widespread populations of the highly polymorphic bumble bee Bombus ephippiatus and its closest relative B. wilmattae. We used DNA sequences from a fragment of cytochrome oxidase I (COI), genotypes for twelve microsatellite markers, and morphometric data from wings to construct a well-supported inference of the divergences among these taxa. We have found complex patterns of genetic isolation and morphological divergence within B. ephippiatus across its geographic range and present evidence that B. ephippiatus comprises multiple independent evolutionary lineages. The pattern of their diversification corresponds to geographic and environmental isolating mechanisms, including the Mexican highlands, the lowlands of the Isthmus of Tehuantepec in southern Mexico, the Nicaraguan Depression, the patchily distributed volcanic ranges in Nuclear Central America and Pleistocene glacial cycles. These results have important implications for the development and distribution of B. ephippiatus as a commercial pollinator in Mexico and Central America.     

Phylogeographic patterns and demographic history of Schiedea globosa (Caryophyllaceae) on the Hawaiian Islands

Geomorphological changes have been demonstrated to have had profound impacts on biodiversity, often leading to demographic expansions and contractions and allopatric divergence of taxa. We examined DNA sequence variation at two nuclear and one maternally inherited plastid locus among 10 populations of Schiedea globosa on the Hawaiian Islands to assess the primary factors shaping genetic structure, phylogeographic patterns, and the importance of geographic isolation to population divergence. Schiedea globosa has characteristics that may promote gene flow, including wind pollination and rafting of plants in ocean currents. However, we detected significant differentiation among populations on all islands except Hawaii, with the maternally inherited plastid locus having the greatest genetic structure (F(ST) = 0.81). Migration rates across all loci are less than one migrant per generation. We found evidence of growth in several populations and on the islands of Molokai and Maui, which supports population expansion associated with the formation of Maui Nui during the last glacial maximum. Similar to data for many other Hawaiian taxa, these data suggest S. globosa originated on Oahu and subsequently colonized Molokai, Maui, and Hawaii in progression. Given the high level of genetic structure, allopatric divergence will likely contribute to further divergence of populations.     

Historical Divergence and Gene Flow in the Genus Zea

Gene flow plays a fundamental role in plant evolutionary history, yet its role in population divergence—and ultimately speciation—remains poorly understood. We investigated gene flow and the modalities of divergence in the domesticate Zea mays ssp. mays and three wild Zea taxa using sequence polymorphism data from 26 nuclear loci. We described diversity across loci and assessed evidence for adaptive and purifying selection at nonsynonymous sites. For each of three divergence events in the history of these taxa, we used approximate Bayesian simulation to estimate population sizes and divergence times and explicitly compare among alternative models of divergence. Our estimates of divergence times are surprisingly consistent with previous data from other markers and suggest rapid diversification of lineages within Zea in the last ~150,000 years. We found widespread evidence of historical gene flow, including evidence for divergence in the face of gene flow. We speculate that cultivated maize may serve as a bridge for gene flow among otherwise allopatric wild taxa.     

a molecular phylogeny of apiaceae subfamily Apioideae: evidence from nuclear ribosomal DNA internal transcribed spacer sequences

Phylogenetic relationships among 40 New World and Old World members of Apiaceae subfamily Apioideae, representing seven of the eight tribes and eight of the ten subtribes commonly recognized in the subfamily, were inferred from nucleotide sequence variation in the internal transcribed spacer (ITS) regions of 18-26S nuclear ribosomal DNA. Although the sequences are alignable, with only 11% of sites excluded from the analyses because of alignment ambiguity, divergence values in pairwise comparisons of unambiguous positions among all taxa were high and ranged from 0.5 to 33.2% of nucleotides in ITS 1 and from 0 to 33.2% of nucleotides in ITS 2. Average sequence divergence across both spacer regions was 18.4% of nucleotides. Phylogenies derived from ITS sequences estimated using neighbor-joining analysis of substitution rates, and maximum likelihood and parsimony methods give trees of essentially similar topology and indicate that: (1) there is little support for any existing system of classification of the subfamily that is based largely on morphological and anatomical features of the mericarp; (2) there is a major phylogenetic division within the subfamily, with one clade comprising the genus Smyrnium and those taxa belonging to Drude's tribes Dauceae, Scandiceae, and Laserpitieae and the other clade comprising all other examined taxa; and (3) the genera Arracacia, Coaxana, Coulterophytum, Enantiophylla, Myrrhidendron, Prionosciadium, and Rhodosciadium, all endemic to Mexico and Central America, comprise a clade but their relationships to other New World taxa are equivocal. A phylogeny derived from parsimony analysis of chloroplast DNA rpoC1 intron sequences is consistent with, but considerably less resolved than, relationships derived from these ITS regions. This study affirms that ITS sequences are useful for phylogenetic inference among closely related members of Apioideae but, owing to high rates of nucleotide substitution, are less useful in resolving relationships among the more ancestral nodes of the phylogeny.     

Genetic and plasmid diversity within natural populations of Pseudomonas syringae with various exposures to copper and streptomycin bactericides.

We examined the genetic and plasmid diversity within natural populations of Pseudomonas syringae isolated from three ornamental pear nurseries in eastern Oklahoma. The bactericide spray regimen differed at each nursery; copper and streptomycin, only copper, and no bactericides were applied at nurseries I, II, and III respectively. Resistance to copper (Cur) and resistance to streptomycin (Smr) were determined for 1,938 isolates of P. syringae; isolates from nurseries I and II were generally Cur Sms; whereas most isolates from nursery III were Cus Sms. The plasmid profiles of 362 isolates were determined, and six, one, seven, and four plasmid profiles were obtained for Cur, Smr, Cur Smr, and Cus Sms isolates, respectively. All Smr plasmids contained sequences homologous to the strA and strB Smr genes from broad-host-range plasmid RSF1010 and were associated with Smr transposon Tn5393. Plasmids were placed into two groups on the basis of hybridization to the oriV and par sequences from pOSU900, a cryptic plasmid in P. syringae pv. syringae. A total of 100 randomly chosen P. syringae isolates from nurseries I and III were analyzed for genetic diversity by using the arbitrarily primed PCR (AP-PCR) technique. An analysis of chromosomal genotypes by AP-PCR revealed a high degree of genetic diversity among the isolates, and the results of this analysis indicated that the isolates could be clustered into two distinct groups. The plasmid profiles were specific to isolates belonging to particular AP-PCR groups. Within each AP-PCR group, identical plasmid profiles were produced by isolates that had different chromosomal genotypes, implying that plasmid transfer has played an important role in the dissemination of Cur and Smr within the populations studied.     

Rates of nucleotide substitution in Cornaceae (Cornales)-Pattern of variation and underlying causal factors

Identifying causes of genetic divergence is a central goal in evolutionary biology. Although rates of nucleotide substitution vary among taxa and among genes, the causes of this variation tend to be poorly understood. In the present study, we examined the rate and pattern of molecular evolution for five DNA regions over a phylogeny of Cornus, the single genus of Cornaceae. To identify evolutionary mechanisms underlying the molecular variation, we employed Bayesian methods to estimate divergence times and to infer how absolute rates of synonymous and nonsynonymous substitutions and their ratios change over time. We found that the rates vary among genes, lineages, and through time, and differences in mutation rates, selection type and intensity, and possibly genetic drift all contributed to the variation of substitution rates observed among the major lineages of Cornus. We applied independent contrast analysis to explore whether speciation rates are linked to rates of molecular evolution. The results showed no relationships for individual genes, but suggested a possible localized link between species richness and rate of nonsynonymous nucleotide substitution for the combined cpDNA regions. Furthermore, we detected a positive correlation between rates of molecular evolution and morphological change in Cornus. This was particularly pronounced in the dwarf dogwood lineage, in which genome-wide acceleration in both molecular and morphological evolution has likely occurred.     

DOES RECOMBINATION CONSTRAIN NEUTRAL DIVERGENCE AMONG BACTERIAL TAXA?

A coalescence model for predicting the fate of neutral divergence among closely related taxa distinguishable as separate DNA sequence clusters is presented here. The model simulates iteratively the positive feedback between sequence divergence and sexual isolation among taxa, where increases in sequence divergence result in reduced recombination, and reduced recombination results in increased sequence divergence. Iteration of this feedback is continued until sequence divergence either converges on a steady state or reaches a runaway process. The eventual outcome of sequence divergence was shown to depend on four estimable population-genetic parameters: the expected intrataxon sequence diversity, the baseline rate of intertaxon recombination, the sensitivity of the recombination rate to sequence divergence, and the neutral mutation rate. The model can be used to determine whether neutral divergence among actual taxa is destined to stop at an equilibrium level, or whether neutral divergence will reach a runaway process. Application of the model to the group of taxa containing Bacillus subtilis and its closest relatives showed these taxa to be on a trajectory of unbounded neutral divergence from one another.     

A mitochondrial-DNA-based phylogeny for some evolutionary-genetic model species of Colias butterflies (Lepidoptera, Pieridae)

We study the phylogenetic relationships among some North American Colias ("sulfur") butterflies, using mitochondrial gene sequences (ribosomal RNA, cytochrome oxidase I+II) totaling about 20% of the mitochondrial genome. We find that (1) the lowland species complex shows a branching order different from earlier views; (2) several montane and northern taxa may be more distinct than in earlier views; (3) one morphologically conservative Holarctic assemblage, C. hecla, is differentiated at the molecular-genetic level into at least three taxa which occupy distinct positions in the phylogeny and are sisters to diverse other taxa. These conclusions, constituting phylogenetic hypotheses, are supported by parsimony, maximum-likelihood, and Bayesian reconstruction algorithms. They are tested formally, by interior branch tests and paired-site tests, against alternative hypotheses derived from conventional species and subspecies naming combinations. In all cases our hypotheses are supported by these tests and the conventional alternatives are rejected. The "barcoding" subset of cytochrome oxidase I sequence identifies only some of the taxa supported by our full data set. Comparison of genetic divergence values among Colias taxa with those among related Pierid butterflies suggests that species radiations within Colias are comparatively younger. This emerging Colias phylogeny facilitates comparisons of genetic polymorphism and other adaptive mechanisms among taxa, thereby connecting micro- and macro-evolutionary processes.     

Genetic diversity in Swiss cheese starter cultures assessed by pulsed field gel electrophoresis and arbitrarily primed PCR

AIMS: To assess intraspecific genetic heterogeneity among commercial Swiss cheese starter culture strains of Lactobacillus helveticus, Streptococcus thermophilus and Propionibacterium freudenreichii and to compare the efficacy of two genetic typing methods. METHODS AND RESULTS: Two genetic typing methods, pulsed field gel electrophoresis (PFGE) and arbitrarily primed PCR (AP-PCR), were used. Nine Strep. thermophilus strains revealed eight PFGE and five AP-PCR genotypes. Seventeen Lactobacillus strains yielded 16 and five genotypes by PFGE and AP-PCR, respectively. Eleven Propionibacterium strains yielded 10 PFGE genotypes. Cluster analysis of PFGE profiles generated similarity coefficients for Strep. thermophilus, Lact. helveticus and Prop. freudenreichii strains of 29.5%, 60.3%, and 30.5%, respectively. Milk acidification rates for Strep. thermophilus and Lact. helveticus were determined. CONCLUSIONS: Pulsed field gel electrophoresis is more discriminatory than AP-PCR. The Lact. helveticus group is more homogeneous than the other species examined. Strains with > 87% similarity by PFGE consistently had the same acidification rate and AP-PCR profile. SIGNIFICANCE AND IMPACT OF THE STUDY: Bacterial strains sold for Swiss cheese manufacture in the United States are genetically diverse. Clustering of genetically related bacteria may be useful in identifying new strains with industrially relevant traits.