Cutaneous Mucormycosis in Tornado Survivors

A total of 18 suspected cases of cutaneous mucormycosis were identified in survivors of the May 2011 tornado in Joplin, Missouri. Apophysomyces trapeziformis was identified in 13 of the patients with microbiologically or histologically proven infection by DNA sequencing. Apophysomyces are classically associated with necrotizing skin and soft-tissue infections following traumatic inoculation of the fungal spores. Although cases of Apophysomyces infection were previously reported in survivors of tsunamis and volcanic eruptions, this is believed to be the first reported case series of cutaneous mucormycosis in tornado survivors.     

Molecular phylogenetic diversity of the emerging mucoralean fungus Apophysomyces: Proposal of three new species

BackgroundApophysomyces is a monotypic genus belonging to the order Mucorales. The species Apophysomyces elegans has been reported to cause severe infections in immunocompromised and immunocompetent people. In a previous study of Álvarez et al.³ [J Clin Microbiol 2009;47:1650–6], we demonstrated a high variability among the 5.8S rRNA gene sequences of clinical strains of A. elegans.Material and methodsWe performed a polyphasic study based on the analysis of the sequences of the histone 3 gene, the internal transcribed spacer region of the rDNA gene, and domains D1 and D2 of the 28S rRNA gene, as well as by evaluation of some relevant morphological and physiological characteristics of a set of clinical and environmental strains of A. elegans.Results and conclusionsWe have demonstrated that A. elegans is a complex of species. We propose as new species Apophysomyces ossiformis, characterised by bone-shaped sporangiospores, Apophysomyces trapeziformis, with trapezoid-shaped sporangiospores, and Apophysomyces variabilis, with variable-shaped sporangiospores. These species failed to assimilate esculin, whereas A. elegans was able to assimilate that glycoside. Amphotericin B and posaconazole are the most active in vitro drugs against Apophysomyces.     

Genome-Wide SNP-Genotyping Array to Study the Evolution of the Human Pathogen Vibrio vulnificus Biotype 3

Vibrio vulnificus is an aquatic bacterium and an important human pathogen. Strains of V. vulnificus are classified into three different biotypes. The newly emerged biotype 3 has been found to be clonal and restricted to Israel. In the family Vibrionaceae, horizontal gene transfer is the main mechanism responsible for the emergence of new pathogen groups. To better understand the evolution of the bacterium, and in particular to trace the evolution of biotype 3, we performed genome-wide SNP genotyping of 254 clinical and environmental V. vulnificus isolates with worldwide distribution recovered over a 30-year period, representing all phylogeny groups. A custom single-nucleotide polymorphism (SNP) array implemented on the Illumina GoldenGate platform was developed based on 570 SNPs randomly distributed throughout the genome. In general, the genotyping results divided the V. vulnificus species into three main phylogenetic lineages and an additional subgroup, clade B, consisting of environmental and clinical isolates from Israel. Data analysis suggested that 69% of biotype 3 SNPs are similar to SNPs from clade B, indicating that biotype 3 and clade B have a common ancestor. The rest of the biotype 3 SNPs were scattered along the biotype 3 genome, probably representing multiple chromosomal segments that may have been horizontally inserted into the clade B recipient core genome from other phylogroups or bacterial species sharing the same ecological niche. Results emphasize the continuous evolution of V. vulnificus and support the emergence of new pathogenic groups within this species as a recurrent phenomenon. Our findings contribute to a broader understanding of the evolution of this human pathogen.     

Single-Nucleotide Polymorphism Phylotyping of Escherichia coli

We describe a rapid and easily automated phylogenetic grouping technique based on analysis of bacterial genome single-nucleotide polymorphisms (SNPs). We selected 13 SNPs derived from a complete sequence analysis of 11 essential genes previously used for multilocus sequence typing (MLST) of 30 Escherichia coli strains representing the genetic diversity of the species. The 13 SNPs were localized in five genes, trpA, trpB, putP, icdA, and polB, and were selected to allow recovery of the main phylogenetic groups (groups A, B1, E, D, and B2) and subgroups of the species. In the first step, we validated the SNP approach in silico by extracting SNP data from the complete sequences of the five genes for a panel of 65 pathogenic strains belonging to different E. coli pathovars, which were previously analyzed by MLST. In the second step, we determined these SNPs by dideoxy single-base extension of unlabeled oligonucleotide primers for a collection of 183 commensal and extraintestinal clinical E. coli isolates and compared the SNP phylotyping method to previous well-established typing methods. This SNP phylotyping method proved to be consistent with the other methods for assigning phylogenetic groups to the different E. coli strains. In contrast to the other typing methods, such as multilocus enzyme electrophoresis, ribotyping, or PCR phylotyping using the presence/absence of three genomic DNA fragments, the SNP typing method described here is derived from a solid phylogenetic analysis, and the results obtained by this method are more meaningful. Our results indicate that similar approaches may be used for a wide variety of bacterial species.     

Effect of a glucose impulse on the CcpA regulon in Staphylococcus aureus

Background

A low genetic diversity in Francisella tularensis has been documented. Current DNA based genotyping methods for typing F. tularensis offer a limited and varying degree of subspecies, clade and strain level discrimination power. Whole genome sequencing is the most accurate and reliable method to identify, type and determine phylogenetic relationships among strains of a species. However, lower cost typing schemes are necessary in order to enable typing of hundreds or even thousands of isolates.

Results

We have generated a high-resolution phylogenetic tree from 40 Francisella isolates, including 13 F. tularensis subspecies holarctica (type B) strains, 26 F. tularensis subsp. tularensis (type A) strains and a single F. novicida strain. The tree was generated from global multi-strain single nucleotide polymorphism (SNP) data collected using a set of six Affymetrix GeneChip^® resequencing arrays with the non-repetitive portion of LVS (type B) as the reference sequence complemented with unique sequences of SCHU S4 (type A). Global SNP based phylogenetic clustering was able to resolve all non-related strains. The phylogenetic tree was used to guide the selection of informative SNPs specific to major nodes in the tree for development of a genotyping assay for identification of F. tularensis subspecies and clades. We designed and validated an assay that uses these SNPs to accurately genotype 39 additional F. tularensis strains as type A (A1, A2, A1a or A1b) or type B (B1 or B2).

Conclusion

Whole-genome SNP based clustering was shown to accurately identify SNPs for differentiation of F. tularensis subspecies and clades, emphasizing the potential power and utility of this methodology for selecting SNPs for typing of F. tularensis to the strain level. Additionally, whole genome sequence based SNP information gained from a representative population of strains may be used to perform evolutionary or phylogenetic comparisons of strains, or selection of unique strains for whole-genome sequencing projects.     

Comparative Genome Structure,Secondary Metabolite,and Effector Coding Capacity across Cochliobolus Pathogens

The genomes of five Cochliobolus heterostrophus strains, two Cochliobolus sativus strains, three additional Cochliobolus species (Cochliobolus victoriae, Cochliobolus carbonum, Cochliobolus miyabeanus), and closely related Setosphaeria turcica were sequenced at the Joint Genome Institute (JGI). The datasets were used to identify SNPs between strains and species, unique genomic regions, core secondary metabolism genes, and small secreted protein (SSP) candidate effector encoding genes with a view towards pinpointing structural elements and gene content associated with specificity of these closely related fungi to different cereal hosts. Whole-genome alignment shows that three to five percent of each genome differs between strains of the same species, while a quarter of each genome differs between species. On average, SNP counts among field isolates of the same C. heterostrophus species are more than 25× higher than those between inbred lines and 50× lower than SNPs between Cochliobolus species. The suites of nonribosomal peptide synthetase (NRPS), polyketide synthase (PKS), and SSP–encoding genes are astoundingly diverse among species but remarkably conserved among isolates of the same species, whether inbred or field strains, except for defining examples that map to unique genomic regions. Functional analysis of several strain-unique PKSs and NRPSs reveal a strong correlation with a role in virulence.     

Whole genome sequence analysis of Cryptococcus gattii from the Pacific Northwest reveals unexpected diversity

A recent emergence of Cryptococcus gattii in the Pacific Northwest involves strains that fall into three primarily clonal molecular subtypes: VGIIa, VGIIb and VGIIc. Multilocus sequence typing (MLST) and variable number tandem repeat analysis appear to identify little diversity within these molecular subtypes. Given the apparent expansion of these subtypes into new geographic areas and their ability to cause disease in immunocompetent individuals, differentiation of isolates belonging to these subtypes could be very important from a public health perspective. We used whole genome sequence typing (WGST) to perform fine-scale phylogenetic analysis on 20 C. gattii isolates, 18 of which are from the VGII molecular type largely responsible for the Pacific Northwest emergence. Analysis both including and excluding (289,586 SNPs and 56,845 SNPs, respectively) molecular types VGI and VGIII isolates resulted in phylogenetic reconstructions consistent, for the most part, with MLST analysis but with far greater resolution among isolates. The WGST analysis presented here resulted in identification of over 100 SNPs among eight VGIIc isolates as well as unique genotypes for each of the VGIIa, VGIIb and VGIIc isolates. Similar levels of genetic diversity were found within each of the molecular subtype isolates, despite the fact that the VGIIb clade is thought to have emerged much earlier. The analysis presented here is the first multi-genome WGST study to focus on the C. gattii molecular subtypes involved in the Pacific Northwest emergence and describes the tools that will further our understanding of this emerging pathogen.     

SYMBIODINIUM (PYRRHOPHYTA) GENOME SIZES (DNA CONTENT) ARE SMALLEST AMONG DINOFLAGELLATES1

Using flow cytometric analysis of fluorescence, we measured the genome sizes of 18 cultured “free‐living” species and 29 Symbiodinium spp. isolates cultured from stony corals, gorgonians, anemones, jellyfish, and giant clams. Genome size directly correlated with cell size, as documented previously for most eukaryotic cell lines. Among the smallest of dinoflagellates, Symbiodinium spp. (6–15 μm) possessed the lowest DNA content that we measured (1.5–4.8 pg·cell^?1). Bloom‐forming or potentially harmful species in the genera Alexandrium, Karenia, Pfiesteria, and Prorocentrum possessed genomes approximately 2 to 50 times larger in size. A phylogenetic analysis indicated that genome/cell size has apparently increased and decreased repeatedly during the evolution of dinoflagellates. In contrast, genome sizes were relatively consistent across distantly and closely related Symbiodinium spp. This may be the product of intracellular host habitats imposing strong selective pressures that have restricted symbiont size.     

Identifying genetic diversity of avirulence genes in Leptosphaeria maculans using whole genome sequencing

Next generation sequencing technology allows rapid re-sequencing of individuals, as well as the discovery of single nucleotide polymorphisms (SNPs), for genomic diversity and evolutionary analyses. By sequencing two isolates of the fungal plant pathogen Leptosphaeria maculans, the causal agent of blackleg disease in Brassica crops, we have generated a resource of over 76 million sequence reads aligned to the reference genome. We identified over 21,000 SNPs with an overall SNP frequency of one SNP every 2,065 bp. Sequence validation of a selection of these SNPs in additional isolates collected throughout Australia indicates a high degree of polymorphism in the Australian population. In preliminary phylogenetic analysis, isolates from Western Australia clustered together and those collected from Brassica juncea stubble were identical. These SNPs provide a novel marker resource to study the genetic diversity of this pathogen. We demonstrate that re-sequencing provides a method of validating previously characterised SNPs and analysing differences in important genes, such as the disease related avirulence genes of L. maculans. Understanding the genetic characteristics of this devastating pathogen is vital in developing long-term solutions to managing blackleg disease in Brassica crops.     

Highly accurate-single chromosomal complete genomes using IonTorrent and MinION sequencing of clinical pathogens

Oxford Nanopore MinION sequencing technology has been gaining immense importance in identification of pathogen and antimicrobial resistance, though with 10–15% error rate. Short read technologies generates high accurate genome but with multiple fragments of genome. This study proposes a novel workflow to reduce the indels resulted from MinION long read sequencing by overlaying short read sequences from IonTorrent in the clinical isolates. Best of both techniques were employed which generated highly accurate-single chromosomal microbial genomes with increase in completeness of genomes from 44.5%, 30% and 43% to 98.6%, 98.6% and 96.6% for P. aeruginosa, A. veronii and B. pertussis respectively. To the best of our knowledge, this is the first study to generate a hybrid of IonTorrent and MinION reads to obtain single chromosomal genomes. This would enable to precisely infer both structural arrangement of genes and SNP based analysis for phylogenetic information.     

Full genome SNP-based phylogenetic analysis reveals the origin and global spread of Brucella melitensis

Background

Brucellosis is an important zoonotic disease that affects both humans and animals. We sequenced the full genome and characterised the genetic diversity of two Brucella melitensis isolates from Malaysia and the Philippines. In addition, we performed a comparative whole-genome single nucleotide polymorphism (SNP) analysis of B. melitensis strains collected from around the world, to investigate the potential origin and the history of the global spread of B. melitensis.

Results

Single sequencing runs of each genome resulted in draft genome sequences of MY1483/09 and Phil1136/12, which covered 99.85% and 99.92% of the complete genome sequences, respectively. The B. melitensis genome sequences, and two B. abortus strains used as the outgroup strains, yielded a total of 13,728 SNP sites. Phylogenetic analysis using whole-genome SNPs and geographical distribution of the isolates revealed spatial clustering of the B. melitensis isolates into five genotypes, I, II, III, IV and V. The Mediterranean strains, identified as genotype I, occupied the basal node of the phylogenetic tree, suggesting that B. melitensis may have originated from the Mediterranean regions. All of the Asian B. melitensis strains clustered into genotype II with the SEA strains, including the two isolates sequenced in this study, forming a distinct clade denoted here as genotype IId. Genotypes III, IV and V of B. melitensis demonstrated a restricted geographical distribution, with genotype III representing the African lineage, genotype IV representing the European lineage and genotype V representing the American lineage.

Conclusion

We showed that SNPs retrieved from the B. melitensis draft full genomes were sufficient to resolve the interspecies relationships between B. melitensis strains and to discriminate between the vaccine and endemic strains. Phylogeographic reconstruction of the history of B. melitensis global spread at a finer scale by using whole-genome SNP analyses supported the origin of all B. melitensis strains from the Mediterranean region. The possible global distribution of B. melitensis following the ancient trade routes was also consistent with whole-genome SNP phylogeny. The whole genome SNP phylogenetics analysis, hence is a powerful tool for intraspecies discrimination of closely related species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1294-x) contains supplementary material, which is available to authorized users.     

Possible Origins of CTnBST, a Conjugative Transposon Found Recently in a Human Colonic Bacteroides Strain

A previous survey of Bacteroides isolates suggested that the ermB gene entered Bacteroides spp. recently. Previously, ermB had been found almost exclusively in gram-positive bacteria. In one Bacteroides strain, ermB was located on 100-kb conjugative transposon (CTn) CTnBST. To assess the possible origin of this CTn, we obtained the full DNA sequence of CTnBST and used this information to investigate its possible origins. Over one-half of CTnBST had high sequence identity to a putative CTn found in the genome of Bacteroides fragilis YCH46. This included the ends of the CTn and genes involved in integration, transfer, and excision. However, the region around the ermB gene contained genes that appeared to originate from gram-positive organisms. In particular, a 7-kb segment containing the ermB gene was 100% identical to an ermB region found in the genome of the gram-positive bacterium Arcanobacterium pyogenes. A screen of Bacteroides isolates whose DNA cross-hybridized with a CTnBST probe revealed that several isolates did not carry the 7-kb region, implying that the acquisition of this region may be more recent than the acquisition of the entire CTnBST element by Bacteroides spp. We have also identified other Bacteroides isolates that carry a slightly modified 7-kb region but have no other traces of CTnBST. Thus, it is possible that this 7-kb region could itself be part of a mobile element that has inserted in a Bacteroides CTn. Our results show that CTnBST is a hybrid element which has acquired a portion of its coding region from gram-positive bacteria but which may originally have come from Bacteroides spp. or some related species.     

<Emphasis Type="Italic">Phytophthora colocasiae</Emphasis> from Vietnam,China, Hawaii and Nepal: intra- and inter-genomic variations in ploidy and a long-lived,diploid Hawaiian lineage.

Phytophthora colocasiae is an important pathogen of taro and is widely distributed. Our goal was to develop whole genome sequence and single nucleotide polymorphism (SNP) markers to characterize historical and current populations from Hawaii (2010 and 2016, HA), historical isolates from Vietnam and China (2010, VN and CH) and current isolates from Nepal (2016, NEP). Seven isolates (VN = 2, CH = 1, HA = 1, NEP = 3) were sequenced (NCBI BioProject PRJNA378784) and compared using the reference genome of the closely related vegetable pathogen P. capsici. Genome-wide SNP analysis using 27,537 markers revealed genomes of diploid, triploid, tetraploid and higher ploidy. Ploidy varied within and between populations, with HA being primarily diploid, CH primarily triploid, VN containing diploid and triploid isolates, and NEP having predominantly triploid, tetraploid and higher ploidy. A total of 37 SNP markers were genotyped in 89 samples (grown in culture or directly from infected tissue) using targeted-sequencing. Analyses indicate a single clonal lineage dominated populations in HA from 2010 to 2016 and targeted-sequencing was useful to estimate ploidy. The implications for adaptation and evolution of P. colocasiae are discussed, as well as consequences for selection and breeding of resistant taro cultivars.     

A Real-Time PCR Array for Hierarchical Identification of Francisella Isolates

A robust, rapid and flexible real-time PCR assay for hierarchical genetic typing of clinical and environmental isolates of Francisella is presented. Typing markers were found by multiple genome and gene comparisons, from which 23 canonical single nucleotide polymorphisms (canSNPs) and 11 canonical insertion-deletion mutations (canINDELs) were selected to provide phylogenetic guidelines for classification from genus to isolate level. The specificity of the developed assay, which uses 68 wells of a 96-well real-time PCR format with a detection limit of 100 pg DNA, was assessed using 62 Francisella isolates of diverse genetic and geographical origins. It was then successfully used for typing 14 F. tularensis subsp. holarctica isolates obtained from tularemia patients in Sweden in 2008 and five more genetically diverse Francisella isolates of global origins. When applied to human ulcer specimens for direct pathogen detection the results were incomplete due to scarcity of DNA, but sufficient markers were identified to detect fine-resolution differences among F. tularensis subsp. holarctica isolates causing infection in the patients. In contrast to other real-time PCR assays for Francisella, which are typically designed for specific detection of a species, subspecies, or strain, this type of assay can be easily tailored to provide appropriate phylogenetic and/or geographical resolution to meet the objectives of the analysis.     

Genomic Comparative Study of Bovine Mastitis Escherichia coli

Escherichia coli, one of the main causative agents of bovine mastitis, is responsible for significant losses on dairy farms. In order to better understand the pathogenicity of E. coli mastitis, an accurate characterization of E. coli strains isolated from mastitis cases is required. By using phylogenetic analyses and whole genome comparison of 5 currently available mastitis E. coli genome sequences, we searched for genotypic traits specific for mastitis isolates. Our data confirm that there is a bias in the distribution of mastitis isolates in the different phylogenetic groups of the E. coli species, with the majority of strains belonging to phylogenetic groups A and B1. An interesting feature is that clustering of strains based on their accessory genome is very similar to that obtained using the core genome. This finding illustrates the fact that phenotypic properties of strains from different phylogroups are likely to be different. As a consequence, it is possible that different strategies could be used by mastitis isolates of different phylogroups to trigger mastitis. Our results indicate that mastitis E. coli isolates analyzed in this study carry very few of the virulence genes described in other pathogenic E. coli strains. A more detailed analysis of the presence/absence of genes involved in LPS synthesis, iron acquisition and type 6 secretion systems did not uncover specific properties of mastitis isolates. Altogether, these results indicate that mastitis E. coli isolates are rather characterized by a lack of bona fide currently described virulence genes.     

An evaluation of the species and subspecies of the genus Salmonella with whole genome sequence data: Proposal of type strains and epithets for novel S. enterica subspecies VII,VIII, IX,X and XI

Species and subspecies within the Salmonella genus have been defined for public health purposes by biochemical properties; however, reference laboratories have increasingly adopted sequence-based, and especially whole genome sequence (WGS), methods for surveillance and routine identification. This leads to potential disparities in subspecies definitions, routine typing, and the ability to detect novel subspecies. A large-scale analysis of WGS data from the routine sequencing of clinical isolates was employed to define and characterise Salmonella subspecies population structure, demonstrating that the Salmonella species and subspecies were genetically distinct, including those previously identified through phylogenetic approaches, namely: S. enterica subspecies londinensis (VII), subspecies brasiliensis (VIII), subspecies hibernicus (IX) and subspecies essexiensis (X). The analysis also identified an additional novel subspecies, reptilium (XI). Further, these analyses indicated that S. enterica subspecies arizonae (IIIa) isolates were divergent from the other S. enterica subspecies, which clustered together and, on the basis of ANI analysis, subspecies IIIa was sufficiently distinct to be classified as a separate species, S. arizonae. Multiple phylogenetic and statistical approaches generated congruent results, suggesting that the proposed species and subspecies structure was sufficiently biologically robust for routine application. Biochemical analyses demonstrated that not all subspecies were distinguishable by these means and that biochemical approaches did not capture the genomic diversity of the genus. We recommend the adoption of standardised genomic definitions of species and subspecies and a genome sequence-based approach to routine typing for the identification and definition of novel subspecies.     

A Brevibacillus sp. antagonistic to mycotoxigenic Fusarium spp.

Antagonistic microbes were isolated from soils to control mycotoxin contamination of cereals by limiting the growth of mycotoxigenic Fusarium species. In total, 341 bacterial isolates were examined for antifungal activity against eight mycotoxigenic Fusarium species using dual culture assays. The screening identified 11 isolates that inhibited mycelial growth of all Fusarium species tested. The culture filtrates of 2 of the 11 isolates completely inhibited germination of conidia up to 21 days of incubation. These two isolates exhibited identical activity toward the fungi tested and were identified as Brevibacillus spp. based on 16S rRNA sequence homology. The most closely related species based on phylogenetic analysis was Brevibacillus reuszeri. Additional dual culturing using further fungal species showed that the antagonistic Brevibacillus inhibited the growth of most Fusarium species tested (39 of 46 species), two Epicoccum spp., one Alternaria sp., three Aspergillus spp. (3 of 11), and three Penicillium spp. (3 of 8). The in vivo assay was performed to test the efficacy of antagonistic Brevibacillus isolates on maize ears and revealed that the application of microbes suppressed ear rot (ANOVA, p = 0.0020). This Brevibacillus sp. may be an antagonist of the majority of Fusarium species, including mycotoxigenic species.     

Using Neisseria meningitidis genomic diversity to inform outbreak strain identification

Meningococcal disease is a life-threatening illness caused by the human-restricted bacterium Neisseria meningitidis. Outbreaks in the USA involve at least two cases in an organization or community caused by the same serogroup within three months. Genome comparisons, including phylogenetic analysis and quantification of genome distances can provide confirmatory evidence of pathogen transmission during an outbreak. Interpreting genome distances depends on understanding their distribution both among isolates from outbreaks and among those not from outbreaks. Here, we identify outbreak strains based on phylogenetic relationships among 141 N. meningitidis isolates collected from 28 outbreaks in the USA during 2010–2017 and 1516 non-outbreak isolates collected through contemporaneous meningococcal surveillance. We show that genome distance thresholds based on the maximum SNPs and allele distances among isolates in the phylogenetically defined outbreak strains are sufficient to separate most pairs of non-outbreak isolates into separate strains. Non-outbreak isolate pairs that could not be distinguished from each other based on genetic distances were concentrated in the clonal complexes CC11, CC103, and CC32. Within each of these clonal complexes, phylodynamic analysis identified a group of isolates with extremely low diversity, collected over several years and multiple states. Clusters of isolates with low genetic diversity could indicate increased pathogen transmission, potentially resulting in local outbreaks or nationwide clonal expansions.     

Evidence for the existence of two new members of the family Chlamydiaceae and proposal of Chlamydia avium sp. nov. and Chlamydia gallinacea sp. nov.

The family Chlamydiaceae with the recombined single genus Chlamydia currently comprises nine species, all of which are obligate intracellular organisms distinguished by a unique biphasic developmental cycle. Anecdotal evidence from epidemiological surveys in flocks of poultry, pigeons and psittacine birds have indicated the presence of non-classified chlamydial strains, some of which may act as pathogens. In the present study, phylogenetic analysis of ribosomal RNA and ompA genes, as well as multi-locus sequence analysis of 11 field isolates were conducted. All independent analyses assigned the strains into two different clades of monophyletic origin corresponding to pigeon and psittacine strains or poultry isolates, respectively. Comparative genome analysis involving the type strains of currently accepted Chlamydiaceae species and the designated type strains representing the two new clades confirmed that the latter could be classified into two different species as their average nucleotide identity (ANI) values were always below 94%, both with the closest relative species and between themselves.     

Diversity of Pseudomonas Genomes,Including Populus-Associated Isolates,as Revealed by Comparative Genome Analysis

The Pseudomonas genus contains a metabolically versatile group of organisms that are known to occupy numerous ecological niches, including the rhizosphere and endosphere of many plants. Their diversity influences the phylogenetic diversity and heterogeneity of these communities. On the basis of average amino acid identity, comparative genome analysis of >1,000 Pseudomonas genomes, including 21 Pseudomonas strains isolated from the roots of native Populus deltoides (eastern cottonwood) trees resulted in consistent and robust genomic clusters with phylogenetic homogeneity. All Pseudomonas aeruginosa genomes clustered together, and these were clearly distinct from other Pseudomonas species groups on the basis of pangenome and core genome analyses. In contrast, the genomes of Pseudomonas fluorescens were organized into 20 distinct genomic clusters, representing enormous diversity and heterogeneity. Most of our 21 Populus-associated isolates formed three distinct subgroups within the major P. fluorescens group, supported by pathway profile analysis, while two isolates were more closely related to Pseudomonas chlororaphis and Pseudomonas putida. Genes specific to Populus-associated subgroups were identified. Genes specific to subgroup 1 include several sensory systems that act in two-component signal transduction, a TonB-dependent receptor, and a phosphorelay sensor. Genes specific to subgroup 2 contain hypothetical genes, and genes specific to subgroup 3 were annotated with hydrolase activity. This study justifies the need to sequence multiple isolates, especially from P. fluorescens, which displays the most genetic variation, in order to study functional capabilities from a pangenomic perspective. This information will prove useful when choosing Pseudomonas strains for use to promote growth and increase disease resistance in plants.