PIPS: pathogenicity island prediction software

Pan-Genome Plasticity and Virulence Factors: A Natural Treasure Trove for Acinetobacter baumannii

Acinetobacter baumannii is a Gram-negative pathogen responsible for a variety of community- and hospital-acquired infections. It is recognized as a life-threatening pathogen among hospitalized individuals and, in particular, immunocompromised patients in many countries. A. baumannii, as a member of the ESKAPE group, encompasses high genomic plasticity and simultaneously is predisposed to receive and exchange the mobile genetic elements (MGEs) through horizontal genetic transfer (HGT). Indeed, A. baumannii is a treasure trove that contains a high number of virulence factors. In accordance with these unique pathogenic characteristics of A. baumannii, the authors aim to discuss the natural treasure trove of pan-genome and virulence factors pertaining to this bacterial monster and try to highlight the reasons why this bacterium is a great concern in the global public health system.

Genomic islands and their role in fitness traits of two key sepsis-causing bacterial pathogens

To survive and establish a niche for themselves, bacteria constantly evolve. Toward that, they not only insert point mutations and promote illegitimate recombinations within their genomes but also insert pieces of 'foreign' deoxyribonucleic acid, which are commonly referred to as 'genomic islands' (GEIs). The GEIs come in several forms, structures and types, often providing a fitness advantage to the harboring bacterium. In pathogenic bacteria, some GEIs may enhance virulence, thus altering disease burden, morbidity and mortality. Hence, delineating (i) the GEIs framework, (ii) their encoded functions, (iii) the triggers that help them move, (iv) the mechanisms they exploit to move among bacteria and (v) identification of their natural reservoirs will aid in superior tackling of several bacterial diseases, including sepsis. Given the vast array of comparative genomics data, in this short review, we provide an overview of the GEIs, their types and the compositions therein, especially highlighting GEIs harbored by two important pathogens, viz. Acinetobacter baumannii and Klebsiella pneumoniae, which prominently trigger sepsis in low- and middle-income countries. Our efforts help shed some light on the challenges these pathogens pose when equipped with GEIs. We hope that this review will provoke intense research into understanding GEIs, the cues that drive their mobility across bacteria and the ways and means to prevent their transfer, especially across pathogenic bacteria.

Horizontal transfer of probable chicken-pathogenicity chromosomal islands between Staphylococcus aureus and Staphylococcus agnetis

Staphylococcus agnetis is an emerging pathogen in chickens but has been most commonly isolated from sub-clinical mastitis in bovines. Previous whole-genome analyses for known virulence genes failed to identify determinants for the switch from mild ductal infections in cattle to severe infections in poultry. We now report identification of a family of 15 kbp, 17-19 gene mobile genetic elements (MGEs) specific to chicken osteomyelitis and dermatitis isolates of S. agnetis. These MGEs can be present in multiple copies per genome. The MGE has been vectored on a Staphylococcus phage that separately lysogenized two S. agnetis osteomyelitis strains. The S. agnetis genome from a broiler breeder case of ulcerative dermatitis contains 2 orthologs of this MGE, not associated with a prophage. BLASTn and phylogenetic analyses show that there are closely related intact MGEs found in genomes of S. aureus. The genome from a 1980s isolate from chickens in Ireland contains 3 copies of this MGE. More recent chicken isolates descended from that genome (Poland 2009, Oklahoma 2010, and Arkansas 2018) contain 2 to 4 related copies. Many of the genes of this MGE can be identified in disparate regions of the genomes of other chicken isolates of S. aureus. BLAST searches of the NCBI databases detect no similar MGEs outside of S. aureus and S. agnetis. These MGEs encode no proteins related to those produced by Staphylococcus aureus Pathogenicity Islands, which have been associated with the transition of S. aureus from human to chicken hosts. Other than mobilization functions, most of the genes in these new MGEs annotate as hypothetical proteins. The MGEs we describe appear to represent a new family of Chromosomal Islands (CIs) shared amongst S. agnetis and S. aureus. Further work is needed to understand the role of these CIs/MGEs in pathogenesis. Analysis of horizontal transfer of genetic elements between isolates and species of Staphylococci provides clues to evolution of host-pathogen interactions as well as revealing critical determinants for animal welfare and human diseases.

Comparative genomic analysis of Canadian Mycoplasma bovis strains isolated from Bison and Cattle

Mycoplasma bovis (M. bovis) in cattle causes pneumonia, arthritis, otitis media, and mastitis. In addition, multiple outbreaks have been recorded in North American bison. The genomic data on Canadian M. bovis in bison and cattle to date is limited. Whole-genome sequencing (WGS) was used to assess the degree of genome conservation across four Canadian M. bovis strains recovered from bison and cattle. Whole-genome sequences of four M. bovis isolates (Mb1, Mb160, Mb300, Mb304) and the PG45 reference genome were utilized to identify the M. bovis genomic similarity, whole-genome single nucleotide polymorphism (WGS-SNP), virulence determinants, and genomic islands. The pan-genome analysis showed that M. bovis encodes a minimum of 971 genes, while the core genome contained 637 genes. Comparative genomics revealed limited diversity in gene content between bison and cattle isolates. Whole-genome SNP analysis showed that the four M. bovis isolates differed from each other and to PG45. A total of 40 putative virulence genes associated with adhesion, colonization, and destruction of tissues were found in the bison and cattle isolates using the virulence factors database (VFDB). These putative virulence factors were equally distributed among isolates. Genomic Islands (GIs) ranging from 4 to 9 and associated with transposases, restriction-modification, ribosomal hypothetical proteins, variable surface lipoproteins, and unknowns were also identified. Overall, the genomic characterization of these isolates may provide new insights into the mechanisms of pathogenicity in M. bovis.

Performance assessment of genomic island prediction tools with an improved version of Design-Island

Genomic Islands (GIs) play an important role in the evolution and adaptation of prokaryotes. The origin and extent of ecological diversity of prokaryotes can be analyzed by comparing GIs across closely or distantly related prokaryotes. Understanding the importance of GI and to study the bacterial evolution, several GI prediction tools have been generated. An unsupervised method, Design-Island, was developed to identify GIs using Monte-Carlo statistical test on randomly selected segments of a chromosome. Here, in the present study Design-Island was modified with the incorporation of majority voting, multiple hypothesis testing correction. The performance of the modified version, Design-Island-II was tested and compared with the existing GI prediction tools. The performance assessment and benchmarking of the GI prediction tools require experimentally validated dataset, which is lacking. So, different datasets, generated or taken from literature were utilized to compare the sensitivity (SN), specificity (SP), precision (PPV) and accuracy (AC) of Design-Island-II. It showed substantial enhancement in term of SN, SP, PPV and AC, and significantly reduced the computation time of the algorithm. The performance of Design-Island-II has also been compared with several GI prediction tools using curated dataset of putative horizontally transferred genes. Design-Island-II showed the highest sensitivity and F1 score, comparable specificity, precision and accuracy in comparison to the other available methods. IslandViewer4 and Islander outperformed all the available methods in terms of AC and PPV respectively. Our study suggested Design-Island-II, IslandViewer4 and GIHunter among the top performing GI prediction tools considering both sensitivity and specificity of the methods.

In silico Prediction of New Drug Candidates Against the Multidrug-Resistant and Potentially Zoonotic Fish Pathogen Serotype III Streptococcus agalactiae

Streptococcus agalactiae is an invasive multi-host pathogen that causes invasive diseases mainly in newborns, elderly, and individuals with underlying health complications. In fish, S. agalactiae causes streptococcosis, which is characterized by septicemia and neurological signs, and leads to great economic losses to the fish farming industry worldwide. These bacteria can be classified into different serotypes based on capsular antigens, and into different sequence types (ST) based on multilocus sequence typing (MLST). In 2015, serotype III ST283 was identified to be associated with a foodborne invasive disease in non-pregnant immunocompetent humans in Singapore, and the infection was related to raw fish consumption. In addition, a serotype III strain isolated from tilapia in Brazil has been reported to be resistant to five antibiotic classes. This specific serotype can serve as a reservoir of resistance genes and pose a serious threat to public health. Thus, new approaches for the control and treatment of S. agalactiae infections are needed. In the present study, 24 S. agalactiae serotype III complete genomes, isolated from human and fish hosts, were compared. The core genome was identified, and, using bioinformatics tools and subtractive criteria, five proteins were identified as potential drug targets. Furthermore, 5,008 drug-like natural compounds were virtually screened against the identified targets. The ligands with the best binding properties are suggested for further in vitro and in vivo analysis.

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.

A panel of recombinant proteins for the serodiagnosis of caseous lymphadenitis in goats and sheep

Caseous lymphadenitis (CLA) is a small ruminant disease characterized by the development of granulomatous lesions in superficial and internal lymph nodes, as well as in some organs, and causes significant economic losses worldwide. The aetiological agent of CLA is the bacterium Corynebacterium pseudotuberculosis; however, the commercially available diagnostic tools present problems with regard to specificity, which can lead to false-negative results. This study aimed to develop an indirect enzyme-linked immunosorbent assay (ELISA) for the detection of specific immunoglobulins in goats and sheep using recombinant C. pseudotuberculosis PLD, CP40, PknG, DtxR and Grx proteins. For validation of the ELISAs, 130 goat serum samples and 160 sheep serum samples were used. The best ELISA for goats was developed using a combination of PLD and CP40 as antigens at a 1:1 ratio, which presented 96.9% sensitivity and 98.4% specificity. The most effective ELISA for sheep presented 91% sensitivity and 98.7% specificity when recombinant PLD alone was used as the antigen. These ELISAs can be used as highly accurate tools in epidemiological surveys and for the serodiagnosis of C. pseudotuberculosis infection in goats and sheep.

Reverse vaccinology and subtractive genomics reveal new therapeutic targets against Mycoplasma pneumoniae: a causative agent of pneumonia

Pneumonia is an infectious disease caused by bacteria, viruses or fungi that results in millions of deaths globally. Despite the existence of prophylactic methods against some of the major pathogens of the disease, there is no efficient prophylaxis against atypical agents such as Mycoplasma pneumoniae, a bacterium associated with cases of community-acquired pneumonia. Because of the morphological peculiarity of M. pneumoniae, which leads to an increased resistance to antibiotics, studies that prospectively investigate the development of vaccines and drug targets appear to be one of the best ways forward. Hence, in this paper, bioinformatics tools were used for vaccine and pharmacological prediction. We conducted comparative genomic analysis on the genomes of 88 M. pneumoniae strains, as opposed to a reverse vaccinology analysis, in relation to the capacity of M. pneumoniae proteins to bind to the major histocompatibility complex, revealing seven targets with immunogenic potential. Predictive cytoplasmic proteins were tested as potential drug targets by studying their structures in relation to other proteins, metabolic pathways and molecular anchorage, which identified five possible drug targets. These findings are a valuable addition to the development of vaccines and the selection of new in vivo drug targets that may contribute to further elucidating the molecular basis of M. pneumoniae-host interactions.

Comparative Analysis of Genomic Island Prediction Tools

Tools for genomic island prediction use strategies for genomic comparison analysis and sequence composition analysis. The goal of comparative analysis is to identify unique regions in the genomes of related organisms, whereas sequence composition analysis evaluates and relates the composition of specific regions with other regions in the genome. The goal of this study was to qualitatively and quantitatively evaluate extant genomic island predictors. We chose tools reported to produce significant results using sequence composition prediction, comparative genomics, and hybrid genomics methods. To maintain diversity, the tools were applied to eight complete genomes of organisms with distinct characteristics and belonging to different families. Escherichia coli CFT073 was used as a control and considered as the gold standard because its islands were previously curated in vitro. The results of predictions with the gold standard were manually curated, and the content and characteristics of each predicted island were analyzed. For other organisms, we created GenBank (GBK) files using Artemis software for each predicted island. We copied only the amino acid sequences from the coding sequence and constructed a multi-FASTA file for each predictor. We used BLASTp to compare all results and generate hits to evaluate similarities and differences among the predictions. Comparison of the results with the gold standard revealed that GIPSy produced the best results, covering ~91% of the composition and regions of the islands, followed by Alien Hunter (81%), IslandViewer (47.8%), Predict Bias (31%), GI Hunter (17%), and Zisland Explorer (16%). The tools with the best results in the analyzes of the set of organisms were the same ones that presented better performance in the tests with the gold standard.

Rapidly evolving changes and gene loss associated with host switching in Corynebacterium pseudotuberculosis

Phylogenomics and genome scale positive selection analyses were performed on 29 Corynebacterium pseudotuberculosis genomes that were isolated from different hosts, including representatives of the Ovis and Equi biovars. A total of 27 genes were identified as undergoing adaptive changes. An analysis of the clades within this species and these biovars, the genes specific to each branch, and the genes responding to selective pressure show clear differences, indicating that adaptation and specialization is occurring in different clades. These changes are often correlated with the isolation host but could indicate responses to some undetermined factor in the respective niches. The fact that some of these more-rapidly evolving genes have homology to known virulence factors, antimicrobial resistance genes and drug targets shows that this type of analysis could be used to identify novel targets, and that these could be used as a way to control this pathogen.

Transcriptome analysis of Corynebacterium pseudotuberculosis biovar Equi in two conditions of the environmental stress

Corynebacterium pseudotuberculosis has been widely studied in an effort to understand its biological evolution. Transcriptomics has revealed possible candidates for virulence and pathogenicity factors of strain 1002 (biovar Ovis). Because C. pseudotuberculosis is classified into two biovars, Ovis and Equi, it was interesting to assess the transcriptional profile of biovar Equi strain 258, the causative agent of ulcerative lymphangitis. The genome of this strain was re-sequenced; the reassembly was completed using optical mapping technology, and the sequence was subsequently re-annotated. Two growth conditions that occur during the host infection process were simulated for the transcriptome: the osmotic and acid medium. Genes that may be associated with the microorganism's resilience under unfavorable conditions were identified through RNAseq, including genes present in pathogenicity islands. The RT-qPCR was performed to confirm the results in biological triplicate for each condition for some genes. The results extend our knowledge of the factors associated with the spread and persistence of C. pseudotuberculosis during the infection process and suggest possible avenues for studies related to the development of vaccines, diagnosis, and therapies that might help minimize damage to agribusinesses.

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.

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.

Heat shock stress: Profile of differential expression in Corynebacterium pseudotuberculosis biovar Equi

Transcriptome studies on Corynebacterium pseudotuberculosis have recently contributed to the understanding about this microorganism's survival mechanisms in various hostile conditions. The gene expression profile of the C. pseudotuberculosis strain 1002 (Ovis biovar), has revealed genes that are possible candidates responsible for its maintenance in adverse environments, such as those found in the host. In another strain of this bacterium, 258 (Equi biovar), a high temperature condition was simulated, in order to verify which genes are responsible for promoting the persistence of the bacterium in these conditions, since it tolerates temperatures higher than 40°C, despite being a mesophilic bacterium. It was possible to generate a list of genes using RNAseq technology that possibly contribute to the survival of the bacteria in this hostile environment. A total of 562 genes were considered as differentially expressed, then, after the fold-change cutoff, 113 were considered induced and 114 repressed, resulting in a total of 227 genes. Therefore, hypothetical proteins presented a fold change above 6, and genes characteristically in control for this type of stress, such as hspR, grpE, and dnaK, presented a fold change above 3. The clpB gene, a chaperone, drew attention due to presenting a fold change above 3 and located in a pathogenicity island. These genes may contribute towards efficient solutions to the effects caused by ulcerative lymphangitis in equines, thus attenuating the damage it causes to agribusiness.

Identification of genomic islands in Chilean Piscirickettsia salmonis strains and analysis of gene expression involved in virulence

Piscirickettsia salmonis, an agent of Piscirickettsiosis, is the cause of major losses in the Chilean salmon industry. We identified, characterized and bioinformatically analysed genomic islands in field strains of P. Salmonis, using the bioinformatic software PIPS, that uses the characteristics of the islands of pathogenicity to identify them. We analysed nine partially sequenced genomes in different new field strains, and compared them with the LF-89 (Type strain) genome, selecting a genomic island present in all of them. We then evaluated the relative expression of three genes present in that island. From the obtained results, we conclude that the expression of the tcf gene is directly proportional to the cytopathogenicity in vitro of the bacteria; the product of the dnsa gene could contribute to its pathogenicity, but would be potentiated by one or more factors. The product of the gene liso is necessary for the virulence process and could have functions in early stages of infection. Regarding the strains, the IBM-040 strain showed a significant increase in the expression of all the genes in the study. Contrarily, LF-89 only presented a significant increase in expression of the gene liso, which correlates with the cytopathogenicity in vitro observed in the SHK-1 cells.

Quantitative Proteomic Analysis Reveals Changes in the Benchmark Corynebacterium pseudotuberculosis Biovar Equi Exoproteome after Passage in a Murine Host

Corynebacterium pseudotuberculosis biovar equi is the etiologic agent of ulcerative lymphangitis. To investigate proteins that could be related to the virulence of this pathogen, we combined an experimental passage process using a murine model and high-throughput proteomics with a mass spectrometry, data-independent acquisition (LC-MS^E) approach to identify and quantify the proteins released into the supernatants of strain 258_equi. To our knowledge, this approach allowed characterization of the exoproteome of a C. pseudotuberculosis equi strain for the first time. Interestingly, the recovery of this strain from infected mouse spleens induced a change in its virulence potential, and it became more virulent in a second infection challenge. Proteomic screening performed from culture supernatant of the control and recovered conditions revealed 104 proteins that were differentially expressed between the two conditions. In this context, proteomic analysis of the recovered condition detected the induction of proteins involved in bacterial pathogenesis, mainly related to iron uptake. In addition, KEGG enrichment analysis showed that ABC transporters, bacterial secretion systems and protein export pathways were significantly altered in the recovered condition. These findings show that secretion and secreted proteins are key elements in the virulence and adaptation of C. pseudotuberculosis. Collectively, bacterial pathogenesis-related proteins were identified that contribute to the processes of adherence, intracellular growth and evasion of the immune system. Moreover, this study enhances our understanding of the factors that may influence the pathogenesis of C. pseudotuberculosis.

LifeStyle-Specific-Islands (LiSSI): Integrated Bioinformatics Platform for Genomic Island Analysis

Distinct bacteria are able to cope with highly diverse lifestyles; for instance, they can be free living or host-associated. Thus, these organisms must possess a large and varied genomic arsenal to withstand different environmental conditions. To facilitate the identification of genomic features that might influence bacterial adaptation to a specific niche, we introduce LifeStyle-Specific-Islands (LiSSI). LiSSI combines evolutionary sequence analysis with statistical learning (Random Forest with feature selection, model tuning and robustness analysis). In summary, our strategy aims to identify conserved consecutive homology sequences (islands) in genomes and to identify the most discriminant islands for each lifestyle.

Intriguing Interaction of Bacteriophage-Host Association: An Understanding in the Era of Omics

Innovations in next-generation sequencing technology have introduced new avenues in microbial studies through “omics” approaches. This technology has considerably augmented the knowledge of the microbial world without isolation prior to their identification. With an enormous volume of bacterial “omics” data, considerable attempts have been recently invested to improve an insight into virosphere. The interplay between bacteriophages and their host has created a significant influence on the biogeochemical cycles, microbial diversity, and bacterial population regulation. This review highlights various concepts such as genomics, transcriptomics, proteomics, and metabolomics to infer the phylogenetic affiliation and function of bacteriophages and their impact on diverse microbial communities. Omics technologies illuminate the role of bacteriophage in an environment, the influences of phage proteins on the bacterial host and provide information about the genes important for interaction with bacteria. These investigations will reveal some of bio-molecules and biomarkers of the novel phage which demand to be unveiled.

A shift in the virulence potential of Corynebacterium pseudotuberculosis biovar ovis after passage in a murine host demonstrated through comparative proteomics

Background

Corynebacterium pseudotuberculosis biovar ovis, a facultative intracellular pathogen, is the etiologic agent of caseous lymphadenitis in small ruminants. During the infection process, C. pseudotuberculosis changes its gene expression to resist different types of stresses and to evade the immune system of the host. However, factors contributing to the infectious process of this pathogen are still poorly documented. To better understand the C. pseudotuberculosis infection process and to identify potential factors which could be involved in its virulence, experimental infection was carried out in a murine model using the strain 1002_ovis and followed by a comparative proteomic analysis of the strain before and after passage.

Results

The experimental infection assays revealed that strain 1002_ovis exhibits low virulence potential. However, the strain recovered from the spleen of infected mice and used in a new infection challenge showed a dramatic change in its virulence potential. Label-free proteomic analysis of the culture supernatants of strain 1002_ovis before and after passage in mice revealed that 118 proteins were differentially expressed. The proteome exclusive to the recovered strain contained important virulence factors such as CP40 proteinase and phospholipase D exotoxin, the major virulence factor of C. pseudotuberculosis. Also, the proteome from recovered condition revealed different classes of proteins involved in detoxification processes, pathogenesis and export pathways, indicating the presence of distinct mechanisms that could contribute in the infectious process of this pathogen.

Conclusions

This study shows that C. pseudotuberculosis modifies its proteomic profile in the laboratory versus infection conditions and adapts to the host context during the infection process. The screening proteomic performed us enable identify known virulence factors, as well as potential proteins that could be related to virulence this pathogen. These results enhance our understanding of the factors that might influence in the virulence of C. pseudotuberculosis.

Electronic supplementary material

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

Intraclonal Genome Stability of the Metallo-β-lactamase SPM-1-producing Pseudomonas aeruginosa ST277, an Endemic Clone Disseminated in Brazilian Hospitals

Carbapenems represent the mainstay therapy for the treatment of serious P. aeruginosa infections. However, the emergence of carbapenem resistance has jeopardized the clinical use of this important class of compounds. The production of SPM-1 metallo-β-lactamase has been the most common mechanism of carbapenem resistance identified in P. aeruginosa isolated from Brazilian medical centers. Interestingly, a single SPM-1-producing P. aeruginosa clone belonging to the ST277 has been widely spread within the Brazilian territory. In the current study, we performed a next-generation sequencing of six SPM-1-producing P. aeruginosa ST277 isolates. The core genome contains 5899 coding genes relative to the reference strain P. aeruginosa PAO1. A total of 26 genomic islands were detected in these isolates. We identified remarkable elements inside these genomic islands, such as copies of the bla_SPM−1 gene conferring resistance to carbapenems and a type I-C CRISPR-Cas system, which is involved in protection of the chromosome against foreign DNA. In addition, we identified single nucleotide polymorphisms causing amino acid changes in antimicrobial resistance and virulence-related genes. Together, these factors could contribute to the marked resistance and persistence of the SPM-1-producing P. aeruginosa ST277 clone. A comparison of the SPM-1-producing P. aeruginosa ST277 genomes showed that their core genome has a high level nucleotide similarity and synteny conservation. The variability observed was mainly due to acquisition of genomic islands carrying several antibiotic resistance genes.

Deciphering the virulence factors of the opportunistic pathogen Mycobacterium colombiense

Mycobacterium avium complex (MAC) contains clinically important nontuberculous mycobacteria worldwide and is the second largest medical complex in the Mycobacterium genus after the Mycobacterium tuberculosis complex. MAC comprises several species that are closely phylogenetically related but diverse regarding their host preference, course of disease, virulence and immune response. In this study we provided immunologic and virulence-related insights into the M. colombiense genome as a model of an opportunistic pathogen in the MAC. By using bioinformatic tools we found that M. colombiense has deletions in the genes involved in p-HBA/PDIM/PGL, PLC, SL-1 and HspX production, and loss of the ESX-1 locus. This information not only sheds light on our understanding the virulence mechanisms used by opportunistic MAC pathogens but also has great potential for the designing of species-specific diagnostic tools.

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.

Comparative Genome Analysis Provides Insights into the Pathogenicity of Flavobacterium psychrophilum

Flavobacterium psychrophilum is a fish pathogen in salmonid aquaculture worldwide that causes cold water disease (CWD) and rainbow trout fry syndrome (RTFS). Comparative genome analyses of 11 F. psychrophilum isolates representing temporally and geographically distant populations were used to describe the F. psychrophilum pan-genome and to examine virulence factors, prophages, CRISPR arrays, and genomic islands present in the genomes. Analysis of the genomic DNA sequences were complemented with selected phenotypic characteristics of the strains. The pan genome analysis showed that F. psychrophilum could hold at least 3373 genes, while the core genome contained 1743 genes. On average, 67 new genes were detected for every new genome added to the analysis, indicating that F. psychrophilum possesses an open pan genome. The putative virulence factors were equally distributed among isolates, independent of geographic location, year of isolation and source of isolates. Only one prophage-related sequence was found which corresponded to the previously described prophage 6H, and appeared in 5 out of 11 isolates. CRISPR array analysis revealed two different loci with dissimilar spacer content, which only matched one sequence in the database, the temperate bacteriophage 6H. Genomic Islands (GIs) were identified in F. psychrophilum isolates 950106-1/1 and CSF 259–93, associated with toxins and antibiotic resistance. Finally, phenotypic characterization revealed a high degree of similarity among the strains with respect to biofilm formation and secretion of extracellular enzymes. Global scale dispersion of virulence factors in the genomes and the abilities for biofilm formation, hemolytic activity and secretion of extracellular enzymes among the strains suggested that F. psychrophilum isolates have a similar mode of action on adhesion, colonization and destruction of fish tissues across large spatial and temporal scales of occurrence. Overall, the genomic characterization and phenotypic properties may provide new insights to the mechanisms of pathogenicity in F. psychrophilum.

Whole-Genome Sequence of Corynebacterium pseudotuberculosis 262 Biovar equi Isolated from Cow Milk

We report the complete genome sequence ofCorynebacterium pseudotuberculosis262, isolated from a bovine host.C. pseudotuberculosisis an etiological agent of diseases with medical and veterinary relevance. The genome contains 2,325,749 bp, 52.8% G+C content, 2,022 coding sequences (CDS), 50 pseudogenes, 48 tRNAs, and 12 rRNAs.

Inside the Pan-genome - Methods and Software Overview

The number of genomes that have been deposited in databases has increased exponentially after the advent of Next-Generation Sequencing (NGS), which produces high-throughput sequence data; this circumstance has demanded the development of new bioinformatics software and the creation of new areas, such as comparative genomics. In comparative genomics, the genetic content of an organism is compared against other organisms, which helps in the prediction of gene function and coding region sequences, identification of evolutionary events and determination of phylogenetic relationships. However, expanding comparative genomics to a large number of related bacteria, we can infer their lifestyles, gene repertoires and minimal genome size. In this context, a powerful approach called Pan-genome has been initiated and developed. This approach involves the genomic comparison of different strains of the same species, or even genus. Its main goal is to establish the total number of non-redundant genes that are present in a determined dataset. Pan-genome consists of three parts: core genome; accessory or dispensable genome; and species-specific or strain-specific genes. Furthermore, pan-genome is considered to be “open” as long as new genes are added significantly to the total repertoire for each new additional genome and “closed” when the newly added genomes cannot be inferred to significantly increase the total repertoire of the genes. To perform all of the required calculations, a substantial amount of software has been developed, based on orthologous and paralogous gene identification.

Whole-Genome Sequence of Corynebacterium pseudotuberculosis Strain 226, Isolated from the Abscess of a Goat in California

Corynebacterium pseudotuberculosis is the etiological agent of a caseous lymphadenitis disease. Herein, we present the first complete genome sequencing of C. pseudotuberculosis strain 226, isolated from an abscess of the sub-iliac lymph node of a goat from California (USA). The genome contains 2,138 coding sequences (CDSs), 12 rRNAs, 49 tRNAs, and 72 pseudogenes.

Comparative genome analysis of Weissella ceti, an emerging pathogen of farm-raised rainbow trout

Background

The genus Weissella belongs to the lactic acid bacteria and includes 18 currently identified species, predominantly isolated from fermented food but rarely from cases of bacteremia in animals. Recently, a new species, designated Weissella ceti, has been correlated with hemorrhagic illness in farm-raised rainbow trout in China, Brazil, and the USA, with high transmission and mortality rates during outbreaks. Although W. ceti is an important emerging veterinary pathogen, little is known about its genomic features or virulence mechanisms. To better understand these and to characterize the species, we have previously sequenced the genomes of W. ceti strains WS08, WS74, and WS105, isolated from different rainbow trout farms in Brazil and displaying different pulsed-field gel electrophoresis patterns. Here, we present a comparative analysis of the three previously sequenced genomes of W. ceti strains from Brazil along with W. ceti NC36 from the USA and those of other Weissella species.

Results

Phylogenomic and orthology-based analyses both showed a high-similarity in the genetic structure of these W. ceti strains. This structure is corroborated by the highly syntenic order of their genes and the neutral evolution inferred from Tajima’s D. A whole-genome multilocus sequence typing analysis distinguished strains WS08 and NC36 from strains WS74 and WS105. We predicted 10 putative genomic islands (GEI), among which PAIs 3a and 3b are phage sequences that occur only in WS105 and WS74, respectively, whereas PAI 1 is species specific.

Conclusions

We identified several genes putatively involved in the basic processes of bacterial physiology and pathogenesis, including survival in aquatic environment, adherence in the host, spread inside the host, resistance to immune-system-mediated stresses, and antibiotic resistance. These data provide new insights in the molecular epidemiology and host adaptation for this emerging pathogen in aquaculture.

Electronic supplementary material

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

Genome informatics and vaccine targets in Corynebacterium urealyticum using two whole genomes, comparative genomics, and reverse vaccinology

BACKGROUND

Corynebacterium urealyticum is an opportunistic pathogen that normally lives on skin and mucous membranes in humans. This high Gram-positive bacteria can cause acute or encrusted cystitis, encrusted pyelitis, and pyelonephritis in immunocompromised patients. The bacteria is multi-drug resistant, and knowledge about the genes that contribute to its virulence is very limited. Two complete genome sequences were used in this comparative genomic study: C. urealyticum DSM 7109 and C. urealyticum DSM 7111.

RESULTS

We used comparative genomics strategies to compare the two strains, DSM 7109 and DSM 7111, and to analyze their metabolic pathways, genome plasticity, and to predict putative antigenic targets. The genomes of these two strains together encode 2,115 non-redundant coding sequences, 1,823 of which are common to both genomes. We identified 188 strain-specific genes in DSM 7109 and 104 strain-specific genes in DSM 7111. The high number of strain-specific genes may be a result of horizontal gene transfer triggered by the large number of transposons in the genomes of these two strains. Screening for virulence factors revealed the presence of the spaDEF operon that encodes pili forming proteins. Therefore, spaDEF may play a pivotal role in facilitating the adhesion of the pathogen to the host tissue. Application of the reverse vaccinology method revealed 19 putative antigenic proteins that may be used in future studies as candidate drug or vaccine targets.

CONCLUSIONS

The genome features and the presence of virulence factors in genomic islands in the two strains of C. urealyticum provide insights in the lifestyle of this opportunistic pathogen and may be useful in developing future therapeutic strategies.

An integrated structural proteomics approach along the druggable genome of Corynebacterium pseudotuberculosis species for putative druggable targets

Background

The bacterium Corynebacterium pseudotuberculosis (Cp) causes caseous lymphadenitis (CLA), mastitis, ulcerative lymphangitis, and oedema in a number of hosts, comprising ruminants, thereby intimidating economic and dairy industries worldwide. So far there is no effective drug or vaccine available against Cp. Previously, a pan-genomic analysis was performed for both biovar equi and biovar ovis and a Pathogenicity Islands (PAIS) analysis within the strains highlighted a large set of proteins that could be relevant therapeutic targets for controlling the onset of CLA. In the present work, a structural druggability analysis pipeline was accomplished along 15 previously sequenced Cp strains from both biovar equi and biovar ovis.

Methods and results

We computed the whole modelome of a reference strain Cp1002 (NCBI Accession: NC_017300.1) and then the homology models of proteins, of 14 different Cp strains, with high identity (≥ 85%) to the reference strain were also done. Druggability score of all proteins pockets was calculated and only those targets that have a highly druggable (HD) pocket in all strains were kept, a set of 58 proteins. Finally, this information was merged with the previous PAIS analysis giving two possible highly relevant targets to conduct drug discovery projects. Also, off-targeting information against host organisms, including Homo sapiens and a further analysis for protein essentiality provided a final set of 31 druggable, essential and non-host homologous targets, tabulated in table S4, additional file 1. Out of 31 globally druggable targets, 9 targets have already been reported in other pathogenic microorganisms, 3 of them (3-isopropylmalate dehydratase small subunit, 50S ribosomal protein L30, Chromosomal replication initiator protein DnaA) in C. pseudotuberculosis.

Conclusion

Overall we provide valuable information of possible targets against C. pseudotuberculosis where some of these targets have already been reported in other microorganisms for drug discovery projects, also discarding targets that might be physiologically relevant but are not amenable for drug binding. We propose that the constructed in silico dataset might serve as a guidance for the scientific community to have a better understanding while selecting putative therapeutic protein candidates as druggable ones as effective measures against C. pseudotuberculosis.

Single-cell genomics-based analysis of virus-host interactions in marine surface bacterioplankton

Viral infections dynamically alter the composition and metabolic potential of marine microbial communities and the evolutionary trajectories of host populations with resulting feedback on biogeochemical cycles. It is quite possible that all microbial populations in the ocean are impacted by viral infections. Our knowledge of virus–host relationships, however, has been limited to a minute fraction of cultivated host groups. Here, we utilized single-cell sequencing to obtain genomic blueprints of viruses inside or attached to individual bacterial and archaeal cells captured in their native environment, circumventing the need for host and virus cultivation. A combination of comparative genomics, metagenomic fragment recruitment, sequence anomalies and irregularities in sequence coverage depth and genome recovery were utilized to detect viruses and to decipher modes of virus–host interactions. Members of all three tailed phage families were identified in 20 out of 58 phylogenetically and geographically diverse single amplified genomes (SAGs) of marine bacteria and archaea. At least four phage–host interactions had the characteristics of late lytic infections, all of which were found in metabolically active cells. One virus had genetic potential for lysogeny. Our findings include first known viruses of Thaumarchaeota, Marinimicrobia, Verrucomicrobia and Gammaproteobacteria clusters SAR86 and SAR92. Viruses were also found in SAGs of Alphaproteobacteria and Bacteroidetes. A high fragment recruitment of viral metagenomic reads confirmed that most of the SAG-associated viruses are abundant in the ocean. Our study demonstrates that single-cell genomics, in conjunction with sequence-based computational tools, enable in situ, cultivation-independent insights into host–virus interactions in complex microbial communities.

Pan-genome analysis of human gastric pathogen H. pylori: comparative genomics and pathogenomics approaches to identify regions associated with pathogenicity and prediction of potential core therapeutic targets

Helicobacter pylori is a human gastric pathogen implicated as the major cause of peptic ulcer and second leading cause of gastric cancer (~70%) around the world. Conversely, an increased resistance to antibiotics and hindrances in the development of vaccines against H. pylori are observed. Pan-genome analyses of the global representative H. pylori isolates consisting of 39 complete genomes are presented in this paper. Phylogenetic analyses have revealed close relationships among geographically diverse strains of H. pylori. The conservation among these genomes was further analyzed by pan-genome approach; the predicted conserved gene families (1,193) constitute ~77% of the average H. pylori genome and 45% of the global gene repertoire of the species. Reverse vaccinology strategies have been adopted to identify and narrow down the potential core-immunogenic candidates. Total of 28 nonhost homolog proteins were characterized as universal therapeutic targets against H. pylori based on their functional annotation and protein-protein interaction. Finally, pathogenomics and genome plasticity analysis revealed 3 highly conserved and 2 highly variable putative pathogenicity islands in all of the H. pylori genomes been analyzed.

Label-free proteomic analysis to confirm the predicted proteome of Corynebacterium pseudotuberculosis under nitrosative stress mediated by nitric oxide

Background

Corynebacterium pseudotuberculosis biovar ovis is a facultative intracellular pathogen, and the etiological agent of caseous lymphadenitis in small ruminants. During the infection process, the bacterium is subjected to several stress conditions, including nitrosative stress, which is caused by nitric oxide (NO). In silico analysis of the genome of C. pseudotuberculosis ovis 1002 predicted several genes that could influence the resistance of this pathogen to nitrosative stress. Here, we applied high-throughput proteomics using high definition mass spectrometry to characterize the functional genome of C. pseudotuberculosis ovis 1002 in the presence of NO-donor Diethylenetriamine/nitric oxide adduct (DETA/NO), with the aim of identifying proteins involved in nitrosative stress resistance.

Results

We characterized 835 proteins, representing approximately 41% of the predicted proteome of C. pseudotuberculosis ovis 1002, following exposure to nitrosative stress. In total, 102 proteins were exclusive to the proteome of DETA/NO-induced cells, and a further 58 proteins were differentially regulated between the DETA/NO and control conditions. An interactomic analysis of the differential proteome of C. pseudotuberculosis in response to nitrosative stress was also performed. Our proteomic data set suggested the activation of both a general stress response and a specific nitrosative stress response, as well as changes in proteins involved in cellular metabolism, detoxification, transcriptional regulation, and DNA synthesis and repair.

Conclusions

Our proteomic analysis validated previously-determined in silico data for C. pseudotuberculosis ovis 1002. In addition, proteomic screening performed in the presence of NO enabled the identification of a set of factors that can influence the resistance and survival of C. pseudotuberculosis during exposure to nitrosative stress.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1065) contains supplementary material, which is available to authorized users.

PAIDB v2.0: exploration and analysis of pathogenicity and resistance islands

Pathogenicity is a complex multifactorial process confounded by the concerted activity of genetic regions associated with virulence and/or resistance determinants. Pathogenicity islands (PAIs) and resistance islands (REIs) are key to the evolution of pathogens and appear to play complimentary roles in the process of bacterial infection. While PAIs promote disease development, REIs give a fitness advantage to the host against multiple antimicrobial agents. The Pathogenicity Island Database (PAIDB, http://www.paidb.re.kr) has been the only database dedicated to providing comprehensive information on all reported PAIs and candidate PAIs in prokaryotic genomes. In this study, we present PAIDB v2.0, whose functionality is extended to incorporate REIs. PAIDB v2.0 contains 223 types of PAIs with 1331 accessions, and 88 types of REIs with 108 accessions. With an improved detection scheme, 2673 prokaryotic genomes were analyzed to locate candidate PAIs and REIs. With additional quantitative and qualitative advancements in database content and detection accuracy, PAIDB will continue to facilitate pathogenomic studies of both pathogenic and non-pathogenic organisms.

Fine-tuned characterization of Staphylococcus aureus Newbould 305, a strain associated with mild and chronic mastitis in bovines

S. aureus is a major aetiological agent of ruminant mastitis worldwide. The chronic nature of S. aureus mastitis makes it difficult to cure and prone to resurgence. In order to identify the bacterial factors involved in this chronicity, Newbould 305 (N305), a strain that can reproducibly induce mild and chronic mastitis in an experimental setting, was characterized in depth. We employed genomic and proteomic techniques combined with phenotype characterization, in order to comprehensively analyse N305. The results were compared with data obtained on S. aureus RF122, a strain representative of the major clone involved in severe bovine mastitis worldwide. Five mobile genetic elements were identified in the N305 genome as carrying virulence factors which correlated with phenotypic features such as cytotoxicity, mammary epithelial cell invasion or host-adaptation. In particular, the presence and characteristics of surface exposed proteins correlated well with the greater adhesion and internalization capacities of N305 in bovine mammary epithelial cells. N305 also displayed less diversity of toxin genes but secreted larger quantities of these toxins, associated with a higher cytotoxicity potential. Our data are consistent with the invasiveness and host-adaptation features which contribute to the chronicity of S. aureus mastitis. Mobile genetic elements, exoproteins and surface exposed proteins constitute good targets for further research to explore the underlying mechanisms related to mastitis chronicity.

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.

Value of a newly sequenced bacterial genome

Next-generation sequencing (NGS) technologies have made high-throughput sequencing available to medium- and small-size laboratories, culminating in a tidal wave of genomic information. The quantity of sequenced bacterial genomes has not only brought excitement to the field of genomics but also heightened expectations that NGS would boost antibacterial discovery and vaccine development. Although many possible drug and vaccine targets have been discovered, the success rate of genome-based analysis has remained below expectations. Furthermore, NGS has had consequences for genome quality, resulting in an exponential increase in draft (partial data) genome deposits in public databases. If no further interests are expressed for a particular bacterial genome, it is more likely that the sequencing of its genome will be limited to a draft stage, and the painstaking tasks of completing the sequencing of its genome and annotation will not be undertaken. It is important to know what is lost when we settle for a draft genome and to determine the “scientific value” of a newly sequenced genome. This review addresses the expected impact of newly sequenced genomes on antibacterial discovery and vaccinology. Also, it discusses the factors that could be leading to the increase in the number of draft deposits and the consequent loss of relevant biological information.

Identifying pathogenicity islands in bacterial pathogenomics using computational approaches

High-throughput sequencing technologies have made it possible to study bacteria through analyzing their genome sequences. For instance, comparative genome sequence analyses can reveal the phenomenon such as gene loss, gene gain, or gene exchange in a genome. By analyzing pathogenic bacterial genomes, we can discover that pathogenic genomic regions in many pathogenic bacteria are horizontally transferred from other bacteria, and these regions are also known as pathogenicity islands (PAIs). PAIs have some detectable properties, such as having different genomic signatures than the rest of the host genomes, and containing mobility genes so that they can be integrated into the host genome. In this review, we will discuss various pathogenicity island-associated features and current computational approaches for the identification of PAIs. Existing pathogenicity island databases and related computational resources will also be discussed, so that researchers may find it to be useful for the studies of bacterial evolution and pathogenicity mechanisms.

Progression of 'OMICS' methodologies for understanding the pathogenicity of Corynebacterium pseudotuberculosis: the Brazilian experience

Since the first successful attempt at sequencing the Corynebacterium pseudotuberculosis genome, large amounts of genomic, transcriptomic and proteomic data have been generated. C. pseudotuberculosis is an interesting bacterium due to its great zoonotic potential and because it causes considerable economic losses worldwide. Furthermore, different strains of C. pseudotuberculosis are capable of causing various diseases in different hosts. Currently, we seek information about the phylogenetic relationships between different strains of C. pseudotuberculosis isolates from different hosts across the world and to employ these data to develop tools to diagnose and eradicate the diseases these strains cause. In this review, we present the latest findings on C. pseudotuberculosis that have been obtained with the most advanced techniques for sequencing and genomic organization. We also discuss the development of in silico tools for processing these data to prompt a better understanding of this pathogen.

Draft Genome Sequence of Pseudomonas plecoglossicida Strain NB2011, the Causative Agent of White Nodules in Large Yellow Croaker (Larimichthys crocea)

We describe the draft genome sequence of Pseudomonas plecoglossicida strain NB2011, the causative agent of white nodules in cultured large yellow croaker (Larimichthys crocea) in China. The draft genome sequence of the bacterium consists of 5.41 million bp, with a G+C content of 62.8%. A total of 4,952 genes were identified.

Complete genome analysis of three Acinetobacter baumannii clinical isolates in China for insight into the diversification of drug resistance elements

Background

The emergence and rapid spreading of multidrug-resistant Acinetobacter baumannii strains has become a major health threat worldwide. To better understand the genetic recombination related with the acquisition of drug-resistant elements during bacterial infection, we performed complete genome analysis on three newly isolated multidrug-resistant A. baumannii strains from Beijing using next-generation sequencing technology.

Methodologies/Principal Findings

Whole genome comparison revealed that all 3 strains share some common drug resistant elements including carbapenem-resistant bla_OXA-23 and tetracycline (tet) resistance islands, but the genome structures are diversified among strains. Various genomic islands intersperse on the genome with transposons and insertions, reflecting the recombination flexibility during the acquisition of the resistant elements. The blood-isolated BJAB07104 and ascites-isolated BJAB0868 exhibit high similarity on their genome structure with most of the global clone II strains, suggesting these two strains belong to the dominant outbreak strains prevalent worldwide. A large resistance island (RI) of about 121-kb, carrying a cluster of resistance-related genes, was inserted into the ATPase gene on BJAB07104 and BJAB0868 genomes. A 78-kb insertion element carrying tra-locus and bla_OXA-23 island, can be either inserted into one of the tniB gene in the 121-kb RI on the chromosome, or transformed to conjugative plasmid in the two BJAB strains. The third strains of this study, BJAB0715, which was isolated from spinal fluid, exhibit much more divergence compared with above two strains. It harbors multiple drug-resistance elements including a truncated AbaR-22-like RI on its genome. One of the unique features of this strain is that it carries both bla_OXA-23 and bla_OXA-58 genes on its genome. Besides, an Acinetobacter lwoffii adeABC efflux element was found inserted into the ATPase position in BJAB0715.

Conclusions

Our comparative analysis on currently completed Acinetobacter baumannii genomes revealed extensive and dynamic genome organizations, which may facilitate the bacteria to acquire drug-resistance elements into their genomes.

An overview of the domestication and impact of the Salmonella mobilome

Salmonella spp. are accountable for a large fraction of the global infectious disease burden, with most of their infections being food- or water-borne. The phenotypic features and adaptive potential of Salmonella spp. appear to be driven to a large extent by mobile or laterally acquired genetic elements. A better understanding of the conduct and diversification of these important pathogens consequently requires a more profound insight into the different mechanisms by which these pivotal elements establish themselves in the cell and affect its behavior. This review, therefore, provides an overview of the physiological impact and domestication of the Salmonella mobilome.

The pan-genome of the animal pathogen Corynebacterium pseudotuberculosis reveals differences in genome plasticity between the biovar ovis and equi strains

Corynebacterium pseudotuberculosis is a facultative intracellular pathogen and the causative agent of several infectious and contagious chronic diseases, including caseous lymphadenitis, ulcerative lymphangitis, mastitis, and edematous skin disease, in a broad spectrum of hosts. In addition, Corynebacterium pseudotuberculosis infections pose a rising worldwide economic problem in ruminants. The complete genome sequences of 15 C. pseudotuberculosis strains isolated from different hosts and countries were comparatively analyzed using a pan-genomic strategy. Phylogenomic, pan-genomic, core genomic, and singleton analyses revealed close relationships among pathogenic corynebacteria, the clonal-like behavior of C. pseudotuberculosis and slow increases in the sizes of pan-genomes. According to extrapolations based on the pan-genomes, core genomes and singletons, the C. pseudotuberculosis biovar ovis shows a more clonal-like behavior than the C. pseudotuberculosis biovar equi. Most of the variable genes of the biovar ovis strains were acquired in a block through horizontal gene transfer and are highly conserved, whereas the biovar equi strains contain great variability, both intra- and inter-biovar, in the 16 detected pathogenicity islands (PAIs). With respect to the gene content of the PAIs, the most interesting finding is the high similarity of the pilus genes in the biovar ovis strains compared with the great variability of these genes in the biovar equi strains. Concluding, the polymerization of complete pilus structures in biovar ovis could be responsible for a remarkable ability of these strains to spread throughout host tissues and penetrate cells to live intracellularly, in contrast with the biovar equi, which rarely attacks visceral organs. Intracellularly, the biovar ovis strains are expected to have less contact with other organisms than the biovar equi strains, thereby explaining the significant clonal-like behavior of the biovar ovis strains.

Tips and tricks for the assembly of a Corynebacterium pseudotuberculosis genome using a semiconductor sequencer

New sequencing platforms have enabled rapid decoding of complete prokaryotic genomes at relatively low cost. The Ion Torrent platform is an example of these technologies, characterized by lower coverage, generating challenges for the genome assembly. One particular problem is the lack of genomes that enable reference-based assembly, such as the one used in the present study, Corynebacterium pseudotuberculosis biovar equi, which causes high economic losses in the US equine industry. The quality treatment strategy incorporated into the assembly pipeline enabled a 16-fold greater use of the sequencing data obtained compared with traditional quality filter approaches. Data preprocessing prior to the de novo assembly enabled the use of known methodologies in the next-generation sequencing data assembly. Moreover, manual curation was proved to be essential for ensuring a quality assembly, which was validated by comparative genomics with other species of the genus Corynebacterium. The present study presents a modus operandi that enables a greater and better use of data obtained from semiconductor sequencing for obtaining the complete genome from a prokaryotic microorganism, C. pseudotuberculosis, which is not a traditional biological model such as Escherichia coli.

Genome sequence of Corynebacterium pseudotuberculosis biovar equi strain 258 and prediction of antigenic targets to improve biotechnological vaccine production

Corynebacterium pseudotuberculosis is the causative agent of several veterinary diseases in a broad range of economically important hosts, which can vary from caseous lymphadenitis in sheep and goats (biovar ovis) to ulcerative lymphangitis in cattle and horses (biovar equi). Existing vaccines against C. pseudotuberculosis are mainly intended for small ruminants and, even in these hosts, they still present remarkable limitations. In this study, we present the complete genome sequence of C. pseudotuberculosis biovar equi strain 258, isolated from a horse with ulcerative lymphangitis. The genome has a total size of 2,314,404 bp and contains 2088 predicted protein-coding regions. Using in silico analysis, eleven pathogenicity islands were detected in the genome sequence of C. pseudotuberculosis 258. The application of a reverse vaccinology strategy identified 49 putative antigenic proteins, which can be used as candidate vaccine targets in future works.

Contribution of phage-derived genomic islands to the virulence of facultative bacterial pathogens

Facultative pathogens have extremely dynamic pan-genomes, to a large extent derived from bacteriophages and other mobile elements. We developed a simple approach to identify phage-derived genomic islands and apply it to show that pathogens from diverse bacterial genera are significantly enriched in clustered phage-derived genes compared with related benign strains. These findings show that genome expansion by integration of prophages containing virulence factors is a major route of evolution of facultative bacterial pathogens.

Campylobacter fetus subspecies: comparative genomics and prediction of potential virulence targets

The genus Campylobacter contains pathogens causing a wide range of diseases, targeting both humans and animals. Among them, the Campylobacter fetus subspecies fetus and venerealis deserve special attention, as they are the etiological agents of human bacterial gastroenteritis and bovine genital campylobacteriosis, respectively. We compare the whole genomes of both subspecies to get insights into genomic architecture, phylogenetic relationships, genome conservation and core virulence factors. Pan-genomic approach was applied to identify the core- and pan-genome for both C. fetus subspecies and members of the genus. The C. fetus subspecies conserved (76%) proteome were then analyzed for their subcellular localization and protein functions in biological processes. Furthermore, with pathogenomic strategies, unique candidate regions in the genomes and several potential core-virulence factors were identified. The potential candidate factors identified for attenuation and/or subunit vaccine development against C. fetus subspecies contain: nucleoside diphosphate kinase (Ndk), type IV secretion systems (T4SS), outer membrane proteins (OMP), substrate binding proteins CjaA and CjaC, surface array proteins, sap gene, and cytolethal distending toxin (CDT). Significantly, many of those genes were found in genomic regions with signals of horizontal gene transfer and, therefore, predicted as putative pathogenicity islands. We found CRISPR loci and dam genes in an island specific for C. fetus subsp. fetus, and T4SS and sap genes in an island specific for C. fetus subsp. venerealis. The genomic variations and potential core and unique virulence factors characterized in this study would lead to better insight into the species virulence and to more efficient use of the candidates for antibiotic, drug and vaccine development.

Pangenomic study of Corynebacterium diphtheriae that provides insights into the genomic diversity of pathogenic isolates from cases of classical diphtheria, endocarditis, and pneumonia

Corynebacterium diphtheriae is one of the most prominent human pathogens and the causative agent of the communicable disease diphtheria. The genomes of 12 strains isolated from patients with classical diphtheria, endocarditis, and pneumonia were completely sequenced and annotated. Including the genome of C. diphtheriae NCTC 13129, we herewith present a comprehensive comparative analysis of 13 strains and the first characterization of the pangenome of the species C. diphtheriae. Comparative genomics showed extensive synteny and revealed a core genome consisting of 1,632 conserved genes. The pangenome currently comprises 4,786 protein-coding regions and increases at an average of 65 unique genes per newly sequenced strain. Analysis of prophages carrying the diphtheria toxin gene tox revealed that the toxoid vaccine producer C. diphtheriae Park-Williams no. 8 has been lysogenized by two copies of the ω(tox)(+) phage, whereas C. diphtheriae 31A harbors a hitherto-unknown tox(+) corynephage. DNA binding sites of the tox-controlling regulator DtxR were detected by genome-wide motif searches. Comparative content analysis showed that the DtxR regulons exhibit marked differences due to gene gain, gene loss, partial gene deletion, and DtxR binding site depletion. Most predicted pathogenicity islands of C. diphtheriae revealed characteristics of horizontal gene transfer. The majority of these islands encode subunits of adhesive pili, which can play important roles in adhesion of C. diphtheriae to different host tissues. All sequenced isolates contain at least two pilus gene clusters. It appears that variation in the distributed genome is a common strategy of C. diphtheriae to establish differences in host-pathogen interactions.