Impact Factor 4.019
2017 JCR, Clarivate Analytics 2018

The world's most-cited Microbiology journal

Original Research ARTICLE

Front. Microbiol., 26 October 2018 | https://doi.org/10.3389/fmicb.2018.02569

Organic Particles: Heterogeneous Hubs for Microbial Interactions in Aquatic Ecosystems

  • 1Leibniz-Institute of Freshwater Ecology and Inland Fisheries, Stechlin, Germany
  • 2Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany

The dynamics and activities of microbes colonizing organic particles (hereafter particles) greatly determine the efficiency of the aquatic carbon pump. Current understanding is that particle composition, structure and surface properties, determined mostly by the forming organisms and organic matter, dictate initial microbial colonization and the subsequent rapid succession events taking place as organic matter lability and nutrient content change with microbial degradation. We applied a transcriptomic approach to assess the role of stochastic events on initial microbial colonization of particles. Furthermore, we asked whether gene expression corroborates rapid changes in carbon-quality. Commonly used size fractionated filtration averages thousands of particles of different sizes, sources, and ages. To overcome this drawback, we used replicate samples consisting each of 3–4 particles of identical source and age and further evaluated the consequences of averaging 10–1000s of particles. Using flow-through rolling tanks we conducted long-term experiments at near in situ conditions minimizing the biasing effects of closed incubation approaches often referred to as “the bottle-effect.” In our open flow-through rolling tank system, however, active microbial communities were highly heterogeneous despite an identical particle source, suggesting random initial colonization. Contrasting previous reports using closed incubation systems, expression of carbon utilization genes didn’t change after 1 week of incubation. Consequently, we suggest that in nature, changes in particle-associated community related to carbon availability are much slower (days to weeks) due to constant supply of labile, easily degradable organic matter. Initial, random particle colonization seems to be subsequently altered by multiple organismic interactions shaping microbial community interactions and functional dynamics. Comparative analysis of thousands particles pooled togethers as well as pooled samples suggests that mechanistic studies of microbial dynamics should be done on single particles. The observed microbial heterogeneity and inter-organismic interactions may have important implications for evolution and biogeochemistry in aquatic systems.

Introduction

Organic particles (hereafter particles) and their surrounding solutes plume are hotspots of microbial activity in aquatic systems (Kiørboe, 2001; Simon et al., 2002). They are usually derived from several sources such as live and dead phyto- and zooplankton, fecal-pellets or allochthonous organic matter (Jackson, 1980; Alldredge and Silver, 1988; Grossart and Simon, 1993; Simon et al., 2002). Carbon and nutrients on particles often exceed background levels by 2–4 orders of magnitude (Prézelin and Alldredge, 1983; Grossart and Simon, 1998b) and thus particles are heavily colonized by microorganisms (McCarthy and Goldman, 1979; Shanks and Trent, 1979; Smriga et al., 2016) including heterotrophic bacteria (Alldredge and Youngbluth, 1985; Alldredge et al., 1986) and fungi (Wurzbacher et al., 2014). Consequently, microbial dynamics play an essential role in organic matter degradation (Alldredge and Silver, 1988; Simon et al., 2002a), aggregation and sinking and thus for the biological (De La Rocha and Passow, 2007) and microbial (Jiao and Azam, 2011; Legendre et al., 2015) carbon pump efficiency.

Microorganisms consume and solubilize particles into dissolved organic matter by a suite of ectohydrolases (Smith et al., 1992; Grossart and Simon, 1998a), but not all microorganisms possess the enzymatic machinery for all organic matter types (Kiørboe et al., 2002). Thus, microbial organic matter consumption over time will greatly affect the particle’s carbon quantity and quality, increasing the semi-labile and refractory fractions over the labile, rapidly bioavailable carbon. This will lead to pronounced temporal shifts in particle-associated community composition since initial communities will be outcompeted by more specialized ones (Kiørboe et al., 2003; Datta et al., 2016). Therefore, changes in particle composition will shape dynamics and composition of the associated communities (Grossart, 1999; Eiler and Bertilsson, 2004; Grossart et al., 2005; LeCleir et al., 2014). In addition, particle associated communities are rarely homogenous because multiple antagonistic reactions between different populations occur (Bidle and Azam, 2001; Long and Azam, 2001; Grossart et al., 2004). Although particles have been suggested to serve as microbial refuges from phages (Grossart and Tang, 2010; Tang et al., 2011), they have also been shown to harbor large numbers of viral particles and activities (Proctor and Fuhrman, 1992; Peduzzi and Weinbauer, 1993; Müller-Niklas and Schuster, 1994; Riemann and Grossart, 2008). Until now, relatively few studies have followed temporal changes in bacterial community composition on particles (Grossart and Simon, 1998a; Schweitzer et al., 2001; Grossart et al., 2003a; Kiørboe et al., 2003; LeCleir et al., 2014; Datta et al., 2016), and even fewer studies tried to link it to activity on particles (Schweitzer et al., 2001; Kong et al., 2013; Moran et al., 2013; Smith et al., 2013; D’Ambrosio et al., 2014). Most of these studies focus on single-point measurements and thus do not provide any temporal resolution. Those that do, however, have used closed rolling tank systems (Grossart and Simon, 1998a; Schweitzer et al., 2001; Datta et al., 2016), and are likely to be affected by a series of phenomena summed under the term “the bottle effect” (Zobell and Anderson, 1936). These include a rapid shift in microbial community composition and an increase in heterotrophy vs. autotrophy (Ferguson et al., 1984; Calvo-Díaz et al., 2011; Baltar et al., 2012).

To overcome this “bottle effect” and allow for constant exchange of nutrients and microorganisms from the environment, we connected flow-through rolling tanks (Ionescu et al., 2015; Figure 1) directly to the meso-oligotrophic Lake Stechlin (northeastern Germany). The tanks were inoculated with pre-formed particles (see section “Materials and Methods”) of an axenic diatom (Navicula sp.) previously isolated from the lake. Particles were collected at several defined time points for respiration and photosynthesis measurements (Ionescu et al., 2015), bacterial community composition, metagenomic, and metatranscriptomic analyses of replicates of 3–4 pooled particles (hereafter single particles). These were further compared to hundreds to thousands of pooled particles (hereafter pooled particles; Figure 2) by filtering the ambient water in the flow-through reactor. Last the single-particles and pooled-particle samples were compared to the results obtained when pooling RNA from two biological replicates. First results (Ionescu et al., 2015) show that under near in situ conditions, and as long as light is available, photosynthesis remains active for up to 9 days, twice as much as in a closed rolling tank. Accordingly, given continuous production of photosynthesis products during these 9 days as well as an expected continuous supply of fresh labile organic matter from the lake, respiration on particles remains high under near in situ conditions for the whole duration of the experiment. This is in clear contrast to a closed incubation system, which excludes the supply of fresh, labile organic matter and thus reveals a rapid decrease in respiration with time (Ionescu et al., 2015).

FIGURE 1
www.frontiersin.org

FIGURE 1. Schematics of the flow-through rolling tank system (A) with live diatom particles during an experiment (B). The conceptual experimental setup is shown in panel (C). The incubation chamber of the flow-through rolling tank is 9 cm in diameter, 30 cm long and has a volume of ca. 1.9 L. Extensive technical details are given in Ionescu et al. (2015).

FIGURE 2
www.frontiersin.org

FIGURE 2. The sampling procedure (A) included collection of 3–4 particles (referred to, when pooled, as single particles; sPs) from parallel incubating rolling tanks (marked here as R1 and R2 for replicates 1 and 2) after which the remaining water, pooling lake particles (referred to as pooled particles; pPs) and microscopic fragments of inoculum particles were filtered on a 5 μm pore size filter. RNA was extracted and mRNA was enriched from each fraction collected. Pooled samples consist of equi-volume mixing of mRNA of two single particle or two pooled lake particles samples. The sampling schemes (B) depict the start and end of incubation times in April 2013. Samples were collected at the end of the incubation period. A subset of the functional annotation results is shown, demonstrating the different results obtain from single sample vs. pooled sample analysis in both single particles and pooled lake particles (C). In subsequent figures the first and second starting replicates of each incubation type (long or short) are named R1–R4 for replicates 1–4.

We hypothesize that: (i) Community composition as well as transcriptomic profile of pooled particles differ from those of individual particles; (ii) In the photic zone, changes in carbon-quality, as can be interpreted by transcriptomic profile of photosynthetic particles are slow due to active photosynthesis and constant supply of fresh, labile organic matter thus minimizing carbon-quality-based changes in the particle-associated microbial community composition; (iii) Initial microbial colonization is random and; (iv) early changes in microbial communities on particles are mainly due to dynamics in microbial interactions.

Our data calls for modifications of current sampling and processing strategies of particles. Individual particle studies are required to address the degree of heterogeneity in microbial dynamics on particles of similar (e.g., homogenous blooms) or diverse organic matter sources.

Materials and Methods

Algal Culture

An axenic culture of the freshwater diatom Navicula sp., isolated from Lake Stechlin, was grown in 900 mL of sterile Z-medium (Staub, 1961) in 10 sterile, pre-combusted 1 L Duran bottles (Schott). The culture was grown under a 12 h dark-light cycle at 15°C and 30 μmol photons m-2 s-1 for 2 weeks. Afterward, bottles with axenic diatoms were placed on a roller table at half the light intensity for particle formation prior to in situ inoculation with freshwater bacteria from Lake Stechlin.

Experimental Setup

To conduct long-term experiments with a constant supply of fresh OM, nutrients, and bacteria from the habitat while avoiding the “bottle effect” we used a newly designed flow-through rolling tank system (Ionescu et al., 2015; Figure 1A). Shortly, this device consists of an inner tube in which the continuous sinking of particles is simulated and an external casing through which lake water flows in and out of the inner tube without disturbing the particles (Figure 1B). A detailed description necessary for construction of the flow-through rolling tank is available in Ionescu et al. (2015).

Freshly formed axenic diatom particles were cautiously transferred to the prefilled inner core of the rolling tank device carefully avoiding air bubbles. The axenic nature of the diatom culture renders fresh particles rather fragile (Gärdes et al., 2011). Thus, since some particles may break apart and exit the system, the flow-through rolling tank was inoculated with ca. 100 particles 1–2 mm in diameter each. The device was then sealed and allowed to roll for up to an hour prior to onset of lake water flow. Lake water was pre-filtered through a 100 μm pore size sieve to prevent large objects including mesozooplankton from entering the experimental setup. Flow was maintained at ca. 300 ml h-1 and rotation speed was set at 2 RPM keeping the particles suspended and without contact with the walls of the flow-through rolling tank, mimicking a sinking speed of ca. 7 m d-1.

Sampling Procedure (Figure 2A)

Each rolling tank was opened, and 3–4 individual Navicula particles (referred to as single particles when pooled) were carefully sampled using a cut-end syringe and avoiding collection of surrounding water. Collected particles were rapidly gathered and placed on a 5 μm pore size polycarbonate filter (47 mm in diameter, Sartorius, Germany) to remove excess water and free-living bacteria. Subsequently, the filters were immediately placed in Z6-buffer (8 M guanidinium-HCl, 20 mM MES, 20 mM EDTA [pH 7.0], and 0.7% [v/v] 2-mercaptoethanol) to deactivate RNAses and denature all proteins. To collect the lake particles that include particle from the lake and microscopic particles from the Navicula innoculum, the total water volume of each rolling tank (1.8 L) was filtered through 5 μm pore size filters. Filters were changed each time the filtration speed slowed down to avoid clogging and collection of free bacteria on the pooled lake particles (pooled lake particles). All filters were immediately placed in 700 μL of Z6 buffer and placed on ice for further analyses.

A total of eight individual rolling tanks were set up for the experiment forming four biological replicate experiments for short and long incubation each as depicted in Figure 2B. Initially, four reactors were run, two of which were sampled after 24 h and two after 8 days. Two additional reactors were inoculated and run for 7 days once the first set had been sampled. A second two-reactor set of one-day incubations was inoculated on the 8th day of the experiment thus testing for specific colonization of particles, given an 8-day time gap.

DNA and RNA Extraction

Nucleic acids were extracted as described for DNA in Ionescu et al. (2012), except for the last desalting step. Shortly, cells were lysed by incubation for 15 min. at 95°C in a Z6-buffer followed by 15 min. incubation at 65°C with phenol (0.5:1 phenol:sample, v:v; Phenol pH 8 and 4.3 for DNA and RNA, respectively) after which the aqueous phase was extracted with chloroform. After an additional phenol-chloroform (1:0.5:0.5, sample:phenol:chloroform, v:v:v) extraction, the samples were extracted twice with chloroform and precipitated over night at -20°C using 2-propanol (1:1). The DNA was cleaned in 70% ethanol and dissolved in Diethylpyrocarbonate (DEPC) treated water. Following the extraction, RNA, samples were treated with Turbo DNA Free (Ambion) according to the manufacturer’s instructions.

mRNA Enrichment and cDNA Synthesis

Since ribosomal RNA is the most common transcript in RNA extractions we have enriched the mRNA by removing all rRNA transcripts following the procedure described in Stewart et al. (2010). Shortly, we generated long 16S and 23S rRNA probes by obtaining sample specific rRNA gene PCR amplicons which were later transcribed to biotin labeled sample-specific RNA probes using the MegaScript T7 kit (Ambion) and biotin labeled C and U (Roche). To achieve maximum rRNA removal, all probes from all samples were pooled and equal aliquots were added to the RNA samples. Following an incubation of 5 min. at 70°C, the hybridized rRNA molecules were extracted with streptavidin coated magnetic beads. The mRNA enriched RNA was concentrated by using the Easy RNA kit (Qiagen). First strand cDNA was synthesized using the Superscript III kit (Invitrogen) according to the manufacturer’s instructions.

Pooled samples consist of an equal volume of mRNA-enriched RNA from a pair of individual samples of either single particles or pooled lake particles. Equal volumes were chosen over equal RNA concentrations to better represent particles of similar sizes colonized at different microbial densities with additionally different level of activity.

Sequencing

Shotgun Sequencing From DNA

Metagenome sequencing steps included fragmentation of genomic DNA, ligation to sequencing adapters and purification. Following the amplification and denaturation steps, libraries were pooled and sequenced. We used 50 ng of DNA from each sample to prepare the libraries using Nextera DNA Sample Preparation Kit (Illumina). Library insert size was determined by Experion Automated Electrophoresis Station (Bio-Rad). The insert size of the libraries ranged from 300 to 1400-bp. The pooled library (12 pM) was loaded to a 600 Cycles v3 Reagent cartridge (Illumina) and the sequencing was performed on a Miseq (Illumina) sequencer. Sequencing was carried out at Molecular Research Laboratories (Mr. DNA), Shallowater, Texas.

cDNA Sequencing

cDNA sequencing steps include purification and generation of blunt end cDNA followed by ligation to sequencing adapters, amplification, denaturation, and sequencing. We used 250 ng of double strand cDNA from each sample to prepare the libraries using TruSeq RNA sample preparation kits (Illumina). The pooled library (10 pM) was loaded to a 600 Cycles v3 Reagent cartridge (Illumina) and the sequencing was performed on Miseq (illumina). cDNA sequencing was carried out twice, the first-time biological duplicates were pooled and the second time each sample was sequenced individually. The results of the first sequencing are provided as Supplementary Information. Sequencing was carried out at Molecular Research Laboratories (Mr. DNA), Shallowater, Texas.

Sequence Quality Control

Sequence quality was assessed using the FastQC program (Babraham Institute) after which each sequence file was cleaned from low quality sequences and traces of the sequencing method with the Clip tool from the Nesoni package (Victoria bioinformatics) using the following parameters: the first 12 and last 3 bases were trimmed and sequences shorter than 50 nucleotides were removed; from the remaining data, sequences with an overall quality score below 20 were removed as well.

Sequence Analysis

The metagenome was analyzed on the MG-RAST server (Meyer et al., 2008).

Transcriptome

All paired end and single end libraries from the transcriptome analysis generated during the sequence quality control procedure were concatenated into three study-wide data files: left, right and non-paired reads, to maximize the chances of transcript assembly.

Transcripts were assembled and analyzed using the Trinotate pipeline1. This includes transcript assembly with the Trinity pipeline (Grabherr et al., 2011b), open reading frame prediction using TransDecoder (part of the Trinity software package), protein identification using BlastX and BlastP (Altschul et al., 1990) against the Uniprot protein database (The uniProt Consortium, Bateman et al., 2017), protein domain prediction using HMMER (V3.1; Finn et al., 2011) against the Pfam database (Finn et al., 2014). To obtain additional annotation data, all transcripts were run through the KEGG KAAS system (Moriya et al., 2007). The basic analysis pipeline is amended with data analysis against functional databases (CAZY; Lombard et al., 2014), and structural databases (SCOP; Fox et al., 2015 and SUPERFAMILY; Wilson et al., 2009) to maximize annotation of unknown proteins. A prediction of polyketide synthetases and non-ribosomal peptides synthetases was done using a locally installed version of ANTISMASH (Weber et al., 2015) and NaPDoS (Ziemert et al., 2012).

Differential expression analysis was conducted using edgeR (Robinson et al., 2010) (v 3.14) using the FPKM normalization mode (Wagner et al., 2012). The 16 samples (eight in each group: single particles or pooled lake particles) were compared grouped by samples (generating duplicates per condition) or by duration of incubation (short vs. long; generating four replicates per condition). Benjamini-Hochberg false discovery rate test and Bonferoni corrected p-values were used to test for statistically significant difference in transcriptomic profile.

The presence of CRISPRs (Clustered Regularly Interspaced Short Palindromic Repeats) was tested using CRISPRfinder (Grissa et al., 2007). Spacer-sequences from candidate CRISPRs were than compared to the entire transcriptome data as well as to the assembled transcripts using Bowtie2 and BLAST, respectively.

Reconstruction of Bacterial Community

Since rRNA was removed from the samples, microbial community composition was reconstructed from the raw reads as annotated on the MG-RAST server (Meyer et al., 2008) against the M5NR database.

Species Co-occurrence Networks

Species–species association networks were constructed using CoNet (Faust and Raes, 2016) as implemented in Cytoscape 3.3.0 (Shannon et al., 2003) software. This program was chosen as it allows the merging of networks calculated via several association methods of which we chose: Kendall, Pearson and Spearman correlations, as well as Steinhaus and variance log ratio similarities. The latter results in values between 0 and 1, whereas the rest range between -1 and 1. Only edges confirmed by at least two methods each with a score either smaller than -0.6 (negative association) or larger than 0.6 (positive association) were retained in the final network. The resulting network was visualized and graphically organized using Cytoscape 3.3.0 (Shannon et al., 2003).

Statistical Analysis

Principle component analysis was done on gene and transcript abundance matrices produced by the Trinotate1 pipeline as well as on the function abundance matrix generated by the MG-RAST (Wilke et al., 2016) pipeline. The analysis was done using Primer software version 6 supplemented with the Permanova package. Permanova analysis was conducted on the same data using duration of incubation (short for 1 day or long for 7–8 days) and sample type (single particles or pooled particles). A similar division was used to conduct analysis of similarity (ANOSIM), forming four groups with four member-replicates each. Prior to all analyses the data was square-root transformed and standardized by total reads.

Results

Community Composition

The Active Bacterial Community (mRNA) on Single Particles Was Heterogeneous and Differed Between Biological Replicates and Incubation Time (Figure 3)

FIGURE 3
www.frontiersin.org

FIGURE 3. Community composition of single particles (3–4 individual particles pooled) and pooled lake particles according to MG-RAST metagenomics annotation. (A) The distribution of Archaea, Bacteria, Eukarya, and Viruses; (B) Bacterial families with a sequence frequency of more than 5% in at least one sample; (C) Eukaryotic families with a sequence frequency of more than 5% in at least one sample; (D) Virus families with a sequence frequency of more than 5% in at least one sample. The lower stacked bars in each panel represent individual samples (Replicate reactors R1–R4) while the upper ones the sequencing results of poled RNA from the pair of sample below. (E) Venn diagrams showing the taxonomic overlap between single particle samples (R1–R4) and pooled RNA extracts (R1 + R2 or R3 + R4) at the family level including Bacteria, Archaea, Eukaryotes, and Viruses. The diagrams do not account for abundance but only presence absence. The number of taxonomic entities in each group is shown in each circle of the diagram. The shading color of the text in the legend of the Venn diagram depicts the color of the individual sample circles. The color of squares represents the overlapping between the different samples.

Community data was deduced from the taxonomic annotation of coding RNA (i.e., protein coding), therefore, it is unlikely that it was affected by our rRNA removal protocol. Short-term incubations (1 day), revealed that reads from single particles (3–4 pooled individual particle; see Figure 2 and section “Materials and Methods”) consisted almost exclusively of bacterial sequences almost randomly distributed between the samples. After long-term incubations (7 and 8 days) sequence frequencies of eukaryotes (mostly fungi) and viruses increased (Figure 3A). Whereas sequence frequencies of Bacteria on pooled particles (Figure 2) were initially low, they increased in the long-term incubations. Average sequence frequencies of viruses on pooled lake particles remained similar over time (Figure 3A).

Microbial Composition of Single Particles Differed From That of Pooled Duplicate Samples Across All Domains (Figures 3A–D)

Initially, identified taxa (families) from single particles represented a small subset of those identified in the pooled samples. In the longer incubations, however, the overlap in present taxa (i.e., presence absence similarity) increased (Figure 3E), but sequence frequencies still played a major role for the observed differences between single and pooled samples.

Total Bacterial Communities (DNA) of Both Single Particles and Pooled Lake Particles Were Similar, but Greatly Differed From Active Communities (RNA) (Supplementary Figure S1)

While the total bacterial community on single particles and pooled lake particles was dominated by Caulobacteraceae, Commamonadaceae, and Pseudomonadaceae, these families were either absent or significantly less dominant in the respective active communities (Figure 3B). The dominant active families consisted of Flavobacteriaceae, Pseudomonadaceae, Cytophagaceae, Enterobacteraceae, Moraxellaceae, Sphingomonadaceae, and Bacteroidaceae (Figure 2B). Bacterial communities on pooled lake particles, particularly of long-term incubations, were more similar to each other than single particles (Figure 3B). By using rRNA gene amplicon (i.e., DNA) sequencing, Ionescu et al. (2015) showed that the total microbial community in flow-through rolling tanks was rather stable for a long time (>1 week).

The active eukaryotic community from the long-term incubations of single particles and all pooled lake particles was dominated by Plectosphaerallaceae, a family of fungi (previously identified and isolated from Lake Stechlin; Wurzbacher et al., 2016), and several dinoflagellate families. The eukaryotic community, except for single particles in short-term incubations, was more homogenous than for bacteria (Figure 2C). Pooled lake particles had higher sequence frequencies of eukaryotes (Figure 2A) indicating more diverse source particles. Similar to the active bacterial communities, eukaryotes are highly alike between the total single particles and pooled lake particles communities (DNA), but greatly differed from the active community (RNA).

Viruses made up only a small portion of the reads from the short-term incubated single particles and consisted of Myoviridae, Siphoviridae, and unclassified viruses (Figure 2D). However, with the increase in viral sequence frequency in the long-term incubations, these families were replaced almost entirely by sequences of unclassifiable viral origin and a few Circoviridae. Viruses on pooled lake particles remained mostly unclassified, except for one long term incubation revealing mostly Myoviridae, Siphoviridae, and Phycodnaviridae. Interestingly, these virus families were not abundant on the single particles from the same incubation.

Functional Annotation

The Trinotate assembly pipeline1, as applied to the pooled sequences of all samples, (Supplementary Data Sheets 1–3) resulted in 48,787 transcripts larger than 90 nt (30 amino acids), located on 25,049 contigs. Of these, 22,585 transcripts could not be assigned to any potential function by any of the analysis tools used (Supplementary Tables S1, S2).

One of our major aims was to test whether a rapid change in particle carbon quality also occurs under near in situ conditions. This should result in noticeable changes in the particle transcriptome profile and possible in shifts in microbial communities from generalists to more specialized communities. The particle aging period employed here (7–8 days) is on par with other studies reporting changes in expressed gene profiles (Datta et al., 2016; 6 days). Our results suggest grouping of multiple particles masks heterogeneity among particles (Figure 3). Thus, all transcript combinations were analyzed for differential functional gene expression (four short-term single particles vs. four long-term single particles incubations and two short-term pooled particles vs. two long-term pooled particles incubations). These analyses were conducted on annotated transcripts, translated proteins as well as on the raw reads as annotated by MG-RAST (Wilke et al., 2016). In all cases, the graphical representation suggests that the genes expressed on the single particles in the long-term incubations bare more similarity to the long-term incubation in the lake. This makes sense since with time, given the nature of the open system, there would be lake particles adhering to the experimental inoculum as well as the likely establishment of an equilibrium with the lake community. Nevertheless no significant difference in the expression of any gene between the short- and long-term incubations as well as between single and pooled particles occurred when using Bonferroni corrected p-values (p = 1), Benjamini–Hochberg false discovery rate test (FDR = 0.9–1), or Permanova analysis (Figure 4).

FIGURE 4
www.frontiersin.org

FIGURE 4. Principle component analysis of raw reads annotation using the MG-RAST pipeline (A), assembled genes (B), and assembled transcripts (C) from single and pooled particles (upper panels) or single particles alone (lower panel) grouped by the length of the incubation period where “short” refers to 1 day and “long” to 7/8 days incubations.

To test whether changes in expression occur only in the subset of genes specific to enzymatic degradation and substrate utilization these were extracted from the MG-RAST annotation (Figure 5B). A similar approach was applied to sugar, nitrogen, phosphorus, and iron metabolism (Supplementary Figures S2S5). These analyses showed a large heterogeneity between single particles with no significant difference between young (1 day) or older (7–8 days) particles (ANOSIM p => 0.2 in all cases).

FIGURE 5
www.frontiersin.org

FIGURE 5. A subset of functional modules involved in substrate degradation and utilization obtained from single particles (3–4 individual particles pooled) and pooled lake particles transcriptomes (A) The lower stacked bars represent individual samples while the upper ones the sequencing results of poled RNA from the pair of sample below. The shown modules are Level-3 MG-RAST subsystem hierarchy. Venn diagrams (B) show the overlap between all the recognized modules including the subset in panel A, between single particle samples (R1–R4) and pooled RNA extracts (R1 + R2 or R3 + R4). The number of modules identified in the transcriptome of each sample is given in each circle of the diagram. The shading color of the text in the legend of the Venn diagram depicts the color of the individual sample circles. The color of squares represents the overlapping between the different samples.

Initially, as for the colonizing communities, the functional modules identified on single particles represented a small subset of the pooled single and lake particles (Figure 5A). At the end of the experiment, however, there was a much higher overlap between single particles and pooled single and lake particles.

Considering the ongoing photosynthetic activity of diatoms forming the particles (Ionescu et al., 2015) and the continuous supply of fresh, labile organic matter (both dissolved and particulate) pumped through the vessel from the lake, there are two main explanations for the absence in differences between expression profiles of particle associated communities on fresh vs. aged particles (Figure 5B). First, communities of aged particles still consist of initial colonizers. Second and more likely, competition for resources and other inter-organismic interactions led to a community replacement by others with similar functions. Although our dataset is limited to testing the first possibility as individual particles couldn’t be sampled at multiple time points and particle numbers were insufficient for a proper statistic, we searched for indications of inter-organismic interactions.

Possible Evidence for Antagonistic Reactions

Anti-phage Activity

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) sequences consist of several similar stem-and-loop palindromic sequences separated by spacers of different sequence (Sorek et al., 2008). To identify potential CRISPR sequences, we analyzed the metagenomes of single particles and pooled lake particles and revealed 735 putative CRISPRs including 1066 different spacers (Supplementary Table S3; Supplementary Data Sheet 4). Among the metagenomic CRSIPRs 222 had also matches in the assembled transcriptome. Additionally, we found the expression of different CRISPR associated genes including Cas1, Cas2, Cas3, and Cas6.

A second defense mechanism employed by bacteria against foreign viral DNA is the use of restriction endonucleases (restriction enzymes; Weinbauer, 2004). Overall multiple Type I, II, IIS, and III endonucleases were identified in the Trinity (Grabherr et al., 2011a) assembled transcripts and the MG-RAST analysis.

Interactions Between Other Organisms

In addition to viruses, further interactions between organisms such as antagonistic interactions between bacteria or prokaryotic and eukaryotic organisms affect microbial community composition (Grossart et al., 2004). Such interactions are evident in the expression of protein secretion system type II, III, IV (Supplementary Data Sheet 1), and VII (Supplementary Data Sheet 2), genes involved in antibiotic production and defense mechanisms as well as the expression of toxin-antitoxin systems (Supplementary Data Sheets 1, 2). Several variants of polyketide synthases and non-ribosomal peptide synthetases were found as well (Supplementary Table S4).

Discussion

Our study using an open, flow-through particle incubation system that significantly minimizes or removes the “bottle effect” portrays the important role of particles as surfaces and hotspots for microbial colonization, metabolic activity, and interactions in the water column. First, we show a large heterogeneity between microbial communities attached to single particles of identical source and age. Second, we show a high heterogeneity in specific activities on these particles. Last, expression of genes known to be involved in antagonistic reactions between organisms including viral and anti-viral activity highlight the importance of organismic interactions for microbial community dynamics on particles. These findings highlight particles as hotspots of biogeochemical organic matter cycling in the water column with profound consequences for the microbial carbon pump (Jiao and Azam, 2011) and overall carbon pump efficiency (Legendre et al., 2015). Furthermore, differences between single particles, pooled lake particles and pooled samples highlight pitfalls to mechanistically understand temporal dynamics of microbial particle degradation when averaging samples. Averaging multiple particles of different ages and sources as is the case in nature can only provide in inventory of organisms and processes occurring on particles. As we show, each individual particle contributes only a subset of this information providing specific information on the microbial assemblages and their particle-degrading activity.

Particle type as a proxy for carbon quality has been suggested to play a significant role in composition and activity of the attached bacterial community (Grossart et al., 2005; Teeling et al., 2012). Datta et al. (2016) by incubating particles in a closed system demonstrated that microbial colonization and succession of total communities (DNA based) on artificial particles of the same type are reproducible. Our study, minimizing the “bottle-effect,” demonstrates that, under near in situ conditions, particle associated communities on single particles of identical source have a high inter-particle heterogeneity. This heterogeneity between biological replicates indicates that the nature of the initial particle colonizing community is not a sole function of particle type and composition, but also of other – perhaps partially stochastic events (Cordero and Datta, 2016; Datta et al., 2016), as discussed below. This suggests that in any given environment some members of the residing microbial community will be able to colonize particles. When extrapolating to the natural aquatic environment, where particles of multiple origins are simultaneously formed, it is likely that these will be rapidly and differently colonized regardless of the surrounding microbial community.

Under near in situ conditions, the initial degradation process does not seem to select for specific microbial groups. On aged particles (7–8 days), microbial communities are as heterogeneous as on “fresh” particles (1 day). The higher similarity between pooled lake particles suggests that heterogeneity between individual particles is masked by averaging over large volumes (i.e., ∼2 L in this case) encompassing hundreds of particles of variable type and age. A previous study comparing, using microscopy, microbial colonization on several thousands of individual particles collected from a single environment – similarly reported a large heterogeneity between particles (Bižić-Ionescu et al., 2015).

Conventional molecular studies on particles have mainly used size differential filtration to obtain various particle size classes (Padilla et al., 2015). Thus, particles are usually treated as an integral of various entities with unknown history and chemical quality ignoring almost completely inter-particle heterogeneity. Our study demonstrates for both microbial community composition and activity, as reflected by transcriptomic profiles, that the sum of two samples, either of single particles or of larger volumes (pooled samples of already pooled lake particles) differs from single particles, in particular during early colonization. This is probably the result of inter-particle heterogeneity in microbial community composition, activity, and colonization density. Furthermore, numerous pooled particles mask the composition and activity of specific subsets of the community associated with individual particles. Therefore, sample or particle pooling does not allow for an understanding of microbial interactions involved in the degradation of a particle of a specific age and/ or source. As particles age, they seem to be more homogenous both in community composition and functionality. Therefore, samples consisting of fresh and aged particles are even more problematic for understanding the initial steps of particle colonization and degradation. To better understand bacterial colonization and succession on particles in relation to their specific organic matter quality, either the use of model particles and organisms is required (Enke et al., 2018) or metagenomic and metatranscriptomic studies of individual single particles with well-defined organic matter composition are required. Findings from individual particles of similar age and composition are likely to reveal functional redundancy among individual microorganisms as well as between different assemblages. Additionally, information on microbial assemblages from individual particles will show if certain microorganisms are exclusively mutualistic and perhaps actively recruit one another from the water column, while others always exclude one another.

Potential Determinants of Particle-Associated Microbial Community Composition

High Colonizer Patchiness

Stochasticity of attachment to surfaces, based on the “first comes, first served” concept may play a major role in determining the composition of the initial particle associated bacterial community. Bacterial distribution in water is generally patchy (Cordero and Datta, 2016). The presence of particles may further affect this patchiness by releasing dissolved organic matter and nutrients and thus attracting different chemotactic bacteria (Kiørboe et al., 2002, 2004; Stocker, 2012). Bacterial populations differ in their affinity to the substrates released from particles, and consequently different particles may attract and be colonized subsequently by different bacteria.

Colonizer–Particle Interactions Differ Between Individual Particles

Microbial utilization of particle resources is unlikely to occur without any competitive interactions between the numerous potential groups of colonizers (Grossart et al., 2003a,b). Studies have shown in vitro antagonistic behavior between isolates obtained from particles (Long and Azam, 2001) and competitive particle colonization (Grossart et al., 2003b). Interactions can occur bi-directionally between different populations: (a) populations of eukaryotes such as particle-forming algae and colonizing fungi or among populations of colonizing fungi; (b) bacteria and eukaryotes such as particle-forming algae and colonizing bacteria (Amin et al., 2012) or between colonizing bacteria and fungi; (c) multiple colonizing bacterial populations including protozoan grazers (Long and Azam, 2001; Grossart et al., 2003b), and (d) colonizing organisms and viruses (Riemann and Grossart, 2008). However, solely analyses of single particles have the potential to identify antagonistic interactions between various colonizing organisms and thus can discriminate the role of particle age and quality from other factors. Accordingly, our data support the presence of competitive microbial interactions, e.g., the attempt of particle-associated bacteria to inhibit the particle-forming diatom.

Bacterial – host antagonism is indicated by the expression of type II, III, IV, and VII protein secretion systems (Supplementary Data Sheets 1, 2). Furthermore, evidence for toxin production in our transcriptome alongside with various resistance mechanisms (Supplementary Data Sheets 1, 2) indicate a plethora of antagonistic interactions. For example, genes involved in the synthesis of phenazine, a toxic Pseudomonas pigment (Gibson et al., 2009), were also found among the particle transcripts. This pigment, both in natural and modified forms, was shown to be toxic for several eukaryotic organisms (Cezairliyan et al., 2013).

Multiple Bacteria–Bacteria Interactions

We detected the expression of a few genes directly involved in the synthesis of known antibiotic compounds such as mitomycins and carbapenem as well as several genes predicted to be involved in antibiotic production. A search for polyketide synthetases and non-ribosomal peptide synthetases showed a few matches as well (Supplementary Table S4). On the other hand, numerous genes such as efflux pumps, transporters, and hydrolases associated with resistance to antibiotics such as polymixin, erythromycin, spectinomycin, and streptomycin occur on particles. This supports our hypothesis of an antagonistic and hence partly competitive behavior between different particle colonizing microorganisms (e.g., bacteria and fungi).

The calculation of co-occurrence networks from the single particles associated communities (Supplementary Figure S6) showed that there is no overlap in species associations between short- and long-term incubations. Similarly, no overlap is observed between the single particles and the pooled lake particles networks (data not shown). Hence, while each network in its own suggests that some species co-occur while other are mutually absent, these interactions seem to be highly particle specific.

Intense Bacteria–Virus Interactions

Bacterial and eukaryotic viruses are strong regulators of microbial communities (Weinbauer, 2004) including those associated to particles. The strong increase in viral sequences on single particles in the long-term incubations indicates that as initial populations settle on the particles they become subject to viral attacks. This is further evidenced by the expression of CRISPRs, the bacterial “immune system.” The constant presence of viruses on the readily colonized pooled lake particles further strengthens the significant, yet underestimated role viruses play on particles. Our microscopic analysis of the incubated particles has detected cells with viral like entities on them (Supplementary Figure S7).

No Changes in Microbial Carbon Usage During the Incubation Period

A sequential change in carbon quality is commonly believed to control bacterial successions on particles during a few days of incubation. Nevertheless, indications for such changes come from studies conducted in fully or partially closed systems. In such systems, the bacterial community is strongly altered by the “bottle effect” leaning toward heterotrophy (Sherr et al., 1999; Calvo-Díaz et al., 2011), and naturally occurring DOM is not continuously available. Contrary to these studies, Ionescu et al. (2015) have shown that employing an open (flow-through) experimental system as used in this study, the decrease in algal activity due to bacterial degradation occurs at least 4 days later than in a closed system. Furthermore, respiration on particles incubated in the open system was maintained at a high level for a long period suggesting a labile carbon source was readily available. Although heterogeneity in particle-associated communities was high, we could not detect any (statistically significant) change in the expression of genes, particularly those related to carbon usage. For natural systems, we thus propose that labile carbon is readily available during particle sinking through the upper water layers. Particle sinking velocities range from as low as 1 cm d-1 (Estrum-Yousef et al., 2000) to up to at least 250 m d-1 (McDonnell and Buesseler, 2010). Hence, some particles will remain in the surface layers of aquatic bodies for rather long periods. However, for a substantial fraction of particles 7–8 days (the duration of our experiment) is too long to stay in the water column of a shallow water body, and surely not in the photic zone of even the open ocean. Studies employing image analysis (Ratmeyer and Wefer, 1996; Nowald et al., 2009), echo sounding systems (Jackson and Checkley, 2011) or sediment traps (McDonnell and Buesseler, 2010) to follow particles throughout the water column, did not find a conclusive decrease in particle size with depth. Thus, while organic matter from photosynthetic activity is no longer available, labile organic matter can reach particles from the dissolved organic matter pool of the surrounding water or from adherence of additional particulate OM.

Hotspots of Antagonism as a Source of Genetic Variability

The high number of genes involved in inter-organisms interactions and the strong evidence for viral presence and activity support particles being a hotspot for lateral gene transfer (Riemann and Grossart, 2008; Walsh et al., 2013). The phage-mediated transfer of genetic material between different bacteria in aquatic environments is thought to be limited by host abundance in the water column (Weinbauer, 2004). This would probably not be the case on the limited surface of densely colonized particles, where close vicinity and high number of possible vectors may result in numerous transduction events. Our transcriptomes include both Type II and Type IV pili secretion systems needed for the uptake of naked DNA, as well as the recA and recR recombination proteins. Additionally, we found the expression of the comEA genes which encode for a protein essential for DNA uptake (Provvedi and Dubnau, 1999). Thus, in addition to gene exchange mediated by virus between (usually) closely related organisms, uptake of naked DNA can also increase genetic exchange between even phylogenetically and functionally distant microorganisms.

Conclusion

Our study demonstrates that single particles represent complex micro-niches in the water column which are randomly colonized by a diverse microbial assemblage from the surrounding water (Figure 6A). We propose that competitive and antagonistic activities on particles, between all prokaryotic and eukaryotic entities play an essential, yet underestimated role in determining dynamics of particle-associated community composition (Figure 6B). Depending on particle sinking rate, the available labile organic matter can be continuously replenished by internal and external processes. Thus, carbon driven changes in microbial community composition and activity, as derived from experiments using closed incubation devices, are absent or less pronounced in the upper (euphotic) water layers (Figure 6C). It is obvious that different dynamics in organic matter composition of particles incubated in closed or open systems result in different conclusions on microbial community dynamics both in relation to community composition and activity. Thus it is important to better understand and model microbial carbon pump efficiency as current global climate change will result in an increased retention time of particles in the surface layers of water bodies (Legendre et al., 2015). Continuous inter-organism activities on particles render them generators of genetic diversity in aquatic ecosystems (Figure 6D) with increased horizontal gene transfer between particle-associated organisms as compared to their free-living counterparts (Arias-Andres et al., 2018). Studies averaging over multiple particles are valuable in providing a broad catalog of organisms and functions related to particles but are restricted in providing information on specific particle associated consortia and their functions. Therefore, while “All roads may lead to Rome” individual members and physiological processes may be deduced only from focused studies using single particles.

FIGURE 6
www.frontiersin.org

FIGURE 6. Conceptual summary of events determining the particle associated microbial community composition, from the stochastic attachment of present bacteria in the water column (A) through competitive and viral interactions (B) to carbon-quality driven changes in community (C). The sequence of events occurring on the particle renders it a generator of genetic diversity through transduction, conjugation, and uptake of naked DNA (D). Microbial attachment is likely to occur rapidly within hours from particle formation. For as long as the particle received internal and external supply of dissolved organic matter, interactions between organisms will determine the community compositions. Depending of the particle’s sinking rate this is a matter of hours to days. Once the supply of fresh organic matter is depleted and the particle’s labile organic matter is consumed carbon-quality driven community shifts will occur, a process that will take place within days from particle formation.

Author Contributions

MB-I designed the experiments, analyzed the data, and wrote the paper. DI designed the experiments, analyzed the data, and wrote the paper. H-PG wrote the paper.

Funding

This study was funded by the Leibniz Society and all authors were supported by the German Science Foundation (DFG) projects Aquameth and MicroPrime (GR 1540/21-1 and /23-1, respectively). Further support came from the German Academic Exchange Service DAAD (stipend for MB-I).

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.

Acknowledgments

We acknowledge the technical support of Solvig Pinnow.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02569/full#supplementary-material

FIGURE S1 | Composition of the total microbial community (DNA) of single particles and pooled lake particles as obtained from a metagenomic analysis of a sample collected after 2 days of incubation.

FIGURE S2 | A subset of functional modules involved in sugar utilization obtained from the singe particles and pooled lake particles transcriptomes. The shown modules are Level-3 MG-RAST subsystem hierarchy. Abundance refers to percent of reads out of total attributed to sugar utilization.

FIGURE S3 | A subset of functional modules involved in nitrogen metabolism and utilization obtained from the single particles and pooled lake particles transcriptomes. The shown modules are Level-3 MG-RAST subsystem hierarchy. Abundance refers to percent of reads out of total attributed to nitrogen metabolism. Note that subsystem “ZZ gjo need homes” is a valid subsystem name which includes genes with no other subsystem affiliation. In this case, however, the function of genes is associated with nitrogen. The fourth replicate of long incubation single particles (R4) resulted in no transcripts related to nitrogen metabolism.

FIGURE S4 | A subset of functional modules involved in phosphorus metabolism and utilization obtained from the single particles and pooled lake particles transcriptomes. The shown modules are Level-3 MG-RAST subsystem hierarchy. Abundance refers to percent of reads out of total attributed to phosphorus metabolism. The fourth replicate of long incubation single particles (R4) resulted in no transcripts related to phosphorus metabolism.

FIGURE S5 | A subset of functional modules involved in iron metabolism and utilization obtained from the single particles and pooled lake particles transcriptomes. The shown modules are Level-3 MG-RAST subsystem hierarchy. Abundance refers to percent of reads out of total attributed to iron metabolism. The second replicate of short incubation single particles (R2) resulted in no transcripts related to iron metabolism.

FIGURE S6 | Co-occurrence networks of communities associated with single particles at the family level from short- and long-term incubations as calculated using CoNet (Faust and Raes, 2016). Depicted nodes are supported by at least 3 correlation or similarity methods. Node size represents maximal abundance (percent of reads attributed to taxa). Flavobacteriaceae, Moraxellaceae, and Pseudomonadaceae are the only common species between these networks, however, without common associations to other species.

FIGURE S7 | Scanning electron micrograph of an individual particle from Day 1 (A) or Day 5 (B). Phage-like particles were observed on cells resembling Caulobacter (C) as well as on other bacteria (D). Red and white arrows represent potential exit sites of phage particles and attached phages, respectively.

TABLE S1 | Statistical overview of the transcriptome annotation of all individual samples using the extended Trinotate pipeline.

TABLE S2 | Statistical overview of the transcriptome annotation of pooled duplicate samples using the Trinotate pipeline.

TABLE S3 | List of transcripts in which CRISPR spacers were identified. Spacers location and coverage details are provided as well. Spacer sequences are given in Supplementary Data Sheet S4.

TABLE S4 | Transcripts containing matches to polyketide synthethases and non-ribosomal peptide synthetases.

Footnote

  1. ^https://trinotate.github.io/

References

Alldredge, A., Cole, J., and Caron, D. (1986). Production of heterotrophic bacteria inhabiting macroscopic organic aggregates (marine snow) from surface waters. Limnol. Oceanogr. 31, 68–78. doi: 10.4319/lo.1986.31.1.0068

CrossRef Full Text | Google Scholar

Alldredge, A., and Silver, M. (1988). Characteristics, dynamics and significance of marine snow. Prog. Oceanogr. 20, 41–82. doi: 10.1016/0079-6611(88)90053-5

CrossRef Full Text | Google Scholar

Alldredge, A., and Youngbluth, M. (1985). The significance of macroscopic aggregates (marine snow) as sites for heterotrophic bacterial production in the mesopelagic zone of the subtropical Atlantic. Deep Sea Res. Part A 32, 1445–1456. doi: 10.1016/0198-0149(85)90096-2

CrossRef Full Text | Google Scholar

Altschul, S., Gish, W., and Miller, W. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2

CrossRef Full Text | Google Scholar

Amin, S. A., Parker, M. S., and Armbrust, E. V. (2012). Interactions between diatoms and bacteria. Microbiol. Mol. Biol. Rev. 76, 667–684. doi: 10.1128/MMBR.00007-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Arias-Andres, M., Klümper, U., Rojas-Jimenez, K., and Grossart, H.-P. (2018). Microplastic pollution increases gene exchange in aquatic ecosystems. Environ. Pollut. 237, 253–261. doi: 10.1016/J.ENVPOL.2018.02.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Baltar, F., Lindh, M. V., Parparov, A., Berman, T., and Pinhassi, J. (2012). Prokaryotic community structure and respiration during long-term incubations. Microbiologyopen 1, 214–224. doi: 10.1002/mbo3.25

PubMed Abstract | CrossRef Full Text | Google Scholar

Bateman, A., Martin, M. J., O’Donovan, C., Magrane, M., Alpi, E., Antunes, R., et al. (2017). UniProt: the universal protein knowledgebase. Nucleic Acids Res. 45, D158–D169. doi: 10.1093/nar/gkw1099

PubMed Abstract | CrossRef Full Text | Google Scholar

Bidle, K. D., and Azam, F. (2001). Bacterial control of silicon regeneration from diatom detritus: significance of bacterial ectohydrolases and species identity. Limnol. Oceanogr. 46, 1606–1623. doi: 10.4319/lo.2001.46.7.1606

CrossRef Full Text | Google Scholar

Bižić-Ionescu, M., Zeder, M., Ionescu, D., Orliæ, S., Fuchs, B. M., Grossart, H.-P., et al. (2015). Comparison of bacterial communities on limnic versus coastal marine particles reveals profound differences in colonization. Environ. Microbiol. 17, 3500–3514. doi: 10.1111/1462-2920.12466

PubMed Abstract | CrossRef Full Text | Google Scholar

Calvo-Díaz, A., Díaz-Pérez, L., Suárez, L. Á., Morán, X. A. G., Teira, E., and Marañón, E. (2011). Decrease in the autotrophic-to-heterotrophic biomass ratio of picoplankton in oligotrophic marine waters due to bottle enclosure. Appl. Environ. Microbiol. 77, 5739–5746. doi: 10.1128/AEM.00066-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Cezairliyan, B., Vinayavekhin, N., Grenfell-Lee, D., Yuen, G. J., Saghatelian, A., and Ausubel, F. M. (2013). Identification of Pseudomonas aeruginosa phenazines that kill Caenorhabditis elegans. PLoS Pathog. 9:e1003101. doi: 10.1371/journal.ppat.1003101

PubMed Abstract | CrossRef Full Text | Google Scholar

Cordero, O. X., and Datta, M. S. (2016). Microbial interactions and community assembly at microscales. Curr. Opin. Microbiol. 31, 227–234. doi: 10.1016/j.mib.2016.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

D’Ambrosio, L., Ziervogel, K., MacGregor, B., Teske, A., and Arnosti, C. (2014). Composition and enzymatic function of particle-associated and free-living bacteria: a coastal/offshore comparison. ISME J. 8, 2167–2179. doi: 10.1038/ismej.2014.67

PubMed Abstract | CrossRef Full Text | Google Scholar

Datta, M. S., Sliwerska, E., Gore, J., Polz, M. F., Cordero, O. X., Alldredge, A. L., et al. (2016). Microbial interactions lead to rapid micro-scale successions on model marine particles. Nat. Commun. 7:11965. doi: 10.1038/ncomms11965

PubMed Abstract | CrossRef Full Text | Google Scholar

De La Rocha, C. L., and Passow, U. (2007). Factors influencing the sinking of POC and the efficiency of the biological carbon pump. Deep Sea Res. Part II Top. Stud. Oceanogr. 54, 639–658. doi: 10.1016/j.dsr2.2007.01.004

CrossRef Full Text | Google Scholar

Eiler, A., and Bertilsson, S. (2004). Composition of freshwater bacterial communities associated with cyanobacterial blooms in four Swedish lakes. Environ. Microbiol. 6, 1228–1243. doi: 10.1111/j.1462-2920.2004.00657.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Enke, T. N., Leventhal, G. E., Metzger, M., Saavedra, J. T., and Cordero, O. X. (2018). Microscale ecology regulates particulate organic matter turnover in model marine microbial communities. Nat. Commun. 91:2743. doi: 10.1038/s41467-018-05159-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Estrum-Yousef, S. R., Feuerpfeil, P., Schubert, H., and Schumann, R. (2000). Quality of particulate matter and its potential sinking rates in pelagic samples of inshore and Bodden waters of the Southern Baltic Sea. Int. Rev. Hydrobiol. 85, 341–357. doi: 10.1002/(SICI)1522-2632(200004)85:2/3<341::AID-IROH341>3.0.CO;2-H

CrossRef Full Text | Google Scholar

Faust, K., and Raes, J. (2016). CoNet app: inference of biological association networks using Cytoscape. F1000Res. 5:1519. doi: 10.12688/f1000research.9050.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferguson, R. L., Buckley, E. N., and Palumbo, A. V. (1984). Response of marine bacterioplankton to differential filtration and confinement. Appl. Environ. Microbiol. 47, 49–55.

PubMed Abstract | Google Scholar

Finn, R., Clements, J., and Eddy, S. (2011). HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 39, W29–W37. doi: 10.1093/nar/gkr367

PubMed Abstract | CrossRef Full Text | Google Scholar

Finn, R. D., Bateman, A., Clements, J., Coggill, P., Eberhardt, R. Y., Eddy, S. R., et al. (2014). Pfam: the protein families database. Nucleic Acids Res. 42, D222–D230. doi: 10.1093/nar/gkt1223

PubMed Abstract | CrossRef Full Text | Google Scholar

Fox, N. K., Brenner, S. E., and Chandonia, J.-M. (2015). The value of protein structure classification information-Surveying the scientific literature. Proteins 83, 2025–2038. doi: 10.1002/prot.24915

PubMed Abstract | CrossRef Full Text | Google Scholar

Gärdes, A., Iversen, M. H., Grossart, H.-P., Passow, U., and Ullrich, M. S. (2011). Diatom-associated bacteria are required for aggregation of Thalassiosira weissflogii. ISME J. 5, 436–445. doi: 10.1038/ismej.2010.145

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibson, J., Sood, A., and Hogan, D. A. (2009). Pseudomonas aeruginosa-Candida albicans interactions: localization and fungal toxicity of a phenazine derivative. Appl. Environ. Microbiol. 75, 504–513. doi: 10.1128/AEM.01037-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

Grissa, I., Vergnaud, G., and Pourcel, C. (2007). CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 35, W52–W57. doi: 10.1093/nar/gkm360

PubMed Abstract | CrossRef Full Text | Google Scholar

Grossart, H.-P. (1999). Interactions between marine bacteria and axenic diatoms (Cylindrotheca fusiformis, Nitzschia laevis, and Thalassiosira weissflogii) incubated under various conditions in the lab. Aquat. Microb. Ecol. 19, 1–11. doi: 10.3354/ame019001

CrossRef Full Text | Google Scholar

Grossart, H. P., Hietanen, S., and Ploug, H. (2003a). Microbial dynamics on diatom aggregates in Øresund, Denmark. Mar. Ecol. Prog. Ser. 249, 69–78. doi: 10.3354/meps249069

CrossRef Full Text | Google Scholar

Grossart, H.-P., Kiørboe, T., Tang, K., and Ploug, H. (2003b). Bacterial colonization of particles: growth and interactions. Appl. Environ. Microbiol. 69, 3500–3509. doi: 10.1128/AEM.69.6.3500-3509.2003

CrossRef Full Text | Google Scholar

Grossart, H.-P., Levold, F., Allgaier, M., Simon, M., and Brinkhoff, T. (2005). Marine diatom species harbour distinct bacterial communities. Environ. Microbiol. 7, 860–873. doi: 10.1111/j.1462-2920.2005.00759.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Grossart, H.-P., Schlingloff, A., Bernhard, M., Simon, M., and Brinkhoff, T. (2004). Antagonistic activity of bacteria isolated from organic aggregates of the German Wadden Sea. FEMS Microbiol. Ecol. 47, 387–396. doi: 10.1016/S0168-6496(03)00305-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Grossart, H.-P., and Simon, M. (1993). Limnetic macroscopic organic aggregates (Lake Snow): occurrence, characteristics, and microbial dynamics in lake constance. Limnol. Oceanogr. 38, 532–546. doi: 10.4319/lo.1993.38.3.0532

CrossRef Full Text | Google Scholar

Grossart, H.-P., and Simon, M. (1998a). Bacterial colonization and microbial decomposition of limnetic organic aggregates (lake snow). Aquat. Microb. Ecol. 15, 127–140. doi: 10.3354/ame015127

CrossRef Full Text | Google Scholar

Grossart, H.-P., and Simon, M. (1998b). Significance of limnetic organic aggregates (lake snow) for the sinking flux of particulate organic matter in a large lake. Aquat. Microb. Ecol. 15, 115–125. doi: 10.3354/ame015115

CrossRef Full Text | Google Scholar

Grossart, H.-P., and Tang, K. W. (2010). www.aquaticmicrobial.net. Commun. Integr. Biol. 3, 491–494. doi: 10.4161/cib.3.6.12975

PubMed Abstract | CrossRef Full Text | Google Scholar

Ionescu, D., Bizic-Ionescu, M., Khalili, A., Malekmohammadi, R., Morad, M. R., de Beer, D., et al. (2015). A new tool for long-term studies of POM-bacteria interactions: overcoming the century-old Bottle Effect. Sci. Rep. 5:14706. doi: 10.1038/srep14706

PubMed Abstract | CrossRef Full Text | Google Scholar

Ionescu, D., Siebert, C., Polerecky, L., Munwes, Y. Y., Lott, C., Häusler, S., et al. (2012). Microbial and chemical characterization of underwater fresh water springs in the Dead Sea. PLoS One 7:e38319. doi: 10.1371/journal.pone.0038319

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, G. (1980). Phytoplankton growth and zooplankton grazing in oligotrophic oceans. Nature 284, 439–441. doi: 10.1038/284439a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, G. A., and Checkley, D. M. (2011). Particle size distributions in the upper 100m water column and their implications for animal feeding in the plankton. Deep Sea Res. Part I Oceanogr. Res. Pap. 58, 283–297. doi: 10.1016/j.dsr.2010.12.008

CrossRef Full Text | Google Scholar

Jiao, N., and Azam, F. (2011). “Microbial carbon pump and its significance for carbon sequestration in the ocean,” in Microbial Carbon Pump in the Ocean, eds N. Jiao, F. Azam, and S. Sanders (Washington, DC: AAAS), 43–45.

Google Scholar

Kiørboe, T. (2001). Formation and fate of marine snow: small-scale processes with large-scale implications. Sci. Mar. 65, 57–71. doi: 10.3989/scimar.2001.65s257

CrossRef Full Text | Google Scholar

Kiørboe, T., Grossart, H., Ploug, H., Tang, K., and Auer, B. (2004). Particle-associated flagellates: swimming patterns, colonization rates, and grazing on attached bacteria. Aquat. Microb. Ecol. 35, 141–152. doi: 10.3354/ame035141

CrossRef Full Text | Google Scholar

Kiørboe, T., Grossart, H.-P., Ploug, H., and Tang, K. (2002). Mechanisms and rates of bacterial colonization of sinking aggregates. Appl. Environ. Microbiol. 68, 3996–4006. doi: 10.1128/AEM.68.8.3996-4006.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiørboe, T., Tang, K., Grossart, H. P., and Ploug, H. (2003). Dynamics of microbial communities on marine snow aggregates: colonization, growth, detachment, and grazing mortality of attached bacteria. Appl. Environ. Microbiol. 69, 3036–3047. doi: 10.1128/AEM.69.6.3036-3047.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, L., Jing, H., Kataoka, T., Buchwald, C., and Liu, H. (2013). Diversity and spatial distribution of hydrazine oxidoreductase (hzo) gene in the oxygen minimum zone off Costa Rica. PLoS One 8:e78275. doi: 10.1371/journal.pone.0078275

PubMed Abstract | CrossRef Full Text | Google Scholar

LeCleir, G. R., DeBruyn, J. M., Maas, E. W., Boyd, P. W., and Wilhelm, S. W. (2014). Temporal changes in particle-associated microbial communities after interception by nonlethal sediment traps. FEMS Microbiol. Ecol. 87, 153–163. doi: 10.1111/1574-6941.12213

PubMed Abstract | CrossRef Full Text | Google Scholar

Legendre, L., Rivkin, R. B., Weinbauer, M. G., Guidi, L., and Uitz, J. (2015). The microbial carbon pump concept: potential biogeochemical significance in the globally changing ocean. Prog. Oceanogr. 134, 432–450. doi: 10.1016/j.pocean.2015.01.008

CrossRef Full Text | Google Scholar

Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., and Henrissat, B. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 42, D490–D495. doi: 10.1093/nar/gkt1178

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, R. A., and Azam, F. (2001). Antagonistic interactions among marine pelagic bacteria. Appl. Environ. Microbiol. 67, 4975–4983. doi: 10.1128/AEM.67.11.4975-4983.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

McCarthy, J., and Goldman, J. (1979). Nitrogenous nutrition of marine phytoplankton in nutrient-depleted waters. Science 203, 670–672. doi: 10.1126/science.203.4381.670

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonnell, A. M. P., and Buesseler, K. O. (2010). Variability in the average sinking velocity of marine particles. Limnol. Oceanogr. 55, 2085–2096. doi: 10.4319/lo.2010.55.5.2085

CrossRef Full Text | Google Scholar

Meyer, F., Paarmann, D., and D’Souza, M. (2008). The metagenomics RAST server–a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi: 10.1186/1471-2105-9-386

PubMed Abstract | CrossRef Full Text | Google Scholar

Moran, M. A., Satinsky, B., Gifford, S. M., Luo, H., Rivers, A., Chan, L.-K., et al. (2013). Sizing up metatranscriptomics. ISME J. 7, 237–243. doi: 10.1038/ismej.2012.94

PubMed Abstract | CrossRef Full Text | Google Scholar

Moriya, Y., Itoh, M., Okuda, S., Yoshizawa, A. C., and Kanehisa, M. (2007). KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 35, W182–W185. doi: 10.1093/nar/gkm321

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller-Niklas, G., and Schuster, S. (1994). Organic content and bacterial metabolism in amorphous aggregations of the northern Adriatic Sea. Limnol. Oceanogr. 39, 58–68. doi: 10.4319/lo.1994.39.1.0058

CrossRef Full Text | Google Scholar

Nowald, N., Fischer, G., Ratmeyer, V., Iversen, M., Reuter, C., and Wefer, G. (2009). “In-situ sinking speed measurements of marine snow aggregates acquired with a settling chamber mounted to the Cherokee ROV,” in Proceedings of the OCEANS 2009-EUROPE, (Bremen: IEEE), 1–6. doi: 10.1109/OCEANSE.2009.5278186

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Peduzzi, P., and Weinbauer, M. G. (1993). Effect of concentrating the virus-rich 2-200-nm size fraction of seawater on the formation of algal floes (marine snow). Limnol. Oceanogr. 38, 1562–1565. doi: 10.4319/lo.1993.38.7.1562

CrossRef Full Text | Google Scholar

Prézelin, B., and Alldredge, A. (1983). Primary production of marine snow during and after an upwelling event. Limnol. Oceanogr. 28, 1156–1167. doi: 10.4319/lo.1983.28.6.1156

CrossRef Full Text | Google Scholar

Proctor, L., and Fuhrman, J. (1992). Mortality of marine-bacteria in response to enrichments of the virus size fraction from seawater. Mar. Ecol. Prog. Ser. 87, 283–293. doi: 10.3354/meps087283

CrossRef Full Text | Google Scholar

Provvedi, R., and Dubnau, D. (1999). ComEA is a DNA receptor for transformation of competent Bacillus subtilis. Mol. Microbiol. 31, 271–280. doi: 10.1046/j.1365-2958.1999.01170.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ratmeyer, V., and Wefer, G. (1996). A high resolution camera system (ParCa) for imaging particles in the ocean: system design and results from profiles and a three-month deployment. J. Mar. Res. 54, 589–603. doi: 10.1357/0022240963213565

CrossRef Full Text | Google Scholar

Riemann, L., and Grossart, H.-P. (2008). Elevated lytic phage production as a consequence of particle colonization by a marine Flavobacterium (Cellulophaga sp.). Microb. Ecol. 56, 505–512. doi: 10.1007/s00248-008-9369-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Schweitzer, B., Huber, I., Amann, R., Ludwig, W., and Simon, M. (2001). Alpha-and Betaproteobacteria control the consumption and release of amino acids on lake snow aggregates. Appl. Environ. Microbiol. 67, 632–645. doi: 10.1128/AEM.67.2.632-645.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Shanks, A. L., and Trent, D. (1979). Marine snow: microscale nutrient. Limnol. Oceanogr. 24, 850–854. doi: 10.4319/lo.1979.24.5.0850

CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Sherr, E., Sherr, B., and Sigmon, C. (1999). Activity of marine bacteria under incubated and in situ conditions. Aquat. Microb. Ecol. 20, 213–223. doi: 10.3354/ame020213

CrossRef Full Text | Google Scholar

Simon, M., Grossart, H.-P., Schweitzer, B., and Ploug, H. (2002). Microbial ecology of organic aggregates in aquatic ecosystems. Aquat. Microb. Ecol. 28, 175–211. doi: 10.3354/ame028175

CrossRef Full Text | Google Scholar

Smith, D. C., Simon, M., Alldredge, A. L., and Azam, F. (1992). Intense hydrolythic activity on marine aggregates and implications for rapid particle dissolution. Nature 359, 139–142. doi: 10.1038/359139a0

CrossRef Full Text | Google Scholar

Smith, M. W., Zeigler Allen, L., Allen, A. E., Herfort, L., and Simon, H. M. (2013). Contrasting genomic properties of free-living and particle-attached microbial assemblages within a coastal ecosystem. Front. Microbiol. 4:120. doi: 10.3389/fmicb.2013.00120

PubMed Abstract | CrossRef Full Text | Google Scholar

Smriga, S., Fernandez, V. I., Mitchell, J. G., and Stocker, R. (2016). Chemotaxis toward phytoplankton drives organic matter partitioning among marine bacteria. Proc. Natl. Acad. Sci. U.S.A. 113, 1576–1581. doi: 10.1073/pnas.1512307113

PubMed Abstract | CrossRef Full Text | Google Scholar

Sorek, R., Kunin, V., and Hugenholtz, P. (2008). CRISPR–a widespread system that provides acquired resistance against phages in bacteria and Archaea. Nat. Rev. Microbiol. 6, 181–186. doi: 10.1038/nrmicro1793

PubMed Abstract | CrossRef Full Text | Google Scholar

Staub, R. (1961). Study of nutrition physiology and autecology of planktic blue-green alga Oscillatoria rubescens DC. Schweiz. Z. Hydrol. 23, 82–198.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Stocker, R. (2012). Marine microbes see a sea of gradients. Science 338, 628–633. doi: 10.1126/science.1208929

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, K. W., Dziallas, C., and Grossart, H.-P. (2011). Zooplankton and aggregates as refuge for aquatic bacteria: protection from UV, heat and ozone stresses used for water treatment. Environ. Microbiol. 13, 378–390. doi: 10.1111/j.1462-2920.2010.02335.x

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, G. P., Kin, K., and Lynch, V. J. (2012). Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 131, 281–285. doi: 10.1007/s12064-012-0162-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Walsh, D. A., Lafontaine, J., and Grossart, H.-P. (2013). “On the eco-evolutionary relationships of fresh and salt water bacteria,” in Lateral Gene Transfer in Evolution, ed. U. Gophna (New York, NY: Springer), 443.

Google Scholar

Weber, T., Blin, K., Duddela, S., Krug, D., Kim, H. U., Bruccoleri, R., et al. (2015). antiSMASH 3.0—a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res. 43, W237–W243. doi: 10.1093/nar/gkv437

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinbauer, M. G. (2004). Ecology of prokaryotic viruses. FEMS Microbiol. Rev. 28, 127–181. doi: 10.1016/j.femsre.2003.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilke, A., Bischof, J., Gerlach, W., Glass, E., Harrison, T., Keegan, K. P., et al. (2016). The MG-RAST metagenomics database and portal in 2015. Nucleic Acids Res. 44, D590–D594. doi: 10.1093/nar/gkv1322

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, D., Pethica, R., Zhou, Y., Talbot, C., Vogel, C., Madera, M., et al. (2009). SUPERFAMILY–sophisticated comparative genomics, data mining, visualization and phylogeny. Nucleic Acids Res. 37, D380–D386. doi: 10.1093/nar/gkn762

PubMed Abstract | CrossRef Full Text | Google Scholar

Wurzbacher, C., Rösel, S., Rychła, A., and Grossart, H.-P. (2014). Importance of saprotrophic freshwater fungi for pollen degradation. PLoS One 9:e94643. doi: 10.1371/journal.pone.0094643

PubMed Abstract | CrossRef Full Text | Google Scholar

Wurzbacher, C., Warthmann, N., Bourne, E., Attermeyer, K., Allgaier, M., Powell, J. R., et al. (2016). High habitat-specificity in fungal communities in oligo-mesotrophic, temperate Lake Stechlin (North-East Germany). MycoKeys 16, 17–44. doi: 10.3897/mycokeys.16.9646

CrossRef Full Text | Google Scholar

Ziemert, N., Podell, S., Penn, K., Badger, J. H., Allen, E., and Jensen, P. R. (2012). The natural product domain seeker NaPDoS: a phylogeny based bioinformatic tool to classify secondary metabolite gene diversity. PLoS One 7:e34064. doi: 10.1371/journal.pone.0034064

PubMed Abstract | CrossRef Full Text | Google Scholar

Zobell, C., and Anderson, D. (1936). Observations on the multiplication of bacteria in different volumes of stored sea water and the influence of oxygen tension and solid surfaces. Biol. Bull. 71, 324–342. doi: 10.2307/1537438

CrossRef Full Text | Google Scholar

Keywords: particle-associated bacteria, microbial communities, inter- and intra-species interactions, antagonism, phage, transcriptome

Citation: Bižić-Ionescu M, Ionescu D and Grossart H-P (2018) Organic Particles: Heterogeneous Hubs for Microbial Interactions in Aquatic Ecosystems. Front. Microbiol. 9:2569. doi: 10.3389/fmicb.2018.02569

Received: 13 March 2018; Accepted: 08 October 2018;
Published: 26 October 2018.

Edited by:

Anne Michelle Wood, University of Oregon, United States

Reviewed by:

Brandy Marie Toner, University of Minnesota Twin Cities, United States
David A. Walsh, Concordia University, Canada

Copyright © 2018 Bižić-Ionescu, Ionescu and Grossart. 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) and the copyright owner(s) 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: Mina Bižić-Ionescu, mbizic@igb-berlin.de Hans-Peter Grossart, hgrossart@igb-berlin.de