Nicholas H. Bergman

Interactive metagenomic visualization in a Web browser

BackgroundA critical output of metagenomic studies is the estimation of abundances of taxonomical or functional groups. The inherent uncertainty in assignments to these groups makes it important to consider both their hierarchical contexts and their prediction confidence. The current tools for visualizing metagenomic data, however, omit or distort quantitative hierarchical relationships and lack the facility for displaying secondary variables.ResultsHere we present Krona, a new visualization tool that allows intuitive exploration of relative abundances and confidences within the complex hierarchies of metagenomic classifications. Krona combines a variant of radial, space-filling displays with parametric coloring and interactive polar-coordinate zooming. The HTML5 and JavaScript implementation enables fully interactive charts that can be explored with any modern Web browser, without the need for installed software or plug-ins. This Web-based architecture also allows each chart to be an independent document, making them easy to share via e-mail or post to a standard Web server. To illustrate Kronas utility, we describe its application to various metagenomic data sets and its compatibility with popular metagenomic analysis tools.ConclusionsKrona is both a powerful metagenomic visualization tool and a demonstration of the potential of HTML5 for highly accessible bioinformatic visualizations. Its rich and interactive displays facilitate more informed interpretations of metagenomic analyses, while its implementation as a browser-based application makes it extremely portable and easily adopted into existing analysis packages. Both the Krona rendering code and conversion tools are freely available under a BSD open-source license, and available from: http://krona.sourceforge.net.

Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation

Long-read single-molecule sequencing has revolutionized de novo genome assembly and enabled the automated reconstruction of reference-quality genomes. However, given the relatively high error rates of such technologies, efficient and accurate assembly of large repeats and closely related haplotypes remains challenging. We address these issues with Canu, a successor of Celera Assembler that is specifically designed for noisy single-molecule sequences. Canu introduces support for nanopore sequencing, halves depth-of-coverage requirements, and improves assembly continuity while simultaneously reducing runtime by an order of magnitude on large genomes versus Celera Assembler 8.2. These advances result from new overlapping and assembly algorithms, including an adaptive overlapping strategy based on tf-idf weighted MinHash and a sparse assembly graph construction that avoids collapsing diverged repeats and haplotypes. We demonstrate that Canu can reliably assemble complete microbial genomes and near-complete eukaryotic chromosomes using either Pacific Biosciences (PacBio) or Oxford Nanopore technologies and achieves a contig NG50 of >21 Mbp on both human and Drosophila melanogaster PacBio data sets. For assembly structures that cannot be linearly represented, Canu provides graph-based assembly outputs in graphical fragment assembly (GFA) format for analysis or integration with complementary phasing and scaffolding techniques. The combination of such highly resolved assembly graphs with long-range scaffolding information promises the complete and automated assembly of complex genomes.

Reducing assembly complexity of microbial genomes with single-molecule sequencing.

BackgroundThe short reads output by first- and second-generation DNA sequencing instruments cannot completely reconstruct microbial chromosomes. Therefore, most genomes have been left unfinished due to the significant resources required to manually close gaps in draft assemblies. Third-generation, single-molecule sequencing addresses this problem by greatly increasing sequencing read length, which simplifies the assembly problem.ResultsTo measure the benefit of single-molecule sequencing on microbial genome assembly, we sequenced and assembled the genomes of six bacteria and analyzed the repeat complexity of 2,267 complete bacteria and archaea. Our results indicate that the majority of known bacterial and archaeal genomes can be assembled without gaps, at finished-grade quality, using a single PacBio RS sequencing library. These single-library assemblies are also more accurate than typical short-read assemblies and hybrid assemblies of short and long reads.ConclusionsAutomated assembly of long, single-molecule sequencing data reduces the cost of microbial finishing to

Structure and Complexity of a Bacterial Transcriptome

1,000 for most genomes, and future advances in this technology are expected to drive the cost lower. This is expected to increase the number of completed genomes, improve the quality of microbial genome databases, and enable high-fidelity, population-scale studies of pan-genomes and chromosomal organization.

Formation and Composition of the Bacillus anthracis Endospore

Although gene expression has been studied in bacteria for decades, many aspects of the bacterial transcriptome remain poorly understood. Transcript structure, operon linkages, and information on absolute abundance all provide valuable insights into gene function and regulation, but none has ever been determined on a genome-wide scale for any bacterium. Indeed, these aspects of the prokaryotic transcriptome have been explored on a large scale in only a few instances, and consequently little is known about the absolute composition of the mRNA population within a bacterial cell. Here we report the use of a high-throughput sequencing-based approach in assembling the first comprehensive, single-nucleotide resolution view of a bacterial transcriptome. We sampled the Bacillus anthracis transcriptome under a variety of growth conditions and showed that the data provide an accurate and high-resolution map of transcript start sites and operon structure throughout the genome. Further, the sequence data identified previously nonannotated regions with significant transcriptional activity and enhanced the accuracy of existing genome annotations. Finally, our data provide estimates of absolute transcript abundance and suggest that there is significant transcriptional heterogeneity within a clonal, synchronized bacterial population. Overall, our results offer an unprecedented view of gene expression and regulation in a bacterial cell.

Mash: fast genome and metagenome distance estimation using MinHash

The endospores of Bacillus anthracis are the infectious particles of anthrax. Spores are dormant bacterial morphotypes able to withstand harsh environments for decades, which contributes to their ability to be formulated and dispersed as a biological weapon. We monitored gene expression in B. anthracis during growth and sporulation using full genome DNA microarrays and matched the results against a comprehensive analysis of the mature anthrax spore proteome. A large portion (approximately 36%) of the B. anthracis genome is regulated in a growth phase-dependent manner, and this regulation is marked by five distinct waves of gene expression as cells proceed from exponential growth through sporulation. The identities of more than 750 proteins present in the spore were determined by multidimensional chromatography and tandem mass spectrometry. Comparison of data sets revealed that while the genes responsible for assembly and maturation of the spore are tightly regulated in discrete stages, many of the components ultimately found in the spore are expressed throughout and even before sporulation, suggesting that gene expression during sporulation may be mainly related to the physical construction of the spore, rather than synthesis of eventual spore content. The spore also contains an assortment of specialized, but not obviously related, metabolic and protective proteins. These findings contribute to our understanding of spore formation and function and will be useful in the detection, prevention, and early treatment of anthrax. This study also highlights the complementary nature of genomic and proteomic analyses and the benefits of combining these approaches in a single study.

Transcriptional Profiling of the Bacillus anthracis Life Cycle In Vitro and an Implied Model for Regulation of Spore Formation

Mash extends the MinHash dimensionality-reduction technique to include a pairwise mutation distance and P value significance test, enabling the efficient clustering and search of massive sequence collections. Mash reduces large sequences and sequence sets to small, representative sketches, from which global mutation distances can be rapidly estimated. We demonstrate several use cases, including the clustering of all 54,118 NCBI RefSeq genomes in 33 CPU h; real-time database search using assembled or unassembled Illumina, Pacific Biosciences, and Oxford Nanopore data; and the scalable clustering of hundreds of metagenomic samples by composition. Mash is freely released under a BSD license (https://github.com/marbl/mash).

Biosynthetic Analysis of the Petrobactin Siderophore Pathway from Bacillus anthracis

The life cycle of Bacillus anthracis includes both vegetative and endospore morphologies which alternate based on nutrient availability, and there is considerable evidence indicating that the ability of this organism to cause anthrax depends on its ability to progress through this life cycle in a regulated manner. Here we report the use of a custom B. anthracis GeneChip in defining the gene expression patterns that occur throughout the entire life cycle in vitro. Nearly 5,000 genes were expressed in five distinct waves of transcription as the bacteria progressed from germination through sporulation, and we identified a specific set of functions represented within each wave. We also used these data to define the temporal expression of the spore proteome, and in doing so we have demonstrated that much of the spores protein content is not synthesized de novo during sporulation but rather is packaged from preexisting stocks. We explored several potential mechanisms by which the cell could control which proteins are packaged into the developing spore, and our analyses were most consistent with a model in which B. anthracis regulates the composition of the spore proteome based on protein stability. This study is by far the most comprehensive survey yet of the B. anthracis life cycle and serves as a useful resource in defining the growth-phase-dependent expression patterns of each gene. Additionally, the data and accompanying bioinformatics analyses suggest a model for sporulation that has broad implications for B. anthracis biology and offer new possibilities for microbial forensics and detection.

Transcriptional Profiling of Bacillus anthracis during Infection of Host Macrophages

The asbABCDEF gene cluster from Bacillus anthracis is responsible for biosynthesis of petrobactin, a catecholate siderophore that functions in both iron acquisition and virulence in a murine model of anthrax. We initiated studies to determine the biosynthetic details of petrobactin assembly based on mutational analysis of the asb operon, identification of accumulated intermediates, and addition of exogenous siderophores to asb mutant strains. As a starting point, in-frame deletions of each of the genes in the asb locus (asbABCDEF) were constructed. The individual mutations resulted in complete abrogation of petrobactin biosynthesis when strains were grown on iron-depleted medium. However, in vitro analysis showed that each asb mutant grew to a very limited extent as vegetative cells in iron-depleted medium. In contrast, none of the B. anthracis asb mutant strains were able to outgrow from spores under the same culture conditions. Provision of exogenous petrobactin was able to rescue the growth defect in each asb mutant strain. Taken together, these data provide compelling evidence that AsbA performs the penultimate step in the biosynthesis of petrobactin, involving condensation of 3,4-dihydroxybenzoyl spermidine with citrate to form 3,4-dihydroxybenzoyl spermidinyl citrate. As a final step, the data reveal that AsbB catalyzes condensation of a second molecule of 3,4-dihydroxybenzoyl spermidine with 3,4-dihydroxybenzoyl spermidinyl citrate to form the mature siderophore. This work sets the stage for detailed biochemical studies with this unique acyl carrier protein-dependent, nonribosomal peptide synthetase-independent biosynthetic system.

Bacillus anthracis genome organization in light of whole transcriptome sequencing

ABSTRACT The interaction between Bacillus anthracis and the mammalian phagocyte is one of the central stages in the progression of inhalational anthrax, and it is commonly believed that the host cell plays a key role in facilitating germination and dissemination of inhaled B. anthracis spores. Given this, a detailed definition of the survival strategies used by B. anthracis within the phagocyte is critical for our understanding of anthrax. In this study, we report the first genome-wide analysis of B. anthracis gene expression during infection of host phagocytes. We developed a technique for specific isolation of bacterial RNA from within infected murine macrophages, and we used custom B. anthracis microarrays to characterize the expression patterns occurring within intracellular bacteria throughout infection of the host phagocyte. We found that B. anthracis adapts very quickly to the intracellular environment, and our analyses identified metabolic pathways that appear to be important to the bacterium during intracellular growth, as well as individual genes that show significant induction in vivo. We used quantitative reverse transcription-PCR to verify that the expression trends that we observed by microarray analysis were valid, and we chose one gene (GBAA1941, encoding a putative transcriptional regulator) for further characterization. A deletion strain missing this gene showed no phenotype in vitro but was significantly attenuated in a mouse model of inhalational anthrax, suggesting that the microarray data described here provide not only the first comprehensive view of how B. anthracis survives within the host cell but also a number of promising leads for further research in anthrax.