Role of Geobacter sulfurreducens Outer Surface c-Type Cytochromes in Reduction of Soil Humic Acid and Anthraquinone-2,6-Disulfonate

Deleting individual genes for outer surface c-type cytochromes in Geobacter sulfurreducens partially inhibited the reduction of humic substances and anthraquinone-2,6,-disulfonate. Complete inhibition was obtained only when five of these genes were simultaneously deleted, suggesting that diverse outer surface cytochromes can contribute to the reduction of humic substances and other extracellular quinones.Humic substances can play an important role in the reduction of Fe(III), and possibly other metals, in sedimentary environments (6, 34). Diverse dissimilatory Fe(III)-reducing microorganisms (3, 5, 7, 9, 11, 19-22, 25) can transfer electrons onto the quinone moieties of humic substances (38) or the model compound anthraquinone-2,6-disulfonate (AQDS). Reduced humic substances or AQDS abiotically reduces Fe(III) to Fe(II), regenerating the quinone. Electron shuttling in this manner can greatly increase the rate of electron transfer to insoluble Fe(III) oxides, presumably because soluble quinone-containing molecules are more accessible for microbial reduction than insoluble Fe(III) oxides (19, 22). Thus, catalytic amounts of humic substances have the potential to dramatically influence rates of Fe(III) reduction in soils and sediments and can promote more rapid degradation of organic contaminants coupled to Fe(III) reduction (1, 2, 4, 10, 24).To our knowledge, the mechanisms by which Fe(III)-reducing microorganisms transfer electrons to humic substances have not been investigated previously for any microorganism. However, reduction of AQDS has been studied using Shewanella oneidensis (17, 40). Disruption of the gene for MtrB, an outer membrane protein required for proper localization of outer membrane cytochromes (31), inhibited reduction of AQDS, as did disruption of the gene for the outer membrane c-type cytochrome, MtrC (17). However, in each case inhibition was incomplete, and it was suggested that there was a possibility of some periplasmic reduction (17), which would be consistent with the ability of AQDS to enter the cell (40).The mechanisms for electron transfer to humic substances in Geobacter species are of interest because molecular studies have frequently demonstrated that Geobacter species are the predominant Fe(III)-reducing microorganisms in sedimentary environments in which Fe(III) reduction is an important process (references 20, 32, and 42 and references therein). Geobacter sulfurreducens has routinely been used for investigations of the physiology of Geobacter species because of the availability of its genome sequence (29), a genetic system (8), and a genome-scale metabolic model (26) has made it possible to take a systems biology approach to understanding the growth of this organism in sedimentary environments (23).     

Alignment of the c-Type Cytochrome OmcS along Pili of Geobacter sulfurreducens

Immunogold localization revealed that OmcS, a cytochrome that is required for Fe(III) oxide reduction by Geobacter sulfurreducens, was localized along the pili. The apparent spacing between OmcS molecules suggests that OmcS facilitates electron transfer from pili to Fe(III) oxides rather than promoting electron conduction along the length of the pili.There are multiple competing/complementary models for extracellular electron transfer in Fe(III)- and electrode-reducing microorganisms (8, 18, 20, 44). Which mechanisms prevail in different microorganisms or environmental conditions may greatly influence which microorganisms compete most successfully in sedimentary environments or on the surfaces of electrodes and can impact practical decisions on the best strategies to promote Fe(III) reduction for bioremediation applications (18, 19) or to enhance the power output of microbial fuel cells (18, 21).The three most commonly considered mechanisms for electron transfer to extracellular electron acceptors are (i) direct contact between redox-active proteins on the outer surfaces of the cells and the electron acceptor, (ii) electron transfer via soluble electron shuttling molecules, and (iii) the conduction of electrons along pili or other filamentous structures. Evidence for the first mechanism includes the necessity for direct cell-Fe(III) oxide contact in Geobacter species (34) and the finding that intensively studied Fe(III)- and electrode-reducing microorganisms, such as Geobacter sulfurreducens and Shewanella oneidensis MR-1, display redox-active proteins on their outer cell surfaces that could have access to extracellular electron acceptors (1, 2, 12, 15, 27, 28, 31-33). Deletion of the genes for these proteins often inhibits Fe(III) reduction (1, 4, 7, 15, 17, 28, 40) and electron transfer to electrodes (5, 7, 11, 33). In some instances, these proteins have been purified and shown to have the capacity to reduce Fe(III) and other potential electron acceptors in vitro (10, 13, 29, 38, 42, 43, 48, 49).Evidence for the second mechanism includes the ability of some microorganisms to reduce Fe(III) that they cannot directly contact, which can be associated with the accumulation of soluble substances that can promote electron shuttling (17, 22, 26, 35, 36, 47). In microbial fuel cell studies, an abundance of planktonic cells and/or the loss of current-producing capacity when the medium is replaced is consistent with the presence of an electron shuttle (3, 14, 26). Furthermore, a soluble electron shuttle is the most likely explanation for the electrochemical signatures of some microorganisms growing on an electrode surface (26, 46).Evidence for the third mechanism is more circumstantial (19). Filaments that have conductive properties have been identified in Shewanella (7) and Geobacter (41) species. To date, conductance has been measured only across the diameter of the filaments, not along the length. The evidence that the conductive filaments were involved in extracellular electron transfer in Shewanella was the finding that deletion of the genes for the c-type cytochromes OmcA and MtrC, which are necessary for extracellular electron transfer, resulted in nonconductive filaments, suggesting that the cytochromes were associated with the filaments (7). However, subsequent studies specifically designed to localize these cytochromes revealed that, although the cytochromes were extracellular, they were attached to the cells or in the exopolymeric matrix and not aligned along the pili (24, 25, 30, 40, 43). Subsequent reviews of electron transfer to Fe(III) in Shewanella oneidensis (44, 45) appear to have dropped the nanowire concept and focused on the first and second mechanisms.Geobacter sulfurreducens has a number of c-type cytochromes (15, 28) and multicopper proteins (12, 27) that have been demonstrated or proposed to be on the outer cell surface and are essential for extracellular electron transfer. Immunolocalization and proteolysis studies demonstrated that the cytochrome OmcB, which is essential for optimal Fe(III) reduction (15) and highly expressed during growth on electrodes (33), is embedded in the outer membrane (39), whereas the multicopper protein OmpB, which is also required for Fe(III) oxide reduction (27), is exposed on the outer cell surface (39).OmcS is one of the most abundant cytochromes that can readily be sheared from the outer surfaces of G. sulfurreducens cells (28). It is essential for the reduction of Fe(III) oxide (28) and for electron transfer to electrodes under some conditions (11). Therefore, the localization of this important protein was further investigated.     

Abundance of Novel and Diverse tfdA-Like Genes,Encoding Putative Phenoxyalkanoic Acid Herbicide-Degrading Dioxygenases,in Soil

Phenoxyalkanoic acid (PAA) herbicides are widely used in agriculture. Biotic degradation of such herbicides occurs in soils and is initiated by α-ketoglutarate- and Fe²⁺-dependent dioxygenases encoded by tfdA-like genes (i.e., tfdA and tfdAα). Novel primers and quantitative kinetic PCR (qPCR) assays were developed to analyze the diversity and abundance of tfdA-like genes in soil. Five primer sets targeting tfdA-like genes were designed and evaluated. Primer sets 3 to 5 specifically amplified tfdA-like genes from soil, and a total of 437 sequences were retrieved. Coverages of gene libraries were 62 to 100%, up to 122 genotypes were detected, and up to 389 genotypes were predicted to occur in the gene libraries as indicated by the richness estimator Chao1. Phylogenetic analysis of in silico-translated tfdA-like genes indicated that soil tfdA-like genes were related to those of group 2 and 3 Bradyrhizobium spp., Sphingomonas spp., and uncultured soil bacteria. Soil-derived tfdA-like genes were assigned to 11 clusters, 4 of which were composed of novel sequences from this study, indicating that soil harbors novel and diverse tfdA-like genes. Correlation analysis of 16S rRNA and tfdA-like gene similarity indicated that any two bacteria with D > 20% of group 2 tfdA-like gene-derived protein sequences belong to different species. Thus, data indicate that the soil analyzed harbors at least 48 novel bacterial species containing group 2 tfdA-like genes. Novel qPCR assays were established to quantify such new tfdA-like genes. Copy numbers of tfdA-like genes were 1.0 × 10⁶ to 65 × 10⁶ per gram (dry weight) soil in four different soils, indicating that hitherto-unknown, diverse tfdA-like genes are abundant in soils.Phenoxyalkanoic acid (PAA) herbicides such as MCPA (4-chloro-2-methyl-phenoxyacetic acid) and 2,4-D (2,4-dichlorophenoxyacetic acid) are widely used to control broad-leaf weeds in agricultural as well as nonagricultural areas (19, 77). Degradation occurs primarily under oxic conditions in soil, and microorganisms play a key role in the degradation of such herbicides in soil (62, 64). Although relatively rapidly degraded in soil (32, 45), both MCPA and 2,4-D are potential groundwater contaminants (10, 56, 70), accentuating the importance of bacterial PAA herbicide-degrading bacteria in soils (e.g., references 3, 5, 6, 20, 41, 59, and 78).Degradation can occur cometabolically or be associated with energy conservation (15, 54). The first step in the degradation of 2,4-D and MCPA is initiated by the product of cadAB or tfdA-like genes (29, 30, 35, 67), which constitutes an α-ketoglutarate (α-KG)- and Fe²⁺-dependent dioxygenase. TfdA removes the acetate side chain of 2,4-D and MCPA to produce 2,4-dichlorophenol and 4-chloro-2-methylphenol, respectively, and glyoxylate while oxidizing α-ketoglutarate to CO₂ and succinate (16, 17).Organisms capable of PAA herbicide degradation are phylogenetically diverse and belong to the Alpha-, Beta-, and Gammproteobacteria and the Bacteroidetes/Chlorobi group (e.g., references 2, 14, 29-34, 39, 60, 68, and 71). These bacteria harbor tfdA-like genes (i.e., tfdA or tfdAα) and are categorized into three groups on an evolutionary and physiological basis (34). The first group consists of beta- and gammaproteobacteria and can be further divided into three distinct classes based on their tfdA genes (30, 46). Class I tfdA genes are closely related to those of Cupriavidus necator JMP134 (formerly Ralstonia eutropha). Class II tfdA genes consist of those of Burkholderia sp. strain RASC and a few strains that are 76% identical to class I tfdA genes. Class III tfdA genes are 77% identical to class I and 80% identical to class II tfdA genes and linked to MCPA degradation in soil (3). The second group consists of alphaproteobacteria, which are closely related to Bradyrhizobium spp. with tfdAα genes having 60% identity to tfdA of group 1 (18, 29, 34). The third group also harbors the tfdAα genes and consists of Sphingomonas spp. within the alphaproteobacteria (30).Diverse PAA herbicide degraders of all three groups were identified in soil by cultivation-dependent studies (32, 34, 41, 78). Besides CadAB, TfdA and certain TfdAα proteins catalyze the conversion of PAA herbicides (29, 30, 35). All groups of tfdA-like genes are potentially linked to the degradation of PAA herbicides, although alternative primary functions of group 2 and 3 TfdAs have been proposed (30, 35). However, recent cultivation-independent studies focused on 16S rRNA genes or solely on group 1 tfdA sequences in soil (e.g., references 3-5, 13, and 41). Whether group 2 and 3 tfdA-like genes are also quantitatively linked to the degradation of PAA herbicides in soils is unknown. Thus, tools to target a broad range of tfdA-like genes are needed to resolve such an issue. Primers used to assess the diversity of tfdA-like sequences used in previous studies were based on the alignment of approximately 50% or less of available sequences to date (3, 20, 29, 32, 39, 47, 58, 73). Primers specifically targeting all major groups of tfdA-like genes to assess and quantify a broad diversity of potential PAA degraders in soil are unavailable. Thus, the objectives of this study were (i) to develop primers specific for all three groups of tfdA-like genes, (ii) to establish quantitative kinetic PCR (qPCR) assays based on such primers for different soil samples, and (iii) to assess the diversity and abundance of tfdA-like genes in soil.     

Detection,Distribution, and Organohalogen Compound Discovery Implications of the Reduced Flavin Adenine Dinucleotide-Dependent Halogenase Gene in Major Filamentous Actinomycete Taxonomic Groups

Halogenases have been shown to play a significant role in biosynthesis and introducing the bioactivity of many halogenated secondary metabolites. In this study, 54 reduced flavin adenine dinucleotide (FADH₂)-dependent halogenase gene-positive strains were identified after the PCR screening of a large collection of 228 reference strains encompassing all major families and genera of filamentous actinomycetes. The wide distribution of this gene was observed to extend to some rare lineages with higher occurrences and large sequence diversity. Subsequent phylogenetic analyses revealed that strains containing highly homologous halogenases tended to produce halometabolites with similar structures, and halogenase genes are likely to propagate by horizontal gene transfer as well as vertical inheritance within actinomycetes. Higher percentages of halogenase gene-positive strains than those of halogenase gene-negative ones contained polyketide synthase genes and/or nonribosomal peptide synthetase genes or displayed antimicrobial activities in the tests applied, indicating their genetic and physiological potentials for producing secondary metabolites. The robustness of this halogenase gene screening strategy for the discovery of particular biosynthetic gene clusters in rare actinomycetes besides streptomycetes was further supported by genome-walking analysis. The described distribution and phylogenetic implications of the FADH₂-dependent halogenase gene present a guide for strain selection in the search for novel organohalogen compounds from actinomycetes.It is well known that actinomycetes, notably filamentous actinomycetes, have a remarkable capacity to produce bioactive molecules for drug development (4, 6). However, novel technologies are demanded for the discovery of new bioactive secondary metabolites from these microbes to meet the urgent medical need for drug candidates (5, 9, 31).Genome mining recently has been used to search for new drug leads (7, 20, 42, 51). Based on the hypothesis that secondary metabolites with similar structures are biosynthesized by gene clusters that harbor certain homologous genes, such homologous genes could serve as suitable markers for distinct natural-product gene clusters (26, 51). A wide range of structurally diverse bioactive compounds are synthesized by polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) systems in actinomycetes, therefore much attention has been given to revealing a previously unrecognized biosynthetic potential of actinomycetes through the genome mining of these genes (2, 3, 22). However, the broad distribution of PKS and NRPS genes and their high numbers even in a single actinomycete complicate their use (2, 3). To rationally exploit the genetic potential of actinomycetes, more and more special genes, such as tailoring enzyme genes, are being utilized for this sequence-guided genetic screening strategy (20, 38).Tailoring enzymes, which are responsible for the introduction and generation of diversity and bioactivity in several structural classes during or after NRPS, PKS, or NRPS/PKS assembly lines, usually include acyltransferases, aminotransferases, cyclases, glycosyltransferases, halogenases, ketoreductases, methyltransferases, and oxygenases (36, 45). Halogenation, an important feature for the bioactivity of a large number of distinct natural products (16, 18, 30), frequently is introduced by one type of halogenase, called reduced flavin adenine dinucleotide (FADH₂)-dependent (or flavin-dependent) halogenase (10, 12, 35). More than 4,000 halometabolites have been discovered (15), including commercially important antibiotics such as chloramphenicol, vancomycin, and teicoplanin (43).Previous investigations of FADH₂-dependent halogenase genes were focused largely on related gene clusters in the genera Amycolatopsis (33, 44, 53) and Streptomyces (8, 10, 21, 27, 32, 34, 47-49) and also on those in the genera Actinoplanes (25), Actinosynnema (50), Micromonospora (1), and Nonomuraea (39); however, none of these studies has led to the rest of the major families and genera of actinomycetes. In addition, there is evidence that FADH₂-dependent halogenase genes of streptomycetes usually exist in halometabolite biosynthetic gene clusters (20), but we lack knowledge of such genes and clusters in other actinomycetes.In the present study, we show that the distribution of the FADH₂-dependent halogenase gene in filamentous actinomycetes does indeed correlate with the potential for halometabolite production based on other genetic or physiological factors. We also showed that genome walking near the halogenase gene locus could be employed to identify closely linked gene clusters that likely encode pathways for organohalogen compound production in actinomycetes other than streptomycetes.     

Cell-to-Cell Spread of the RNA Interference Response Suppresses Semliki Forest Virus (SFV) Infection of Mosquito Cell Cultures and Cannot Be Antagonized by SFV

In their vertebrate hosts, arboviruses such as Semliki Forest virus (SFV) (Togaviridae) generally counteract innate defenses and trigger cell death. In contrast, in mosquito cells, following an early phase of efficient virus production, a persistent infection with low levels of virus production is established. Whether arboviruses counteract RNA interference (RNAi), which provides an important antiviral defense system in mosquitoes, is an important question. Here we show that in Aedes albopictus-derived mosquito cells, SFV cannot prevent the establishment of an antiviral RNAi response or prevent the spread of protective antiviral double-stranded RNA/small interfering RNA (siRNA) from cell to cell, which can inhibit the replication of incoming virus. The expression of tombusvirus siRNA-binding protein p19 by SFV strongly enhanced virus spread between cultured cells rather than virus replication in initially infected cells. Our results indicate that the spread of the RNAi signal contributes to limiting virus dissemination.In animals, RNA interference (RNAi) was first described for Caenorhabditis elegans (27). The production or introduction of double-stranded RNA (dsRNA) in cells leads to the degradation of mRNAs containing homologous sequences by sequence-specific cleavage of mRNAs. Central to RNAi is the production of 21- to 26-nucleotide small interfering RNAs (siRNAs) from dsRNA and the assembly of an RNA-induced silencing complex (RISC), followed by the degradation of the target mRNA (23, 84). RNAi is a known antiviral strategy of plants (3, 53) and insects (21, 39, 51). Study of Drosophila melanogaster in particular has given important insights into RNAi responses against pathogenic viruses and viral RNAi inhibitors (31, 54, 83, 86, 91). RNAi is well characterized for Drosophila, and orthologs of antiviral RNAi genes have been found in Aedes and Culex spp. (13, 63).Arboviruses, or arthropod-borne viruses, are RNA viruses mainly of the families Bunyaviridae, Flaviviridae, and Togaviridae. The genus Alphavirus within the family Togaviridae contains several mosquito-borne pathogens: arboviruses such as Chikungunya virus (16) and equine encephalitis viruses (88). Replication of the prototype Sindbis virus and Semliki Forest virus (SFV) is well understood (44, 71, 74, 79). Their genome consists of a positive-stranded RNA with a 5′ cap and a 3′ poly(A) tail. The 5′ two-thirds encodes the nonstructural polyprotein P1234, which is cleaved into four replicase proteins, nsP1 to nsP4 (47, 58, 60). The structural polyprotein is encoded in the 3′ one-third of the genome and cleaved into capsid and glycoproteins after translation from a subgenomic mRNA (79). Cytoplasmic replication complexes are associated with cellular membranes (71). Viruses mature by budding at the plasma membrane (35).In nature, arboviruses are spread by arthropod vectors (predominantly mosquitoes, ticks, flies, and midges) to vertebrate hosts (87). Little is known about how arthropod cells react to arbovirus infection. In mosquito cell cultures, an acute phase with efficient virus production is generally followed by the establishment of a persistent infection with low levels of virus production (9). This is fundamentally different from the cytolytic events following arbovirus interactions with mammalian cells and pathogenic insect viruses with insect cells. Alphaviruses encode host response antagonists for mammalian cells (2, 7, 34, 38).RNAi has been described for mosquitoes (56) and, when induced before infection, antagonizes arboviruses and their replicons (1, 4, 14, 15, 29, 30, 32, 42, 64, 65). RNAi is also functional in various mosquito cell lines (1, 8, 43, 49, 52). In the absence of RNAi, alphavirus and flavivirus replication and/or dissemination is enhanced in both mosquitoes and Drosophila (14, 17, 31, 45, 72). RNAi inhibitors weakly enhance SFV replicon replication in tick and mosquito cells (5, 33), posing the questions of how, when, and where RNAi interferes with alphavirus infection in mosquito cells.Here we use an A. albopictus-derived mosquito cell line to study RNAi responses to SFV. Using reporter-based assays, we demonstrate that SFV cannot avoid or efficiently inhibit the establishment of an RNAi response. We also demonstrate that the RNAi signal can spread between mosquito cells. SFV cannot inhibit cell-to-cell spread of the RNAi signal, and spread of the virus-induced RNAi signal (dsRNA/siRNA) can inhibit the replication of incoming SFV in neighboring cells. Furthermore, we show that SFV expression of a siRNA-binding protein increases levels of virus replication mainly by enhancing virus spread between cells rather than replication in initially infected cells. Taken together, these findings suggest a novel mechanism, cell-to-cell spread of antiviral dsRNA/siRNA, by which RNAi limits SFV dissemination in mosquito cells.     

Conjugative Plasmid Transfer and Adhesion Dynamics in an Escherichia coli Biofilm

A conjugative plasmid from the catheter-associated urinary tract infection strain Escherichia coli MS2027 was sequenced and annotated. This 42,644-bp plasmid, designated pMAS2027, contains 58 putative genes and is most closely related to plasmids belonging to incompatibility group X (IncX1). Plasmid pMAS2027 encodes two important virulence factors: type 3 fimbriae and a type IV secretion (T4S) system. Type 3 fimbriae, recently found to be functionally expressed in E. coli, played an important role in biofilm formation. Biofilm formation by E. coli MS2027 was specifically due to expression of type 3 fimbriae and not the T4S system. The T4S system, however, accounted for the conjugative ability of pMAS2027 and enabled a non-biofilm-forming strain to grow as part of a mixed biofilm following acquisition of this plasmid. Thus, the importance of conjugation as a mechanism to spread biofilm determinants was demonstrated. Conjugation may represent an important mechanism by which type 3 fimbria genes are transferred among the Enterobacteriaceae that cause device-related infections in nosocomial settings.Bacterial biofilms are complex communities of bacterial cells living in close association with a surface (17). Bacterial cells in these protected environments are often resistant to multiple factors, including antimicrobials, changes in the pH, oxygen radicals, and host immune defenses (19, 38). Biofilm formation is a property of many bacterial species, and a range of molecular mechanisms that facilitate this process have been described (2, 3, 11, 14, 16, 29, 33, 34). Often, the ability to form a biofilm is dependent on the production of adhesins on the bacterial cell surface. In Escherichia coli, biofilm formation is enhanced by the production of certain types of fimbriae (e.g., type 1 fimbriae, type 3 fimbriae, F1C, F9, curli, and conjugative pili) (14, 23, 25, 29, 33, 39, 46), cell surface adhesins (e.g., autotransporter proteins such as antigen 43, AidA, TibA, EhaA, and UpaG) (21, 34, 35, 40, 43), and flagella (22, 45).The close proximity of bacterial cells in biofilms creates an environment conducive for the exchange of genetic material. Indeed, plasmid-mediated conjugation in monospecific and mixed E. coli biofilms has been demonstrated (6, 18, 24, 31). The F plasmid represents the best-characterized conjugative system for biofilm formation by E. coli. The F pilus mediates adhesion to abiotic surfaces and stabilizes the biofilm structure through cell-cell interactions (16, 30). Many other conjugative plasmids also contribute directly to biofilm formation upon derepression of the conjugative function (16).One example of a conjugative system employed by gram-negative Enterobacteriaceae is the type 4 secretion (T4S) system. The T4S system is a multisubunit structure that spans the cell envelope and contains a secretion channel often linked to a pilus or other surface filament or protein (8). The Agrobacterium tumefaciens VirB-VirD4 system is the archetypical T4S system and is encoded by 11 genes in the virB operon and one gene (virD4) in the virD operon (7, 8). Genes with strong homology to genes in the virB operon have also been identified on other conjugative plasmids. For example, the pilX1 to pilX11 genes on the E. coli R6K IncX plasmid and the virB1 to virB11 genes are highly conserved at the nucleotide level (28).We recently described identification and characterization of the mrk genes encoding type 3 fimbriae in a uropathogenic strain of E. coli isolated from a patient with a nosocomial catheter-associated urinary tract infection (CAUTI) (29). The mrk genes were located on a conjugative plasmid (pMAS2027) and were strongly associated with biofilm formation. In this study we determined the entire sequence of plasmid pMAS2027 and revealed the presence of conjugative transfer genes homologous to the pilX1 to pilX11 genes of E. coli R6K (in addition to the mrk genes). We show here that biofilm formation is driven primarily by type 3 fimbriae and that the T4S apparatus is unable to mediate biofilm growth in the absence of the mrk genes. Finally, we demonstrate that conjugative transfer of pMAS2027 within a mixed biofilm confers biofilm formation properties on recipient cells due to acquisition of the type 3 fimbria-encoding mrk genes.     

A New Borrelia Species Defined by Multilocus Sequence Analysis of Housekeeping Genes

Analysis of Lyme borreliosis (LB) spirochetes, using a novel multilocus sequence analysis scheme, revealed that OspA serotype 4 strains (a rodent-associated ecotype) of Borrelia garinii were sufficiently genetically distinct from bird-associated B. garinii strains to deserve species status. We suggest that OspA serotype 4 strains be raised to species status and named Borrelia bavariensis sp. nov. The rooted phylogenetic trees provide novel insights into the evolutionary history of LB spirochetes.Multilocus sequence typing (MLST) and multilocus sequence analysis (MLSA) have been shown to be powerful and pragmatic molecular methods for typing large numbers of microbial strains for population genetics studies, delineation of species, and assignment of strains to defined bacterial species (4, 13, 27, 40, 44). To date, MLST/MLSA schemes have been applied only to a few vector-borne microbial populations (1, 6, 30, 37, 40, 41, 47).Lyme borreliosis (LB) spirochetes comprise a diverse group of zoonotic bacteria which are transmitted among vertebrate hosts by ixodid (hard) ticks. The most common agents of human LB are Borrelia burgdorferi (sensu stricto), Borrelia afzelii, Borrelia garinii, Borrelia lusitaniae, and Borrelia spielmanii (7, 8, 12, 35). To date, 15 species have been named within the group of LB spirochetes (6, 31, 32, 37, 38, 41). While several of these LB species have been delineated using whole DNA-DNA hybridization (3, 20, 33), most ecological or epidemiological studies have been using single loci (5, 9-11, 29, 34, 36, 38, 42, 51, 53). Although some of these loci have been convenient for species assignment of strains or to address particular epidemiological questions, they may be unsuitable to resolve evolutionary relationships among LB species, because it is not possible to define any outgroup. For example, both the 5S-23S intergenic spacer (5S-23S IGS) and the gene encoding the outer surface protein A (ospA) are present only in LB spirochete genomes (36, 43). The advantage of using appropriate housekeeping genes of LB group spirochetes is that phylogenetic trees can be rooted with sequences of relapsing fever spirochetes. This renders the data amenable to detailed evolutionary studies of LB spirochetes.LB group spirochetes differ remarkably in their patterns and levels of host association, which are likely to affect their population structures (22, 24, 46, 48). Of the three main Eurasian Borrelia species, B. afzelii is adapted to rodents, whereas B. valaisiana and most strains of B. garinii are maintained by birds (12, 15, 16, 23, 26, 45). However, B. garinii OspA serotype 4 strains in Europe have been shown to be transmitted by rodents (17, 18) and, therefore, constitute a distinct ecotype within B. garinii. These strains have also been associated with high pathogenicity in humans, and their finer-scale geographical distribution seems highly focal (10, 34, 52, 53).In this study, we analyzed the intra- and interspecific phylogenetic relationships of B. burgdorferi, B. afzelii, B. garinii, B. valaisiana, B. lusitaniae, B. bissettii, and B. spielmanii by means of a novel MLSA scheme based on chromosomal housekeeping genes (30, 48).