Recent physiological and ecological studies have challenged the long-held belief that microbial metabolism of molecular hydrogen (H 2 ) is a niche process. To gain a broader insight into the importance of microbial H 2 metabolism, we comprehensively surveyed the genomic and metagenomic distribution of hydrogenases, the reversible enzymes that catalyse the oxidation and evolution of H 2 . The protein sequences of 3286 non-redundant putative hydrogenases were curated from publicly available databases. These metalloenzymes were classified into multiple groups based on (1) amino acid sequence phylogeny, (2) metal-binding motifs, (3) predicted genetic organisation and (4) reported biochemical characteristics. Four groups (22 subgroups) of [NiFe]-hydrogenase, three groups (6 subtypes) of [FeFe]-hydrogenases and a small group of [Fe]-hydrogenases were identified. We predict that this hydrogenase diversity supports H 2 -based respiration, fermentation and carbon fixation processes in both oxic and anoxic environments, in addition to various H 2 -sensing, electron-bifurcation and energy-conversion mechanisms. Hydrogenase-encoding genes were identified in 51 bacterial and archaeal phyla, suggesting strong pressure for both vertical and lateral acquisition. Furthermore, hydrogenase genes could be recovered from diverse terrestrial, aquatic and host-associated metagenomes in varying proportions, indicating a broad ecological distribution and utilisation. Oxygen content (pO 2 ) appears to be a central factor driving the phylum-and ecosystem-level distribution of these genes. In addition to compounding evidence that H 2 was the first electron donor for life, our analysis suggests that the great diversification of hydrogenases has enabled H 2 metabolism to sustain the growth or survival of microorganisms in a wide range of ecosystems to the present day. This work also provides a comprehensive expanded system for classifying hydrogenases and identifies new prospects for investigating H 2 metabolism.

CRISPR–Cas systems provide bacteria with adaptive immunity against foreign nucleic acids by acquiring short, invader-derived sequences called spacers. Here, we use high-throughput sequencing to analyse millions of spacer acquisition events in wild-type populations of Pectobacterium atrosepticum. Plasmids not previously encountered, or plasmids that had escaped CRISPR–Cas targeting via point mutation, are used to provoke naive or primed spacer acquisition, respectively. The origin, location and order of spacer acquisition show that spacer selection through priming initiates near the site of CRISPR–Cas recognition (the protospacer), but on the displaced strand, and is consistent with 3′–5′ translocation of the Cas1:Cas2-3 acquisition machinery. Newly acquired spacers determine the location and strand specificity of subsequent spacers and demonstrate that interference-driven spacer acquisition (‘targeted acquisition') is a major contributor to adaptation in type I-F CRISPR–Cas systems. Finally, we show that acquisition of self-targeting spacers is occurring at a constant rate in wild-type cells and can be triggered by foreign DNA with similarity to the bacterial chromosome.

BackgroundCRISPR (clustered regularly interspaced short palindromic repeats) RNAs provide the specificity for noncoding RNA-guided adaptive immune defence systems in prokaryotes. CRISPR arrays consist of repeat sequences separated by specific spacer sequences. CRISPR arrays have previously been identified in a large proportion of prokaryotic genomes. However, currently available detection algorithms do not utilise recently discovered features regarding CRISPR loci.ResultsWe have developed a new approach to automatically detect, predict and interactively refine CRISPR arrays. It is available as a web program and command line from bioanalysis.otago.ac.nz/CRISPRDetect. CRISPRDetect discovers putative arrays, extends the array by detecting additional variant repeats, corrects the direction of arrays, refines the repeat/spacer boundaries, and annotates different types of sequence variations (e.g. insertion/deletion) in near identical repeats. Due to these features, CRISPRDetect has significant advantages when compared to existing identification tools. As well as further support for small medium and large repeats, CRISPRDetect identified a class of arrays with ‘extra-large’ repeats in bacteria (repeats 44–50 nt). The CRISPRDetect output is integrated with other analysis tools. Notably, the predicted spacers can be directly utilised by CRISPRTarget to predict targets.ConclusionCRISPRDetect enables more accurate detection of arrays and spacers and its gff output is suitable for inclusion in genome annotation pipelines and visualisation. It has been used to analyse all complete bacterial and archaeal reference genomes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2627-0) contains supplementary material, which is available to authorized users.

The bacterial and archaeal CRISPR/Cas adaptive immune system targets specific protospacer nucleotide sequences in invading organisms. This requires base pairing between processed CRISPR RNA and the target protospacer. For type I and II CRISPR/Cas systems, protospacer adjacent motifs (PAM) are essential for target recognition, and for type III, mismatches in the flanking sequences are important in the antiviral response. In this study, we examine the properties of each class of CRISPR. We use this information to provide a tool (CRISPRTarget) that predicts the most likely targets of CRISPR RNAs (http://bioanalysis.otago.ac.nz/CRISPRTarget). This can be used to discover targets in newly sequenced genomic or metagenomic data. To test its utility, we discover features and targets of well-characterized Streptococcus thermophilus and Sulfolobus solfataricus type II and III CRISPR/Cas systems. Finally, in Pectobacterium species, we identify new CRISPR targets and propose a model of temperate phage exposure and subsequent inhibition by the type I CRISPR/Cas systems.

Microbial molecular hydrogen (H 2 ) cycling is central to metabolic homeostasis and microbial composition in the human gastrointestinal tract. Molecular H 2 is produced as an endproduct of carbohydrate fermentation and is reoxidised primarily by sulfate-reduction, acetogenesis, and methanogenesis. However, the enzymatic basis for these processes is incompletely understood and the hydrogenases responsible have not been investigated. In this work, we surveyed the genomic and metagenomic distribution of hydrogenase-encoding genes in the human colon to infer dominant mechanisms of H 2 cycling. The data demonstrate that 70% of gastrointestinal microbial species listed in the Human Microbiome Project encode the genetic capacity to metabolise H 2 . A wide variety of anaerobicallyadapted hydrogenases were present, with [FeFe]-hydrogenases predominant. We subsequently analyzed the hydrogenase gene content of stools from 20 healthy human subjects. The hydrogenase gene content of all samples was overwhelmingly dominated by fermentative and electron-bifurcating [FeFe]-hydrogenases emerging from the Bacteroidetes and Firmicutes. This study supports that H 2 metabolism in the human gut is driven by fermentative H 2 production and interspecies H 2 transfer. However, it suggests that electron-bifurcation rather than respiration is the dominant mechanism of H 2 reoxidation in the human colon, generating reduced ferredoxin to sustain carbon-fixation (e.g. acetogenesis) and respiration (via the Rnf complex). This work provides the first comprehensive bioinformatic insight into the mechanisms of H 2 metabolism in the human colon.