Score-based prediction of genomic islands in prokaryotic genomes using hidden Markov models

Comparative Genomics of Degradative Novosphingobium Strains With Special Reference to Microcystin-Degrading Novosphingobium sp. THN1

Bacteria in genus Novosphingobium associated with biodegradation of substrates are prevalent in environments such as lakes, soil, sea, wood and sediments. To better understand the characteristics linked to their wide distribution and metabolic versatility, we report the whole genome sequence of Novosphingobium sp. THN1, a microcystin-degrading strain previously isolated by Jiang et al. (2011) from cyanobacteria-blooming water samples from Lake Taihu, China. We performed a genomic comparison analysis of Novosphingobium sp. THN1 with 21 other degradative Novosphingobium strains downloaded from GenBank. Phylogenetic trees were constructed using 16S rRNA genes, core genes, protein-coding sequences, and average nucleotide identity of whole genomes. Orthologous protein analysis showed that the 22 genomes contained 674 core genes and each strain contained a high proportion of distributed genes that are shared by a subset of strains. Inspection of their genomic plasticity revealed a high number of insertion sequence elements and genomic islands that were distributed on both chromosomes and plasmids. We also compared the predicted functional profiles of the Novosphingobium protein-coding genes. The flexible genes and all protein-coding genes produced the same heatmap clusters. The COG annotations were used to generate a dendrogram correlated with the compounds degraded. Furthermore, the metabolic profiles predicted from KEGG pathways showed that the majority of genes involved in central carbon metabolism, nitrogen, phosphate, sulfate metabolism, energy metabolism and cell mobility (above 62.5%) are located on chromosomes. Whereas, a great many of genes involved in degradation pathways (21-50%) are located on plasmids. The abundance and distribution of aromatics-degradative mono- and dioxygenases varied among 22 Novosphingoibum strains. Comparative analysis of the microcystin-degrading mlr gene cluster provided evidence for horizontal acquisition of this cluster. The Novosphingobium sp. THN1 genome sequence contained all the functional genes crucial for microcystin degradation and the mlr gene cluster shared high sequence similarity (≥85%) with the sequences of other microcystin-degrading genera isolated from cyanobacteria-blooming water. Our results indicate that Novosphingobium species have high genomic and functional plasticity, rearranging their genomes according to environment variations and shaping their metabolic profiles by the substrates they are exposed to, to better adapt to their environments.

Involvement of β-Carbonic Anhydrase Genes in Bacterial Genomic Islands and Their Horizontal Transfer to Protists

Genomic islands (GIs) are a type of mobile genetic element (MGE) that are present in bacterial chromosomes. They consist of a cluster of genes that produce proteins that contribute to a variety of functions, including, but not limited to, the regulation of cell metabolism, antimicrobial resistance, pathogenicity, virulence, and resistance to heavy metals. The genes carried in MGEs can be used as a trait reservoir in times of adversity. Transfer of genes using MGEs, occurring outside reproduction, is called horizontal gene transfer (HGT). Previous data have shown that numerous HGT events have occurred through endosymbiosis between prokaryotes and eukaryotes. β-Carbonic anhydrase (β-CA) enzymes play a critical role in the biochemical pathways of many prokaryotes and eukaryotes. We previously suggested the horizontal transfer of β-CA genes from plasmids of some prokaryotic endosymbionts to their protozoan hosts. In this study, we set out to identify β-CA genes that might have been transferred between prokaryotic and protist species through HGT in GIs. Therefore, we investigated prokaryotic chromosomes containing β-CA-encoding GIs and utilized multiple bioinformatics tools to reveal the distinct movements of β-CA genes among a wide variety of organisms. Our results identify the presence of β-CA genes in GIs of several medically and industrially relevant bacterial species, and phylogenetic analyses reveal multiple cases of likely horizontal transfer of β-CA genes from GIs of ancestral prokaryotes to protists.IMPORTANCE The evolutionary process is mediated by mobile genetic elements (MGEs), such as genomic islands (GIs). A gene or set of genes in the GIs is exchanged between and within various species through horizontal gene transfer (HGT). Based on the crucial role that GIs can play in bacterial survival and proliferation, they were introduced as environment- and pathogen-associated factors. Carbonic anhydrases (CAs) are involved in many critical biochemical pathways, such as the regulation of pH homeostasis and electrolyte transfer. Among the six evolutionary families of CAs, β-CA gene sequences are present in many bacterial species, which can be horizontally transferred to protists during evolution. This study shows the involvement of bacterial β-CA gene sequences in the GIs and suggests their horizontal transfer to protists during evolution.

Deciphering pathogenicity and antibiotic resistance islands in methicillin-resistant Staphylococcus aureus genomes

Staphylococcus aureus is a versatile pathogen that is capable of causing infections in both humans and animals. It can cause furuncles, septicaemia, pneumonia and endocarditis. Adaptation of S. aureus to the modern hospital environment has been facilitated, in part, by the horizontal acquisition of drug resistance genes, such as mecA gene that imparts resistance to methicillin. Horizontal acquisitions of islands of genes harbouring virulence and antibiotic resistance genes have made S. aureus resistant to commonly used antibiotics. To decipher genomic islands (GIs) in 22 hospital- and 9 community-associated methicillin-resistant S. aureus strains and classify a subset of GIs carrying virulence and resistance genes as pathogenicity and resistance islands respectively, we applied a host of methods for localizing genomic islands in prokaryotic genomes. Surprisingly, none of the frequently used GI prediction methods could perform well in delineating the resistance islands in the S. aureus genomes. Rather, a gene clustering procedure exploiting biases in codon usage for identifying horizontally transferred genes outperformed the current methods for GI detection, in particular in identifying the known islands in S. aureus including the SCCmec island that harbours the mecA resistance gene. The gene clustering approach also identified novel, as yet unreported islands, with many of these found to harbour virulence and/or resistance genes. These as yet unexplored islands may provide valuable information on the evolution of drug resistance in S. aureus.

GI-Cluster: Detecting genomic islands via consensus clustering on multiple features

The accurate detection of genomic islands (GIs) in microbial genomes is important for both evolutionary study and medical research, because GIs may promote genome evolution and contain genes involved in pathogenesis. Various computational methods have been developed to predict GIs over the years. However, most of them cannot make full use of GI-associated features to achieve desirable performance. Additionally, many methods cannot be directly applied to newly sequenced genomes. We develop a new method called GI-Cluster, which provides an effective way to integrate multiple GI-related features via consensus clustering. GI-Cluster does not require training datasets or existing genome annotations, but it can still achieve comparable or better performance than supervised learning methods in comprehensive evaluations. Moreover, GI-Cluster is widely applicable, either to complete and incomplete genomes or to initial GI predictions from other programs. GI-Cluster also provides plots to visualize the distribution of predicted GIs and related features. GI-Cluster is available at https://github.com/icelu/GI_Cluster.

Genus-wide comparison of Pseudovibrio bacterial genomes reveal diverse adaptations to different marine invertebrate hosts

Bacteria belonging to the genus Pseudovibrio have been frequently found in association with a wide variety of marine eukaryotic invertebrate hosts, indicative of their versatile and symbiotic lifestyle. A recent comparison of the sponge-associated Pseudovibrio genomes has shed light on the mechanisms influencing a successful symbiotic association with sponges. In contrast, the genomic architecture of Pseudovibrio bacteria associated with other marine hosts has received less attention. Here, we performed genus-wide comparative analyses of 18 Pseudovibrio isolated from sponges, coral, tunicates, flatworm, and seawater. The analyses revealed a certain degree of commonality among the majority of sponge- and coral-associated bacteria. Isolates from other marine invertebrate host, tunicates, exhibited a genetic repertoire for cold adaptation and specific metabolic abilities including mucin degradation in the Antarctic tunicate-associated bacterium Pseudovibrio sp. Tun.PHSC04_5.I4. Reductive genome evolution was simultaneously detected in the flatworm-associated bacteria and the sponge-associated bacterium P. axinellae AD2, through the loss of major secretion systems (type III/VI) and virulence/symbioses factors such as proteins involved in adhesion and attachment to the host. Our study also unraveled the presence of a CRISPR-Cas system in P. stylochi UST20140214-052 a flatworm-associated bacterium possibly suggesting the role of CRISPR-based adaptive immune system against the invading virus particles. Detection of mobile elements and genomic islands (GIs) in all bacterial members highlighted the role of horizontal gene transfer for the acquisition of novel genetic features, likely enhancing the bacterial ecological fitness. These findings are insightful to understand the role of genome diversity in Pseudovibrio as an evolutionary strategy to increase their colonizing success across a wide range of marine eukaryotic hosts.

A Gene Homologous to rRNA Methylase Genes Confers Erythromycin and Clindamycin Resistance in Bifidobacterium breve

Bifidobacteria are mutualistic intestinal bacteria, and their presence in the human gut has been associated with health-promoting activities. The presence of antibiotic resistance genes in this genus is controversial, since, although bifidobacteria are nonpathogenic microorganisms, they could serve as reservoirs of resistance determinants for intestinal pathogens. However, until now, few antibiotic resistance determinants have been functionally characterized in this genus. In this work, we show that Bifidobacterium breve CECT7263 displays atypical resistance to erythromycin and clindamycin. In order to delimit the genomic region responsible for the observed resistance phenotype, a library of genomic DNA was constructed and a fragment of 5.8 kb containing a gene homologous to rRNA methylase genes was able to confer erythromycin resistance in Escherichia coli This genomic region seems to be very uncommon, and homologs of the gene have been detected in only one strain of Bifidobacterium longum and two other strains of B. breve In this context, analysis of shotgun metagenomics data sets revealed that the gene is also uncommon in the microbiomes of adults and infants. The structural gene and its upstream region were cloned into a B. breve-sensitive strain, which became resistant after acquiring the genetic material. In vitro conjugation experiments did not allow us to detect gene transfer to other recipients. Nevertheless, prediction of genes potentially acquired through horizontal gene transfer events revealed that the gene is located in a putative genomic island.IMPORTANCEBifidobacterium breve is a very common human intestinal bacterium. Often described as a pioneer microorganism in the establishment of early-life intestinal microbiota, its presence has been associated with several beneficial effects for the host, including immune stimulation and protection against infections. Therefore, some strains of this species are considered probiotics. In relation to this, because probiotic bacteria are used for human and animal consumption, one of the safety concerns over these bacteria is the presence of antibiotic resistance genes, since the human gut is a densely populated habitat that could favor the transfer of genetic material to potential pathogens. In this study, we analyzed the genetic basis responsible for the erythromycin and clindamycin resistance phenotype of B. breve CECT7263. We were able to identify and characterize a novel gene homologous to rRNA methylase genes which confers erythromycin and clindamycin resistance. This gene seems to be very uncommon in other bifidobacteria and in the gut microbiomes of both adults and infants. Even though conjugation experiments showed the absence of transferability under in vitro conditions, it has been predicted to be located in a putative genomic island recently acquired by specific bifidobacterial strains.

Pan-Genome Analysis Links the Hereditary Variation of Leptospirillum ferriphilum With Its Evolutionary Adaptation

Niche adaptation has long been recognized to drive intra-species differentiation and speciation, yet knowledge about its relatedness with hereditary variation of microbial genomes is relatively limited. Using Leptospirillum ferriphilum species as a case study, we present a detailed analysis of genomic features of five recognized strains. Genome-to-genome distance calculation preliminarily determined the roles of spatial distance and environmental heterogeneity that potentially contribute to intra-species variation within L. ferriphilum species at the genome level. Mathematical models were further constructed to extrapolate the expansion of L. ferriphilum genomes (an 'open' pan-genome), indicating the emergence of novel genes with new sequenced genomes. The identification of diverse mobile genetic elements (MGEs) (such as transposases, integrases, and phage-associated genes) revealed the prevalence of horizontal gene transfer events, which is an important evolutionary mechanism that provides avenues for the recruitment of novel functionalities and further for the genetic divergence of microbial genomes. Comprehensive analysis also demonstrated that the genome reduction by gene loss in a broad sense might contribute to the observed diversification. We thus inferred a plausible explanation to address this observation: the community-dependent adaptation that potentially economizes the limiting resources of the entire community. Now that the introduction of new genes is accompanied by a parallel abandonment of some other ones, our results provide snapshots on the biological fitness cost of environmental adaptation within the L. ferriphilum genomes. In short, our genome-wide analyses bridge the relation between genetic variation of L. ferriphilum with its evolutionary adaptation.

Resolving the complete genome of Kuenenia stuttgartiensis from a membrane bioreactor enrichment using Single-Molecule Real-Time sequencing

Anaerobic ammonium-oxidizing (anammox) bacteria are a group of strictly anaerobic chemolithoautotrophic microorganisms. They are capable of oxidizing ammonium to nitrogen gas using nitrite as a terminal electron acceptor, thereby facilitating the release of fixed nitrogen into the atmosphere. The anammox process is thought to exert a profound impact on the global nitrogen cycle and has been harnessed as an environment-friendly method for nitrogen removal from wastewater. In this study, we present the first closed genome sequence of an anammox bacterium, Kuenenia stuttgartiensis MBR1. It was obtained through Single-Molecule Real-Time (SMRT) sequencing of an enrichment culture constituting a mixture of at least two highly similar Kuenenia strains. The genome of the novel MBR1 strain is different from the previously reported Kuenenia KUST reference genome as it contains numerous structural variations and unique genomic regions. We find new proteins, such as a type 3b (sulf)hydrogenase and an additional copy of the hydrazine synthase gene cluster. Moreover, multiple copies of ammonium transporters and proteins regulating nitrogen uptake were identified, suggesting functional differences in metabolism. This assembly, including the genome-wide methylation profile, provides a new foundation for comparative and functional studies aiming to elucidate the biochemical and metabolic processes of these organisms.

Comparative Genomics of Environmental and Clinical Burkholderia cenocepacia Strains Closely Related to the Highly Transmissible Epidemic ET12 Lineage

The Burkholderia cenocepacia epidemic ET12 lineage belongs to the genomovar IIIA including the reference strain J2315, a highly transmissible epidemic B. cenocepacia lineage. Members of this lineage are able to cause lung infections in immunocompromised and cystic fibrosis patients. In this study, we describe the genome of F01, an environmental B. cenocepacia strain isolated from soil in Burkina Faso that is, to our knowledge, the most closely related strain to this epidemic lineage. A comparative genomic analysis was performed on this new isolate, in association with five clinical and one environmental B. cenocepacia strains whose genomes were previously sequenced. Antibiotic resistances, virulence phenotype, and genomic contents were compared and discussed with an emphasis on virulent and antibiotic determinants. Surprisingly, no significant differences in antibiotic resistance and virulence were found between clinical and environmental strains, while the most important genomic differences were related to the number of prophages identified in their genomes. The ET12 lineage strains showed a noticeable greater number of prophages (partial or full-length), especially compared to the phylogenetically related environmental F01 strain (i.e., 5–6 and 3 prophages, respectively). Data obtained suggest possible involvements of prophages in the clinical success of opportunistic pathogens.

Comparative genomic analysis reveals the evolution and environmental adaptation strategies of vibrios

BACKGROUND

Vibrios are among the most diverse and ecologically important marine bacteria, which have evolved many characteristics and lifestyles to occupy various niches. The relationship between genome features and environmental adaptation strategies is an essential part for understanding the ecological functions of vibrios in the marine system. The advent of complete genome sequencing technology has provided an important method of examining the genetic characteristics of vibrios on the genomic level.

RESULTS

Two Vibrio genomes were sequenced and found to occupy many unique orthologues families which absent from the previously genes pool of the complete genomes of vibrios. Comparative genomics analysis found vibrios encompass a steady core-genome and tremendous pan-genome with substantial gene gain and horizontal gene transfer events in the evolutionary history. Evolutionary analysis based on the core-genome tree suggested that V. fischeri emerged ~ 385 million years ago, along with the occurrence of cephalopods and the flourish of fish. The relatively large genomes, the high number of 16S rRNA gene copies, and the presence of R-M systems and CRISPR system help vibrios live in various marine environments. Chitin-degrading related genes are carried in nearly all the Vibrio genomes. The number of chitinase genes in vibrios has been extremely expanded compared to which in the most recent ancestor of the genus. The chitinase A genes were estimated to have evolved along with the genus, and have undergone significant purifying selective force to conserve the ancestral state.

CONCLUSIONS

Vibrios have experienced extremely genome expansion events during their evolutionary history, allowing them to develop various functions to spread globally. Despite their close phylogenetic relationships, vibrios were found to have a tremendous pan-genome with a steady core-genome, which indicates the highly plastic genome of the genus. Additionally, the existence of various chitin-degrading related genes and the expansion of chitinase A in the genus demonstrate the importance of the chitin utilization for vibrios. Defensive systems in the Vibrio genomes may protect them from the invasion of external DNA. These genomic features investigated here provide a better knowledge of how the evolutionary process has forged Vibrio genomes to occupy various niches.

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.

Whole-genome sequencing reveals the mechanisms for evolution of streptomycin resistance in Lactobacillus plantarum

In this research, we investigated the evolution of streptomycin resistance in Lactobacillus plantarum ATCC14917, which was passaged in medium containing a gradually increasing concentration of streptomycin. After 25 d, the minimum inhibitory concentration (MIC) of L. plantarum ATCC14917 had reached 131,072 µg/mL, which was 8,192-fold higher than the MIC of the original parent isolate. The highly resistant L. plantarum ATCC14917 isolate was then passaged in antibiotic-free medium to determine the stability of resistance. The MIC value of the L. plantarum ATCC14917 isolate decreased to 2,048 µg/mL after 35 d but remained constant thereafter, indicating that resistance was irreversible even in the absence of selection pressure. Whole-genome sequencing of parent isolates, control isolates, and isolates following passage was used to study the resistance mechanism of L. plantarum ATCC14917 to streptomycin and adaptation in the presence and absence of selection pressure. Five mutated genes (single nucleotide polymorphisms and structural variants) were verified in highly resistant L. plantarum ATCC14917 isolates, which were related to ribosomal protein S12, LPXTG-motif cell wall anchor domain protein, LrgA family protein, Ser/Thr phosphatase family protein, and a hypothetical protein that may correlate with resistance to streptomycin. After passage in streptomycin-free medium, only the mutant gene encoding ribosomal protein S12 remained; the other 4 mutant genes had reverted to the wild type as found in the parent isolate. Although the MIC value of L. plantarum ATCC14917 was reduced in the absence of selection pressure, it remained 128-fold higher than the MIC value of the parent isolate, indicating that ribosomal protein S12 may play an important role in streptomycin resistance. Using the mobile elements database, we demonstrated that streptomycin resistance-related genes in L. plantarum ATCC14917 were not located on mobile elements. This research offers a way of combining laboratory evolution techniques and whole-genome sequencing for evaluating antibiotic resistance in probiotics.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

A Practical Guide for Comparative Genomics of Mobile Genetic Elements in Prokaryotic Genomes

Mobile genetic elements (MGEs) are an important feature of prokaryote genomes but are seldom well annotated and, consequently, are often underestimated. MGEs include transposons (Tn), insertion sequences (ISs), prophages, genomic islands (GEIs), integrons, and integrative and conjugative elements (ICEs). They are intimately involved in genome evolution and promote phenomena such as genomic expansion and rearrangement, emergence of virulence and pathogenicity, and symbiosis. In spite of the annotation bottleneck, there are so far at least 75 different programs and databases dedicated to prokaryotic MGE analysis and annotation, and this number is rapidly growing. Here, we present a practical guide to explore, compare, and visualize prokaryote MGEs using a combination of available software and databases tailored to small scale genome analyses. This protocol can be coupled with expert MGE annotation and exploited for evolutionary and comparative genomic analyses.

Adaptation in toxic environments: comparative genomics of loci carrying antibiotic resistance genes derived from acid mine drainage waters

Several studies have suggested the existence of a close relationship between antibiotic-resistant phenotypes and resistance to other toxic compounds such as heavy metals, which involve co-resistance or cross-resistance mechanisms. A metagenomic library was previously constructed in Escherichia coli with DNA extracted from the bacterial community inhabiting an acid mine drainage (AMD) site highly contaminated with heavy metals. Here, we conducted a search for genes involved in antibiotic resistance using this previously constructed library. In particular, resistance to antibiotics was observed among five clones carrying four different loci originating from CARN5 and CARN2, two genomes reconstructed from the metagenomic data. Among the three CARN2 loci, two carry genes homologous to those previously proposed to be involved in antibiotic resistance. The third CARN2 locus carries a gene encoding a membrane transporter with an unknown function and was found to confer bacterial resistance to rifampicin, gentamycin, and kanamycin. The genome of Thiomonas delicata DSM 16361 and Thiomonas sp. X19 were sequenced in this study. Homologs of genes carried on these three CARN2 loci were found in these genomes, two of these loci were found in genomic islands. Together, these findings confirm that AMD environments contaminated with several toxic metals also constitute habitats for bacteria that function as reservoirs for antibiotic resistance genes.

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/.

The genomic and transcriptomic basis of the potential of Lactobacillus plantarum A6 to improve the nutritional quality of a cereal based fermented food

The objective of this work was to investigate the nutritional potential of Lactobacillus plantarum A6 in a food matrix using next generation sequencing. To this end, we characterized the genome of the A6 strain for a complete overview of its potential. We then compared its transcriptome when grown in a food matrix made from pearl millet to and its transcriptome when cultivated in a laboratory medium. Genomic comparison of the strain L. plantarum A6 with the strains WCFS1, ST-III, JDM1 and ATCC14917 led to the identification of five regions of genomic plasticity. More specifically, 362 coding sequences, mostly annotated as coding for proteins of unknown functions, were specific to L. plantarum A6. A total of 1201 genes were significantly differentially expressed in laboratory medium and food matrix. Among them, 821 genes were up-regulated in the food matrix compared to the laboratory medium, representing 23% of whole genomic objects. In the laboratory medium, the expression of 380 genes, representing 11% of the all genomic objects was at least double than in the food matrix. Genes encoding important functions for the nutritional quality of the food were identified. Considering its efficiency as an amylolytic strain, we investigated all genes involved in carbohydrate metabolism, paying particular attention to starch metabolism. An extracellular alpha amylase, a neopullulanase and maltodextrin transporters were identified, all of which were highly expressed in the food matrix. In addition, genes involved in alpha-galactoside metabolism were identified but only two of them were induced in food matrix than in laboratory medium. This may be because alpha galactosides were already eliminated during soaking. Different biosynthetic pathways involved in the synthesis of vitamin B (folate, riboflavin, and cobalamin) were identified. They allowed the identification of a potential of vitamin synthesis, which should be confirmed through biochemical analysis in further work. Surprisingly, some genes involved in metabolism and bioaccessibility of iron were identified. They were related directly to the use of transport of iron, or indirectly to metabolism of polyphenols, responsible of iron chelation in the food. The combination of genomics and transcriptomics not only revealed previously undocumented nutritional properties of L. plantarum A6, but also documented the behaviour of this bacterium in food.

Genomic Analysis of Clavibacter michiganensis Reveals Insight Into Virulence Strategies and Genetic Diversity of a Gram-Positive Bacterial Pathogen

Clavibacter michiganensis subsp. michiganensis is a gram-positive bacterial pathogen that proliferates in the xylem vessels of tomato, causing bacterial canker disease. In this study, we sequenced and assembled genomes of 11 C. michiganensis subsp. michiganensis strains isolated from infected tomato fields in California as well as five Clavibacter strains that colonize tomato endophytically but are not pathogenic in this host. The analysis of the C. michiganensis subsp. michiganensis genomes supported the monophyletic nature of this pathogen but revealed genetic diversity among strains, consistent with multiple introduction events. Two tomato endophytes that clustered phylogenetically with C. michiganensis strains capable of infecting wheat and pepper and were also able to cause disease in these plants. Plasmid profiles of the California strains were variable and supported the essential role of the pCM1-like plasmid and the CelA cellulase in virulence, whereas the absence of the pCM2-like plasmid in some pathogenic C. michiganensis subsp. michiganensis strains revealed it is not essential. A large number of secreted C. michiganensis subsp. michiganensis proteins were carbohydrate-active enzymes (CAZymes). Glycome profiling revealed that C. michiganensis subsp. michiganensis but not endophytic Clavibacter strains is able to extensively alter tomato cell-wall composition. Two secreted CAZymes found in all C. michiganensis subsp. michiganensis strains, CelA and PelA1, enhanced pathogenicity on tomato. Collectively, these results provide a deeper understanding of C. michiganensis subsp. michiganensis diversity and virulence strategies.

Genomic and Transcriptomic Analysis of Growth-Supporting Dehalogenation of Chlorinated Methanes in Methylobacterium

Bacterial adaptation to growth with toxic halogenated chemicals was explored in the context of methylotrophic metabolism of Methylobacterium extorquens, by comparing strains CM4 and DM4, which show robust growth with chloromethane and dichloromethane, respectively. Dehalogenation of chlorinated methanes initiates growth-supporting degradation, with intracellular release of protons and chloride ions in both cases. The core, variable and strain-specific genomes of strains CM4 and DM4 were defined by comparison with genomes of non-dechlorinating strains. In terms of gene content, adaptation toward dehalogenation appears limited, strains CM4 and DM4 sharing between 75 and 85% of their genome with other strains of M. extorquens. Transcript abundance in cultures of strain CM4 grown with chloromethane and of strain DM4 grown with dichloromethane was compared to growth with methanol as a reference C₁ growth substrate. Previously identified strain-specific dehalogenase-encoding genes were the most transcribed with chlorinated methanes, alongside other genes encoded by genomic islands (GEIs) and plasmids involved in growth with chlorinated compounds as carbon and energy source. None of the 163 genes shared by strains CM4 and DM4 but not by other strains of M. extorquens showed higher transcript abundance in cells grown with chlorinated methanes. Among the several thousand genes of the M. extorquens core genome, 12 genes were only differentially abundant in either strain CM4 or strain DM4. Of these, 2 genes of known function were detected, for the membrane-bound proton translocating pyrophosphatase HppA and the housekeeping molecular chaperone protein DegP. This indicates that the adaptive response common to chloromethane and dichloromethane is limited at the transcriptional level, and involves aspects of the general stress response as well as of a dehalogenation-specific response to intracellular hydrochloric acid production. Core genes only differentially abundant in either strain CM4 or strain DM4 total 13 and 58 CDS, respectively. Taken together, the obtained results suggest different transcriptional responses of chloromethane- and dichloromethane-degrading M. extorquens strains to dehalogenative metabolism, and substrate- and pathway-specific modes of growth optimization with chlorinated methanes.

The Complete Genome Sequence of the Fish Pathogen Tenacibaculum maritimum Provides Insights into Virulence Mechanisms

Tenacibaculum maritimum is a devastating bacterial pathogen of wild and farmed marine fish with a broad host range and a worldwide distribution. We report here the complete genome sequence of the T. maritimum type strain NCIMB 2154^T. The genome consists of a 3,435,971-base pair circular chromosome with 2,866 predicted protein-coding genes. Genes encoding the biosynthesis of exopolysaccharides, the type IX secretion system, iron uptake systems, adhesins, hemolysins, proteases, and glycoside hydrolases were identified. They are likely involved in the virulence process including immune escape, invasion, colonization, destruction of host tissues, and nutrient scavenging. Among the predicted virulence factors, type IX secretion-mediated and cell-surface exposed proteins were identified including an atypical sialidase, a sphingomyelinase and a chondroitin AC lyase which activities were demonstrated in vitro.

Comparative Analysis of the Flavobacterium columnare Genomovar I and II Genomes

Columnaris disease caused by Gram-negative rod Flavobacterium columnare is one of the most common diseases of catfish. F. columnare is also a common problem in other cultured fish species worldwide. F. columnare has three major genomovars; we have sequenced a representative strain from genomovar I (ATCC 49512, which is avirulent in catfish) and genomovar II (94-081, which is highly pathogenic in catfish). Here, we present a comparative analysis of the two genomes. Interestingly, F. columnare ATCC 49512 and 94-081 meet criteria to be considered different species based on the Average Nucleotide Identity (90.71% similar) and DNA–DNA Hybridization (42.6% similar). Genome alignment indicated the two genomes have a large number of rearrangements. However, function-based comparative genomics analysis indicated that the two strains have similar functional capabilities with 2,263 conserved orthologous clusters; strain ATCC 49512 has 290 unique orthologous clusters while strain 94-081 has 391. Both strains carry type I secretion system, type VI secretion system, and type IX secretion system. The two genomes also have similar CRISPR capacities. The F. columnare strain ATCC 49512 genome contains a higher number of insertion sequence families and phage regions, while the F. columnare strain 94-081 genome has more genomic islands and more regulatory gene capacity. Transposon mutagenesis using Tn4351 in pathogenic strain 94-081 yielded six mutants, and experimental infections of fish showed hemolysin and glycine cleavage protein mutants had 15 and 10% mortalities, respectively, while the wild-type strain caused 100% mortalities. Our comparative and mutational analysis yielded important information on classification of genomovars I and II F. columnare as well as potential virulence genes in F. columnare strain 94-081.

Rearranged Biosynthetic Gene Cluster and Synthesis of Hassallidin E in Planktothrix serta PCC 8927

Cyanobacteria produce a wide range of natural products with antifungal bioactivity. The cyclic glycosylated lipopeptides of the hassallidin family have potent antifungal activity and display a great degree of chemical diversity. Here, we report the discovery of a hassallidin biosynthetic gene cluster from the filamentous cyanobacterium Planktothrix serta PCC 8927. The hassallidin gene cluster showed heavy rearrangement and marks of genomic plasticity. Nucleotide bias, differences in GC content, and phylogenetic incongruence suggested the acquisition of the hassallidin biosynthetic gene cluster in Planktothrix serta PCC 8927 by horizontal gene transfer. Chemical analyses by liquid chromatography and mass spectrometry demonstrated that this strain produced hassallidin E, a new glycosylated hassallidin variant. Hassallidin E was the only structural variant produced by Planktothrix serta PCC 8927 in all tested conditions. Further evaluated on human pathogenic fungi, hassallidin E showed an antifungal bioactivity. Hassallidin production levels correlated with nitrogen availability, in the only nitrogen-fixing Planktothrix described so far. Our results provide insights into the distribution and chemical diversity of cyanobacterial antifungal compounds as well as raise questions on their ecological relevance.

'ARMAN' archaea depend on association with euryarchaeal host in culture and in situ

Intriguing, yet uncultured ‘ARMAN’-like archaea are metabolically dependent on other members of the microbial community. It remains uncertain though which hosts they rely upon, and, because of the lack of complete genomes, to what extent. Here, we report the co-culturing of ARMAN-2-related organism, Mia14, with Cuniculiplasma divulgatum PM4 during the isolation of this strain from acidic streamer in Parys Mountain (Isle of Anglesey, UK). Mia14 is highly enriched in the binary culture (ca. 10% genomic reads) and its ungapped 0.95 Mbp genome points at severe voids in central metabolic pathways, indicating dependence on the host, C. divulgatum PM4. Analysis of C. divulgatum isolates from different sites and shotgun sequence data of Parys Mountain samples suggests an extensive genetic exchange between Mia14 and hosts in situ. Within the subset of organisms with high-quality genomic assemblies representing the ‘DPANN’ superphylum, the Mia14 lineage has had the largest gene flux, with dozens of genes gained that are implicated in the host interaction.

In the absence of complete genomes, the metabolic capabilities of uncultured ARMAN-like archaea have been uncertain. Here, Golyshina et al. apply an enrichment culture technique and find that the ungapped genome of the ARMAN-like archaeon Mia14 has lost key metabolic pathways, suggesting dependence on the host archaeon Cuniculiplasma divulgatum.

Comparative Genome Analysis Provides Insights into Both the Lifestyle of Acidithiobacillus ferrivorans Strain CF27 and the Chimeric Nature of the Iron-Oxidizing Acidithiobacilli Genomes

The iron-oxidizing species Acidithiobacillus ferrivorans is one of few acidophiles able to oxidize ferrous iron and reduced inorganic sulfur compounds at low temperatures (<10°C). To complete the genome of At. ferrivorans strain CF27, new sequences were generated, and an update assembly and functional annotation were undertaken, followed by a comparative analysis with other Acidithiobacillus species whose genomes are publically available. The At. ferrivorans CF27 genome comprises a 3,409,655 bp chromosome and a 46,453 bp plasmid. At. ferrivorans CF27 possesses genes allowing its adaptation to cold, metal(loid)-rich environments, as well as others that enable it to sense environmental changes, allowing At. ferrivorans CF27 to escape hostile conditions and to move toward favorable locations. Interestingly, the genome of At. ferrivorans CF27 exhibits a large number of genomic islands (mostly containing genes of unknown function), suggesting that a large number of genes has been acquired by horizontal gene transfer over time. Furthermore, several genes specific to At. ferrivorans CF27 have been identified that could be responsible for the phenotypic differences of this strain compared to other Acidithiobacillus species. Most genes located inside At. ferrivorans CF27-specific gene clusters which have been analyzed were expressed by both ferrous iron-grown and sulfur-attached cells, indicating that they are not pseudogenes and may play a role in both situations. Analysis of the taxonomic composition of genomes of the Acidithiobacillia infers that they are chimeric in nature, supporting the premise that they belong to a particular taxonomic class, distinct to other proteobacterial subgroups.

Comparative genome and transcriptome analysis reveals distinctive surface characteristics and unique physiological potentials of Pseudomonas aeruginosa ATCC 27853

BACKGROUND

Pseudomonas aeruginosa ATCC 27853 was isolated from a hospital blood specimen in 1971 and has been widely used as a model strain to survey antibiotics susceptibilities, biofilm development, and metabolic activities of Pseudomonas spp.. Although four draft genomes of P. aeruginosa ATCC 27853 have been sequenced, the complete genome of this strain is still lacking, hindering a comprehensive understanding of its physiology and functional genome.

RESULTS

Here we sequenced and assembled the complete genome of P. aeruginosa ATCC 27853 using the Pacific Biosciences SMRT (PacBio) technology and Illumina sequencing platform. We found that accessory genes of ATCC 27853 including prophages and genomic islands (GIs) mainly contribute to the difference between P. aeruginosa ATCC 27853 and other P. aeruginosa strains. Seven prophages were identified within the genome of P. aeruginosa ATCC 27853. Of the predicted 25 GIs, three contain genes that encode monoxoygenases, dioxygenases and hydrolases that could be involved in the metabolism of aromatic compounds. Surveying virulence-related genes revealed that a series of genes that encode the B-band O-antigen of LPS are lacking in ATCC 27853. Distinctive SNPs in genes of cellular adhesion proteins such as type IV pili and flagella biosynthesis were also observed in this strain. Colony morphology analysis confirmed an enhanced biofilm formation capability of ATCC 27853 on solid agar surface compared to Pseudomonas aeruginosa PAO1. We then performed transcriptome analysis of ATCC 27853 and PAO1 using RNA-seq and compared the expression of orthologous genes to understand the functional genome and the genomic details underlying the distinctive colony morphogenesis. These analyses revealed an increased expression of genes involved in cellular adhesion and biofilm maturation such as type IV pili, exopolysaccharide and electron transport chain components in ATCC 27853 compared with PAO1. In addition, distinctive expression profiles of the virulence genes lecA, lasB, quorum sensing regulators LasI/R, and the type I, III and VI secretion systems were observed in the two strains.

CONCLUSIONS

The complete genome sequence of P. aeruginosa ATCC 27853 reveals the comprehensive genetic background of the strain, and provides genetic basis for several interesting findings about the functions of surface associated proteins, prophages, and genomic islands. Comparative transcriptome analysis of P. aeruginosa ATCC 27853 and PAO1 revealed several classes of differentially expressed genes in the two strains, underlying the genetic and molecular details of several known and yet to be explored morphological and physiological potentials of P. aeruginosa ATCC 27853.

An Efficient Approach to Explore and Discriminate Anomalous Regions in Bacterial Genomes Based on Maximum Entropy

Recently, there has been an increase in the number of whole bacterial genomes sequenced, mainly due to the advancing of next-generation sequencing technologies. In face of this, there is a need to provide new analytical alternatives that can follow this advance. Given our current knowledge about the genomic plasticity of bacteria and that those genomic regions can uncover important features about this microorganism, our goal was to develop a fast methodology based on maximum entropy (ME) to guide the researcher to regions that could be prioritized during the analysis. This methodology was compared with other available methods. In addition, ME was applied to eight different bacterial genera. The methodology consists of two main steps: processing the nucleotide sequence and ME calculation. We applied ME to Xanthomonas axonopodis pv. citri 306 (XAC) and Xanthomonas campestris pv. campestris ATCC 33913 (XCC), both of which have their anomalous regions well documented. We then compared our results against those from Alien Hunter, HGT-DB, Islander, IslandPath, and SIGI-HMM. ME was shown to be superior in terms of efficiency and analysis duration. Besides, ME only needs the genome sequence in FASTA format as input. The proposed strategy based on ME is able to help in bacterial genome exploration. This is a simple and fast strategy for individual genomes in comparison with other available methods, without relying on previous annotation and alignments. This methodology can also be a new option in the early stages of analysis of newly sequenced bacterial genomes.

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.

No evidence for a bovine mastitis Escherichia coli pathotype

Background

Escherichia coli bovine mastitis is a disease of significant economic importance in the dairy industry. Molecular characterization of mastitis-associated E. coli (MAEC) did not result in the identification of common traits. Nevertheless, a mammary pathogenic E. coli (MPEC) pathotype has been proposed suggesting virulence traits that differentiate MAEC from commensal E. coli. The present study was designed to investigate the MPEC pathotype hypothesis by comparing the genomes of MAEC and commensal bovine E. coli.

Results

We sequenced the genomes of eight E. coli isolated from bovine mastitis cases and six fecal commensal isolates from udder-healthy cows. We analyzed the phylogenetic history of bovine E. coli genomes by supplementing this strain panel with eleven bovine-associated E. coli from public databases. The majority of the isolates originate from phylogroups A and B1, but neither MAEC nor commensal strains could be unambiguously distinguished by phylogenetic lineage. The gene content of both MAEC and commensal strains is highly diverse and dominated by their phylogenetic background. Although individual strains carry some typical E. coli virulence-associated genes, no traits important for pathogenicity could be specifically attributed to MAEC. Instead, both commensal strains and MAEC have very few gene families enriched in either pathotype. Only the aerobactin siderophore gene cluster was enriched in commensal E. coli within our strain panel.

Conclusions

This is the first characterization of a phylogenetically diverse strain panel including several MAEC and commensal isolates. With our comparative genomics approach we could not confirm previous studies that argue for a positive selection of specific traits enabling MAEC to elicit bovine mastitis. Instead, MAEC are facultative and opportunistic pathogens recruited from the highly diverse bovine gastrointestinal microbiota. Virulence-associated genes implicated in mastitis are a by-product of commensalism with the primary function to enhance fitness in the bovine gastrointestinal tract. Therefore, we put the definition of the MPEC pathotype into question and suggest to designate corresponding isolates as MAEC.

Electronic supplementary material

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

Comparative Genomics Unravels the Functional Roles of Co-occurring Acidophilic Bacteria in Bioleaching Heaps

The spatial-temporal distribution of populations in various econiches is thought to be potentially related to individual differences in the utilization of nutrients or other resources, but their functional roles in the microbial communities remain elusive. We compared differentiation in gene repertoire and metabolic profiles, with a focus on the potential functional traits of three commonly recognized members (Acidithiobacillus caldus, Leptospirillum ferriphilum, and Sulfobacillus thermosulfidooxidans) in bioleaching heaps. Comparative genomics revealed that intra-species divergence might be driven by horizontal gene transfer. These co-occurring bacteria shared a few homologous genes, which significantly suggested the genomic differences between these organisms. Notably, relatively more genes assigned to the Clusters of Orthologous Groups category [G] (carbohydrate transport and metabolism) were identified in Sulfobacillus thermosulfidooxidans compared to the two other species, which probably indicated their mixotrophic capabilities that assimilate both organic and inorganic forms of carbon. Further inspection revealed distinctive metabolic capabilities involving carbon assimilation, nitrogen uptake, and iron-sulfur cycling, providing robust evidence for functional differences with respect to nutrient utilization. Therefore, we proposed that the mutual compensation of functionalities among these co-occurring organisms might provide a selective advantage for efficiently utilizing the limited resources in their habitats. Furthermore, it might be favorable to chemoautotrophs' lifestyles to form mutualistic interactions with these heterotrophic and/or mixotrophic acidophiles, whereby the latter could degrade organic compounds to effectively detoxify the environments. Collectively, the findings shed light on the genetic traits and potential metabolic activities of these organisms, and enable us to make some inferences about genomic and functional differences that might allow them to co-exist.

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.

Adaptive Evolution of Extreme Acidophile Sulfobacillus thermosulfidooxidans Potentially Driven by Horizontal Gene Transfer and Gene Loss

Recent phylogenomic analysis has suggested that three strains isolated from different copper mine tailings around the world were taxonomically affiliated with Sulfobacillusthermosulfidooxidans Here, we present a detailed investigation of their genomic features, particularly with respect to metabolic potentials and stress tolerance mechanisms. Comprehensive analysis of the Sulfobacillus genomes identified a core set of essential genes with specialized biological functions in the survival of acidophiles in their habitats, despite differences in their metabolic pathways. The Sulfobacillus strains also showed evidence for stress management, thereby enabling them to efficiently respond to harsh environments. Further analysis of metabolic profiles provided novel insights into the presence of genomic streamlining, highlighting the importance of gene loss as a main mechanism that potentially contributes to cellular economization. Another important evolutionary force, especially in larger genomes, is gene acquisition via horizontal gene transfer (HGT), which might play a crucial role in the recruitment of novel functionalities. Also, a successful integration of genes acquired from archaeal donors appears to be an effective way of enhancing the adaptive capacity to cope with environmental changes. Taken together, the findings of this study significantly expand the spectrum of HGT and genome reduction in shaping the evolutionary history of Sulfobacillus strains.IMPORTANCE Horizontal gene transfer (HGT) and gene loss are recognized as major driving forces that contribute to the adaptive evolution of microbial genomes, although their relative importance remains elusive. The findings of this study suggest that highly frequent gene turnovers within microorganisms via HGT were necessary to incur additional novel functionalities to increase the capacity of acidophiles to adapt to changing environments. Evidence also reveals a fascinating phenomenon of potential cross-kingdom HGT. Furthermore, genome streamlining may be a critical force in driving the evolution of microbial genomes. Taken together, this study provides insights into the importance of both HGT and gene loss in the evolution and diversification of bacterial genomes.

Genotypes Associated with Listeria monocytogenes Isolates Displaying Impaired or Enhanced Tolerances to Cold, Salt, Acid, or Desiccation Stress

The human pathogen Listeria monocytogenes is a large concern in the food industry where its continuous detection in food products has caused a string of recalls in North America and Europe. Most recognized for its ability to grow in foods during refrigerated storage, L. monocytogenes can also tolerate several other food-related stresses with some strains possessing higher levels of tolerances than others. The objective of this study was to use a combination of phenotypic analyses and whole genome sequencing to elucidate potential relationships between L. monocytogenes genotypes and food-related stress tolerance phenotypes. To accomplish this, 166 L. monocytogenes isolates were sequenced and evaluated for their ability to grow in cold (4°C), salt (6% NaCl, 25°C), and acid (pH 5, 25°C) stress conditions as well as survive desiccation (33% RH, 20°C). The results revealed that the stress tolerance of L. monocytogenes is associated with serotype, clonal complex (CC), full length inlA profiles, and the presence of a plasmid which was identified in 55% of isolates. Isolates with full length inlA exhibited significantly (p < 0.001) enhanced cold tolerance relative to those harboring a premature stop codon (PMSC) in this gene. Similarly, isolates possessing a plasmid demonstrated significantly (p = 0.013) enhanced acid tolerance. We also identified nine new L. monocytogenes sequence types, a new inlA PMSC, and several connections between CCs and the presence/absence or variations of specific genetic elements. A whole genome single-nucleotide-variants phylogeny revealed sporadic distribution of tolerant isolates and closely related sensitive and tolerant isolates, highlighting that minor genetic differences can influence the stress tolerance of L. monocytogenes. Specifically, a number of cold and desiccation sensitive isolates contained PMSCs in σB regulator genes (rsbS, rsbU, rsbV). Collectively, the results suggest that knowing the sequence type of an isolate in addition to screening for the presence of full-length inlA and a plasmid, could help food processors and food agency investigators determine why certain isolates might be persisting in a food processing environment. Additionally, increased sequencing of L. monocytogenes isolates in combination with stress tolerance profiling, will enhance the ability to identify genetic elements associated with higher risk strains.

Complete Genome Sequence of Paenibacillus polymyxa YC0573, a Plant Growth-Promoting Rhizobacterium with Antimicrobial Activity

Paenibacillus polymyxa strain YC0573 is a plant growth–promoting rhizobacterium with antimicrobial activity, which was isolated from tobacco rhizosphere. Here, we report the complete genome sequence of P. polymyxa YC0573. Antifungal and antibacterial genes were discovered.

Comparison of biofilm formation and motility processes in arsenic-resistant Thiomonas spp. strains revealed divergent response to arsenite

Bacteria of the genus Thiomonas are found ubiquitously in arsenic contaminated waters such as acid mine drainage (AMD), where they contribute to the precipitation and the natural bioremediation of arsenic. In these environments, these bacteria have developed a large range of resistance strategies among which the capacity to form particular biofilm structures. The biofilm formation is one of the most ubiquitous adaptive response observed in prokaryotes to various stresses, such as those induced in the presence of toxic compounds. This study focused on the process of biofilm formation in three Thiomonas strains (CB1, CB2 and CB3) isolated from the same AMD. The results obtained here show that these bacteria are all capable of forming biofilms, but the architecture and the kinetics of formation of these biofilms differ depending on whether arsenite is present in the environment and from one strain to another. Indeed, two strains favoured biofilm formation, whereas one favoured motility in the presence of arsenite. To identify the underlying mechanisms, the patterns of expression of some genes possibly involved in the process of biofilm formation were investigated in Thiomonas sp. CB2 in the presence and absence of arsenite, using a transcriptomic approach (RNA‐seq). The findings obtained here shed interesting light on how the formation of biofilms, and the motility processes contribute to the adaptation of Thiomonas strains to extreme environments.

Comparative Genomic Analysis of Two Vibrio toranzoniae Strains with Different Virulence Capacity Reveals Clues on Its Pathogenicity for Fish

Vibrio toranzoniae is a Gram-negative bacterium of the Splendidus clade within the Vibrio genus. V. toranzoniae was first isolated from healthy clams in Galicia (Spain) but recently was also identified associated to disease outbreaks of red conger eel in Chile. Experimental challenges showed that the Chilean isolates were able to produce fish mortalities but not the strains isolated from clams. The aim of the present study was to determine the differences at the genomic level between the type strain of the species (CECT 7225^T) and the strain R17, isolated from red conger eel in Chile, which could explain their different virulent capacity. The genome-based comparison showed high homology between both strains but differences were observed in certain gene clusters that include some virulence factors. Among these, we found that iron acquisition systems and capsule synthesis genes were the main differential features between both genomes that could explain the differences in the pathogenicity of the strains. Besides, the studied genomes presented genomic islands and toxins, and the R17 strain presented CRISPR sequences that are absent on the type strain. Taken together, this analysis provided important insights into virulence factors of V. toranzoniae that will lead to a better understanding of the pathogenic process.

Prevalence of Antibiotic Resistance Genes among Human Gut-Derived Bifidobacteria

The microbiota of the human gastrointestinal tract (GIT) may regularly be exposed to antibiotics, which are used to prevent and treat infectious diseases caused by bacteria and fungi. Bacterial communities of the gut retain a reservoir of antibiotic resistance (AR) genes, and antibiotic therapy thus positively selects for those microorganisms that harbor such genetic features, causing microbiota modulation. During the first months following birth, bifidobacteria represent some of the most dominant components of the human gut microbiota, although little is known about their AR gene complement (or resistome). In the current study, we assessed the resistome of the Bifidobacterium genus based on phenotypic and genotypic data of members that represent all currently recognized bifidobacterial (sub)species. Moreover, a comparison between the bifidobacterial resistome and gut metagenome data sets from adults and infants shows that the bifidobacterial community present at the first week following birth possesses a reduced AR arsenal compared to that present in the infant bifidobacterial population in subsequent weeks of the first year of life. Our findings reinforce the concept that the early infant gut microbiota is more susceptible to disturbances by antibiotic treatment than the gut microbiota developed at a later life stage.

IMPORTANCE

The spread of resistance to antibiotics among bacterial communities has represented a major concern since their discovery in the last century. The risk of genetic transfer of resistance genes between microorganisms has been extensively investigated due to its relevance to human health. In contrast, there is only limited information available on antibiotic resistance among human gut commensal microorganisms such as bifidobacteria, which are widely exploited by the food industry as health-promoting microorganisms or probiotic ingredients. In the current study, we explored the occurrence of antibiotic resistance genes in the genomes of bifidobacteria and evaluated their genetic mobility to other human gut commensal microorganisms.

Gene Turnover Contributes to the Evolutionary Adaptation of Acidithiobacillus caldus: Insights from Comparative Genomics

Acidithiobacillus caldus is an extremely acidophilic sulfur-oxidizer with specialized characteristics, such as tolerance to low pH and heavy metal resistance. To gain novel insights into its genetic complexity, we chosen six A. caldus strains for comparative survey. All strains analyzed in this study differ in geographic origins as well as in ecological preferences. Based on phylogenomic analysis, we clustered the six A. caldus strains isolated from various ecological niches into two groups: group 1 strains with smaller genomes and group 2 strains with larger genomes. We found no obvious intraspecific divergence with respect to predicted genes that are related to central metabolism and stress management strategies between these two groups. Although numerous highly homogeneous genes were observed, high genetic diversity was also detected. Preliminary inspection provided a first glimpse of the potential correlation between intraspecific diversity at the genome level and environmental variation, especially geochemical conditions. Evolutionary genetic analyses further showed evidence that the difference in environmental conditions might be a crucial factor to drive the divergent evolution of A. caldus species. We identified a diverse pool of mobile genetic elements including insertion sequences and genomic islands, which suggests a high frequency of genetic exchange in these harsh habitats. Comprehensive analysis revealed that gene gains and losses were both dominant evolutionary forces that directed the genomic diversification of A. caldus species. For instance, horizontal gene transfer and gene duplication events in group 2 strains might contribute to an increase in microbial DNA content and novel functions. Moreover, genomes undergo extensive changes in group 1 strains such as removal of potential non-functional DNA, which results in the formation of compact and streamlined genomes. Taken together, the findings presented herein show highly frequent gene turnover of A. caldus species that inhabit extremely acidic environments, and shed new light on the contribution of gene turnover to the evolutionary adaptation of acidophiles.

Comparative Genome Analysis of Fusobacterium nucleatum

Fusobacterium nucleatum is considered to be a key oral bacterium in recruiting periodontal pathogens into subgingival dental plaque. Currently F. nucleatum can be subdivided into five subspecies. Our previous genome analysis of F. nucleatum W1481 (referred to hereafter as W1481), isolated from an 8-mm periodontal pocket in a patient with chronic periodontitis, suggested the possibility of a new subspecies. To further investigate the biology and relationships of this possible subspecies with other known subspecies, we performed comparative analysis between W1481 and 35 genome sequences represented by the five known Fusobacterium subspecies. Our analyses suggest that W1481 is most likely a new F. nucleatum subspecies, supported by evidence from phylogenetic analyses and maximal unique match indices (MUMi). Interestingly, we found a horizontally transferred W1481-specific genomic island harboring the tripartite ATP-independent (TRAP)-like transporter genes, suggesting this bacterium might have a high-affinity transport system for the C4-dicarboxylates malate, succinate, and fumarate. Moreover, we found virulence genes in the W1481 genome that may provide a strong defense mechanism which might enable it to colonize and survive within the host by evading immune surveillance. This comparative study provides better understanding of F. nucleatum and the basis for future functional work on this important pathogen.

Genomic composition and dynamics among Methanomicrobiales predict adaptation to contrasting environments

Members of the order Methanomicrobiales are abundant, and sometimes dominant, hydrogenotrophic (H₂-CO₂ utilizing) methanoarchaea in a broad range of anoxic habitats. Despite their key roles in greenhouse gas emissions and waste conversion to methane, little is known about the physiological and genomic bases for their widespread distribution and abundance. In this study, we compared the genomes of nine diverse Methanomicrobiales strains, examined their pangenomes, reconstructed gene flow and identified genes putatively mediating their success across different habitats. Most strains slowly increased gene content whereas one, Methanocorpusculum labreanum, evidenced genome downsizing. Peat-dwelling Methanomicrobiales showed adaptations centered on improved transport of scarce inorganic nutrients and likely use H⁺ rather than Na⁺ transmembrane chemiosmotic gradients during energy conservation. In contrast, other Methanomicrobiales show the potential to concurrently use Na⁺ and H⁺ chemiosmotic gradients. Analyses also revealed that the Methanomicrobiales lack a canonical electron bifurcation system (MvhABGD) known to produce low potential electrons in other orders of hydrogenotrophic methanogens. Additional putative differences in anabolic metabolism suggest that the dynamics of interspecies electron transfer from Methanomicrobiales syntrophic partners can also differ considerably. Altogether, these findings suggest profound differences in electron trafficking in the Methanomicrobiales compared with other hydrogenotrophs, and warrant further functional evaluations.

Isolation and Characterization of NDM-Positive Escherichia coli from Municipal Wastewater in Jeddah, Saudi Arabia

The emergence of resistance to last-resort antibiotics is a public health concern of global scale. Besides direct person-to-person propagation, environmental pathways might contribute to the dissemination of antibiotic-resistant bacteria and antibiotic resistance genes (ARGs). Here, we describe the incidence of bla_NDM-1, a gene conferring resistance to carbapenems, in the wastewater of the city of Jeddah, Saudi Arabia, over a 1-year period. bla_NDM-1 was detected at concentrations ranging from 10⁴ to 10⁵ copies/m³ of untreated wastewater during the entire monitoring period. These results indicate the ubiquity and high incidence of bla_NDM-1 in the local wastewater. To track the bacteria carrying bla_NDM-1, we isolated Escherichia coli PI7, a strain of sequence type 101 (ST101), from wastewater around the Hajj event in October 2013. Genome sequencing of this strain revealed an extensive repertoire of ARGs as well as virulence and invasive traits. These traits were further confirmed by antibiotic resistance profiling and in vitro cell internalization in HeLa cell cultures. Given that this strain remains viable even after a certain duration in the sewerage, and that Jeddah lacks a robust sanitary infrastructure to fully capture all generated sewage, the presence of this bacterium in the untreated wastewater represents a potential hazard to the local public health. To the best of our knowledge, this is the first report of a bla_NDM-1-positive E. coli strain isolated from a nonnosocomial environment in Saudi Arabia and may set a priority concern for the need to establish improved surveillance for carbapenem-resistant E. coli in the country and nearby regions.

Identification of Novel Genomic Islands in Liverpool Epidemic Strain of Pseudomonas aeruginosa Using Segmentation and Clustering

Pseudomonas aeruginosa is an opportunistic pathogen implicated in a myriad of infections and a leading pathogen responsible for mortality in patients with cystic fibrosis (CF). Horizontal transfers of genes among the microorganisms living within CF patients have led to highly virulent and multi-drug resistant strains such as the Liverpool epidemic strain of P. aeruginosa, namely the LESB58 strain that has the propensity to acquire virulence and antibiotic resistance genes. Often these genes are acquired in large clusters, referred to as "genomic islands (GIs)." To decipher GIs and understand their contributions to the evolution of virulence and antibiotic resistance in P. aeruginosa LESB58, we utilized a recursive segmentation and clustering procedure, presented here as a genome-mining tool, "GEMINI." GEMINI was validated on experimentally verified islands in the LESB58 strain before examining its potential to decipher novel islands. Of the 6062 genes in P. aeruginosa LESB58, 596 genes were identified to be resident on 20 GIs of which 12 have not been previously reported. Comparative genomics provided evidence in support of our novel predictions. Furthermore, GEMINI unraveled the mosaic structure of islands that are composed of segments of likely different evolutionary origins, and demonstrated its ability to identify potential strain biomarkers. These newly found islands likely have contributed to the hyper-virulence and multidrug resistance of the Liverpool epidemic strain of P. aeruginosa.

Genomic analysis of three Bifidobacterium species isolated from the calf gastrointestinal tract

Ruminant animals contribute significantly to the global value of agriculture and rely on a complex microbial community for efficient digestion. However, little is known of how this microbial-host relationship develops and is maintained. To begin to address this, we have determined the ability of three Bifidobacterium species isolated from the faeces of newborn calves to grow on carbohydrates typical of a newborn ruminant diet. Genome sequences have been determined for these bacteria with analysis of the genomes providing insights into the host association and identification of several genes that may mediate interactions with the ruminant gastrointestinal tract. The present study provides a starting point from which we can define the role of potential beneficial microbes in the nutrition of young ruminants and begin to influence the interactions between the microbiota and the host. The differences observed in genomic content hint at niche partitioning among the bifidobacterial species analysed and the different strategies they employ to successfully adapt to this habitat.

Rhizoxin analogs, orfamide A and chitinase production contribute to the toxicity of Pseudomonas protegens strain Pf-5 to Drosophila melanogaster

Pseudomonas protegens strain Pf-5 is a soil bacterium that was first described for its capacity to suppress plant diseases and has since been shown to be lethal to certain insects. Among these is the common fruit fly Drosophila melanogaster, a well-established model organism for studies evaluating the molecular and cellular basis of the immune response to bacterial challenge. Pf-5 produces the insect toxin FitD, but a ΔfitD mutant of Pf-5 retained full toxicity against D. melanogaster in a noninvasive feeding assay, indicating that FitD is not a major determinant of Pf-5's oral toxicity against this insect. Pf-5 also produces a broad spectrum of exoenzymes and natural products with antibiotic activity, whereas a mutant with a deletion in the global regulatory gene gacA produces none of these exoproducts and also lacks toxicity to D. melanogaster. In this study, we made use of a panel of Pf-5 mutants having single or multiple mutations in the biosynthetic gene clusters for seven natural products and two exoenzymes that are produced by the bacterium under the control of gacA. Our results demonstrate that the production of rhizoxin analogs, orfamide A, and chitinase are required for full oral toxicity of Pf-5 against D. melanogaster, with rhizoxins being the primary determinant.

Characterization of the first cultured representative of Verrucomicrobia subdivision 5 indicates the proposal of a novel phylum

The recently isolated strain L21-Fru-AB^T represents moderately halophilic, obligately anaerobic and saccharolytic bacteria that thrive in the suboxic transition zones of hypersaline microbial mats. Phylogenetic analyses based on 16S rRNA genes, RpoB proteins and gene content indicated that strain L21-Fru-AB^T represents a novel species and genus affiliated with a distinct phylum-level lineage originally designated Verrucomicrobia subdivision 5. A survey of environmental 16S rRNA gene sequences revealed that members of this newly recognized phylum are wide-spread and ecologically important in various anoxic environments ranging from hypersaline sediments to wastewater and the intestine of animals. Characteristic phenotypic traits of the novel strain included the formation of extracellular polymeric substances, a Gram-negative cell wall containing peptidoglycan and the absence of odd-numbered cellular fatty acids. Unusual metabolic features deduced from analysis of the genome sequence were the production of sucrose as osmoprotectant, an atypical glycolytic pathway lacking pyruvate kinase and the synthesis of isoprenoids via mevalonate. On the basis of the analyses of phenotypic, genomic and environmental data, it is proposed that strain L21-Fru-AB^T and related bacteria are specifically adapted to the utilization of sulfated glycopolymers produced in microbial mats or biofilms.

Comparative genome-scale modelling of Staphylococcus aureus strains identifies strain-specific metabolic capabilities linked to pathogenicity

Staphylococcus aureus is a preeminent bacterial pathogen capable of colonizing diverse ecological niches within its human host. We describe here the pangenome of S. aureus based on analysis of genome sequences from 64 strains of S. aureus spanning a range of ecological niches, host types, and antibiotic resistance profiles. Based on this set, S. aureus is expected to have an open pangenome composed of 7,411 genes and a core genome composed of 1,441 genes. Metabolism was highly conserved in this core genome; however, differences were identified in amino acid and nucleotide biosynthesis pathways between the strains. Genome-scale models (GEMs) of metabolism were constructed for the 64 strains of S. aureus These GEMs enabled a systems approach to characterizing the core metabolic and panmetabolic capabilities of the S. aureus species. All models were predicted to be auxotrophic for the vitamins niacin (vitamin B3) and thiamin (vitamin B1), whereas strain-specific auxotrophies were predicted for riboflavin (vitamin B2), guanosine, leucine, methionine, and cysteine, among others. GEMs were used to systematically analyze growth capabilities in more than 300 different growth-supporting environments. The results identified metabolic capabilities linked to pathogenic traits and virulence acquisitions. Such traits can be used to differentiate strains responsible for mild vs. severe infections and preference for hosts (e.g., animals vs. humans). Genome-scale analysis of multiple strains of a species can thus be used to identify metabolic determinants of virulence and increase our understanding of why certain strains of this deadly pathogen have spread rapidly throughout the world.