A Sample-to-Sequence Protocol for Genus Targeted Transcriptomic Profiling: Application to Marine Synechococcus
- 1School of Life Sciences, University of Warwick, Coventry, UK
- 2Warwick Medical School, University of Warwick, Coventry, UK
- 3Department of Biochemistry and Biomolecular Science, Faculty of Science and Engineering, Macquarie University, Sydney, NSW, Australia
- 4Sydney Cytometry, Centenary Institute, University of Sydney, Sydney, NSW, Australia
- 5National Oceanography Centre, Southampton, UK
Recent studies using whole community metagenomic and metatranscriptomic approaches are revealing important new insights into the functional potential and activity of natural marine microbial communities. Here, we complement these approaches by describing a complete ocean sample-to-sequence protocol, specifically designed to target a single bacterial genus for purposes of both DNA and RNA profiling using fluorescence activated cell sorting (FACS). The importance of defining and understanding the effects of a sampling protocol are critical if we are to gain meaningful data from environmental surveys. Rigorous pipeline trials with a cultured isolate, Synechococcus sp. BL107 demonstrate that water filtration has a well-defined but limited impact on the transcriptomic profile of this organism, whilst sample storage and multiple rounds of cell sorting have almost no effect on the resulting RNA sequence profile. Attractively, the means to replicate the sampling strategy is within the budget and expertise of most researchers.
Next generation sequencing (NGS) techniques offer huge potential for the genetic exploration of marine bacteria. Initial Sanger-based genome sequencing of marine isolates (Palenik et al., 2003; Rocap et al., 2003; Giovannoni et al., 2005) has been superseded by NGS approaches capable of sequencing multiple genomes simultaneously (Kashtan et al., 2014; Billerbeck et al., 2016). Corresponding marker gene studies have revealed the spatial partitioning and temporal succession of marine bacterial communities (Fuhrman and Steele, 2008; Gilbert et al., 2012; Teeling et al., 2012; Needham and Fuhrman, 2016), highlighting the diversity and geographic distribution of specific marine microbes (Johnson et al., 2006; Zwirglmaier et al., 2008; Gómez-Pereira et al., 2010; Farrant et al., 2016), exposing the extent and dynamics of the rare microbial biosphere (Campbell et al., 2011; Welch and Huse, 2011). Advances in cell sorting and single cell sequencing circumvent the limitations of culturing to study the genomes of keystone marine species (Zehr et al., 2008; Tripp et al., 2010; Heywood et al., 2011; Worden et al., 2012; Kashtan et al., 2014; Mende et al., 2016) while the greater sequencing depth achieved by the latest NGS platforms unlocks the potential for detailed metagenomic and metatranscriptomic studies of natural microbial communities (Cuvelier et al., 2010; Shi et al., 2010; Batmalle et al., 2014; Satinsky et al., 2014; Aylward et al., 2015; Shilova et al., 2016). Together with large-scale ocean sampling surveys (Nicholls, 2007; Hurwitz and Sullivan, 2013; Armbrust and Palumbi, 2015; Rees et al., 2015; Sunagawa et al., 2015) and well-established ocean time stations (Phillips and Joyce, 2007; Karl and Church, 2014) the means to tap the oceans entire genomic potential is seemingly within our grasp.
Environmental sampling approaches for both metagenomic and metatranscriptomic studies are clearly shaped by the questions researchers wish to ask yet constrained by an understanding of current limits in technology. The task of assessing the genomic potential of the marine picoplankton community has led to much of the early pioneering work going into the development and validation of sampling strategies and methodologies. In all cases, a specific size fraction from the water sample is filtered onto a flat membrane or pleated filter, which is preserved for DNA extraction and sequenced at a later date. The introduction of bias into the profile can be observed with the choice of cell size fractionation (Craig, 1986; Padilla et al., 2015), the method in which the cells are lysed and the manner in which material is amplified–if yield requirements are not met for sequencing (Batmalle et al., 2014). This entire process becomes even more complex when we consider using the same material for combined DNA and RNA studies. Since the half-life stability of bacterial mRNA is on the order of minutes (Steglich et al., 2010) the amount of time taken to collect and concentrate biomass from the environment is a critical factor. However, reducing the time for collection and concentration limits the total yield of biomass, with the consequence that the obtained material must then be amplified to meet DNA/RNA library preparation requirements. Important considerations need to be made, not only in the time and manner in which the material is collected, filtered and preserved but also the handling of the RNA post extraction—all of which have the potential to elicit a transcriptional response and impact on the recovered RNA profile. Considerable effort in recent years has gone into the development of environmental transcriptomic methodologies, with a view to understanding how the sampling process impacts the data produced (Stewart et al., 2010; Ottesen et al., 2011).
While large-scale marine metagenomics projects are beginning to reveal the genomic potential of distinct communities, there are still major hurdles to overcome in understanding functional adaptation. Transcriptomics provides the means to address two of the current major challenges in marine microbiology: (i) to develop an understanding of how the genomic blueprint of an organism is dynamically expressed under a set of environmental conditions, and (ii) to begin to uncover the role of the plethora of hypothetical genes i.e., those with no known function, present in all microbial genomes/metagenomes. Correlating gene expression levels with changing environmental conditions provides a means to group known and unknown genes into functional regulons. However, the ability to assign functional traits to a specific genus or group within the community is underpinned by an ability to map reads accurately to the correct taxonomic group. Community-wide analyses of sequence data may provide an overview of biological functions in organisms that dominate DNA but fail to assign genes and associated traits to keystone species involved in critical environmental processes, such as primary production and nitrogen fixation (Jardillier et al., 2010). This is of particular importance when attempting to dissect the functional potential of genes from marine microbes and to link molecular adaptations with the eco-physiological limits of species and ecotypes.
Genome sequencing of both laboratory isolates and single cells has built a mosaic picture of core genes common to many, and auxiliary genes specific to a few. The auxiliary gene sets present some of the most interesting, perhaps niche-defining features of the genomes. However, their presence on mobile genomic islands makes the taxonomic assignment of reads for some genes potentially problematic since their phylogenetic origins may lie in an entirely different genus (Coleman et al., 2006; Dufresne et al., 2008; Penn et al., 2009). With this in mind, the selective isolation of a single genus from an environmental sample, made possible by fluorescence activated cell sorting (FACS), could lead to a more accurate assignment of transcriptomic reads. Targeting a specific planktonic group via FACS for metagenomic sequencing has already been shown to improve specific gene coverage, particularly important if the group under study is a minority player within the whole community (Cuvelier et al., 2010; Thompson et al., 2013; Batmalle et al., 2014). However, there remain significant challenges in both RNA recovery and amplification bias before FACS targeted transcriptomics can be realized. To address these challenges, we conducted a series of experiments designed to validate a complete ocean-to-sequence protocol with a proof of concept application on marine Synechococcus.
Cyanobacteria of the genus Synechococcus are a ubiquitous component of the marine photosynthetic community, significantly contributing to primary productivity (Jardillier et al., 2010; Flombaum et al., 2013). Whilst numerically less abundant than their sister taxon Prochlorococcus, they have a wider geographic distribution occupying diverse marine habitats from large mid-ocean oligotrophic gyres to nutrient-enriched coastal waters, upwelling systems and polar regions (Li et al., 1983; Partensky et al., 1999; Saito et al., 2005; Cottrell and Kirchman, 2009). Through genome sequencing, phylogenetic marker amplicon analysis and metagenomic recruitment studies, a clear picture of the genus is emerging where the existence of multiple ecotypes inhabiting distinct ocean niches prevails (Dufresne et al., 2008; Zwirglmaier et al., 2008; Farrant et al., 2016; Sohm et al., 2016). Factors controlling Synechococcus abundance and activity have been explored in the form of in situ rate measurements and candidate marker genes representing key processes (Fuller et al., 2005; Lindell et al., 2005; Grob et al., 2015). The obvious limitation of this approach is that it draws on known expression analyses of cultured representatives and relies on the assumption that cells perform in a similar manner in the environment. A targeted in situ transcriptomic approach would ultimately aim to identify clade-specific gene sets under dynamic expression, to interpret what genomic factors are most important in controlling both geographical and temporal distribution patterns and contributing to their activity in global nutrient cycles.
The key challenges for a sampling protocol are to process sufficient sample volume over a short timeframe to yield enough material for amplification-free sequencing without significantly altering the in situ transcriptional profile of the target cells. In the nutrient deplete oligotrophic gyres, Synechococcus cell numbers can fall to as low as 1 × 103 cells ml−1 (Partensky et al., 1999; Zubkov et al., 2007; Zwirglmaier et al., 2008) which can represent < 1% of the cells in the entire community. Whereas, 1–20 l of seawater would provide adequate material for a whole microbial community-level RNA extraction, far larger (>10-fold), volumes of water would be required to yield sufficient flow sorted Synechococcus cells for the extraction of both DNA and RNA. In addition, the combined effects of both concentration and cell sorting on a transcription profile are poorly described and need to be considered carefully in order to draw meaningful biological information from expression data.
To address these issues we designed a new high-speed filtration and cell concentration apparatus using accessible and relatively low-cost components and validated the protocol for transcript recovery from a simulated open ocean filtration and sorting experiment against that of a control using a single reference strain of Synechococcus (BL107). The transcription profile of BL107 was examined at each stage of the protocol and compared with a control profile where no processing had occurred. We were able to demonstrate that (i) all major changes to the transcript profile occur within the filtration phase and not the cell sorting phase (ii) changes to the transcript profile are limited and well-defined, and (iii) rapid cell filtration can be achieved via a low cost, easily maintained filtration rig that does not require specialist operation.
Materials and Methods
High Speed Water Filtration (SWiFT) Apparatus Design and Operation
The SWiFT cell capture rig (Figure 1) is designed to take seawater directly from source. Seawater undergoes a two-phase filtration: first passing through a 35 μm cartridge filter housed within a standard consumer grade under-counter water filter unit before passing through a Memteq High Volume Celltrap (HVCT) 0.2 μm filter unit. Filter units are connected via silicone tubing attached with jubilee clips to the barb connectors. Water is driven through the apparatus via a Watson Marlow 620SN/R cased peristaltic pump. An in-line pressure gage allows the user to monitor pressure increase as cells build up within the HVCT, avoiding HVCT overloading. The SWiFT rig is adaptable for use with alternative peristaltic pump models, subject to their liters/min water pumping capacity. HVCT units are pre-packed sterile. Following a single round of filtration and cell recovery the unit is disposed and not re-used. Before entering the HVCT water enters a large transparent housing unit which, once filled with flow through water, acts as a bubble trap drawing water from the bottom while inflow occurs at the top. Avoidance of air entering the HVCT is a critical step as once a bacterial film begins to line the HVCT filaments, bubbles are unable to pass through with the resulting accumulation of bubbles blocking flow through the HVCT over time.
Figure 1. SWiFT cell capture rig technical specifications. A Watson Marlow 620SN/R encased pump draws source water through the apparatus, water flow indicated by blue arrows. Clear water filtration unit (A) 10″ houses the disposable supaplete II 35 micron pre-filter (Amazon filters) while (B) 20″ acts as a bubble trap once air is purged from the chamber via the adjustable valve. High volume Celltrap units connect to the main body via luerlok connection to perform a dead-end filtration. (C) Cells are recovered from the disconnected CT by drawing filtrate back through the filaments with a 50 ml luer lock syringe (VWR).
Synechococcus sp. BL107 Sampling
A mid-exponential phase culture of Synechococcus sp. BL107 grown in artificial seawater medium (ASW; Wyman et al., 1985) at 22°C with 10 μmol photons m−2 s−1 continuous light was divided in two. Half of the culture was centrifuged at 4754 × g for 4 min, and three RNA extractions performed immediately. The remainder of the culture was diluted in ASW to represent a concentration of Synechococcus more akin to those found in natural environments of the North Atlantic (ca. 3 × 104 cells ml−1; Zubkov et al., 1998) and filtered using the SWiFT rig (Figure 1) for 15 min. Cells were recovered in suspension from the HVCT unit by drawing back through the trap 10 ml filtrate with a syringe followed immediately by flash freezing in liquid nitrogen. Total filtration time inclusive of water transfer to rig and cell freezing amounted to 20 min. Frozen samples were stored at −80°C for 16 weeks to replicate typical on-ship sample transit times following the Atlantic Meridional Transect ocean sampling program.
Flow Cytometric Sorting of Synechococcus sp. BL107
Synechococcus sp. BL107 cells were recovered from the frozen filtrate by FACS using a BD Influx™ Cell Sorter (BD Biosciences). Before each round of sorting both sheath line and sample line were cleaned with two rounds of FACSclean and FACSrinse each taking 30 min to ensure purity of both the sample line and sheath fluid. Filter sterilized (0.22 μm) ASW without phosphate/nitrate, and trace metals was used as sheath fluid. A 70 μm nozzle was used in combination with a fluidic pressure of 65 psi during sorting. For the single round of sorting, 0.5 ml of frozen filtrate was defrosted on ice for 10 min in the presence of potassium citrate (30 mM final concentration) and SYBRGold nucleic acid stain, final concentration 10−2 concentration of the commercial stock (Invitrogen; Marie et al., 1997). Cells were excited using a 488 nm, 200 mW laser with a standard filter setup. The BD Influx™ was set to one-drop pure sort mode (highest sort purity for this instrument) and Synechococcus sp. BL107 cells gated and separated based on their green (530 nm) and orange (580 nm) fluorescence. Cells leaving the sort stream were diverted into a custom collection chamber designed to support a Falcon 5 ml polypropylene round bottom flow tube within a bath of dry ice and ethanol (Figure 2). Cells were instantly frozen on contact with the collection tube preventing prolonged thaw times. Each thawed and stained aliquot of cells took 10 min for sort completion. A total of 2.5 × 107 cells were sorted for each of the three single sort samples. For the two consecutive rounds of cell sorting, the first round served as a basic enrichment. Cells were not stained, with sort gating and selection solely based on their phycoerythrin auto-fluorescence (580 nm). In the second round of flow sorting cells were stained and sorted as described above. A total of 1.5 × 107 cells were flow-sorted for the two double sort samples. Sort purity was regularly assessed as described by Zubkov et al. (2007). Sorted material was 99% enriched with recovery ranging from 85 to 95%.
Figure 2. Custom BD Influx flow tube holder cross-section. The tube holder was designed to support a two-way sort program. Flow tubes are suspended within a hollow chamber designed to bring them into direct contact with a dry ice ethanol bath. Sorted cells freeze on contact with inserted flow tubes.
Flow Cytometric Counts of Synechococcus sp. BL107
Synechococcus sp. BL107 cell counts for the purpose of pre-filter retention testing were performed on a FACScan flow cytometer (BD Biosciences). Cells were excited using a 488 nm 20 mW laser with detection based on 580 nm auto-fluorescence properties (Olson et al., 1993). To calculate cells ml−1, the material recovered through HVCT was diluted 1000-fold in ASW medium spiked with a known concentration of fluorescent beads (Sphero Nile Red Beads 1.7–2.2 μm). From the bead counts the total number of beads detected in a fixed amount of time was calculated as a percentage of the total known bead concentration. Cell concentration could then be calculated based on the assumption that the same percentage of cells pass through the laser compared to beads.
RNA Extraction and Next Generation Sequencing
RNA extraction was based on a variation of the protocol described in Gilbert et al. (2000). Briefly, frozen sorted cells were filtered onto a 13 mm 0.2 μm Supor® filter (Pall) and re-suspended in 0.4 ml solution A (0.3 M sucrose, 0.01 M sodium acetate), 0.4 ml solution B (0.01 M sodium acetate pH 4.2 2% SDS), 0.8 ml acidic phenol (pH 4.3) then bead beaten in Lysing Matrix E tubes (MP Biomedicals) for 45 s, transferred to ice for 45 s and repeated. The aqueous phase was then purified using the RNA Clean & Concentrator™ −5 kit (Zymo Research, CA) with a final elution volume of 15 μl. Samples were then treated with Turbo DNase (Ambion) for 30 min followed by a second round of purification with the RNA Clean & Concentrator™ −5 kit. Total RNA quality and yield were assessed using Bioanalyzer RNA 6000 nano and pico chips (Agilent) and QuantiFluor RNA System (Promega) analysis, respectively. Total RNA was prepared for sequencing using the ScriptSeq v2 (Epicentre) library workflow for the HiSeq2000 2 × 100 bp reads (V3 chemistry). Samples were sequenced to various depths by the NERC Biomolecular Analysis Facility (NBAF) located at the Centre for Genomic Research (University of Liverpool).
Analysis of RNAseq Read Data
Illumina reads (ArrayExpress accession E-MTAB-4814) were quality trimmed and mapped to the Cyanorak curated Synechococcus sp. BL107 genome (www.sb-roscoff.fr/cyanorak/) using CLC Genomics Workbench 8.0.1 (Qiagen). Briefly, reads were quality trimmed before global read mapping to the Synechococcus sp. BL107 genome with the following altered mapping stringencies: mismatch penalty 3, deletion cost 3, length fraction 0.9, similarity fraction 0.9. Counting of both paired and singlet reads that mapped to protein coding regions only gave total reads per gene. The reads mapping over a gene boundary were counted toward the total expression value for the gene they best mapped to. Gene counts for each sample were then normalized and compared to the control (RNA extracted from the material that had not been filtered or sorted) using the DESeq2 package in R (Love et al., 2014). The four copies of the psbA gene were excluded from the analysis due to their high degree of sequence similarity (95%) at the nucleotide level. Combined with a short transcript read length and the nature by which reads are assigned when they map equally well to multiple regions, it is not possible to accurately distinguish changes in transcript ratios between the four copies (Garczarek et al., 2008).
SWiFT: An Effective Filtration Rig for the Rapid Recovery of Marine Picoplankton
SWiFT (Figure 1) is a unique three-stage water filtration unit designed to concentrate and recover picoplankton in liquid from large volumes of seawater in a short space of time. Seawater is propelled through the system via a peristaltic pump, undergoing an initial filtration to remove large biomass >35 μm, a second stage which removes air bubbles from the system and finally a third stage in which cells (>0.2 μm) are captured via a parallel 0.22 μm diameter pore size hollow fiber filter (CTHV400, Celltrap™, MEMTEQ Ventures Ltd. UK).
Initial cell sorting trials with Synechococcus sp. BL107 revealed a minimum of 1.5 × 107 cells were required for sufficient RNA recovery to meet minimum total RNA yields for the ScriptSeq v2 (Epicentre) library preparation without the need for amplification. Based on the lowest measurements for open ocean Synechococcus cell numbers 1 × 103 cells ml−1 (Partensky et al., 1999; Zubkov et al., 2007; Zwirglmaier et al., 2008) and an approximation of 93% cell recovery (Table 1) we estimated that 200 l of seawater would need to be processed within a 20 min time frame, inclusive of water transfer times and cell freezing.
While transcripts turning over at the timescale of a few minutes are not captured by most current sampling devices i.e., except for autosamplers (Jacquet et al., 1998), a 20 min window was deemed effective for capturing the majority of transcripts with half-lives in the midrange based on a previous study of the RNA half-life from Prochlorococcus (Steglich et al., 2010). Flat membrane filters commonly used for bacterial community recovery were not able to handle the biomass load or speed of filtration. A significant advantage of the hollow fiber filter is that concentrated cells can be efficiently recovered in ~20 s in liquid, using a syringe, whereas recovering the same biomass from a large flat filter would be considerably slower with a lower efficiency.
Celltraps that consist of parallel hollow fiber filter tubing designed for dead end filtration (Figure 1C) have been routinely used for cell recovery in liquid (Gómez-Pereira et al., 2013). While the CT40 model is still restrictive with a 2–20 l seawater capacity, the new high volume CT400, with an estimated equivalent surface area of 100,000 mm2 is now able to meet both volume and time restrictions processing up to 300 l within 20 min when used in conjunction with a 620SN/R pump (Watson Marlow) and a 15.9 mm bore silicone tubing.
We used a pre-filtration step to remove particles larger than 35 μm through a large capacity cartridge pleat filter (Amazon Filters Ltd). Industrial pleat filters are graded to prevent passage of particles above a target size range. However, smaller particle retention rates were unknown. Both pleat and flat membrane filtration trials were carried out using cultured Synechococcus sp. BL107 cells to assess recovery rates. Data showed that the Supapleat II (35 μm pore graded) pleat filter gave comparable cell recovery rates to a Nuclepore flat membrane filtration system (Table 1).
Changes to the Synechococcus Transcription Profile Caused by Filtration were Limited and Well-Defined
A trial protocol was designed to isolate and assess the impact of a 20-min HVCT cell concentration, combined with a frozen sample transit period, followed by the cumulative impact of FACS on the transcriptome of Synechococcus sp. BL107. The RNA yield following extraction was low (< 35 ng), which ruled out the possibility to deplete the rRNA pool by subtractive hybridisation, with the risk of biasing the resultant mRNA pool (Giannoukos et al., 2012). Total RNA was therefore sequenced with no prior amplification on the Illumina Hi-Seq platform (Supplementary Table 1).
Maximum and average total read counts for each of the conditions demonstrated comprehensive coverage of the Synechococcus sp. BL107 transcriptome with 6% of genes (150) on average across all conditions reporting no counts (Table 2). Comparison of the SWiFT filtered with the unfiltered Synechococcus sp. BL107 transcriptome showed 9% (226) of protein encoding genes demonstrated a >2-fold change (q < 0.05) in their expression levels (Figure 3). This is a relatively limited gene set that is being affected, indeed, when compared to nutrient starvation (e.g., iron stress in Synechococcus sp. BL107 causes the differential expression of 15% of genes (Millard, Eriksson and Scanlan unpublished data).
Figure 3. Volcano plot demonstrating significant changes to the Synechococcus sp. BL107 transcriptome post SWiFT filtration. Genes demonstrating < 0.05 q-value and >2-fold change in gene expression are highlighted in blue.
Close to half (48%, 107) of the genes significantly altered had known annotations within the Kyoto Encyclopaedia of Genes and Genomes (KEGG) allowing further categorisation into 15 broad functional groups (Table 3). Five subgroups were found to contain a significant enrichment of genes known to be associated with a particular pathway or functional cluster (Fisher's exact test, Table 4). These groups were: photosynthetic antenna, photosynthetic core, phycobilin metabolism, chaperones, and nitrogen metabolism associated genes. All of the above groups demonstrated a decrease in expression with the exception of the photosynthetic core associated ferredoxin (petF), circadian regulators rpaA, and primary sigmaF2, and all of the chaperone proteins.
Table 3. Functional groups defining the effect of SWiFT filtration on the transcriptome profile of Synechococcus sp. BL107.
Changes to the gene expression profile appear well-defined when put in the context of the filtration process. Transfer of Synechococcus sp. BL107 through the SWiFT rig involves pumping the water through the apparatus at pressure. Water pressure entering the HVCT increases gradually (up to 25 psi) as the bacterial load builds up into a film on the cell trap restricting water flow. This would account for the increase in chaperone and heat shock associated proteins observed in the profile, corresponding to the bacterial response to moderate increases in hydrostatic pressure (Welch et al., 1993). Cells traveling through the filtration rig enter the cell trap filaments, which effectively place them in darkness for a period of up to 20 min depending on the time they spend within the filament core. The movement of cells into the dark filament cores has by far the largest impact on the transcriptome of Synechococcus sp. BL107. Indicators that this transition could be occurring can be seen in the decrease in expression of rpaA and increase in rpaB (Table 3). These two key response regulators are known to function cooperatively as part of the circadian kaiC-associated output pathway (Hanaoka et al., 2012; Markson et al., 2013) and have been shown to demonstrate a similar day to night transcriptome profile in the closely related cyanobacterium Prochlorococcus sp. MED4 (Zinser et al., 2009).
While no change was observed in the kaiABC core circadian oscillating genes, the transition to dark is strongly reflected in the changes in the rest of the transcriptome profile. Most notably down-regulation of the photosynthetic antenna and core photosystem genes occurs in conjunction with the down-regulation of the nitrogen transport and metabolism pathways, all common features of the switch from day to night (Waldbauer et al., 2012). Notable absences from the profile include changes to genes associated with glycogen metabolism, pentose phosphate pathway and the Calvin cycle. These genes together play important roles in orchestrating the cells' shift to the major night time activity of respiratory metabolism (Diamond et al., 2015). Failure to detect all common genes involved in the day to night transition is perhaps to be expected as typical circadian gene levels oscillate over a 24-h period whereas the BL107 cells only experience a dark period up to 20 min during SWiFT filtration.
Flow Sorting has Minimal Effect on the Synechococcus sp. BL107 Transcriptome Profile
Following both a single and double round of FACS only 5 and 6 genes, respectively changed their expression levels (q < 0.05, fold change >2) when compared with the post SWiFT filtration transcriptome profile (Table 5) The fact that changes to the transcription profile are extremely limited provides good evidence that the rapid thaw and cell sorting periods effectively preserve the transcript profile post SWiFT filtration at the point of snap freezing.
Table 5. Genes demonstrating a change in expression post SWiFT filtration following both a single and double round of flow cytometric sorting.
To our knowledge, this is the first study in which sorted cells are frozen immediately after drop deflection from the main sort stream as a result of a custom designed cooling chamber (Figure 2). This effectively reduces total cell thaw time by up to two-thirds compared with a standard 20-min workflow for a 1 ml sample from thaw to final sort collection. Prior to FACS, cells were stored at −80°C for 16 weeks to replicate sample transit times during long cruise sampling transects. While RNA yield comparisons between the control and sorted samples are complicated by the fact that the cells were treated differently prior to RNA extraction, Bioanalyzer analysis demonstrates RNA integrity remains constant following −80°C storage for 16-weeks (Figure 4). The observed limited changes in the transcript profile also suggest RNA degradation during extended periods of freezing is limited, at least in the case of Synechococcus sp. BL107.
Figure 4. Bioanalyzer electrophoretogram output for prokaryotic total RNA fragment analysis. (A) Control samples: cells which underwent immediate RNA extraction after removal from the culture; samples run on a nano 6000 chipset (B) Samples post SWiFT filtration, 16 weeks frozen at −80°C and a single round of flow sorting; samples run on a pico 6000 chipset.
Meeting the Need for an Effective Affordable Filtration Rig
The ability to sample the aquatic microbial community has been in place for over 30 years thanks to remotely operated bottle samplers mounted on observational instruments (James Baker, 1981). These devices monitor conductivity (salinity), temperature and underwater pressure (depth) (CTD) in real time allowing researchers to be selective with respect to the position in the water column sampled. Traditionally, water samples are manually taken from the CTD mounted Niskin bottles or pumped directly onboard the ship. Pico and micro-planktonic community cell capture then occurs via vacuum or gravity filtration onto flat membrane filters or concentrated with tangential flow filtration (Schmidt et al., 1991; Ganesh et al., 2014).
Technological advancements in recent years have moved water filtration apparatus from manned on-ship rigs to automated remote moorings able to remain at sea for months at a time (Scholin et al., 2009; Harvey et al., 2012), allowing access to water bodies otherwise inaccessible (Winslow et al., 2014). Filtration devices are even mounted on long-range autonomous underwater vehicles (AUV) able to descend the water column down to 4000 m (http://www.mbari.org/esp-technical-information/). For cruise transect sampling, robotic filtration devices are now being mounted directly onto the on-ship CTD, negating transit times required to bring deep water to the surface (Lauro et al., 2014).
While impressive in their ability to remotely sample, these advanced deployments come with two significant drawbacks. The first involves scale and recoverability. Most on-board filtration units use either flat membrane or sterivex™ filters (Millipore) with a limited capacity and poor recovery of cells for the purposes of cell sorting. These constraints effectively restrict downstream genetic studies to whole community analysis. The second drawback, which applies to the wider research community, is the build costs and expertise required for operation. At the budget end of the spectrum, development and build costs of robotic samplers is estimated around $200,000 (Lauro et al., 2014) and requires a skilled technician to maintain and operate the equipment. In this study, we have developed an alternative low cost (total build $5000), high-speed cell sampler (Figure 1) that needs no specialist skills to build, operate, and maintain. We have been able to demonstrate that SWiFT cell capture not only provides an effective planktonic cell recovery system but for the first time that post filtration cell sorting does little to either degrade or change the RNA transcript profile of a marine cyanobacterium. Application of this complete ocean sampling to sequence protocol for different phytoplankton groups could become increasingly attractive in the future as new multi-laser cytometers are now able to differentiate natural populations of phytoplankton at the class level (Thompson and van den Engh, 2016).
The Effects of SWiFT Cell Capture and How They Can Be Accommodated Within an Environmental Sampling Context
The effect of filtration on the transcription profile (Table 3) is well-defined when placed in context with the movement of cells into the dark filament cores of the HVCT concentrating device. The down-regulation of photosynthesis genes and up-regulation of genes involved in respiration is congruent with the known circadian cell programmed shift from day to night and warrants important considerations when sampling in an environmental context. CTD Niskin collection bottles are opaque effectively placing any sample collected by them into the dark. The period of time any given sample spends in the bottle is then dependent on both the depth sampled and the number of sample depths completed within the column profile. CTD transit times can take anywhere between 5 and 40 min from the point of sampling to on-deck retrieval. Cell trap filtration devices are also not alone in their reduced opacity, with several of the membrane filter holders on the market completely encapsulating flat membrane filters within a stainless steel or opaque plastic holder. Given the findings of this study, it would perhaps be unwise to perform a diel cycle sampling program using SWiFT cell capture. However, because the impact to the transcriptome is both limited (with 9% of genes affected for our model strain Synechococcus sp. BL107) and well-defined (light-dark transition), this can be accounted for when assessing the transcriptome of environmental samples.
The ability to account for sampling impacts on cellular processes is paramount if we are to derive meaningful data within an environmental context. Taking into account the effects of sampling lends confidence to our ability to detect both temporal and spatial changes in a targeted community profile. Moving toward the next stage of testing with open ocean trials, this genus targeted approach aims to identify clade-specific gene sets, with the opportunity to uncover those actively being expressed across different geographic provinces. This method will help to explore the role that specific gene sets play in defining the ecological distinctness of these lineages, giving new insights into how Synechococcus populates the environment.
FP and DS designed and conceived the project and critically evaluated practical work at all stages. FP conducted all practical work and data analysis. AM contributed to analysis and Interpretation of data. MO contributed to conceptual design and implementation of filtration apparatus. SD contributed to conceptual design of the FACS custom collection chamber. SM contributed to the design of pipeline trials. MZ contributed to field sampling application. FP and DS wrote the manuscript. All authors critically evaluated the manuscript.
The research performed in this manuscript was funded by National Environment Research Council grant NE/I00985X/1.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Authors would like to thank Sydney Cytometry for their co-operation with the development of the custom cell sort collection chamber. SD is an ISAC SRL EL.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.01592
HVCT, High Volume Celltrap; SWiFT, High Speed Water Filtration; NGS, next generation sequencing; CTD, conductivity (salinity), temperature and underwater pressure (depth); ASW, artificial seawater.
Aylward, F. O., Eppley, J. M., Smith, J. M., Chavez, F. P., Scholin, C. A., and DeLong, E. F. (2015). Microbial community transcriptional networks are conserved in three domains at ocean basin scales. Proc. Natl. Acad. Sci. U.S.A. 112, 5443–5448. doi: 10.1073/pnas.1502883112
Batmalle, C. S., Chiang, H.-I., Zhang, K., Lomas, M. W., and Martiny, A. C. (2014). Development and bias assessment of a method for targeted metagenomic sequencing of marine cyanobacteria. Appl. Environ. Microbiol. 80, 1116–1125. doi: 10.1128/AEM.02834-13
Billerbeck, S., Wemheuer, B., Voget, S., Poehlein, A., Giebel, H.-A., Brinkhoff, T., et al. (2016). Biogeography and environmental genomics of the Roseobacter-affiliated pelagic CHAB-I-5 lineage. Nat. Microbiol. 552:16063. doi: 10.1038/nmicrobiol.2016.63
Campbell, B. J., Yu, L., Heidelberg, J. F., and Kirchman, D. L. (2011). Activity of abundant and rare bacteria in a coastal ocean. Proc. Natl. Acad. Sci. U.S.A. 108, 12776–12781. doi: 10.1073/pnas.1101405108
Coleman, M. L., Sullivan, M. B., Martiny, A. C., Steglich, C., Barry, K., Delong, E. F., et al. (2006). Genomic islands and the ecology and evolution of Prochlorococcus. Science 311, 1768–1770. doi: 10.1126/science.1122050
Cuvelier, M. L., Allen, A. E., Monier, A., McCrow, J. P., Messié, M., Tringe, S. G., et al. (2010). Targeted metagenomics and ecology of globally important uncultured eukaryotic phytoplankton. Proc. Natl. Acad. Sci. U.S.A. 107, 14679–14684. doi: 10.1073/pnas.1001665107
Diamond, S., Jun, D., Rubin, B. E., and Golden, S. S. (2015). The circadian oscillator in Synechococcus elongatus controls metabolite partitioning during diurnal growth. Proc. Natl. Acad. Sci. U.S.A. 112, E1916–E1925. doi: 10.1073/pnas.1504576112
Dufresne, A., Ostrowski, M., Scanlan, D. J., Garczarek, L., Mazard, S., Palenik, B., et al. (2008). Unraveling the genomic mosaic of a ubiquitous genus of marine cyanobacteria. Genome Biol. 9:R90. doi: 10.1186/gb-2008-9-5-r90
Farrant, G. K., Doré, H., Cornejo-Castillo, F. M., Partensky, F., Ratin, M., Ostrowski, M., et al. (2016). Delineating ecologically significant taxonomic units from global patterns of marine picocyanobacteria. Proc. Natl. Acad. Sci. U.S.A. 113, E3365–E3374. doi: 10.1073/pnas.1524865113
Flombaum, P., Gallegos, J. L., Gordillo, R. A., Rincón, J., Zabala, L. L., Jiao, N., et al. (2013). Present and future global distributions of the marine Cyanobacteria Prochlorococcus and Synechococcus. Proc. Natl. Acad. Sci. U.S.A. 110, 9824–9829. doi: 10.1073/pnas.1307701110
Fuller, N. J., West, N. J., Marie, D., Yallop, M., Rivlin, T., Post, A. F., et al. (2005). Dynamics of community structure and phosphate status of picocyanobacterial populations in the Gulf of Aqaba, Red Sea. Limnol. Oceanogr. 50, 363–375. doi: 10.4319/lo.2005.50.1.0363
Ganesh, S., Parris, D. J., DeLong, E. F., and Stewart, F. J. (2014). Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone. ISME J. 8, 187–211. doi: 10.1038/ismej.2013.144
Garczarek, L., Dufresne, A., Blot, N., Cockshutt, A. M., Peyrat, A., Campbell, D. A., et al. (2008). Function and evolution of the psbA gene family in marine Synechococcus: Synechococcus sp. WH7803 as a case study. ISME J. 2, 937–953. doi: 10.1038/ismej.2008.46
Giannoukos, G., Ciulla, D. M., Huang, K., Haas, B. J., Izard, J., Levin, J. Z., et al. (2012). Efficient and robust RNA-seq process for cultured bacteria and complex community transcriptomes. Genome Biol. 13:R23. doi: 10.1186/gb-2012-13-3-r23
Gilbert, B., McDonald, I. R., Finch, R., Stafford, G. P., Nielsen, A. K., and Murrell, J. C. (2000). Molecular analysis of the pmo (particulate methane monooxygenase) operons from two type II methanotrophs. Appl. Environ. Microbiol. 66, 966–975. doi: 10.1128/AEM.66.3.966-975.2000
Gilbert, J. A., Steele, J. A., Caporaso, J. G., Steinbrück, L., Reeder, J., Temperton, B., et al. (2012). Defining seasonal marine microbial community dynamics. ISME J. 6, 298–308. doi: 10.1038/ismej.2011.107
Giovannoni, S. J., Tripp, H. J., Givan, S., Podar, M., Vergin, K. L., Baptista, D., et al. (2005). Genome streamlining in a cosmopolitan oceanic bacterium. Science 309, 1242–1245. doi: 10.1126/science.1114057
Gómez-Pereira, P. R., Fuchs, B. M., Alonso, C., Oliver, M. J., van Beusekom, J. E., and Amann, R. (2010). Distinct flavobacterial communities in contrasting water masses of the North Atlantic Ocean. ISME J. 4, 472–487. doi: 10.1038/ismej.2009.142
Gómez-Pereira, P. R., Kennaway, G., Fuchs, B. M., Tarran, G. A., Zubkov, M., and V (2013). Flow cytometric identification of Mamiellales clade II in the Southern Atlantic Ocean. FEMS Microbiol. Ecol. 83, 664–671. doi: 10.1111/1574-6941.12023
Grob, C., Jardillier, L., Hartmann, M., Ostrowski, M., Zubkov, M. V., and Scanlan, D. J. (2015). Cell-specific CO2 fixation rates of two distinct groups of plastidic protists in the Atlantic Ocean remain unchanged after nutrient addition. Environ. Microbiol. Rep. 7, 211–218. doi: 10.1111/1758-2229.12228
Hanaoka, M., Takai, N., Hosokawa, N., Fujiwara, M., Akimoto, Y., Kobori, N., et al. (2012). RpaB, another response regulator operating circadian clock-dependent transcriptional regulation in Synechococcus elongatus PCC 7942. J. Biol. Chem. 287, 26321–26327. doi: 10.1074/jbc.M111.338251
Harvey, J. B. J., Ryan, J. P., Marin, R., Preston, C. M., Alvarado, N., Scholin, C. A., et al. (2012). Robotic sampling, in situ monitoring and molecular detection of marine zooplankton. J. Exp. Mar. Bio. Ecol. 413, 60–70. doi: 10.1016/j.jembe.2011.11.022
Heywood, J. L., Sieracki, M. E., Bellows, W., Poulton, N. J., and Stepanauskas, R. (2011). Capturing diversity of marine heterotrophic protists: one cell at a time. ISME J. 5, 674–684. doi: 10.1038/ismej.2010.155
Hurwitz, B. L., and Sullivan, M. B. (2013). The Pacific Ocean virome (POV): a marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology. PLoS ONE 8:e57355. doi: 10.1371/journal.pone.0057355
Jacquet, S., Lennon, J.-F., and Vaulot, D. (1998). Application of a compact automatic sea water sampler to high frequency picoplankton studies. Aquat. Microb. Ecol. 14, 309–314. doi: 10.3354/ame014309
Jardillier, L., Zubkov, M. V., Pearman, J., and Scanlan, D. J. (2010). Significant CO2 fixation by small prymnesiophytes in the subtropical and tropical northeast Atlantic Ocean. ISME J. 4, 1180–1192. doi: 10.1038/ismej.2010.36
Johnson, Z. I., Zinser, E. R., Coe, A., McNulty, N. P., Woodward, E. M. S., and Chisholm, S. W. (2006). Niche partitioning among Prochlorococcus ecotypes along ocean-scale environmental gradients. Science 311, 1737–1740. doi: 10.1126/science.1118052
Kashtan, N., Roggensack, S. E., Rodrigue, S., Thompson, J. W., Biller, S. J., Coe, A., et al. (2014). Single-cell genomics reveals hundreds of coexisting subpopulations in wild Prochlorococcus. Science 344, 416–420. doi: 10.1126/science.1248575
Lauro, F. M., Senstius, S. J., Cullen, J., Neches, R., Jensen, R. M., Brown, M. V., et al. (2014). The common oceanographer: crowdsourcing the collection of oceanographic data. PLoS Biol. 12:e1001947. doi: 10.1371/journal.pbio.1001947
Lindell, D., Penno, S., Al-Qutob, M., David, E., Rivlin, T., Lazar, B., et al. (2005). Expression of the nitrogen stress response gene ntcA reveals nitrogen sufficient Synechococcus populations in the oligotrophic northern Red Sea. Limnol. Oceanogr. 50, 1932–1944. doi: 10.4319/lo.2005.50.6.1932
Marie, D., Partensky, F., Jacquet, S., and Vaulot, D. (1997). Enumeration and cell cycle analysis of natural populations of marine picoplankton by flow cytometry using the nucleic acid stain SYBR Green, I. Appl. Environ. Microbiol. 63, 186–193.
Markson, J. S., Piechura, J. R., Puszynska, A. M., and O'shea, E. K. (2013). Circadian control of global gene expression by the cyanobacterial master regulator RpaA. Cell 155, 1396–1408. doi: 10.1016/j.cell.2013.11.005
Mende, D. R., Aylward, F. O., Eppley, J. M., Nielsen, T. N., and DeLong, E. F. (2016). Improved environmental genomes via integration of metagenomic and single-cell assemblies. Front. Microbiol. 7:143. doi: 10.3389/fmicb.2016.00143
Olson, R. J., Zettler E., and DuRand, M. (1993). “Phytoplankton analysis using flow cytometry,” in Handbook of Methods in Aquatic Microbial Ecology, eds P. F. Kemp, J. J. Cole, B. F. Sherr, and E. B. Sherr (Boca Raton, FL: CRC Press), 175–186.
Ottesen, E. A., Marin, R. III, Preston, C. M., Young, C. R., Ryan, J. P., Scholin, C. A., et al. (2011). Metatranscriptomic analysis of autonomously collected and preserved marine bacterioplankton. ISME J. 5, 1881–1895. doi: 10.1038/ismej.2011.70
Padilla, C. C., Ganesh, S., Gantt, S., Huhman, A., Parris, D. J., Sarode, N., et al. (2015). Standard filtration practices may significantly distort planktonic microbial diversity estimates. Front. Microbiol. 6:547. doi: 10.3389/fmicb.2015.00547
Partensky, F., Blanchot, J., and Vaulot, D. (1999). Differential distribution and ecology of Prochlorococcus and Synechococcus in oceanic waters : a review. Bull. l'Instit. Océanogr. Monaco. 19, 457–475.
Penn, K., Jenkins, C., Nett, M., Udwary, D. W., Gontang, E. A., McGlinchey, R. P., et al. (2009). Genomic islands link secondary metabolism to functional adaptation in marine Actinobacteria. ISME J. 3, 1193–1203. doi: 10.1038/ismej.2009.58
Rocap, G., Larimer, F. W., Lamerdin, J., Malfatti, S., Chain, P., Ahlgren, N. A., et al. (2003). Genome divergence in two Prochlorococcus ecotypes reflects oceanic niche differentiation. Nature 424, 1042–1047. doi: 10.1038/nature01947
Saito, M. A Rocap, G., and Moffett, J. W. (2005). Production of cobalt binding ligands in a Synechococcus feature at the Costa Rica upwelling dome. Limnol. Oceanogr. 50, 279–290. doi: 10.4319/lo.2005.50.1.0279
Satinsky, B. M., Crump, B. C., Smith, C. B., Sharma, S., Zielinski, B. L., Doherty, M., et al. (2014). Microspatial gene expression patterns in the Amazon River Plume. Proc. Natl. Acad. Sci. U.S.A. 111, 11085–11090. doi: 10.1073/pnas.1402782111
Scholin, C., Everlove, C., Harris, A., Alvarado, N., Birch, J., Greenfield, D., et al. (2009). Remote detection of marine microbes, small invertebrates, harmful algae, and biotoxins using the Environmental Sample Processor (ESP). Oceanography 22, 1580167. doi: 10.5670/oceanog.2009.46
Shi, Y., Tyson, G. W., Eppley, J. M., and Delong, E. F. (2010). Integrated metatranscriptomic and metagenomic analyses of stratified microbial assemblages in the open ocean. ISME J. 5, 999–1013. doi: 10.1038/ismej.2010.189
Shilova, I. N., Robidart, J. C., DeLong, E. F., and Zehr, J. P. (2016). Genetic diversity affects the daily transcriptional oscillations of marine microbial populations. PLoS ONE 11:e0146706. doi: 10.1371/journal.pone.0146706
Sohm, J. A., Ahlgren, N. A., Thomson, Z. J., Williams, C., Moffett, J. W., Saito, M. A., et al. (2016). Co-occurring Synechococcus ecotypes occupy four major oceanic regimes defined by temperature, macronutrients and iron. ISME J. 10, 333–345. doi: 10.1038/ismej.2015.115
Steglich, C., Lindell, D., Futschik, M., Rector, T., Steen, R., and Chisholm, S. W. (2010). Short RNA half-lives in the slow-growing marine cyanobacterium Prochlorococcus. Genome Biol. 11:R54. doi: 10.1186/gb-2010-11-5-r54
Stewart, F. J., Ottesen, E. A., and DeLong, E. F. (2010). Development and quantitative analyses of a universal rRNA-subtraction protocol for microbial metatranscriptomics. ISME J. 4, 896–907. doi: 10.1038/ismej.2010.18
Sunagawa, S., Coelho, L. P., Chaffron, S., Kultima, J. R., Labadie, K., Salazar, G., et al. (2015). Ocean plankton. Structure and function of the global ocean microbiome. Science 348:1261359. doi: 10.1126/science.1261359
Teeling, H., Fuchs, B. M., Becher, D., Klockow, C., Gardebrecht, A., Bennke, C. M., et al. (2012). Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science 336, 608–611. doi: 10.1126/science.1218344
Thompson, A., Bench, S., Carter, B., and Zehr, J. (2013). Coupling FACS and genomic methods for the characterization of uncultivated symbionts. Methods Enzymol. 531, 45–60. doi: 10.1016/B978-0-12-407863-5.00003-4
Thompson, A. W., and van den Engh, G. (2016). A multi-laser flow cytometry method to measure single cell and population-level relative fluorescence action spectra for the targeted study and isolation of phytoplankton in complex assemblages. Limnol. Oceanogr. Methods 14, 39–49. doi: 10.1002/lom3.10068
Tripp, H. J., Bench, S. R., Turk, K. A., Foster, R. A., Desany, B. A., Niazi, F., et al. (2010). Metabolic streamlining in an open-ocean nitrogen-fixing cyanobacterium. Nature 464, 90–94. doi: 10.1038/nature08786
Waldbauer, J. R., Rodrigue, S., Coleman, M. L., and Chisholm, S. W. (2012). Transcriptome and proteome dynamics of a light-dark synchronized bacterial cell cycle. PLoS ONE 7:e43432. doi: 10.1371/journal.pone.0043432
Welch, D. B., and Huse, S. M. (2011). “Microbial diversity in the deep sea and the underexplored ‘rare biosphere’,” in Handbook of Molecular Microbial Ecology II: Metagenomics in Different Habitats, ed F. J. de Bruijn (Hoboken, NJ: Wiley), 243–252. doi: 10.1002/9781118010549.ch24
Winslow, L., Dugan, H., Buelow, H., Cronin, K., Priscu, J., Vesbach, C., et al. (2014). Autonomous year-round sampling and sensing to explore the physical and biological habitability of permanently ice-covered antarctic lakes. Mar. Technol. Soc. J. 48, 8–17. doi: 10.4031/MTSJ.48.5.6
Worden, A. Z., Janouskovec, J., McRose, D., Engman, A., Welsh, R. M., Malfatti, S., et al. (2012). Global distribution of a wild alga revealed by targeted metagenomics. Curr. Biol. 22, R675–R677. doi: 10.1016/j.cub.2012.07.054
Zehr, J. P., Bench, S. R., Carter, B. J., Hewson, I., Niazi, F., Shi, T., et al. (2008). Globally distributed uncultivated oceanic N2-fixing cyanobacteria lack oxygenic photosystem II. Science 322, 1110–1112. doi: 10.1126/science.1165340
Zinser, E. R., Lindell, D., Johnson, Z. I., Futschik, M. E., Steglich, C., Coleman, M. L., et al. (2009). Choreography of the transcriptome, photophysiology, and cell cycle of a minimal photoautotroph, Prochlorococcus. PLoS ONE 4:e5135. doi: 10.1371/journal.pone.0005135
Zubkov, M. V., Mary, I., Woodward, E. M. S., Warwick, P. E., Fuchs, B. M., Scanlan, D. J., et al. (2007). Microbial control of phosphate in the nutrient-depleted North Atlantic subtropical gyre. Environ. Microbiol. 9, 2079–2089. doi: 10.1111/j.1462-2920.2007.01324.x
Zubkov, M. V., Sleigh, M. A., Tarran, G. A., Burkill, P. H., and Leakey, R. J. (1998). Picoplanktonic community structure on an Atlantic transect from 50°N to 50°S. Deep Sea Res. Part I Oceanogr. Res. Pap. 45, 1339–1355. doi: 10.1016/S0967-0637(98)00015-6
Zwirglmaier, K., Jardillier, L., Ostrowski, M., Mazard, S., Garczarek, L., Vaulot, D., et al. (2008). Global phylogeography of marine Synechococcus and Prochlorococcus reveals a distinct partitioning of lineages among oceanic biomes. Environ. Microbiol. 10, 147–161. doi: 10.1111/j.1462-2920.2007.01440.x
Keywords: marine microbiology, Synechococcus, flow cytometry, transcriptomics, fluorescence activated cell sorting, RNA
Citation: Pitt FD, Millard A, Ostrowski M, Dervish S, Mazard S, Paulsen IT, Zubkov MV and Scanlan DJ (2016) A Sample-to-Sequence Protocol for Genus Targeted Transcriptomic Profiling: Application to Marine Synechococcus. Front. Microbiol. 7:1592. doi: 10.3389/fmicb.2016.01592
Received: 08 July 2016; Accepted: 22 September 2016;
Published: 14 October 2016.
Edited by:George S. Bullerjahn, Bowling Green State University, USA
Reviewed by:Lucas Stal, Royal Netherlands Institute for Sea Research (NWO), Netherlands
Katherine Rose Marie Mackey, University of California, Irvine, USA
Copyright © 2016 Pitt, Millard, Ostrowski, Dervish, Mazard, Paulsen, Zubkov and Scanlan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Frances D. Pitt, firstname.lastname@example.org