IslandViewer update: Improved genomic island discovery and visualization

Unveiling genome plasticity and a novel phage in Mycoplasma felis: Genomic investigations of four feline isolates

Mycoplasma felis has been isolated from diseased cats and horses, but to date only a single fully assembled genome of this species, of an isolate from a horse, has been characterized. This study aimed to characterize and compare the completely assembled genomes of four clinical isolates of M. felis from three domestic cats, assembled with the aid of short- and long-read sequencing methods. The completed genomes encoded a median of 759 ORFs (range 743-777) and had a median average nucleotide identity of 98.2 % with the genome of the available equid origin reference strain. Comparative genomic analysis revealed the occurrence of multiple horizontal gene transfer events and significant genome reassortment. This had resulted in the acquisition or loss of numerous genes within the Australian felid isolate genomes, encoding putative proteins involved in DNA transfer, metabolism, DNA replication, host cell interaction and restriction modification systems. Additionally, a novel mycoplasma phage was detected in one Australian felid M. felis isolate by genomic analysis and visualized using cryo-transmission electron microscopy. This study has highlighted the complex genomic dynamics in different host environments. Furthermore, the sequences obtained in this work will enable the development of new diagnostic tools, and identification of future infection control and treatment options for the respiratory disease complex in cats.

Delving into the lifestyle of Sundarban Wetland resident, biofilm producing, halotolerant Salinicoccus roseus: a comparative genomics-based intervention

BACKGROUND

Microbial community played an essential role in ecosystem processes, be it mangrove wetland or other intertidal ecologies. Several enzymatic activities like hydrolases are effective ecological indicators of soil microbial function. So far, little is known on halophilic bacterial contribution and function on a genomic viewpoint of Indian Sundarban Wetland. Considering the above mentioned issues, the aims of this study was to understand the life style, metabolic functionalities and genomic features of the isolated bacterium, Salinicoccus roseus strain RF1H. A comparative genome-based study of S. roseus has not been reported yet. Henceforth, we have considered the inclusion of the intra-species genome comparison of S. roseus to gain insight into the high degree of variation in the genome of strain RF1H among others.

RESULTS

Salinicoccus roseus strain RF1H is a pink-red pigmented, Gram-positive and non-motile cocci. The bacterium exhibited high salt tolerance (up to 15% NaCl), antibiotic resistance, biofilm formation and secretion of extracellular hydrolytic enzymes. The circular genome was approximately 2.62978 Mb in size, encoding 574 predicted genes with GC content 49.5%. Presence of genomic elements (prophages, transposable elements, CRISPR-Cas system) represented bacterial virulence and multidrug-resistance. Furthermore, genes associated with salt tolerance, temperature adaptation and DNA repair system were distributed in 17 genomic islands. Genes related to hydrocarbon degradation manifested metabolic capability of the bacterium for potential biotechnological applications. A comparative pangenome analysis revealed two-component response regulator, modified C4-dicarboxylate transport system and osmotic stress regulated ATP-binding proteins. Presence of genes encoding arginine decarboxylase (ADC) enzyme being involved in biofilm formation was reported from the genome. In silico study revealed the protein is thermostable and made up with ~ 415 amino acids, and hydrophilic in nature. Three motifs appeared to be evolutionary conserved in all Salinicoccus sequences.

CONCLUSION

The first report of whole genome analysis of Salinicoccus roseus strain RF1H provided information of metabolic functionalities, biofilm formation, resistance mechanism and adaptation strategies to thrive in climate-change induced vulnerable spot like Sundarban. Comparative genome analysis highlighted the unique genome content that contributed the strain's adaptability. The biomolecules produced during metabolism are important sources of compounds with potential beneficial applications in pharmaceuticals.

Comparative Genomic Analysis of Shrimp-Pathogenic Vibrio parahaemolyticus LC and Intraspecific Strains with Emphasis on Virulent Factors of Mobile Genetic Elements

Vibrio parahaemolyticus exhibits severe pathogenicity in humans and animals worldwide. In this study, genome sequencing and comparative analyses were conducted for in-depth characterization of the virulence factor (VF) repertoire of V. parahaemolyticus strain LC, which presented significant virulence to shrimp Litopenaeus vannamei. Strain LC, harboring two circular chromosomes and three linear plasmids, demonstrated ≥98.14% average nucleotide identities with 31 publicly available V. parahaemolyticus genomes, including 13, 11, and 7 shrimp-, human-, and non-pathogenic strains, respectively. Phylogeny analysis based on dispensable genes of pan-genome clustered 11 out of 14 shrimp-pathogenic strains and 7 out of 11 clinical strains into two distinct clades, indicating the close association between host-specific pathogenicity and accessory genes. The VFDB database revealed that 150 VFs of LC were mainly associated with the secretion system, adherence, antiphagocytosis, chemotaxis, motility, and iron uptake, whereas no homologs of the typical pathogenic genes pirA, pirB, tdh, and trh were detected. Four genes, mshB, wbfT, wbfU, and wbtI, were identified in both types of pathogenic strains but were absent in non-pathogens. Notably, a unique cluster similar to Yen-Tc, which encodes an insecticidal toxin complex, and diverse toxin-antitoxin (TA) systems, were identified on the mobile genetic elements (MGEs) of LC. Conclusively, in addition to the common VFs, various unique MGE-borne VFs, including the Yen-Tc cluster, TA components, and multiple chromosome-encoded chitinase genes, may contribute to the full spectrum of LC virulence. Moreover, V. parahaemolyticus demonstrates host-specific virulence, which potentially drives the origin and spread of pathogenic factors.

Systematic comparison of genome information processing and boundary recognition tools used for genomic island detection

Genomic islands are fragments of foreign DNA that are found in bacterial and archaeal genomes, and are typically associated with symbiosis or pathogenesis. While numerous genomic island detection methods have been proposed, there has been limited evaluation of the efficiency of the genome information processing and boundary recognition tools. In this study, we conducted a review of the statistical methods involved in genomic signatures, host signature extraction, informative signature selection, divergence measures, and boundary detection steps in genomic island prediction. We compared the performances of these methods on simulated experiments using alien fragments obtained from both artificial and real genomes. Our results indicate that among the nine genomic signatures evaluated, genomic signature frequency and full probability performed the best. However, their performance declined when normalized to their expectations and variances, such as Z-score and composition vector. Based on our experiments of the E. coli genome, we found that the confidence intervals of the window variances achieved the best performance in the signature extraction of the host, with the best confidence interval being 1.5-2 times the standard error. Ordered kurtosis was most effective in selecting informative signatures from a single genome, without requiring prior knowledge from other datasets. Among the three divergence measures evaluated, the two-sample t-test was the most successful, and a non-overlapping window with a small eye window (size 2) was best suited for identifying compositionally distinct regions. Finally, the maximum of the Markovian Jensen-Shannon divergence score, in terms of GC-content bias, was found to make boundary detection faster while maintaining a similar error rate.

Comparative genomics of food-derived probiotic Lactiplantibacillus plantarum K25 reveals its hidden potential, compactness, and efficiency

This study aimed to investigate the intricate genetic makeup of the Lactiplantibacillus plantarum K25 strain by conducting a comprehensive analysis of comparative genomics. The results of our study demonstrate that the genome exhibits a high-level efficiency and compactness, comprising a total of 3,199 genes that encode proteins and a GC content of 43.38%. The present study elucidates the evolutionary lineage of Lactiplantibacillus plantarum strains through an analysis of the degree of gene order conservation and synteny across a range of strains, thereby underscoring their closely interrelated evolutionary trajectories. The identification of various genetic components in the K25 strain, such as bacteriocin gene clusters and prophage regions, highlights its potential utility in diverse domains, such as biotechnology and medicine. The distinctive genetic elements possess the potential to unveil innovative therapeutic and biotechnological remedies in future. This study provides a comprehensive analysis of the L. plantarum K25 strain, revealing its remarkable genomic potential and presenting novel prospects for utilizing its unique genetic features in diverse scientific fields. The present study contributes to the existing literature on Lactiplantibacillus plantarum and sets the stage for prospective investigations and practical implementations that leverage the exceptional genetic characteristics of this adap organism.

Translational informatics for human microbiota: data resources, models and applications

With the rapid development of human intestinal microbiology and diverse microbiome-related studies and investigations, a large amount of data have been generated and accumulated. Meanwhile, different computational and bioinformatics models have been developed for pattern recognition and knowledge discovery using these data. Given the heterogeneity of these resources and models, we aimed to provide a landscape of the data resources, a comparison of the computational models and a summary of the translational informatics applied to microbiota data. We first review the existing databases, knowledge bases, knowledge graphs and standardizations of microbiome data. Then, the high-throughput sequencing techniques for the microbiome and the informatics tools for their analyses are compared. Finally, translational informatics for the microbiome, including biomarker discovery, personalized treatment and smart healthcare for complex diseases, are discussed.

A Maverick-like cluster in the genome of a pathogenic, moderately virulent strain of Gallibacterium anatis, ESV200, a transient biofilm producer

Introduction

Gallibacterium anatis causes gallibacteriosis in birds. These bacteria produce biofilms and secrete several fimbrial appendages as tools to cause disease in animals. G. anatis strains contain up to three types of fimbriae. Complete genome sequencing is the strategy currently used to determine variations in the gene content of G. anatis, although today only the completely circularized genome of G. anatis UMN179 is available.

Methods

The appearance of growth of various strains of G. anatis in liquid culture medium was studied. Biofilm production and how the amount of biofilm was affected by DNase, Proteinase K, and Pronase E enzymes were analyzed. Fimbrial gene expression was performed by protein analysis and qRT-PCR. In an avian model, the pathogenesis generated by the strains G. anatis ESV200 and 12656-12 was investigated. Using bioinformatic tools, the complete genome of G. anatis ESV200 was comparatively studied to search for virulence factors that would help explain the pathogenic behavior of this strain.

Results and Discussion

G. anatis ESV200 strain differs from the 12656-12 strain because it produces a biofilm at 20%. G. anatis ESV200 strain express fimbrial genes and produces biofilm but with a different structure than that observed for strain 12656-12. ESV200 and 12656-12 strains are pathogenic for chickens, although the latter is the most virulent. Here, we show that the complete genome of the ESV200 strain is similar to that of the UNM179 strain. However, these strains have evolved with many structural rearrangements; the most striking chromosomal arrangement is a Maverick-like element present in the ESV200 strain.

Comparative genome analysis of Streptococcus strains to identify virulent genes causing neonatal meningitis

AIM

To determine Streptococcus agalactiae genes responsible for causing neonatal meningitis.

BACKGROUND

Streptococcus agalactiae strain 2603 V/R is causative agent of neonatal meningitis, maternal infection and sepsis in young children. World health organisation reported high burden of new born death caused by this bacterium. Streptococcus agalactiae colonizing epithelial cells of vagina and endothelial cells have high resistance to available antibiotic drugs which makes it essential to determine new drug targets.

OBJECTIVES

To compare the genome of selected strain with the non-pathogenic strains of streptococcus and identify the virulent and antibiotic resistant genes for adaptation in host environment.

METHOD

The whole genome of human pathogen Streptococcus agalactiae strain 2603 V/R was analysed and compared with Streptococcus dysgalactiae strains using visualization and annotation tools. Genomic islands, mobile genetic elements, virulent and resistant genes were studied.

RESULTS

Genetically pathogenic strain is most similar to Streptococcus dysgalactiae subsp. equisimilis strain NCTC 7136. Comparative analysis revealed the importance of capsular polysaccharides and surface proteins responsible for avoiding immune system attachment to host epithelial cells and virulent behaviour. High number of genes coding for antibiotics resistance may provide a competitive advantage for survival of pathogenic Streptococcus agalactiae strain 2603 V/R in its niche.

CONCLUSIONS

The comparative analysis of pathogenic strain Streptococcus agalactiae with non-pathogenic strains of Streptococcus dysgalactiae provided new insights in pathogenicity that could aid in recognization for new regions and genes for development of new drug development strategies considering presence of high number of resistance genes.

Comparative Genomic Analysis of Seven Vibrio alginolyticus Strains Isolated From Shrimp Larviculture Water With Emphasis on Chitin Utilization

The opportunistic pathogen Vibrio alginolyticus is gaining attention because of its disease-causing risks to aquatic animals and humans. In this study, seven Vibrio strains isolated from different shrimp hatcheries in Southeast China were subjected to genome sequencing and subsequent comparative analysis to explore their intricate relationships with shrimp aquaculture. The seven isolates had an average nucleotide identity of ≥ 98.3% with other known V. alginolyticus strains. The species V. alginolyticus had an open pan-genome, with the addition of ≥ 161 novel genes following each new genome for seven isolates and 14 publicly available V. alginolyticus strains. The percentages of core genes of the seven strains were up to 83.1–87.5%, indicating highly conserved functions, such as chitin utilization. Further, a total of 14 core genes involved in the chitin degradation pathway were detected on the seven genomes with a single copy, 12 of which had undergone significant purifying selection (dN/dS < 1). Moreover, the seven strains could utilize chitin as the sole carbon-nitrogen source. In contrast, mobile genetic elements (MGEs) were identified in seven strains, including plasmids, prophages, and genomic islands, which mainly encoded accessory genes annotated as hypothetical proteins. The infection experiment showed that four of the seven strains might be pathogenic because the survival rates of Litopenaeus vannamei postlarvae were significantly reduced (P < 0.05) when compared to the control. However, no obvious correlation was noted between the number of putative virulence factors and toxic effects of the seven strains. Collectively, the persistence of V. alginolyticus in various aquatic environments may be attributed to its high genomic plasticity via the acquisition of novel genes by various MGEs. In view of the strong capability of chitin utilization by diverse vibrios, the timely removal of massive chitin-rich materials thoroughly in shrimp culture systems may be a key strategy to inhibit proliferation of vibrios and subsequent infection of shrimp. In addition, transcontinental transfer of potentially pathogenic V. alginolyticus strains should receive great attention to avoid vibriosis.

Enabling genomic island prediction and comparison in multiple genomes to investigate bacterial evolution and outbreaks

Outbreaks of virulent and/or drug-resistant bacteria have a significant impact on human health and major economic consequences. Genomic islands (GIs; defined as clusters of genes of probable horizontal origin) are of high interest because they disproportionately encode virulence factors, some antimicrobial-resistance (AMR) genes, and other adaptations of medical or environmental interest. While microbial genome sequencing has become rapid and inexpensive, current computational methods for GI analysis are not amenable for rapid, accurate, user-friendly and scalable comparative analysis of sets of related genomes. To help fill this gap, we have developed IslandCompare, an open-source computational pipeline for GI prediction and comparison across several to hundreds of bacterial genomes. A dynamic and interactive visualization strategy displays a bacterial core-genome phylogeny, with bacterial genomes linearly displayed at the phylogenetic tree leaves. Genomes are overlaid with GI predictions and AMR determinants from the Comprehensive Antibiotic Resistance Database (CARD), and regions of similarity between the genomes are also displayed. GI predictions are performed using Sigi-HMM and IslandPath-DIMOB, the two most precise GI prediction tools based on nucleotide composition biases, as well as a novel blast-based consistency step to improve cross-genome prediction consistency. GIs across genomes sharing sequence similarity are grouped into clusters, further aiding comparative analysis and visualization of acquisition and loss of mobile GIs in specific sub-clades. IslandCompare is an open-source software that is containerized for local use, plus available via a user-friendly, web-based interface to allow direct use by bioinformaticians, biologists and clinicians (at https://islandcompare.ca).

Enhanced production of polyhydroxyalkanoates in Pseudomonas putida KT2440 by a combination of genome streamlining and promoter engineering

Polyhydroxyalkanoates (PHAs), a class of bioplastics produced by a variety of microorganisms, have become the ideal alternatives for oil-derived plastics due to their superior physicochemical and material characteristics. Pseudomonas putida KT2440 can produce medium-chain-length PHA (mcl-PHA) from various substrates. In this study, a novel strategy of the large-scale deletion of genomic islands (GIs) coupling with promoter engineering was developed in P. putida KT2440 for constructing the minimal genome cell factories (MGF) capable of efficiently producing mcl-PHA. Firstly, P. putida KTU-U13, a 13 GIs- and upp-deleted mutant derived from the parental strain P. putida KT2440, was used as a starting strain for further deletion of GIs to generate a series of genome-reduced strains. Subsequently, the two minimal genome strains KTU-U24 and KTU-U27, which had a 7.19% and 8.35% reduction relative to the genome size of KT2440 and were advantageous over the strain KTU (KT2440∆upp) and KTU-U13 in several physiological traits such as the maximum specific growth rate, plasmid transformation efficiency, heterologous protein expression capacity and PHA production capacity, were selected as the chassis cells for PHA metabolic engineering. To prevent the formation of the by-product gluconic acid, the glucose dehydrogenase gene was deleted in KTU-U24 and KTU-U27, resulting in KTU-U24∆gcd and KTU-U27∆gcd. To enhance the transcriptional level of PHA synthase genes (phaC) and the supply of the precursor acetyl-CoA, a strong endogenous promoter P46 was inserted into upstream of the phaC operon and pyruvate dehydrogenase gene in the genome of KTU-U24∆gcd and KTU-U27∆gcd, to generate KTU-U24∆gcd-P46CA and KTU-U27∆gcd-P46CA, with the PHA yield of 50.5 wt% and 53.8 wt% (weight percent of PHA in cell dry weight). Finally, KTU-U27∆gcd-P46CA, the most minimal KT2440 chassis currently available, was able to accumulate the PHA to 55.82 wt% in a 5-l fermentor, which is the highest PHA yield obtained with P. putida KT2440 so far. This study suggests that genome streamlining in combination with promoter engineering may be a feasible strategy for the development of the MGF for the efficient production of high value products. Moreover, further streamlining of the P. putida KT2440 genome has great potential to create the optimal chassis for synthetic biology applications.

Identification of the Genes of the Plant Pathogen Pseudomonas syringae MB03 Required for the Nematicidal Activity Against Caenorhabditis elegans Through an Integrated Approach

Nematicidal potential of the common plant pathogen Pseudomonas syringae has been recently identified against Caenorhabditis elegans. The current study was designed to investigate the detailed genetic mechanism of the bacterial pathogenicity by applying comparative genomics, transcriptomics, mutant library screening, and protein expression. Results showed that P. syringae strain MB03 could kill C. elegans in the liquid assay by gut colonization. The genome of P. syringae MB03 was sequenced and comparative analysis including multi locus sequence typing, and genome-to-genome distance placed MB03 in phylogroup II of P. syringae. Furthermore, comparative genomics of MB03 with nematicidal strains of Pseudomonas aeruginosa (PAO1 and PA14) predicted 115 potential virulence factors in MB03. However, genes for previously reported nematicidal metabolites, such as phenazine, pyochelin, and pyrrolnitrin, were found absent in the MB03 genome. Transcriptomics analysis showed that the growth phase of the pathogen considerably affected the expression of virulence factors, as genes for the flagellum, glutamate ABC transporter, phoP/phoQ, fleS/fleR, type VI secretion system, and serralysin were highly up-regulated when stationary phase MB03 cells interacted with C. elegans. Additionally, screening of a transposon insertion mutant library led to the identification of other nematicidal genes such as acnA, gltP, oprD, and zapE. Finally, the nematicidal activity of selected proteins was confirmed by heterologous expression in Escherichia coli.

Genomic Island Prediction via Chi-Square Test and Random Forest Algorithm

Genomic islands are related to microbial adaptation and carry different genomic characteristics from the host. Therefore, many methods have been proposed to detect genomic islands from the rest of the genome by evaluating its sequence composition. Many sequence features have been proposed, but many of them have not been applied to the identification of genomic islands. In this paper, we present a scheme to predict genomic islands using the chi-square test and random forest algorithm. We extract seven kinds of sequence features and select the important features with the chi-square test. All the selected features are then input into the random forest to predict the genome islands. Three experiments and comparison show that the proposed method achieves the best performance. This understanding can be useful to design more powerful method for the genomic island prediction.

Microbial genomic island discovery, visualization and analysis

Horizontal gene transfer (also called lateral gene transfer) is a major mechanism for microbial genome evolution, enabling rapid adaptation and survival in specific niches. Genomic islands (GIs), commonly defined as clusters of bacterial or archaeal genes of probable horizontal origin, are of particular medical, environmental and/or industrial interest, as they disproportionately encode virulence factors and some antimicrobial resistance genes and may harbor entire metabolic pathways that confer a specific adaptation (solvent resistance, symbiosis properties, etc). As large-scale analyses of microbial genomes increases, such as for genomic epidemiology investigations of infectious disease outbreaks in public health, there is increased appreciation of the need to accurately predict and track GIs. Over the past decade, numerous computational tools have been developed to tackle the challenges inherent in accurate GI prediction. We review here the main types of GI prediction methods and discuss their advantages and limitations for a routine analysis of microbial genomes in this era of rapid whole-genome sequencing. An assessment is provided of 20 GI prediction software methods that use sequence-composition bias to identify the GIs, using a reference GI data set from 104 genomes obtained using an independent comparative genomics approach. Finally, we present guidelines to assist researchers in effectively identifying these key genomic regions.

Genomic variation among closely related Vibrio alginolyticus strains is located on mobile genetic elements

BACKGROUND

Species of the genus Vibrio, one of the most diverse bacteria genera, have undergone niche adaptation followed by clonal expansion. Niche adaptation and ultimately the formation of ecotypes and speciation in this genus has been suggested to be mainly driven by horizontal gene transfer (HGT) through mobile genetic elements (MGEs). Our knowledge about the diversity and distribution of Vibrio MGEs is heavily biased towards human pathogens and our understanding of the distribution of core genomic signatures and accessory genes encoded on MGEs within specific Vibrio clades is still incomplete. We used nine different strains of the marine bacterium Vibrio alginolyticus isolated from pipefish in the Kiel-Fjord to perform a multiscale-comparative genomic approach that allowed us to investigate [1] those genomic signatures that characterize a habitat-specific ecotype and [2] the source of genomic variation within this ecotype.

RESULTS

We found that the nine isolates from the Kiel-Fjord have a closed-pangenome and did not differ based on core-genomic signatures. Unique genomic regions and a unique repertoire of MGEs within the Kiel-Fjord isolates suggest that the acquisition of gene-blocks by HGT played an important role in the evolution of this ecotype. Additionally, we found that ~ 90% of the genomic variation among the nine isolates is encoded on MGEs, which supports ongoing theory that accessory genes are predominately located on MGEs and shared by HGT. Lastly, we could show that these nine isolates share a unique virulence and resistance profile which clearly separates them from all other investigated V. alginolyticus strains and suggests that these are habitat-specific genes, required for a successful colonization of the pipefish, the niche of this ecotype.

CONCLUSION

We conclude that all nine V. alginolyticus strains from the Kiel-Fjord belong to a unique ecotype, which we named the Kiel-alginolyticus ecotype. The low sequence variation of the core-genome in combination with the presence of MGE encoded relevant traits, as well as the presence of a suitable niche (here the pipefish), suggest, that this ecotype might have evolved from a clonal expansion following HGT driven niche-adaptation.

2SigFinder: the combined use of small-scale and large-scale statistical testing for genomic island detection from a single genome

BACKGROUND

Genomic islands are associated with microbial adaptations, carrying genomic signatures different from the host. Some methods perform an overall test to identify genomic islands based on their local features. However, regions of different scales will display different genomic features.

RESULTS

We proposed here a novel method "2SigFinder ", the first combined use of small-scale and large-scale statistical testing for genomic island detection. The proposed method was tested by genomic island boundary detection and identification of genomic islands or functional features of real biological data. We also compared the proposed method with the comparative genomics and composition-based approaches. The results indicate that the proposed 2SigFinder is more efficient in identifying genomic islands.

CONCLUSIONS

From real biological data, 2SigFinder identified genomic islands from a single genome and reported robust results across different experiments, without annotated information of genomes or prior knowledge from other datasets. 2SigHunter identified 25 Pathogenicity, 1 tRNA, 2 Virulence and 2 Repeats from 27 Pathogenicity, 1 tRNA, 2 Virulence and 2 Repeats, and detected 101 Phage and 28 HEG out of 130 Phage and 36 HEGs in S. enterica Typhi CT18, which shows that it is more efficient in detecting functional features associated with GIs.

MTGIpick allows robust identification of genomic islands from a single genome

Genomic islands (GIs) that are associated with microbial adaptations and carry sequence patterns different from that of the host are sporadically distributed among closely related species. This bias can dominate the signal of interest in GI detection. However, variations still exist among the segments of the host, although no uniform standard exists regarding the best methods of discriminating GIs from the rest of the genome in terms of compositional bias. In the present work, we proposed a robust software, MTGIpick, which used regions with pattern bias showing multiscale difference levels to identify GIs from the host. MTGIpick can identify GIs from a single genome without annotated information of genomes or prior knowledge from other data sets. When real biological data were used, MTGIpick demonstrated better performance than existing methods, as well as revealed potential GIs with accurate sizes missed by existing methods because of a uniform standard. Software and supplementary are freely available at http://bioinfo.zstu.edu.cn/MTGI or https://github.com/bioinfo0706/MTGIpick.

An efficient genomic signature ranking method for genomic island prediction from a single genome

Genomic islands that are associated with microbial adaptations and carry genomic signatures different from that of the host, and thus many methods have been proposed to select the informative genomic signatures from a range of organisms and discriminate genomic islands from the rest of the genome in terms of these signature biases. However, they are of limited use when closely related genomes are unavailable. In the present work, we proposed a kurtosis-based ranking method to select the informative genomic signatures from a single genome. In simulations with alien fragments from artificial and real genomes, the proposed kurtosis-based ranking method efficiently selected the informative genomic signatures from a single genome, without annotated information of genomes or prior knowledge from other datasets. This understanding can be useful to design more powerful method for genomic island detection.

Shared genome analyses of notable listeriosis outbreaks, highlighting the critical importance of epidemiological evidence, input datasets and interpretation criteria

The persuasiveness of genomic evidence has pressured scientific agencies to supplement or replace well-established methodologies to inform public health and food safety decision-making. This study of 52 epidemiologically defined Listeria monocytogenes isolates, collected between 1981 and 2011, including nine outbreaks, was undertaken (1) to characterize their phylogenetic relationship at finished genome-level resolution, (2) to elucidate the underlying genetic diversity within an endemic subtype, CC8, and (3) to re-evaluate the genetic relationship and epidemiology of a CC8-delimited outbreak in Canada in 2008. Genomes representing Canadian Listeria outbreaks between 1981 and 2010 were closed and manually annotated. Single nucleotide variants (SNVs) and horizontally acquired traits were used to generate phylogenomic models. Phylogenomic relationships were congruent with classical subtyping and epidemiology, except for CC8 outbreaks, wherein the distribution of SNV and prophages revealed multiple co-evolving lineages. Chronophyletic reconstruction of CC8 evolution indicates that prophage-related genetic changes among CC8 strains manifest as PFGE subtype reversions, obscuring the relationship between CC8 isolates, and complicating the public health interpretation of subtyping data, even at maximum genome resolution. The size of the shared genome interrogated did not change the genetic relationship measured between highly related isolates near the tips of the phylogenetic tree, illustrating the robustness of these approaches for routine public health applications where the focus is recent ancestry. The possibility exists for temporally and epidemiologically distinct events to appear related even at maximum genome resolution, highlighting the continued importance of epidemiological evidence.

Genomic and Geographic Context for the Evolution of High-Risk Carbapenem-Resistant Enterobacter cloacae Complex Clones ST171 and ST78

Recent reports have established the escalating threat of carbapenem-resistant Enterobacter cloacae complex (CREC). Here, we demonstrate that CREC has evolved as a highly antibiotic-resistant rather than highly virulent nosocomial pathogen. Applying genomics and Bayesian phylogenetic analyses to a 7-year collection of CREC isolates from a northern Manhattan hospital system and to a large set of publicly available, geographically diverse genomes, we demonstrate clonal spread of a single clone, ST171. We estimate that two major clades of epidemic ST171 diverged prior to 1962, subsequently spreading in parallel from the Northeastern to the Mid-Atlantic and Midwestern United States and demonstrating links to international sites. Acquisition of carbapenem and fluoroquinolone resistance determinants by both clades preceded widespread use of these drugs in the mid-1980s, suggesting that antibiotic pressure contributed substantially to its spread. Despite a unique mobile repertoire, ST171 isolates showed decreased virulence in vitro While a second clone, ST78, substantially contributed to the emergence of CREC, it encompasses diverse carbapenemase-harboring plasmids, including a potentially hypertransmissible IncN plasmid, also present in other sequence types. Rather than heightened virulence, CREC demonstrates lineage-specific, multifactorial adaptations to nosocomial environments coupled with a unique potential to acquire and disseminate carbapenem resistance genes. These findings indicate a need for robust surveillance efforts that are attentive to the potential for local and international spread of high-risk CREC clones.IMPORTANCE Carbapenem-resistant Enterobacter cloacae complex (CREC) has emerged as a formidable nosocomial pathogen. While sporadic acquisition of plasmid-encoded carbapenemases has been implicated as a major driver of CREC, ST171 and ST78 clones demonstrate epidemic potential. However, a lack of reliable genomic references and rigorous statistical analyses has left many gaps in knowledge regarding the phylogenetic context and evolutionary pathways of successful CREC. Our reconstruction of recent ST171 and ST78 evolution represents a significant addition to current understanding of CREC and the directionality of its spread from the Eastern United States to the northern Midwestern United States with links to international collections. Our results indicate that the remarkable ability of E. cloacae to acquire and disseminate cross-class antibiotic resistance rather than virulence determinants, coupled with its ability to adapt under conditions of antibiotic pressure, likely led to the wide dissemination of CREC.

Trajectories and Drivers of Genome Evolution in Surface-Associated Marine Phaeobacter

The extent of genome divergence and the evolutionary events leading to speciation of marine bacteria have mostly been studied for (locally) abundant, free-living groups. The genus Phaeobacter is found on different marine surfaces, seems to occupy geographically disjunct habitats, and is involved in different biotic interactions, and was therefore targeted in the present study. The analysis of the chromosomes of 32 closely related but geographically spread Phaeobacter strains revealed an exceptionally large, highly syntenic core genome. The flexible gene pool is constantly but slightly expanding across all Phaeobacter lineages. The horizontally transferred genes mostly originated from bacteria of the Roseobacter group and horizontal transfer most likely was mediated by gene transfer agents. No evidence for geographic isolation and habitat specificity of the different phylogenomic Phaeobacter clades was detected based on the sources of isolation. In contrast, the functional gene repertoire and physiological traits of different phylogenomic Phaeobacter clades were sufficiently distinct to suggest an adaptation to an associated lifestyle with algae, to additional nutrient sources, or toxic heavy metals. Our study reveals that the evolutionary trajectories of surface-associated marine bacteria can differ significantly from free-living marine bacteria or marine generalists.

Exploring the Genomic Traits of Non-toxigenic Vibrio parahaemolyticus Strains Isolated in Southern Chile

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis worldwide. As reported in other countries, after the rise and fall of the pandemic strain in Chile, other post-pandemic strains have been associated with clinical cases, including strains lacking the major toxins TDH and TRH. Since the presence or absence of tdh and trh genes has been used for diagnostic purposes and as a proxy of the virulence of V. parahaemolyticus isolates, the understanding of virulence in V. parahaemolyticus strains lacking toxins is essential to detect these strains present in water and marine products to avoid possible food-borne infection. In this study, we characterized the genome of four environmental and two clinical non-toxigenic strains (tdh-, trh-, and T3SS2-). Using whole-genome sequencing, phylogenetic, and comparative genome analysis, we identified the core and pan-genome of V. parahaemolyticus of strains of southern Chile. The phylogenetic tree based on the core genome showed low genetic diversity but the analysis of the pan-genome revealed that all strains harbored genomic islands carrying diverse virulence and fitness factors or prophage-like elements that encode toxins like Zot and RTX. Interestingly, the three strains carrying Zot-like toxin have a different sequence, although the alignment showed some conserved areas with the zot sequence found in V. cholerae. In addition, we identified an unexpected diversity in the genetic architecture of the T3SS1 gene cluster and the presence of the T3SS2 gene cluster in a non-pandemic environmental strain. Our study sheds light on the diversity of V. parahaemolyticus strains from the southern Pacific which increases our current knowledge regarding the global diversity of this organism.

xenoGI: reconstructing the history of genomic island insertions in clades of closely related bacteria

Background

Genomic islands play an important role in microbial genome evolution, providing a mechanism for strains to adapt to new ecological conditions. A variety of computational methods, both genome-composition based and comparative, have been developed to identify them. Some of these methods are explicitly designed to work in single strains, while others make use of multiple strains. In general, existing methods do not identify islands in the context of the phylogeny in which they evolved. Even multiple strain approaches are best suited to identifying genomic islands that are present in one strain but absent in others. They do not automatically recognize islands which are shared between some strains in the clade or determine the branch on which these islands inserted within the phylogenetic tree.

Results

We have developed a software package, xenoGI, that identifies genomic islands and maps their origin within a clade of closely related bacteria, determining which branch they inserted on. It takes as input a set of sequenced genomes and a tree specifying their phylogenetic relationships. Making heavy use of synteny information, the package builds gene families in a species-tree-aware way, and then attempts to combine into islands those families whose members are adjacent and whose most recent common ancestor is shared. The package provides a variety of text-based analysis functions, as well as the ability to export genomic islands into formats suitable for viewing in a genome browser. We demonstrate the capabilities of the package with several examples from enteric bacteria, including an examination of the evolution of the acid fitness island in the genus Escherichia. In addition we use output from simulations and a set of known genomic islands from the literature to show that xenoGI can accurately identify genomic islands and place them on a phylogenetic tree.

Conclusions

xenoGI is an effective tool for studying the history of genomic island insertions in a clade of microbes. It identifies genomic islands, and determines which branch they inserted on within the phylogenetic tree for the clade. Such information is valuable because it helps us understand the adaptive path that has produced living species.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2038-0) contains supplementary material, which is available to authorized users.

Transcriptional changes when Myxococcus xanthus preys on Escherichia coli suggest myxobacterial predators are constitutively toxic but regulate their feeding

Predation is a fundamental ecological process, but within most microbial ecosystems the molecular mechanisms of predation remain poorly understood. We investigated transcriptome changes associated with the predation of Escherichia coli by the myxobacterium Myxococcus xanthus using mRNA sequencing. Exposure to pre-killed prey significantly altered expression of 1319 predator genes. However, the transcriptional response to living prey was minimal, with only 12 genes being significantly up-regulated. The genes most induced by prey presence (kdpA and kdpB, members of the kdp regulon) were confirmed by reverse transcriptase quantitative PCR to be regulated by osmotic shock in M. xanthus, suggesting indirect sensing of prey. However, the prey showed extensive transcriptome changes when co-cultured with predator, with 40 % of its genes (1534) showing significant changes in expression. Bacteriolytic M. xanthus culture supernatant and secreted outer membrane vesicles (OMVs) also induced changes in expression of large numbers of prey genes (598 and 461, respectively). Five metabolic pathways were significantly enriched in prey genes up-regulated on exposure to OMVs, supernatant and/or predatory cells, including those for ribosome and lipopolysaccharide production, suggesting that the prey cell wall and protein production are primary targets of the predator's attack. Our data suggest a model of the myxobacterial predatome (genes and proteins associated with predation) in which the predator constitutively produces secretions which disable its prey whilst simultaneously generating a signal that prey is present. That signal then triggers a regulated feeding response in the predator.

Zisland Explorer: detect genomic islands by combining homogeneity and heterogeneity properties

Genomic islands are genomic fragments of alien origin in bacterial and archaeal genomes, usually involved in symbiosis or pathogenesis. In this work, we described Zisland Explorer, a novel tool to predict genomic islands based on the segmental cumulative GC profile. Zisland Explorer was designed with a novel strategy, as well as a combination of the homogeneity and heterogeneity of genomic sequences. While the sequence homogeneity reflects the composition consistence within each island, the heterogeneity measures the composition bias between an island and the core genome. The performance of Zisland Explorer was evaluated on the data sets of 11 different organisms. Our results suggested that the true-positive rate (TPR) of Zisland Explorer was at least 10.3% higher than that of four other widely used tools. On the other hand, the new tool did not lose overall accuracy with the improvement in the TPR and showed better equilibrium among various evaluation indexes. Also, Zisland Explorer showed better accuracy in the prediction of experimental island data. Overall, the tool provides an alternative solution over other tools, which expands the field of island prediction and offers a supplement to increase the performance of the distinct predicting strategy. We have provided a web service as well as a graphical user interface and open-source code across multiple platforms for Zisland Explorer, which is available at http://cefg.uestc.edu.cn/Zisland_Explorer/ or http://tubic.tju.edu.cn/Zisland_Explorer/.

Adaptation of Surface-Associated Bacteria to the Open Ocean: A Genomically Distinct Subpopulation of Phaeobacter gallaeciensis Colonizes Pacific Mesozooplankton

The marine Roseobacter group encompasses numerous species which occupy a large variety of ecological niches. However, members of the genus Phaeobacter are specifically adapted to a surface-associated lifestyle and have so far been found nearly exclusively in disjunct, man-made environments including shellfish and fish aquacultures, as well as harbors. Therefore, the possible natural habitats, dispersal and evolution of Phaeobacter spp. have largely remained obscure. Applying a high-throughput cultivation strategy along a longitudinal Pacific transect, the present study revealed for the first time a widespread natural occurrence of Phaeobacter in the marine pelagial. These bacteria were found to be specifically associated to mesoplankton where they constitute a small but detectable proportion of the bacterial community. The 16S rRNA gene sequences of 18 isolated strains were identical to that of Phaeobacter gallaeciensis DSM26640^T but sequences of internal transcribed spacer and selected genomes revealed that the strains form a distinct clade within P. gallaeciensis. The genomes of the Pacific and the aquaculture strains were highly conserved and had a fraction of the core genome of 89.6%, 80 synteny breakpoints, and differed 2.2% in their nucleotide sequences. Diversification likely occurred through neutral mutations. However, the Pacific strains exclusively contained two active Type I restriction modification systems which is commensurate with a reduced acquisition of mobile elements in the Pacific clade. The Pacific clade of P. gallaeciensis also acquired a second, homolog phosphonate transport system compared to all other P. gallaeciensis. Our data indicate that a previously unknown, distinct clade of P. gallaeciensis acquired a limited number of clade-specific genes that were relevant for its association with mesozooplankton and for colonization of the marine pelagial. The divergence of the Pacific clade most likely was driven by the adaptation to this novel ecological niche rather than by geographic isolation.

A Systematic In-silico Analysis of Helicobacter pylori Pathogenic Islands for Identification of Novel Drug Target Candidates

Background:

Helicobacter pylori is associated with inflammation of different areas, such as the duodenum and stomach, causing gastritis and gastric ulcers leading to lymphoma and cancer. Pathogenic islands are a type of clustered mobile elements ranging from 10-200 Kb contributing to the virulence of the respective pathogen coding for one or more virulence factors. Virulence factors are molecules expressed and secreted by pathogen and are responsible for causing disease in the host. Bacterial genes/virulence factors of the pathogenic islands represent a promising source for identifying novel drug targets.

Objective:

The study aimed at identifying novel drug targets from pathogenic islands in H. pylori.

Material & Methods:

The genome of 23 H. pylori strains were screened for pathogenic islands and bacterial genes/virulence factors to identify drug targets. Protein-protein interactions of drug targets were predicted for identifying interacting partners. Further, host-pathogen interactions of interacting partners were predicted to identify important molecules which are closely associated with gastric cancer.

Results:

Screening the genome of 23 H. pylori strains revealed 642 bacterial genes/virulence factors in 31 pathogenic islands. Further analysis identified 101 genes which were non-homologous to human and essential for the survival of the pathogen, among them 31 are potential drug targets. Protein-protein interactions for 31 drug targets predicted 609 interacting partners. Predicted interacting partners were further subjected to host-pathogen interactions leading to identification of important molecules like TNF receptor associated factor 6, (TRAF6) and MAPKKK7 which are closely associated with gastric cancer.

Conclusion:

These provocative studies enabled us to identify important molecules in H. pylori and their counter interacting molecules in the host leading to gastric cancer and also a pool of novel drug targets for therapeutic intervention of gastric cancer.

Draft Genome Sequence of Klebsiella michiganensis 3T412C, Harboring an Arsenic Resistance Genomic Island, Isolated from Mine Tailings in Peru

An arsenic resistance genomic island in the bacterium Klebsiella michiganensis 3T412C was isolated from mine tailings from Peru. This genomic island confers adaptation to extreme environments with high concentrations of arsenic. Isolate 3T412C contained a complete set of genes involved in resistance to arsenic. This operon is surrounded by putative genes for resistance to other heavy metals.

Isolation and Whole-genome Sequence Analysis of the Imipenem Heteroresistant Acinetobacter baumannii Clinical Isolate HRAB-85

OBJECTIVES

Heteroresistance is a phenomenon in which there are various responses to antibiotics from bacterial cells within the same population. Here, we isolated and characterised an imipenem heteroresistant Acinetobacter baumannii strain (HRAB-85).

METHODS

The genome of strain HRAB-85 was completely sequenced and analysed to understand its antibiotic resistance mechanisms. Population analysis and multilocus sequence typing were performed.

RESULTS

Subpopulations grew in the presence of imipenem at concentrations of up to 64μg/mL, and the strain was found to belong to ST208. The total length of strain HRAB-85 was 4,098,585bp with a GC content of 39.98%. The genome harboured at least four insertion sequences: the common ISAba1, ISAba22, ISAba24, and newly reported ISAba26. Additionally, 19 antibiotic-resistance genes against eight classes of antimicrobial agents were found, and 11 genomic islands (GIs) were identified. Among them, GI3, GI10, and GI11 contained many ISs and antibiotic-resistance determinants.

CONCLUSIONS

The existence of imipenem heteroresistant phenotypes in A. baumannii was substantiated in this hospital, and imipenem pressure, which could induce imipenem-heteroresistant subpopulations, may select for highly resistant strains. The complete genome sequencing and bioinformatics analysis of HRAB-85 could improve our understanding of the epidemiology and resistance mechanisms of carbapenem-heteroresistant A. baumannii.

Characterization of the emerging zoonotic pathogen Arcobacter thereius by whole genome sequencing and comparative genomics

Four Arcobacter species have been associated with human disease, and based on current knowledge, these Gram negative bacteria are considered as potential food and waterborne zoonotic pathogens. At present, only the genome of the species Arcobacter butzleri has been analysed, and still little is known about their physiology and genetics. The species Arcobacter thereius has first been isolated from tissue of aborted piglets, duck and pig faeces, and recently from stool of human patients with enteritis. In the present study, the complete genome and analysis of the A. thereius type strain LMG24486T, as well as the comparative genome analysis with 8 other A. thereius strains are presented. Genome analysis revealed metabolic pathways for the utilization of amino acids, which represent the main source of energy, together with the presence of genes encoding for respiration-associated and chemotaxis proteins. Comparative genome analysis with the A. butzleri type strain RM4018 revealed a large correlation, though also unique features. Furthermore, in silico DDH and ANI based analysis of the nine A. thereius strains disclosed clustering into two closely related genotypes. No discriminatory differences in genome content nor phenotypic behaviour were detected, though recently the species Arcobacter porcinus was proposed to encompass part of the formerly identified Arcobacter thereius strains. The report of the presence of virulence associated genes in A. thereius, the presence of antibiotic resistance genes, verified by in vitro susceptibility testing, as well as other pathogenic related relevant features, support the classification of A. thereius as an emerging pathogen.

Genomic characterization of a local epidemic Pseudomonas aeruginosa reveals specific features of the widespread clone ST395

Pseudomonas aeruginosa is a ubiquitous opportunistic pathogen with several clones being frequently associated with outbreaks in hospital settings. ST395 is among these so-called ‘international’ clones. We aimed here to define the biological features that could have helped the implantation and spread of the clone ST395 in hospital settings. The complete genome of a multidrug resistant index isolate (DHS01) of a large hospital outbreak was analysed. We identified DHS01-specific genetic elements, among which were identified those shared with a panel of six independent ST395 isolates responsible for outbreaks in other hospitals. DHS01 has the fifth largest chromosome of the species (7.1 Mbp), with most of its 1555 accessory genes borne by either genomic islands (GIs, n=48) or integrative and conjugative elements (ICEs, n=5). DHS01 is multidrug resistant mostly due to chromosomal mutations. It displayed signatures of adaptation to chronic infection in part due to the loss of a 131 kbp chromosomal fragment. Four GIs were specific to the clone ST395 and contained genes involved in metabolism (GI-4), in virulence (GI-6) and in resistance to copper (GI-7). GI-7 harboured an array of six copper transporters and was shared with non-pathogenic Pseudomonas sp. retrieved from copper-contaminated environments. Copper resistance was confirmed phenotypically in all other ST395 isolates and possibly accounted for the spreading capability of the clone in hospital outbreaks, where water networks have been incriminated. This suggests that genes transferred from copper-polluted environments may have favoured the implantation and spread of the international clone P. aeruginosa ST395 in hospital settings.

A genomic view of food-related and probiotic Enterococcus strains

The study of enterococcal genomes has grown considerably in recent years. While special attention is paid to comparative genomic analysis among clinical relevant isolates, in this study we performed an exhaustive comparative analysis of enterococcal genomes of food origin and/or with potential to be used as probiotics. Beyond common genetic features, we especially aimed to identify those that are specific to enterococcal strains isolated from a certain food-related source as well as features present in a species-specific manner. Thus, the genome sequences of 25 Enterococcus strains, from 7 different species, were examined and compared. Their phylogenetic relationship was reconstructed based on orthologous proteins and whole genomes. Likewise, markers associated with a successful colonization (bacteriocin genes and genomic islands) and genome plasticity (phages and clustered regularly interspaced short palindromic repeats) were investigated for lifestyle specific genetic features. At the same time, a search for antibiotic resistance genes was carried out, since they are of big concern in the food industry. Finally, it was possible to locate 1617 FIGfam families as a core proteome universally present among the genera and to determine that most of the accessory genes code for hypothetical proteins, providing reasonable hints to support their functional characterization.

Comparative Genomics and Transcriptional Analysis of Flavobacterium columnare Strain ATCC 49512

Flavobacterium columnare is a Gram-negative fish pathogen causing columnaris disease in wild and cultured fish species. Although the pathogen is widespread in aquatic environments and fish worldwide, little is known about biology of F. columnare and mechanisms of columnaris disease pathogenesis. Previously we presented the complete genome sequence of F. columnare strain ATCC 49512. Here we present a comparison of the strain ATCC 49512 genome to four other Flavobacterium genomes. In this analysis, we identified predicted proteins whose functions indicate F. columnare is capable of denitrification, which would enable anaerobic growth in aquatic pond sediments. Anaerobic growth of F. columnare ATCC 49512 with nitrate supplementation was detected experimentally. F. columnare ATCC 49512 had a relatively high number of insertion sequences and genomic islands compared to the other Flavobacterium species, suggesting a larger degree of horizontal gene exchange and genome plasticity. A type VI subtype III secretion system was encoded in F. columnare along with F. johnsoniae and F. branchiophilum. RNA sequencing proved to be a valuable technique to improve annotation quality; 41 novel protein coding regions were identified, 16 of which had a non-traditional start site (TTG, GTG, and CTT). Candidate small noncoding RNAs were also identified. Our results improve our understanding of F. columnare ATCC 49512 biology, and our results support the use of RNA sequencing to improve annotation of bacterial genomes, particularly for type strains.

Comparative pan genome analysis of oral Prevotella species implicated in periodontitis

Prevotella is part of the oral bacterial community implicated in periodontitis. Pan genome analyses of eight oral Prevotella species, P. dentalis, P. enoeca, P. fusca, P. melaninogenica, P. denticola, P. intermedia 17, P. intermedia 17-2 and P. sp. oral taxon 299 are presented in this study. Analysis of the Prevotella pan genome revealed features such as secretion systems, resistance to oxidative stress and clustered regularly interspaced short palindromic repeat (CRISPR)-Cas systems that enable the bacteria to adapt to the oral environment. We identified the presence of type VI secretion system (T6SS) in P. fusca and P. intermedia strains. For some VgrG and Hcp proteins which were not part of the core T6SS loci, we used gene neighborhood analysis and identified putative effector proteins and putative polyimmunity loci in P. fusca and polymorphic toxin systems in P. intermedia strains. Earlier studies have identified the presence of Por secretion system (PorSS) in P. gingivalis, P. melaninogenica and P. intermedia. We noted the presence of their homologs in six other oral Prevotella studied here. We suggest that in Prevotella, PorSS is used to secrete cysteine proteases such as interpain and C-terminal domain containing proteins with a "Por_secre_tail" domain. We identified subtype I-B CRISPR-Cas system in P. enoeca. Putative CRISPR-Cas system subtypes for 37 oral Prevotella and 30 non-oral Prevotella species were also predicted. Further, we performed a BLASTp search of the Prevotella proteins which are also conserved in the red-complex pathogens, against the human proteome to identify potential broad-spectrum drug targets. In summary, the use of a pan genome approach enabled identification of secretion systems and defense mechanisms in Prevotella that confer adaptation to the oral cavity.

A comprehensive phylogenetic analysis of copper transporting P1B ATPases from bacteria of the Rhizobiales order uncovers multiplicity, diversity and novel taxonomic subtypes

The ubiquitous cytoplasmic membrane copper transporting P_1B‐1 and P_1B‐3‐type ATPases pump out Cu⁺ and Cu²⁺, respectively, to prevent cytoplasmic accumulation and avoid toxicity. The presence of five copies of Cu‐ATPases in the symbiotic nitrogen‐fixing bacteria Sinorhizobium meliloti is remarkable; it is the largest number of Cu⁺‐transporters in a bacterial genome reported to date. Since the prevalence of multiple Cu‐ATPases in members of the Rhizobiales order is unknown, we performed an in silico analysis to understand the occurrence, diversity and evolution of Cu⁺‐ATPases in members of the Rhizobiales order. Multiple copies of Cu‐ATPase coding genes (2–8) were detected in 45 of the 53 analyzed genomes. The diversity inferred from a maximum‐likelihood (ML) phylogenetic analysis classified Cu‐ATPases into four monophyletic groups. Each group contained additional subtypes, based on the presence of conserved motifs. This novel phylogeny redefines the current classification, where they are divided into two subtypes (P_1B‐1 and P_1B‐3). Horizontal gene transfer (HGT) as well as the evolutionary dynamic of plasmid‐borne genes may have played an important role in the functional diversification of Cu‐ATPases. Homologous cytoplasmic and periplasmic Cu⁺‐chaperones, CopZ, and CusF, that integrate a CopZ‐CopA‐CusF tripartite efflux system in gamma‐proteobacteria and archeae, were found in 19 of the 53 surveyed genomes of the Rhizobiales. This result strongly suggests a high divergence of CopZ and CusF homologs, or the existence of unexplored proteins involved in cellular copper transport.

Comparative Genome Analysis of the Daptomycin-Resistant Streptococcus anginosus Strain J4206 Associated with Breakthrough Bacteremia

Streptococcus anginosus is a member of the normal oral flora that can become a pathogen causing pyogenic infections in humans. The genome of daptomycin-resistant strain J4206, originally isolated from a patient suffering from breakthrough bacteremia and septic shock at the University of Texas Health Science Center at San Antonio, was determined. The circular genome is 2,001,352 bp long with a GC content of 38.62% and contains multiple mobile genetic elements, including the phage-like chromosomal island SanCI that mediates a mutator phenotype, transposons, and integrative conjugative elements. Daptomycin resistance involves multiple alterations in the cell membrane and cell wall, and unique features were identified in J4206 that may contribute to resistance. A cluster of capsular polysaccharide (CPS) genes for choline metabolism and transport are present that may help neutralize cell surface charges, destabilizing daptomycin binding. Further, unique J4206 genes encoding sortases and LPXTG-target proteins that are involved in cell wall modification were present. The J4206 genome is phylogenetically closely related to the recently reported vancomycin-resistant SA1 strain; however, these genomes differ with SNPs in cardiolipin synthetase, histidine kinase yycG, teichoic acid modification genes, and other genes involved in cell surface modification. Transmission electron microscopy showed that the cell walls of both strains J4206 and SA1 were significantly thicker and more electron dense than daptomycin- and vancomycin-sensitive strain J4211. This comparative genomic study has identified unique genes as well as allelic variants in the J4206 genome that are involved in cell surface modification and thus might contribute to the acquisition of daptomycin resistance.

Comparative genomics and physiology of the butyrate-producing bacterium Intestinimonas butyriciproducens

Intestinimonas is a newly described bacterial genus with representative strains present in the intestinal tract of human and other animals. Despite unique metabolic features including the production of butyrate from both sugars and amino acids, there is to date no data on their diversity, ecology, and physiology. Using a comprehensive phylogenetic approach, Intestinimomas was found to include at least three species that colonize primarily the human and mouse intestine. We focused on the most common and cultivable species of the genus, Intestinimonas butyriciproducens, and performed detailed genomic and physiological comparison of strains SRB521^T and AF211, isolated from the mouse and human gut respectively. The complete 3.3-Mb genomic sequences of both strains were highly similar with 98.8% average nucleotide identity, testifying to their assignment to one single species. However, thorough analysis revealed significant genomic rearrangements, variations in phage-derived sequences, and the presence of new CRISPR sequences in both strains. Moreover, strain AF211 appeared to be more efficient than strain SRB521^T in the conversion of the sugars arabinose and galactose. In conclusion, this study provides genomic and physiological insight into Intestinimonas butyriciproducens, a prevalent butyrate-producing species, differentiating strains that originate from the mouse and human gut.

Two draft genome sequences of Pseudomonas jessenii strains isolated from a copper contaminated site in Denmark

Pseudomonas jessenii C2 and Pseudomonas jessenii H16 were isolated from low-Cu and high-Cu industrially contaminated soil, respectively. P. jessenii H16 displayed significant resistance to copper when compared to P. jessenii C2. Here we describe genome sequences and interesting features of these two strains. The genome of P. jessenii C2 comprised 6,420,113 bp, with 5814 protein-coding genes and 67 RNA genes. P. jessenii H16 comprised 6,807,788 bp, with 5995 protein-coding genes and 70 RNA genes. Of special interest was a specific adaptation to this harsh copper-contaminated environment as P. jessenii H16 contained a novel putative copper resistance genomic island (GI) of around 50,000 bp.

Unusual Legionnaires' outbreak in cool, dry Western Canada: an investigation using genomic epidemiology

An outbreak of Legionnaires' disease occurred in an inner city district in Calgary, Canada. This outbreak spanned a 3-week period in November-December 2012, and a total of eight cases were identified. Four of these cases were critically ill requiring intensive care admission but there was no associated mortality. All cases tested positive for Legionella pneumophila serogroup 1 (LP1) by urinary antigen testing. Five of the eight patients were culture positive for LP1 from respiratory specimens. These isolates were further identified as Knoxville monoclonal subtype and sequence subtype ST222. Whole-genome sequencing revealed that the isolates differed by no more than a single vertically acquired single nucleotide variant, supporting a single point-source outbreak. Hypothesis-based environmental investigation and sampling was conducted; however, a definitive source was not identified. Geomapping of case movements within the affected urban sector revealed a 1·0 km common area of potential exposure, which coincided with multiple active construction sites that used water spray to minimize transient dust. This community point-source Legionnaires' disease outbreak is unique due to its ST222 subtype and occurrence in a relatively dry and cold weather setting in Western Canada. This report suggests community outbreaks of Legionella should not be overlooked as a possibility during late autumn and winter months in the Northern Hemisphere.

Complete genome sequence of Bacillus amyloliquefaciens subsp. plantarum S499, a rhizobacterium that triggers plant defences and inhibits fungal phytopathogens

Bacillus amyloliquefaciens subsp. plantarum S499 is a plant beneficial rhizobacterium with a good antagonistic potential against phytopathogens through the release of active secondary metabolites. Moreover, it can induce systemic resistance in plants by producing considerable amounts of surfactins. The complete genome sequence of B. amyloliquefaciens subsp. plantarum S499 includes a circular chromosome of 3,927,922bp and a plasmid of 8,008bp. A remarkable abundance in genomic regions of putative horizontal origin emerged from the analysis. Furthermore, we highlighted the presence of genes involved in the establishment of interactions with the host plants at the root level and in the competition with other soil-borne microorganisms. More specifically, genes related to the synthesis of amylolysin, amylocyclicin, and butirosin were identified. These antimicrobials were not known before to be part of the antibiotic arsenal of the strain. The information embedded in the genome will support the upcoming studies regarding the application of B. amyloliquefaciens isolates as plant-growth promoters and biocontrol agents.

A Primer on Infectious Disease Bacterial Genomics

The number of large-scale genomics projects is increasing due to the availability of affordable high-throughput sequencing (HTS) technologies. The use of HTS for bacterial infectious disease research is attractive because one whole-genome sequencing (WGS) run can replace multiple assays for bacterial typing, molecular epidemiology investigations, and more in-depth pathogenomic studies. The computational resources and bioinformatics expertise required to accommodate and analyze the large amounts of data pose new challenges for researchers embarking on genomics projects for the first time. Here, we present a comprehensive overview of a bacterial genomics projects from beginning to end, with a particular focus on the planning and computational requirements for HTS data, and provide a general understanding of the analytical concepts to develop a workflow that will meet the objectives and goals of HTS projects.

Complete genome sequences of the Serratia plymuthica strains 3Rp8 and 3Re4-18, two rhizosphere bacteria with antagonistic activity towards fungal phytopathogens and plant growth promoting abilities

The Serratia plymuthica strains 3Rp8 and 3Re4-18 are motile, Gram-negative, non-sporulating bacteria. Strain 3Rp8 was isolated from the rhizosphere of Brassica napus L. and strain 3Re4-18 from the endorhiza of Solanum tuberosum L. Studies have shown in vitro activity against the soil-borne fungi Verticillium dahliae Kleb., Rhizoctonia solani Kühn, and Sclerotinia sclerotiorum. Here, we announce and describe the complete genome sequence of S. plymuthica 3Rp8 consisting of a single circular chromosome of 5.5 Mb that encodes 4954 protein-coding and 108 RNA-only encoding genes and of S. plymuthica 3Re4-18 consisting of a single circular chromosome of 5.4 Mb that encodes 4845 protein-coding and 109 RNA-only encoding genes. The whole genome sequences and annotations are available in NCBI under the locus numbers CP012096 and CP012097, respectively. The genome analyses revealed genes putatively responsible for the promising plant growth promoting and biocontrol properties including predicting factors such as secretion systems, iron scavenging siderophores, chitinases, secreted proteases, glucanases and non-ribosomal peptide synthetases, as well as unique genomic islands.

The early events underlying genome evolution in a localized Sinorhizobium meliloti population

BACKGROUND

Population genetic analyses based on genome-wide sequencing data have been carried out for Sinorhizobium medicae and S. meliloti, two closely related bacterial species forming nitrogen-fixing symbioses with plants of the genus Medicago. However, genome coverage was low or the isolates had a broad geographic distribution, making it difficult to interpret the estimated diversity and to unravel the early events underlying population genetic variations and ecological differentiation.

RESULTS

Here, to gain insight into the early genome level variation and diversification within S. meliloti populations, we first used Illumina paired-end reads technology to sequence a new clone of S. meliloti strain GR4, a highly competitive strain for alfalfa nodulation. The Illumina data and the GR4 genome sequence previously obtained with 454 technology were used to generate a high-quality reference genome sequence. We then used Illumina technology to sequence the genomes of 13 S. meliloti isolates representative of the genomic variation within the GR4-type population, obtained from a single field site with a high degree of coverage. The genome sequences obtained were analyzed to determine nucleotide diversity, divergence times, polymorphism and genomic variation. Similar low levels of nucleotide diversity were observed for the chromosome, pSymB and pSymA replicons. The isolates displayed other types of variation, such as indels, recombination events, genomic island excision and the transposition of mobile elements.

CONCLUSIONS

Our results suggest that the GR4-type population has experienced a process of demographic expansion and behaves as a stable genotypic cluster of genome-wide similarity, with most of the genome following a clonal pattern of evolution. Although some of genetic variation detected within the GR4-type population is probably due to genetic drift, others might be important in diversification and environmental adaptation.

Structural and functional analysis of the finished genome of the recently isolated toxic Anabaena sp. WA102

Background

Very few closed genomes of the cyanobacteria that commonly produce toxic blooms in lakes and reservoirs are available, limiting our understanding of the properties of these organisms. A new anatoxin-a-producing member of the Nostocaceae, Anabaena sp. WA102, was isolated from a freshwater lake in Washington State, USA, in 2013 and maintained in non-axenic culture.

Results

The Anabaena sp. WA102 5.7 Mbp genome assembly has been closed with long-read, single-molecule sequencing and separately a draft genome assembly has been produced with short-read sequencing technology. The closed and draft genome assemblies are compared, showing a correlation between long repeats in the genome and the many gaps in the short-read assembly. Anabaena sp. WA102 encodes anatoxin-a biosynthetic genes, as does its close relative Anabaena sp. AL93 (also introduced in this study). These strains are distinguished by differences in the genes for light-harvesting phycobilins, with Anabaena sp. AL93 possessing a phycoerythrocyanin operon. Biologically relevant structural variants in the Anabaena sp. WA102 genome were detected only by long-read sequencing: a tandem triplication of the anaBCD promoter region in the anatoxin-a synthase gene cluster (not triplicated in Anabaena sp. AL93) and a 5-kbp deletion variant present in two-thirds of the population. The genome has a large number of mobile elements (160). Strikingly, there was no synteny with the genome of its nearest fully assembled relative, Anabaena sp. 90.

Conclusion

Structural and functional genome analyses indicate that Anabaena sp. WA102 has a flexible genome. Genome closure, which can be readily achieved with long-read sequencing, reveals large scale (e.g., gene order) and local structural features that should be considered in understanding genome evolution and function.

Electronic supplementary material

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

Computational methods for predicting genomic islands in microbial genomes

Clusters of genes acquired by lateral gene transfer in microbial genomes, are broadly referred to as genomic islands (GIs). GIs often carry genes important for genome evolution and adaptation to niches, such as genes involved in pathogenesis and antibiotic resistance. Therefore, GI prediction has gradually become an important part of microbial genome analysis. Despite inherent difficulties in identifying GIs, many computational methods have been developed and show good performance. In this mini-review, we first summarize the general challenges in predicting GIs. Then we group existing GI detection methods by their input, briefly describe representative methods in each group, and discuss their advantages as well as limitations. Finally, we look into the potential improvements for better GI prediction.

Complete Genome Sequence of Salmonella enterica Serovar Typhimurium Strain SO3 (Sequence Type 302) Isolated from a Baby with Meningitis in Mexico

The complete genome of ITALIC! Salmonella entericaserovar Typhimurium strain SO3 (sequence type 302), isolated from a fatal meningitis infection in Mexico, was determined using PacBio technology. The chromosome hosts six complete prophages and is predicted to harbor 51 genomic islands, including 13 pathogenicity islands (SPIs). It carries the ITALIC! Salmonellavirulence plasmid (pSTV).

Complete Genome Sequence of Salmonella enterica Serovar Typhimurium Strain SO2 (Sequence Type 302) Isolated from an Asymptomatic Child in Mexico

The complete genome sequence of Salmonella enterica serovar Typhimurium strain SO2, isolated from an asymptomatic child in Mexico, was determined using PacBio single-molecule real-time technology. Strain SO2 has six complete chromosomal prophages, namely, ST104, Gifsy-2, ST64B, Gifsy-1, ELPhiS, and FSL SP-004, and carries a Salmonella virulence plasmid.

Complete Genome Sequence of Salmonella enterica Serovar Typhimurium Strain YU15 (Sequence Type 19) Harboring the Salmonella Genomic Island 1 and Virulence Plasmid pSTV

The complete genome ofSalmonella entericasubsp.entericaserovar Typhimurium sequence type 19 (ST19) strain YU15, isolated in Yucatán, Mexico, from a human baby stool culture, was determined using PacBio technology. The chromosome contains five intact prophages and theSalmonellagenomic island 1 (SGI1). This strain carries theSalmonellavirulence plasmid pSTV.