Targeted hypermutation of putative antigen sensors in multicellular bacteria.

Diversity-generating retroelements (DGRs) are used by bacteria, archaea, and viruses as a targeted mutagenesis tool. Through error-prone reverse transcription, DGRs introduce random mutations at specific genomic loci, enabling rapid evolution of these targeted genes. However, the function and benefits of DGR-diversified proteins in cellular hosts remain elusive. We find that 82% of DGRs from one of the major monophyletic lineages of DGR reverse transcriptases are encoded by multicellular bacteria, which often have two or more DGR loci in their genomes. Using the multicellular purple sulfur bacterium Thiohalocapsa sp. PB-PSB1 as an example, we characterized nine distinct DGR loci capable of generating 10282 different combinations of target proteins. With environmental metagenomes from individual Thiohalocapsa aggregates, we show that most of PB-PSB1's DGR target genes are diversified across its biogeographic range, with spatial heterogeneity in the diversity of each locus. In Thiohalocapsa PB-PSB1 and other bacteria hosting this lineage of cellular DGRs, the diversified target genes are associated with NACHT-domain anti-phage defenses and putative ternary conflict systems previously shown to be enriched in multicellular bacteria. We propose that these DGR-diversified targets act as antigen sensors that confer a form of adaptive immunity to their multicellular consortia, though this remains to be experimentally tested. These findings could have implications for understanding the evolution of multicellularity, as the NACHT-domain anti-phage systems and ternary systems share both domain homology and conceptual similarities with the innate immune and programmed cell death pathways of plants and metazoans.

Horizontal Gene Transfer and CRISPR Targeting Drive Phage-Bacterial Host Interactions and Coevolution in "Pink Berry" Marine Microbial Aggregates.

Bacteriophages (phages), which are viruses that infect bacteria, are the most abundant components of microbial communities and play roles in community dynamics and host evolution. However, the study of phage-host interactions is hindered by a paucity of model systems from natural environments. Here, we investigate phage-host interactions in the "pink berry" consortia, which are naturally occurring, low-diversity, macroscopic bacterial aggregates that are found in the Sippewissett Salt Marsh (Falmouth, MA, USA). We leverage metagenomic sequence data and a comparative genomics approach to identify eight compete phage genomes, infer their bacterial hosts from host-encoded clustered regularly interspaced short palindromic repeats (CRISPRs), and observe the potential evolutionary consequences of these interactions. Seven of the eight phages identified infect known pink berry symbionts, namely, Desulfofustis sp. PB-SRB1, Thiohalocapsa sp. PB-PSB1, and Rhodobacteraceae sp. A2, and they are largely divergent from known viruses. In contrast to the conserved bacterial community structure of pink berries, the distribution of these phages across aggregates is highly variable. Two phages persisted over a period of seven years with high sequence conservation, allowing us to identify gene gain and loss. Increased nucleotide variation in a conserved phage capsid gene that is commonly targeted by host CRISPR systems suggests that CRISPRs may drive phage evolution in pink berries. Finally, we identified a predicted phage lysin gene that was horizontally transferred to its bacterial host, potentially via a transposon intermediary. Taken together, our results demonstrate that pink berry consortia contain diverse and variable phages as well as provide evidence for phage-host coevolution via multiple mechanisms in a natural microbial system. IMPORTANCE Phages, which are viruses that infect bacteria, are important components of all microbial systems, in which they drive the turnover of organic matter by lysing host cells, facilitate horizontal gene transfer (HGT), and coevolve with their bacterial hosts. Bacteria resist phage infection, which is often costly or lethal, through a diversity of mechanisms. One of these mechanisms is CRISPR systems, which encode arrays of phage-derived sequences from past infections to block subsequent infection with related phages. Here, we investigate the bacteria and phage populations from a simple marine microbial community, known as "pink berries", found in salt marshes of Falmouth, Massachusetts, as a model of phage-host coevolution. We identify eight novel phages and characterize a case of putative CRISPR-driven phage evolution as well as an instance of HGT between a phage and its host, together suggesting that phages have large evolutionary impacts in a naturally occurring microbial community.

Host-associated microbiomes drive structure and function of marine ecosystems.

The significance of symbioses between eukaryotic hosts and microbes extends from the organismal to the ecosystem level and underpins the health of Earth's most threatened marine ecosystems. Despite rapid growth in research on host-associated microbes, from individual microbial symbionts to host-associated consortia of significantly relevant taxa, little is known about their interactions with the vast majority of marine host species. We outline research priorities to strengthen our current knowledge of host-microbiome interactions and how they shape marine ecosystems. We argue that such advances in research will help predict responses of species, communities, and ecosystems to stressors driven by human activity and inform future management strategies.

Elevated pCO2 alters marine heterotrophic bacterial community composition and metabolic potential in response to a pulse of phytoplankton organic matter.

Factors that affect the respiration of organic carbon by marine bacteria can alter the extent to which the oceans act as a sink of atmospheric carbon dioxide. We designed seawater dilution experiments to assess the effect of pCO2 enrichment on heterotrophic bacterial community composition and metabolic potential in response to a pulse of phytoplankton-derived organic carbon. Experiments included treatments of elevated (1000 p.p.m.) and low (250 p.p.m.) pCO2 amended with 10 µmol L-1 dissolved organic carbon from Emiliana huxleyi lysates, and were conducted using surface-seawater collected from the South Pacific Subtropical Gyre. To assess differences in community composition and metabolic potential, shotgun metagenomic libraries were sequenced from low and elevated pCO2 treatments collected at the start of the experiment and following exponential growth. Our results indicate bacterial communities changed markedly in response to the organic matter pulse over time and were significantly affected by pCO2 enrichment. Elevated pCO2 also had disproportionate effects on the abundance of sequences related to proton pumps, carbohydrate metabolism, modifications of the phospholipid bilayer, resistance to toxic compounds and conjugative transfer. These results contribute to a growing understanding of the effects of elevated pCO2 on bacteria-mediated carbon cycling during phytoplankton bloom conditions in the marine environment.

Primer selection impacts specific population abundances but not community dynamics in a monthly time-series 16S rRNA gene amplicon analysis of coastal marine bacterioplankton.

Primers targeting the 16S small subunit ribosomal RNA marker gene, used to characterize bacterial and archaeal communities, have recently been re-evaluated for marine planktonic habitats. To investigate whether primer selection affects the ecological interpretation of bacterioplankton populations and community dynamics, amplicon sequencing with four primer sets targeting several hypervariable regions of the 16S rRNA gene was conducted on both mock communities constructed from cloned 16S rRNA genes and a time-series of DNA samples from the temperate coastal Santa Barbara Channel. Ecological interpretations of community structure (delineation of depth and seasonality, correlations with environmental factors) were similar across primer sets, while population dynamics varied. We observed substantial differences in relative abundances of taxa known to be poorly resolved by some primer sets, such as Thaumarchaeota and SAR11, and unexpected taxa including Roseobacter clades. Though the magnitude of relative abundances of common OTUs differed between primer sets, the relative abundances of the OTUs were nonetheless strongly correlated. We do not endorse one primer set but rather enumerate strengths and weaknesses to facilitate selection appropriate to a system or experimental goal. While 16S rRNA gene primer bias suggests caution in assessing quantitative population dynamics, community dynamics appear robust across studies using different primers.

Microscale sulfur cycling in the phototrophic pink berry consortia of the Sippewissett Salt Marsh.

Microbial metabolism is the engine that drives global biogeochemical cycles, yet many key transformations are carried out by microbial consortia over short spatiotemporal scales that elude detection by traditional analytical approaches. We investigate syntrophic sulfur cycling in the 'pink berry' consortia of the Sippewissett Salt Marsh through an integrative study at the microbial scale. The pink berries are macroscopic, photosynthetic microbial aggregates composed primarily of two closely associated species: sulfide-oxidizing purple sulfur bacteria (PB-PSB1) and sulfate-reducing bacteria (PB-SRB1). Using metagenomic sequencing and (34) S-enriched sulfate stable isotope probing coupled with nanoSIMS, we demonstrate interspecies transfer of reduced sulfur metabolites from PB-SRB1 to PB-PSB1. The pink berries catalyse net sulfide oxidation and maintain internal sulfide concentrations of 0-500 µm. Sulfide within the berries, captured on silver wires and analysed using secondary ion mass spectrometer, increased in abundance towards the berry interior, while δ(34) S-sulfide decreased from 6‰ to -31‰ from the exterior to interior of the berry. These values correspond to sulfate-sulfide isotopic fractionations (15-53‰) consistent with either sulfate reduction or a mixture of reductive and oxidative metabolisms. Together this combined metagenomic and high-resolution isotopic analysis demonstrates active sulfur cycling at the microscale within well-structured macroscopic consortia consisting of sulfide-oxidizing anoxygenic phototrophs and sulfate-reducing bacteria.

A workflow for genome-wide mapping of archaeal transcription factors with ChIP-seq.

Deciphering the structure of gene regulatory networks across the tree of life remains one of the major challenges in postgenomic biology. We present a novel ChIP-seq workflow for the archaea using the model organism Halobacterium salinarum sp. NRC-1 and demonstrate its application for mapping the genome-wide binding sites of natively expressed transcription factors. This end-to-end pipeline is the first protocol for ChIP-seq in archaea, with methods and tools for each stage from gene tagging to data analysis and biological discovery. Genome-wide binding sites for transcription factors with many binding sites (TfbD) are identified with sensitivity, while retaining specificity in the identification the smaller regulons (bacteriorhodopsin-activator protein). Chromosomal tagging of target proteins with a compact epitope facilitates a standardized and cost-effective workflow that is compatible with high-throughput immunoprecipitation of natively expressed transcription factors. The Pique package, an open-source bioinformatics method, is presented for identification of binding events. Relative to ChIP-Chip and qPCR, this workflow offers a robust catalog of protein-DNA binding events with improved spatial resolution and significantly decreased cost. While this study focuses on the application of ChIP-seq in H. salinarum sp. NRC-1, our workflow can also be adapted for use in other archaea and bacteria with basic genetic tools.

Horizontal gene transfer and CRISPR targeting drive phage-bacterial host interactions and coevolution in pink berry marine microbial aggregates.

Bacteriophages (phages), viruses that infect bacteria, are the most abundant components of microbial communities and play roles in community dynamics and host evolution. The study of phage-host interactions, however, is made difficult by a paucity of model systems from natural environments and known and cultivable phage-host pairs. Here, we investigate phage-host interactions in the "pink berry" consortia, naturally-occurring, low-diversity, macroscopic aggregates of bacteria found in the Sippewissett Salt Marsh (Falmouth, MA, USA). We leverage metagenomic sequence data and a comparative genomics approach to identify eight compete phage genomes, infer their bacterial hosts from host-encoded clustered regularly interspaced short palindromic repeats (CRISPR), and observe the potential evolutionary consequences of these interactions. Seven of the eight phages identified infect the known pink berry symbionts Desulfofustis sp. PB-SRB1, Thiohalocapsa sp. PB-PSB1, and Rhodobacteraceae sp. A2, and belong to entirely novel viral taxa, except for one genome which represents the second member of the Knuthellervirus genus. We further observed increased nucleotide variation over a region of a conserved phage capsid gene that is commonly targeted by host CRISPR systems, suggesting that CRISPRs may drive phage evolution in pink berries. Finally, we identified a predicted phage lysin gene that was horizontally transferred to its bacterial host, potentially via a transposon intermediary, emphasizing the role of phages in bacterial evolution in pink berries. Taken together, our results demonstrate that pink berry consortia contain diverse and variable phages, and provide evidence for phage-host co-evolution via multiple mechanisms in a natural microbial system. IMPORTANCE: Phages (viruses that infect bacteria) are important components of all microbial systems, where they drive the turnover of organic matter by lysing host cells, facilitate horizontal gene transfer (HGT), and co-evolve with their bacterial hosts. Bacteria resist phage infection, which is often costly or lethal, through a diversity of mechanisms. One of these mechanisms are CRISPR systems, which encode arrays of phage-derived sequences from past infections to block subsequent infection with related phages. Here, we investigate bacteria and phage populations from a simple marine microbial community known as "pink berries" found in salt marshes of Falmouth, Massachusetts, as a model of phage-host co-evolution. We identify eight novel phages, and characterize a case of putative CRISPR-driven phage evolution and an instance of HGT between phage and host, together suggesting that phages have large evolutionary impacts in a naturally-occuring microbial community.

A Monte Carlo-based framework enhances the discovery and interpretation of regulatory sequence motifs.

BACKGROUND: Discovery of functionally significant short, statistically overrepresented subsequence patterns (motifs) in a set of sequences is a challenging problem in bioinformatics. Oftentimes, not all sequences in the set contain a motif. These non-motif-containing sequences complicate the algorithmic discovery of motifs. Filtering the non-motif-containing sequences from the larger set of sequences while simultaneously determining the identity of the motif is, therefore, desirable and a non-trivial problem in motif discovery research. RESULTS: We describe MotifCatcher, a framework that extends the sensitivity of existing motif-finding tools by employing random sampling to effectively remove non-motif-containing sequences from the motif search. We developed two implementations of our algorithm; each built around a commonly used motif-finding tool, and applied our algorithm to three diverse chromatin immunoprecipitation (ChIP) data sets. In each case, the motif finder with the MotifCatcher extension demonstrated improved sensitivity over the motif finder alone. Our approach organizes candidate functionally significant discovered motifs into a tree, which allowed us to make additional insights. In all cases, we were able to support our findings with experimental work from the literature. CONCLUSIONS: Our framework demonstrates that additional processing at the sequence entry level can significantly improve the performance of existing motif-finding tools. For each biological data set tested, we were able to propose novel biological hypotheses supported by experimental work from the literature. Specifically, in Escherichia coli, we suggested binding site motifs for 6 non-traditional LexA protein binding sites; in Saccharomyces cerevisiae, we hypothesize 2 disparate mechanisms for novel binding sites of the Cse4p protein; and in Halobacterium sp. NRC-1, we discoverd subtle differences in a general transcription factor (GTF) binding site motif across several data sets. We suggest that small differences in our discovered motif could confer specificity for one or more homologous GTF proteins. We offer a free implementation of the MotifCatcher software package at http://www.bme.ucdavis.edu/facciotti/resources_data/software/.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity.

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the "pink berries" of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

The Green Berry Consortia of the Sippewissett Salt Marsh: Millimeter-Sized Aggregates of Diazotrophic Unicellular Cyanobacteria.

Microbial interactions driving key biogeochemical fluxes often occur within multispecies consortia that form spatially heterogeneous microenvironments. Here, we describe the "green berry" consortia of the Sippewissett salt marsh (Falmouth, MA, United States): millimeter-sized aggregates dominated by an uncultured, diazotrophic unicellular cyanobacterium of the order Chroococcales (termed GB-CYN1). We show that GB-CYN1 is closely related to Crocosphaera watsonii (UCYN-B) and "Candidatus Atelocyanobacterium thalassa" (UCYN-A), two groups of unicellular diazotrophic cyanobacteria that play an important role in marine primary production. Other green berry consortium members include pennate diatoms and putative heterotrophic bacteria from the Alphaproteobacteria and Bacteroidetes. Tight coupling was observed between photosynthetic oxygen production and heterotrophic respiration. When illuminated, the green berries became supersaturated with oxygen. From the metagenome, we observed that GB-CYN1 encodes photosystem II genes and thus has the metabolic potential for oxygen production unlike UCYN-A. In darkness, respiratory activity rapidly depleted oxygen creating anoxia within the aggregates. Metagenomic data revealed a suite of nitrogen fixation genes encoded by GB-CYN1, and nitrogenase activity was confirmed at the whole-aggregate level by acetylene reduction assays. Metagenome reads homologous to marker genes for denitrification were observed and suggest that heterotrophic denitrifiers might co-occur in the green berries, although the physiology and activity of facultative anaerobes in these aggregates remains uncharacterized. Nitrogen fixation in the surface ocean was long thought to be driven by filamentous cyanobacterial aggregates, though recent work has demonstrated the importance of unicellular diazotrophic cyanobacteria (UCYN) from the order Chroococcales. The green berries serve as a useful contrast to studies of open ocean UCYN and may provide a tractable model system to investigate microbial dynamics within phytoplankton aggregates, a phenomenon of global importance to the flux of particulate organic carbon and nitrogen in surface waters.

Candidatus Frankia Datiscae Dg1, the Actinobacterial Microsymbiont of Datisca glomerata, Expresses the Canonical nod Genes nodABC in Symbiosis with Its Host Plant.

Frankia strains are nitrogen-fixing soil actinobacteria that can form root symbioses with actinorhizal plants. Phylogenetically, symbiotic frankiae can be divided into three clusters, and this division also corresponds to host specificity groups. The strains of cluster II which form symbioses with actinorhizal Rosales and Cucurbitales, thus displaying a broad host range, show suprisingly low genetic diversity and to date can not be cultured. The genome of the first representative of this cluster, Candidatus Frankia datiscae Dg1 (Dg1), a microsymbiont of Datisca glomerata, was recently sequenced. A phylogenetic analysis of 50 different housekeeping genes of Dg1 and three published Frankia genomes showed that cluster II is basal among the symbiotic Frankia clusters. Detailed analysis showed that nodules of D. glomerata, independent of the origin of the inoculum, contain several closely related cluster II Frankia operational taxonomic units. Actinorhizal plants and legumes both belong to the nitrogen-fixing plant clade, and bacterial signaling in both groups involves the common symbiotic pathway also used by arbuscular mycorrhizal fungi. However, so far, no molecules resembling rhizobial Nod factors could be isolated from Frankia cultures. Alone among Frankia genomes available to date, the genome of Dg1 contains the canonical nod genes nodA, nodB and nodC known from rhizobia, and these genes are arranged in two operons which are expressed in D. glomerata nodules. Furthermore, Frankia Dg1 nodC was able to partially complement a Rhizobium leguminosarum A34 nodC::Tn5 mutant. Phylogenetic analysis showed that Dg1 Nod proteins are positioned at the root of both α- and ß-rhizobial NodABC proteins. NodA-like acyl transferases were found across the phylum Actinobacteria, but among Proteobacteria only in nodulators. Taken together, our evidence indicates an Actinobacterial origin of rhizobial Nod factors.

Sequencing of seven haloarchaeal genomes reveals patterns of genomic flux.

We report the sequencing of seven genomes from two haloarchaeal genera, Haloferax and Haloarcula. Ease of cultivation and the existence of well-developed genetic and biochemical tools for several diverse haloarchaeal species make haloarchaea a model group for the study of archaeal biology. The unique physiological properties of these organisms also make them good candidates for novel enzyme discovery for biotechnological applications. Seven genomes were sequenced to ∼20×coverage and assembled to an average of 50 contigs (range 5 scaffolds-168 contigs). Comparisons of protein-coding gene compliments revealed large-scale differences in COG functional group enrichment between these genera. Analysis of genes encoding machinery for DNA metabolism reveals genera-specific expansions of the general transcription factor TATA binding protein as well as a history of extensive duplication and horizontal transfer of the proliferating cell nuclear antigen. Insights gained from this study emphasize the importance of haloarchaea for investigation of archaeal biology.

Evaluation of algorithm performance in ChIP-seq peak detection.

Next-generation DNA sequencing coupled with chromatin immunoprecipitation (ChIP-seq) is revolutionizing our ability to interrogate whole genome protein-DNA interactions. Identification of protein binding sites from ChIP-seq data has required novel computational tools, distinct from those used for the analysis of ChIP-Chip experiments. The growing popularity of ChIP-seq spurred the development of many different analytical programs (at last count, we noted 31 open source methods), each with some purported advantage. Given that the literature is dense and empirical benchmarking challenging, selecting an appropriate method for ChIP-seq analysis has become a daunting task. Herein we compare the performance of eleven different peak calling programs on common empirical, transcription factor datasets and measure their sensitivity, accuracy and usability. Our analysis provides an unbiased critical assessment of available technologies, and should assist researchers in choosing a suitable tool for handling ChIP-seq data.