Inferring the evolutionary history of vibrios by means of multilocus sequence analysis

We performed the first broad study aiming at the reconstruction of the evolutionary history of vibrios by means of multilocus sequence analysis of nine genes. Overall, 14 distinct clades were recognized using the SplitsTree decomposition method. Some of these clades may correspond to families, e.g., the clades Salinivibrio and Photobacteria, while other clades, e.g., Splendidus and Harveyi, correspond to genera. The common ancestor of all vibrios was estimated to have been present 600 million years ago. We can define species of vibrios as groups of strains that share >95% gene sequence similarity and >99.4% amino acid identity based on the eight protein-coding housekeeping genes. The gene sequence data were used to refine the standard online electronic taxonomic scheme for vibrios (http://www.taxvibrio.lncc.br).     

Core gene set as the basis of multilocus sequence analysis of the subclass Actinobacteridae

Comparative genomic sequencing is shedding new light on bacterial identification, taxonomy and phylogeny. An in silico assessment of a core gene set necessary for cellular functioning was made to determine a consensus set of genes that would be useful for the identification, taxonomy and phylogeny of the species belonging to the subclass Actinobacteridae which contained two orders Actinomycetales and Bifidobacteriales. The subclass Actinobacteridae comprised about 85% of the actinobacteria families. The following recommended criteria were used to establish a comprehensive gene set; the gene should (i) be long enough to contain phylogenetically useful information, (ii) not be subject to horizontal gene transfer, (iii) be a single copy (iv) have at least two regions sufficiently conserved that allow the design of amplification and sequencing primers and (v) predict whole-genome relationships. We applied these constraints to 50 different Actinobacteridae genomes and made 1,224 pairwise comparisons of the genome conserved regions and gene fragments obtained by using Sequence VARiability Analysis Program (SVARAP), which allow designing the primers. Following a comparative statistical modeling phase, 3 gene fragments were selected, ychF, rpoB, and secY with R2>0.85. Selected sets of broad range primers were tested from the 3 gene fragments and were demonstrated to be useful for amplification and sequencing of 25 species belonging to 9 genera of Actinobacteridae. The intraspecies similarities were 96.3-100% for ychF, 97.8-100% for rpoB and 96.9-100% for secY among 73 strains belonging to 15 species of the subclass Actinobacteridae compare to 99.4-100% for 16S rRNA. The phylogenetic topology obtained from the combined datasets ychF+rpoB+secY was globally similar to that inferred from the 16S rRNA but with higher confidence. It was concluded that multi-locus sequence analysis using core gene set might represent the first consensus and valid approach for investigating the bacterial identification, phylogeny and taxonomy.     

A multilocus sequence analysis of the genus Xanthomonas

A multilocus sequence analysis (MLSA) of strains representing all validly published Xanthomonas spp. (119 strains) was conducted using four genes; dnaK, fyuA, gyrB and rpoD, a total of 440 sequences. Xanthomonas spp. were divided into two groups similar to those indicated in earlier 16S rDNA comparative analyses, and they possibly represent distinct genera. The analysis clearly differentiated most species that have been established by DNA-DNA reassociation. A similarity matrix of the data indicated clear numerical differences that could form the basis for species differentiation in the future, as an alternative to DNA-DNA reassociation. Some species, X. cynarae, X. gardneri and X. hortorum, formed a single heterogeneous group that is in need of further investigation. X. gardneri appeared to be a synonym of X. cynarae. Recently proposed new species, X. alfalfae, X. citri, X. euvesicatoria, X. fuscans and X. perforans, were not clearly differentiated as species from X. axonopodis, and X. euvesicatoria and X. perforans are very probably synonyms. MLSA offers a powerful tool for further investigation of the classification of Xanthomonas. Based on the dataset produced, the method also offers a relatively simple way of identifying strains as members of known species, or of indicating their status as members of new species.     

Bayesian identification of admixture events using multilocus molecular markers

Bayesian statistical methods for the estimation of hidden genetic structure of populations have gained considerable popularity in the recent years. Utilizing molecular marker data, Bayesian mixture models attempt to identify a hidden population structure by clustering individuals into genetically divergent groups, whereas admixture models target at separating the ancestral sources of the alleles observed in different individuals. We discuss the difficulties involved in the simultaneous estimation of the number of ancestral populations and the levels of admixture in studied individuals' genomes. To resolve this issue, we introduce a computationally efficient method for the identification of admixture events in the population history. Our approach is illustrated by analyses of several challenging real and simulated data sets. The software (baps), implementing the methods introduced here, is freely available at http://www.rni.helsinki.fi/~jic/bapspage.html.     

rpoB sequence analysis as a novel basis for bacterial identification

Comparison of the sequences of conserved genes, most commonly those encoding 16S rRNA, is used for bacterial genotypic identification. Among some taxa, such as the Enterobacteriaceae, variation within this gene does not allow confident species identification. We investigated the usefulness of RNA polymerase beta-subunit encoding gene ( rpoB  ) sequences as an alternative tool for universal bacterial genotypic identification. We generated a database of partial rpoB for 14 Enterobacteriaceae species and then assessed the intra- and interspecies divergence between the rpoB and the 16S rRNA genes by pairwise comparisons. We found that levels of divergence between the rpoB sequences of different strains were markedly higher than those between their 16S rRNA genes. This higher discriminatory power was further confirmed by assigning 20 blindly selected clinical isolates to the correct enteric species on the basis of rpoB sequence comparison. Comparison of rpoB sequences from Enterobacteriaceae was also used as the basis for their phylogenetic analysis and demonstrated the genus Klebsiella to be polyphyletic. The trees obtained with rpoB were more compatible with the currently accepted classification of Enterobacteriaceae than those obtained with 16S rRNA. These data indicate that rpoB is a powerful identification tool, which may be useful for universal bacterial identification.     

Phylogenetic diversity of insecticolous fusaria inferred from multilocus DNA sequence data and their molecular identification via FUSARIUM-ID and Fusarium MLST

We constructed several multilocus DNA sequence datasets to assess the phylogenetic diversity of insecticolous fusaria, especially focusing on those housed at the Agricultural Research Service Collection of Entomopathogenic Fungi (ARSEF), and to aid molecular identifications of unknowns via the FUSARIUM-ID and Fusarium MLST online databases and analysis packages. Analyses of a 190-taxon, two-locus dataset, which included 159 isolates from insects, indicated that: (i) insect-associated fusaria were nested within 10 species complexes spanning the phylogenetic breadth of Fusarium, (ii) novel, putatively unnamed insecticolous species were nested within 8/10 species complexes and (iii) Latin binomials could be applied with confidence to only 18/58 phylogenetically distinct fusaria associated with pest insects. Phylogenetic analyses of an 82-taxon, three-locus dataset nearly fully resolved evolutionary relationships among the 10 clades containing insecticolous fusaria. Multilocus typing of isolates within four species complexes identified surprisingly high genetic diversity in that 63/65 of the fusaria typed represented newly discovered haplotypes. The DNA sequence data, together with corrected ABI sequence chromatograms and alignments, have been uploaded to the following websites dedicated to identifying fusaria: FUSARIUM-ID (http://isolate.fusariumdb.org) at Pennsylvania State University's Department of Plant Pathology and Fusarium MLST (http://www.cbs.knaw.nl/fusarium) at the Centraalbureau voor Schimmelcultures (CBS-KNAW) Fungal Biodiversity Center.     

High-throughput multilocus sequence typing: bringing molecular typing to the next level

Multilocus sequence typing (MLST) is a widely used system for typing microorganisms by sequence analysis of housekeeping genes. The main advantage of MLST in comparison to other typing techniques is the unambiguity and transferability of sequence data. However, a main disadvantage is the high cost of DNA sequencing. Here we introduce a high-throughput MLST (HiMLST) method that employs next-generation sequencing (NGS) technology (Roche 454), to generate large quantities of high-quality MLST data at low costs. The HiMLST protocol consists of two steps. In the first step MLST target genes are amplified by PCR in multi-well plates. During this PCR the amplicons of each bacterial isolate are provided with a unique DNA barcode, the multiplex identifier (MID). In the second step all amplicons are pooled and sequenced in a single NGS-run. The MLST profile of each individual isolate can be retrieved easily using its unique MID. With HiMLST we have profiled 575 isolates of Legionella pneumophila, Staphylococcus aureus, Pseudomonas aeruginosa and Streptococcus pneumoniae in mixed species HiMLST experiments. In conclusion, the introduction of HiMLST paves the way for a broad employment of the MLST as a high-quality and cost-effective method for typing microbial species.     

Phylogenetic analysis of the pathogenic bacteria Spiroplasma penaei based on multilocus sequence analysis

A pathogenic Spiroplasma penaei strain was isolated from the hemolymph of moribund Pacific white shrimp, Penaeus vannamei. The shrimp sample originated from a shrimp farm near Cartagena, Colombia, that was suffering from high mortalities in ponds with very low salinity and high temperatures. This new emerging disease in a marine crustacean in the Americas is described as a systemic infection. The multilocus phylogenetic analysis suggests that S. penaei strain has a terrestrial origin. Evolutionary relationship trees, based on five partial DNA sequences of 16S rDNA, 23S rDNA, 5S rDNA, gyrB, rpoB genes and two complete DNA sequences of 16S-23S rDNA and 23S-5S rDNA intergenic spacer region, were reconstructed using the distance-based Neighboring-Joining (NJ) method with Kimura-2-parameter substitution model. The NJ trees based on all DNA sequences investigated in this study positioned S. penaei in the Citri-Poulsonii clade and corroborate the observations by other investigators using the 16S rDNA gene. Pairwise genetic distance calculation between sequences of spiroplasmas showed S. penaei to be closely related to Spiroplasma insolitum and distantly related to Spiroplasma sp. SHRIMP from China.     

Application of recA and rpoB sequence analysis on phylogeny and molecular identification of Geobacillus species

Aims:  Some Geobacillus species have highly similar 16S rRNA gene sequences, making 16S rDNA sequence analysis-based identification problematic. To overcome this limitation, recA and rpoB sequence analysis was evaluated as an alternative for distinguishing Geobacillus species.
Methods and Results:  The phylogram of 16S rRNA gene sequences inferred from the neighbour-joining method showed that nine clusters of Geobacillus species were characterized with bootstrap values >90%. The recA and rpoB sequences of 10 reference strains in clusters V, VIb and VIc were amplified and sequenced using consensus primers. Alignment of recA sequences in clusters V, VIb and VIc revealed three types of recA genes, consistent with the putative amino acid sequences and in vivo recA splicing analysis. The phylogram constructed from rpoB sequences showed more divergence than that constructed from 16S rRNA gene sequences.
Conclusions:  recA and rpoB sequence analysis differentiated closely-related Geobacillus species and provided direct evidence for reclassifying some species dubiously categorized as Geobacilli . Additionally, this study revealed three types of recA genes in the different Geobacillus species.
Significance and Impact of the Study:  This study highlights the advantage of recA and rpoB sequence analysis to supplement 16S rRNA gene sequence analysis for efficient and convenient determination of Geobacillus species.     

Determination of molecular phylogenetics of Vibrio parahaemolyticus strains by multilocus sequence typing

Vibrio parahaemolyticus is an important human pathogen whose transmission is associated with the consumption of contaminated seafood. There is a growing public health concern due to the emergence of a pandemic strain causing severe outbreaks worldwide. Many questions remain unanswered regarding the evolution and population structure of V. parahaemolyticus. In this work, we describe a multilocus sequence typing (MLST) scheme for V. parahaemolyticus based on the internal fragment sequences of seven housekeeping genes. This MLST scheme was applied to 100 V. parahaemolyticus strains isolated from geographically diverse clinical (n = 37) and environmental (n = 63) sources. The sequences obtained from this work were deposited and are available in a public database (http://pubmlst.org/vparahaemolyticus). Sixty-two unique sequence types were identified, and most (50) were represented by a single isolate, suggesting a high level of genetic diversity. Three major clonal complexes were identified by eBURST analysis. Separate clonal complexes were observed for V. parahaemolyticus isolates originating from the Pacific and Gulf coasts of the United States, while a third clonal complex consisted of strains belonging to the pandemic clonal complex with worldwide distribution. The data reported in this study indicate that V. parahaemolyticus is genetically diverse with a semiclonal population structure and an epidemic structure similar to that of Vibrio cholerae. Genetic diversity in V. parahaemolyticus appears to be driven primarily by frequent recombination rather than mutation, with recombination ratios estimated at 2.5:1 and 8.8:1 by allele and site, respectively. Application of this MLST scheme to more V. parahaemolyticus strains and by different laboratories will facilitate production of a global picture of the epidemiology and evolution of this pathogen.     

Population genetic analysis of Pseudomonas aeruginosa using multilocus sequence typing

To study the population genetic structure of Pseudomonas aeruginosa, we developed a multilocus sequence typing scheme. The sequences of internal fragments of seven housekeeping genes were obtained for 34 P. aeruginosa isolates from patients hospitalized in five different European cities. Twenty-six different allelic profiles were identified. The mean allelic diversity was 0.854 (range: 0.606-0.978), which was about six times greater than the results obtained with the multilocus enzyme electrophoresis method. Linkage disequilibrium was measured with the index of association. An index of 1.95+/-0.24 was calculated when all the strains were considered. This index was 1.76+/-0.27 when only one strain per sequence type was considered. Both results were different from 0, indicating linkage among loci, which means that the population structure of our set of P. aeruginosa isolates is clonal. The clonal structure of the population was also suggested by the congruence of the topology of the different trees obtained from the seven housekeeping genes. These results are in contrast to previous studies, finding a non clonal population structure. Since a small number of isolates was analyzed in this study, there might be a bias of selection which includes the possibility that they belong to widely disseminated epidemic clones. Another possibility is that recombination did not occurred homogeneously throughout the genome of P. aeruginosa, so that part of it has a clonal structure, while the remaining part of the genome is more frequently subject to recombination.     

Proteomics and multilocus sequence analysis confirm intraspecific variability of Vibrio tapetis

Vibrio tapetis is the etiological agent of brown ring disease (BRD) in clams. Phenotypic, antigenic and genetic variability have been demonstrated, with three groups being established associated with host origin. In this work we analyze the variability of representative strains of these three groups, CECT 4600(T) and GR0202RD, isolated from Manila clam and carpet-shell clam, respectively, and HH6087, isolated from halibut, on the basis of the whole proteome analysis by 2D-PAGE and multilocus sequence analysis (MLSA). A quantitative analysis of the proteome match coefficient showed a similarity of 79% between the clam isolates, whereas fish isolate showed similarities lower than 70%. A preliminary mass spectrometry (MS) assay allowed the identification of 27 proteins including 50S ribosomal protein L9, riboflavin synthase β subunit, ribose-phosphate pyrophosphokinase and succinyl-CoA synthase α subunit. The MLSA approach gave similar results, showing a 99.4% similarity of the clam isolates, which was higher than that observed between the fish isolate and either clam strain (98.2%). The topology of the maximum parsimony tree, obtained from 2D-PAGE analysis, and the phylogenetic tree, constructed with the maximum likelihood algorithm from concatenated sequences of 16S rRNA gene and five housekeeping genes (atpA, pyrH, recA, rpoA and rpoD), was very similar, confirming the closer relationship between the two clam isolates.     

Comparison of multilocus variable-number tandem-repeat analysis with multilocus sequence typing and pulsed-field gel electrophoresis for Enterococcus faecalis

Enterococcus faecalis represents recently an important etiological agent of health care-associated infections (HAIs) and there is a need for evaluation and comparison of typing methods available for this microorganism. We tested multilocus VNTR (variable-number tandem repeats) analysis (MLVA) on a well-characterized collection of 153 clinical isolates of E. faecalis, corresponding to 52 multilocus sequence types and 67 pulsed-field gel electrophoresis (PFGE) profiles. MLVA showed high discriminatory power, discerning 111 different types (diversity index equal 98.9%). The concordance MLVA/MLST and MLVA/PFGE was 0.95 and 0.74, respectively. High discriminatory power of MLVA indicates its utility for local epidemiology such as outbreak investigation, and for differentiation of clones defined by other methods.     

Fusarium culmorum is a single phylogenetic species based on multilocus sequence analysis

Fusarium culmorum is a major pathogen of wheat and barley causing head blight and crown rot in cooler temperate climates of Australia, Europe, West Asia and North Africa. To better understand its evolutionary history we partially sequenced single copy nuclear genes encoding translation elongation factor 1-α (TEF), reductase (RED) and phosphate permease (PHO) in 100 F. culmorum isolates with 11 isolates of Fusarium crookwellense, Fusarium graminearum and Fusarium pseudograminearum. Phylogenetic analysis of multilocus sequence (MLS) data using Bayesian inference and maximum parsimony analysis showed that F. culmorum from wheat is a single phylogenetic species with no significant linkage disequilibrium and little or no lineage development along geographic origin. Both MLS and TEF and RED gene sequence analysis separated the four Fusarium species used and delineated three to four groups within the F. culmorum clade. But the PHO gene could not completely resolve isolates into their respective species. Fixation index and gene flow suggest significant genetic exchange between the isolates from distant geographic regions. A lack of strong lineage structure despite the geographic separation of the three collections indicates a frequently recombining species and/or widespread distribution of genotypes due to international trade, tourism and long-range dispersal of macroconidia. Moreover, the two mating type genes were present in equal proportion among the F. culmorum collection used in this study, leaving open the possibility of sexual reproduction.     

A probability matrix for the identification of vibrios

A probability matrix for computer-assisted identification of vibrios has been constructed, based on the API 20E system. Data were gathered from 173 strains representing 31 taxa of vibrios and related organisms, from a variety of sources. The matrix was tested internally by four statistical programs. Program OVERMAT tested the separation and program MOSTTYP the discretion and homogeneity of the taxa. Most of the taxa were satisfactory but a few were less so; reasons for this are discussed. Program CHARSEP and program DIACHAR tested the separation and diagnostic values, respectively, of the characters used. The overall test error was 4–5%. The matrix was assessed externally by its performance in the identification of vibrio-like strains isolated from freshwater. Of 243 wild strains, 79–4% were identified with ten taxa, with a Willcox score of ζ 0–99.     

A probability matrix for the identification of vibrios

A probability matrix for computer-assisted identification of vibrios has been constructed, based on the API 20E system. Data were gathered from 173 strains representing 31 taxa of vibrios and related organisms, from a variety of sources. The matrix was tested internally by four statistical programs. Program OVERMAT tested the separation and program MOSTTYP the discretion and homogeneity of the taxa. Most of the taxa were satisfactory but a few were less so; reasons for this are discussed. Program CHARSEP and program DIACHAR tested the separation and diagnostic values, respectively, of the characters used. The overall test error was 4.5%. The matrix was assessed externally by its performance in the identification of vibrio-like strains isolated from freshwater. Of 243 wild strains, 79.4% were identified with ten taxa, with a Willcox score of greater than or equal to 0.99.     

安徽部分地区动物源弯曲菌的分离鉴定及多位点序列分型

【背景】弯曲菌是一种重要的食源性人兽共患病原菌,革兰氏阴性、微需氧、弯曲螺旋状。【目的】为了解安徽地区弯曲菌流行状况和分子遗传特征,对安徽6个不同地区动物源的弯曲菌进行分离鉴定,并研究分离株分子分型。【方法】通过形态学及培养特性观察、生化试验、PCR方法对菌株进行鉴定。以弯曲菌7个管家基因asp A、gln A、glt A、gly A、pgm、tkt和unc A为目的基因对分离株进行多位点序列分型,并制成遗传进化树。【结果】共分离到42株弯曲菌菌株,源自6个地区的分离株具有较为一致的形态特性和相似的生化特性。多位点序列分型结果显示,本研究中共获得32种ST型,共发现9种新的ST型(8190、8222、8223、8831、8833、8841、8832、8834和8843)和6个新的等位基因(gln A606、gln A607、glt A518、gly A680、pgm863和unc A541)。进化树结果显示,空肠弯曲菌与结肠弯曲菌遗传关系相差甚远,聚集归为两个大群,分别有5个分支和3个分支。【结论】安徽6个地区不同来源的空肠弯曲菌与结肠弯曲菌均有丰富的基因型,且没有明显优势的基因型。从遗传变异的角度来看,空肠弯曲菌复杂多样,结肠弯曲菌相对保守。     

多位点序列和比较基因组学分析对盐单胞菌科菌株JSM 104105的系统发育学研究

【目的】为了精准地判定产铁载体嗜盐新菌株JSM 104105在盐单胞菌科(Halomonadaceae)中的系统发育地位。【方法】采用基于16SrRNA基因、DNA促旋酶B亚基基因(DNAgyraseB subunit gene, gyrB)和核糖核酸聚合酶σ因子RpoD基因(RNA polymerase sigma factor RpoD gene,rpoD)序列的多位点序列分析方法(multilocus sequence analysis, MLSA),对菌株JSM 104105的系统发育地位进行初步分析;采用比较基因组学分析(comparative genomics analysis),分析该菌株与其系统发育关系密切的典型菌株之间的GC含量、平均核苷酸一致性(average nucleotide identity,ANI)和数字DNA-DNA杂交值(digital DNA-DNA hybridization, dDDH),并进行基因组系统发育学分析(phylogenomic analysis),更准确地判定菌株JSM 104105在盐单胞菌科中的系统发育地位。【结果】MLSA分析...     

Biogeography of the ubiquitous marine bacterium Alteromonas macleodii determined by multilocus sequence analysis

Twenty‐three isolates of the widely distributed marine bacteria Alteromonas macleodii have been analysed by multilocus sequence analysis combined with phylogenetic and multivariate statistical analyses. The strains originated from the Pacific Ocean, Mediterranean Sea, English Channel, Black Sea and Thailand. Using the nucleotide sequences of nine loci for each of the 23 isolates, a robust identification was achieved of different clades within the single species. Strains generally clustered with the depth in the water column from which the isolate originated. Strains also showed more recombination with isolates from the same vicinity, suggesting that genetic exchange plays a role in diversification of planktonic marine prokaryotes. This study thus shows for the first time for a large set of isolates of a species of planktonic marine prokaryotes that multilocus sequence analysis overcomes the problems associated with the analysis of individual marker genes or presence of extensive recombination events. It can thus achieve intraspecific identification to the level of genotypes and, by comparison with relevant environmental data, ecotypes.