Additional file 2: Table S1. Flow cytometric counts (cells mL-1) and forward scatter (FSC) of Synechococcus and two photosynthetic eukaryotic populations along with the ratio of the Synechococcus to total photosynthetic eukaryote populations. Table S2. Flow cytometric counts (cells mL-1) of bacterial and viral subpopulations and total populations along with the ratio of total virus populations to total bacterial populations. Table S3. Statistics for changes in flow cytometrically-quantified population abundances (count) of various populations at experimental day 6 for different treatments. Table S4. Measurements of photosynthetic efficiency for each treatment during the 6 days of the experiment. Table S5. (a) Concentrations of all of the metal tested with ICP-MS. Metal concentration is espressed in mg/L. Highlighted in green are the metal for which the concentration was above the limit of the detection for at least one of the treatment or for the seawater control. ICP-MS analysis was performed at the Elemental Analysis Facility at the Southern Cross University in Lismore,NSW (Australia) an accredited NATA facility for the analysis of seawater samples. Methods reference has the code for the NATA standard methodology followed. T0 indicates samples that were collected at the beginning of the experiment, day 6 samples that were collected at the end of the experiment . *PQL= Minimum dectection limit. (b) Concentrations of all of the metal tested with ICP-MS for the preliminary batch of PVC. Metal concentration is espressed in ug/L. Highlighted in green are the metal for which the concentration was above the limit of the detection for at least one of the replicate. ICP-MS analysis was performed at the Elemental Analysis Facility at the Southern Cross University in Lismore,NSW (Australia) an accredited NATA facility for the analysis of seawater samples. Table S6. Statistical test results for community profile groupings based on amplicon analyses. Table S7. Alpha diversity calculations for the bacterial and eukaryotic communities based on amplicon analyses. Table S8. DESeq2 results for the bacterial community. Only results with an alpha < 0.01 were chosen. Table S9. Anova followed by pairwise t-test FDR corrected calculated for the photosynthetic groups identified based on cyanobacteria and plastid sequences from 16S rRNA. Table S10. DESeq2 results for the eukaryotic community. Only results with an alpha < 0.01 were chosen. Table S11. Taxonomic assignment of reads for genes in the bacterial metagenomes and the statistical analyses of the taxonomic groups in each treatment. Table S12. DESeq2 results for PVC10 vs SW on day 6. Only KEGG orthologs with log2Fold change >2 and < -2 are displayed. Table S13. DESeq2 results for PVC1 vs SW on day 6. Only KEGG orthologs with log2Fold change >2 and < -2 are displayed. Table S14. DESeq2 results for ZnH vs SW on day 6. Only KEGG orthologs with log2Fold change >2 and < -2 are displayed. Table S15. DESeq2 results for ZnL vs SW on day 6. Only KEGG orthologs with log2Fold change >2 and < -2 are displayed. Table S16. Pairwise comparison of the abundance of integrase genes in the different treatments. Only significant comparisons (p.adj < 0.05) are reported here. Significance was calculated with an ANOVA followed by a pairwise t-test and the p value FDR adjusted. n1 and n2 represents the number of replicates for treatment group 1 and group 2 respectively. Table S17. (a) Summary information on the assembled metagenome-assembled genomes (bins) and MAGs coverage in the different samples. (b) Pairwise test on the MAGs distribution between the different treatments, p-value was calculate in R with the r-statix package and FDR corrected. Table S18. Pairwise comparison of the coverage of almost complete metabolic pathway modules of the identified MAGs (completion >49%) across the different treatments. Significance was calculated with an ANOVA followed by a pairwise t-test and the pvalue FDR adjusted. n1 and n2 represents the number of replicates for group 1 and group 2 respectively.