Retrospective analysis of Vibrio spp. isolated from marketed crustaceans using multilocus sequence analysis

Objective The phytoplasmas in 16SrI group cause severe diseases of many crops and ecological plants in China. At present, the genetic variation and population structure of the yet-uncultured and very closely related phytoplasma strains in China are still not fully understood. A multilocus sequence analysis (MLSA) scheme was used to elucidate genetic diversity of phytoplasma strains in 16SrI group from different regions in China and their relationships with geographical distribution, and to compare the levels of variation in different housekeeping gene of different phytoplasma strains. The study aimed to provide some reference approach and evidence for finer detection, identification, classification and phylogenetic analysis of different phytoplasma strains in China. Method Ten housekeeping gene (rp, tuf, secA, secY, ipt, dnaK, fusA, gyrB, pyrG and rpoB) fragments combined with 16S rDNA sequences were employed to analyze the genetic variation and phylogenetic relationships of 18 phytoplasma strains infecting chinaberry, lettuce, mulberry, paulownia and periwinkle from ten provinces in China, using five phytoplasma strains, onion yellows phytoplasma (OY-M), aster yellows witches'-broom phytoplasma (AYWB), Candidatus Phytoplasma australiense (CPA), strawberry lethal yellows phytoplasma (SLY) and Candidatus Phytoplasma mali (CPM), of which whole genomes have been sequenced, as references. Sequence polymorphism and variation levels of different gene fragments were analyzed by multiple sequence alignment. Result The nucleotide site polymorphisms of rp, tuf, secA, secY, ipt, dnaK, fusA, gyrB, pyrG and rpoB gene fragments were used to resolve all strains into 15 sequence types (STs), demonstrating extensive genetic diversity among the 16SrI group strain population. All the strains, classified into 16SrI-B and-D subgroup by 16S rDNA analysis, were clustered into one clade and clearly differentiated into discrete subclades by phylogenetic analysis of the concatenated gene sequences. Ten chinaberry witches'-broom strains, which were most closely related to two mulberry dwarf strains and hardly distinguished with 16S rDNA, were definitely split into four distinct clusters and 8 STs apparently relative to their geographical origins. Two lettuce yellows strains in Sanming, Fujian province, China, were more closely related to the onion yellows OY-M strain in Japan than the periwinkle virescence and paulownia witches'-broom strains in China. The levels of variation in dnaK locus were higher than those in 16S rDNA and other genes tested. Conclusion It is suggested that the MLSA should potentially be a useful and reliable approach for phytoplasma identification and differentiation as well as for depth examination of strain diversity and this method can be wildly used to study the genetic variation and evolutionary relationships of phytoplasma strains among various groups or subgroups in the future.

Streptomyces griseus and related species form the biggest but least well-defined clade in the whole Streptomyces 16S rRNA gene tree. Multilocus sequence analysis (MLSA) has shown promising potential for refining Streptomyces systematics. In this investigation, strains of 18 additional S. griseus clade species were analysed and data from a previous pilot study were integrated in a larger MLSA phylogeny. The results demonstrated that MLSA of five housekeeping genes (atpD, gyrB, recA, rpoB and trpB) is better than the previous six-gene scheme, as it provides equally good resolution and stability and is more cost-effective; MLSA using three or four of the genes also shows good resolution and robustness for differentiating most of the strains and is therefore of value for everyday use. MLSA is more suitable for discriminating strains that show >99 % 16S rRNA gene sequence similarity. DNA-DNA hybridization (DDH) between strains with representative MLSA distances revealed a strong correlation between the data of MLSA and DDH. The 70 % DDH value for current species definition corresponds to a five-gene MLSA distance of 0.007, which could be considered as the species cut-off for the S. griseus clade. It is concluded that the MLSA procedure can be a practical, reliable and robust alternative to DDH for the identification and classification of streptomycetes at the species and intraspecies levels. Based on the data from MLSA and DDH, as well as cultural and morphological characteristics, 18 species and three subspecies of the S. griseus clade are considered to be later heterotypic synonyms of 11 genomic species: Streptomyces griseinus and Streptomyces mediolani as synonyms of Streptomyces albovinaceus; Streptomyces praecox as a synonym of Streptomyces anulatus; Streptomyces olivoviridis as a synonym of Streptomyces atroolivaceus; Streptomyces griseobrunneus as a synonym of Streptomyces bacillaris; Streptomyces cavourensis subsp. washingtonensis as a synonym of Streptomyces cyaneofuscatus; Streptomyces acrimycini, Streptomyces baarnensis, Streptomyces caviscabies and Streptomyces flavofuscus as synonyms of Streptomyces fimicarius; Streptomyces flavogriseus as a synonym of Streptomyces flavovirens; Streptomyces erumpens, 'Streptomyces ornatus' and Streptomyces setonii as synonyms of Streptomyces griseus; Streptomyces graminofaciens as a synonym of Streptomyces halstedii; Streptomyces alboviridis, Streptomyces griseus subsp. alpha, Streptomyces griseus subsp. cretosus and Streptomyces luridiscabiei as synonyms of Streptomyces microflavus; and Streptomyces californicus and Streptomyces floridae as synonyms of Streptomyces puniceus.

The symbioses between legumes and nitrogen-fixing rhizobia make the greatest contribution to the global nitrogen input via the process of biological nitrogen fixation (BNF). Bradyrhizobium stands out as the main genus nodulating basal Caesalpinioideae. We performed a polyphasic study with 11 strains isolated from root nodules of Chamaecristafasciculata, an annual multi-functional native legume of the USA. In the 16S rRNA gene phylogeny the strains were clustered in the Bradyrhizobium japonicumsuperclade. The results of analysis of the intergenic transcribed spacer (ITS) indicated less than 89.9 % similarity to other Bradyrhizobium species. Multilocus sequence analysis (MLSA) with four housekeeping genes (glnII, gyrB, recA and rpoB) confirmed the new group, sharing less than 95.2 % nucleotide identity with other species. The MLSA with 10 housekeeping genes (atpD, dnaK, gap, glnII, gltA, gyrB, pnp, recA, rpoB and thrC) indicated Bradyrhizobium daqingense as the closest species. Noteworthy, high genetic diversity among the strains was confirmed in the analyses of ITS, MLSA and BOX-PCR. Average nucleotide identity and digital DNA-DNA hybridization values were below the threshold of described Bradyrhizobium species, of 89.7 and 40 %, respectively. In the nifH and nodC phylogenies, the strains were grouped together, but with an indication of horizontal gene transfer, showing higher similarity to Bradyrhizobium arachidis and Bradyrhizobium forestalis. Other phenotypic, genotypic and symbiotic properties were evaluated, and the results altogether support the description of the CNPSo strains as representatives of the new species Bradyrhizobiumfrederickii sp. nov., with CNPSo 3426T (=USDA 10052T=U686T=CL 20T) as the type strain.

Strain aSej3T was isolated from a root nodule of a Lupinus angustifolius plant growing in Bizerte, Tunisia. 16S rRNA gene analysis placed this strain within the genus Bradyrhizobium. Multilocus sequence analysis (MLSA) including three housekeeping genes (glnII, gyrB and recA) grouped aSej3T together with Bradyrhizobium rifense CTAW71T, Bradyrhizobium cytisi CTAW11T, Bradyrhizobium ganzhouense RITF806T, Bradyrhizobium lupini USDA 3051T and Bradyrhizobium canariense BTA-1T. MLSA with five housekeeping genes (dnaK, glnII, gyrB, recA and rpoB) revealed that this strain shares less than 93.5 % nucleotide identity with other type strains. Genome sequencing and inspection revealed a genome size of 8.83 Mbp with a G+C content of 62.8 mol%. Genome-wide average nucleotide identity and digital DNA-DNA hybridization values were below 87.5 and 36.2 %, respectively, when compared to described Bradyrhizobium species. Strain aSej3T nodulated L. angustifolius plants under axenic conditions and its nodC gene clustered within the genistearum symbiovar. Altogether, the phylogenetic data and the chemotaxonomic characteristics of this strain support that aSej3T represents a new species for which we propose the name Bradyrhizobium hipponense sp. nov. with the type strain aSej3T (=DSM 108913T=LMG 31020T).

Systematics can provide a fundamental framework for understanding the relationships and diversification of organisms. Multilocus sequence analysis (MLSA) has shown great promise for an elaborate taxonomic grouping of streptomycete diversity. To evaluate the practical significance of MLSA as a valuable systematic tool for streptomycetes, we examined six endophytic Streptomyces griseus isolates and two S. griseus reference strains possessing obvious antagonistic activities and identical 16S rRNA gene sequences, using both housekeeping genes and secondary metabolic genes. All the eight strains contained PKS-I and NRPS genes, but not PKS-II genes, and showed similar diversity in both the MLSA phylogeny based on five housekeeping genes (atpD, gyrB, recA, rpoB and trpB) and fingerprinting of KS-AT genes. We also inferred a phylogeny based on concatenated amino acid sequences of representative KS-AT genes from the strains, which displayed a topology correlated well with those of housekeeping-gene MLSA and KS-AT fingerprinting. The good congruence observed between phylogenies based on the different datasets verified that the MLSA scheme provided robust resolution at intraspecific level and could predict the overall diversity of secondary metabolic potential within a Streptomyces species, despite somewhat of a discrepancy with antimicrobial data. It is therefore feasible to apply MLSA to dissecting natural diversity of streptomycetes for a better understanding of their evolution and ecology, as well as for facilitating their bioprospecting.

Members of the genus Pantoea are Gram-negative bacteria isolated from various environments. Taxonomic affiliation based on multilocus sequence analysis (MLSA) is used routinely for inferring accurate phylogeny and identification of bacterial species and genera. Partial sequences of five housekeeping genes (fusA, gyrB, leuS, rpoB, and pyrG) were extracted from 206 draft or complete genomes of Pantoea strains publicly available in databases and analyzed together with the representative sequences of the 25 validly published Pantoea type strains to verify and assess their phylogenetic assignations. Of a total of 159 strains assigned to species level, 11.3% of the non-type strains were incorrectly assigned within suitable Pantoea species. The highest proportion of misidentified strains was recorded in Pantoea vagans, 8 out of 15 (53.3%) inaccurate assignations at the species level. One probable reason for this incorrect classification could be the method previously used for strain identification. Forty-seven (22.8%) genome sequences were from strains identified at the genus level only (Pantoea sp.). A combination of MLSA, average nucleotide identities [ANI and MuMmer-based ANI (ANIm)], tetranucleotide usage pattern (TETRA), and genome-based DNA-DNA hybridization (gDDH) data was used to accurately assign 25 of the 47 strains to validly published Pantoea species, while 17 strains could be assigned as putative novel species within the genus Pantoea. Four genomes designed as Pantoea sp. were identified as Mixta calida. Positive and significant correlation coefficients were computed between MLSA and all the indices derived from whole-genome sequences being proposed for species delimitation. gDDH exhibited the best correlation with MLSA while TETRA was the worst. Accurate species-level identification is key to a better understanding of bacterial diversity and evolution. The MLSA scheme used here could be instrumental to determine the correct taxonomic status of new whole-genome sequenced Pantoea strains, especially non-type strains, before depositing into public databases.

Multi-locus Sequence Analysis

Taxonomic evaluation of the Streptomyces hygroscopicus clade using multilocus sequence analysis and DNA–DNA hybridization, validating the MLSA scheme for systematics of the whole genus

There is a need for easy, practical, reliable and robust techniques for the identification and classification of bacterial isolates to the species level as alternatives to 16S rRNA gene sequence analysis and DNA-DNA hybridization. Here, we demonstrate that multilocus sequence analysis (MLSA) of housekeeping genes is a valuable alternative technique. An MLSA study of 10 housekeeping genes (atpD, dnaK, gap, glnA, gltA, gyrB, pnp, recA, rpoB and thrC) was performed on 34 representatives of the genus Ensifer. Genetic analysis and comparison with 16S and 23S rRNA gene sequences demonstrated clear species boundaries and a higher discrimination potential for all housekeeping genes. Comparison of housekeeping gene sequence data with DNA-DNA reassociation data revealed good correlation at the intraspecies level, but indicated that housekeeping gene sequencing is superior to DNA-DNA hybridization for the assessment of genetic relatedness between Ensifer species. Our MLSA data, confirmed by DNA-DNA hybridizations, support the suggestion that Ensifer xinjiangensis is a later heterotypic synonym of Ensifer fredii.

Vigna minima is a climbing annual plant widely distributed in barren wilderness, grass land, and shrub bush of China and other countries such as Japan. However, the rhizobia nodulating with this plant has never been systematically studied. In order to reveal the biodiversity of nodulating rhizobia symbiosis with V. minima, a total of 874 rhizobium isolates were obtained from root nodules of the plant spread in 11 sampling sites of Shandong Peninsula, China, and they were designated as 41 haplotypes in the genus Bradyrhizobium based upon recA sequence analyses. By multilocus sequence analysis (MLSA) of five housekeeping genes (dnaK, glnII, gyrB, recA, and rpoB), the 41 strains representing different recA haplotypes were classified into nine defined species and nine novel genospecies. Bradyrhizobium elkanii, Bradyrhizobium ferriligni, and Bradyrhizobium pachyrhizi were the predominant and universally distributed groups. The phylogeny of symbiotic genes of nodC and nifH showed similar topology and phylogenetic relationships, in which all the representative strains were classified into two clades grouped with strains nodulating with Vigna spp., demonstrating that Vigna spp. shared common nodulating groups in the natural environment. All the representative strains formed nodules with V. minima in a nodulation test performed in green house conditions. The correlation between V. minima nodulating rhizobia and soil characteristics analyzed by CANOCO indicates that available nitrogen, total nitrogen, and organic carbon in the soil samples were the main factors affecting the distribution of rhizobia isolated in this study. This study systematically uncovered the biodiversity and distribution characteristics of V. minima nodulating rhizobia for the first time, which provided novel information for the formation of the corresponding rhizobium community.

The bacterium Pseudomonas anguilliseptica has in recent years emerged as a serious threat to production of lumpfish in Norway. Little is known about the population structure of this bacterium despite its association with disease in a wide range of different fish species throughout the world. The phylogenetic relationships between 53 isolates, primarily derived from diseased lumpfish, but including a number of reference strains from diverse geographical origins and fish species, were reconstructed by Multi-Locus Sequence Analysis (MLSA) using nine housekeeping genes (rpoB, atpD, gyrB, rpoD, ileS, aroE, carA, glnS and recA). MLSA revealed a high degree of relatedness between the studied isolates, altough the seven genotypes identified formed three main phylogenetic lineages. While four genotypes were identified amongst Norwegian lumpfish isolates, a single genotype dominated, irrespective of geographic origin. This suggests the existence of a dominant genotype associated with disease in production of lumpfish in Norwegian aquaculture. Elucidation of the population structure of the bacterium has provided valuable information for potential future vaccine development.

A novel Vibrio strain (CAIM 722T=SW9T=DSM 24596T) was isolated in 2003 from water of a shrimp (Penaeus vannamei) culture pond located in Los Mochis, Sinaloa, Mexico, and taxonomically characterized using a polyphasic approach. The 16S rRNA gene sequence clustered within those of the genus Vibrio, showing high similarity to the type strains of the Porteresiae clade. Multilocus sequence analysis using eight housekeeping genes (ftsZ, gapA, gyrB, mreB, pyrH, recA, rpoA, topA and 16S rRNA) and phylogenetic analysis with 139 single-copy genes showed that the strain forms an independent branch. Whole genome sequencing and genomic analyses (average nucleotide identity, OrthoANI, average amino acid identity and in silico DNA-DNA hybridization) produced values well below the thresholds for species delineation with all methods tested. In addition, a phenotypic characterization was performed to support the description and differentiation of the novel strain from related taxa. The results obtained demonstrate that the strain represent a novel species, for which the name Vibrio eleionomae sp. nov. is proposed.

Two bacterial strains were isolated from the hepatopancreas of a cultured shrimp (Penaeus vannamei) in Sinaloa, México. Their partial 16S rRNA gene sequences clustered within those of the genus Photobacterium, showing high similarity to the type strains of Photobacterium angustum and Photobacterium leiognathi, were 87.1% and 97.5%, respectively. Multilocus sequence analysis using eight housekeeping genes (ftsZ, gapA, gyrB, mreB, pyrH, recA, rpoA, topA and 16S rRNA) and phylogenetic analysis with 139 single-copy genes showed that the new strains form an independent branch whole genome sequencing and genomic analyses (average nucleotide identity, average amino acid identity, and in silico DNA-DNA hybridization) produced values well below the thresholds for species delineation with all methods tested. In addition, a phenotypic characterization was performed to support the description and differentiation of the novel strains from related taxa. The results obtained demonstrate that the two strains represent a novel species for which the name Photobacterium lucens sp. nov. is proposed.

Aeschynomene indica is a semiaquatic legume that forms both stem and root nodules with rhizobia. Some A. indica rhizobia (AIRs) have been reported to nodulate the host using a Nod factor-independent pathway and possess photosynthetic abilities. To investigate the diversity and community structure of AIRs in China, a total of 300 rhizobial isolates were acquired from the root and stem nodules of A. indica grown at 4 sites in Shandong Peninsula, China. Nineteen representative strains were selected according to their recA phylogeny. With further classification in comparison with reference strains, 10 Bradyrhizobium genospecies were defined based on the 16S rRNA gene phylogeny and multilocus sequence analysis (MLSA) of housekeeping genes (HKGs) recA, atpD, glnII, dnaK, gyrB, and rpoB In addition, 6 genospecies were found only in China. No nodulation gene (nodA, nodB, nodC, or nodZ) was detected in the AIRs isolates by PCR amplification and Southern blotting. Phylogenetic analysis of nifH and the photosynthesis-related gene pufLM revealed their common origins. All representative strains formed root nodules, but only 9 representative strains for 4 genospecies formed stem nodules on A. indica, indicating that the stem nodulation process of A. indica is limited to some strains. The nucleotide diversity and recombination events of the HKGs, as well as nifH and pufLM genes, showed that mutation contributes more than recombination in evolution. The distribution of dominant AIR genospecies was mainly affected by available nitrogen, organic carbon, total nitrogen, and pH. Our study helps to characterize the diversity and evolution of AIRs.IMPORTANCEAeschynomene indica rhizobia (AIRs) can form both root and stem nodules via Nod factor-independent processes, which distinguishes them from other rhizobia. This study systematically uncovered the diversity and community composition of A. indica rhizobia distributed in eastern China. Our results reclassified all the A. indica rhizobia across the world and represent a useful contribution to evaluating the diversity and distribution of the symbiont. The presence of novel genospecies specifically distributed in China enriched the A. indica rhizobia resources and provided insight into the geographic distribution of rhizobia. The phylogenetic relationship between nifH and pufLM of A. indica rhizobia across the world provides insight into the evolution of their nitrogen fixation and photosynthetic abilities.

Multilocus sequence analysis (MLSA) protocols have been developed for species circumscription for many taxa. However, at present, no studies based on MLSA have been performed within any moderately halophilic bacterial group. To test the usefulness of MLSA with these kinds of micro-organisms, the family Halomonadaceae, which includes mainly halophilic bacteria, was chosen as a model. This family comprises ten genera with validly published names and 85 species of environmental, biotechnological and clinical interest. In some cases, the phylogenetic relationships between members of this family, based on 16S rRNA gene sequence comparisons, are not clear and a deep phylogenetic analysis using several housekeeping genes seemed appropriate. Here, MLSA was applied using the 16S rRNA, 23S rRNA, atpA, gyrB, rpoD and secA genes for species of the family Halomonadaceae. Phylogenetic trees based on the individual and concatenated gene sequences revealed that the family Halomonadaceae formed a monophyletic group of micro-organisms within the order Oceanospirillales. With the exception of the genera Halomonas and Modicisalibacter, all other genera within this family were phylogenetically coherent. Five of the six studied genes (16S rRNA, 23S rRNA, gyrB, rpoD and secA) showed a consistent evolutionary history. However, the results obtained with the atpA gene were different; thus, this gene may not be considered useful as an individual gene phylogenetic marker within this family. The phylogenetic methods produced variable results, with those generated from the maximum-likelihood and neighbour-joining algorithms being more similar than those obtained by maximum-parsimony methods. Horizontal gene transfer (HGT) plays an important evolutionary role in the family Halomonadaceae; however, the impact of recombination events in the phylogenetic analysis was minimized by concatenating the six loci, which agreed with the current taxonomic scheme for this family. Finally, the findings of this study also indicated that the 16S rRNA, gyrB and rpoD genes were the most suitable genes for future taxonomic studies using MLSA within the family Halomonadaceae.