Single-copy Core Genes Research Articles

Genomic insight into the origins and evolution of symbiosis genes in Phaseolus vulgaris microsymbionts

BackgroundPhaseolus vulgaris (common bean) microsymbionts belonging to the bacterial genera Rhizobium, Bradyrhizobium, and Ensifer (Sinorhizobium) have been isolated across the globe. Individual symbiosis genes (e.g., nodC) of these rhizobia can be different within each genus and among distinct genera. Little information is available about the symbiotic structure of indigenous Rhizobium strains nodulating introduced bean plants or the emergence of a symbiotic ability to associate with bean plants in Bradyrhizobium and Ensifer strains. Here, we sequenced the genomes of 29 representative bean microsymbionts (21 Rhizobium, four Ensifer, and four Bradyrhizobium) and compared them with closely related reference strains to estimate the origins of symbiosis genes among these Chinese bean microsymbionts.ResultsComparative genomics demonstrated horizontal gene transfer exclusively at the plasmid level, leading to expanded diversity of bean-nodulating Rhizobium strains. Analysis of vertically transferred genes uncovered 191 (out of the 2654) single-copy core genes with phylogenies strictly consistent with the taxonomic status of bacterial species, but none were found on symbiosis plasmids. A common symbiotic region was wholly conserved within the Rhizobium genus yet different from those of the other two genera. A single strain of Ensifer and two Bradyrhizobium strains shared similar gene content with soybean microsymbionts in both chromosomes and symbiotic regions.ConclusionsThe 19 native bean Rhizobium microsymbionts were assigned to four defined species and six putative novel species. The symbiosis genes of R. phaseoli, R. sophoriradicis, and R. esperanzae strains that originated from Mexican bean-nodulating strains were possibly introduced alongside bean seeds. R. anhuiense strains displayed distinct host ranges, indicating transition into bean microsymbionts. Among the six putative novel species exclusive to China, horizontal transfer of symbiosis genes suggested symbiosis with other indigenous legumes and loss of originally symbiotic regions or non-symbionts before the introduction of common bean into China. Genome data for Ensifer and Bradyrhizobium strains indicated symbiotic compatibility between microsymbionts of common bean and other hosts such as soybean.

Genomic and phenotypic characterization of Burkholderia isolates from the potable water system of the International Space Station.

The opportunistic pathogens Burkholderia cepacia and Burkholderia contaminans, both genomovars of the Burkholderia cepacia complex (BCC), are frequently cultured from the potable water dispenser (PWD) of the International Space Station (ISS). Here, we sequenced the genomes and conducted phenotypic assays to characterize these Burkholderia isolates. All recovered isolates of the two species fall within monophyletic clades based on phylogenomic trees of conserved single-copy core genes. Within species, the ISS-derived isolates all demonstrate greater than 99% average nucleotide identity (with 95–99% of genomes aligning) and share around 90% of the identified gene clusters from a pangenomic analysis–suggesting that the two groups are each composed of highly similar genomic lineages and their members may have all stemmed from the same two founding populations. The differences that can be observed between the recovered isolates at the pangenomic level are primarily located within putative plasmids. Phenotypically, macrophage intracellularization and lysis occurred at generally similar rates between all ISS-derived isolates, as well as with their respective type-terrestrial strain references. All ISS-derived isolates exhibited antibiotic sensitivity similar to that of the terrestrial reference strains, and minimal differences between isolates were observed. With a few exceptions, biofilm formation rates were generally consistent across each species. And lastly, though isolation date does not necessarily provide any insight into how long a given isolate had been aboard the ISS, none of the assayed physiology correlated with either date of isolation or distances based on nucleotide variation. Overall, we find that while the populations of Burkholderia present in the ISS PWS each maintain virulence, they are likely are not more virulent than those that might be encountered on planet and remain susceptible to clinically used antibiotics.

Comparative Genomics of Pandoraea, a Genus Enriched in Xenobiotic Biodegradation and Metabolism.

Comparative analysis of partial gyrB, recA, and gltB gene sequences of 84 Pandoraea reference strains and field isolates revealed several clusters that included no taxonomic reference strains. The gyrB, recA, and gltB phylogenetic trees were used to select 27 strains for whole-genome sequence analysis and for a comparative genomics study that also included 41 publicly available Pandoraea genome sequences. The phylogenomic analyses included a Genome BLAST Distance Phylogeny approach to calculate pairwise digital DNA–DNA hybridization values and their confidence intervals, average nucleotide identity analyses using the OrthoANIu algorithm, and a whole-genome phylogeny reconstruction based on 107 single-copy core genes using bcgTree. These analyses, along with subsequent chemotaxonomic and traditional phenotypic analyses, revealed the presence of 17 novel Pandoraea species among the strains analyzed, and allowed the identification of several unclassified Pandoraea strains reported in the literature. The genus Pandoraea has an open pan genome that includes many orthogroups in the ‘Xenobiotics biodegradation and metabolism’ KEGG pathway, which likely explains the enrichment of these species in polluted soils and participation in the biodegradation of complex organic substances. We propose to formally classify the 17 novel Pandoraea species as P. anapnoica sp. nov. (type strain LMG 31117T = CCUG 73385T), P. anhela sp. nov. (type strain LMG 31108T = CCUG 73386T), P. aquatica sp. nov. (type strain LMG 31011T = CCUG 73384T), P. bronchicola sp. nov. (type strain LMG 20603T = ATCC BAA-110T), P. capi sp. nov. (type strain LMG 20602T = ATCC BAA-109T), P. captiosa sp. nov. (type strain LMG 31118T = CCUG 73387T), P. cepalis sp. nov. (type strain LMG 31106T = CCUG 39680T), P. commovens sp. nov. (type strain LMG 31010T = CCUG 73378T), P. communis sp. nov. (type strain LMG 31110T = CCUG 73383T), P. eparura sp. nov. (type strain LMG 31012T = CCUG 73380T), P. horticolens sp. nov. (type strain LMG 31112T = CCUG 73379T), P. iniqua sp. nov. (type strain LMG 31009T = CCUG 73377T), P. morbifera sp. nov. (type strain LMG 31116T = CCUG 73389T), P. nosoerga sp. nov. (type strain LMG 31109T = CCUG 73390T), P. pneumonica sp. nov. (type strain LMG 31114T = CCUG 73388T), P. soli sp. nov. (type strain LMG 31014T = CCUG 73382T), and P. terrigena sp. nov. (type strain LMG 31013T = CCUG 73381T).

Comparative genomic analysis of the Hafnia genus reveals an explicit evolutionary relationship between the species alvei and paralvei and provides insights into pathogenicity

BackgroundThe Hafnia genus is an opportunistic pathogen that has been implicated in both nosocomial and community-acquired infections. Although Hafnia is fairly often isolated from clinical material, its taxonomy has remained an unsolved riddle, and the involvement and importance of Hafnia in human disease is also uncertain. Here, we used comparative genomic analysis to define the taxonomy of Hafnia, identify species-specific genes that may be the result of ecological and pathogenic specialization, and reveal virulence-related genetic profiles that may contribute to pathogenesis.ResultsOne complete genome sequence and 19 draft genome sequences for Hafnia strains were generated and combined with 27 publicly available genomes. We provided high-resolution typing methods by constructing phylogeny and population structure based on single-copy core genes in combination with whole genome average nucleotide identity to identify two distant Hafnia species (alvei and paralvei) and one mislabeled strain. The open pan-genome and the presence of numerous mobile genetic elements reveal that Hafnia has undergone massive gene rearrangements. Presence of species-specific core genomes associated with metabolism and transport suggests the putative niche differentiation between alvei and paralvei. We also identified possession of diverse virulence-related profiles in both Hafnia species., including the macromolecular secretion system, virulence, and antimicrobial resistance. In the macromolecular system, T1SS, Flagellum 1, Tad pilus and T6SS-1 were conserved in Hafnia, whereas T4SS, T5SS, and other T6SSs exhibited the evolution of diversity. The virulence factors in Hafnia are related to adherence, toxin, iron uptake, stress adaptation, and efflux pump. The identified resistance genes are associated with aminoglycoside, beta-lactam, bacitracin, cationic antimicrobial peptide, fluoroquinolone, and rifampin. These virulence-related profiles identified at the genomic level provide insights into Hafnia pathogenesis and the differentiation between alvei and paralvei.ConclusionsOur research using core genome phylogeny and comparative genomics analysis of a larger collection of strains provides a comprehensive view of the taxonomy and species-specific traits between Hafnia species. Deciphering the genome of Hafnia strains possessing a reservoir of macromolecular secretion systems, virulence factors, and resistance genes related to pathogenicity may provide insights into addressing its numerous infections and devising strategies to combat the pathogen.

Genome and transcriptome characterization of the glycoengineered Nicotiana benthamiana line \u0394XT/FT

BackgroundThe allotetraploid tobacco species Nicotiana benthamiana native to Australia has become a popular host for recombinant protein production. Although its usage grows every year, little is known on this plant’s genomic and transcriptomic features. Most N. benthamiana accessions currently used in research lack proper documentation of their breeding history and provenance. One of these, the glycoengineered N. benthamiana line ΔXT/FT is increasingly used for the production of biopharmaceutical proteins.ResultsBased on an existing draft assembly of the N. benthamiana genome we predict 50,516 protein –encoding genes (62,216 transcripts) supported by expression data derived from 2.35 billion mRNA-seq reads. Using single-copy core genes we show high completeness of the predicted gene set. We functionally annotate more than two thirds of the gene set through sequence homology to genes from other Nicotiana species. We demonstrate that the expression profiles from leaf tissue of ΔXT/FT and its wild type progenitor only show minimal differences. We identify the transgene insertion sites in ΔXT/FT and show that one of the transgenes was inserted inside another predicted gene that most likely lost its function upon insertion. Based on publicly available mRNA-seq data, we confirm that the N. benthamiana accessions used by different research institutions most likely derive from a single source.ConclusionsThis work provides gene annotation of the N. benthamiana genome, a genomic and transcriptomic characterization of a transgenic N. benthamiana line in comparison to its wild-type progenitor, and sheds light onto the relatedness of N. benthamiana accessions that are used in laboratories around the world.

Characterization of novel Gluconobacter species from fruits and fermented food products: Gluconobacter cadivus sp. nov., Gluconobacter vitians sp. nov. and Gluconobacter potus sp. nov.

Strains LMG 1744T, LMG 1745, LMG 31484T, LMG 1764T and R-71646 were isolated from rotting fruits and fermented food products. A phylogenomic analysis based on 107 single-copy core genes revealed that they grouped in a Gluconobacter lineage comprising Gluconobacter oxydans, Gluconobacter roseus, Gluconobacter sphaericus, Gluconobacter kanchanaburiensis, Gluconobacter albidus, Gluconobacter cerevisiae, Gluconobacter kondonii and Gluconobacter aidae. OrthoANIu and digital DNA hybridization analyses demonstrated that these five strains represented three novel Gluconobacter species, which could be differentiated from the type strains of closely related Gluconobacter species by multiple phenotypic characteristics. We therefore propose to classify strains LMG 1744T and LMG 1745 in the novel species Gluconobacter cadivus sp. nov., with LMG 1744T (=CECT 30141T) as the type strain; to classify strain LMG 31484T as the novel species Gluconobacter vitians sp. nov., with LMG 31484T (=CECT 30132T) as the type strain; and to classify strains LMG 1764T and R-71646 in the novel species Gluconobacter potus sp. nov., with LMG 1764T (=CECT 30140T) as the type strain.

Pan-Genomic Study of Mycobacterium tuberculosis Reflecting the Primary/Secondary Genes, Generality/Individuality, and the Interconversion Through Copy Number Variations.

Tuberculosis (TB) has surpassed HIV as the leading infectious disease killer worldwide since 2014. The main pathogen, Mycobacterium tuberculosis (Mtb), contains ~4,000 genes that account for ~90% of the genome. However, it is still unclear which of these genes are primary/secondary, which are responsible for generality/individuality, and which interconvert during evolution. Here we utilized a pan-genomic analysis of 36 Mtb genomes to address these questions. We identified 3,679 Mtb core (i.e., primary) genes, determining their phenotypic generality (e.g., virulence, slow growth, dormancy). We also observed 1,122 dispensable and 964 strain-specific secondary genes, reflecting partially shared and lineage-/strain-specific individualities. Among which, five L2 lineage-specific genes might be related to the increased virulence of the L2 lineage. Notably, we discovered 28 Mtb “Super Core Genes” (SCGs: more than a copy in at least 90% strains), which might be of increased importance, and reflected the “super phenotype generality.” Most SCGs encode PE/PPE, virulence factors, antigens, and transposases, and have been verified as playing crucial roles in Mtb pathogenicity. Further investigation of the 28 SCGs demonstrated the interconversion among SCGs, single-copy core, dispensable, and strain-specific genes through copy number variations (CNVs) during evolution; different mutations on different copies highlight the delicate adaptive-evolution regulation amongst Mtb lineages. This reflects that the importance of genes varied through CNVs, which might be driven by selective pressure from environment/host-adaptation. In addition, compared with Mycobacterium bovis (Mbo), Mtb possesses 48 specific single core genes that partially reflect the differences between Mtb and Mbo individuality.

Genetic and Phylogenetic Characteristics of Pasteurella multocida Isolates From Different Host Species.

Pasteurella multocida is a leading cause of respiratory diseases in many host species. To understand the genetic characteristics of P. multocida strains isolated from different host species, we sequenced the genomic DNA of P. multocida isolated from pigs and analyzed the genetic characteristics of strains from avian species, bovine species, pigs, and rabbits using whole genome sequence (WGS) data. Our results found that a capsular: lipopolysaccharide (LPS): multilocus sequence typing (MLST) genotype A: L1: ST129 (43.75%) was predominant in avian P. multocida; while genotypes B: L2: ST122 (60.00%) and A: L3: ST79 (30.00%) were predominate in bovine P. multocida; genotype D: L6: ST50 (37.50%) in porcine P. multocida; and genotype A: L3: ST9 (76.47%) in rabbit P. multocida. Comparative genomic analysis of P. multocida from different host species found that there are no genes in the P. multocida genome that are specific to any type of host. Phylogenetic analysis using either whole-genome single nucleotide polymorphisms (SNPs) or the set of SNPs present in all single-copy core genes across genomes showed that P. multocida strains with the same LPS genotype and MLST genotype were clustered together, suggesting the combining both the LPS and MLST typing schemes better explained the topology seen in the P. multocida phylogeny.

Burkholderia insecticola sp. nov., a gut symbiotic bacterium of the bean bug Riptortus pedestris.

A Gram-negative, aerobic, rod-shaped, non-spore-forming, motile bacterium, designated strain RPE64T, was isolated from the gut symbiotic organ of the bean bug Riptortus pedestris, collected in Tsukuba, Japan, in 2007. 16S rRNA gene sequencing showed that this strain belongs to the Burkholderia glathei clade, exhibiting the highest sequence similarity to Burkholderia peredens LMG 29314T (100 %), Burkholderia turbans LMG 29316T (99.52 %) and Burkholderia ptereochthonis LMG 29326T (99.04 %). Phylogenomic analyses based on 107 single-copy core genes and Genome blast Distance Phylogeny confirmed B. peredens LMG 29314T, B. ptereochthonis LMG 29326T and several uncultivated, endophytic Burkholderia species as its nearest phylogenetic neighbours. Digital DNA-DNA hybridization experiments unambiguously demonstrated that strain RPE64T represents a novel species in this lineage. The G+C content of its genome was 63.2 mol%. The isoprenoid quinone was ubiquinone 8 and the predominant fatty acid components were C16 : 0, C18 : 1ω7c and C17 : 0 cyclo. The absence of nitrate reduction and the capacity to grow at pH 8 clearly differentiated strain RPE64T from related Burkholderia species. Based on these genotypic and phenotypic characteristics, strain RPE64T is classified as representing a novel species of the genus Burkholderia, for which the name Burkholderia insecticola sp. nov. is proposed. The type strain is RPE64T (=NCIMB 15023T=JCM 31142T).

A comprehensive understanding of the biocontrol potential of Bacillus velezensis LM2303 against Fusarium head blight.

Fusarium head blight (FHB) mainly caused by F. graminearum, always brings serious damage to wheat production worldwide. In this study, we found that strain LM2303 had strong antagonist activity against F. graminearum and significantly reduced disease severity of FHB with the control efficiency of 72.3% under field conditions. To gain a comprehensive understanding of the biocontrol potential of strain LM2303 against FHB, an integrated approach of genome mining and chemical analysis was employed. The whole genome of strain LM2303 was obtained and analyzed, showing the largest number of genes/gene clusters associated with biocontrol functions as compared with the known biocontrol strains (FZB42, M75, CAU B946). And strain LM2303 was accurately determined as a member of the B. velezensis clade using the phylogenomic analysis of single-copy core genes. Through genome mining, 13 biosynthetic gene clusters(BGCs) encoding secondary metabolites with biocontrol functions were identified, which were further confirmed through chemical analyses such as UHPLC-ESI-MS, including three antifungal metabolites (fengycin B, iturin A, and surfactin A), eight antibacterial metabolites (surfactin A, butirosin, plantazolicin and hydrolyzed plantazolicin, kijanimicin, bacilysin, difficidin, bacillaene A and bacillaene B, 7-o-malonyl macrolactin A and 7-o-succinyl macrolactin A), the siderophore bacillibactin, molybdenum cofactor and teichuronic acid. In addition, genes/gene clusters involved in plant colonization, plant growth promotion and induced systemic resistance were also found and analyzed, along with the corresponding metabolites. Finally, four different mechanisms of strain LM2303 involved in the biocontrol of FHB were putatively obtained. This work provides better insights into a mechanistic understanding of strain LM2303 in control of FHB, reinforcing the higher potential of this strain as a powerful biocontrol strain agent (BCA) for FHB control. The results also provide scientific reference and comparison for other biocontrol strains.

Emended description of Salinivibrio proteolyticus, including Salinivibrio costicola subsp. vallismortis and five new isolates.

We carried out a comparative taxonomic study of Salinivibrio proteolyticus and Salinivibrio costicola subsp. vallismortis, as well as of five halophilic strains (IB574, IB872, PR5, PR919 and PR932), isolated from salterns in Spain and Puerto Rico that were closely related to these bacteria. Multilocus sequence analysis of concatenated gyrB, recA, rpoA and rpoD housekeeping genes showed that they constituted a single cluster separate from the other species and subspecies of Salinivibrio. Experimental and in silico DNA-DNA hybridization studies indicated that they are members of the same species, with relatedness of 100-74 % and 97.8-70.0 %, respectively. The average nucleotide identity (ANI) determined for these strains was 99.7-95.6 % for ANIb and 99.7-95.7 % for OrthoANI. However, the ANI values for S. costicolasubsp.vallismortis DSM 8285T with respect to S. costicolasubsp.costicola DSM 11403T and S. costicolasubsp.alcaliphilus DSM 16359T were 78.7 and 78.9 % (ANIb) and 79.4 and 79.4 % (OrthoANI), respectively. The phylogenomic tree based on 1072 concatenated orthologous single-copy core genes confirmed that S. proteolyticus, S. costicolasubsp.vallismortis and the five new isolates constitute a coherent single phylogroup, separated from the other species and subspecies of Salinivibrio. All these data indicate that S. costicolasubsp.vallismortis is a heterotypic synonym of S. proteolyticus and we propose an emended description of this species.

Genomic insight into the taxonomy of Rhizobium genospecies that nodulate Phaseolus vulgaris

Due to the wide cultivation of bean (Phaseolus vulgaris L.), rhizobia associated with this plant have been isolated from many different geographical regions. In order to investigate the species diversity of bean rhizobia, comparative genome sequence analysis was performed in the present study for 69 Rhizobium strains mainly isolated from root nodules of bean and clover (Trifolium spp.). Based on genome average nucleotide identity, digital DNA:DNA hybridization, and phylogenetic analysis of 1,458 single-copy core genes, these strains were classified into 28 clusters, consistent with their species definition based on multilocus sequence analysis (MLSA) of atpD, glnII, and recA. The bean rhizobia were found in 16 defined species and nine putative novel species; in addition, 35 strains previously described as Rhizobium etli, Rhizobium phaseoli, Rhizobium vallis, Rhizobium gallicum, Rhizobium leguminosarum and Rhizobium spp. should be renamed. The phylogenetic patterns of symbiotic genes nodC and nifH were highly host-specific and inconsistent with the genomic phylogeny. Multiple symbiovars (sv.) within the Rhizobium species were found as a common feature: sv. phaseoli, sv. trifolii and sv. viciae in Rhizobium anhuiense; sv. phaseoli and sv. mimosae in Rhizobium sophoriradicis/R. etli/Rhizobium sp. III; sv. phaseoli and sv. trifolii in Rhizobium hidalgonense/Rhizobium acidisoli; sv. phaseoli and sv. viciae in R. leguminosarum/Rhizobium sp. IX; sv. trifolii and sv. viciae in Rhizobium laguerreae. Thus, genomic comparison revealed great species diversity in bean rhizobia, corrected the species definition of some previously misnamed strains, and demonstrated the MLSA a valuable and simple method for defining Rhizobium species.

Linking pangenomes and metagenomes: the Prochlorococcus metapangenome.

Pangenomes offer detailed characterizations of core and accessory genes found in a set of closely related microbial genomes, generally by clustering genes based on sequence homology. In comparison, metagenomes facilitate highly resolved investigations of the relative distribution of microbial genomes and individual genes across environments through read recruitment analyses. Combining these complementary approaches can yield unique insights into the functional basis of microbial niche partitioning and fitness, however, advanced software solutions are lacking. Here we present an integrated analysis and visualization strategy that provides an interactive and reproducible framework to generate pangenomes and to study them in conjunction with metagenomes. To investigate its utility, we applied this strategy to a Prochlorococcus pangenome in the context of a large-scale marine metagenomic survey. The resulting Prochlorococcus metapangenome revealed remarkable differential abundance patterns between very closely related isolates that belonged to the same phylogenetic cluster and that differed by only a small number of gene clusters in the pangenome. While the relationships between these genomes based on gene clusters correlated with their environmental distribution patterns, phylogenetic analyses using marker genes or concatenated single-copy core genes did not recapitulate these patterns. The metapangenome also revealed a small set of core genes that mostly occurred in hypervariable genomic islands of the Prochlorococcus populations, which systematically lacked read recruitment from surface ocean metagenomes. Notably, these core gene clusters were all linked to sugar metabolism, suggesting potential benefits to Prochlorococcus from a high sequence diversity of sugar metabolism genes. The rapidly growing number of microbial genomes and increasing availability of environmental metagenomes provide new opportunities to investigate the functioning and the ecology of microbial populations, and metapangenomes can provide unique insights for any taxon and biome for which genomic and sufficiently deep metagenomic data are available.

Leptospira venezuelensis sp. nov., a new member of the intermediate group isolated from rodents, cattle and humans.

Three strains, CLM-U50T, CLM-R50 and IVIC-Bov1, belonging to the genus Leptospira, were isolated in Venezuela from a patient with leptospirosis, a domestic rat (Rattus norvegicus) and a cow (Bos taurus), respectively. The initial characterisation of these strains based on the rrs gene (16S rRNA) suggested their designation as a novel species within the 'intermediates' group of the genus Leptospira. Further phylogenomic characterisation based on single copy core genes was consistent with their separation into a novel species. The average nucleotide identity between these three strains was >99 %, but below 89 % with respect to any previously described leptospiral species, also supporting their designation as a novel species. Given this evidence, these three isolates were considered to represent a novel species, for which the name Leptospiravenezuelensis sp. nov. is proposed, with CLM-U50T (=CIP 111407T=DSM 105752T) as the type strain.

Genome-Wide Analysis of the First Sequenced Mycoplasma capricolum subsp. capripneumoniae Strain M1601.

Mycoplasma capricolum subsp. capripneumoniae (Mccp) is a common pathogen of goats that causes contagious caprine pleuropneumonia. We closed the gap and corrected rRNA operons in the draft genome of Mccp M1601: a strain isolated from an infected goat in a farm in Gansu, China. The genome size of M1601 is 1,016,707 bp with a GC content of 23.67%. We identified 915 genes (occupying 90.27% of the genome), of which 713 are protein-coding genes (excluding 163 pseudogenes). No genomic islands and complete insertion sequences were found in the genome. Putative determinants associated with the organism’s virulence were analyzed, and 26 genes (including one adhesion protein gene, two capsule synthesis gene clusters, two lipoproteins, hemolysin A, ClpB, and proteins involved in pyruvate metabolism and cation transport) were potential virulence factors. In addition, two transporter systems (ATP-binding cassette [ABC] transporters and phosphotransferase) and two secretion systems (Sec and signal recognition particle [SRP] pathways) were observed in the Mccp genome. Genome synteny analysis reveals a good collinear relationship between M1601 and Mccp type strain F38. Phylogenetic analysis based on 11 single-copy core genes of 31 Mycoplasma strains revealed good collinearity between M1601 and Mycoplasma capricolum subsp. capricolum (Mcc) and close relationship among Mycoplasma mycoides cluster strains. Our genome-wide analysis of Mccp M1601 provides helpful information on the pathogenic mechanisms and genetics of Mccp.

BcgTree: automatized phylogenetic tree building from bacterial core genomes.

The need for multi-gene analyses in scientific fields such as phylogenetics and DNA barcoding has increased in recent years. In particular, these approaches are increasingly important for differentiating bacterial species, where reliance on the standard 16S rDNA marker can result in poor resolution. Additionally, the assembly of bacterial genomes has become a standard task due to advances in next-generation sequencing technologies. We created a bioinformatic pipeline, bcgTree, which uses assembled bacterial genomes either from databases or own sequencing results from the user to reconstruct their phylogenetic history. The pipeline automatically extracts 107 essential single-copy core genes, found in a majority of bacteria, using hidden Markov models and performs a partitioned maximum-likelihood analysis. Here, we describe the workflow of bcgTree and, as a proof-of-concept, its usefulness in resolving the phylogeny of 293 publically available bacterial strains of the genus Lactobacillus. We also evaluate its performance in both low- and high-level taxonomy test sets. The tool is freely available at github ( https://github.com/iimog/bcgTree ) and our institutional homepage ( http://www.dna-analytics.biozentrum.uni-wuerzburg.de ).

The Dynamics of Genetic Interactions between Vibrio metoecus and Vibrio cholerae, Two Close Relatives Co-Occurring in the Environment.

Vibrio metoecus is the closest relative of Vibrio cholerae, the causative agent of the potent diarrheal disease cholera. Although the pathogenic potential of this new species is yet to be studied in depth, it has been co-isolated with V. cholerae in coastal waters and found in clinical specimens in the United States. We used these two organisms to investigate the genetic interaction between closely related species in their natural environment. The genomes of 20 V. cholerae and 4 V. metoecus strains isolated from a brackish coastal pond on the US east coast, as well as 4 clinical V. metoecus strains were sequenced and compared with reference strains. Whole genome comparison shows 86–87% average nucleotide identity (ANI) in their core genes between the two species. On the other hand, the chromosomal integron, which occupies approximately 3% of their genomes, shows higher conservation in ANI between species than any other region of their genomes. The ANI of 93–94% observed in this region is not significantly greater within than between species, meaning that it does not follow species boundaries. Vibrio metoecus does not encode toxigenic V. cholerae major virulence factors, the cholera toxin and toxin-coregulated pilus. However, some of the pathogenicity islands found in pandemic V. cholerae were either present in the common ancestor it shares with V. metoecus, or acquired by clinical and environmental V. metoecus in partial fragments. The virulence factors of V. cholerae are therefore both more ancient and more widespread than previously believed. There is high interspecies recombination in the core genome, which has been detected in 24% of the single-copy core genes, including genes involved in pathogenicity. Vibrio metoecus was six times more often the recipient of DNA from V. cholerae as it was the donor, indicating a strong bias in the direction of gene transfer in the environment.

Bacillus Cereus: A Leukemia-Specific, Neuroinvasive Pathogen?

BackgroundBacillus cereus has been associated with central nervous system (CNS) infections in immunocompromised patients.1 We report 5 cases of neuroinvasive B.cereus infections (4 fatal) in patients with acute myelogenous leukemia (AML) at Brigham and Women’s Hospital (BWH) in Boston from May 2013 to February 2014.MethodsClinical Investigation: Cases were defined as patients with positive blood or CNS cultures for B.cereus and/or positive PCR tests from CNS specimens and radiological evidence of acute CNS disease >48 hours after admission. Infection control practices were observed. A case-control analysis of medications and exposures was conducted; odds ratios and p-values were calculated via Fisher’s exact test.Genomic Analysis: Genomic DNA was extracted from patient-derived and environmental B. cereus isolates. Sequencing libraries were prepared and paired-end whole genome sequencing (WGS) was performed. Sequencing reads were computationally assembled using several existing tools (SPAdes, AllPaths-LG, VELVET, and Newbler). Draft genome annotation was performed using RAST and Prodigal. All single copy core genes were aligned using MUSCLE and a phylogenetic tree was generated using FastTree and RAxML.ResultsClinical Investigation: All 5 cases were AML patients who received chemotherapy. All had neutropenic fevers early in their course that resolved with empiric antibiotics. Fever recurred > 1 week into hospitalization (median onset 17 days, range 9-20) with headaches, neurological changes, and abnormal CNS imaging. Three had documented B. cereus bacteremia with bacilli identified on CNS pathological stains; two had positive CNS cultures. Four patients expired. Autopsies were notable for intestinal ulcers with Gram-positive bacilli in 2 patients.No infection control breaches were observed. Of 43 environmental samples, 4 tested positive for B.cereus, including a kitchen cart where bananas were stored, the shelf of a blanket warmer, and an air sample near a construction site. On case-control analysis, no medications were associated with cases. The only significant food exposures were bananas (OR 9.3, p=0.04) and cranberry juice (OR 9.6, p=0.04).Genomic Analysis:WGS data was obtained for all B. cereus isolates from BWH (n=4) and swabs of environmental screening efforts from BWH (n=4). Two patient isolates (7299 and 6254) were taken forward for pilot analysis. A comparison of various computational assembly tools demonstrated that the SPAdes assembler exhibited superior performance. SPAdes was thus used for draft genome assembly for the remaining isolates. Phylogenetic analysis of isolates 7299 and 6254 revealed that the organisms were genetically distinct from one another (Fig. 1). Genomic analysis and multilocus sequence typing (performed by WGS and independently by PCR amplification/Sanger sequencing) confirmed that all strains were genetically distinct from one another, and suggested that the environmental strains were relatively divergent from the pathogenic strains studied.ConclusionsNeuroinvasive B.cereus infection should be suspected in febrile AML patients with acute neurological symptoms, even those with negative blood cultures.The clinical investigation suggests the outbreak may have been attributable to gastrointestinal inoculation, possibly secondary to bananas. Of note, the B. cereus strain cultured from the suspected banana source was not identical to the pathogenic strains. Multiple changes in practice were instituted as a consequence of this outbreak including enhanced environmental cleanings, temporary modification of empiric antibiotic recommendations for recurrent neutropenic fever (vancomycin, cefepime, and ciprofloxacin) and the placement of all AML patients on diets that excluded fresh fruits, vegetables, and rice. Six months have elapsed since the last case.WGS demonstrated that the case series was not the result of a single, hyper-virulent B. cereus organism. Rather, the infections appear to have been attributable to distinct B. cereus strains that fall into two distinct genotypic clusters. Ongoing comparative microbial genomic and biological efforts are underway and have identified candidate bacterial and host factors that may induce heightened virulence in this patient population.1. Inoue et al. Leukemia & lymphoma. 2010;51(5):860-9. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Comparative Genomic Analysis of N2-Fixing and Non-N2-Fixing Paenibacillus spp.: Organization, Evolution and Expression of the Nitrogen Fixation Genes

We provide here a comparative genome analysis of 31 strains within the genus Paenibacillus including 11 new genomic sequences of N2-fixing strains. The heterogeneity of the 31 genomes (15 N2-fixing and 16 non-N2-fixing Paenibacillus strains) was reflected in the large size of the shell genome, which makes up approximately 65.2% of the genes in pan genome. Large numbers of transposable elements might be related to the heterogeneity. We discovered that a minimal and compact nif cluster comprising nine genes nifB, nifH, nifD, nifK, nifE, nifN, nifX, hesA and nifV encoding Mo-nitrogenase is conserved in the 15 N2-fixing strains. The nif cluster is under control of a σ70-depedent promoter and possesses a GlnR/TnrA-binding site in the promoter. Suf system encoding [Fe–S] cluster is highly conserved in N2-fixing and non-N2-fixing strains. Furthermore, we demonstrate that the nif cluster enabled Escherichia coli JM109 to fix nitrogen. Phylogeny of the concatenated NifHDK sequences indicates that Paenibacillus and Frankia are sister groups. Phylogeny of the concatenated 275 single-copy core genes suggests that the ancestral Paenibacillus did not fix nitrogen. The N2-fixing Paenibacillus strains were generated by acquiring the nif cluster via horizontal gene transfer (HGT) from a source related to Frankia. During the history of evolution, the nif cluster was lost, producing some non-N2-fixing strains, and vnf encoding V-nitrogenase or anf encoding Fe-nitrogenase was acquired, causing further diversification of some strains. In addition, some N2-fixing strains have additional nif and nif-like genes which may result from gene duplications. The evolution of nitrogen fixation in Paenibacillus involves a mix of gain, loss, HGT and duplication of nif/anf/vnf genes. This study not only reveals the organization and distribution of nitrogen fixation genes in Paenibacillus, but also provides insight into the complex evolutionary history of nitrogen fixation.