Comprehensive single-PCR 16S and 18S rRNA community analysis validated with mock communities, and estimation of sequencing bias against 18S.

Universal primers for SSU rRNA genes allow profiling of natural communities by simultaneously amplifying templates from Bacteria, Archaea, and Eukaryota in a single PCR reaction. Despite the potential to show relative abundance for all rRNA genes, universal primers are rarely used, due to various concerns including amplicon length variation and its effect on bioinformatic pipelines. We thus developed 16S and 18S rRNA mock communities and a bioinformatic pipeline to validate this approach. Using these mocks, we show that universal primers (515Y/926R) outperformed eukaryote-specific V4 primers in observed versus expected abundance correlations (slope = 0.88 vs. 0.67-0.79), and mock community members with single mismatches to the primer were strongly underestimated (threefold to eightfold). Using field samples, both primers yielded similar 18S beta-diversity patterns (Mantel test, p < 0.001) but differences in relative proportions of many rarer taxa. To test for length biases, we mixed mock communities (16S + 18S) before PCR and found a twofold underestimation of 18S sequences due to sequencing bias. Correcting for the twofold underestimation, we estimate that, in Southern California field samples (1.2-80 µm), there were averages of 35% 18S, 28% chloroplast 16S, and 37% prokaryote 16S rRNA genes. These data demonstrate the potential for universal primers to generate comprehensive microbiome profiles.

Dynamic marine viral infections and major contribution to photosynthetic processes shown by spatiotemporal picoplankton metatranscriptomes.

Viruses provide top-down control on microbial communities, yet their direct study in natural environments was hindered by culture limitations. The advance of bioinformatics enables cultivation-independent study of viruses. Many studies assemble new viral genomes and study viral diversity using marker genes from free viruses. Here we use cellular metatranscriptomics to study active community-wide viral infections. Recruitment to viral contigs allows tracking infection dynamics over time and space. Our assemblies represent viral populations, but appear biased towards low diversity viral taxa. Tracking relatives of published T4-like cyanophages and pelagiphages reveals high genomic continuity. We determine potential hosts by matching dynamics of infection with abundance of particular microbial taxa. Finally, we quantify the relative contribution of cyanobacteria and viruses to photosystem-II psbA (reaction center) expression in our study sites. We show sometimes >50% of all cyanobacterial+viral psbA expression is of viral origin, highlighting the contribution of viruses to photosynthesis and oxygen production.

Dynamics and interactions of highly resolved marine plankton via automated high-frequency sampling.

Short timescale observations are valuable for understanding microbial ecological processes. We assessed dynamics in relative abundance and potential activities by sequencing the small sub-unit ribosomal RNA gene (rRNA gene) and rRNA molecules (rRNA) of Bacteria, Archaea, and Eukaryota once to twice daily between March 2014 and May 2014 from the surface ocean off Catalina Island, California. Typically Ostreococcus, Braarudosphaera, Teleaulax, and Synechococcus dominated phytoplankton sequences (including chloroplasts) while SAR11, Sulfitobacter, and Fluviicola dominated non-phytoplankton Bacteria and Archaea. We observed short-lived increases of diatoms, mostly Pseudo-nitzschia and Chaetoceros, with quickly responding Bacteria and Archaea including Flavobacteriaceae (Polaribacter & Formosa), Roseovarius, and Euryarchaeota (MGII), notably the exact amplicon sequence variants we observed responding similarly to another diatom bloom nearby, 3 years prior. We observed correlations representing known interactions among abundant phytoplankton rRNA sequences, demonstrating the biogeochemical and ecological relevance of such interactions: (1) The kleptochloroplastidic ciliate Mesodinium 18S rRNA gene sequences and a single Teleaulax taxon (via 16S rRNA gene sequences) were correlated (Spearman r = 0.83) yet uncorrelated to a Teleaulax 18S rRNA gene OTU, or any other taxon (consistent with a kleptochloroplastidic or karyokleptic relationship) and (2) the photosynthetic prymnesiophyte Braarudosphaera bigelowii and two strains of diazotrophic cyanobacterium UCYN-A were correlated and each taxon was also correlated to other taxa, including B. bigelowii to a verrucomicrobium and a dictyochophyte phytoplankter (all r > 0.8). We also report strong correlations (r > 0.7) between various ciliates, bacteria, and phytoplankton, suggesting interactions via currently unknown mechanisms. These data reiterate the utility of high-frequency time series to show rapid microbial reactions to stimuli, and provide new information about in situ dynamics of previously recognized and hypothesized interactions.

Microbial Mat Compositional and Functional Sensitivity to Environmental Disturbance.

The ability of ecosystems to adapt to environmental perturbations depends on the duration and intensity of change and the overall biological diversity of the system. While studies have indicated that rare microbial taxa may provide a biological reservoir that supports long-term ecosystem stability, how this dynamic population is influenced by environmental parameters remains unclear. In this study, a microbial mat ecosystem located on San Salvador Island, The Bahamas was used as a model to examine how environmental disturbance affects the protein synthesis potential (PSP) of rare and abundant archaeal and bacterial communities and how these changes impact potential biogeochemical processes. This ecosystem experienced a large shift in salinity (230 to 65 g kg-1) during 2011-2012 following the landfall of Hurricane Irene on San Salvador Island. High throughput sequencing and analysis of 16S rRNA and rRNA genes from samples before and after the pulse disturbance showed significant changes in the diversity and PSP of abundant and rare taxa, suggesting overall compositional and functional sensitivity to environmental change. In both archaeal and bacterial communities, while the majority of taxa showed low PSP across conditions, the overall community PSP increased post-disturbance, with significant shifts occurring among abundant and rare taxa across and within phyla. Broadly, following the post-disturbance reduction in salinity, taxa within Halobacteria decreased while those within Crenarchaeota, Thaumarchaeota, Thermoplasmata, Cyanobacteria, and Proteobacteria, increased in abundance and PSP. Quantitative PCR of genes and transcripts involved in nitrogen and sulfur cycling showed concomitant shifts in biogeochemical cycling potential. Post-disturbance conditions increased the expression of genes involved in N-fixation, nitrification, denitrification, and sulfate reduction. Together, our findings show complex community adaptation to environmental change and help elucidate factors connecting disturbance, biodiversity, and ecosystem function that may enhance ecosystem models.

Long-term operation of microbial electrosynthesis systems improves acetate production by autotrophic microbiomes.

Microbial electrosynthesis is the biocathode-driven production of chemicals from CO2 and has the promise to be a sustainable, carbon-consuming technology. To date, microbial electrosynthesis of acetate, the first step in order to generate liquid fuels from CO2, has been characterized by low rates and yields. To improve performance, a previously established acetogenic biocathode was operated in semi-batch mode at a poised potential of -590 mV vs SHE for over 150 days beyond its initial development. Rates of acetate production reached a maximum of 17.25 mM day(-1) (1.04 g L(-1) d(-1)) with accumulation to 175 mM (10.5 g L(-1)) over 20 days. Hydrogen was also produced at high rates by the biocathode, reaching 100 mM d(-1) (0.2 g L(-1) d(-1)) and a total accumulation of 1164 mM (2.4 g L(-1)) over 20 days. Phylogenetic analysis of the active electrosynthetic microbiome revealed a similar community structure to what was observed during an earlier stage of development of the electroacetogenic microbiome. Acetobacterium spp. dominated the active microbial population on the cathodes. Also prevalent were Sulfurospirillum spp. and an unclassified Rhodobacteraceae. Taken together, these results demonstrate the stability, resilience, and improved performance of electrosynthetic biocathodes following long-term operation. Furthermore, sustained product formation at faster rates by a carbon-capturing microbiome is a key milestone addressed in this study that advances microbial electrosynthesis systems toward commercialization.

Microbial phylogenetic profiling with the Pacific Biosciences sequencing platform.

High-throughput sequencing of 16S rRNA gene amplicons has revolutionized the capacity and depth of microbial community profiling. Several sequencing platforms are available, but most phylogenetic studies are performed on the 454-pyrosequencing platform because its longer reads can give finer phylogenetic resolution. The Pacific Biosciences (PacBio) sequencing platform is significantly less expensive per run, does not rely on amplification for library generation, and generates reads that are, on average, four times longer than those from 454 (C2 chemistry), but the resulting high error rates appear to preclude its use in phylogenetic profiling. Recently, however, the PacBio platform was used to characterize four electrosynthetic microbiomes to the genus-level for less than USD 1,000 through the use of PacBio's circular consensus sequence technology. Here, we describe in greater detail: 1) the output from successful 16S rRNA gene amplicon profiling with PacBio, 2) how the analysis was contingent upon several alterations to standard bioinformatic quality control workflows, and 3) the advantages and disadvantages of using the PacBio platform for community profiling.

Electrosynthesis of commodity chemicals by an autotrophic microbial community.

A microbial community originating from brewery waste produced methane, acetate, and hydrogen when selected on a granular graphite cathode poised at -590 mV versus the standard hydrogen electrode (SHE) with CO(2) as the only carbon source. This is the first report on the simultaneous electrosynthesis of these commodity chemicals and the first description of electroacetogenesis by a microbial community. Deep sequencing of the active community 16S rRNA revealed a dynamic microbial community composed of an invariant Archaea population of Methanobacterium spp. and a shifting Bacteria population. Acetobacterium spp. were the most abundant Bacteria on the cathode when acetogenesis dominated. Methane was generally the dominant product with rates increasing from <1 to 7 mM day(-1) (per cathode liquid volume) and was concomitantly produced with acetate and hydrogen. Acetogenesis increased to >4 mM day(-1) (accumulated to 28.5 mM over 12 days), and methanogenesis ceased following the addition of 2-bromoethanesulfonic acid. Traces of hydrogen accumulated during initial selection and subsequently accelerated to >11 mM day(-1) (versus 0.045 mM day(-1) abiotic production). The hypothesis of electrosynthetic biocatalysis occurring at the microbe-electrode interface was supported by a catalytic wave (midpoint potential of -460 mV versus SHE) in cyclic voltammetry scans of the biocathode, the lack of redox active components in the medium, and the generation of comparatively high amounts of products (even after medium exchange). In addition, the volumetric production rates of these three commodity chemicals are marked improvements for electrosynthesis, advancing the process toward economic feasibility.

Characterization and quantitation of a novel ß-lactamase gene found in a wastewater treatment facility and the surrounding coastal ecosystem.

Wastewater treatment plants (WWTPs) are engineered structures that collect, concentrate, and treat human waste, ultimately releasing treated wastewater into local environments. While WWTPs efficiently remove most biosolids, it has been shown that many antibiotics and antibiotic-resistant bacteria can survive the treatment process. To determine how WWTPs influence the concentration and dissemination of antibiotic-resistant genes into the environment, a functional metagenomic approach was used to identify a novel antibiotic resistance gene within a WWTP, and quantitative PCR (qPCR) was used to determine gene copy numbers within the facility and the local coastal ecosystem. From the WWTP metagenomic library, the fosmid insert contained in one highly resistant clone (MIC, ≈ 416 µg ml(-1) ampicillin) was sequenced and annotated, revealing 33 putative genes, including a 927-bp gene that is 42% identical to a functionally characterized ß-lactamase from Staphylococcus aureus PC1. Isolation and subcloning of this gene, referred to as bla(M-1), conferred ampicillin resistance to its Escherichia coli host. When normalized to volume, qPCR showed increased concentrations of bla(M-1) during initial treatment stages but 2-fold-decreased concentrations during the final treatment stage. The concentration ng(-1) DNA increased throughout the WWTP process from influent to effluent, suggesting that bla(M-1) makes up a significant proportion of the overall genetic material being released into the coastal ecosystem. Average discharge was estimated to be 3.9 × 10(14) copies of the bla(M-1) gene released daily into this coastal ecosystem. Furthermore, the gene was observed in all sampled coastal water and sediment samples surrounding the facility. Our results suggest that WWTPs may be a pathway for the dissemination of novel antibiotic resistance genes into the environment.

Extraction of high molecular weight DNA from microbial mats.

Successful and accurate analysis and interpretation of metagenomic data is dependent upon the efficient extraction of high-quality, high molecular weight (HMW) community DNA. However, environmental mat samples often pose difficulties to obtaining large concentrations of high-quality, HMW DNA. Hypersaline microbial mats contain high amounts of extracellular polymeric substances (EPS)1 and salts that may inhibit downstream applications of extracted DNA. Direct and harsh methods are often used in DNA extraction from refractory samples. These methods are typically used because the EPS in mats, an adhesive matrix, binds DNA during direct lysis. As a result of harsher extraction methods, DNA becomes fragmented into small sizes. The DNA thus becomes inappropriate for large-insert vector cloning. In order to circumvent these limitations, we report an improved methodology to extract HMW DNA of good quality and quantity from hypersaline microbial mats. We employed an indirect method involving the separation of microbial cells from the background mat matrix through blending and differential centrifugation. A combination of mechanical and chemical procedures was used to extract and purify DNA from the extracted microbial cells. Our protocol yields approximately 2 µg of HMW DNA (35-50 kb) per gram of mat sample, with an A(260/280) ratio of 1.6. Furthermore, amplification of 16S rRNA genes suggests that the protocol is able to minimize or eliminate any inhibitory effects of contaminants. Our results provide an appropriate methodology for the extraction of HMW DNA from microbial mats for functional metagenomic studies and may be applicable to other environmental samples from which DNA extraction is challenging.

Extraction of high molecular weight DNA from microbial mats.

Due to the presence of inhibitors such as extracellular polymeric substances (EPSs) and salts, most microbial mat studies have relied on harsh methods of direct DNA extraction that result in DNA fragments too small for large-insert vector cloning. High molecular weight (HMW) DNA is crucial in functional metagenomic studies, because large fragments present greater access to genes of interest. Here we report improved methodologies for extracting HMW DNA from EPS-rich hypersaline microbial mats. The protocol uses a combination of microbial cell separation with mechanical and chemical methods for DNA extraction and purification followed by precipitation with polyethylene glycol (PEG). The protocol yields >2 µg HMW DNA (>48 kb) per gram of mat sample, with A260:280 ratios >1.7. In addition, 16S rRNA gene analysis using denaturing gradient gel electrophoresis and pyrosequencing showed that this protocol extracts representative DNA from microbial mat communities and results in higher overall calculated diversity indices compared with three other standard methods of DNA extraction. Our results show the importance of validating the DNA extraction methods used in metagenomic studies to ensure optimal recovery of microbial richness.