Jochen Blom

EDGAR 2.0: an enhanced software platform for comparative gene content analyses

The rapidly increasing availability of microbial genome sequences has led to a growing demand for bioinformatics software tools that support the functional analysis based on the comparison of closely related genomes. By utilizing comparative approaches on gene level it is possible to gain insights into the core genes which represent the set of shared features for a set of organisms under study. Vice versa singleton genes can be identified to elucidate the specific properties of an individual genome. Since initial publication, the EDGAR platform has become one of the most established software tools in the field of comparative genomics. Over the last years, the software has been continuously improved and a large number of new analysis features have been added. For the new version, EDGAR 2.0, the gene orthology estimation approach was newly designed and completely re-implemented. Among other new features, EDGAR 2.0 provides extended phylogenetic analysis features like AAI (Average Amino Acid Identity) and ANI (Average Nucleotide Identity) matrices, genome set size statistics and modernized visualizations like interactive synteny plots or Venn diagrams. Thereby, the software supports a quick and user-friendly survey of evolutionary relationships between microbial genomes and simplifies the process of obtaining new biological insights into their differential gene content. All features are offered to the scientific community via a web-based and therefore platform-independent user interface, which allows easy browsing of precomputed datasets. The web server is accessible at http://edgar.computational.bio.

Whole Genome Sequencing versus Traditional Genotyping for Investigation of a Mycobacterium tuberculosis Outbreak: A Longitudinal Molecular Epidemiological Study

In an outbreak investigation of Mycobacterium tuberculosis comparing whole genome sequencing (WGS) with traditional genotyping, Stefan Niemann and colleagues found that classical genotyping falsely clustered some strains, and WGS better reflected contact tracing.

Whole-genome comparison of disease and carriage strains provides insights into virulence evolution in Neisseria meningitidis

Neisseria meningitidis is a leading cause of infectious childhood mortality worldwide. Most research efforts have hitherto focused on disease isolates belonging to only a few hypervirulent clonal lineages. However, up to 10% of the healthy human population is temporarily colonized by genetically diverse strains mostly with little or no pathogenic potential. Currently, little is known about the biology of carriage strains and their evolutionary relationship with disease isolates. The expression of a polysaccharide capsule is the only trait that has been convincingly linked to the pathogenic potential of N. meningitidis. To gain insight into the evolution of virulence traits in this species, whole-genome sequences of three meningococcal carriage isolates were obtained. Gene content comparisons with the available genome sequences from three disease isolates indicate that there is no core pathogenome in N. meningitidis. A comparison of the chromosome structure suggests that a filamentous prophage has mediated large chromosomal rearrangements and the translocation of some candidate virulence genes. Interspecific comparison of the available Neisseria genome sequences and dot blot hybridizations further indicate that the insertion sequence IS1655 is restricted only to N. meningitidis; its low sequence diversity is an indicator of an evolutionarily recent population bottleneck. A genome-based phylogenetic reconstruction provides evidence that N. meningitidis has emerged as an unencapsulated human commensal from a common ancestor with Neisseria gonorrhoeae and Neisseria lactamica and consecutively acquired the genes responsible for capsule synthesis via horizontal gene transfer.

Relationship of Bacillus amyloliquefaciens clades associated with strains DSM 7T and FZB42T: a proposal for Bacillus amyloliquefaciens subsp. amyloliquefaciens subsp. nov. and Bacillus amyloliquefaciens subsp. plantarum subsp. nov. based on complete genome sequence comparisons

The whole-genome-sequenced rhizobacterium Bacillus amyloliquefaciens FZB42(T) (Chen et al., 2007) and other plant-associated strains of the genus Bacillus described as belonging to the species Bacillus amyloliquefaciens or Bacillus subtilis are used commercially to promote the growth and improve the health of crop plants. Previous investigations revealed that a group of strains represented a distinct ecotype related to B. amyloliquefaciens; however, the exact taxonomic position of this group remains elusive (Reva et al., 2004). In the present study, we demonstrated the ability of a group of Bacillus strains closely related to strain FZB42(T) to colonize Arabidopsis roots. On the basis of their phenotypic traits, the strains were similar to Bacillus amyloliquefaciens DSM 7(T) but differed considerably from this type strain in the DNA sequences of genes encoding 16S rRNA, gyrase subunit A (gyrA) and histidine kinase (cheA). Phylogenetic analysis performed with partial 16S rRNA, gyrA and cheA gene sequences revealed that the plant-associated strains of the genus Bacillus, including strain FZB42(T), formed a lineage, which could be distinguished from the cluster of strains closely related to B. amyloliquefaciens DSM 7(T). DNA-DNA hybridizations (DDH) performed with genomic DNA from strains DSM 7(T) and FZB42(T) yielded relatedness values of 63.7-71.2 %. Several methods of genomic analysis, such as direct whole-genome comparison, digital DDH and microarray-based comparative genomichybridization (M-CGH) were used as complementary tests. The group of plant-associated strains could be distinguished from strain DSM 7(T) and the type strain of B. subtilis by differences in the potential to synthesize non-ribosomal lipopeptides and polyketides. Based on the differences found in the marker gene sequences and the whole genomes of these strains, we propose two novel subspecies, designated B. amyloliquefaciens subsp. plantarum subsp. nov., with the type strain FZB42(T) ( = DSM 23117(T) = BGSC 10A6(T)), and B. amyloliquefaciens subsp. amyloliquefaciens subsp. nov., with the type strain DSM 7(T)( = ATCC 23350(T) = Fukumoto Strain F(T)), for plant-associated and non-plant-associated representatives, respecitvely. This is in agreement with results of DDH and M-CGH tests and the MALDI-TOF MS of cellular components, all of which suggested that the ecovars represent two different subspecies.

Complete genome sequence of the fire blight pathogen Erwinia amylovora CFBP 1430 and comparison to other Erwinia spp.

Fire blight, caused by the enterobacterium Erwinia amylovora, is a devastating disease of rosaceous plants that has global economic importance for apple and pear production and trade. The complete genome of E. amylovora CFBP 1430 was sequenced, annotated, and compared with the genomes of other Erwinia spp. Several singleton and shared features of the E. amylovora CFBP 1430 genome were identified that offer a first view into evolutionary aspects within the genus Erwinia. Comparative genomics identified or clarified virulence and fitness determinants and secretion systems. Novel insights revealed in the genome of E. amylovora CFBP 1430 hold potential for exploitation to improve the design of more effective fire blight control strategies.

Complete genome sequencing of Agrobacterium sp. H13-3, the former Rhizobium lupini H13-3, reveals a tripartite genome consisting of a circular and a linear chromosome and an accessory plasmid but lacking a tumor-inducing Ti-plasmid

Agrobacterium sp. H13-3, formerly known as Rhizobium lupini H13-3, is a soil bacterium that was isolated from the rhizosphere of Lupinus luteus. The isolate has been established as a model system for studying novel features of flagellum structure, motility and chemotaxis within the family Rhizobiaceae. The complete genome sequence of Agrobacterium sp. H13-3 has been established and the genome structure and phylogenetic assignment of the organism was analysed. For de novo sequencing of the Agrobacterium sp. H13-3 genome, a combined strategy comprising 454-pyrosequencing on the Genome Sequencer FLX platform and PCR-based amplicon sequencing for gap closure was applied. The finished genome consists of three replicons and comprises 5,573,770 bases. Based on phylogenetic analyses, the isolate could be assigned to the genus Agrobacterium biovar I and represents a genomic species G1 strain within this biovariety. The highly conserved circular chromosome (2.82 Mb) of Agrobacterium sp. H13-3 mainly encodes housekeeping functions characteristic for an aerobic, heterotrophic bacterium. Agrobacterium sp. H13-3 is a motile bacterium driven by the rotation of several complex flagella. Its behaviour towards external stimuli is regulated by a large chemotaxis regulon and a total of 17 chemoreceptors. Comparable to the genome of Agrobacterium tumefaciens C58, Agrobacterium sp. H13-3 possesses a linear chromosome (2.15 Mb) that is related to its reference replicon and features chromosomal and plasmid-like properties. The accessory plasmid pAspH13-3a (0.6 Mb) is only distantly related to the plasmid pAtC58 of A. tumefaciens C58 and shows a mosaic structure. A tumor-inducing Ti-plasmid is missing in the sequenced strain H13-3 indicating that it is a non-virulent isolate.

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.

Exact and complete short-read alignment to microbial genomes using Graphics Processing Unit programming

MOTIVATION The introduction of next-generation sequencing techniques and especially the high-throughput systems Solexa (Illumina Inc.) and SOLiD (ABI) made the mapping of short reads to reference sequences a standard application in modern bioinformatics. Short-read alignment is needed for reference based re-sequencing of complete genomes as well as for gene expression analysis based on transcriptome sequencing. Several approaches were developed during the last years allowing for a fast alignment of short sequences to a given template. Methods available to date use heuristic techniques to gain a speedup of the alignments, thereby missing possible alignment positions. Furthermore, most approaches return only one best hit for every query sequence, thus losing the potentially valuable information of alternative alignment positions with identical scores. RESULTS We developed SARUMAN (Semiglobal Alignment of short Reads Using CUDA and NeedleMAN-Wunsch), a mapping approach that returns all possible alignment positions of a read in a reference sequence under a given error threshold, together with one optimal alignment for each of these positions. Alignments are computed in parallel on graphics hardware, facilitating an considerable speedup of this normally time-consuming step. Combining our filter algorithm with CUDA-accelerated alignments, we were able to align reads to microbial genomes in time comparable or even faster than all published approaches, while still providing an exact, complete and optimal result. At the same time, SARUMAN runs on every standard Linux PC with a CUDA-compatible graphics accelerator. AVAILABILITY http://www.cebitec.uni-bielefeld.de/brf/saruman/saruman.html.

Comparative genomics of Lactobacillus.

The genus Lactobacillus includes a diverse group of bacteria consisting of many species that are associated with fermentations of plants, meat or milk. In addition, various lactobacilli are natural inhabitants of the intestinal tract of humans and other animals. Finally, several Lactobacillus strains are marketed as probiotics as their consumption can confer a health benefit to host. Presently, 154 Lactobacillus species are known and a growing fraction of these are subject to draft genome sequencing. However, complete genome sequences are needed to provide a platform for detailed genomic comparisons. Therefore, we selected a total of 20 genomes of various Lactobacillus strains for which complete genomic sequences have been reported. These genomes had sizes varying from 1.8 to 3.3 Mb and other characteristic features, such as G+C content that ranged from 33% to 51%. The Lactobacillus pan genome was found to consist of approximately 14 000 protein‐encoding genes while all 20 genomes shared a total of 383 sets of orthologous genes that defined the Lactobacillus core genome (LCG). Based on advanced phylogeny of the proteins encoded by this LCG, we grouped the 20 strains into three main groups and defined core group genes present in all genomes of a single group, signature group genes shared in all genomes of one group but absent in all other Lactobacillus genomes, and Group‐specific ORFans present in core group genes of one group and absent in all other complete genomes. The latter are of specific value in defining the different groups of genomes. The study provides a platform for present individual comparisons as well as future analysis of new Lactobacillus genomes.

Genome sequence of B. amyloliquefaciens type strain DSM7(T) reveals differences to plant-associated B. amyloliquefaciens FZB42.

The complete genome sequence of Bacillus amyloliquefaciens type strain DSM7(T) is presented. A comparative analysis between the genome sequences of the plant associated strain FZB42 (Chen et al., 2007) with the genome of B. amyloliquefaciens DSM7(T) revealed obvious differences in the variable part of the genomes, whilst the core genomes were found to be very similar. The strains FZB42 and DSM7(T) have in common 3345 genes (CDS) in their core genomes; whilst 547 and 344 CDS were found to be unique in DSM7(T) and FZB42, respectively. The core genome shared by both strains exhibited 97.89% identity on amino acid level. The number of genes representing the core genome of the strains FZB42, DSM7(T), and Bacillus subtilis DSM10(T) was calculated as being 3098 and their identity was 92.25%. The 3,980,199 bp genome of DSM7(T) contains numerous genomic islands (GI) detected by different methods. Many of them were located in vicinity of tRNA, glnA, and glmS gene copies. In contrast to FZB42, but similar to B. subtilis DSM10(T), the GI were enriched in prophage sequences and often harbored transposases, integrases and recombinases. Compared to FZB42, B. amyloliquefaciens DSM7(T) possessed a reduced potential to non-ribosomally synthesize secondary metabolites with antibacterial and/or antifungal action. B. amyloliquefaciens DSM7(T) did not produce the polyketides difficidin and macrolactin and was impaired in its ability to produce lipopeptides other than surfactin. Differences established within the variable part of the genomes, justify our proposal to discriminate the plant-associated ecotype represented by FZB42 from the group of type strain related B. amyloliquefaciens soil bacteria.