Environmental Factors Support the Formation of Specific Bacterial Assemblages on Microplastics

While the global distribution of microplastics (MP) in the marine environment is currently being critically evaluated, the potential role of MP as a vector for distinct microbial assemblages or even pathogenic bacteria is hardly understood. To gain a deeper understanding, we investigated how different in situ conditions contribute to the composition and specificity of MP-associated bacterial communities in relation to communities on natural particles. Polystyrene (PS), polyethylene (PE), and wooden pellets were incubated for 2 weeks along an environmental gradient, ranging from marine (coastal Baltic Sea) to freshwater (waste water treatment plant, WWTP) conditions. The associated assemblages as well as the water communities were investigated applying high-throughput 16S rRNA gene sequencing. Our setup allowed for the first time to determine MP-dependent and -independent assemblage factors as subject to different environmental conditions in one system. Most importantly, plastic-specific assemblages were found to develop solely under certain conditions, such as lower nutrient concentration and higher salinity, while the bacterial genus Erythrobacter, known for the ability to utilize polycyclic aromatic hydrocarbons (PAH), was found specifically on MP across a broader section of the gradient. We discovered no enrichment of potential pathogens on PE or PS; however, the abundant colonization of MP in a WWTP by certain bacteria commonly associated with antibiotic resistance suggests MP as a possible hotspot for horizontal gene transfer. Taken together, our study clarifies that the surrounding environment prevailingly shapes the biofilm communities, but that MP-specific assemblage factors exist. These findings point to the ecological significance of specific MP-promoted bacterial populations in aquatic environments and particularly in plastic accumulation zones.

Introduction

Microplastics (MP) can be found worldwide in aquatic ecosystems, reaching from populated to pristine, from marine to freshwater regions. Specific MP accumulation zones emerge due to oceanographic currents, such as in the subtropical gyres (Cózar et al., 2014) or, due to anthropogenic point sources, in industrial agglomeration areas (Mani et al., 2015). The ecological implications of this pollution in the different environmental systems can so far not be estimated.

The interaction of MP and microbial communities in marine and freshwater is especially underrepresented in the scientific literature. It is, however, well-established that plastic biofilms contain several “usual suspects,” which are abundant in biofilms on many surfaces in aquatic ecosystems. In the marine environment, these include typical initial colonizers of the family Rhodobacteraceae (Dang et al., 2008; Zettler et al., 2013) and the algae group Bacillariophyta (Carson et al., 2013; Oberbeckmann et al., 2014; Reisser et al., 2014). Members of the Campylobacteraceae were found to be abundant on MP sampled in an urban river in Chicago (USA) (McCormick et al., 2014) as well as on MP incubated in coastal sediment from the Humber Estuary (UK) (Harrison et al., 2014). Following one assumption by Maso et al. (2003), that a harmful algae bloom species was dispersed in the Mediterranean while hitchhiking on macroplastic debris, the question emerged whether also MP can serve as a vector for harmful microorganisms. In 2013, members of the potentially pathogenic genus Vibrio were detected to be abundant on a polypropylene (PP) particle from the North Atlantic (Zettler et al., 2013). Since then, MP from WWTPs and MP passing the guts of aquatic organisms have been hypothesized to play a particular role as vector for pathogenic microorganisms (McCormick et al., 2014; Oberbeckmann et al., 2015). Studies on Arenicola marina and Mytilus edulis, however, did not reveal a distinct enrichment of potential pathogens on PS (Kesy et al., 2016, 2017).

Pioneering studies in the field of MP biofilms reveal that biogeography plays an important role for the composition of MP colonizing communities (Oberbeckmann et al., 2014; Amaral-Zettler et al., 2015). Not much, however, is known about the plastic-colonizing processes on a local to regional scale, and a detailed understanding of the factors which are shaping the composition, activity, and potential specificity of MP-associated assemblages is lacking. No investigation of MP-colonization comprised marine as well as freshwater environments, and only very limited previous research enabled a direct comparisons of MP-attached assemblages with those colonizing natural particles, e.g., marine snow.

Our study addressed these knowledge gaps to be able to reliably evaluate the potential of MP to be a vector of distinct microbial assemblages or even pathogens and to better understand the characteristics and dynamics of microbial MP colonizing communities. Our hypothesis was that environmental factors play a major role for the colonization of natural and synthetic surfaces alike, but that a specific plastic-associated microbiome would develop, different from the microbiomes associated with natural particles. We designed an in situ incubation experiment, where PE, PS, and wooden pellets were exposed for 2 weeks at 7 different stations in the coastal Baltic Sea, in the estuary of the river Warnow (Rostock, Germany), and in the effluent basins of a WWTP. Corresponding water communities (free-living and attached to naturally occurring particulate material) were sampled as well. The locations covered different environmental conditions along a salinity and nutrient gradient and included different anthropogenic impacts.

Experimental Procedures

Please see Additional file 1 for more detailed information.

Experimental Design and Setup

HDPE (ExxonMobil^TM HDPE HTA 108, density 0.961 g/cm³) and PS (BASF Polystyrole 143E, density 1.04 g/cm³) pellets (both ø 3 mm) were incubated for 14 days (19./20.08.−03./04.09.2014) at five stations in the estuary of the river Warnow and in the Baltic Sea, Germany (Figure 1). Station 1 was in Heiligendamm (54.146 N 11.843 E), station 2 in the Alter Strom in Warnemünde (54.181 N, 12.087 E), station 3 in the Marine Science Center in Hohe Düne (54.183 N, 12.095 E), station 4 in the Unterwarnow close to a WWTP discharge (54.106 N, 12.096 E), and station 5 in the Unterwarnow close to the city center of Rostock (54.097 N, 12.151 E). As non-plastic control, wooden pellets were incubated. Nine incubators were built in-house and installed in the water column at each station (n = 3 per material). Each incubator was covered with a 500 μm pore-sized gauze to keep the pellets inside while assuring sufficient flow-though. Considering the different densities of the materials, 54 g of PS, 50 g of HDPE, or 35 g of wood were added to the incubators. Recovered plastic and wood pellets were rinsed thrice with sterile station water, shock frozen with liquid nitrogen, and stored at −80°C until further processing. To determine whether plastic and wood biofilm communities are unique from surrounding microbial communities, background water samples were collected in triplicate at the end of the 2-week exposure experiment from every station. The >3 μm (1 l, “particle-attached”) and the 3–0.2 μm (300–500 ml, “free-living”) fractions were collected by serial filtration. Filters were shock frozen and stored at −80°C until further analyses.

FIGURE 1

Figure 1. Map section of Northern Europe, with box highlighting sampling area (Top), and individual stations 1-5 (Bottom) near Rostock, Germany, from the coastal Baltic Sea to the river Warnow.

The experiment was repeated as described above at two stations in a WWTP (28.4.−12.05.2015, name and location remain anonymous as requested by operator). One incubation (station 6) took place in a basin following three treatment steps (a mechanical pre-treatment, an activated sludge, a biological treatment), the other incubation (station 7) following an additional biological filtration, just before the WWTP outflow. Most German WWTPs have a similar setup, releasing the treated water after these 3 or 4 stages.

During both incubation experiments, environmental metadata were collected in triplicates at the beginning and end of the incubation (Table 2). Temperature and salinity were measured with a portable meter (HQ40d, Hach), while ${\text{PO}}_4^{3-}$, ${\text{NO}}_3^{-}$, ${\text{NO}}_2^{-}$, ${\text{NH}}_4^{+}$, and SiO₂ were determined using a colorimetric method (Grasshoff et al., 1999).

DNA Isolation

DNA was extracted from incubated plastic pellets (20 pellets per sample), wood pellets (1/2 teaspoon), and filters. The procedure involved a treatment with SLS buffer [20% w/v SLS], bead beating and extraction with 1 vol. phenol/chloroform, followed by precipitation using 70% ethanol. Pre-incubated control pellets were subjected to the same procedure.

Amplification and Library Generation

Dual-indexing was used to generate a barcoded MiSeq library of tag sequences for each sample of extracted DNA (Kozich et al., 2013). Primers targeting the 16S V4 region were used: 515F (5′ GTGCCAGCMGCCGCGGTAA 3′; (Giovannoni et al., 1988; Caporaso et al., 2011) and 806R (5′ GGACTACHVGGGTWTCTAAT 3′; (Caporaso et al., 2011). Due to low DNA concentrations, a nested PCR approach [with primer pair 27F: 5′ AGAGTTTGATCCTGGCTCAG 3′ (Lane, 1991) and 1492R: 5′ GGTTACCTTGTTACGACTT 3′ Lane, 1991] was applied for samples from station 1–5, while direct amplification could be applied for the samples from station 6 and 7. Sterile water control samples revealed consistently negative results in the nested as well as in the direct amplification approach. Likewise, pre-incubated wood and plastic pellet samples did not lead to a successful amplification using neither the nested nor the direct PCR approach.

Sequence Data Processing

Quality filtering of reads (adapter, barcode, primer clipping, permitted length = 250–275 bp, max. number of ambiguous bases per sequence = 0, max. number of homopolymers per sequence = 8), taxonomy assignment (Wang classification, reference database SSURef_123_SILVA, required bootstrap value ≥ 85%) and picking of operational taxonomic units, OTUs (label = 0.03), were carried out using Mothur (Schloss et al., 2009). Chloroplasts, mitochondria, eukaryotes and unknown sequences were removed. OTUs with a total abundance of ≤ 3 were excluded from downstream analyses. All raw sequence files, including sequencing controls, are available from the NCBI Short Read Archive (SRA) database (BioProject PRJNA338729; BioSamples SAMN05567827-SAMN05567879, SAMN06883266, SAMN06883268).

Alpha and Beta Diversity and Statistical Analysis

Chao1 (species richness; Chao, 1984), Shannon H (diversity; Shannon, 1948) and Pielou (evenness; Pielou, 1969) indices were calculated based on a subset of 500 reads per samples (100 iterations). Alpha diversity analyses were carried out with RStudio (v 3.0.2 RStudio Team, 2015) and the vegan package (Oksanen et al., 2015). Using SigmaPlot (v 10, Systat Software GmbH), results of alpha diversity analyses were visualized, and it was tested for significant differences between richness, diversity, and evenness between the three datasets (t-test).

For beta diversity analyses, samples with less than 10,000 reads were excluded. Relative abundances of OTUs were calculated and square root transformed. The dataset “Baltic” consisted of samples from stations 2 and 3, the “Warnow” dataset of samples from stations 4 and 5, and the “WWTP” dataset of samples from stations 6 and 7. This segmentation of the data was based on an hierarchical clustering of the environmental factors (Figure 3). Station 1 was omitted for the dataset segmentation to assure a similar number of samples per dataset. The “plastic” dataset contained all data on PE and PS associated communities. In our study, the factor “environment” refers to the different sample locations (1–7), the factor “source” refers to the different sample types (PE, PS, wood, free-living, and particle attached water fraction).

Permutational multivariate analysis of variance (PERMANOVA) (Anderson, 2001) was applied to test for significant differences between environments and sources based on Bray-Curtis similarity matrices. P-values were obtained using type III sums of squares and 999 permutations. When analyzing solely subsets of the whole dataset, Monte Carlo permutations were calculated due to the lower number of possible permutations. To test whether significant (p < 0.05) PERMANOVA results were based on location or dispersion effects, homogeneity of dispersion (PERMDISP) was applied based on calculated distances to centroids (Anderson and Walsh, 2013). To disclose false discovery rates, Benjamini-Hochberg corrections (Benjamini and Hochberg, 1995) for multiple testing was applied using RStudio v 3.3.1 (RStudio Team, 2015). To visualize the patterns of samples regarding environment and source, PCO was performed (Gower, 1966).

The environmental data were normalized and used to construct a Euclidian distance matrix that was subjected to hierarchical clustering (group average), omitting SiO₂ due to a missing value. CLUSTER, PERMANOVA, PERMDISP, and PCO were calculated with the software PRIMER 6 including the PERMANOVA+ add-on package (PRIMER-E, UK; Clarke and Gorley, 2006; Anderson et al., 2008).

To identify taxa that discriminated (i) the plastic-attached from the other communities, and (ii) the PS from the PE associated communities, the linear discriminant analysis effect size method (LEfSe; Segata et al., 2011) was used. This was performed with the LEfSe online tool in the Galaxy framework, using all default settings for data formatting and LDA effect size. Instead of using a square root transformation, non-transformed relative abundances were used and a LEfSe-internal per-sample normalization was performed. Factor “source” was set as class, and the strategy “all-against-all” was applied. Solely OTUs with a mean (across the whole dataset) relative abundance of ≥0.1% were considered.

A phylogenetic tree was built with the 175 OTUs displaying a mean relative abundance of ≥0.1%. The tree was constructed using mothur, and visualized using the online tool iTOL version 3.2.4 (Letunic and Bork, 2016).

Quantitative PCR (qPCR)

Standard DNA for Vibrio-targeting qPCR originated from an in house culture of V. aestuarianus and was extracted using the QIAGEN^® QIAamp^® DNA Mini Kit. A 114 bp long DNA fragment was amplified using the primer pair 567F (GGCGTAAAGCGCATGCAGGT) and 680R (GAAATTCTACCCCCCTCTACAG) (Thompson et al., 2004), applying an annealing temperature of 60.2°C. Cloning of the purified PCR products into the vector took place with the StrataClone Blunt PCR Cloning Kit (Agilent Technologies, USA) and extraction of the vector from colonies with the QIAGEN^® Plasmid Mini Kit (QIAGEN, DE). The samples were analyzed in three qPCR runs, each run including three standard samples. The mean values of the run conditions are listed in the following: slope = 3.155 ± 0.062, PCR efficiency = 107.53 ± 2.97, γ intercept = 45.609 ± 1.440, R² = 0.998±0.002, Ct cut-off = 30.7.

Results

The Environment Determines the Substrate Specificity of Microbial Colonization

To verify the overall statistical relevance of the factors “environment” and “source,” PERMANOVA was carried out, and the community patterns were visualized using principal coordinates analysis (PCO, Figure 2).

FIGURE 2

Figure 2. PCO displaying variation in microbial community composition in the different sample sources water (free-living fraction=water_0.2; particle-attached fraction=water_3), wood, PE, PS. Labels refer to sampling stations 1-7. Plot is based on a Bray-Curtis similarity matrix of square root transformed relative abundances of OTUs (16S rRNA gene data).

The factor “environment” differentiated the communities significantly (global p = 0.001; 20 of 21 pairwise p < 0.05, Additional file 2). The factor “source” (PE, PS, wood, water free-living, and particle-attached) also differentiated the communities significantly (global p = 0.001; 9 of 10 pairwise p < 0.05, Table 1). Solely the communities on PE and PS did not differ significantly. Benjamini-Hochberg corrections supported the significant p-values from PERMANOVA. All pairwise tests comparing samples from station 6 or 7 with samples from station 1 to 5 displayed significant PERMDISP results. Likewise, pairwise tests involving the particle-attached water fraction (“3”) displayed significant or near-significant PERMDISP results.

TABLE 1

Table 1. Results of global and pairwise PERMANOVA, Benjamini-Hochberg correction, and corresponding PERMDISP of microbial communities from the different sample sources PE, PS, wood (L), particle-attached (3), and free-living water fraction (0.2).

PCO of all samples displayed a gradient, which was formed by salinity, temperature and, negatively correlating, nutrient concentrations (Figure 2). The samples roughly line up along this gradient according to station. However, the first two axes solely amount to 38% of the variation within the communities. Samples from the coastal Baltic and Warnow form two joint clusters, one comprising the plastic and wood-associated communities, the other the water communities. All samples from the WWTP form another cluster.

Hierarchical clustering of the environmental data indicated similar environmental conditions at stations 1 to 3 (coastal Baltic Sea) and at stations 4 and 5 (Warnow estuary) (Figure 3). Station 6 and 7 (WWTP) cluster away from the other stations, while displaying different conditions between each other. Station 1, 2, and 3 featured marine/brackish conditions (salinity 13.5–14.7) and lowest nutrient concentrations (e.g., 0.09–3.84 μmol ${\text{NO}}_3^{-}$/l) (Figure 1, Table 2). Station 4 and 5 displayed lower salinities (6.5–10.0, brackish conditions) and increased nutrient concentrations (e.g., 6.03–48.76 μmol ${\text{NO}}_3^{-}$/l). Station 6 and 7 had freshwater (salinity of < 0.5) and high nutrient levels (e.g., 701.80–751.35 μmol ${\text{NO}}_3^{-}$/l). The water temperature in the sampling period was similar at stations 1–5 (17.0–18.3°C), while lower at stations 6 and 7 (14.7–15.2°C).

FIGURE 3

Figure 3. Hierarchical cluster (group average) based on Euclidian distances of environmental data (temperature, salinity, ${\text{PO}}_4^{3-}$, ${\text{NO}}_3^{-}$, ${\text{NO}}_2^{-}$, ${\text{NH}}_4^{+}$). Cluster displays similarity of environmental conditions at stations 1 to 7.

TABLE 2

Table 2. Temperature, salinity and nutrient concentrations (${\text{PO}}_4^{3-}$, ${\text{NO}}_3^{-}$, ${\text{NO}}_2^{-}$, ${\text{NH}}_4^{+}$, and SiO₂) of the water at stations 1–7 at the start and end dates of the in situ experiment.

Based on the clustering of the environmental data, the sequencing data were partitioned into “Baltic” (station 2 and 3), “Warnow” (station 4 and 5), and “WWTP” (station 6 and 7) datasets. Pairwise PERMANOVA (Monte Carlo) revealed a significant separation between all sample sources in the “Baltic” dataset (Additional file 3). This includes the communities attached to plastic compared to wood (p = 0.001) as well as the communities colonizing PE compared to PS (p = 0.047). The “Warnow” dataset displayed no significant separation between the communities attached to PE compared to PS, and to wood compared to PS. All other sample sources, however, differed significantly. In the “WWTP” dataset, no significant differences were found between either of the attached communities. Solely the free-living and attached communities differed significantly. The different levels of substrate specificity of the microbial communities in the three datasets are visualized in Figure 4. All significant PERMANOVA results were supported by Benjamini-Hochberg corrections, except the differentiation between communities associated with wood and PE in the “Warnow” dataset. Dispersion effects could be detected in some cases (Additional file 3).

FIGURE 4

Figure 4. PCO displaying variation in microbial community composition in the different sample sources water (free-living fraction=water_0.2; particle-attached fraction=water_3), wood, PE, PS. Labels refer to sampling stations. Circles indicate 20 (straight line) and 40% (dashed line) of similarity. Plots are based on Bray-Curtis similarity matrices of square root-transformed relative abundances of OTUs (16S rRNA gene data), displaying the “Baltic” (A), the “Warnow” (B), the “WWTP” (C) dataset. Also given are P-values of pairwise PERMANOVA (Monte Carlo) from the tests (1) wood vs. PE, (2) wood vs. PS, (3) PE vs. PS, with significant values (p < 0.05) in bold.

Based on alpha diversity analyses, assemblages from the “WWTP” dataset displayed significantly higher richness, diversity and evenness compared to the other stations (Figure 5).

FIGURE 5

Figure 5. Box plots illustrating the Chao1, Shannon H and Pielou E indices of datasets “Baltic,” “Warnow,” “WWTP”; *indicates significant differences (p < 0.01) between datasets.

Partitioning of MP-Associated Microbial Communities

Some families displayed high relative abundances at all stations in the PE and PS-associated as well as other communities: Flavobacteriaceae (esp. the genus Flavobacterium), Rhodobacteraceae, Methylophilaceae (Methylotenera), and Planctomycetaceae (Planctomyces, Pirellula; decreasing from station 1 to 7) (Figure 6 Additional file 4). The family Hyphomonadaceae (Figure 6 OTU 155-158), in particular the genus Hyphomonas, was found abundantly at stations 1–4 and 7, being consistently more abundant on both plastic surfaces (up to 6.3% relative abundance per sample) than on natural surfaces. Within the family Planctomycetaceae, one OTU assigned to the genus Blastopirellula was very abundant solely on plastic (up to 6.8%; Figure 6 OTU 96). At the Baltic and Warnow stations, the family Erythrobacteraceae (genus Erythrobacter, Figure 6 OTU 152) was abundant on PS and particularly on PE (up to 4.1%). One OTU (Figure 6 OTU 149), classified as Sphingopyxis (Sphingomonadaceae), displayed high relative abundances solely in the plastic-assemblages, particularly at station 7 (up to 8%).

FIGURE 6

Figure 6. Phylogenetic tree of the 175 OTUs (16S rRNA data based) displaying an average relative abundance of ≥0.1% across the whole dataset. The colors of the branches refer to different bacterial phyla. The outer rings represent different sample sources: particle-attached water community (PAW), wood-, PE-, and PS-associated assemblages. The colors in the rings indicate the mean relative abundance of the OTUs in the individual sample sources across all stations. Scale bar represents 10 nucleotide substitutions per 100 nucleotides. Families with ≥ 3 OTUs are listed below. For more taxonomic details of OTUs see Additional file 4. 1–36: Flavobacteriaceae; 41–45: Cryomorphaceae; 52–58, 60: Saprospiraceae; 68, 69, 71, 73–75: Verrucomicrobiaceae; 78–80: Microbacteriaceae; 82–84: Cyanobacteria Family I; 86–96: Planctomycetaceae; 106, 107, 118, 119: Cellvibrionaceae; 122–130: Comamonadaceae; 131, 134–136: Rhodocyclaceae; 144–146: Campylobacteraceae; 150, 152, 153: Erythrobacteraceae; 155–158: Hyphomonadaceae; 159–175: Rhodobacteraceae.

LEfSe analysis (considering solely OTUs with a mean relative abundance ≥0.1%) confirmed the significant differential abundance in the plastic-assemblages of the above described genera Hyphomonas, Erythrobacter, Blastopirellula, and Sphingopyxis (Additional file 5). Other taxa were also significantly more abundant on plastic, despite not showing overall high relative abundances. These include the family Phyllobacteriaceae and the genus Tenacibaculum (Flavobacteriaceae). Likewise, the communities were tested toward their preference for PE or PS (Additional file 6). Most differential taxa between polymers were found in the “Baltic” (N = 37) dataset, and more differential abundant taxa were detected for PS- (N = 34) than PE- (N = 12) assemblages. Members of the families Erythrobacteraceae were in higher abundances associated with PE, while members of the families Verrucomicrobiaceae were significantly more abundant on PS.

Occurrence of MP-Associated Potentially Pathogenic Bacteria

No increased relative abundances of potential pathogens on plastic compared to other sample types were detected. The genus Vibrio (Vibrionaceae) was found in plastic-attached communities, but LEfSe revealed significant higher abundance of this genus in wood-associated samples (LDA score = 4.7). The 16S rRNA gene barcode data were supported by qPCR. Of 65 samples tested toward the occurrence of Vibrio sp., 27 gave positive results, with members of the genus Vibrio being found at all tested stations and in all tested sources (Additional file 7). The particle-attached water communities displayed the most positive Vibrio sp. samples (N = 10), while the highest Vibrio abundances were detected for wood samples (1.30*10⁷-2.16*10⁸ gene copies/ng DNA, station 2). Solely two positive samples each were found for the plastic assemblages, with Vibrio concentrations of up to 6.95*10⁴ copies/ng at PS at station 4. Wood and water communities from station 4, however, displayed higher Vibrio abundances (1.16–3.28*10⁵ copies/ng).

Some OTUs from the plastic-associated sequences were assigned to members of the family Enterobacteriaceae. The highest relative abundance of Enterobacteriaceae in plastic samples was found on PS at station 4 and 6 (0.13 and 0.12%), with higher abundances of this family in water samples at the same stations (up to 1.7%).

Discussion

MP as Specific Vector Relevant under Certain Environmental Conditions

The present study sheds light on the relevance of MP, in comparison with natural particles, as a vector for specific bacterial communities under different environmental conditions. Our field experiment revealed that plastic-associated assemblages show spatial variation, in accordance with previous studies (Hoellein et al., 2014; Oberbeckmann et al., 2014; Amaral-Zettler et al., 2015). Here, we now demonstrate that we can find a MP-specific microbiome in marine waters, after McCormick and coauthors established this for freshwater habitats already in 2014. Our study reveals for the first time that the degree of specificity depends on the ambient environmental conditions, and that we can refer to the term “plastisphere” (Zettler et al., 2013) solely in certain environmental scenarios. The scenario “lower nutrient levels, higher salinity, coastal system” led to substrate specificity in our study. In areas with lower nutrient concentrations (Baltic stations) we could not only observe a complete differentiation between plastic, wood and particle-attached water communities but also significant differences between the PE- and PS-assemblages. At the stations with higher nutrient concentrations, especially within the WWTP, we could not detect major differences between the communities associated with different substrates. The fungal community from the same experimental setup showed similar patterns, as published recently by Kettner et al. (2017). We detected, using PERMDISP analyses, dispersion factors to play a role in the differentiation between some of the sample groups. Linear discriminant analysis and visualization via PCO, however, support that location factors play at least an equivalent role.

Rummel et al. (2017) summarized substrate parameters relevant for microbial colonization, e.g., surface roughness and hydrophobicity. Considering the significant differences between PE, PS, and wood biofilms at the coastal stations, these factors assumingly did play a role. The surface parameters, however, were overall secondary to the environmental influence. Potentially, surface parameters play mainly a role in early colonization stages, with spatial factors increasing in importance over time. Alternatively, environmental factors have a stronger influence on shaping MP assemblages starting with the initial attachment, e.g., by determining the conditioning film. In order to clarify these early attachment processes, future experiments on MP biofilms might include the colonization phase of the first few hours to days. In all aquatic systems biofilm formation, architecture and stability depend on complex interactions with biotic and abiotic factors (Gerbersdorf and Wieprecht, 2015). A study by De Tender et al. (2015) suggested an influence of factors as salinity, temperature and oxygen content, on microbial colonization of macroplastics. The affirmation of these complex influences on MP biofilms can now be provided with our data. Several parameters including nutrient concentrations, salinity, and temperature contributed to shaping the microbial communities on PE and PS.

Low nutrient levels trigger the attachment to surfaces in many bacterial species, while under high nutrient levels biofilm formation seems less advantageous (Stanley and Lazazzera, 2004). Zobell postulated already in 1943 that biofilm formation plays a major role in nutrient-low waters, and mentioned also plastic to be surface active (Zobell, 1943). Over 70 years later, we confirm his findings with regard to MP. Presumably, the microbial colonization at the coastal stations with lower nutrient availability was a directed process, specific to each surface type (e.g., hard vs. soft substratum), and involving microorganisms with a more specialized metabolism. The more nutrients that are available, the quicker a conditioning film and primary biofilm can develop, the faster a secondary, less substrate specific, biofilm can be established. Datta et al. (2016) investigated micro-scale successions on marine particles, and described three phases of colonization. Specifically, the third phase is dominated by secondary consumers and displays the highest biofilm diversity. The term “Fast-but-inefficient metabolism strategy” for bacteria was coined by Polz and Cordero (2016) based on a study by Roller et al. (2016): under high nutrient flux, ribosomal RNA operon numbers increase, followed by higher biomass production but lower carbon use efficiency. Valuable metabolic intermediates are excreted, which in a biofilm can be used by other colonizers. This way, an interactive biofilm network can develop, leading to high richness and diversity. In agreement with this, in our study bacterial richness and diversity was higher for the Warnow and especially the nutrient-rich WWTP compared to the Baltic stations.

MP is ubiquitously present in the marine ecosystem with particularly high concentrations in the subtropical ocean gyres (Law et al., 2010; Eriksen et al., 2013; Cózar et al., 2014). Our findings suggest that the significance of MP in these accumulation zones could be even higher than believed to date. Besides plastic, these waters harbor a rich microbial diversity, such as the oligotrophic Sargasso Sea (Venter et al., 2004). Leaving aside the issue of attached potential pathogenic or invasive microorganisms, merely the formation of non-natural, plastic-specific biofilms in such a plastic accumulation zone might have implications for the sensitive ecosystem. The microbial communities in the oceans form complex interactions with each other and with other organisms and play key roles in vital ecosystem processes such as primary production and nutrient cycling (Falkowski et al., 2008). Potentially, these natural synergies are being altered by MP biofilms under certain environmental conditions, including low nutrient availability.

Hyphomonadaceae and Erythrobacteraceae Are Abundant Specifically in MP Biofilms

On the family level, PE and PS colonizing communities consisted of many typical aquatic colonizers, such as members of the family Rhodobacteraceae (Dang et al., 2008). Several families abundant on the plastic after exposure, have previously been found associated with various hard substrata, plastic surfaces just being one of them (Hoellein et al., 2014; Oberbeckmann et al., 2016). Our in situ experiment, however, identified several taxa attracted specifically by plastic surfaces, especially at the coastal Baltic stations. Two families were abundant exclusively on MP across several stations, namely Hyphomonadaceae (mainly Hyphomonas) and Erythrobacteraceae (mainly Erythrobacter). Zettler et al. (2013) found OTUs affiliated with these families to be abundantly associated with, some of them even unique to, plastic fragments floating in the Sargasso Sea. Plastic fragments in the North Pacific Gyre were abundantly colonized by members of the family Hyphomonadaceae and to a slightly lesser extent also Erythrobacteraceae (Bryant et al., 2016). Members of both families were also part of the “core microbiome,” defined by De Tender et al. (2017) for biofilms on PE sheets and ropes incubated on the seafloor in the Belgian part of the North Sea. This conformity with our findings indicates that the ecological conclusions from our Baltic Sea study apply to other ecosystems as well.

Hyphomonadaceae are prosthecate bacteria and produce the polysaccharide holdfast, which enables them to attach firmly to surfaces as primary colonizers. The prosthecae allow them to take up nutrients in a wide surrounding area (Abraham and Rohde, 2014). The outer membrane of Hyphomonadaceae features transport systems that actively take up compounds too large for passive diffusion. The efficient use of nutrients and the ability to attach firmly to even smooth surfaces makes members of the Hyphomonadaceae prime candidates to colonize MP.

Erythrobacteraceae on the other hand, stand out due to their production of carotenoids and bacteriochlorophyll a (Tonon et al., 2014), which is beneficial for a life on surface-floating plastic debris. Interestingly, this family has a high potential for biotechnological purposes, due to the possession of important hydrolases. There is very strong evidence for the ability of members of the genus Erythrobacter to degrade PAH (Gao et al., 2015; Zhuang et al., 2015). Members of Erythrobacter might be able to degrade PAH associated with plastic. MP sorb PAH from the surrounding water, and especially PE has a high capacity to accumulate PAH (Adams et al., 2007; Teuten et al., 2009; Rochman et al., 2013). Therefore, the metabolization of MP-associated PAH by members of Erythrobacter is a likely explanation for the significantly higher abundance of Erythrobacter on PS and particularly PE compared to natural surfaces in our study.

The genus Tenacibaculum (Flavobacteriaceae) has consistently been reported to be abundant in plastic colonizing communities, across polymers, and sampling areas (Zettler et al., 2013; Bryant et al., 2016; Oberbeckmann et al., 2016; De Tender et al., 2017). While not being among the most abundant taxa on MP from our experiment, Tenacibaculum displayed significantly higher abundance on PE and PS than in other sample types, and should be investigated further regarding its role in MP biofilms.

Colonization of PE and PS by Potentially Pathogenic Microorganisms

Potentially pathogenic Vibrio spp. have been recently in the spotlight, following a study detecting their high relative abundances (>20%) on a PP particle from the North Atlantic (Zettler et al., 2013). Subsequently, members of the genus Vibrio were also found on plastic debris at a Scottish beach (Quilliam et al., 2014) and on floating MP in the North and Baltic Sea (Kirstein et al., 2016). Bryant et al. (2016), however, did not confirm high abundances of Vibrionaceae on plastic debris from the North Pacific Gyre. Our study revealed that Vibrio bacteria did colonize PE and PS in the Warnow and coastal Baltic Sea, but in lower concentrations than they colonized the surface of wood and other natural particles in the water column. The enrichment of Vibrio sp. in aquatic environments is common on wooden debris, plants and other natural surfaces (Takemura et al., 2014), but might play a minor role with regard to the colonization of MP.

It has also been previously hypothesized that MP in WWTPs might take up sewage-related microorganisms and import them into the environment via the treated wastewater effluent (McCormick et al., 2014; Oberbeckmann et al., 2015). Our station 4, located in the Warnow estuary close to a WWTP runoff, featured very high nutrient levels, indicating a strong impact by the treated wastewater. Indeed, elevated levels of the family Enterobacteriaceae were detected on PE and PS at station 4 as well as within the WWTP (station 7). This family comprises several potentially gut-related and pathogenic genera (Octavia and Lan, 2014). The relative abundances of Enterobacteriaceae on MP, however, remained low (<0.5%) in relation to the levels in the water communities outside and within the WWTP. Despite the low numbers of potential pathogens associated with MP, the possible long-term effect of MP colonization should be considered. Plastic fragments represent more stable surfaces than wood or most other natural particles, which might lead to more durable biofilms. Already Hoellein et al. (2014) presumed that plastic biofilms might stay intact longer and transport biofilms further compared to natural surfaces.

One genus significantly more abundant on plastic than in other sample types within the WWTP was Sphingopyxis (Sphingomonadaceae). Sphingomonadaceae are known to be a reservoir for antibiotic-resistance (Vaz-Moreira et al., 2011; Iredell et al., 2016). More specifically, the genus Sphingopyxis has been described as secondary colonizer of PE drinking water pipes (Douterelo et al., 2014), and a draft genome analysis revealed the presence of a class I integron associated with two antibiotic resistance genes (Gomez-Alvarez et al., 2016). Class I integrons are often embedded in plasmids or transposons, enabling their lateral transfer (Gillings et al., 2008), and have been suggested to be a proxy for anthropogenic pollution (Gillings et al., 2015). The potential for MP to be a hotspot for the transport and transfer of antibiotic resistance has not been investigated so far, but there is urgent need to do so.

Conclusions

Our study highlights the role of MP for specific microbial communities as subject to environmental conditions, including nutrient levels. Consequently, this emphasizes the relevance of MP accumulation zones, such as the oligotrophic Sargasso Sea. Most potentially pathogenic taxa were not enriched on MP; however, the presence of certain bacteria commonly associated with antibiotic resistance on MP in a WWTP suggests that small plastic particles may be hotspots for horizontal gene transfer. In comparison with natural surfaces, MP attracted significantly higher numbers of members of the families Hyphomonadaceae and Erythrobacteraceae, with the latter potentially being associated with PAH degradation. Future studies should examine, whether our short-term observations can be extended to longer time frames and to other important marine areas.

Author Contributions

SO and ML designed the experiment and analyzed the data. SO and BK performed the experiments and measurements. All authors were involved in data interpretation, drafting of the manuscript, and read and approved the final manuscript.

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

The authors thank Jana Normann (University Rostock), Christian Burmeister, Heike Benterbusch, Lars Möller, and Stephanie Mothes (all IOW) for technical assistance. The authors thank the workshop team of IOW, Alexander Hentzsch (IOW), Jennifer Rudolph (DSMZ), and Therese Kettner (IGB) for sampling and setup support, and Gunnar Gerdts (AWI) for statistical advice. The authors acknowledge access to sampling stations by Bootsterassen Warnemünde, Marine Science Center Hohe Düne (Prof. Dr. G. Dehnhardt, Dr. S. Wieskotten) and the WWTP. The work was funded by the Leibniz Association in the scope of the project MikrOMIK (SAW-2014-IOW-2) to ML. Purchase of the Illumina MiSeq was kindly supported by the EU-EFRE (European Funds for Regional Development) program and funds from the University Medicine Rostock. We thank the editor and two reviewers for helpful comments on this manuscript.

Supplementary Material

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

References

Abraham, W. R., and Rohde, M. (2014). “The family Hyphomonadaceae,” in The Prokaryotes: Alphaproteobacteria and Betaproteobacteria, eds E. Rosenberg, E. F. Delong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin; Heidelberg: Springer), 283–299.

Google Scholar

Adams, R. G., Lohmann, R., Fernandez, L. A., Macfarlane, J. K., and Gschwend, P. M. (2007). Polyethylene devices: passive samplers for measuring dissolved hydrophobic organic compounds in aquatic environments. Environ. Sci. Technol. 41, 1317–1323. doi: 10.1021/es0621593

PubMed Abstract | CrossRef Full Text | Google Scholar

Amaral-Zettler, L. A., Zettler, E. R., Slikas, B., Boyd, G. D., Melvin, D. W., Morrall, C. E., et al. (2015). The biogeography of the plastisphere: implications for policy. Front. Ecol. Environ. 13, 541–546. doi: 10.1890/150017

CrossRef Full Text | Google Scholar

Anderson, M., Gorley, R., and Clarke, K. (2008). “PERMANOVA+ for PRIMER: Guide to software and statistical methods,” in PRIMER-E (Plymouth: PRIMER-E).

Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral. Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x

CrossRef Full Text | Google Scholar

Anderson, M. J., and Walsh, D. C. I. (2013). PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: what null hypothesis are you testing? Ecol. Monogr. 83, 557–574. doi: 10.1890/12-2010.1

CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300.

Google Scholar

Bryant, J. A., Clemente, T. M., Viviani, D. A., Fong, A. A., Thomas, K. A., Kemp, P., et al. (2016). Diversity and activity of communities inhabiting plastic debris in the North Pacific Gyre. mSystems 1:e00024-16. doi: 10.1128/mSystems.00024-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U.S.A. 108, 4516–4522. doi: 10.1073/pnas.1000080107

PubMed Abstract | CrossRef Full Text | Google Scholar

Carson, H. S., Nerheim, M. S., Carroll, K. A., and Eriksen, M. (2013). The plastic-associated microorganisms of the North Pacific Gyre. Mar. Pollut. Bull. 75, 126–132. doi: 10.1016/j.marpolbul.2013.07.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Chao, A. (1984). Non-parametric estimation of the number of classes in a population. Scand. J. Stat. 11, 265–270.

Google Scholar

Clarke, K. R., and Gorley. R. N. (2006). PRIMER v6: User Manual/Tutorial (Plymouth Routines in Multivariate Ecological Research). Plymouth: PRIMER-E

Cózar, A., Echevarría, F., González-Gordillo, J. I., Irigoien, X., Ubeda, B., Hernández-León, S., et al. (2014). Plastic debris in the open ocean. Proc. Natl. Acad. Sci. U.S.A. 111, 10239–10244. doi: 10.1073/pnas.1314705111

PubMed Abstract | CrossRef Full Text | Google Scholar

Dang, H., Li, T., Chen, M., and Huang, G. (2008). Cross-Ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters. Appl. Environ. Microbiol. 74, 52–60. doi: 10.1128/AEM.01400-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Datta, M. S., Sliwerska, E., Gore, J., Polz, M. F., and Cordero, O. X. (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 Tender, C. A., Devriese, L. I., Haegeman, A., Maes, S., Ruttink, T., and Dawyndt, P. (2015). Bacterial community profiling of plastic litter in the Belgian part of the North Sea. Environ. Sci. Technol. 49, 9629–9638. doi: 10.1021/acs.est.5b01093

PubMed Abstract | CrossRef Full Text | Google Scholar

De Tender, C., Devriese, L. I., Haegeman, A., Maes, S., Vangeyte, J., Cattrijsse, A., et al. (2017). Temporal dynamics of bacterial and fungal colonization on plastic debris in the North Sea. Environ. Sci. Technol. 51, 7350–7360. doi: 10.1021/acs.est.7b00697

PubMed Abstract | CrossRef Full Text | Google Scholar

Douterelo, I., Sharpe, R., and Boxall, J. (2014). Bacterial community dynamics during the early stages of biofilm formation in a chlorinated experimental drinking water distribution system: implications for drinking water discolouration. J. Appl. Microbiol. 117, 286–301. doi: 10.1111/jam.12516

PubMed Abstract | CrossRef Full Text | Google Scholar

Eriksen, M., Maximenko, N., Thiel, M., Cummins, A., Lattin, G., Wilson, S., et al. (2013). Plastic pollution in the South Pacific subtropical gyre. Mar. Pollut. Bull. 68, 71–76. doi: 10.1016/j.marpolbul.2012.12.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Falkowski, P. G., Fenchel, T., and Delong, E. F. (2008). The microbial engines that drive Earth's biogeochemical cycles. Science 320, 1034–1039. doi: 10.1126/science.1153213

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, X., Gao, W., Cui, Z., Han, B., Yang, P., Sun, C., et al. (2015). Biodiversity and degradation potential of oil-degrading bacteria isolated from deep-sea sediments of South Mid-Atlantic Ridge. Mar. Pollut. Bull. 97, 373–380. doi: 10.1016/j.marpolbul.2015.05.065

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerbersdorf, S. U., and Wieprecht, S. (2015). Biostabilization of cohesive sediments: revisiting the role of abiotic conditions, physiology and diversity of microbes, polymeric secretion, and biofilm architecture. Geobiology 13, 68–97. doi: 10.1111/gbi.12115

PubMed Abstract | CrossRef Full Text | Google Scholar

Gillings, M., Boucher, Y., Labbate, M., Holmes, A., Krishnan, S., Holley, M., et al. (2008). The evolution of class 1 integrons and the rise of antibiotic resistance. J. Bacteriol. 190, 5095–5100. doi: 10.1128/JB.00152-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Gillings, M. R., Gaze, W. H., Pruden, A., Smalla, K., Tiedje, J. M., and Zhu, Y. G. (2015). Using the class 1 integron-integrase gene as a proxy for anthropogenic pollution. ISME J. 9, 1269–1279. doi: 10.1038/ismej.2014.226

PubMed Abstract | CrossRef Full Text | Google Scholar

Giovannoni, S. J., Delong, E. F., Olsen, G. J., and Pace, N. R. (1988). Phylogenetic group-specific oligodeoxynucleotide probes for identification of single microbial cells. J. Bacteriol. 170, 720–726. doi: 10.1128/jb.170.2.720-726.1988

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomez-Alvarez, V., Pfaller, S., and Revetta, R. P. (2016). Draft genome sequence of two sphingopyxis sp. strains, dominant members of the bacterial community associated with a drinking water distribution system simulator. Genome Announc. 4, e00183–e00116. doi: 10.1128/genomeA.00183-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Gower, J. (1966). Some distance properties of latent root and vector methods used in multivariate analysis. Biometrika 53, 325–338. doi: 10.1093/biomet/53.3-4.325

CrossRef Full Text | Google Scholar

Grasshoff, K., Kremling, K., and Ehrhardt, M. (1999). Methods of Seawater Analysis, 3rd Edn. Weinheim: WILEY-VCH.

Harrison, J. P., Schratzberger, M., Sapp, M., and Osborn, A. M. (2014). Rapid bacterial colonization of low-density polyethylene microplastics in coastal sediment microcosms. BMC Microbiol. 14:232. doi: 10.1186/s12866-014-0232-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoellein, T., Rojas, M., Pink, A., Gasior, J., and Kelly, J. (2014). Anthropogenic litter in urban freshwater ecosystems: distribution and microbial interactions. PLoS ONE 9:e98485. doi: 10.1371/journal.pone.0098485

PubMed Abstract | CrossRef Full Text | Google Scholar

Iredell, J., Brown, J., and Tagg, K. (2016). Antibiotic resistance in Enterobacteriaceae: mechanisms and clinical implications. BMJ 352:h6420 doi: 10.1136/bmj.h6420

PubMed Abstract | CrossRef Full Text | Google Scholar

Kesy, K., Hentzsch, A., Klaeger, F., Oberbeckmann, S., Mothes, S., and Labrenz, M. (2017). Fate and stability of polyamide-associated bacterial assemblages after their passage through the digestive tract of the blue mussel Mytilus edulis. Mar. Pollut. Bull. 125, 132–138. doi: 10.1016/j.marpolbul.2017.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Kesy, K., Oberbeckmann, S., Müller, F., and Labrenz, M. (2016). Polystyrene influences bacterial assemblages in Arenicola marina-populated aquatic environments in vitro. Environ. Pollut. 219, 219–227. doi: 10.1016/j.envpol.2016.10.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Kettner, M. T., Rojas-Jimenez, K., Oberbeckmann, S., Labrenz, M., and Grossart, H.-P. (2017). Microplastics alter composition of fungal communities in aquatic ecosystems. Environ. Microbiol. 19, 4447–4459. doi: 10.1111/1462-2920.1389

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirstein, I. V., Kirmizi, S., Wichels, A., Garin-Fernandez, A., Erler, R., Löder, M., et al. (2016). Dangerous hitchhikers? Evidence for potentially pathogenic Vibrio spp. on microplastic particles. Mar. Environ. Res. 120, 1–8. doi: 10.1016/j.marenvres.2016.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/AEM.01043-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Lane, D. J. (1991). “16S/23S rRNA Sequencing,” in Nucleic Acid Techniques in Bacterial Systematics, eds E. Stackebrandt and M. Goodfellow (New York, NY: John Wiley and Sons), 115–175.

Google Scholar

Law, K. L., Morét-Ferguson, S., Maximenko, N. A., Proskurowski, G., Peacock, E. E., Hafner, J., et al. (2010). Plastic accumulation in the North Atlantic subtropical gyre. Science 329, 1185–1188. doi: 10.1126/science.1192321

PubMed Abstract | CrossRef Full Text | Google Scholar

Letunic, I., and Bork, P. (2016). Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245. doi: 10.1093/nar/gkw290

PubMed Abstract | CrossRef Full Text | Google Scholar

Mani, T., Hauk, A., Walter, U., and Burkhardt-Holm, P. (2015). Microplastics profile along the Rhine River. Sci. Rep. 5:17988. doi: 10.1038/srep17988

PubMed Abstract | CrossRef Full Text | Google Scholar

Maso, M., Garces, E., Pages, F., and Camp, J. (2003). Drifting plastic debris as a potential vector for dispersing Harmful Algal Bloom (HAB) species. Sci. Mar. 67, 107–111. doi: 10.3989/scimar.2003.67n1107

CrossRef Full Text | Google Scholar

McCormick, A., Hoellein, T. J., Mason, S. A., Schluep, J., and Kelly, J. J. (2014). Microplastic is an abundant and distinct microbial habitat in an urban river. Environ. Sci. Technol. 48, 11863–11871. doi: 10.1021/es503610r

PubMed Abstract | CrossRef Full Text | Google Scholar

Oberbeckmann, S., Löder, M. G. J., and Labrenz, M. (2015). Marine microplastic-associated biofilms – a review. Environ. Chem. 12, 551–562. doi: 10.1071/EN15069

CrossRef Full Text | Google Scholar

Oberbeckmann, S., Loeder, M. G., Gerdts, G., and Osborn, A. M. (2014). Spatial and seasonal variation in diversity and structure of microbial biofilms on marine plastics in Northern European waters. FEMS Microbiol. Ecol. 90, 478–492. doi: 10.1111/1574-6941

PubMed Abstract | CrossRef Full Text | Google Scholar

Oberbeckmann, S., Osborn, A. M., and Duhaime, M. B. (2016). Microbes on a bottle: substrate, season and geography influence community composition of microbes colonizing marine plastic debris. PLoS ONE 11:e0159289. doi: 10.1371/journal.pone.0159289

PubMed Abstract | CrossRef Full Text | Google Scholar

Octavia, S., and Lan, R. (2014). “The family enterobacteriaceae,” in The Prokaryotes: Gammaproteobacteria, eds E. Rosenberg, E. F. Delong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin; Heidelberg: Springer), 225–286.

Google Scholar

Oksanen, J. F., Blanchet, G., Friendly, M., Kindt, R., Legendre, P., Minchin, P. R., et al. (2015). vegan: Community Ecology Package. R package version 2.2–1. Available online at: http://CRAN.R-project.org/package=vegan

Pielou, E. C. (1969). An Introduction to Mathematical Ecology. New York, NY: Wiley Interscience; John Wiley & Sons.

Google Scholar

Polz, M. F., and Cordero, O. X. (2016). Bacterial evolution: genomics of metabolic trade-offs. Nat. Microbiol. 1:16181. doi: 10.1038/nmicrobiol.2016.181

PubMed Abstract | CrossRef Full Text | Google Scholar

Quilliam, R. S., Jamieson, J., and Oliver, D. M. (2014). Seaweeds and plastic debris can influence the survival of faecal indicator organisms in beach environments. Mar. Pollut. Bull. 84, 201–207. doi: 10.1016/j.marpolbul.2014.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Reisser, J., Shaw, J., Hallegraeff, G., Proietti, M., Barnes, D. K., Thums, M., et al. (2014). Millimeter-sized marine plastics: a new pelagic habitat for microorganisms and invertebrates. PLoS ONE 9:e100289. doi: 10.1371/journal.pone.0100289

PubMed Abstract | CrossRef Full Text | Google Scholar

Rochman, C. M., Hoh, E., Hentschel, B. T., and Kaye, S. (2013). Long-Term Field measurement of sorption of organic contaminants to five types of plastic pellets: implications for plastic marine debris. Environ. Sci. Technol. 47, 1646–1654. doi: 10.1021/es303700s

PubMed Abstract | CrossRef Full Text | Google Scholar

Roller, B. R., Stoddard, S. F., and Schmidt, T. M. (2016). Exploiting rRNA operon copy number to investigate bacterial reproductive strategies. Nat. Microbiol. 1:6160. doi: 10.1038/nmicrobiol.2016.160

PubMed Abstract | CrossRef Full Text

RStudio Team (2015). RStudio: Integrated Development for R. Boston, MA: RStudio, Inc. Available online at: http://www.rstudio.com/

Rummel, C. D., Jahnke, A., Gorokhova, E., Kühnel, D., and Schmitt-Jansen, M. (2017). Impacts of biofilm formation on the fate and potential effects of microplastic in the aquatic environment. Environ. Sci. Technol. Lett. 4, 258–267. doi: 10.1021/acs.estlett.7b00164

CrossRef Full Text | Google Scholar

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, C. E. (1948). A mathematical theory of communication. Bell. Syst. Tech. J. 27, 379–423. doi: 10.1002/j.1538-7305.1948.tb01338.x

CrossRef Full Text | Google Scholar

Stanley, N. R., and Lazazzera, B. A. (2004). Environmental signals and regulatory pathways that influence biofilm formation. Mol. Microbiol. 52, 917–924. doi: 10.1111/j.1365-2958.2004.04036.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Takemura, A. F., Chien, D. M., and Polz, M. F. (2014). Associations and dynamics of Vibrionaceae in the environment, from the genus to the population level. Front. Microbiol. 5:38. doi: 10.3389/fmicb.2014.00038

PubMed Abstract | CrossRef Full Text | Google Scholar

Teuten, E. L., Saquing, J. M., Knappe, D. R., Barlaz, M. A., Jonsson, S., Björn, A., et al. (2009). Transport and release of chemicals from plastics to the environment and to wildlife. Philos. Trans. R. Soc. B. 364, 2027–2045. doi: 10.1098/rstb.2008.0284

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, J. R., Randa, M. A., Marcelino, L. A., Tomita-Mitchell, A., Lim, E., and Polz, M. F. (2004). Diversity and dynamics of a North Atlantic coastal Vibrio community. Appl. Environ. Microbiol. 70, 4103–4110. doi: 10.1128/AEM.70.7.4103-4110.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Tonon, L. A. C., Moreira, A. P. B., and Thompson, F. (2014). “The Family Erythrobacteraceae,” in The Prokaryotes: Alphaproteobacteria and Betaproteobacteria, eds E. Rosenberg, E. F. Delong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin; Heidelberg: Springer), 213–235.

Google Scholar

Vaz-Moreira, I., Nunes, O. C., and Manaia, C. M. (2011). Diversity and antibiotic resistance patterns of Sphingomonadaceae isolates from drinking water. Appl. Environ. Microbiol. 77, 5697–5706. doi: 10.1128/AEM.00579-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Venter, J. C., Remington, K., Heidelberg, J. F., Halpern, A. L., Rusch, D., Eisen, J. A., et al. (2004). Environmental genome shotgun sequencing of the Sargasso Sea. Science 304, 66–74. doi: 10.1126/science.1093857

PubMed Abstract | CrossRef Full Text | Google Scholar

Zettler, E. R., Mincer, T. J., and Amaral-Zettler, L. A. (2013). Life in the “Plastisphere”: Microbial communities on plastic marine debris. Environ. Sci. Technol. 47, 7137–7146. doi: 10.1021/es401288x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhuang, L., Liu, Y., Wang, L., Wang, W., and Shao, Z. (2015). Erythrobacter atlanticus sp nov., a bacterium from ocean sediment able to degrade polycyclic aromatic hydrocarbons. Int. J. Syst. Evol. Microbiol. 65, 3714–3719. doi: 10.1099/ijsem.0.000481

PubMed Abstract | CrossRef Full Text | Google Scholar

Zobell, C. (1943). The Effect of solid surfaces upon bacterial activity. J. Bacteriol. 46, 39–56.

PubMed Abstract | Google Scholar