Comparative genomics of 274 Vibrio cholerae genomes reveals mobile functions structuring three niche dimensions.

BACKGROUND: Vibrio cholerae is a globally dispersed pathogen that has evolved with humans for centuries, but also includes non-pathogenic environmental strains. Here, we identify the genomic variability underlying this remarkable persistence across the three major niche dimensions space, time, and habitat. RESULTS: Taking an innovative approach of genome-wide association applicable to microbial genomes (GWAS-M), we classify 274 complete V. cholerae genomes by niche, including 39 newly sequenced for this study with the Ion Torrent DNA-sequencing platform. Niche metadata were collected for each strain and analyzed together with comprehensive annotations of genetic and genomic attributes, including point mutations (single-nucleotide polymorphisms, SNPs), protein families, functions and prophages. CONCLUSIONS: Our analysis revealed that genomic variations, in particular mobile functions including phages, prophages, transposable elements, and plasmids underlie the metadata structuring in each of the three niche dimensions. This underscores the role of phages and mobile elements as the most rapidly evolving elements in bacterial genomes, creating local endemicity (space), leading to temporal divergence (time), and allowing the invasion of new habitats. Together, we take a data-driven approach for comparative functional genomics that exploits high-volume genome sequencing and annotation, in conjunction with novel statistical and machine learning analyses to identify connections between genotype and phenotype on a genome-wide scale.

Clinical insights from metagenomic analysis of sputum samples from patients with cystic fibrosis.

As DNA sequencing becomes faster and cheaper, genomics-based approaches are being explored for their use in personalized diagnoses and treatments. Here, we provide a proof of principle for disease monitoring using personal metagenomic sequencing and traditional clinical microbiology by focusing on three adults with cystic fibrosis (CF). The CF lung is a dynamic environment that hosts a complex ecosystem composed of bacteria, viruses, and fungi that can vary in space and time. Not surprisingly, the microbiome data from the induced sputum samples we collected revealed a significant amount of species diversity not seen in routine clinical laboratory cultures. The relative abundances of several species changed as clinical treatment was altered, enabling the identification of the climax and attack communities that were proposed in an earlier work. All patient microbiomes encoded a diversity of mechanisms to resist antibiotics, consistent with the characteristics of multidrug-resistant microbial communities that are commonly observed in CF patients. The metabolic potentials of these communities differed by the health status and recovery route of each patient. Thus, this pilot study provides an example of how metagenomic data might be used with clinical assessments for the development of treatments tailored to individual patients.

Metagenomic detection of phage-encoded platelet-binding factors in the human oral cavity.

The human oropharynx is a reservoir for many potential pathogens, including streptococcal species that cause endocarditis. Although oropharyngeal microbes have been well described, viral communities are essentially uncharacterized. We conducted a metagenomic study to determine the composition of oropharyngeal DNA viral communities (both phage and eukaryotic viruses) in healthy individuals and to evaluate oropharyngeal swabs as a rapid method for viral detection. Viral DNA was extracted from 19 pooled oropharyngeal swabs and sequenced. Viral communities consisted almost exclusively of phage, and complete genomes of several phage were recovered, including Escherichia coli phage T3, Propionibacterium acnes phage PA6, and Streptococcus mitis phage SM1. Phage relative abundances changed dramatically depending on whether samples were chloroform treated or filtered to remove microbial contamination. pblA and pblB genes of phage SM1 were detected in the metagenomes. pblA and pblB mediate the attachment of S. mitis to platelets and play a significant role in S. mitis virulence in the endocardium, but have never previously been detected in the oral cavity. These genes were also identified in salivary metagenomes from three individuals at three time points and in individual saliva samples by PCR. Additionally, we demonstrate that phage SM1 can be induced by commonly ingested substances. Our results indicate that the oral cavity is a reservoir for pblA and pblB genes and for phage SM1 itself. Further studies will determine the association between pblA and pblB genes in the oral cavity and the risk of endocarditis.

Distinct lineage of vesiculovirus from big brown bats, United States.

We identified a novel rhabdovirus, American bat vesiculovirus, from postmortem tissue samples from 120 rabies-negative big brown bats with a history of human contact. Five percent of the tested bats were infected with this virus. The extent of zoonotic exposure and possible health effects in humans from this virus are unknown.

Identification and removal of ribosomal RNA sequences from metatranscriptomes.

SUMMARY: Here, we present riboPicker, a robust framework for the rapid, automated identification and removal of ribosomal RNA sequences from metatranscriptomic datasets. The results can be exported for subsequent analysis, and the databases used for the web-based version are updated on a regular basis. riboPicker categorizes rRNA-like sequences and provides graphical visualizations and tabular outputs of ribosomal coverage, alignment results and taxonomic classifications. AVAILABILITY AND IMPLEMENTATION: This open-source application was implemented in Perl and can be used as stand-alone version or accessed online through a user-friendly web interface. The source code, user help and additional information is available at http://ribopicker.sourceforge.net/. CONTACT: rschmied@sciences.sdsu.edu; rschmied@sciences.sdsu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SEQanswers: an open access community for collaboratively decoding genomes.

SUMMARY: The affordability of high-throughput sequencing has created an unprecedented surge in the use of genomic data in basic, translational and clinical research. The rapid evolution of sequencing technology, coupled with its broad adoption across biology and medicine, necessitates fast, collaborative interdisciplinary discussion. SEQanswers provides a real-time knowledge-sharing resource to address this need, covering experimental and computational aspects of sequencing and sequence analysis. Developers of popular analysis tools are among the >4000 active members, and ~40 peer-reviewed publications have referenced SEQanswers. AVAILABILITY: The SEQanswers community is freely accessible at http://SEQanswers.com/

Reference-independent comparative metagenomics using cross-assembly: crAss.

MOTIVATION: Metagenomes are often characterized by high levels of unknown sequences. Reads derived from known microorganisms can easily be identified and analyzed using fast homology search algorithms and a suitable reference database, but the unknown sequences are often ignored in further analyses, biasing conclusions. Nevertheless, it is possible to use more data in a comparative metagenomic analysis by creating a cross-assembly of all reads, i.e. a single assembly of reads from different samples. Comparative metagenomics studies the interrelationships between metagenomes from different samples. Using an assembly algorithm is a fast and intuitive way to link (partially) homologous reads without requiring a database of reference sequences. RESULTS: Here, we introduce crAss, a novel bioinformatic tool that enables fast simple analysis of cross-assembly files, yielding distances between all metagenomic sample pairs and an insightful image displaying the similarities.

PhiSiGns: an online tool to identify signature genes in phages and design PCR primers for examining phage diversity.

BACKGROUND: Phages (viruses that infect bacteria) have gained significant attention because of their abundance, diversity and important ecological roles. However, the lack of a universal gene shared by all phages presents a challenge for phage identification and characterization, especially in environmental samples where it is difficult to culture phage-host systems. Homologous conserved genes (or "signature genes") present in groups of closely-related phages can be used to explore phage diversity and define evolutionary relationships amongst these phages. Bioinformatic approaches are needed to identify candidate signature genes and design PCR primers to amplify those genes from environmental samples; however, there is currently no existing computational tool that biologists can use for this purpose. RESULTS: Here we present PhiSiGns, a web-based and standalone application that performs a pairwise comparison of each gene present in user-selected phage genomes, identifies signature genes, generates alignments of these genes, and designs potential PCR primer pairs. PhiSiGns is available at (http://www.phantome.org/phisigns/; http://phisigns.sourceforge.net/) with a link to the source code. Here we describe the specifications of PhiSiGns and demonstrate its application with a case study. CONCLUSIONS: PhiSiGns provides phage biologists with a user-friendly tool to identify signature genes and design PCR primers to amplify related genes from uncultured phages in environmental samples. This bioinformatics tool will facilitate the development of novel signature genes for use as molecular markers in studies of phage diversity, phylogeny, and evolution.

Case studies of the spatial heterogeneity of DNA viruses in the cystic fibrosis lung.

Microbial communities in the lungs of patients with cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) have been shown to be spatially heterogeneous. Viral communities may also vary spatially, leading to localized viral populations and infections. Here, we characterized viral communities from multiple areas of the lungs of two patients with late-stage CF using metagenomics, that is, the explanted lungs from a transplant patient and lungs acquired postmortem. All regions harbored eukaryotic viruses that may infect the human host, notably herpesviruses, anelloviruses, and papillomaviruses. In the highly diseased apical lobes of explant lungs, viral diversity was extremely low, and only eukaryotic viruses were present. The absence of phage suggests that CF-associated microbial biofilms may escape top-down controls by phage predation. The phages present in other lobes of explant lungs and in all lobes of postmortem lungs comprised distinct communities, and encoded genes for clinically important microbial phenotypes, including small colony variants and antibiotic resistance. Based on the these observations, we postulate that viral communities in CF lungs are spatially distinct and contribute to CF pathology by augmenting the metabolic potential of resident microbes, as well as by directly damaging lung tissue via carcinomas and herpesviral outbreaks.

Quality control and preprocessing of metagenomic datasets.

SUMMARY: Here, we present PRINSEQ for easy and rapid quality control and data preprocessing of genomic and metagenomic datasets. Summary statistics of FASTA (and QUAL) or FASTQ files are generated in tabular and graphical form and sequences can be filtered, reformatted and trimmed by a variety of options to improve downstream analysis. AVAILABILITY AND IMPLEMENTATION: This open-source application was implemented in Perl and can be used as a stand alone version or accessed online through a user-friendly web interface. The source code, user help and additional information are available at http://prinseq.sourceforge.net/.

Microfluidic PCR combined with pyrosequencing for identification of allelic variants with phenotypic associations among targeted Salmonella genes.

A novel targeted massive parallel sequencing approach identified genetic variation in eight known or predicted fimbrial adhesins for 46 Salmonella strains. The results highlight associations between specific adhesin alleles, host species, and antimicrobial resistance. The differentiation of allelic variants has potential applications for diagnostic microbiology and epidemiological investigations.

Coastal bacterioplankton community diversity along a latitudinal gradient in Latin America by means of V6 tag pyrosequencing.

The bacterioplankton diversity of coastal waters along a latitudinal gradient between Puerto Rico and Argentina was analyzed using a total of 134,197 high-quality sequences from the V6 hypervariable region of the small-subunit ribosomal RNA gene (16S rRNA) (mean length of 60 nt). Most of the OTUs were identified into Proteobacteria, Bacteriodetes, Cyanobacteria, and Actinobacteria, corresponding to approx. 80% of the total number of sequences. The number of OTUs corresponding to species varied between 937 and 1946 in the seven locations. Proteobacteria appeared at high frequency in the seven locations. An enrichment of Cyanobacteria was observed in Puerto Rico, whereas an enrichment of Bacteroidetes was detected in the Argentinian shelf and Uruguayan coastal lagoons. The highest number of sequences of Actinobacteria and Acidobacteria were obtained in the Amazon estuary mouth. The rarefaction curves and Good coverage estimator for species diversity suggested a significant coverage, with values ranging between 92 and 97% for Good coverage. Conserved taxa corresponded to aprox. 52% of all sequences. This study suggests that human-contaminated environments may influence bacterioplankton diversity.

TagCleaner: Identification and removal of tag sequences from genomic and metagenomic datasets.

BACKGROUND: Sequencing metagenomes that were pre-amplified with primer-based methods requires the removal of the additional tag sequences from the datasets. The sequenced reads can contain deletions or insertions due to sequencing limitations, and the primer sequence may contain ambiguous bases. Furthermore, the tag sequence may be unavailable or incorrectly reported. Because of the potential for downstream inaccuracies introduced by unwanted sequence contaminations, it is important to use reliable tools for pre-processing sequence data. RESULTS: TagCleaner is a web application developed to automatically identify and remove known or unknown tag sequences allowing insertions and deletions in the dataset. TagCleaner is designed to filter the trimmed reads for duplicates, short reads, and reads with high rates of ambiguous sequences. An additional screening for and splitting of fragment-to-fragment concatenations that gave rise to artificial concatenated sequences can increase the quality of the dataset. Users may modify the different filter parameters according to their own preferences. CONCLUSIONS: TagCleaner is a publicly available web application that is able to automatically detect and efficiently remove tag sequences from metagenomic datasets. It is easily configurable and provides a user-friendly interface. The interactive web interface facilitates export functionality for subsequent data processing, and is available at http://edwards.sdsu.edu/tagcleaner.

The GAAS metagenomic tool and its estimations of viral and microbial average genome size in four major biomes.

Metagenomic studies characterize both the composition and diversity of uncultured viral and microbial communities. BLAST-based comparisons have typically been used for such analyses; however, sampling biases, high percentages of unknown sequences, and the use of arbitrary thresholds to find significant similarities can decrease the accuracy and validity of estimates. Here, we present Genome relative Abundance and Average Size (GAAS), a complete software package that provides improved estimates of community composition and average genome length for metagenomes in both textual and graphical formats. GAAS implements a novel methodology to control for sampling bias via length normalization, to adjust for multiple BLAST similarities by similarity weighting, and to select significant similarities using relative alignment lengths. In benchmark tests, the GAAS method was robust to both high percentages of unknown sequences and to variations in metagenomic sequence read lengths. Re-analysis of the Sargasso Sea virome using GAAS indicated that standard methodologies for metagenomic analysis may dramatically underestimate the abundance and importance of organisms with small genomes in environmental systems. Using GAAS, we conducted a meta-analysis of microbial and viral average genome lengths in over 150 metagenomes from four biomes to determine whether genome lengths vary consistently between and within biomes, and between microbial and viral communities from the same environment. Significant differences between biomes and within aquatic sub-biomes (oceans, hypersaline systems, freshwater, and microbialites) suggested that average genome length is a fundamental property of environments driven by factors at the sub-biome level. The behavior of paired viral and microbial metagenomes from the same environment indicated that microbial and viral average genome sizes are independent of each other, but indicative of community responses to stressors and environmental conditions.

Corrigendum: Allelic variation contributes to bacterial host specificity.

This corrects the article DOI: 10.1038/ncomms9754.

Genomic Comparison of the Closely-Related Salmonella enterica Serovars Enteritidis, Dublin and Gallinarum.

The Salmonella enterica serovars Enteritidis, Dublin, and Gallinarum are closely related but differ in virulence and host range. To identify the genetic elements responsible for these differences and to better understand how these serovars are evolving, we sequenced the genomes of Enteritidis strain LK5 and Dublin strain SARB12 and compared these genomes to the publicly available Enteritidis P125109, Dublin CT 02021853 and Dublin SD3246 genome sequences. We also compared the publicly available Gallinarum genome sequences from biotype Gallinarum 287/91 and Pullorum RKS5078. Using bioinformatic approaches, we identified single nucleotide polymorphisms, insertions, deletions, and differences in prophage and pseudogene content between strains belonging to the same serovar. Through our analysis we also identified several prophage cargo genes and pseudogenes that affect virulence and may contribute to a host-specific, systemic lifestyle. These results strongly argue that the Enteritidis, Dublin and Gallinarum serovars of Salmonella enterica evolve by acquiring new genes through horizontal gene transfer, followed by the formation of pseudogenes. The loss of genes necessary for a gastrointestinal lifestyle ultimately leads to a systemic lifestyle and niche exclusion in the host-specific serovars.

Allelic variation contributes to bacterial host specificity.

Understanding the molecular parameters that regulate cross-species transmission and host adaptation of potential pathogens is crucial to control emerging infectious disease. Although microbial pathotype diversity is conventionally associated with gene gain or loss, the role of pathoadaptive nonsynonymous single-nucleotide polymorphisms (nsSNPs) has not been systematically evaluated. Here, our genome-wide analysis of core genes within Salmonella enterica serovar Typhimurium genomes reveals a high degree of allelic variation in surface-exposed molecules, including adhesins that promote host colonization. Subsequent multinomial logistic regression, MultiPhen and Random Forest analyses of known/suspected adhesins from 580 independent Typhimurium isolates identifies distinct host-specific nsSNP signatures. Moreover, population and functional analyses of host-associated nsSNPs for FimH, the type 1 fimbrial adhesin, highlights the role of key allelic residues in host-specific adherence in vitro. Together, our data provide the first concrete evidence that functional differences between allelic variants of bacterial proteins likely contribute to pathoadaption to diverse hosts.

Breath gas metabolites and bacterial metagenomes from cystic fibrosis airways indicate active pH neutral 2,3-butanedione fermentation.

The airways of cystic fibrosis (CF) patients are chronically colonized by patient-specific polymicrobial communities. The conditions and nutrients available in CF lungs affect the physiology and composition of the colonizing microbes. Recent work in bioreactors has shown that the fermentation product 2,3-butanediol mediates cross-feeding between some fermenting bacteria and Pseudomonas aeruginosa, and that this mechanism increases bacterial current production. To examine bacterial fermentation in the respiratory tract, breath gas metabolites were measured and several metagenomes were sequenced from CF and non-CF volunteers. 2,3-butanedione was produced in nearly all respiratory tracts. Elevated levels in one patient decreased during antibiotic treatment, and breath concentrations varied between CF patients at the same time point. Some patients had high enough levels of 2,3-butanedione to irreversibly damage lung tissue. Antibiotic therapy likely dictates the activities of 2,3-butanedione-producing microbes, which suggests a need for further study with larger sample size. Sputum microbiomes were dominated by P. aeruginosa, Streptococcus spp. and Rothia mucilaginosa, and revealed the potential for 2,3-butanedione biosynthesis. Genes encoding 2,3-butanedione biosynthesis were disproportionately abundant in Streptococcus spp, whereas genes for consumption of butanedione pathway products were encoded by P. aeruginosa and R. mucilaginosa. We propose a model where low oxygen conditions in CF lung lead to fermentation and a decrease in pH, triggering 2,3-butanedione fermentation to avoid lethal acidification. We hypothesize that this may also increase phenazine production by P. aeruginosa, increasing reactive oxygen species and providing additional electron acceptors to CF microbes.

Species-specific viromes in the ancestral holobiont Hydra.

Recent evidence showing host specificity of colonizing bacteria supports the view that multicellular organisms are holobionts comprised of the macroscopic host in synergistic interdependence with a heterogeneous and host-specific microbial community. Whereas host-bacteria interactions have been extensively investigated, comparatively little is known about host-virus interactions and viral contribution to the holobiont. We sought to determine the viral communities associating with different Hydra species, whether these viral communities were altered with environmental stress, and whether these viruses affect the Hydra-associated holobiont. Here we show that each species of Hydra harbors a diverse host-associated virome. Primary viral families associated with Hydra are Myoviridae, Siphoviridae, Inoviridae, and Herpesviridae. Most Hydra-associated viruses are bacteriophages, a reflection of their involvement in the holobiont. Changes in environmental conditions alter the associated virome, increase viral diversity, and affect the metabolism of the holobiont. The specificity and dynamics of the virome point to potential viral involvement in regulating microbial associations in the Hydra holobiont. While viruses are generally regarded as pathogenic agents, our study suggests an evolutionary conserved ability of viruses to function as holobiont regulators and, therefore, constitutes an emerging paradigm shift in host-microbe interactions.

Characterization of Amoeboaphelidium protococcarum, an algal parasite new to the cryptomycota isolated from an outdoor algal pond used for the production of biofuel.

Mass culture of algae for the production of biofuels is a developing technology designed to offset the depletion of fossil fuel reserves. However, large scale culture of algae in open ponds can be challenging because of incidences of infestation with algal parasites. Without knowledge of the identity of the specific parasite and how to control these pests, algal-based biofuel production will be limited. We have characterized a eukaryotic parasite of Scenedesmus dimorphus growing in outdoor ponds used for biofuel production. We demonstrated that as the genomic DNA of parasite FD01 increases, the concentration of S. dimorphus cells decreases; consequently, this is a highly destructive pathogen. Techniques for culture of the parasite and host were developed, and the endoparasite was identified as the Aphelidea, Amoeboaphelidium protococcarum. Phylogenetic analysis of ribosomal sequences revealed that parasite FD01 placed within the recently described Cryptomycota, a poorly known phylum based on two species of Rozella and environmental samples. Transmission electron microscopy demonstrated that aplanospores of the parasite produced filose pseudopodia, which contained fine fibers the diameter of actin microfilaments. Multiple lipid globules clustered and were associated with microbodies, mitochondria and a membrane cisternae, an arrangement characteristic of the microbody-lipid globule complex of chytrid zoospores. After encystment and attachment to the host cells, the parasite injected its protoplast into the host between the host cell wall and plasma membrane. At maturity the unwalled parasite occupied the entire host cell. After cleavage of the protoplast into aplanospores, a vacuole and lipids remained in the host cell. Amoeboaphelidium protococcarum isolate FD01 is characteristic of the original description of this species and is different from strain X-5 recently characterized. Our results help put a face on the Cryptomycota, revealing that the phylum is more diverse than previously understood and include some of the Aphelidea as well as Rozella species and potentially Microsporidia.