FIGfams: yet another set of protein families

Taxonomic distribution of metabolic functions in bacteria associated with Trichodesmium consortia

IMPORTANCE

Colonies of the cyanobacteria Trichodesmium act as a biological hotspot for the usage and recycling of key resources such as C, N, P, and Fe within an otherwise oligotrophic environment. While Trichodesmium colonies are known to interact and support a unique community of algae and particle-associated microbes, our understanding of the taxa that populate these colonies and the gene functions they encode is still limited. Characterizing the taxa and adaptive strategies that influence consortium physiology and its concomitant biogeochemistry is critical in a future ocean predicted to have increasingly resource-depleted regions.

Analyzing Prokaryotic Transcriptomics in the Light of Genome Data with the MicroScope Platform

Large-scale genome sequencing and the increasingly massive use of high-throughput approaches produce a vast amount of new information that completely transforms our understanding of thousands of microbial species occurring in our environment. However, despite the development of powerful bioinformatics approaches, full interpretation of the content of these genomes remains a difficult task. To address this challenge, the MicroScope platform has been developed. It is an integrated Web platform for management, annotation, comparative analysis, and visualization of microbial genomes ( https://mage.genoscope.cns.fr/microscope ). Launched in 2005, the platform has been under continuous development and provides analyzes for complete and ongoing genome projects together with metabolic network reconstruction and transcriptomic experiments allowing users to improve the understanding of gene functions. MicroScope platform is widely used by microbiologists from academia and industry all around the world for collaborative studies and expert annotation. It enables collaborative work in a rich comparative genomic context and improves community-based curation efforts. Here, we describe the protocol to follow for the integration and analysis of transcriptomics data in the Microscope platform. The chapter reviews each key step from the experimental design to the analysis and interpretation of the experiment data and results. The integration of transcriptomics data gives a dynamic view of the genome by allowing the users to improve the understanding of gene functions by interpreting them in the light of regulatory cell processes. Moreover, they can also contribute to the refinement of genome annotation through the discovery of new genes and help to fill metabolic gaps.

A shotgun approach to explore the bacterial diversity and a brief insight into the glycoside hydrolases of Samiti lake located in the Eastern Himalayas

BACKGROUND

The Himalayas have always been an enigma and, being biodiversity hotspots, are considered extremely important from an ecological point of view. Recent advances in studies regarding high-altitude lakes have garnered relevant importance as these habitats could harbor potential psychrophilic and psychrotrophic microbes with bio-prospective applications. Contemplating the above scenario, the present study has been undertaken to understand the diversity and the functional capacities of the microbes thriving in this lake.

RESULTS

In our present study on Samiti Lake, the abundance of Proteobacteria as the major phylum was seen in both the soil and water samples. Incase of the ABSLW (water) and ABS1 (soil) sample, 148,066 and 239,754 predicted genes, were taken for functional analysis. The KEGG analysis showed that ABSLW and ABS1 had 122,911 and 160,268, genes assigned to KO terms respectively. Whereas in case of COG functional analysis, 104,334 and 130,191 genes were assigned to different COG classes for ABSLW and ABS1 respectively. Further, on studying the glycoside hydrolases, an abundance of GH13, GH2, GH3, GH43, and GH23 in both the soil and water samples were seen.

CONCLUSION

Our study has provided a comprehensive report about the bacterial diversity and functional capacities of microbes thriving in Samiti Lake.  It has also thrown some light on the occurrence of glycoside hydrolases in this region, as they have numerous biotechnological applications in different sectors.

Integrative functional analysis uncovers metabolic differences between Candida species

Candida species are a dominant constituent of the human mycobiome and associated with the development of several diseases. Understanding the Candida species metabolism could provide key insights into their ability to cause pathogenesis. Here, we have developed the BioFung database, providing an efficient annotation of protein-encoding genes. Along, with BioFung, using carbohydrate-active enzyme (CAZymes) analysis, we have uncovered core and accessory features across Candida species demonstrating plasticity, adaption to the environment and acquired features. We show a greater importance of amino acid metabolism, as functional analysis revealed that all Candida species can employ amino acid metabolism. However, metabolomics revealed that only a specific cluster of species (AGAu species—C. albicans, C. glabrata and C. auris) utilised amino acid metabolism including arginine, cysteine, and methionine metabolism potentially improving their competitive fitness in pathogenesis. We further identified critical metabolic pathways in the AGAu cluster with biomarkers and anti-fungal target potential in the CAZyme profile, polyamine, choline and fatty acid biosynthesis pathways. This study, combining genomic analysis, and validation with gene expression and metabolomics, highlights the metabolic diversity with AGAu species that underlies their remarkable ability to dominate they mycobiome and cause disease.

Metabolic differences between Candida species are uncovered using the BioFung database alongside genomic and metabolic analysis.

New Multidrug Efflux Systems in a Microcystin-Degrading Bacterium Blastomonas fulva and Its Genomic Feature

A microcystin-degrading bacterial strain, Blastomonas fulva T2, was isolated from the culture of a microalgae Microcystis. The strain B. fulva T2 is Gram-stain-negative, non-motile, aerobic, non-spore-forming and phototrophic. The cells of B. fulva T2 are able to grow in ranges of temperature from 15 to 37 °C, with a pH of 6 to 8 and a salinity of 0 to 1% NaCl. Here, we sequenced the complete genome of B. fulva T2, aiming to better understand the evolutionary biology and the function of the genus Blastomonas at the molecular level. The complete genome of B. fulva T2 contained a circular chromosome (3,977,381 bp) with 64.3% GC content and a sizable plasmid (145.829 bp) with 60.7% GC content which comprises about 3.5% of the total genetic content. A total of 3842 coding genes, including 46 tRNAs and 6 rRNAs, were predicted in the genome. The genome contains genes for glycolysis, citric acid cycle, Entner–Doudoroff pathways, photoreaction center and bacteriochlorophylla synthesis. A 7.9 K gene cluster containing mlrA, mlrB, mlrC and mlrD_1,2,3,4 of microcystin-degrading enzymes was identified. Notably, eight different efflux pumps categorized into RND, ABC and MFS types have been identified in the genome of strain T2. Our findings should provide new insights of the alternative reaction pathway as well as the enzymes which mediated the degradation of microcystin by bacteria, as well as the evolution, architectures, chemical mechanisms and physiological roles of the new bacterial multidrug efflux system.

MCRL: using a reference library to compress a metagenome into a non-redundant list of sequences, considering viruses as a case study

MOTIVATION

Metagenomes offer a glimpse into the total genomic diversity contained within a sample. Currently, however, there is no straightforward way to obtain a non-redundant list of all putative homologs of a set of reference sequences present in a metagenome.

RESULTS

To address this problem, we developed a novel clustering approach called 'metagenomic clustering by reference library' (MCRL), where a reference library containing a set of reference genes is clustered with respect to an assembled metagenome. According to our proposed approach, reference genes homologous to similar sets of metagenomic sequences, termed 'signatures', are iteratively clustered in a greedy fashion, retaining at each step the reference genes yielding the lowest E values, and terminating when signatures of remaining reference genes have a minimal overlap. The outcome of this computation is a non-redundant list of reference genes homologous to minimally overlapping sets of contigs, representing potential candidates for gene families present in the metagenome. Unlike metagenomic clustering methods, there is no need for contigs to overlap to be associated with a cluster, enabling MCRL to draw on more information encoded in the metagenome when computing tentative gene families. We demonstrate how MCRL can be used to extract candidate viral gene families from an oral metagenome and an oral virome that otherwise could not be determined using standard approaches. We evaluate the sensitivity, accuracy and robustness of our proposed method for the viral case study and compare it with existing analysis approaches.

AVAILABILITY AND IMPLEMENTATION

https://github.com/a-tadmor/MCRL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Accessory Genome Dynamics and Structural Variation of Shigella from Persistent Infections

Shigellosis is a diarrheal disease caused mainly by Shigella flexneri and Shigella sonnei Infection is thought to be largely self-limiting, with short- to medium-term and serotype-specific immunity provided following clearance. However, cases of men who have sex with men (MSM)-associated shigellosis have been reported where Shigella of the same serotype were serially sampled from individuals between 1 and 1,862 days apart, possibly due to persistent carriage or reinfection with the same serotype. Here, we investigate the accessory genome dynamics of MSM-associated S. flexneri and S. sonnei isolates serially sampled from individual patients at various days apart to shed light on the adaptation of these important pathogens during infection. We find that pairs likely associated with persistent infection/carriage and with a smaller single nucleotide polymorphism (SNP) distance, demonstrated significantly less variation in accessory genome content than pairs likely associated with reinfection, and with a greater SNP distance. We observed antimicrobial resistance acquisition during Shigella carriage, including the gain of an extended-spectrum beta-lactamase gene during carriage. Finally, we explored large chromosomal structural variations and rearrangements in seven (five chronic and two reinfection associated) pairs of S. flexneri 3a isolates from an MSM-associated epidemic sublineage, which revealed variations at several common regions across isolate pairs, mediated by insertion sequence elements and comprising a distinct predicted functional profile. This study provides insight on the variation of accessory genome dynamics and large structural genomic changes in Shigella during persistent infection/carriage. In addition, we have also created a complete reference genome and biobanked isolate of the globally important pathogen, S. flexneri 3a.IMPORTANCEShigella spp. are Gram-negative bacteria that are the etiological agent of shigellosis, the second most common cause of diarrheal illness among children under the age of five in low-income countries. In high-income countries, shigellosis is also a sexually transmissible disease among men who have sex with men. Within the latter setting, we have captured prolonged and/or recurrent infection with shigellae of the same serotype, challenging the belief that Shigella infection is short lived and providing an early opportunity to study the evolution of the pathogen over the course of infection. Using this recently emerged transmission scenario, we comprehensively characterize the genomic changes that occur over the course of individual infection with Shigella and uncover a distinct functional profile of variable genomic regions, findings that have relevance for other Enterobacteriaceae.

Three Novel Clostridia Isolates Produce n-Caproate and iso-Butyrate from Lactate: Comparative Genomics of Chain-Elongating Bacteria

The platform chemicals n-caproate and iso-butyrate can be produced by anaerobic fermentation from agro-industrial residues in a process known as microbial chain elongation. Few lactate-consuming chain-elongating species have been isolated and knowledge on their shared genetic features is still limited. Recently we isolated three novel clostridial strains (BL-3, BL-4, and BL-6) that convert lactate to n-caproate and iso-butyrate. Here, we analyzed the genetic background of lactate-based chain elongation in these isolates and other chain-elongating species by comparative genomics. The three strains produced n-caproate, n-butyrate, iso-butyrate, and acetate from lactate, with the highest proportions of n-caproate (18%) for BL-6 and of iso-butyrate (23%) for BL-4 in batch cultivation at pH 5.5. They show high genomic heterogeneity and a relatively small core-genome size. The genomes contain highly conserved genes involved in lactate oxidation, reverse β-oxidation, hydrogen formation and either of two types of energy conservation systems (Rnf and Ech). Including genomes of another eleven experimentally validated chain-elongating strains, we found that the chain elongation-specific core-genome encodes the pathways for reverse β-oxidation, hydrogen formation and energy conservation, while displaying substantial genome heterogeneity. Metabolic features of these isolates are important for biotechnological applications in n-caproate and iso-butyrate production.

An assessment of genome annotation coverage across the bacterial tree of life

Although gene-finding in bacterial genomes is relatively straightforward, the automated assignment of gene function is still challenging, resulting in a vast quantity of hypothetical sequences of unknown function. But how prevalent are hypothetical sequences across bacteria, what proportion of genes in different bacterial genomes remain unannotated, and what factors affect annotation completeness? To address these questions, we surveyed over 27 000 bacterial genomes from the Genome Taxonomy Database, and measured genome annotation completeness as a function of annotation method, taxonomy, genome size, 'research bias' and publication date. Our analysis revealed that 52 and 79 % of the average bacterial proteome could be functionally annotated based on protein and domain-based homology searches, respectively. Annotation coverage using protein homology search varied significantly from as low as 14 % in some species to as high as 98 % in others. We found that taxonomy is a major factor influencing annotation completeness, with distinct trends observed across the microbial tree (e.g. the lowest level of completeness was found in the Patescibacteria lineage). Most lineages showed a significant association between genome size and annotation incompleteness, likely reflecting a greater degree of uncharacterized sequences in 'accessory' proteomes than in 'core' proteomes. Finally, research bias, as measured by publication volume, was also an important factor influencing genome annotation completeness, with early model organisms showing high completeness levels relative to other genomes in their own taxonomic lineages. Our work highlights the disparity in annotation coverage across the bacterial tree of life and emphasizes a need for more experimental characterization of accessory proteomes as well as understudied lineages.

Genetics and regulation of nitrogen fixation in Paenibacillus brasilensis PB24

Biological nitrogen fixation (BNF), performed by diazotrophic prokaryotes, is responsible for reducing dinitrogen (N₂) present in the biosphere into biologically available forms of nitrogen. Paenibacillus brasilensis PB24 is a diazotrophic Gram-positive bacterium and is considered ecologically and industrially important because it is able to produce antimicrobial substances and 2,3-butanediol. However, the genetics and regulation of its nitrogen fixing (nif) genes have never been assessed so far. Therefore, the present study aimed to (i) identify the structural and regulatory genes related to BNF in the PB24 genome, (ii) perform comparative genomics analysis of the nif operon among different Paenibacillus species and (iii) study the expression of these genes in the presence and absence of NH₄. Strain PB24 showed a nif operon composed of nine genes (nifBHDKENXhesAV), with a conserved synteny (with small variations) among the Paenibacillus species evaluated. BNF regulatory genes, glnK and amtB (encoding GlnK signal transduction protein and AmtB transmembrane protein, respectively) and glnR and glnA genes (encoding the transcription factor GlnR and glutamine synthetase) were found in the PB24 genome. Primers were designed for qPCR amplification of the nitrogenase structural (nifH, nifD and nifK) and regulatory (glnA and amtB) BNF genes. The structural gene expression in PB24 was up- and downregulated in the absence and presence of NH₄, respectively. The gene expression levels indicated a GlnR-mediated repression of genes associated with ammonium import (amtBglnK) and BNF (nif genes). Additionally, the regulatory mechanism of GlnR in P. brasilensis PB24 differed from the other Paenibacillus evaluated, considering the different distribution of binding sites recognized by GlnR.

Construction of a complete set of Neisseria meningitidis mutants and its use for the phenotypic profiling of this human pathogen

The bacterium Neisseria meningitidis causes life-threatening meningitis and sepsis. Here, we construct a complete collection of defined mutants in protein-coding genes of this organism, identifying all genes that are essential under laboratory conditions. The collection, named NeMeSys 2.0, consists of individual mutants in 1584 non-essential genes. We identify 391 essential genes, which are associated with basic functions such as expression and preservation of genome information, cell membrane structure and function, and metabolism. We use this collection to shed light on the functions of diverse genes, including a gene encoding a member of a previously unrecognised class of histidinol-phosphatases; a set of 20 genes required for type IV pili function; and several conditionally essential genes encoding antitoxins and/or immunity proteins. We expect that NeMeSys 2.0 will facilitate the phenotypic profiling of a major human bacterial pathogen.

Integrated Metabolic Modeling, Culturing, and Transcriptomics Explain Enhanced Virulence of Vibrio cholerae during Coinfection with Enterotoxigenic Escherichia coli

Most studies proposing new strategies to manage and treat infections have been largely focused on identifying druggable targets that can inhibit a pathogen's growth when it is the single cause of infection. In vivo, however, infections can be caused by multiple species. This is important to take into account when attempting to develop or use current antibacterials since their efficacy can change significantly between single infections and coinfections. In this study, we used genome-scale metabolic models (GEMs) to interrogate the growth capabilities of Vibrio cholerae in single infections and coinfections with enterotoxigenic E. coli (ETEC), which cooccur in a large fraction of diarrheagenic patients. Coinfection model predictions showed that V. cholerae growth capabilities are enhanced in the presence of ETEC relative to V. cholerae single infection, through cross-fed metabolites made available to V. cholerae by ETEC. In vitro, cocultures of the two enteric pathogens further confirmed model predictions showing an increased growth of V. cholerae in coculture relative to V. cholerae single cultures while ETEC growth was suppressed. Dual RNAseq analysis of the cocultures also confirmed that the transcriptome of V. cholerae was distinct during coinfection compared to single-infection scenarios where processes related to metabolism were significantly perturbed. Further, in silico gene-knockout simulations uncovered discrepancies in gene essentiality for V. cholerae growth between single infections and coinfections. Integrative model-guided analysis thus identified druggable targets that would be critical for V. cholerae growth in both single infections and coinfections; thus, designing inhibitors against those targets would provide a broader spectrum of coverage against cholera infections.

Gene essentiality is altered during polymicrobial infections. Nevertheless, most studies rely on single-species infections to assess pathogen gene essentiality. Here, we use genome-scale metabolic models (GEMs) to explore the effect of coinfection of the diarrheagenic pathogen Vibrio cholerae with another enteric pathogen, enterotoxigenic Escherichia coli (ETEC). Model predictions showed that V. cholerae metabolic capabilities were increased due to ample cross-feeding opportunities enabled by ETEC. This is in line with increased severity of cholera symptoms known to occur in patients with dual infections by the two pathogens. In vitro coculture systems confirmed that V. cholerae growth is enhanced in cocultures relative to single cultures. Further, expression levels of several V. cholerae metabolic genes were significantly perturbed as shown by dual RNA sequencing (RNAseq) analysis of its cocultures with different ETEC strains. A decrease in ETEC growth was also observed, probably mediated by nonmetabolic factors. Single gene essentiality analysis predicted conditionally independent genes that are essential for the pathogen’s growth in both single-infection and coinfection scenarios. Our results reveal growth differences that are of relevance to drug targeting and efficiency in polymicrobial infections.

IMPORTANCE Most studies proposing new strategies to manage and treat infections have been largely focused on identifying druggable targets that can inhibit a pathogen's growth when it is the single cause of infection. In vivo, however, infections can be caused by multiple species. This is important to take into account when attempting to develop or use current antibacterials since their efficacy can change significantly between single infections and coinfections. In this study, we used genome-scale metabolic models (GEMs) to interrogate the growth capabilities of Vibrio cholerae in single infections and coinfections with enterotoxigenic E. coli (ETEC), which cooccur in a large fraction of diarrheagenic patients. Coinfection model predictions showed that V. cholerae growth capabilities are enhanced in the presence of ETEC relative to V. cholerae single infection, through cross-fed metabolites made available to V. cholerae by ETEC. In vitro, cocultures of the two enteric pathogens further confirmed model predictions showing an increased growth of V. cholerae in coculture relative to V. cholerae single cultures while ETEC growth was suppressed. Dual RNAseq analysis of the cocultures also confirmed that the transcriptome of V. cholerae was distinct during coinfection compared to single-infection scenarios where processes related to metabolism were significantly perturbed. Further, in silico gene-knockout simulations uncovered discrepancies in gene essentiality for V. cholerae growth between single infections and coinfections. Integrative model-guided analysis thus identified druggable targets that would be critical for V. cholerae growth in both single infections and coinfections; thus, designing inhibitors against those targets would provide a broader spectrum of coverage against cholera infections.

Comparative genomics of the genus Roseburia reveals divergent biosynthetic pathways that may influence colonic competition among species

Roseburia species are important denizens of the human gut microbiome that ferment complex polysaccharides to butyrate as a terminal fermentation product, which influences human physiology and serves as an energy source for colonocytes. Previous comparative genomics analyses of the genus Roseburia have examined polysaccharide degradation genes. Here, we characterize the core and pangenomes of the genus Roseburia with respect to central carbon and energy metabolism, as well as biosynthesis of amino acids and B vitamins using orthology-based methods, uncovering significant differences among species in their biosynthetic capacities. Variation in gene content among Roseburia species and strains was most significant for cofactor biosynthesis. Unlike all other species of Roseburia that we analysed, Roseburia inulinivorans strains lacked biosynthetic genes for riboflavin or pantothenate but possessed folate biosynthesis genes. Differences in gene content for B vitamin synthesis were matched with differences in putative salvage and synthesis strategies among species. For example, we observed extended biotin salvage capabilities in R. intestinalis strains, which further suggest that B vitamin acquisition strategies may impact fitness in the gut ecosystem. As differences in the functional potential to synthesize components of biomass (e.g. amino acids, vitamins) can drive interspecies interactions, variation in auxotrophies of the Roseburia spp. genomes may influence in vivo gut ecology. This study serves to advance our understanding of the potential metabolic interactions that influence the ecology of Roseburia spp. and, ultimately, may provide a basis for rational strategies to manipulate the abundances of these species.

Similar genomic patterns of clinical infective endocarditis and oral isolates of Streptococcus sanguinis and Streptococcus gordonii

Streptococcus gordonii and Streptococcus sanguinis belong to the Mitis group streptococci, which mostly are commensals in the human oral cavity. Though they are oral commensals, they can escape their niche and cause infective endocarditis, a severe infection with high mortality. Several virulence factors important for the development of infective endocarditis have been described in these two species. However, the background for how the commensal bacteria, in some cases, become pathogenic is still not known. To gain a greater understanding of the mechanisms of the pathogenic potential, we performed a comparative analysis of 38 blood culture strains, S. sanguinis (n = 20) and S. gordonii (n = 18) from patients with verified infective endocarditis, along with 21 publicly available oral isolates from healthy individuals, S. sanguinis (n = 12) and S. gordonii (n = 9). Using whole genome sequencing data of the 59 streptococci genomes, functional profiles were constructed, using protein domain predictions based on the translated genes. These functional profiles were used for clustering, phylogenetics and machine learning. A clear separation could be made between the two species. No clear differences between oral isolates and clinical infective endocarditis isolates were found in any of the 675 translated core-genes. Additionally, random forest-based machine learning and clustering of the pan-genome data as well as amino acid variations in the core-genome could not separate the clinical and oral isolates. A total of 151 different virulence genes was identified in the 59 genomes. Among these homologs of genes important for adhesion and evasion of the immune system were found in all of the strains. Based on the functional profiles and virulence gene content of the genomes, we believe that all analysed strains had the ability to become pathogenic.

Greedy de novo motif discovery to construct motif repositories for bacterial proteomes

BACKGROUND

Bacterial surfaces are complex systems, constructed from membranes, peptidoglycan and, importantly, proteins. The proteins play crucial roles as critical regulators of how the bacterium interacts with and survive in its environment. A full catalog of the motifs in protein families and their relative conservation grade is a prerequisite to target the protein-protein interaction that bacterial surface protein makes to host proteins.

RESULTS

In this paper, we propose a greedy approach to identify conserved motifs in large sequence families iteratively. Each iteration discovers a motif de novo and masks all occurrences of that motif. Remaining unmasked sequences are subjected to the next round of motif detection until no more significant motifs can be found. We demonstrate the utility of the method through the construction of a proteome-wide motif repository for Group A Streptococcus (GAS), a significant human pathogen. GAS produce numerous surface proteins that interact with over 100 human plasma proteins, helping the bacteria to evade the host immune response. We used the repository to find that proteins part of the bacterial surface has motif architectures that differ from intracellular proteins.

CONCLUSIONS

We elucidate that the M protein, a coiled-coil homodimer that extends over 500 A from the cell wall, has a motif architecture that differs between various GAS strains. As the M protein is known to bind a variety of different plasma proteins, the results indicate that the different motif architectures are responsible for the quantitative differences of plasma proteins that various strains bind. The speed and applicability of the method enable its application to all major human pathogens.

Genome properties in 2019: a new companion database to InterPro for the inference of complete functional attributes

Automatic annotation of protein function is routinely applied to newly sequenced genomes. While this provides a fine-grained view of an organism's functional protein repertoire, proteins, more commonly function in a coordinated manner, such as in pathways or multimeric complexes. Genome Properties (GPs) define such functional entities as a series of steps, originally described by either TIGRFAMs or Pfam entries. To increase the scope of coverage, we have migrated GPs to function as a companion resource utilizing InterPro entries. Having introduced GPs-specific versioned releases, we provide software and data via a GitHub repository, and have developed a new web interface to GPs (available at https://www.ebi.ac.uk/interpro/genomeproperties). In addition to exploring each of the 1286 GPs, the website contains GPs pre-calculated for a representative set of proteomes; these results can be used to profile GPs phylogenetically via an interactive viewer. Users can upload novel data to the viewer for comparison with the pre-calculated results. Over the last year, we have added ∼700 new GPs, increasing the coverage of eukaryotic systems, as well as increasing general coverage through automatic generation of GPs from related resources. All data are freely available via the website and the GitHub repository.

Metabolic specialization of denitrifiers in permeable sediments controls N2 O emissions

Coastal oceans receive large amounts of anthropogenic fixed nitrogen (N), most of which is denitrified in the sediment before reaching the open ocean. Sandy sediments, which are common in coastal regions, seem to play an important role in catalysing this N-loss. Permeable sediments are characterized by advective porewater transport, which supplies high fluxes of organic matter into the sediment, but also leads to fluctuations in oxygen and nitrate concentrations. Little is known about how the denitrifying communities in these sediments are adapted to such fluctuations. Our combined results indicate that denitrification in eutrophied sandy sediments from the world's largest tidal flat system, the Wadden Sea, is carried out by different groups of microorganisms. This segregation leads to the formation of N₂ O which is advectively transported to the overlying waters and thereby emitted to the atmosphere. At the same time, the production of N₂ O within the sediment supports a subset of Flavobacteriia which appear to be specialized on N₂ O reduction. If the mechanisms shown here are active in other coastal zones, then denitrification in eutrophied sandy sediments may substantially contribute to current marine N₂ O emissions.

The Complete Genome and Physiological Analysis of the Microbialite-Dwelling Agrococcus pavilionensis sp. nov; Reveals Genetic Promiscuity and Predicted Adaptations to Environmental Stress

Members of the bacterial genus Agrococcus are globally distributed and found across environments so highly diverse that they include forests, deserts, and coal mines, as well as in potatoes and cheese. Despite how widely Agrococcus occurs, the extent of its physiology, genomes, and potential roles in the environment are poorly understood. Here we use whole-genome analysis, chemotaxonomic markers, morphology, and 16S rRNA gene phylogeny to describe a new isolate of the genus Agrococcus from freshwater microbialites in Pavilion Lake, British Columbia, Canada. We characterize this isolate as a new species Agrococcus pavilionensis strain RW1 and provide the first complete genome from a member of the genus Agrococcus. The A. pavilionensis genome consists of one chromosome (2,627,177 bp) as well as two plasmids (HC-CG1 1,427 bp, and LC-RRW783 31,795 bp). The genome reveals considerable genetic promiscuity via mobile elements, including a prophage and plasmids involved in integration, transposition, and heavy-metal stress. A. pavilionensis strain RW1 differs from other members of the Agrococcus genus by having a novel phospholipid fatty acid iso-C15:1Δ⁴, β-galactosidase activity and amygdalin utilization. Carotenoid biosynthesis is predicted by genomic metabolic reconstruction, which explains the characteristic yellow pigmentation of A. pavilionensis. Metabolic reconstructions of strain RW1 genome predicts a pathway for releasing ammonia via ammonification amino acids, which could increase the saturation index leading to carbonate precipitation. Our genomic analyses suggest signatures of environmental adaption to the relatively cold and oligotrophic conditions of Pavilion Lake microbialites. A. pavilionensis strain RW1 in modern microbialites has an ecological significance in Pavilion Lake microbialites, which include potential roles in heavy-metal cycling and carbonate precipitation (e.g., ammonification of amino acids and filamentation which many trap carbonate minerals).

Protection induced by a Francisella tularensis subunit vaccine delivered by glucan particles

Francisella tularensis is an intracellular pathogen causing the disease tularemia, and an organism of concern to biodefence. There is no licensed vaccine available. Subunit approaches have failed to induce protection, which requires both humoral and cellular immune memory responses, and have been hampered by a lack of understanding as to which antigens are immunoprotective. We undertook a preliminary in silico analysis to identify candidate protein antigens. These antigens were then recombinantly expressed and encapsulated into glucan particles (GPs), purified Saccharomyces cerevisiae cell walls composed primarily of β-1,3-glucans. Immunological profiling in the mouse was used to down-selection to seven lead antigens: FTT1043 (Mip), IglC, FTT0814, FTT0438, FTT0071 (GltA), FTT0289, FTT0890 (PilA) prior to transitioning their evaluation to a Fischer 344 rat model for efficacy evaluation. F344 rats were vaccinated with the GP protein antigens co-delivered with GP-loaded with Francisella LPS. Measurement of cell mediated immune responses and computational epitope analysis allowed down-selection to three promising candidates: FTT0438, FTT1043 and FTT0814. Of these, a GP vaccine delivering Francisella LPS and the FTT0814 protein was able to induce protection in rats against an aerosol challenge of F. tularensis SchuS4, and reduced organ colonisation and clinical signs below that which immunisation with a GP-LPS alone vaccine provided. This is the first report of a protein supplementing protection induced by LPS in a Francisella vaccine. This paves the way for developing an effective, safe subunit vaccine for the prevention of inhalational tularemia, and validates the GP platform for vaccine delivery where complex immune responses are required for prevention of infections by intracellular pathogens.

Environmentally-driven gene content convergence and the Bacillus phylogeny

Background

Members of the Bacillus genus have been isolated from a variety of environments. However, the relationship between potential metabolism and the niche from which bacteria of this genus have been isolated has not been extensively studied. The existence of a monophyletic aquatic Bacillus group, composed of members isolated from both marine and fresh water has been proposed. Here, we present a phylogenetic/phylogenomic analysis to investigate the potential relationship between the environment from which group members have been isolated and their evolutionary origin. We also carried out hierarchical clustering based on functional content to test for potential environmental effects on the genetic content of these bacteria.

Results

The phylogenetic reconstruction showed that Bacillus strains classified as aquatic have evolutionary origins in different lineages. Although we observed the presence of a clade consisting exclusively of aquatic Bacillus, it is not comprised of the same strains previously reported. In contrast to phylogeny, clustering based on the functional categories of the encoded proteomes resulted in groups more compatible with the environments from which the organisms were isolated. This evidence suggests a detectable environmental influence on bacterial genetic content, despite their different evolutionary origins.

Conclusion

Our results suggest that aquatic Bacillus species have polyphyletic origins, but exhibit convergence at the gene content level.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1261-7) contains supplementary material, which is available to authorized users.

Phylogeny-corrected identification of microbial gene families relevant to human gut colonization

The mechanisms by which different microbes colonize the healthy human gut versus other body sites, the gut in disease states, or other environments remain largely unknown. Identifying microbial genes influencing fitness in the gut could lead to new ways to engineer probiotics or disrupt pathogenesis. We approach this problem by measuring the statistical association between a species having a gene and the probability that the species is present in the gut microbiome. The challenge is that closely related species tend to be jointly present or absent in the microbiome and also share many genes, only a subset of which are involved in gut adaptation. We show that this phylogenetic correlation indeed leads to many false discoveries and propose phylogenetic linear regression as a powerful solution. To apply this method across the bacterial tree of life, where most species have not been experimentally phenotyped, we use metagenomes from hundreds of people to quantify each species' prevalence in and specificity for the gut microbiome. This analysis reveals thousands of genes potentially involved in adaptation to the gut across species, including many novel candidates as well as processes known to contribute to fitness of gut bacteria, such as acid tolerance in Bacteroidetes and sporulation in Firmicutes. We also find microbial genes associated with a preference for the gut over other body sites, which are significantly enriched for genes linked to fitness in an in vivo competition experiment. Finally, we identify gene families associated with higher prevalence in patients with Crohn's disease, including Proteobacterial genes involved in conjugation and fimbria regulation, processes previously linked to inflammation. These gene targets may represent new avenues for modulating host colonization and disease. Our strategy of combining metagenomics with phylogenetic modeling is general and can be used to identify genes associated with adaptation to any environment.

The Human Virome Protein Cluster Database (HVPC): A Human Viral Metagenomic Database for Diversity and Function Annotation

Human virome, including those of bacteria (bacteriophages) have received an increasing attention recently, owing to the rapid developments in human microbiome research and the awareness of the far-reaching influence of microbiomes on health and disease. Nevertheless, human viromes are still underrepresented in literature making viruses a virtually untapped resource of diversity, functional and physiological information. Here we present the human virome protein cluster database as an effort to improve functional annotation and characterization of human viromes. The database was built out of hundreds of virome datasets from six different body sites. We also show the utility of this database through its use for the characterization of three bronchoalveolar lavage (BAL) viromes from one healthy control in addition to one moderate and one severe chronic obstructive pulmonary disease (COPD) patients. The use of the database allowed for a better functional annotation, which were otherwise poorly characterized when limited to annotation using sequences from full-length viral genomes. In addition, our BAL samples gave a first insight into viral communities of COPD patients and confirm a state of dysbiosis for viruses that increases with disease progression. Moreover, they shed light on the potential role of phages in the horizontal gene transfer of bacterial virulence factors, a phenomenon that highlights a possible contribution of phages to etiopathology.

In vivo growth and genomic characterization of rickettsia-like organisms isolated from farmed Chinook salmon (Oncorhynchus tshawytscha) in New Zealand

A rickettsia-like organism, designated NZ-RLO2, was isolated from Chinook salmon (Oncorhynchus tshawytscha) farmed in the South Island, New Zealand. In vivo growth showed NZ-RLO2 was able to grow in CHSE-214, EPC, BHK-21, C6/36 and Sf21 cell lines, while Piscirickettsia salmonis LF-89T grew in all but BHK-21 and Sf21. NZ-RLO2 grew optimally in EPC at 15°C, CHSE-214 and EPC at 18°C. The growth of LF-89 T was optimal at 15°C, 18°C and 22°C in CHSE-24, but appeared less efficient in EPC cells at all temperatures. Pan-genome comparison of predicted proteomes shows that available Chilean strains of P. salmonis grouped into two clusters (p-value = 94%). NZ-RLO2 was genetically different from previously described NZ-RLO1, and both strains grouped separately from the Chilean strains in one of the two clusters (p-value = 88%), but were closely related to each other. TaqMan and Sybr Green real-time PCR targeting RNA polymerase (rpoB) and DNA primase (dnaG), respectively, were developed to detect NZ-RLO2. This study indicates that the New Zealand strains showed a closer genetic relationship to one of the Chilean P. salmonis clusters; however, more Piscirickettsia genomes from wider geographical regions and diverse hosts are needed to better understand the classification within this genus.

Characterization of Antimicrobial-Producing Beneficial Bacteria Isolated from Huanglongbing Escape Citrus Trees

The microbiome associated with crop plants has a strong impact on their health and productivity. Candidatus Liberibacter asiaticus (Las), the bacterial pathogen responsible for Huanglongbing (HLB) disease, lives inside the phloem of citrus plants including the root system. It has been suggested that Las negatively affects citrus microbiome. On the other hand, members of citrus microbiome also influence the interaction between Las and citrus. Here, we report the isolation and characterization of multiple putative beneficial bacteria from healthy citrus rhizosphere. Firstly, six bacterial strains showing antibacterial activity against two bacteria closely related to Las: Agrobacterium tumefaciens and Sinorhizobium meliloti were selected. Among them, Burkholderia metallica strain A53 and Burkholderia territorii strain A63 are within the β-proteobacteria class, whereas Pseudomonas granadensis strain 100 and Pseudomonas geniculata strain 95 are within the γ-proteobacteria class. Additionally, two gram-positive bacteria Rhodococcus jialingiae strain 108 and Bacillus pumilus strain 104 were also identified. Secondly, antimicrobial activity against three fungal pathogens: Alternaria alternata, Colletotrichum acutatum, Phyllosticta citricarpa, and two oomycetes: Phytophthora nicotianae and Phytophthora palmivora. Four bacterial strains Burkholderia territorii A63, Burkholderia metallica A53, Pseudomonas geniculata 95, and Bacillus pumilus 104 were shown to have antagonistic activity against the citrus root pathogen Phytophthora nicotianae based on dual culture antagonist assays and compartmentalized petri dish assays. The four selected bacteria were sequenced. Genes involved in phosphate solubilization, siderophore production and iron acquisition, volatile organic compound production, osmoprotection and osmotic tolerance, phytohormone production, antagonism, and nutrient competition were predicted and discussed related to the beneficial traits.

A sensitive short read homology search tool for paired-end read sequencing data

Background

Homology search is still a significant step in functional analysis for genomic data. Profile Hidden Markov Model-based homology search has been widely used in protein domain analysis in many different species. In particular, with the fast accumulation of transcriptomic data of non-model species and metagenomic data, profile homology search is widely adopted in integrated pipelines for functional analysis. While the state-of-the-art tool HMMER has achieved high sensitivity and accuracy in domain annotation, the sensitivity of HMMER on short reads declines rapidly. The low sensitivity on short read homology search can lead to inaccurate domain composition and abundance computation. Our experimental results showed that half of the reads were missed by HMMER for a RNA-Seq dataset. Thus, there is a need for better methods to improve the homology search performance for short reads.

Results

We introduce a profile homology search tool named Short-Pair that is designed for short paired-end reads. By using an approximate Bayesian approach employing distribution of fragment lengths and alignment scores, Short-Pair can retrieve the missing end and determine true domains. In particular, Short-Pair increases the accuracy in aligning short reads that are part of remote homologs. We applied Short-Pair to a RNA-Seq dataset and a metagenomic dataset and quantified its sensitivity and accuracy on homology search. The experimental results show that Short-Pair can achieve better overall performance than the state-of-the-art methodology of profile homology search.

Conclusions

Short-Pair is best used for next-generation sequencing (NGS) data that lack reference genomes. It provides a complementary paired-end read homology search tool to HMMER. The source code is freely available at https://sourceforge.net/projects/short-pair/.

MEGGASENSE - The Metagenome/Genome Annotated Sequence Natural Language Search Engine: A Platform for the Construction of Sequence Data Warehouses

The MEGGASENSE platform constructs relational databases of DNA or protein sequences. The default functional analysis uses 14 106 hidden Markov model (HMM) profiles based on sequences in the KEGG database. The Solr search engine allows sophisticated queries and a BLAST search function is also incorporated. These standard capabilities were used to generate the SCATT database from the predicted proteome of Streptomyces cattleya. The implementation of a specialised metagenome database (AMYLOMICS) for bioprospecting of carbohydrate-modifying enzymes is described. In addition to standard assembly of reads, a novel 'functional' assembly was developed, in which screening of reads with the HMM profiles occurs before the assembly. The AMYLOMICS database incorporates additional HMM profiles for carbohydrate-modifying enzymes and it is illustrated how the combination of HMM and BLAST analyses helps identify interesting genes. A variety of different proteome and metagenome databases have been generated by MEGGASENSE.

The In-Feed Antibiotic Carbadox Induces Phage Gene Transcription in the Swine Gut Microbiome

Carbadox is a quinoxaline-di-N-oxide antibiotic fed to over 40% of young pigs in the United States that has been shown to induce phage DNA transduction in vitro; however, the effects of carbadox on swine microbiome functions are poorly understood. We investigated the in vivo longitudinal effects of carbadox on swine gut microbial gene expression (fecal metatranscriptome) and phage population dynamics (fecal dsDNA viromes). Microbial metagenome, transcriptome, and virome sequences were annotated for taxonomic inference and gene function by using FIGfam (isofunctional homolog sequences) and SEED subsystems databases. When the beta diversities of microbial FIGfam annotations were compared, the control and carbadox communities were distinct 2 days after carbadox introduction. This effect was driven by carbadox-associated lower expression of FIGfams (n = 66) related to microbial respiration, carbohydrate utilization, and RNA metabolism (q < 0.1), suggesting bacteriostatic or bactericidal effects within certain populations. Interestingly, carbadox treatment caused greater expression of FIGfams related to all stages of the phage lytic cycle 2 days following the introduction of carbadox (q ≤0.07), suggesting the carbadox-mediated induction of prophages and phage DNA recombination. These effects were diminished by 7 days of continuous carbadox in the feed, suggesting an acute impact. Additionally, the viromes included a few genes that encoded resistance to tetracycline, aminoglycoside, and beta-lactam antibiotics but these did not change in frequency over time or with treatment. The results show decreased bacterial growth and metabolism, prophage induction, and potential transduction of bacterial fitness genes in swine gut bacterial communities as a result of carbadox administration.

FDA regulations on agricultural antibiotic use have focused on antibiotics that are important for human medicine. Carbadox is an antibiotic not used in humans but frequently used on U.S. pig farms. It is important to study possible side effects of carbadox use because it has been shown to promote bacterial evolution, which could indirectly impact antibiotic resistance in bacteria of clinical importance. Interestingly, the present study shows greater prophage gene expression in feces from carbadox-fed animals than in feces from nonmedicated animals 2 days after the initiation of in-feed carbadox treatment. Importantly, the phage genetic material isolated in this study contained genes that could provide resistance to antibiotics that are important in human medicine, indicating that human-relevant antibiotic resistance genes are mobile between bacteria via phages. This study highlights the collateral effects of antibiotics and demonstrates the need to consider diverse antibiotic effects whenever antibiotics are being used or new regulations are considered.

Improve homology search sensitivity of PacBio data by correcting frameshifts

MOTIVATION

Single-molecule, real-time sequencing (SMRT) developed by Pacific BioSciences produces longer reads than secondary generation sequencing technologies such as Illumina. The long read length enables PacBio sequencing to close gaps in genome assembly, reveal structural variations, and identify gene isoforms with higher accuracy in transcriptomic sequencing. However, PacBio data has high sequencing error rate and most of the errors are insertion or deletion errors. During alignment-based homology search, insertion or deletion errors in genes will cause frameshifts and may only lead to marginal alignment scores and short alignments. As a result, it is hard to distinguish true alignments from random alignments and the ambiguity will incur errors in structural and functional annotation. Existing frameshift correction tools are designed for data with much lower error rate and are not optimized for PacBio data. As an increasing number of groups are using SMRT, there is an urgent need for dedicated homology search tools for PacBio data.

RESULTS

In this work, we introduce Frame-Pro, a profile homology search tool for PacBio reads. Our tool corrects sequencing errors and also outputs the profile alignments of the corrected sequences against characterized protein families. We applied our tool to both simulated and real PacBio data. The results showed that our method enables more sensitive homology search, especially for PacBio data sets of low sequencing coverage. In addition, we can correct more errors when comparing with a popular error correction tool that does not rely on hybrid sequencing.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available at https://sourceforge.net/projects/frame-pro/

CONTACT

yannisun@msu.edu.

Improvements to PATRIC, the all-bacterial Bioinformatics Database and Analysis Resource Center

The Pathosystems Resource Integration Center (PATRIC) is the bacterial Bioinformatics Resource Center (https://www.patricbrc.org). Recent changes to PATRIC include a redesign of the web interface and some new services that provide users with a platform that takes them from raw reads to an integrated analysis experience. The redesigned interface allows researchers direct access to tools and data, and the emphasis has changed to user-created genome-groups, with detailed summaries and views of the data that researchers have selected. Perhaps the biggest change has been the enhanced capability for researchers to analyze their private data and compare it to the available public data. Researchers can assemble their raw sequence reads and annotate the contigs using RASTtk. PATRIC also provides services for RNA-Seq, variation, model reconstruction and differential expression analysis, all delivered through an updated private workspace. Private data can be compared by ‘virtual integration’ to any of PATRIC's public data. The number of genomes available for comparison in PATRIC has expanded to over 80 000, with a special emphasis on genomes with antimicrobial resistance data. PATRIC uses this data to improve both subsystem annotation and k-mer classification, and tags new genomes as having signatures that indicate susceptibility or resistance to specific antibiotics.

Metagenomics: Probing pollutant fate in natural and engineered ecosystems

Polluted environments are a reservoir of microbial species able to degrade or to convert pollutants to harmless compounds. The proper management of microbial resources requires a comprehensive characterization of their genetic pool to assess the fate of contaminants and increase the efficiency of bioremediation processes. Metagenomics offers appropriate tools to describe microbial communities in their whole complexity without lab-based cultivation of individual strains. After a decade of use of metagenomics to study microbiomes, the scientific community has made significant progress in this field. In this review, we survey the main steps of metagenomics applied to environments contaminated with organic compounds or heavy metals. We emphasize technical solutions proposed to overcome encountered obstacles. We then compare two metagenomic approaches, i.e. library-based targeted metagenomics and direct sequencing of metagenomes. In the former, environmental DNA is cloned inside a host, and then clones of interest are selected based on (i) their expression of biodegradative functions or (ii) sequence homology with probes and primers designed from relevant, already known sequences. The highest score for the discovery of novel genes and degradation pathways has been achieved so far by functional screening of large clone libraries. On the other hand, direct sequencing of metagenomes without a cloning step has been more often applied to polluted environments for characterization of the taxonomic and functional composition of microbial communities and their dynamics. In this case, the analysis has focused on 16S rRNA genes and marker genes of biodegradation. Advances in next generation sequencing and in bioinformatic analysis of sequencing data have opened up new opportunities for assessing the potential of biodegradation by microbes, but annotation of collected genes is still hampered by a limited number of available reference sequences in databases. Although metagenomics is still facing technical and computational challenges, our review of the recent literature highlights its value as an aid to efficiently monitor the clean-up of contaminated environments and develop successful strategies to mitigate the impact of pollutants on ecosystems.

An integrated metagenomics pipeline for strain profiling reveals novel patterns of bacterial transmission and biogeography

We present the Metagenomic Intra-species Diversity Analysis System (MIDAS), which is an integrated computational pipeline for quantifying bacterial species abundance and strain-level genomic variation, including gene content and single-nucleotide polymorphisms (SNPs), from shotgun metagenomes. Our method leverages a database of more than 30,000 bacterial reference genomes that we clustered into species groups. These cover the majority of abundant species in the human microbiome but only a small proportion of microbes in other environments, including soil and seawater. We applied MIDAS to stool metagenomes from 98 Swedish mothers and their infants over one year and used rare SNPs to track strains between hosts. Using this approach, we found that although species compositions of mothers and infants converged over time, strain-level similarity diverged. Specifically, early colonizing bacteria were often transmitted from an infant’s mother, while late colonizing bacteria were often transmitted from other sources in the environment and were enriched for spore-formation genes. We also applied MIDAS to 198 globally distributed marine metagenomes and used gene content to show that many prevalent bacterial species have population structure that correlates with geographic location. Strain-level genetic variants present in metagenomes clearly reveal extensive structure and dynamics that are obscured when data are analyzed at a coarser taxonomic resolution.

Comparative Genome Analysis of Fusobacterium nucleatum

Fusobacterium nucleatum is considered to be a key oral bacterium in recruiting periodontal pathogens into subgingival dental plaque. Currently F. nucleatum can be subdivided into five subspecies. Our previous genome analysis of F. nucleatum W1481 (referred to hereafter as W1481), isolated from an 8-mm periodontal pocket in a patient with chronic periodontitis, suggested the possibility of a new subspecies. To further investigate the biology and relationships of this possible subspecies with other known subspecies, we performed comparative analysis between W1481 and 35 genome sequences represented by the five known Fusobacterium subspecies. Our analyses suggest that W1481 is most likely a new F. nucleatum subspecies, supported by evidence from phylogenetic analyses and maximal unique match indices (MUMi). Interestingly, we found a horizontally transferred W1481-specific genomic island harboring the tripartite ATP-independent (TRAP)-like transporter genes, suggesting this bacterium might have a high-affinity transport system for the C4-dicarboxylates malate, succinate, and fumarate. Moreover, we found virulence genes in the W1481 genome that may provide a strong defense mechanism which might enable it to colonize and survive within the host by evading immune surveillance. This comparative study provides better understanding of F. nucleatum and the basis for future functional work on this important pathogen.

Functional Annotations of Paralogs: A Blessing and a Curse

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

A Primer on Infectious Disease Bacterial Genomics

The number of large-scale genomics projects is increasing due to the availability of affordable high-throughput sequencing (HTS) technologies. The use of HTS for bacterial infectious disease research is attractive because one whole-genome sequencing (WGS) run can replace multiple assays for bacterial typing, molecular epidemiology investigations, and more in-depth pathogenomic studies. The computational resources and bioinformatics expertise required to accommodate and analyze the large amounts of data pose new challenges for researchers embarking on genomics projects for the first time. Here, we present a comprehensive overview of a bacterial genomics projects from beginning to end, with a particular focus on the planning and computational requirements for HTS data, and provide a general understanding of the analytical concepts to develop a workflow that will meet the objectives and goals of HTS projects.

HoloVir: A Workflow for Investigating the Diversity and Function of Viruses in Invertebrate Holobionts

Abundant bioinformatics resources are available for the study of complex microbial metagenomes, however their utility in viral metagenomics is limited. HoloVir is a robust and flexible data analysis pipeline that provides an optimized and validated workflow for taxonomic and functional characterization of viral metagenomes derived from invertebrate holobionts. Simulated viral metagenomes comprising varying levels of viral diversity and abundance were used to determine the optimal assembly and gene prediction strategy, and multiple sequence assembly methods and gene prediction tools were tested in order to optimize our analysis workflow. HoloVir performs pairwise comparisons of single read and predicted gene datasets against the viral RefSeq database to assign taxonomy and additional comparison to phage-specific and cellular markers is undertaken to support the taxonomic assignments and identify potential cellular contamination. Broad functional classification of the predicted genes is provided by assignment of COG microbial functional category classifications using EggNOG and higher resolution functional analysis is achieved by searching for enrichment of specific Swiss-Prot keywords within the viral metagenome. Application of HoloVir to viral metagenomes from the coral Pocillopora damicornis and the sponge Rhopaloeides odorabile demonstrated that HoloVir provides a valuable tool to characterize holobiont viral communities across species, environments, or experiments.

Capturing One of the Human Gut Microbiome's Most Wanted: Reconstructing the Genome of a Novel Butyrate-Producing, Clostridial Scavenger from Metagenomic Sequence Data

The role of the microbiome in health and disease is attracting great attention, yet we still know little about some of the most prevalent microorganisms inside our bodies. Several years ago, Human Microbiome Project (HMP) researchers generated a list of "most wanted" taxa: bacteria both prevalent among healthy volunteers and distantly related to any sequenced organisms. Unfortunately, the challenge of assembling high-quality genomes from a tangle of metagenomic reads has slowed progress in learning about these uncultured bacteria. Here, we describe how recent advances in sequencing and analysis allowed us to assemble "most wanted" genomes from metagenomic data collected from four stool samples. Using a combination of both de novo and guided assembly methods, we assembled and binned over 100 genomes from an initial data set of over 1,300 Gbp. One of these genome bins, which met HMP's criteria for a "most wanted" taxa, contained three essentially complete genomes belonging to a previously uncultivated species. This species is most closely related to Eubacterium desmolans and the clostridial cluster IV/Clostridium leptum subgroup species Butyricicoccus pullicaecorum (71-76% average nucleotide identity). Gene function analysis indicates that the species is an obligate anaerobe, forms spores, and produces the anti-inflammatory short-chain fatty acids acetate and butyrate. It also appears to take up metabolically costly molecules such as cobalamin, methionine, and branch-chained amino acids from the environment, and to lack virulence genes. Thus, the evidence is consistent with a secondary degrader that occupies a host-dependent, nutrient-scavenging niche within the gut; its ability to produce butyrate, which is thought to play an anti-inflammatory role, makes it intriguing for the study of diseases such as colon cancer and inflammatory bowel disease. In conclusion, we have assembled essentially complete genomes from stool metagenomic data, yielding valuable information about uncultured organisms' metabolic and ecologic niches, factors that may be required to successfully culture these bacteria, and their role in maintaining health and causing disease.

Genomic insights into members of the candidate phylum Hyd24-12 common in mesophilic anaerobic digesters

Members of the candidate phylum Hyd24-12 are globally distributed, but no genomic information or knowledge about their morphology, physiology or ecology is available. In this study, members of the Hyd24-12 lineage were shown to be present and abundant in full-scale mesophilic anaerobic digesters at Danish wastewater treatment facilities. In some samples, a member of the Hyd24-12 lineage was one of the most abundant genus-level bacterial taxa, accounting for up to 8% of the bacterial biomass. Three closely related and near-complete genomes were retrieved using metagenome sequencing of full-scale anaerobic digesters. Genome annotation and metabolic reconstruction showed that they are Gram-negative bacteria likely involved in acidogenesis, producing acetate and hydrogen from fermentation of sugars, and may play a role in the cycling of sulphur in the digesters. Fluorescence in situ hybridization revealed single rod-shaped cells dispersed within the flocs. The genomic information forms a foundation for a more detailed understanding of their role in anaerobic digestion and provides the first insight into a hitherto undescribed branch in the tree of life.

Genome Sequence of a Clinical Staphylococcus aureus Strain from a Prosthetic Joint Infection

Here, we report the genome sequence of Staphylococcus aureus LYO-S2, an isolate with sequence type (ST) 45 that was isolated in 2001 from a prosthetic joint infection.

PATtyFams: Protein Families for the Microbial Genomes in the PATRIC Database

The ability to build accurate protein families is a fundamental operation in bioinformatics that influences comparative analyses, genome annotation, and metabolic modeling. For several years we have been maintaining protein families for all microbial genomes in the PATRIC database (Pathosystems Resource Integration Center, patricbrc.org) in order to drive many of the comparative analysis tools that are available through the PATRIC website. However, due to the burgeoning number of genomes, traditional approaches for generating protein families are becoming prohibitive. In this report, we describe a new approach for generating protein families, which we call PATtyFams. This method uses the k-mer-based function assignments available through RAST (Rapid Annotation using Subsystem Technology) to rapidly guide family formation, and then differentiates the function-based groups into families using a Markov Cluster algorithm (MCL). This new approach for generating protein families is rapid, scalable and has properties that are consistent with alignment-based methods.

MycoCAP - Mycobacterium Comparative Analysis Platform

Mycobacterium spp. are renowned for being the causative agent of diseases like leprosy, Buruli ulcer and tuberculosis in human beings. With more and more mycobacterial genomes being sequenced, any knowledge generated from comparative genomic analysis would provide better insights into the biology, evolution, phylogeny and pathogenicity of this genus, thus helping in better management of diseases caused by Mycobacterium spp.With this motivation, we constructed MycoCAP, a new comparative analysis platform dedicated to the important genus Mycobacterium. This platform currently provides information of 2108 genome sequences of at least 55 Mycobacterium spp. A number of intuitive web-based tools have been integrated in MycoCAP particularly for comparative analysis including the PGC tool for comparison between two genomes, PathoProT for comparing the virulence genes among the Mycobacterium strains and the SuperClassification tool for the phylogenic classification of the Mycobacterium strains and a specialized classification system for strains of Mycobacterium abscessus. We hope the broad range of functions and easy-to-use tools provided in MycoCAP makes it an invaluable analysis platform to speed up the research discovery on mycobacteria for researchers. Database URL: http://mycobacterium.um.edu.my

Automated and Accurate Estimation of Gene Family Abundance from Shotgun Metagenomes

Shotgun metagenomic DNA sequencing is a widely applicable tool for characterizing the functions that are encoded by microbial communities. Several bioinformatic tools can be used to functionally annotate metagenomes, allowing researchers to draw inferences about the functional potential of the community and to identify putative functional biomarkers. However, little is known about how decisions made during annotation affect the reliability of the results. Here, we use statistical simulations to rigorously assess how to optimize annotation accuracy and speed, given parameters of the input data like read length and library size. We identify best practices in metagenome annotation and use them to guide the development of the Shotgun Metagenome Annotation Pipeline (ShotMAP). ShotMAP is an analytically flexible, end-to-end annotation pipeline that can be implemented either on a local computer or a cloud compute cluster. We use ShotMAP to assess how different annotation databases impact the interpretation of how marine metagenome and metatranscriptome functional capacity changes across seasons. We also apply ShotMAP to data obtained from a clinical microbiome investigation of inflammatory bowel disease. This analysis finds that gut microbiota collected from Crohn’s disease patients are functionally distinct from gut microbiota collected from either ulcerative colitis patients or healthy controls, with differential abundance of metabolic pathways related to host-microbiome interactions that may serve as putative biomarkers of disease.

Microbial communities perform a wide variety of functions, from marine photosynthesis to aiding digestion in the human gut. Shotgun “metagenomic” sequencing can be used to sample millions of short DNA sequences from such communities directly, without needing to first culture its constituents in the laboratory. Using these data, researchers can survey which functions are encoded by mapping these short sequences to known protein families and pathways. Several tools for this annotation already exist. But, annotation is a multi-step process that includes identification of genes in a metagenome and determination of the type of protein each gene encodes. We currently know little about how different choices of parameters during annotation influences the final results. In this work, we systematically test how several key decisions affect the accuracy and speed of annotation, and based on these results, develop new software for annotation, which we named ShotMAP. We then use ShotMAP to functionally characterize marine communities and gut communities in a clinical cohort of inflammatory bowel disease. We find several functions are differentially represented in the gut microbiome of Crohn’s disease patients, which could be candidates for biomarkers and could also offer insight into the pathophysiology of Crohn’s. ShotMAP is freely available (https://github.com/sharpton/shotmap).

Genome reconstructions indicate the partitioning of ecological functions inside a phytoplankton bloom in the Amundsen Sea, Antarctica

Antarctica polynyas support intense phytoplankton blooms, impacting their environment by a substantial depletion of inorganic carbon and nutrients. These blooms are dominated by the colony-forming haptophyte Phaeocystis antarctica and they are accompanied by a distinct bacterial population. Yet, the ecological role these bacteria may play in P. antarctica blooms awaits elucidation of their functional gene pool and of the geochemical activities they support. Here, we report on a metagenome (~160 million reads) analysis of the microbial community associated with a P. antarctica bloom event in the Amundsen Sea polynya (West Antarctica). Genomes of the most abundant Bacteroidetes and Proteobacteria populations have been reconstructed and a network analysis indicates a strong functional partitioning of these bacterial taxa. Three of them (SAR92, and members of the Oceanospirillaceae and Cryomorphaceae) are found in close association with P. antarctica colonies. Distinct features of their carbohydrate, nitrogen, sulfur and iron metabolisms may serve to support mutualistic relationships with P. antarctica. The SAR92 genome indicates a specialization in the degradation of fatty acids and dimethylsulfoniopropionate (compounds released by P. antarctica) into dimethyl sulfide, an aerosol precursor. The Oceanospirillaceae genome carries genes that may enhance algal physiology (cobalamin synthesis). Finally, the Cryomorphaceae genome is enriched in genes that function in cell or colony invasion. A novel pico-eukaryote, Micromonas related genome (19.6 Mb, ~94% completion) was also recovered. It contains the gene for an anti-freeze protein, which is lacking in Micromonas at lower latitudes. These draft genomes are representative for abundant microbial taxa across the Southern Ocean surface.

Average genome size estimation improves comparative metagenomics and sheds light on the functional ecology of the human microbiome

Average genome size is an important, yet often overlooked, property of microbial communities. We developed MicrobeCensus to rapidly and accurately estimate average genome size from shotgun metagenomic data and applied our tool to 1,352 human microbiome samples. We found that average genome size differs significantly within and between body sites and tracks with major functional and taxonomic differences. In the gut, average genome size is positively correlated with the abundance of Bacteroides and genes related to carbohydrate metabolism. Importantly, we found that average genome size variation can bias comparative analyses, and that normalization improves detection of differentially abundant genes.

Insights on the emergence of Mycobacterium tuberculosis from the analysis of Mycobacterium kansasii

By phylogenetic analysis, Mycobacterium kansasii is closely related to Mycobacterium tuberculosis. Yet, although both organisms cause pulmonary disease, M. tuberculosis is a global health menace, whereas M. kansasii is an opportunistic pathogen. To illuminate the differences between these organisms, we have sequenced the genome of M. kansasii ATCC 12478 and its plasmid (pMK12478) and conducted side-by-side in vitro and in vivo investigations of these two organisms. The M. kansasii genome is 6,432,277 bp, more than 2 Mb longer than that of M. tuberculosis H37Rv, and the plasmid contains 144,951 bp. Pairwise comparisons reveal conserved and discordant genes and genomic regions. A notable example of genomic conservation is the virulence locus ESX-1, which is intact and functional in the low-virulence M. kansasii, potentially mediating phagosomal disruption. Differences between these organisms include a decreased predicted metabolic capacity, an increased proportion of toxin–antitoxin genes, and the acquisition of M. tuberculosis-specific genes in the pathogen since their common ancestor. Consistent with their distinct epidemiologic profiles, following infection of C57BL/6 mice, M. kansasii counts increased by less than 10-fold over 6 weeks, whereas M. tuberculosis counts increased by over 10,000-fold in just 3 weeks. Together, these data suggest that M. kansasii can serve as an image of the environmental ancestor of M. tuberculosis before its emergence as a professional pathogen, and can be used as a model organism to study the switch from an environmental opportunistic pathogen to a professional host-restricted pathogen.