- Open Access
A multi-OMIC characterisation of biodegradation and microbial community succession within the PET plastisphere
Microbiome volume 9, Article number: 141 (2021)
Plastics now pollute marine environments across the globe. On entering these environments, plastics are rapidly colonised by a diverse community of microorganisms termed the plastisphere. Members of the plastisphere have a myriad of diverse functions typically found in any biofilm but, additionally, a number of marine plastisphere studies have claimed the presence of plastic-biodegrading organisms, although with little mechanistic verification. Here, we obtained a microbial community from marine plastic debris and analysed the community succession across 6 weeks of incubation with different polyethylene terephthalate (PET) products as the sole carbon source, and further characterised the mechanisms involved in PET degradation by two bacterial isolates from the plastisphere.
We found that all communities differed significantly from the inoculum and were dominated by Gammaproteobacteria, i.e. Alteromonadaceae and Thalassospiraceae at early time points, Alcanivoraceae at later time points and Vibrionaceae throughout. The large number of encoded enzymes involved in PET degradation found in predicted metagenomes and the observation of polymer oxidation by FTIR analyses both suggested PET degradation was occurring. However, we were unable to detect intermediates of PET hydrolysis with metabolomic analyses, which may be attributed to their rapid depletion by the complex community. To further confirm the PET biodegrading potential within the plastisphere of marine plastic debris, we used a combined proteogenomic and metabolomic approach to characterise amorphous PET degradation by two novel marine isolates, Thioclava sp. BHET1 and Bacillus sp. BHET2. The identification of PET hydrolytic intermediates by metabolomics confirmed that both isolates were able to degrade PET. High-throughput proteomics revealed that whilst Thioclava sp. BHET1 used the degradation pathway identified in terrestrial environment counterparts, these were absent in Bacillus sp. BHET2, indicating that either the enzymes used by this bacterium share little homology with those characterised previously, or that this bacterium uses a novel pathway for PET degradation.
Overall, the results of our multi-OMIC characterisation of PET degradation provide a significant step forwards in our understanding of marine plastic degradation by bacterial isolates and communities and evidences the biodegrading potential extant in the plastisphere of marine plastic debris.
Plastics are both ubiquitous and problematic in the marine environment . Since the mass-production of plastic began almost 70 years ago, annual production has reached hundreds of millions of tonnes . There are a plethora of routes for plastic to end up in the ocean, such as mismanagement of waste , lost or discarded fishing gear , fibres released during the washing of clothes  or microplastics in cosmetic products . Current estimates put plastic input into the oceans at 19 to 23 million tonnes every year  and, although their ultimate fate and durability is currently unknown , some suggest a persistence of hundreds of years  or fragmentation rates as low as 1-5% per year . Unfortunately, the properties that make plastics so widely used also underlie the reasons that they are so difficult to break down [11, 12].
Whilst a large number of microbes have been reported to degrade different types of plastics, currently very few studies identify the mechanisms and enzymes involved. One particularly noteworthy exception to this is in the degradation of polyethylene terephthalate (PET), where a number of PET hydrolases, termed PETases, have been identified . The PETase that has garnered the most attention to date is that of the bacterium Ideonella sakaiensis, a terrestrial Betaproteobacterium isolated from outside a plastic bottle factory [14, 15]. This PETase is different from other similar esterases as it exhibits higher hydrolytic activity on PET than other substrates and is active at lower temperatures. The isolation of I. sakaiensis and its PETase allowed for further characterisation [16, 17] and engineering  of the enzyme, as well as for homologues to be searched for in environmental metagenomes , thus gaining a broader understanding of how widespread the ability to degrade PET is. Danso et al.  were also successful in expressing a PETase originating from a marine metagenome in the laboratory, although, to our knowledge, currently no marine microbes with this ability have been isolated, so the conditions necessary for their growth and the metabolic pathways they use to catabolise PET degradation subproducts are not yet known.
The organisms found colonizing plastics in the ocean, termed the ‘plastisphere’ , are clearly distinct from microbial assemblages found in the surrounding water  and can differ from those colonizing natural surfaces . The plastisphere may be specific to a particular polymer type [23,24,25,26], location [22, 26,27,28,29,30,31,32] or season [30, 31], but the largest factors shaping plastisphere communities are (i) the environment they are incubated in; (ii) the methodology used for collection and sequencing and (iii) the incubation time used . Microbial communities that colonise surfaces and polymers—including plastics—in the marine environment are complex and are known to go through distinct stages of community succession [34, 35]. This means that if the time at which a microbial biofilm community is being studied is not right, then a community efficient at degrading that compound may not be identified , which could impact the differences found between those microbial communities attached to plastics and those attached to glass [25, 31, 37] or natural particles [22, 38]. Previous studies characterising succession on marine plastics have also found that plastic-specific communities are only found at earlier stages of colonisation and that these communities tend to converge at later time points as the biofilm matures [12, 34, 39, 40].
In order to determine the fate of plastics in the oceans, it is important not only to characterise the real biodegrading potential extant within the plastisphere but also to understand the microbial community dynamics that may be driving this degradation. To address this current gap in knowledge, we took two approaches: (i) to study the dynamics of a microbial community obtained from marine plastic debris when exposed to PET (and its derivatives) over 6 weeks and (ii) to determine the PET biodegrading potential within the plastisphere, as well as in two fully characterised isolates, Thioclava sp. BHET1 and Bacillus sp. BHET2, both able to break down amorphous PET and its derivatives. The proteogenomic and metabolomic analysis of these marine isolates allowed the identification of potential PETases and catabolic pathways involved in PET degradation in marine ecosystems.
Microbial community succession on PET
The microbial community obtained from beached plastics diverged over time in all treatments when incubated with (i) no additional carbon (control), (ii) amorphous PET films, (iii) PET powder, (iv) weathered PET powder or (v) PET monomer BHET (bis(2-hydroxy ethyl) terephthalate (Fig. 1)). Whilst PET powder had a highly crystalline conformation, amorphous PET and BHET monomers were expected to be more accessible for microbial biodegradation. Interestingly, although the assessment of microbial growth is problematic in such settings , DNA extraction yields significantly increased over time for the weathered and non-weathered PET powder as well as the amorphous PET film treatments when compared with the no-carbon control (Figure S1). Bacterial community structure was assessed via 16S rRNA gene sequencing obtained from the inoculum as well as from all five treatments across 6 weeks of incubation (i.e. days 1, 3, 7, 14, 21, 30, 42), including separate analysis of planktonic and biofilm communities grown with the amorphous PET films. Only two samples (replicate 2 from day 42 amorphous PET biofilm and replicate 1 from day 42 PET powder) as well as procedural controls (i.e. DNA extraction and PCR negative controls) were removed due to low numbers of reads (<1000). All other samples had a minimum of 4000 reads and a mean of 19,000 reads per sample, with a total of 18,114 Amplicon Sequence Variants (ASVs) being detected across all samples.
Microbial community differences between treatments
Gammaproteobacteria dominated all samples (>60% relative abundance; Figure S2), largely due to the dominance of Vibrionaceae at all time points along with Alteromonadaceae (e.g. ASV2, maximum abundance 25%) and Thalassospiraceae (e.g. ASV18, maximum abundance 9%) during early stages and Alcanivoraceae (e.g. ASV8, maximum abundance 15%) at later stages in all treatments but BHET (Fig. 1 and Figure S2). All communities were both significantly different (Fig. 1a and Table S1; ANOSIM R=0.885, p=0.001) and less diverse than (Figure S3) the inoculum, but the amorphous PET biofilm and BHET treatments showed the most remarkable differences from all other communities (ANOSIM R=0.709, p=0.001 and R=0.446; p=0.001, respectively; Fig. 1a and Table S1). This may come as a consequence of a higher availability of substrate when compared with the other treatments. The principal response curve redundancy analysis (PRC)  was used to identify the ASVs that drove these community differences (i.e. had weights higher or lower than one) or contributed to similarities (i.e. had weights close to one) between treatments (Fig. 1b). Interestingly, many of the ASVs that were contributing towards the differences between the amorphous PET biofilm and no carbon control treatments were not present in the BHET treatment or were present only in very low relative abundances and vice versa (Fig. 1c; Table S2). The planktonic communities surrounding the amorphous PET films became more similar to the no carbon control communities over time, possibly due to a reduction in their access to the PET substrate as mature biofilms developed on the material. On the other hand, the PET and weathered PET powder communities slowly diverged from the no carbon control communities over time, suggesting a possible increase in crystalline-PET degradation and substrate availability as also supported by the increase of PETases within the community, as shown below.
Community succession induced by PET-like substrates
The clear distinctness of the microbial communities growing in the presence of BHET and amorphous PET films suggests the presence of an available substrate (i.e. the BHET monomer and available PET chains from low crystallinity/amorphous PET films) that may have selected for distinct biodegrading microbes. To further identify these microbial groups, ASVs were defined as early, middle or late colonisers depending on whether they peaked in abundance on days 1-7, 14-30 or 42, respectively (Fig. 2). This was carried out separately for each treatment and ASVs were only included if they were above 0.5% abundance in at least one time point for that treatment. Overall, this analysis identified 77 ASVs, of which some showed a clear early (n=24), middle (n=15) or late colonisation pattern (n=2), and also revealed a number of ASVs (n=42) that were prevalent in only one condition. Interestingly, the ASVs that drove community divergence in amorphous PET biofilm (PRC analysis contribution to effect >2) or BHET treatments (contribution to effect <0.9; see Fig. 1b) compared with the no carbon control were predominantly early and middle colonisers, although these ASVs were also present in other treatments (Fig. 2; Table S3). Other ASVs that were abundant but that did not contribute to community divergence in the PRC analysis in Fig. 1b (i.e. with weights close to one), were generally middle or late colonisers, or varied between treatments (Fig. 2). This may suggest that the readily available substrates to both the amorphous PET biofilm and BHET microbial communities, that initially exert a selection for organisms that are capable of degrading them, may be depleted after 1 week of incubation. After substrate depletion, the community experiences a succession similar to other treatments in which the substrate is less available. Curiously, the thermal weathering of PET did not produce an apparent increase in PET-derived substrate availability. PET is a highly thermostable polymeric material and, therefore, only small (but significant; two independent samples T test p<0.05) chemical variations were observed by Fourier-transform infrared spectroscopy (FTIR) after 9 months of thermal weathering (incubation at 80 °C; Figure S4). Specifically, there were small increases in the ratios between the reference wavenumber 1410 cm−1 (I1410) and wavenumbers corresponding to C=O (I1711/I1410) and C-O (I1240/I1410) carboxylic acid, C-O (I1090/I1410) ester and C-H (I725/I1410) aromatic bonds (Figure S4). This low thermal oxidation and generation of oligomeric PET by weathering may explain the high similarity observed between weathered and non-weathered PET powder exposed communities (Table S1 and Fig. 1) and colonisation dynamics in these treatments (Fig. 2). This contrasts with other plastic materials, e.g. polyethylene, that release large amounts of carbon when thermo-oxidised, which is known to induce the growth of a distinct microbial community during early stages of colonisation [40, 42].
Isolation of PET-degrading microbes from the marine plastisphere
The isolation of PET-degrading microbes from plastispheres was carried out to further confirm the biodegrading potential extant on marine plastic debris and characterise the metabolic pathways involved. Microbial enrichments using a mix of PET powder and BHET led to the isolation of two bacteria that grew on agar plates with BHET as the sole carbon source. These were identified through partial sequencing of their 16S rRNA genes as Thioclava dalianensis (99% identity) and Bacillus aquimaris (99% identity) (named hereafter as Thioclava sp. BHET1 and Bacillus sp. BHET2) and were selected for further proteogenomic and metabolomic characterisation. Their genome sequences revealed that Thioclava sp. BHET1 had a genome size of 7.66 Mb, with 7568 coding sequences and a GC content of 63.26%, whilst Bacillus sp. BHET2 had a genome size of 4.23 Mb with 4368 coding sequences and a GC content of 40.97% (Table S4).
Distribution of isolates Thioclava sp. BHET1 and Bacillus sp. BHET2 amongst marine plastispheres and the global ocean
We searched for sequences matching both isolates in our PET community incubations as well as in all marine plastisphere community samples included in our recent meta-analysis  and in the global oceanic survey Tara [43,44,45]. As expected, we found sequences matching both isolates within our PET community incubations, being highest in abundance within the inoculum (0.42%) and the amorphous PET biofilm (0.08% on day 7) for Thioclava sp. BHET1 (Fig. 3a) and Bacillus sp. BHET2 (Fig. 3d), respectively. Interestingly, when exploring the abundance of both isolates in available plastisphere studies from around the globe, they were both predominantly found on plastics obtained from the North-East Atlantic Ocean and the Mediterranean Sea (Fig. 3b and e), with Thioclava sp. BHET1 reaching a maximum relative abundance of 4.4% (>99% identity) in a plastic sample of unknown type (Fig. 3c) and Bacillus sp. BHET1 reaching a maximum relative abundance of 6.8% (>99% identity) in a PVC sample (Fig. 3f). In the Tara Oceans dataset, where planktonic marine communities were surveyed around the world, there were no sequences above 99% identity with either isolate, although sequences sharing >97% identity with both isolates were widely distributed (Figure S5). Sequences matching (>97% identity) Thioclava sp. BHET1 and Bacillus sp. BHET2 reached maximum abundances of 2.91% and 0.03%, respectively, in surface waters, and 1.20% and 0.03%, respectively, in deep waters (Figure S5). Hence, whilst both isolates seem to prevail in marine plastispheres, they are likely rare within the planktonic microbiome. Thioclava sp. BHET1 and Bacillus sp. BHET2 are certainly not the only taxa performing PET degradation within the microbial community and as shown further below, where key PET biodegradation intermediates accumulate in the milieu of these cultures, PET metabolisation is likely a task more efficiently performed by a consortium of microbes.
PET degradation by the isolates Thioclava sp. BHET1 and Bacillus sp. BHET2
PET is known to be degraded through an initial hydrolysis (by a PET hydrolase, or PETase) to PET oligomers, BHET, mono(2-hydroxy ethyl) terephthalate (MHET), terephthalic acid and ethylene glycol . MHET may then be acted upon by a MHET hydrolase, producing ethylene glycol and terephthalic acid, although PETases may also exhibit hydrolytic activity on BHET and MHET [15, 17]. Ethylene glycol metabolism usually takes place either via conversion to acetaldehyde and acetate [46, 47] or via the formation of glyoxylate, whilst terephthalic acid is usually metabolised to protocatechuate via dioxygenases . A search of the genomes of both isolates was carried out for potential PETases involved in PET hydrolysis using a Hidden Markov Model (HMM) constructed from known PETase sequences (Table S5), as in Danso et al. , revealing seven enzymes that were above the inclusion threshold in Thioclava sp. BHET1 and four in Bacillus sp. BHET2 (Table S6). Interestingly, though, whilst the genome of Thioclava sp. BHET1 encoded canonical catabolic pathways for PET intermediate degradation (e.g. terephthalic acid degradation), these were not found in Bacillus sp. BHET2 (Table S4).
A comprehensive proteomic (i.e. of the bacterial isolates; Tables S7 and S8) and metabolomic analysis (i.e. of both the bacterial isolates and microbial communities; Tables S9 and S10) was performed to further identify the mechanisms used by marine microbes to breakdown PET and its intermediates.
Proteomic analysis of PET degradation by two marine isolates
The two new marine isolates, Thioclava sp. BHET1 and Bacillus sp. BHET2, were incubated with the labile substrate fructose (Figure S6), amorphous PET films, BHET and terephthalic acid for a full cellular- and extracellular-proteomic analysis that would shed light on the pathways and enzymes induced by PET substrates.
In Thioclava sp. BHET1, the proteomic data suggests the esterase with coding sequence (CDS) 0051 as the enzyme involved in PET depolymerisation, which was slightly increased in the presence of PET and its derivatives (i.e. BHET and terephthalic acid) relative to the fructose positive control (Fig. 4). This enzyme was highlighted as a possible PETase by our HMM analysis (Table S6), is predicted to be secreted and contains an alpha/beta hydrolase fold domain which is characteristic of this kind of esterase. No enzymes that were similar to the MHETase of I. sakaiensis could be identified in the Thioclava sp. BHET1 genome; however, the carboxylesterase (CDS 1741) is likely also capable of this hydrolysis and was 2.1-, 3.5- and 2.3-fold more abundant in the PET, BHET and terephthalic acid cellular proteomes, respectively. There were also several tripartite ATP-independent periplasmic (TRAP) transporters that were upregulated in all treatments when compared with the positive control (Table S11), and that could be involved in the transport of PET degradation products, i.e. ethylene glycol and terephthalic acid, into the cell.
As expected, enzymes involved in ethylene glycol catabolism (encoded by the CDSs 0884, 0999 and 0538; Fig. 4) were more abundantly detected in Thioclava sp. BHET1 in the BHET treatment. Particularly, the acetaldehyde dehydrogenase (encoded by CDS 0999)—necessary for the conversion of acetaldehyde to acetate—was highly abundant in the cellular proteomes, representing 3.33 and 1.42% relative abundance in the BHET and PET treatments, respectively.
The other PET biodegradation product, terephthalic acid, is usually converted to protocatechuate for its degradation. In this case, the Thioclava sp. BHET1 proteins that were annotated with these functions (i.e. proteins tphA1, tphA2, tphA3 and tphB encoded by the gene cluster 2867-2870; Table S4) were not detected in the proteomes (Table S8). The conversion of terephthalic acid to protocatechuate is more likely catalysed by other terephthalate dioxygenase orthologues, i.e. 1142-1143, that were particularly induced by the presence of terephthalic acid (Fig. 4). Protocatechuate is expected to be funnelled into the β-ketoadipate pathway via 3-oxoadipate-enol-lactone although, again, the expected enzyme (i.e. protocatechuate 3,4-dioxygenase made by subunits pcaG and pcaH and encoded by genes 1124-1125) was not detected or not particularly induced by the presence of PET degradation products. Interestingly, though, the incredibly high induction of a catechol 1,2-dioxygenase (i.e. catA encoded by 1817; over 8000× increased in the terephthalic acid treatment versus the control) and a muconate cycloisomerase (i.e. catB encoded by 1816; 57× increased) may suggest that protocatechuate may be degraded via catechol, as previously hypothesised by Hara et al.  (Fig. 4). Nevertheless, the lack of a specific protocatechuate decarboxylase in Thioclava sp. BHET1 raises the question of whether the catA and catB-like enzymes may directly attack protocatechuate. Curiously, PET and PET sub-products BHET and terephthalic acid also seemed to co-induce all enzymes involved in the anaerobic degradation of phenylacetate (i.e. via phenylacetyl-CoA and other intermediates to acetyl-CoA ; encoded by the gene cluster 1572-1593; Table S8) as well as some of the enzymes for aerobic phenylacetate degradation (i.e. via homogentisate and other intermediates to fumarate and acetoacetate ; gene cluster 2777-2782).
The pathway used by Bacillus sp. BHET2 to metabolise PET products could not be determined by our proteogenomic analysis. For Bacillus sp. BHET2, the genome annotations by Prokka and BlastKOALA (Blast, Basic Local Alignment Search Tool; KEGG, Kyoto Encyclopedia of Genes and Genomes; KOALA, KEGG Orthology And Links Annotation)  as well as subsequent local BLAST searches with known terephthalic acid and protocatechuate degradation proteins did not reveal any proteins with significant homology. However, PET treatments did seem to induce a large number of proteins that are usually involved in the degradation of xenobiotics, such as Cytochrome C oxidases and monooxygenases (Table S12) that were upregulated in the PET, BHET and terephthalic acid treatments compared with the control (i.e. with fructose). This, along with the metabolomic analyses shown below (Fig. 4), suggests that these PET compounds are being degraded, but possibly using enzymes that share little homology with those previously described or via a pathway that is currently unknown and needs further characterisation.
PET biodegrading potential in the predicted metagenomes from the communities
Having identified canonical and alternative pathways for PET degradation in Thioclava sp. BHET1, we used PICRUSt2  to determine their predicted abundance in the communities (Fig. 5). The default database used for PICRUSt2 uses the KEGG ortholog annotations for 20,000 genomes contained within the Joint Genome Institute (JGI) Integrated Microbial Genome (IMG) database . The version used did not contain the KEGG ortholog for PETase (i.e. K21104) and, hence, we used the PETase HMM constructed above to determine the abundance of PETases in the JGI genomes (full details are given in the methods section). Using the default HMM E value cut-off (0.01), 416 of the JGI genomes were predicted to contain at least one PETase, which was reduced to 370 genomes by choosing a more stringent cut-off (1 × 10−4; Table S13A). Weighted Nearest Sequenced Taxon Indices (NSTI) ranged between 0.03 (in planktonic amorphous PET samples) and 0.1 (in the inoculum) and had a median of 0.05 (Table S14), indicating that these predictions are expected to be of acceptable accuracy , although we do note that these are not real metagenomes and the results therefore need to be taken with caution.
The predicted metagenome revealed that 156 ASVs (of the 18,114 total ASVs) potentially encoded for a PETase-like enzyme, 130 of which were highest in abundance in the PET or weathered PET powder treatments. We examined the closest JGI genome matches to each of these ASVs (that contained PETases) and found that they were all similar to one of seven genomes: (i) Thalassolituus oleivorans R6-15 (Oceanospirillales; 2 copies; 5 ASVs); (ii) Lentzea violacea DSM 44796 (Actinomycetales; 3 copies; 5 ASVs); (iii) Lentzea flaviverrucosa CGMCC 4.578 (Actinomycetales; 2 copies; 1 ASV, highest in abundance in a DNA extraction control); (iv) Pseudomonas aestusnigri VGX014 (Pseudomonadales; 2 copies; 127 ASVs, 123 of which were highest in abundance in the PET or weathered PET powder treatments); (v) Loktanella atrilutea DSM 29326 (Rhodobacterales; 1 copy; 16 ASVs); (vi) Plantactinospora sp. CNZ320 (Micromonosporales; 1 copy; 1 ASV, highest in abundance in a DNA extraction control); and (vii) Oleibacter sp. HI0075 (Oceanospirillales; 5 copies; 1 ASV, highest in abundance in an amorphous PET biofilm sample) . Each of these PETases was further verified manually through a National Centre for Biotechnology Information (NCBI) conserved domain search ; the top domains found in these predicted PETases were the Abhydrolase super family (11 PETases), PldB super family (3 PETases) or DAP2 super family (2 PETases), typically found in alpha/beta hydrolases, lysophospholipases or dipeptidyl aminopeptidases, respectively (Table S13B) . Of the 156 ASVs that were predicted to contain PETases, only three of these ASVs were above 0.5% abundance at any time point (all most similar to Pseudomonas aestusnigri VGX014): ASV20 (Pseudomonas, 2 copies, maximum abundance 7.4%, NSTI 0.002), ASV43 (Azomonas, 2 copies, maximum abundance 3.6%, NSTI 0.035) and ASV98 (Pseudomonas, 2 copies, maximum abundance 0.9%, NSTI 0.025; taxonomically classified by NCBI BLAST). Each of these three ASVs were middle or late colonisers and, curiously, were only abundant in the PET and weathered PET powder treatments (Figs. 2 and 5). Hence, whilst known PETase-like enzymes were identified in <0.5% of microbial community members in all other treatments, both PET and weathered PET powders showed a remarkable abundance of bacteria that encode one of these enzymes, reaching 20-25% of the microbial community by the end of the incubation (i.e. 1 gene copy per every 4 or 5 bacteria; Fig. 5). Only three confirmed MHETases (i.e. responsible for the conversion of MHET to terephthalic acid) are currently known  and the initial conversion of ethylene glycol to glyoxylate is catalysed by dehydrogenases with broad specificity and, hence, these genes were not included in this analysis.
Enzymes involved in the conversion of terephthalic acid to protocatechuate (i.e. terephthalate 1,2-dioxygenases) were predicted by using a HMM of known terephthalate degradation genes as done above for PETases. The genes tphA2 and tphA3 showed a general decrease in abundance over time in the PET and weathered PET powder treatments, as well as in the no carbon control treatment (Fig. 5). These enzymes are only useful after an initial conversion of the PET, or BHET, to terephthalic acid, and we had therefore expected that in the PET treatments the pattern of their abundance would follow that of the PETases, i.e. would increase in abundance over time. It is possible that the rate at which PET is being hydrolysed is too slow to exert an effect on the abundance of genes for terephthalic acid degradation. It is interesting, though, to note the high abundance of tphB in the biofilm on amorphous PET and in the presence of BHET, possibly because these were the treatments where terephthalic acid was most available (Fig. 5). Interestingly, the alternative pathway detected by proteomics in Thioclava sp. BHET1 for terephthalic acid degradation, i.e. genes 1142-1143 (Fig. 4), followed a similar abundance pattern as tphB (Figure S7) and may well be worth further biochemical characterization to confirm the hypothesised function given in this study.
The abundance of the genes involved in catalysing protocatechuate towards the β-ketoadipate pathway, i.e. the genes pcaG and pcaH, were remarkably abundant in the BHET treatment (Fig. 5), as this substrate may be more readily available than PET. We also explored the abundance within the communities of the alternative pathway suggested for protocatechuate degradation via catechol, i.e. catechol dioxygenase genes catA and catB, because of its strong induction in the proteome of Thioclava sp. BHET1. In this case, the abundance of both catA and catB decreased over time in almost all treatments when compared with their abundance in the community inoculum (Figure S7). We also analysed the abundance of phenylacetate degradation genes, a pathway that seemed to be co-regulated by the presence of PET sub-products in Thioclava sp. BHET1, observing a consistent increase in abundance of all genes paaABCDEGJKZ in the amorphous PET biofilm treatment (Figure S7).
Whilst Bacillus sp. BHET2 did not encode for any of the known enzymes for terephthalate biodegradation, there were a considerable number of oxidases and monooxygenases upregulated in its proteome when exposed to PET substrates. Despite that these are very generic enzymes, we analysed the abundance of mono- and dioxygenases in each one of the communities and found, on average, more than one dioxygenase and monooxygenase gene copy per bacterium in the predicted metagenomes (Figure S7). Hence, no analysis of the distribution of the biodegradation pathway of Bacillus sp. BHET2 could be made within the communities.
Metabolomic assessment of the degradation of PET by isolates and communities
The detection of PET degradation intermediates and the build-up of these metabolites in the culture supernatant are the clearest evidence of PET breakdown and, furthermore, flags bottlenecks where the biodegrading potential of the bacteria may be less efficient (Fig. 4). Non-targeted metabolomics were carried out on the supernatants of our newly isolated marine strains Thioclava sp. BHET1 and Bacillus sp. BHET2 as well as the characterised terrestrial PET-degrader I. sakaiensis [15, 18, 58], when incubated with amorphous PET films, BHET, terephthalic acid and fructose. Substrates with no inoculum were included as negative controls. We also performed a non-targeted metabolomic analysis on each of the community incubations on day 42 in order to identify products of PET degradation.
The metabolomic analyses confirmed that all three bacterial isolates as well as the BHET-grown microbial community (i.e. the community incubated with BHET for 42 days) were able to hydrolyse BHET as they significantly accumulated MHET (Thioclava sp. BHET1 and the BHET community), terephthalic acid (Bacillus sp. BHET2) or both (I. sakaiensis), when compared with control incubations with no bacteria (Fig. 4 and Tables S9 and S10). Furthermore, two potential oxidised derivatives of BHET (C12H12O6 and C12H12O7) also accumulated significantly in the Thioclava sp. BHET1 and I. sakaiensis incubations, i.e. C12H12O6 and C12H12O7, the first of which was also identified in the incubations with the BHET-grown microbial communities (Fig. 4, Tables S9 and S10). The generation of these oxidised derivatives of BHET seems to occur only when this substrate is in excess as they were not detected in the PET treatments, and could be carried out by oxidases, oxidoreductases or dehydrogenases, of which there are many detected in the Thioclava sp. BHET1 proteome. Most interestingly, during the incubations with PET, BHET accumulated in the supernatants of Thioclava sp. BHET1 and Bacillus sp. BHET2 (i.e. almost 22-fold and 2.6-fold higher than in controls, respectively; Fig. 4 and Table S9). In incubations with I. sakaiensis, BHET also accumulated in the PET treatments (1.88-fold higher than in controls) although this accumulation was not statistically significant.
Curiously, no PET degradation sub-products were observed in the community incubations where polymeric PET was present. PET sub-products may have not been observed because the rate at which PET hydrolysis occurs is lower than the assimilation of these oligomers by the microbial community and, therefore, there is no measurable build-up of these metabolic intermediates. This may also explain the strong differences between the biofilm and planktonic microbial community in the amorphous PET film condition (Figs. 1 and 2), where degradation intermediates may be rapidly consumed close to the surface of the plastics and are not accessible to the planktonic community. Hence, due to the lack of PET degradation intermediates, it is not surprising that only the supernatants of BHET treatments grouped separately from supernatants from all other conditions on the metabolomic nMDS plot (Figure S8). The high accumulation of the degradation intermediate BHET in cultures of Thioclava sp. BHET1 exposed to PET suggests that, whilst it is capable of PET hydrolysis, it is not as efficient at using the BHET as the community is. No degradation intermediates were detected when each of the three bacteria were incubated with terephthalic acid suggesting (i) that the toxicity of terephthalic acid limits degradation of this compound when no other carbon source is present, as we previously found for phthalic acid ; or (ii) the degradation pathway of terephthalic acid has no bottleneck that produces an accumulation of detectable levels of the intermediate in the culture supernatant.
PET surface oxidation
Given the distinct proteomic response in the presence of PET and the metabolomic detection of PET hydrolysis products (Fig. 4), we performed an additional experiment to determine the modifications to the amorphous PET surface after incubation with both our marine isolates (i.e. Thioclava sp. BHET1 and Bacillus sp. BHET2), the microbial community (i.e. the community used to inoculate the PET community succession experiment) and control incubations with no microbial inoculum. In this additional experiment, we incubated the amorphous PET for 5 months and then measured the absorbance at key oxidation peaks using Fourier-transform infrared spectroscopy (FTIR). The increase in absorbance measured by FTIR at these key oxidation peaks is indicative of PET polymer chain hydrolysis, thus exposing more functional groups and specifically leading to (i) an increase in the number of C=O and C-O carboxylic acid end groups (increased I1711 and I1240, respectively); (ii) C-H bending of the aromatic ring (increased I725) and (iii) an increase in the number of C-O ester end groups (increased I1090) [60, 61]. Indexes of peak variation were normalised using the invariable peak I1410 that corresponds to ring C-H bending and ring C-C stretching, as in Donelli et al. . Both isolates significantly oxidised the amorphous PET surface (Fig. 6), but the truly remarkable increase in the oxidation produced by Bacillus sp. BHET2 would be in accordance with a non-specific oxidation carried out by the large number of cytochrome C oxidases and oxygenases detected by proteomics. This may also explain that, whilst Thioclava sp. BHET1 generated a large accumulation of the degradation intermediate BHET, Bacillus sp. BHET2 may produce a diversity of oligomeric intermediates other than BHET. The incubation with the community produced a slight increase in PET surface oxidation although this was not statistically significant (Fig. 6).
The microbial community succession across 6 weeks of incubation with different types of PET substrates as well as a multi-OMIC analysis of two new marine isolates, Thioclava sp. BHET1 and Bacillus sp. BHET2, has provided a comprehensive overview of the plastic biodegrading potential extant in marine plastispheres. We were able to confirm via metabolomic and proteogenomic analyses that Thioclava sp. BHET1 degrades PET through an initial hydrolysis into monomers, mainly BHET. Although a PETase-like candidate was detected by comparative genomics and proteomics in Thioclava sp. BHET1, it is worth noting that this enzyme shares much lower homology with previous PETases identified by Danso et al.  and, hence, further biochemical testing is required to confirm its PET hydrolytic activity. The ester branches of BHET (and MHET) are removed generating terephthalic acid and ethylene glycol, both of which are transported inside the cell and catabolised as indicated in Fig. 4. The proteomic analysis of Thioclava sp. BHET1 revealed some alternative enzymes induced by the presence of PET products, and that may be involved in their biodegradation, which were not anticipated by the genomic screening. The induction of these enzymes (Fig. 4), as well as the co-induction of the phenylacetate catabolic pathway, requires further work to confirm whether they catalyse these reactions or whether these are just co-induced by a common intermediate. Oxygenases are usually involved in the degradation of most aromatic hydrocarbons, such as terephthalic acid and protocatechuate as well as a potential further intermediate, catechol (Fig. 4). These oxygenases tend to be abundant in the marine environment  as well as on marine plastics [21, 63, 64] and, as expected, we found they were abundant in the PICRUSt2 predicted metagenomes from the communities (Figure S7). This suggests that whilst enzymes capable of the hydrolysis of PET, PETases, are generally present in relatively low abundances , the enzymes involved in the metabolisation of PET degradation sub-products may be more widely distributed in the marine environment. Whilst we tend to assign enzymes to very specific substrates, it is not new that some oxygenases may display a broader substrate specificity , especially those with larger catalytic pockets . The enzymes identified in this study are, hence, excellent candidates for further substrate specificity evaluation. The catabolic pathways for terephthalic acid degradation by Bacillus sp. BHET2 were not determined, although the FTIR and metabolomics results suggest that the amorphous PET and BHET were being degraded (Figs. 4 and 6). The large number of cytochrome C oxidases and undefined monooxygenases detected in the proteome of this strain suggests that this bacterium uses a yet uncharacterised pathway for PET degradation.
Both Bacillus and Thioclava spp. have previously been found to be colonisers of plastics in the marine environment [26, 67,68,69,70] and we confirmed that both were present in marine plastisphere biofilms (Fig. 3). Thioclava spp. have also been found to be potent degraders of crude oil , whilst the degradative ability of Bacillus spp. has previously been reported in the terrestrial environment for PET  and also for PE in the marine environment . In Sudhakar et al. , growth was higher on PE that was thermally pre-treated, and it was also suggested that an initial abiotic oxidation step was necessary for Bacillus spp. to be capable of PE degradation. As we recently discussed in detail , abiotic degradation is likely a prerequisite for the biodegradation of many plastics, particularly for plastics that do not contain heteroatoms in their backbone, such as PE and PP . The products of this abiotic degradation for PE, PP, PS and PET are low-molecular weight compounds with oxidised end groups  that share structural similarity with BHET, MHET and terephthalic acid. Curiously, although there was an increase in absorbance ratios relating to an increase in carboxylic acid end groups in FTIR spectra after thermal pre-treatment of PET powder (weathered PET powder; Figure S4), this did not lead to large differences between the microbial communities growing on the weathered and non-weathered PET powder (Figs. 1, 2 and S7). Furthermore, the crystallinity of the PET determines the ease with which degradation may take place , meaning that the high crystallinity PET powder (>40%) used here is more difficult to degrade due to the lack of chain accessibility for microbial attack. PET is also highly thermo-stable , meaning that the effect of the thermal treatment that we use for artificial weathering may be less effective than for other materials, such as we previously observed for polyethylene . Most interestingly, Bacillus sp. BHET2 produced a strong oxidation of the amorphous PET surface (Fig. 6) and, hence, gives evidence that microbial production of extracellular reactive oxygen species may be a potent initiator of recalcitrant polymer degradation.
Typically, on particle surfaces in the marine environment, microbial community succession leads to an initial dominance of particle substrate degraders before these are overtaken by cheaters and cross-feeders, which stands true for both natural particles and polymers [35, 36] as well as marine plastics . Typically, this succession occurs at the taxonomic level of order, but in our 6-week succession experiment, the Gammaproteobacteria make up 70-80% of the community at all time points and succession was only observed at lower taxonomic levels (Figs. 1 and 2). The hydrocarbonoclastic Alteromonadaceae and Thalassospiraceae were abundant at early stages whilst the Alcanivoraceae were abundant at later incubation stages and the Vibrionaceae dominated throughout. Interestingly, one ASV (ASV1), with 99% 16S rRNA gene identity with both Vibrio parahaemolyticus and Vibrio alginolyticus, made up approximately 20-40% of all communities. The inoculum used here was removed from plastics collected at a beach in Porthcawl, Wales, which is close to a wastewater treatment plant and this is one possible explanation for the dominance of Vibrio spp., both in the inoculum and in all samples at all time points (Figs. 1, 2 and S2), emphasising the persistence of such potential pathogenic microbes once they colonise plastic surfaces. Several previous studies have found Vibrio spp. in high abundances on some marine plastics [20, 38] and others have investigated the presence of Vibrio spp. on marine plastics as potentially pathogenic microbial hitchhikers [23, 77,78,79]. Nevertheless, as we [12, 33] and others  have previously highlighted, Vibrio spp. are not exclusively found on plastics; they are well-known marine biofilm-forming microbes  and a consortia of Vibrio spp. was previously reported to degrade a polyvinyl alcohol-linear low-density polyethylene plastic blend . We could not determine whether Vibrio spp. were capable of PET degradation, mainly because no Vibrio-like organism was found amongst our PET biodegrading isolates. Furthermore, the PICRUSt2 predicted metagenome showed that neither ASV1 nor any other Vibrio spp. ASV that contributed more than 0.5% relative abundance possessed any of the genes involved in PET, terephthalic acid or protocatechuate degradation that are shown in Fig. 5, suggesting that the abundant Vibrio spp. are more likely cross-feeders, possibly utilising the organic matter generated by other members of the microbial community [82, 83] rather than PET or any of its immediate breakdown products.
Several previous studies have noted that the effects of substrate type are only likely to be seen at early stages of biofilm formation [35, 84, 85]. The same seemed to be true here, where those species with large contributions to differences between treatments were predominantly early colonisers of their corresponding treatments (defined as those that peaked in abundance between days 1 and 7; Figs. 1 and 2). We also found that genes involved in protocatechuate degradation peaked in abundance on days 3 and 7 for the BHET and days 21 and 30 for the amorphous PET biofilm treatments, fitting with the hypothesis that BHET, being a much more labile substrate than PET, would induce an earlier selection of microbes capable of biodegrading the aromatic ring. Intriguingly, though, potential PETase homologues were only abundantly identified in the PET and weathered PET powder treatments, and these increased over time, reaching their maximum abundance on day 42 of incubation (Fig. 5). This suggests that for highly crystalline PET, degradation by the microbial community had not yet reached its peak within 6 weeks of incubation, which is potentially the reason that the PET and weathered PET powder communities did not differ from the no carbon control communities until the end of the incubations (Fig. 1). We also did not detect intermediates of PET degradation in these communities through our metabolomic analyses, although this could be due to the reduced rates of intermediate production from such a recalcitrant polymer and, thus, these are consumed more rapidly than they are produced. More research—perhaps using isotopically labelled plastics [12, 86]—should therefore be carried out to conclusively determine the biodegradation of high crystallinity PET and track its transfer across the complex microbial community.
Here, we have characterised the first marine PET-degrading bacterial isolates, Thioclava sp. BHET1 and Bacillus sp. BHET2, through a proteogenomic and metabolomic approach. Whilst Thioclava sp. BHET1 revealed an interesting and unanticipated array of enzymes to process PET sub-products, Bacillus sp. BHET2 showed an even more intriguing and uncharacterised range of PET-degrading enzymes. The further characterisation of these enzymes will allow for the search of these enzymes in environmental metagenomes and help assess the real biodegrading potential extant in the plastisphere. We have also characterised the microbial community succession of a plastisphere exposed to different PET products across 6 weeks of incubation and provided some evidence that the peak of PET degradation may occur at different time points depending on the recalcitrance of the substrate (i.e. crystallinity) and the accessibility of the substrate to the microbes. The results of our multi-OMIC analyses highlight the potential for PET degradation that exists in the marine plastisphere and flag the need for a more mechanistic characterisation of the biodegrading potential in different environments.
Materials and methods
Three types of polyethylene terephthalate (PET) were used for microbial incubations: (i) PET powder (crystallinity >40%; particles <300 μm; ES306031 Goodfellow, UK); (ii) amorphous PET films (250 μm thickness; ES301445 Goodfellow, UK); and (iii) weathered PET powder (artificially weathered through incubation at 80 °C for 9 months; Figure S4). The manufacturing intermediates/PET breakdown products bis(2-hydroxy ethyl) terephthalate (BHET), mono(2-hydroxy ethyl) terephthalate (MHET) and terephthalic acid (Sigma Aldrich, UK) were also used.
All microbial cultures were performed in Bushnell-Haas mineral medium  with 3% NaCl (hereafter referred to as mineral medium), unless otherwise stated, and all incubations were carried out at 30 °C in the dark with constant shaking at 150 rpm (liquid cultures only).
The microbial community used as an inoculum for laboratory incubations was obtained from bulk marine plastic debris collected from Porthcawl beach (Wales, UK) in July 2018. The debris were washed with sterile seawater (autoclaved 20 min at 121 °C) and sonicated to detach the biofilm from the plastic. Briefly, the plastics were placed into 50 mL falcon tubes containing 50 mL mineral medium and subject to sonication in a Branson 2510 Ultrasonic water bath for 10 min followed by 30 s vortexing. After removing the plastics, the detached biofilm was concentrated by centrifugation (4000×g for 5 min). The cell pellet was resuspended in 15 mL and 500 μL were used as the inoculum for each of five treatments, all in independent biological replicates (n=15). The five treatments were (i) no carbon source (control); (ii) 10 films of amorphous PET (1.5 × 0.5 cm); (iii) 1% (w/v) PET powder; (iv) 1% (w/v) weathered PET powder and (v) 1% (w/v) BHET. Inoculated (n=15) and non-inoculated control flasks (n=15) were grown in 75 cm2 tissue culture flasks containing 50 mL mineral medium supplemented with one of the four substrates or no additional carbon source (control). Aliquots of 1.5 mL were collected from each flask on days 1, 3, 7, 14, 21, 30 and 42 of incubation and cells were pelleted by centrifugation at 14,000×g for 5 min. Additionally, one film of amorphous PET was also taken at each time point for community analysis. Cell pellets and films were stored in Buffer AL (Qiagen, UK) at −20 °C until further DNA extraction and community analysis.
DNA extraction and amplicon sequencing
Comparison of DNA extraction method efficiency
In order to determine the optimal method for DNA extraction from biofilms, four DNA extraction methods were tested: (i) DNeasy Plant Mini Kit (Qiagen) with the addition of an initial bead-beating step as detailed in ; (ii) DNeasy Blood and Tissue Kit (Qiagen) with modifications as in (a); (iii) DNeasy Power Biofilm Kit (Qiagen), following the manufacturer’s instructions and (iv) GeneJET Genomic DNA Purification Kit (Thermofisher Scientific) following the manufacturer’s instructions for gram-positive bacteria with an additional bead beating step, as in (i). For procedural purposes, all four kits were tested using 60-day colonised amorphous PET films from a preliminary experiment using the same setup as above (no other measurements were taken from these preliminary PET-incubated communities) and mineral medium as procedural blanks. A Qubit® HS DNA kit (Life Technologies Corporation) was used for DNA quantification. The GeneJet Genomic DNA Purification Kit was found to yield the highest concentrations of DNA (Table S15) and was therefore used for all DNA extractions of microbial communities in this study. A Qubit® HS DNA kit was used for DNA quantification, after which samples were diluted to equalise the concentrations across all samples.
The V4-V5 region of the 16S rRNA gene was amplified using primers 515F-Y and 926R  after which PCR indexing, amplicon purification and normalisation was carried out as in . Pooled libraries were additionally quantified using the NEBNext Library Quant Kit for Illumina (New England Biolabs, UK) and diluted to 4 nM. Libraries were denatured using 0.2N NaOH and MiSeq amplicon sequencing was carried out using the MiSeq Reagent Kit v3 (600 cycles; Illumina, UK), following the manufacturer’s instructions for a 14 pM library with 2% phiX as an internal reference. Reads were demultiplexed using Illumina BaseSpace.
Microbial community structure determination and statistical analysis
Sequencing data were processed following the DADA2 (version 1.8.0) pipeline (Callahan et al., 2016b ) in R (version 3.6.1), as in Wright et al. . PICRUSt2 artificial metagenome predictions  were carried out in QIIME2 using the additional packages castor , HMMER , EPA-NG  and gappa . In order to include known genes for PET degradation that are not included in the default PICRUSt2 predictions, all prokaryotic genomes were downloaded from the Joint Genome Institute (JGI) Integrated Microbial Genome (IMG) database  and were filtered using a custom Python script to create a database containing only those genomes used in PICRUSt2. This included 14,286 of the 20,000 PICRUSt2 genomes (these can be found in the Figshare file at ); the others were not available from the JGI genome database. Following the methods of Danso et al. , an alignment of known PET hydrolases (Table S5) was constructed using the Clustal Omega programme  on the European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI) multiple sequence alignment server . This alignment was used for the construction of a Hidden Markov Model (HMM ;) which was used to search the 14,286 genomes for PETase homologues. This was repeated for other genes involved in PET degradation, such as terephthalate dioxygenases (Table S5) and the counts of these genes within the genomes were added to the default PICRUSt2 database. Whilst HMM searches of the JGI genomes used the default E value cut-off of 0.01, we further verified the presence of PETases within our predicted metagenomes by (i) computing a distance matrix  from the tree output by PICRUSt2 to find the PETase-containing organism most closely related (by 16S rRNA sequence) to each of the ASVs predicted to contain a PETase; and (ii) obtaining the predicted PETase sequences from these organisms and performing manual NCBI conserved domain searches . Principal response curves  were calculated using the vegan package in R  and all other analyses were carried out using custom Python (version 3.7.6) scripts (https://github.com/R-Wright-1/PET-plastisphere).
Enrichment, isolation and characterisation of bacteria capable of PET degradation
Tissue culture flasks (25 cm2) containing 25 mL mineral medium supplemented with 0.005% (w/v) yeast extract (Merck KGaA, Germany), 0.1% (w/v) BHET and PET powder, were inoculated with 1 mL of a microbial community obtained from bulk marine plastic debris as indicated above. Cultures were incubated at 30 °C in the dark with shaking at 150 rpm. When growth was visible after 4 weeks (assessed as a change in turbidity), 1 mL was used to repeat this enrichment step, incubating for a further 2 weeks until turbidity was visible.
Isolation and genome sequencing
The microbial community from the second enrichment (100 μL) was spread on replicate solid agar plates made with supplemented Bushnell-Haas mineral medium, as above (i.e. 1% (w/v) BHET and 0.005% (w/v) yeast extract), containing 1.5% (w/v) agar and incubated for 3 weeks at 30 °C in the dark. Morphologically distinct colonies were picked and streaked onto fresh plates until isolates were obtained. The identification of isolates was carried out by partial sequencing of the 16S rRNA gene (GATC BioTech, Germany) using primers 27F and 1492R  after DNA extraction using the DNeasy Plant Mini Kit (Qiagen) with modifications (as above) and purification using the QIAquick PCR purification kit (Qiagen). Two fast growing isolates were sent to MicrobesNG (Birmingham, UK) for whole genome sequencing and were used for further characterisation of their ability to degrade PET. Assembled genomes were annotated using Prokka  and Blast KEGG Orthology and Links Annotation (BlastKOALA)  (Table S4). A PETase Hidden Markov Model (HMM) was constructed, as above, and used to search for PETase homologues in the genomes.
Growth of both isolates was tested in marine broth 2216 (BD Difco™) or mineral medium supplemented with 0.1% (w/v) glucose, fructose, succinate, glycerol, pyruvate or N-acetyl-D-glucosamine in order to define a suitable labile substrate for the control condition. Their growth on each substrate was measured over 3 days by absorbance (600 nm) measurements taken every 30 min on a Synergy HTX microplate reader. Material for cellular and exo-proteomics was generated by growing each isolate in 100 mL glass Erlenmeyer flasks containing 40 mL mineral medium, supplemented with either fructose (0.1% w/v), BHET (0.1% w/v), terephthalic acid (0.02% w/v) or three 0.5 × 0.75 cm amorphous PET films (Goodfellow, UK) as sole sources of carbon and energy. Samples (1 mL) were taken on days 0, 1, 3, 7 and 14 to monitor growth through absorbance (600 nm) measurements. When there was visible growth, or the incubation time was 2 weeks (whichever was sooner), cultures were centrifuged (4000×g for 15 min). Cellular pellets were immediately stored at −20 °C for further proteomic analysis. Culture supernatants (~40 ml), used for exoproteomic analysis, were further filtered through 0.2 μm pore PTFE filters (Millex-LG; Millipore, Germany) prior to freezing. Before freezing, 1.5 mL of the filtered supernatant was separated in a different vial for metabolomic analysis. Ideonella sakaiensis 201-F6T was obtained from the National Institute of Technology and Evaluation Biological Resource Center (strain NBRC 110686T). I. sakaiensis cultures (in mineral medium with no additional NaCl) and no inoculum controls were incubated in parallel but were only used for absorbance measurements and metabolomics and not proteomics.
Global distribution of our isolates in marine plastispheres and planktonic samples
In order to determine the distribution of the two isolates, several searches were carried out: (i) local Blast  searches using the 16S rRNA gene of each isolate against the community succession amplicon sequencing data; (ii) local Blast searches using the 16S rRNA gene of each isolate against the 16S rRNA gene amplicon sequencing data from all marine plastisphere samples included in our recent plastisphere meta-analysis , i.e. 124,319 ASVs in 1185 samples; (iii) local Blast searches using the 16S rRNA gene of each isolate against 16S rRNA fragments from the Tara oceans metagenomes (mitags) [43,44,45] and (iv) MetaQUAST [102, 103] to determine the coverage for each of the isolates’ genomes within the assembled Tara metagenomes [43,44,45], downloaded from . For local Blast searches, all matches with above 90% identity were initially kept and these were then further filtered to keep only those with above 95, 97 or 99% identity, depending on the comparison being made. For plotting, all ASVs with above the chosen identity threshold (i.e. 95, 97 or 99%) were summed and for the map plots in Fig. 3 only the sample with the maximum sum of relative abundance was plotted for each location (5 × 5 latitude/longitude area). All custom scripts used for carrying out these searches as well as analysing and plotting these results are at https://github.com/R-Wright-1/PET-plastisphere.
Proteome preparation and shotgun analysis
Exoproteomes were concentrated using a trichloroacetic acid precipitation protocol as previously described . Exoproteome precipitates and cell pellets were then dissolved in 1 × LDS loading buffer (Invitrogen, USA) and further processed following the methods in . Tryptic digested proteomes were analysed by nanoLC-ESI MS/MS using an Ultimate 2000 LC system (Dionex-LC Packings) coupled with an Orbitrap Fusion mass spectrometer (Thermo Scientific, USA) using a 60 min LC separation (exoproteomes) or a 120 min LC separation (cellular proteomes) on a 25 cm column and settings as described in . Compiled mass spectra were identified and quantified using an in-house database comprised of the coding sequences (CDS) of each sequenced bacterium in MaxQuant version 220.127.116.11 . Comparative proteomic analysis between samples was performed using custom Python scripts (version 3.6.8 with modules numpy, os, csv and math; https://github.com/R-Wright-1/PET-plastisphere) written to carry out two-sample Student’s T tests for significance and calculate fold changes. Conserved domain searches  were carried out for manual curation of the functions assigned to all key proteins identified and further mapped onto KEGG degradation pathways [108, 109]. Peptides and protein groups are shown in Table S7 whilst the results of proteomic analyses including relative abundance and fold changes are shown in Table S8.
Metabolomics for the identification of PET degradation intermediates (LC-MS)
The untargeted metabolomic analyses of isolate and community supernatants were carried out using an UHPLC system (Ultimate 3000; Thermo Fisher Scientific, Waltham, MA, USA) coupled to a Q-Exactive Hybrid Quadrupole-Orbitrap mass spectrometer (ThermoFisher Scientific) operating with a heated electrospray interface (HESI). Sample preparation, LC-MS conditions and data processing were as previously described . The following pure compounds were used as internal standards: terephthalic acid (retention time 13.16 min; ∆ 1.20 ppm), MHET (retention time 13.81 min; ∆ 2.33 ppm) and BHET (retention time 14.08 min; ∆ 5.12 ppm). Full details of all detected compounds are in Tables S9 and S10 for the isolates and communities, respectively.
PET biodegradation by isolates and microbial communities
In parallel with the plastisphere succession experiment, microbial communities (obtained as above), isolates (Thioclava sp. BHET1 and Bacillus sp. BHET2) or controls with no inoculum were incubated in 25 cm2 tissue culture flasks containing 25 mL mineral medium with carbon sources as follows (n=24 total): (i) no additional carbon (control) or (ii) five films of 1 × 3 cm amorphous PET. Flasks were incubated for 5 months, after which the amorphous PET pieces were removed from each incubation. Flasks were topped up with sterile water as necessary to keep the volume at 25 mL.
Fourier transform infrared spectroscopy (FTIR) as an indication of PET biodegradation
Biofilms colonising amorphous PET films were removed as described in Erni-Cassola et al. . Briefly, colonised and non-colonised control films were soaked in 15% H2O2 and incubated at 60 °C for 90 min with shaking at 100 rpm, after which plastics were placed into fresh H2O2 and incubated overnight at 60 °C. Films were then thoroughly washed three times with MilliQ water and dried overnight at 60 °C. Following the removal of the biofilm, PET films were analysed by FTIR. Procedural controls that were not incubated with the microbial inoculums were included in the analysis. Three measurements were taken from each PET film using an Agilent Technologies Cary 630 FTIR spectrometer. FTIR spectra were smoothed and baseline normalised using the hyperSpec package in R . Ratios between the absorbance peak height at 1410 cm−1 and the absorbance peak heights at 725, 1090, 1240 and 1711 cm−1 for the C-H aromatic, C-O ester, C-O and C=O carboxylic acid bonds (I725/I1410, I1090/I1410, I1240/I1410 and I1711/I1410, respectively) [60, 61] and two sample T tests between the ratios of microbial incubations and control incubations with no inoculum were calculated using a custom python script (https://github.com/R-Wright-1/PET-plastisphere).
Availability of data and materials
All scripts used for analyses and adding additional genes to the PICRUSt2 database can be found at https://github.com/R-Wright-1/PET-Plastiphere. The JGI IMG genomes that were searched using the Hidden Markov Models are in the Figshare file at  (https://doi.org/10.6084/m9.figshare.12233192). MiSeq sequencing data are deposited in the NCBI Sequence Read Archive (SRA) under Bioproject accession number PRJNA544783. The complete genome sequences of Thioclava sp. BHET1 and Bacillus sp. BHET2 are deposited in the GenBank database under the accession numbers PRJNA544734 and PRJNA525098, respectively.
Law KL. Plastics in the marine environment. Ann Rev Mar Sci. 2017;9(1):205–29. https://doi.org/10.1146/annurev-marine-010816-060409.
PlasticsEurope. Plastics – the Facts 2018. 2018.
Geyer R, Jambeck JR, Law KL. Production, use, and fate of all plastics ever made. Sci Adv. 2017;3(7):e1700782. https://doi.org/10.1126/sciadv.1700782.
Gilman E. Status of international monitoring and management of abandoned, lost and discarded fishing gear and ghost fishing. Mar Policy. 2015;60:225–39. https://doi.org/10.1016/j.marpol.2015.06.016.
Napper IE, Thompson RC. Release of synthetic microplastic plastic fibres from domestic washing machines: effects of fabric type and washing conditions. Mar Pollut Bull. 2016;112(1-2):39–45. https://doi.org/10.1016/j.marpolbul.2016.09.025.
Napper IE, Bakir A, Rowland SJ, Thompson RC. Characterisation, quantity and sorptive properties of microplastics extracted from cosmetics. Mar Pollut Bull. 2015;99(1-2):178–85. https://doi.org/10.1016/j.marpolbul.2015.07.029.
Borrelle SB, Ringama J, Law KL, Monnahan CC, Lebreton L, McGivern A, et al. Predicted growth in plastic waste exceeds efforts to mitigate plastic pollution. Science (80- ). 2020;1518:1515–8.
Ward CP, Reddy CM. Opinion: we need better data about the environmental persistence of plastic goods. Proc Natl Acad Sci. 2020;117:202008009.
Avio CG, Gorbi S, Regoli F. Plastics and microplastics in the oceans: from emerging pollutants to emerged threat. Mar Environ Res. 2016;128:2–11.
Gerritse J, Leslie HA, de Tender CA, Devriese LI, Vethaak AD. Fragmentation of plastic objects in a laboratory seawater microcosm. Sci Rep. 2020;10:1–16.
Andrady AL. The plastic in microplastics: a review. Mar Pollut Bull. 2017;119(1):12–22. https://doi.org/10.1016/j.marpolbul.2017.01.082.
Wright RJ, Erni-Cassola G, Zadjelovic V, Latva M, Christie-Oleza J. Marine plastic debris – a new surface for microbial colonization. Environ Sci Technol. 2020;4(19):11657–72.
Danso D, Chow J, Streit WR. Plastics: Microbial Degradation, Environmental and Biotechnological Perspectives. Appl Environ Microbiol. 2019;85(19):e01095-19.
Tanasupawat S, Takehana T, Yoshida S, Hiraga K, Oda K. Ideonella sakaiensis sp. nov., isolated from a microbial consortium that degrades PET. Int J Syst Evol Microbiol. 2016;66(8):2813–8. https://doi.org/10.1099/ijsem.0.001058.
Yoshida S, Hiraga K, Takehana T, Taniguchi I, Yamaji H, Maeda Y, et al. A bacterium that degrades and assimilates poly(ethylene terephthalate). Science (80- ). 2016;351:1196–9.
Han X, Liu W, Huang J, Ma J, Zheng Y, Ko T, et al. Structural insight into catalytic mechanism of PET. Nat Commun. 2017;8:1–6.
Joo S, Cho IJ, Seo H, Son HF, Sagong H-Y, Shin TJ, et al. Structural insight into molecular mechanism of poly(ethylene terephthalate) degradation. Nat Commun. 2018;9:1–12.
Austin HP, Allen MD, Donohoe BS, Rorrer NA, Kearns FL, Silveira RL, et al. Characterization and engineering of a plastic-degrading aromatic polyesterase. Proc Natl Acad Sci. 2018;115(19):E4350–7. https://doi.org/10.1073/pnas.1718804115.
Danso D, Schmeisser C, Chow J, Zimmermann W, Wei R, Leggewie C, et al. New insights into the function and global distribution of polyethylene terephthalate (PET) degrading bacteria and enzymes in marine and terrestrial metagenomes. Appl Environ Microbiol. 2018;84:e02773–17.
Zettler ER, Mincer TJ, Amaral-Zettler LA. Life in the “Plastisphere”: microbial communities on plastic marine debris. Environ Sci Technol. 2013;47(13):7137–46. https://doi.org/10.1021/es401288x.
Bryant JA, Clemente TM, Viviani DA, Fong AA, Thomas KA, Kemp P, et al. Diversity and activity of communities inhabiting plastic debris in the North Pacific Gyre. mSystems. 2016;1:e00024–16.
Oberbeckmann S, Kreikemeyer B, Labrenz M. Environmental factors support the formation of specific bacterial assemblages on microplastics. Front Microbiol. 2018;8:1–12.
Frère L, Maignien L, Chalopin M, Huvet A, Rinnert E, Morrison H, et al. Microplastic bacterial communities in the Bay of Brest: influence of polymer type and size. Environ Pollut. 2018;242(Pt A):614–25. https://doi.org/10.1016/j.envpol.2018.07.023.
Kirstein IV, Wichels A, Gullans E, Krohne G, Gerdts G. The plastisphere – uncovering tightly attached plastic “specific” microorganisms. PLoS One. 2019;14:1–17.
Kirstein IV, Wichels A, Krohne G, Gerdts G. Mature biofilm communities on synthetic polymers in seawater - specific or general? Mar Environ Res. 2018;142:147–54. https://doi.org/10.1016/j.marenvres.2018.09.028.
Carson HS, Nerheim MS, Carroll KA, Eriksen M. The plastic-associated microorganisms of the North Pacific Gyre. Mar Pollut Bull. 2013;75(1-2):126–32. https://doi.org/10.1016/j.marpolbul.2013.07.054.
Jiang P, Zhao S, Zhu L, Li D. Microplastic-associated bacterial assemblages in the intertidal zone of the Yangtze Estuary. Sci Total Environ. 2018;624:48–54. https://doi.org/10.1016/j.scitotenv.2017.12.105.
Kettner MT, Rojas-Jimenez K, Oberbeckmann S, Labrenz M, Grossart H-P. Microplastics alter composition of fungal communities in aquatic ecosystems. Environ Microbiol. 2017;19(11):4447–59. https://doi.org/10.1111/1462-2920.13891.
Muthukrishnan T, Khaburi M Al, Abed RMM. Fouling microbial communities on plastics compared with wood and steel: Are they substrate- or location-specific? Microb Ecol. 2019;78:361–74.
Oberbeckmann S, Loeder MGJ, Gerdts G, Osborn AM. Spatial and seasonal variation in diversity and structure of microbial biofilms on marine plastics in Northern European waters. FEMS Microbiol Ecol. 2014;90(2):478–92. https://doi.org/10.1111/1574-6941.12409.
Oberbeckmann S, Osborn AM, Duhaime MB. Microbes on a bottle: substrate, season and geography influence community composition of microbes colonizing marine plastic debris. PLoS One. 2016;11:1–24.
Amaral-Zettler LA, Zettler ER, Slikas B, Boyd GD, Melvin DW, Morrall CE, et al. The biogeography of the plastisphere: implications for policy. Front Ecol Environ. 2015;13(10):541–6. https://doi.org/10.1890/150017.
Wright RJ, Langille MGI, Walker TR. Food or just a free ride? A meta-analysis reveals the global diversity of the Plastisphere. ISME J. 2021;15:789–806.
Pollet T, Berdjeb L, Garnier CC, Durrieu GG, Le Poupon C, Misson B, et al. Prokaryotic community successions and interactions in marine biofilms: the key role of Flavobacteriia. FEMS Microb Ecol. 2018;94:1–13.
Datta MS, Sliwerska E, Gore J, Polz MF, Cordero OX. Microbial interactions lead to rapid micro-scale successions on model marine particles. Nat Commun. 2016;7:1–7.
Wright RJ, Gibson MI, Christie-Oleza JA. Understanding microbial community dynamics to improve optimal microbiome selection. Microbiome. 2019;7:1–14.
Woodall LC, Jungblut AD, Hopkins K, Id AH, Robinson F, Gwinnett C, et al. Deep-sea anthropogenic macrodebris harbours rich and diverse communities of bacteria and archaea. PLoS One. 2018;13(11):e0206220.
De Tender CA, Devriese LI, Haegeman A, Maes S, Ruttink T, Dawyndt P. Bacterial community profiling of plastic litter in the Belgian part of the North Sea. Environ Sci Technol. 2015;49(16):9629–38. https://doi.org/10.1021/acs.est.5b01093.
Dang H, Li T, Chen M, Huang G. Cross-ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters. Appl Environ Microbiol. 2008;74(1):52–60. https://doi.org/10.1128/AEM.01400-07.
Erni-Cassola G, Wright RJ, Gibson MI, Christie-Oleza JA. Early colonization of weathered polyethylene by distinct bacteria in marine coastal seawater. Microb Ecol. Microbial Ecol. 2019;79:517–26.
van den Brink PJ, den Besten PJ. bij de Vaate A, ter Braak CJF. Principal response curves technique for the analysis of multivariate biomonitoring time series. Environ Monit Assess. 2009;152(1-4):271–81. https://doi.org/10.1007/s10661-008-0314-6.
Romera-Castillo C, Pinto M, Langer TM, Álvarez-Salgado XA, Herndl GJ. Dissolved organic carbon leaching from plastics stimulates microbial activity in the ocean. Nat Commun. 2018;9:1–7.
Sunagawa S, Coelho LP, Chaffron S, Kultima JR, Labadie K, Salazar G, et al. Structure and function of the global ocean microbiome. Science (80- ). 2015;348:1–10.
Logares R, Sunagawa S, Salazar G, Cornejo-Castillo FM, Ferrera I, Sarmento H, et al. Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities. Environ Microbiol. 2014;16(9):2659–71. https://doi.org/10.1111/1462-2920.12250.
Tully BJ, Graham ED, Heidelberg JF. The reconstruction of 2,631 draft metagenome-assembled genomes from the global oceans. Sci Data. 2018;5:1–8.
Wiegant WM, de Bont JA. A new route for ethylene glycol metabolism in Mycobacterium E44. J Gen Microbiol. 1980;120:325–31.
Trifunovic D, Schuchmann K, Muller V. Ethylene glycol metabolism in the acetogen Acetobacterium woodii. J Bacteriol. 2016;198(7):1058–65. https://doi.org/10.1128/JB.00942-15.
de Lara MS, Abdulmutalib U, Gonzalez J, Kim J, Smith AA, Faulon J-L, et al. Genes for a circular and sustainable bio-PET economy. Genes (Basel). 2019;10:1–15.
Hara H, Eltis LD, Davies JE, Mohn WW. Transcriptomic analysis reveals a bifurcated terephthalate degradation pathway in Rhodococcus sp. strain RHA1. J Bacteriol. 2007;189(5):1641–7. https://doi.org/10.1128/JB.01322-06.
Teufel R, Mascaraque V, Ismail W, Voss M, Perera J, Eisenreich W, et al. Bacterial phenylalanine and phenylacetate catabolic pathway revealed. Proc Natl Acad Sci U S A. 2010;107(32):14390–5. https://doi.org/10.1073/pnas.1005399107.
Arias-Barrau E, Olivera ER, Luengo JM, Fernández C, Galán B, García JL, et al. The homogentisate pathway: a central catabolic pathway involved in the degradation of L-phenylalanine, L-tyrosine, and 3-hydroxyphenylacetate in Pseudomonas putida. J Bacteriol. 2004;186(15):5062–77. https://doi.org/10.1128/JB.186.15.5062-5077.2004.
Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428(4):726–31. https://doi.org/10.1016/j.jmb.2015.11.006.
Douglas GM, Maffei VJ, Zaneveld JR, Yurgel SN, Brown JR, Taylor CM, et al. PICRUSt2 for prediction of metagenome functions. Nat Biotechnol. 2020;38:685–8.
Chen IMA, Chu K, Palaniappan K, Pillay M, Ratner A, Huang J, et al. IMG/M v.5.0: an integrated data management and comparative analysis system for microbial genomes and microbiomes. Nucleic Acids Res. 2019;47(D1):D666–77. https://doi.org/10.1093/nar/gky901.
Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:1814–21.
Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S, Chitsaz F, Geer LY, et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015;43:222–6.
Knott BC, Erickson E, Allen MD, Gado JE, Graham R, Kearns FL, et al. Characterization and engineering of a two-enzyme system for plastics depolymerization. Proc Natl Acad Sci. 2020;117:202006753.
Tanasupawat S, Takehana T, Yoshida S, Hiraga K, Oda K. Ideonella sakaiensis sp. nov., isolated from a microbial consortium that degrades poly(ethylene terephthalate). Int J Syst Evol Microbiol. 2016;66(8):2813–8. https://doi.org/10.1099/ijsem.0.001058.
Wright RJ, Bosch R, Gibson MI, Christie-Oleza JA. Plasticizer degradation by marine bacterial isolates: a proteogenomic and metabolomic characterization. Environ Sci Technol. 2020;54(4):2244–56. https://doi.org/10.1021/acs.est.9b05228.
Gok A. Degradation of poly(ethylene-terephthalate) under accelerated weathering exposures. Case West Reserv Univ. 2015. https://etd.ohiolink.edu/apexprod/rws_olink/r/1501/10?p10_etd_subid=110243&clear=10.
Donelli I, Freddi G, Nierstrasz VA, Taddei P. Surface structure and properties of poly-(ethylene terephthalate) hydrolyzed by alkali and cutinase. Polym Degrad Stab. 2010;95(9):1542–50. https://doi.org/10.1016/j.polymdegradstab.2010.06.011.
Quero GM, Luna GM. Surfing and dining on the “plastisphere”: microbial life on plastic marine debris. Adv Oceanogr Limnol. 2017;8:199–207.
Didier D, Anne M, Alexandra TH. Plastics in the North Atlantic garbage patch: a boat-microbe for hitchhikers and plastic degraders. Sci Total Environ. 2017;599–600:1222–32.
Dussud C, Meistertzheim AL, Conan P, Pujo-pay M, George M, Fabre P, et al. Evidence of niche partitioning among bacteria living on plastics, organic particles and surrounding seawaters. Environ Pollut. 2018;236:807–16. https://doi.org/10.1016/j.envpol.2017.12.027.
Wackett LP. Broad specificity microbial enzymes. Microb Biotechnol. 2015;8(1):188–9. https://doi.org/10.1111/1751-7915.12270.
Musumeci MA, Loviso CL, Lozada M, Ferreira FV, Dionisi HM. Substrate specificities of aromatic ring-hydroxylating oxygenases of an uncultured gammaproteobacterium from chronically-polluted subantarctic sediments. Int Biodeterior Biodegrad. 2019;137:127–36. https://doi.org/10.1016/j.ibiod.2018.12.005.
Syranidou E, Karkanorachaki K, Amorotti F, Avgeropoulos A, Kolvenbach B, Zhou N, et al. Biodegradation of mixture of plastic films by tailored marine consortia. J Hazard Mater. 2019;375:33–42. https://doi.org/10.1016/j.jhazmat.2019.04.078.
Syranidou E, Karkanorachaki K, Amorotti F, Repouskou E, Kroll K, Kolvenbach B, et al. Development of tailored indigenous marine consortia for the degradation of naturally weathered polyethylene films. PLoS One. 2017;12:1–21.
Syranidou E, Karkanorachaki K, Amorotti F, Franchini M, Repouskou E, Kaliva M, et al. Biodegradation of weathered polystyrene films in seawater microcosms. Sci Rep. 2017;7:1–12.
Curren E, Leong SCY. Profiles of bacterial assemblages from microplastics of tropical coastal environments. Sci Total Environ. 2019;655:313–20. https://doi.org/10.1016/j.scitotenv.2018.11.250.
Wang W, Zhong R, Shan D, Shao Z. Indigenous oil-degrading bacteria in crude oil-contaminated seawater of the Yellow sea, China. Appl Microbiol Biotechnol. 2014;98(16):7253–69. https://doi.org/10.1007/s00253-014-5817-1.
Nakkabi A, Elmoualij N, Sadiki M, Ibnsouda Koraichi S, Fahim M. Biodegradation of poly (ethylene terephthalate) by Bacillus Subtilis. Int J Recent Adv Multidiscip Res. 2015;2:1060–2.
Sudhakar M, Doble M, Murthy PS, Venkatesan R. Marine microbe-mediated biodegradation of low- and high-density polyethylenes. Int Biodeterior Biodegradation. 2008;61(3):203–13. https://doi.org/10.1016/j.ibiod.2007.07.011.
Gewert B, Plassmann MM, MacLeod M. Pathways for degradation of plastic polymers floating in the marine environment. Environ Sci Process Impacts. 2015;17(9):1513–21. https://doi.org/10.1039/C5EM00207A.
Gewert B, Plassmann M, Sandblom O, Macleod M. Identification of chain scission products released to water by plastic exposed to ultraviolet light. Environ Sci Technol Lett. Am Chem Soc. 2018;5:272–6.
Freire MTDA, Damant AP, Castle L, Reyes FGR. Thermal stability of polyethylene terephthalate (PET): oligomer distribution and formation of volatiles. Packag Technol Sci. 1999;12(1):29–36.
Kirstein IV, Kirmizi S, Wichels A, Garin-Fernandez A, Erler R, Löder M, et al. Dangerous hitchhikers? Evidence for potentially pathogenic Vibrio spp. on microplastic particles. Mar Environ Res. 2016;120:1–8. https://doi.org/10.1016/j.marenvres.2016.07.004.
Quilliam RS, Jamieson J, Oliver DM. Seaweeds and plastic debris can influence the survival of faecal indicator organisms in beach environments. Mar Pollut Bull. 2014;84(1-2):201–7. https://doi.org/10.1016/j.marpolbul.2014.05.011.
Keswani A, Oliver DM, Gutierrez T, Quilliam RS. Microbial hitchhikers on marine plastic debris: human exposure risks at bathing waters and beach environments. Mar Environ Res. 2016;118:10–9. https://doi.org/10.1016/j.marenvres.2016.04.006.
Huq A, Small EB, West PA, Huq MI, Rahman R, Colwell RR. Ecological relationships between Vibrio cholerae and planktonic crustacean copepods. Appl Environ Microbiol. 1983;45(1):275–83. https://doi.org/10.1128/AEM.45.1.275-283.1983.
Raghul SS, Bhat SG, Chandrasekaran M, Francis V, Thachil ET. Biodegradation of polyvinyl alcohol-low linear density polyethylene-blended plastic film by consortium of marine benthic vibrios. Int J Environ Sci Technol. 2014;11(7):1827–34. https://doi.org/10.1007/s13762-013-0335-8.
Noriega-Ortega BE, Wienhausen G, Mentges A, Dittmar T, Simon M, Niggemann J. Does the chemodiversity of bacterial exometabolomes sustain the chemodiversity of marine dissolved organic matter? Front Microbiol. 2019;10:1–13.
Christie-Oleza JA, Sousoni D, Lloyd M, Armengaud J, Scanlan DJ. Nutrient recycling facilitates long-term stability of marine microbial phototroph-heterotroph interactions. Nat Microbiol. 2017;2:1–10.
Zhang WP, Wang Y, Tian RM, Bougouffa S, Yang B, Cao HL, et al. Species sorting during biofilm assembly by artificial substrates deployed in a cold seep system. Sci Rep. 2014;4:1–7.
Ogonowski M, Motiei A, Ininbergs K, Hell E, Gerdes Z, Udekwu KI, et al. Evidence for selective bacterial community structuring on microplastics. Environ Microbiol. 2018;20(8):2796–808. https://doi.org/10.1111/1462-2920.14120.
Zumstein MT, Narayan R, Kohler H-PE, McNeill K, Sander M. Dos and do nots when assessing the biodegradation of plastics. Environ Sci Technol. 2019;53(17):9967–9. https://doi.org/10.1021/acs.est.9b04513.
Bushnell LD, Haas HF. The utilization of certain hydrocarbons by microorganisms. J Bacteriol. 1941;41(5):653–73. https://doi.org/10.1128/JB.41.5.653-673.1941.
Parada AE, Needham DM, Fuhrman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol. 2016;18(5):1403–14. https://doi.org/10.1111/1462-2920.13023.
Callahan BJ, Sankaran K, Fukuyama JA, McMurdie PJ, Holmes SP. Bioconductor workflow for microbiome data analysis: from raw reads to community analyses. F1000Research. 2016;5:1–49.
Louca S, Doebeli M. Efficient comparative phylogenetics on large trees. Bioinformatics. 2018;34(6):1053–5. https://doi.org/10.1093/bioinformatics/btx701.
Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011;7(10):e1002195. https://doi.org/10.1371/journal.pcbi.1002195.
Barbera P, Kozlov AM, Czech L, Morel B, Darriba D, Flouri T, et al. EPA-ng: massively parallel evolutionary placement of genetic sequences. Syst Biol. 2019;68(2):365–9. https://doi.org/10.1093/sysbio/syy054.
Czech L, Stamatakis A. Scalable methods for analyzing and visualizing phylogenetic placement of metagenomic samples. PLoS One. 2019;14(5):e0217050.
Wright RJ. JGI PICRUSt2 genomes [Internet]. Figshare. Available from: https://doi.org/10.6084/m9.figshare.12233192.
Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7(1):539. https://doi.org/10.1038/msb.2011.75.
Madeira F, Park YM, Lee J, Buso N, Gur T, Madhusoodanan N, et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019;47(W1):W636–41. https://doi.org/10.1093/nar/gkz268.
Paradis E, Schliep K. Ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics. 2019;35(3):526–8. https://doi.org/10.1093/bioinformatics/bty633.
Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14(6):927–30. https://doi.org/10.1111/j.1654-1103.2003.tb02228.x.
Lane DJ. 16S/23S rRNA Sequencing. In: Stackebrandt E, Goodfellow M, editors. Nucleic Acid Tech Bact Syst; 1991. p. 115–75.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9. https://doi.org/10.1093/bioinformatics/btu153.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10. https://doi.org/10.1016/S0022-2836(05)80360-2.
Mikheenko A, Prjibelski A, Saveliev V, Antipov D, Gurevich A. Versatile genome assembly evaluation with QUAST-LG. Bioinformatics. 2018;34(13):i142–50. https://doi.org/10.1093/bioinformatics/bty266.
Mikheenko A, Saveliev V, Gurevich A. MetaQUAST: Evaluation of metagenome assemblies. Bioinformatics. 2016;32(7):1088–90. https://doi.org/10.1093/bioinformatics/btv697.
Eren AM. TARA-RAW-ASSEMBLIES-1000nt. [Internet]. Figshare. Available from: https://doi.org/10.6084/m9.figshare.4902920.v1.
Christie-Oleza JA, Armengaud J. In-depth analysis of exoproteomes from marine bacteria by shotgun liquid chromatography-tandem mass spectrometry: the Ruegeria pomeroyi DSS-3 case-study. Mar Drugs. 2010;8(8):2223–39. https://doi.org/10.3390/md8082223.
Christie-Oleza JA, Scanlan DJ, Armengaud J. “You produce while I clean up”, a strategy revealed by exoproteomics during Synechococcus-Roseobacter interactions. Proteomics. 2015;15(20):3454–62. https://doi.org/10.1002/pmic.201400562.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26(12):1367–72. https://doi.org/10.1038/nbt.1511.
Marchler-Bauer A, Bo Y, Han L, He J, Lanczycki CJ, Lu S, et al. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017;45:200–3.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61. https://doi.org/10.1093/nar/gkw1092.
Erni-Cassola G, Gibson MI, Thompson RC, Christie-Oleza JA. Lost, but found with Nile Red: a novel method for detecting and quantifying small microplastics (1 mm to 20 μm) in environmental samples. Environ Sci Technol. 2017;51(23):13641–8. https://doi.org/10.1021/acs.est.7b04512.
Beleites C, Sergo V. hyperSpec: a package to handle hyperspectral data sets in R; 2020.
We thank Andrew Millard and Slawomir Michniewski for assistance with running the MiSeq, Maria del Mar Aguilo-Ferretjans for help with proteomics sample preparation, Alex Baker and Christopher Stubbs for assistance with the μFTIR, and the Christie-Oleza and Langille groups (in particular Vinko Zadjelovic, Audam Chhun, Gabriel Erni Cassola and Gavin Douglas) for helpful discussions throughout the project. We also acknowledge technical support from the WPH Proteomic Facility at the University of Warwick and thank the Warwick Institute for Synthetic Biology for access to the GeneVac. Technical assistance in processing the metabolomics samples by LC-Orbitrap-MS/MS provided by Gabriel Martorell and Rosa Gomila, both from the Scientific and Technical Center (SCT) of the University of the Balearic Islands, is also greatly appreciated.
RW was supported by a Midlands Integrative Biosciences Training Partnership PhD scholarship via grant BB/M01116X/1, JC-O by the Natural Environment Research Council Independent Research Fellowship NE/K009044/1 and Ramón y Cajal contract RYC-2017-22452 (funded by the Ministry of Science, Innovation and Universities, the National Agency of Research, and the European Social Fund), RB by the MINECO project CTM2015-70180-R (FEDER cofunding) and MIG by European Research Council grant 638631. This work was supported by the Agencia Estatal de Investigación (project PID2019-109509RB-I00 / AEI / 10.13039/501100011033).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
VideoID for videobyte was corrected.
DNA yields from all PET succession experiment mesocosms and the inoculum. The DNA concentration in all negative extraction controls was too low to measure aside from day 30 (0.04 ng μL-1). Where DNA yields were significantly (two sample T-test p<0.05) higher than the no carbon control these days are marked with an asterisk. Figure S2. Relative abundances of taxa within all samples, with each bar representing the mean of three biological replicates. Each row shows taxa grouped to a different taxonomic level (shown on y label) and other represents all that were present at below 0.5% relative abundance at that level. ASVs are classified to species level where possible: ASV1 Vibrio, ASV2 Alteromonas, ASV3 Bacillus, ASV4 Vibrio, ASV5 Sediminibacterium salmoneum, ASV6 Pseudoalteromonas, ASV7 Sunxiuqinia, ASV8 Alcanivorax, ASV9 Thiobacimonas profunda, ASV10 Thalassospira lucentensis, ASV11 Methylophaga, ASV12 Pseudoalteromonas, ASV13 Vibrio, ASV14 Thalassospira, ASV15 Alcanivorax, ASV16 Thalassospira, ASV17 Vibrio alginolyticus, ASV18 Halomonas, ASV19 Catenococcus, ASV20 Pseudomonas, ASV21 Pseudoalteromonas, ASV22 Thalassospira, ASV23 Shewanella, ASV26 Exiguobacterium, ASV28 Roseivirga, ASV31 Alteromonas, ASV34 Oricola cellulosilytica, ASV36 Lysobacter maris, ASV39 Catenococcus, ASV40 Exiguobacterium, ASV42 Tistlia, ASV43 Azomonas, ASV45 Rhizobiaceae, ASV46 Maritimibacter, ASV48 Catenococcus, ASV49 Halomonadaceae, ASV53 Vibrionaceae, ASV54 Oricola, ASV55 Catenococcus, ASV56 Sphingomonas, ASV63 Vibrionaceae, ASV64 Parvibaculum, ASV65 Catenococcus, ASV70 Catenococcus, ASV73 Catenococcus, ASV75 Vibrionaceae, ASV76 Aestuariibacter aggregatus, ASV78 Catenococcus, ASV85 Hyphobacterium, ASV145 Tenacibaculum litoreum, ASV153 Mesoflavibacter zeaxanthinifaciens, ASV203 Tenacibaculum. Figure S3. Diversity for all samples across 42 days of incubation. Showing Simpsons index of diversity (top) and species richness (bottom). Figure S4. Fourier transform infrared spectra of PET powder and weathered PET powder before incubation with communities or isolates. All wavelengths are shown in (A) and ratios between the wavenumbers at 1410 and 1711, 1240, 725 and 1090 cm-1 are shown in (B), while (C), (D), (E) and (F) show 1000-600, 1200-1000, 1600-1200 and 1800-1600 cm-1, respectively. Each line shows the mean absorbance for three technical replicates for each of three biological replicates (i.e. nine total measurements) per treatment while bars and error bars represent means and standard deviations for biological replicates. Asterisks denote significant differences (two independent samples T-test, p<0.05) between the wavenumber ratios before and after thermal weathering of PET. Dashed lines indicate the wavenumbers used for ratio calculations. Figure S5. Distribution of Thioclava sp. BHET1 (left) and Bacillus sp. BHET2 (right) in surface (top) or deep-chlorophyll maximum (bottom) waters samples by the Tara oceans expedition [2,3,4]. Sequences within the Tara oceans mitags dataset that shared above 90, 95 or 97% identity with each of the Thioclava sp. BHET1 and Bacillus sp. BHET2 16S rRNA genes were identified and the relative abundance of all matches were summed to give the abundances shown here. We also calculated the coverage for each of Thioclava sp. BHET1 and Bacillus sp. BHET2 in the assembled Tara oceans metagenomes. These were co-assembled for each ocean and the same values are therefore plotted for all stations within each ocean. For both relative abundance and coverage, purple indicates 0% and green indicates 3% or higher. Figure S6. Growth of the isolates Thioclava sp. BHET1 (A-C) and Bacillus sp. BHET2 (D-F) on a range of common growth substrates across three days of incubation. Panels show biological replicates. Figure S7. Predicted abundance of genes that are potentially involved in PET degradation in PICRUSt2-assembled predicted metagenomes for all communities over time. Genes that are in the standard PET degradation pathway (i.e. shown in Fig. 4) are shown in red in (A), while those that are predicted to be involved in PET degradation based on the proteomics results of isolates (Fig. 3) are shown in blue. The abundance and taxonomic contributions to each KEGG ortholog is shown in (B). Note that K14037 was not found in the PICRUSt2 predicted metagenome and is therefore not shown here. Also shown is all monooxygenases and all dioxygenases. Taxonomic contributions shown are scaled by the relative abundance of each taxon as well as the number of gene copies possessed by that taxon. All taxa with a total contribution below 0.5% are grouped to Other. See Table S2 for individual ASVs and Table S12 for NSTI values for all treatments at all time points. Figure S7. Predicted abundance of genes that are potentially involved in PET degradation in PICRUSt2-assembled predicted metagenomes for all communities over time. Genes that are in the standard PET degradation pathway (i.e. shown in Fig. 4) are shown in red in (A), while those that are predicted to be involved in PET degradation based on the proteomics results of isolates (Fig. 3) are shown in blue. The abundance and taxonomic contributions to each KEGG ortholog is shown in (B). Note that K14037 was not found in the PICRUSt2 predicted metagenome and is therefore not shown here. Also shown is all monooxygenases and all dioxygenases. Taxonomic contributions shown are scaled by the relative abundance of each taxon as well as the number of gene copies possessed by that taxon. All taxa with a total contribution below 0.5% are grouped to Other. See Table S2 for individual ASVs and Table S12 for NSTI values for all treatments at all time points. Figure S8. nMDS plot showing Bray-Curtis distance between metabolomic analyses of culture supernatants of the community incubations used for MiSeq on day 42. Supernatants from incubations with no inoculum are shown with crosses while supernatants from incubations with the microbial community are shown with circles. Table S1. Results of PERMANOVA and ANOSIM tests for statistical significance (using Bray-Curtis distance) on all succession experiment samples. ANOSIM results are mentioned in the text as these values are more conservative. Table S2. ASVs identified by the PRC analysis. This includes ASV classifications using DADA2 and BLAST, PRC species weights, the closest representative whole genomes (from the NCBI database, where this was >97% similarity), whether these genomes potentially contain PETases and MHETases, PICRUSt2 nearest sequenced taxon indices (NSTI), the KEGG orthologs present in the genomes according to PICRUSt2 and other relevant information. Table S3. Analysis of early, middle or late colonisers, showing the day on which that ASV was most abundant in that treatment. Only ASVs that were above 0.5% in abundance in at least one time point in that treatment were included. Table S4. Genomic analysis of Thioclava sp. BHET1 and Bacillus sp. BHET2 (separate excel file). Table S5. Sequences that were used to construct the Hidden Markov Models (HMMs) for PETase, pcaG, pcaH, tphA2, tphA3 and tphB. Table S6. Potential PETases found in the genomes of Thioclava sp. BHET1 and Bacillus sp. BHET2 using a Hidden Markov Model (HMM) constructed with known PETases. Table S7. Peptides and protein groups for cellular and extracellular proteomics performed on Thioclava sp. BHET1 and Bacillus sp. BHET2 growing with fructose, TPA, BHET and PET (separate excel file). Table S8. Proteomic analysis of cellular and extracellular proteomics performed on Thioclava sp. BHET1 and Bacillus sp. BHET2 growing with fructose, TPA, BHET and PET (separate excel file). Table S9. Metabolomic analysis performed on Thioclava sp. BHET1 and Bacillus sp. BHET2 growing with fructose, TPA, BHET and PET (separate excel file). Table S10. Metabolomic analysis performed on microbial communities after incubation with BHET, amorphous PET, PET powder and weathered PET powder (separate excel file). Table S11. Proteins in the Thioclava sp. BHET1 cellular proteome that are potentially related to PET, BHET and TPA degradation, including relative abundance within the proteome and fold change when compared with the positive control. Table S12. Proteins that are potentially involved in xenobiotics degradation that were upregulated in one or more treatments in the Bacillus sp. BHET2 cellular proteome, including relative abundance within the proteome and fold change when compared with the positive control. Table S13. Details of the PETases found within the PICRUSt2 artificial metagenome and predicted to be in ASVs (separate excel file). Table S14. Nearest Sequenced Taxon Indices (NSTI) for all samples included in the PICRUSt2 analysis. Table S15. Comparison of different kits for DNA extraction from plastic pieces incubated with microbial communities.
About this article
Cite this article
Wright, R.J., Bosch, R., Langille, M.G.I. et al. A multi-OMIC characterisation of biodegradation and microbial community succession within the PET plastisphere. Microbiome 9, 141 (2021). https://doi.org/10.1186/s40168-021-01054-5
- Polyethylene terephthalate
- Plastic biodegradation
- Microbial community succession