ArrayOme: a program for estimating the sizes of microarray-visualized bacterial genomes

VRprofile2: detection of antibiotic resistance-associated mobilome in bacterial pathogens

VRprofile2 is an updated pipeline that rapidly identifies diverse mobile genetic elements in bacterial genome sequences. Compared with the previous version, three major improvements were made. First, the user-friendly visualization could aid users in investigating the antibiotic resistance gene cassettes in conjunction with various mobile elements in the multiple resistance region with mosaic structure. VRprofile2 could compare the predicted mobile elements to the collected known mobile elements with similar architecture. A new mobilome indicator was proposed to give an overall estimation of the mobilome size in individual bacterial genomes. Second, the relationship between antibiotic resistance genes, mobile elements, and host strains would be efficiently examined with the aid of predicted strain's sequence typing, the incompatibility group and the transferability of plasmids. Finally, the updated back-end database, MobilomeDB2, now collected nearly a thousand active mobile elements retrieved from literature or based on prediction. The pre-computed results of the antibiotic resistance gene-carrying mobile elements of >5500 ESKAPEE genomes were also provided. We expect that VRprofile2 will provide better support for researchers interested in bacterial mobile elements and the dissemination of antibiotic resistance. VRprofile2 is freely available to all users without any login requirement at https://tool2-mml.sjtu.edu.cn/VRprofile.

mGenomeSubtractor: a web-based tool for parallel in silico subtractive hybridization analysis of multiple bacterial genomes

mGenomeSubtractor performs an mpiBLAST-based comparison of reference bacterial genomes against multiple user-selected genomes for investigation of strain variable accessory regions. With parallel computing architecture, mGenomeSubtractor is able to run rapid BLAST searches of the segmented reference genome against multiple subject genomes at the DNA or amino acid level within a minute. In addition to comparison of protein coding sequences, the highly flexible sliding window-based genome fragmentation approach offered can be used to identify short unique sequences within or between genes. mGenomeSubtractor provides powerful schematic outputs for exploration of identified core and accessory regions, including searches against databases of mobile genetic elements, virulence factors or bacterial essential genes, examination of G+C content and binucleotide distribution bias, and integrated primer design tools. mGenomeSubtractor also allows for the ready definition of species-specific gene pools based on available genomes. Pan-genomic arrays can be easily developed using the efficient oligonucleotide design tool. This simple high-throughput in silico ‘subtractive hybridization’ analytical tool will support the rapidly escalating number of comparative bacterial genomics studies aimed at defining genomic biomarkers of evolutionary lineage, phenotype, pathotype, environmental adaptation and/or disease-association of diverse bacterial species. mGenomeSubtractor is freely available to all users without any login requirement at: http://bioinfo-mml.sjtu.edu.cn/mGS/.

Translational genomics to develop a Salmonella enterica serovar Paratyphi A multiplex polymerase chain reaction assay

The use of pathogen genome sequence data for the control and management of infections remains an ongoing challenge. We describe a broadly applicable, web-enabled approach that can be used to develop bacteria-specific polymerase chain reaction (PCR) assays. Salmonella enterica Paratyphi A has emerged as a major cause of enteric fever in Asia. Culture-based diagnosis is slow and frequently negative in patients with suspected typhoid and paratyphoid fever, potentially compromising patient management and public health. We used the MobilomeFINDER web-server to perform in silico subtractive hybridization, thus identifying 43 protein-coding sequences (CDSs) that were present in two Paratyphi A strains but not in other sequenced Salmonella genomes. After exclusion of 29 CDSs found to be variably present in Paratyphi A strains by microarray hybridization and grouping of remaining CDSs by genomic location, four dispersed targets (stkF, spa2473, spa2539, hsdM) were used to develop a highly discriminatory multiplex PCR assay. All 52 Paratyphi A strains within the diverse panel investigated produced one of two pathognomonic four-band signatures. Given rapid and ongoing expansion of DNA and comparative genomics databases, our universally accessible web-server-supported do-it-yourself approach offers the potential to contribute significantly to the rapid development of species-, serovar-, or pathotype-specific PCR assays targeting pre-existing and emerging bacterial pathogens.

Analysis of a genomic island housing genes for DNA S-modification system in Streptomyces lividans 66 and its counterparts in other distantly related bacteria

The complete sequence (92 770 bp) of a genomic island (GI) named SLG from Streptomyces lividans 66, encoding a novel DNA S-modification system (dnd), was determined. Its overall G+C content was 67.8%, lower than those of three sequenced Streptomyces genomes. Among 85 predicted open reading frames (ORFs) in SLG, 22 ORFs showed little homology with previously known proteins. SLG displays a mosaic structure composed of four modules, indicative of multiple recombination events in its formation. Spontaneous excision and circularization of SLG was observed, and the excision rate appeared to be induced at least fivefold by MNNG exposure. Using constructed mini-islands of SLG, we demonstrated that Slg01, a P4-like integrase, was sufficient to promote SLG integration, excision and circularization. Eleven counterpart dnd clusters, which also mapped to GIs in 10 chromosomes and a plasmid, were found in taxonomically unrelated bacterial species from various geographic niches. Additionally, c. 10% of actinomycetes were found to possess a dnd cluster in a survey involving 74 strains. Comparison of dnd clusters in the 12 bacteria strongly suggests that these dnd-bearing elements might have evolved from a common ancestor similar to plasmid-originated chromosome II of Pseudoalteromonas haloplanktis TAC125.

MobilomeFINDER: web-based tools for in silico and experimental discovery of bacterial genomic islands

MobilomeFINDER (http://mml.sjtu.edu.cn/MobilomeFINDER) is an interactive online tool that facilitates bacterial genomic island or 'mobile genome' (mobilome) discovery; it integrates the ArrayOme and tRNAcc software packages. ArrayOme utilizes a microarray-derived comparative genomic hybridization input data set to generate 'inferred contigs' produced by merging adjacent genes classified as 'present'. Collectively these 'fragments' represent a hypothetical 'microarray-visualized genome (MVG)'. ArrayOme permits recognition of discordances between physical genome and MVG sizes, thereby enabling identification of strains rich in microarray-elusive novel genes. Individual tRNAcc tools facilitate automated identification of genomic islands by comparative analysis of the contents and contexts of tRNA sites and other integration hotspots in closely related sequenced genomes. Accessory tools facilitate design of hotspot-flanking primers for in silico and/or wet-science-based interrogation of cognate loci in unsequenced strains and analysis of islands for features suggestive of foreign origins; island-specific and genome-contextual features are tabulated and represented in schematic and graphical forms. To date we have used MobilomeFINDER to analyse several Enterobacteriaceae, Pseudomonas aeruginosa and Streptococcus suis genomes. MobilomeFINDER enables high-throughput island identification and characterization through increased exploitation of emerging sequence data and PCR-based profiling of unsequenced test strains; subsequent targeted yeast recombination-based capture permits full-length sequencing and detailed functional studies of novel genomic islands.

Prediction of effective genome size in metagenomic samples

A novel computational approach shows a link between genome size and habitat from analysis of environmental metagenomic DNA reads.

We introduce a novel computational approach to predict effective genome size (EGS; a measure that includes multiple plasmid copies, inserted sequences, and associated phages and viruses) from short sequencing reads of environmental genomics (or metagenomics) projects. We observe considerable EGS differences between environments and link this with ecologic complexity as well as species composition (for instance, the presence of eukaryotes). For example, we estimate EGS in a complex, organism-dense farm soil sample at about 6.3 megabases (Mb) whereas that of the bacteria therein is only 4.7 Mb; for bacteria in a nutrient-poor, organism-sparse ocean surface water sample, EGS is as low as 1.6 Mb. The method also permits evaluation of completion status and assembly bias in single-genome sequencing projects.

Survey of genomic diversity among Enterococcus faecalis strains by microarray-based comparative genomic hybridization

We have compared nine Enterococcus faecalis strains with E. faecalis V583 by comparative genomic hybridization using microarrays (CGH). The strains used in this study (the "test" strains) originated from various environments. CGH is a powerful and promising tool for obtaining novel information on genome diversity in bacteria. By CGH, one obtains clues about which genes are present or divergent in the strains, compared to a reference strain (here, V583). The information obtained by CGH is important from both ecological and systematic points of view. CGH of E. faecalis showed considerable diversity in gene content: Compared to V583, the percentage of divergent genes in the test strains varied from 15% to 23%, and 154 genes were divergent in all strains. The main variation was found in regions corresponding to exogenously acquired or mobile DNA in V583. Antibiotic resistance genes, virulence factors, and integrated plasmid genes dominated among the divergent genes. The strains examined showed various contents of genes corresponding to the pTEF1, pTEF2, and pTEF3 genes in V583. The extensive transport and metabolic capabilities of V583 appeared similar in the test strains; CGH indicated that the ability to transport and metabolize various carbohydrates was similar in the test strains (verified by API 50 CH assays). The contents of genes related to stress tolerance appeared similar in V583 and the nine test strains, supporting the view of E. faecalis as an organism able to resist harsh conditions.

Development of a diagnostic test for the Midlands 1 cystic fibrosis epidemic strain of Pseudomonas aeruginosa

In a previous study of isolates from cystic fibrosis (CF) patients in England and Wales, the Midlands 1 strain of Pseudomonas aeruginosa was identified as the second most common clone, representing 10% of isolates and found in nearly one-third of all CF centres [Scott, F. W. & Pitt, T. L. (2004). J Med Microbiol 53, 609-615]. Using suppression subtractive hybridization, 54 sequences were identified as present in a Midlands 1 strain but were absent from strain PAO1. The distribution of 14 of these sequences amongst representatives of Midlands 1, other CF epidemic strains and unrelated P. aeruginosa CF isolates was determined using PCR assays. Using these data, a PCR-based test was developed that was specific for the Midlands 1 clone, which was confirmed using dot-blot hybridization. By applying the test to CF isolates from a CF centre in Liverpool, a Midlands 1 clone was identified. The identity was confirmed using typing by PFGE. The PCR test should facilitate a greater understanding of the distribution of the Midlands 1 strain in the UK and elsewhere.