Corrigendum: Allelic variation contributes to bacterial host specificity.

This corrects the article DOI: 10.1038/ncomms9754.

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.

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.

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.

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.

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.

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.

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.

Combining de novo and reference-guided assembly with scaffold_builder.

Genome sequencing has become routine, however genome assembly still remains a challenge despite the computational advances in the last decade. In particular, the abundance of repeat elements in genomes makes it difficult to assemble them into a single complete sequence. Identical repeats shorter than the average read length can generally be assembled without issue. However, longer repeats such as ribosomal RNA operons cannot be accurately assembled using existing tools. The application Scaffold_builder was designed to generate scaffolds - super contigs of sequences joined by N-bases - based on the similarity to a closely related reference sequence. This is independent of mate-pair information and can be used complementarily for genome assembly, e.g. when mate-pairs are not available or have already been exploited. Scaffold_builder was evaluated using simulated pyrosequencing reads of the bacterial genomes Escherichia coli 042, Lactobacillus salivarius UCC118 and Salmonella enterica subsp. enterica serovar Typhi str. P-stx-12. Moreover, we sequenced two genomes from Salmonella enterica serovar Typhimurium LT2 G455 and Salmonella enterica serovar Typhimurium SDT1291 and show that Scaffold_builder decreases the number of contig sequences by 53% while more than doubling their average length. Scaffold_builder is written in Python and is available at http://edwards.sdsu.edu/scaffold_builder. A web-based implementation is additionally provided to allow users to submit a reference genome and a set of contigs to be scaffolded.

Mechanistic model of Rothia mucilaginosa adaptation toward persistence in the CF lung, based on a genome reconstructed from metagenomic data.

The impaired mucociliary clearance in individuals with Cystic Fibrosis (CF) enables opportunistic pathogens to colonize CF lungs. Here we show that Rothia mucilaginosa is a common CF opportunist that was present in 83% of our patient cohort, almost as prevalent as Pseudomonas aeruginosa (89%). Sequencing of lung microbial metagenomes identified unique R. mucilaginosa strains in each patient, presumably due to evolution within the lung. The de novo assembly of a near-complete R. mucilaginosa (CF1E) genome illuminated a number of potential physiological adaptations to the CF lung, including antibiotic resistance, utilization of extracellular lactate, and modification of the type I restriction-modification system. Metabolic characteristics predicted from the metagenomes suggested R. mucilaginosa have adapted to live within the microaerophilic surface of the mucus layer in CF lungs. The results also highlight the remarkable evolutionary and ecological similarities of many CF pathogens; further examination of these similarities has the potential to guide patient care and treatment.

Multivariate analysis of functional metagenomes.

Metagenomics is a primary tool for the description of microbial and viral communities. The sheer magnitude of the data generated in each metagenome makes identifying key differences in the function and taxonomy between communities difficult to elucidate. Here we discuss the application of seven different data mining and statistical analyses by comparing and contrasting the metabolic functions of 212 microbial metagenomes within and between 10 environments. Not all approaches are appropriate for all questions, and researchers should decide which approach addresses their questions. This work demonstrated the use of each approach: for example, random forests provided a robust and enlightening description of both the clustering of metagenomes and the metabolic processes that were important in separating microbial communities from different environments. All analyses identified that the presence of phage genes within the microbial community was a predictor of whether the microbial community was host-associated or free-living. Several analyses identified the subtle differences that occur with environments, such as those seen in different regions of the marine environment.

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.

Metagenomics and metatranscriptomics: windows on CF-associated viral and microbial communities.

BACKGROUND: Samples collected from CF patient airways often contain large amounts of host-derived nucleic acids that interfere with recovery and purification of microbial and viral nucleic acids. This study describes metagenomic and metatranscriptomic methods that address these issues. METHODS: Microbial and viral metagenomes, and microbial metatranscriptomes, were successfully prepared from sputum samples from five adult CF patients. RESULTS: Contaminating host DNA was dramatically reduced in the metagenomes. Each CF patient presented a unique microbiome; in some Pseudomonas aeruginosa was replaced by other opportunistic bacteria. Even though the taxonomic composition of the microbiomes is very different, the metabolic potentials encoded by the community are very similar. The viral communities were dominated by phages that infect major CF pathogens. The metatranscriptomes reveal differential expression of encoded metabolic potential with changing health status. CONCLUSIONS: Microbial and viral metagenomics combined with microbial transcriptomics characterize the dynamic polymicrobial communities found in CF airways, revealing both the taxa present and their current metabolic activities. These approaches can facilitate the development of individualized treatment plans and novel therapeutic approaches.

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.

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.

Abrolhos bank reef health evaluated by means of water quality, microbial diversity, benthic cover, and fish biomass data.

The health of the coral reefs of the Abrolhos Bank (Southwestern Atlantic) was characterized with a holistic approach using measurements of four ecosystem components: (i) inorganic and organic nutrient concentrations, [1] fish biomass, [1] macroalgal and coral cover and (iv) microbial community composition and abundance. The possible benefits of protection from fishing were particularly evaluated by comparing sites with varying levels of protection. Two reefs within the well-enforced no-take area of the National Marine Park of Abrolhos (Parcel dos Abrolhos and California) were compared with two unprotected coastal reefs (Sebastião Gomes and Pedra de Leste) and one legally protected but poorly enforced coastal reef (the "paper park" of Timbebas Reef). The fish biomass was lower and the fleshy macroalgal cover was higher in the unprotected reefs compared with the protected areas. The unprotected and protected reefs had similar seawater chemistry. Lower vibrio CFU counts were observed in the fully protected area of California Reef. Metagenome analysis showed that the unprotected reefs had a higher abundance of archaeal and viral sequences and more bacterial pathogens, while the protected reefs had a higher abundance of genes related to photosynthesis. Similar to other reef systems in the world, there was evidence that reductions in the biomass of herbivorous fishes and the consequent increase in macroalgal cover in the Abrolhos Bank may be affecting microbial diversity and abundance. Through the integration of different types of ecological data, the present study showed that protection from fishing may lead to greater reef health. The data presented herein suggest that protected coral reefs have higher microbial diversity, with the most degraded reef (Sebastião Gomes) showing a marked reduction in microbial species richness. It is concluded that ecological conditions in unprotected reefs may promote the growth and rapid evolution of opportunistic microbial pathogens.

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/

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.

Spatial distribution of microbial communities in the cystic fibrosis lung.

Cystic fibrosis (CF) is a common fatal genetic disorder with mortality most often resulting from microbial infections of the lungs. Culture-independent studies of CF-associated microbial communities have indicated that microbial diversity in the CF airways is much higher than suggested by culturing alone. However, these studies have relied on indirect methods to sample the CF lung such as expectorated sputum and bronchoalveolar lavage (BAL). Here, we characterize the diversity of microbial communities in tissue sections from anatomically distinct regions of the CF lung using barcoded 16S amplicon pyrosequencing. Microbial communities differed significantly between different areas of the lungs, and few taxa were common to microbial communities in all anatomical regions surveyed. Our results indicate that CF lung infections are not only polymicrobial, but also spatially heterogeneous suggesting that treatment regimes tailored to dominant populations in sputum or BAL samples may be ineffective against infections in some areas of the lung.