A marine heatwave drives significant shifts in pelagic microbiology.

Marine heatwaves (MHWs) cause disruption to marine ecosystems, deleteriously impacting macroflora and fauna. However, effects on microorganisms are relatively unknown despite ocean temperature being a major determinant of assemblage structure. Using data from thousands of Southern Hemisphere samples, we reveal that during an "unprecedented" 2015/16 Tasman Sea MHW, temperatures approached or surpassed the upper thermal boundary of many endemic taxa. Temperate microbial assemblages underwent a profound transition to niche states aligned with sites over 1000 km equatorward, adapting to higher temperatures and lower nutrient conditions bought on by the MHW. MHW conditions also modulate seasonal patterns of microbial diversity and support novel assemblage compositions. The most significant affects of MHWs on microbial assemblages occurred during warmer months, when temperatures exceeded the upper climatological bounds. Trends in microbial response across several MHWs in different locations suggest these are emergent properties of temperate ocean warming, which may facilitate monitoring, prediction and adaptation efforts.

The Microbiological Drivers of Temporally Dynamic Dimethylsulfoniopropionate Cycling Processes in Australian Coastal Shelf Waters.

The organic sulfur compounds dimethylsulfoniopropionate (DMSP) and dimethyl sulfoxide (DMSO) play major roles in the marine microbial food web and have substantial climatic importance as sources and sinks of dimethyl sulfide (DMS). Seasonal shifts in the abundance and diversity of the phytoplankton and bacteria that cycle DMSP are likely to impact marine DMS (O) (P) concentrations, but the dynamic nature of these microbial interactions is still poorly resolved. Here, we examined the relationships between microbial community dynamics with DMS (O) (P) concentrations during a 2-year oceanographic time series conducted on the east Australian coast. Heterogenous temporal patterns were apparent in chlorophyll a (chl a) and DMSP concentrations, but the relationship between these parameters varied over time, suggesting the phytoplankton and bacterial community composition were affecting the net DMSP concentrations through differential DMSP production and degradation. Significant increases in DMSP were regularly measured in spring blooms dominated by predicted high DMSP-producing lineages of phytoplankton (Heterocapsa, Prorocentrum, Alexandrium, and Micromonas), while spring blooms that were dominated by predicted low DMSP-producing phytoplankton (Thalassiosira) demonstrated negligible increases in DMSP concentrations. During elevated DMSP concentrations, a significant increase in the relative abundance of the key copiotrophic bacterial lineage Rhodobacterales was accompanied by a three-fold increase in the gene, encoding the first step of DMSP demethylation (dmdA). Significant temporal shifts in DMS concentrations were measured and were significantly correlated with both fractions (0.2-2 µm and >2 µm) of microbial DMSP lyase activity. Seasonal increases of the bacterial DMSP biosynthesis gene (dsyB) and the bacterial DMS oxidation gene (tmm) occurred during the spring-summer and coincided with peaks in DMSP and DMSO concentration, respectively. These findings, along with significant positive relationships between dsyB gene abundance and DMSP, and tmm gene abundance with DMSO, reinforce the significant role planktonic bacteria play in producing DMSP and DMSO in ocean surface waters. Our results highlight the highly dynamic nature and myriad of microbial interactions that govern sulfur cycling in coastal shelf waters and further underpin the importance of microbial ecology in mediating important marine biogeochemical processes.

Diatom Biogeography, Temporal Dynamics, and Links to Bacterioplankton across Seven Oceanographic Time-Series Sites Spanning the Australian Continent.

Diatom communities significantly influence ocean primary productivity and carbon cycling, but their spatial and temporal dynamics are highly heterogeneous and are governed by a complex diverse suite of abiotic and biotic factors. We examined the seasonal and biogeographical dynamics of diatom communities in Australian coastal waters using amplicon sequencing data (18S-16S rRNA gene) derived from a network of oceanographic time-series spanning the Australian continent. We demonstrate that diatom community composition in this region displays significant biogeography, with each site harbouring distinct community structures. Temperature and nutrients were identified as the key environmental contributors to differences in diatom communities at all sites, collectively explaining 21% of the variability observed in diatoms assemblages. However, specific groups of bacteria previously implicated in mutualistic ecological interactions with diatoms (Rhodobacteraceae, Flavobacteriaceae and Alteromonadaceae) also explained a further 4% of the spatial dynamics observed in diatom community structure. We also demonstrate that the two most temperate sites (Port Hacking and Maria Island) exhibited strong seasonality in diatom community and that at these sites, winter diatom communities co-occurred with higher proportion of Alteromonadaceae. In addition, we identified significant co-occurrence between specific diatom and bacterial amplicon sequence variants (ASVs), with members of the Roseobacter and Flavobacteria clades strongly correlated with some of the most abundant diatom genera (Skeletonema, Thalassiosira, and Cylindrotheca). We propose that some of these co-occurrences might be indicative of ecologically important interactions between diatoms and bacteria. Our analyses reveal that in addition to physico-chemical conditions (i.e., temperature, nutrients), the relative abundance of specific groups of bacteria appear to play an important role in shaping the spatial and temporal dynamics of marine diatom communities.

Biogeography of Southern Ocean prokaryotes: a comparison of the Indian and Pacific sectors.

We investigated the Southern Ocean (SO) prokaryote community structure via zero-radius operational taxonomic unit (zOTU) libraries generated from 16S rRNA gene sequencing of 223 full water column profiles. Samples reveal the prokaryote diversity trend between discrete water masses across multiple depths and latitudes in Indian (71-99°E, summer) and Pacific (170-174°W, autumn-winter) sectors of the SO. At higher taxonomic levels (phylum-family) we observed water masses to harbour distinct communities across both sectors, but observed sectorial variations at lower taxonomic levels (genus-zOTU) and relative abundance shifts for key taxa such as Flavobacteria, SAR324/Marinimicrobia, Nitrosopumilus and Nitrosopelagicus at both epi- and bathy-abyssopelagic water masses. Common surface bacteria were abundant in several deep-water masses and vice-versa suggesting connectivity between surface and deep-water microbial assemblages. Bacteria from same-sector Antarctic Bottom Water samples showed patchy, high beta-diversity which did not correlate well with measured environmental parameters or geographical distance. Unconventional depth distribution patterns were observed for key archaeal groups: Crenarchaeota was found across all depths in the water column and persistent high relative abundances of common epipelagic archaeon Nitrosopelagicus was observed in deep-water masses. Our findings reveal substantial regional variability of SO prokaryote assemblages that we argue should be considered in wide-scale SO ecosystem microbial modelling.

Biogeographical and seasonal dynamics of the marine Roseobacter community and ecological links to DMSP-producing phytoplankton.

Ecological interactions between marine bacteria and phytoplankton play a pivotal role in governing the ocean's major biogeochemical cycles. Among these, members of the marine Roseobacter Group (MRG) can establish mutualistic relationships with phytoplankton that are, in part, maintained by exchanges of the organosulfur compound, dimethylsulfoniopropionate (DMSP). Yet most of what is known about these interactions has been derived from culture-based laboratory studies. To investigate temporal and spatial co-occurrence patterns between members of the MRG and DMSP-producing phytoplankton we analysed 16S and 18S rRNA gene amplicon sequence variants (ASVs) derived from 5 years of monthly samples from seven environmentally distinct Australian oceanographic time-series. The MRG and DMSP-producer communities often displayed contemporaneous seasonality, which was greater in subtropical and temperate environments compared to tropical environments. The relative abundance of both groups varied latitudinally, displaying a poleward increase, peaking (MRG at 33% of total bacteria, DMSP producers at 42% of eukaryotic phototrophs) during recurrent spring-summer phytoplankton blooms in the most temperate site (Maria Island, Tasmania). Network analysis identified 20,140 significant positive correlations between MRG ASVs and DMSP producers and revealed that MRGs exhibit significantly stronger correlations to high DMSP producers relative to other DMSP-degrading bacteria (Pelagibacter, SAR86 and Actinobacteria). By utilising the power of a continental network of oceanographic time-series, this study provides in situ confirmation of interactions found in laboratory studies and demonstrates that the ecological dynamics of an important group of marine bacteria are shaped by the production of an abundant and biogeochemically significant organosulfur compound.

Metabolic pathways inferred from a bacterial marker gene illuminate ecological changes across South Pacific frontal boundaries.

Global oceanographic monitoring initiatives originally measured abiotic essential ocean variables but are currently incorporating biological and metagenomic sampling programs. There is, however, a large knowledge gap on how to infer bacterial functions, the information sought by biogeochemists, ecologists, and modelers, from the bacterial taxonomic information (produced by bacterial marker gene surveys). Here, we provide a correlative understanding of how a bacterial marker gene (16S rRNA) can be used to infer latitudinal trends for metabolic pathways in global monitoring campaigns. From a transect spanning 7000 km in the South Pacific Ocean we infer ten metabolic pathways from 16S rRNA gene sequences and 11 corresponding metagenome samples, which relate to metabolic processes of primary productivity, temperature-regulated thermodynamic effects, coping strategies for nutrient limitation, energy metabolism, and organic matter degradation. This study demonstrates that low-cost, high-throughput bacterial marker gene data, can be used to infer shifts in the metabolic strategies at the community scale.

Temperate southern Australian coastal waters are characterised by surprisingly high rates of nitrogen fixation and diversity of diazotrophs.

Biological dinitrogen (N2) fixation is one mechanism by which specific microorganisms (diazotrophs) can ameliorate nitrogen (N) limitation. Historically, rates of N2 fixation were believed to be limited outside of the low nutrient tropical and subtropical open ocean; however, emerging evidence suggests that N2 fixation is also a significant process within temperate coastal waters. Using a combination of amplicon sequencing, targeting the nitrogenase reductase gene (nifH), quantitative nifH PCR, and 15N2 stable isotope tracer experiments, we investigated spatial patterns of diazotroph assemblage structure and N2 fixation rates within the temperate coastal waters of southern Australia during Austral autumn and summer. Relative to previous studies in open ocean environments, including tropical northern Australia, and tropical and temperate estuaries, our results indicate that high rates of N2 fixation (10-64 nmol L-1 d-1) can occur within the large inverse estuary Spencer Gulf, while comparatively low rates of N2 fixation (2 nmol L-1 d-1) were observed in the adjacent continental shelf waters. Across the dataset, low concentrations of NO3/NO2 were significantly correlated with the highest N2 fixation rates, suggesting that N2 fixation could be an important source of new N in the region as dissolved inorganic N concentrations are typically limiting. Overall, the underlying diazotrophic community was dominated by nifH sequences from Cluster 1 unicellular cyanobacteria of the UCYN-A clade, as well as non-cyanobacterial diazotrophs related to Pseudomonas stutzeri, and Cluster 3 sulfate-reducing deltaproteobacteria. Diazotroph community composition was significantly influenced by salinity and SiO4 concentrations, reflecting the transition from UCYN-A-dominated assemblages in the continental shelf waters, to Cluster 3-dominated assemblages in the hypersaline waters of the inverse estuary. Diverse, transitional diazotrophic communities, comprised of a mixture of UCYN-A and putative heterotrophic bacteria, were observed at the mouth and southern edge of Spencer Gulf, where the highest N2 fixation rates were observed. In contrast to observations in other environments, no seasonal patterns in N2 fixation rates and diazotroph community structure were apparent. Collectively, our findings are consistent with the emerging view that N2 fixation within temperate coastal waters is a previously overlooked dynamic and potentially important component of the marine N cycle.

Microbial tropicalization driven by a strengthening western ocean boundary current.

Western boundary currents (WBCs) redistribute heat and oligotrophic seawater from the tropics to temperate latitudes, with several displaying substantial climate change-driven intensification over the last century. Strengthening WBCs have been implicated in the poleward range expansion of marine macroflora and fauna, however, the impacts on the structure and function of temperate microbial communities are largely unknown. Here we show that the major subtropical WBC of the South Pacific Ocean, the East Australian Current (EAC), transports microbial assemblages that maintain tropical and oligotrophic (k-strategist) signatures, to seasonally displace more copiotrophic (r-strategist) temperate microbial populations within temperate latitudes of the Tasman Sea. We identified specific characteristics of EAC microbial assemblages compared with non-EAC assemblages, including strain transitions within the SAR11 clade, enrichment of Prochlorococcus, predicted smaller genome sizes and shifts in the importance of several functional genes, including those associated with cyanobacterial photosynthesis, secondary metabolism and fatty acid and lipid transport. At a temperate time-series site in the Tasman Sea, we observed significant reductions in standing stocks of total carbon and chlorophyll a, and a shift towards smaller phytoplankton and carnivorous copepods, associated with the seasonal impact of the EAC microbial assemblage. In light of the substantial shifts in microbial assemblage structure and function associated with the EAC, we conclude that climate-driven expansions of WBCs will expand the range of tropical oligotrophic microbes, and potentially profoundly impact the trophic status of temperate waters.

Investigating the Diversity of Marine Bacteriophage in Contrasting Water Masses Associated with the East Australian Current (EAC) System.

Virus- and bacteriophage-induced mortality can have a significant impact on marine productivity and alter the flux of nutrients in marine microbial food-webs. Viral mediated horizontal gene transfer can also influence host fitness and community composition. However, there are very few studies of marine viral diversity in the Southern Hemisphere, which hampers our ability to fully understand the complex interplay of biotic and abiotic factors that shape microbial communities. We carried out the first genetic study of bacteriophage communities within a dynamic western boundary current (WBC) system, the east Australian current (EAC). Virus DNA sequences were extracted from 63 assembled metagenomes and six metaviromes obtained from various depths at 24 different locations. More than 1700 bacteriophage genomic fragments (>9 kbps) were recovered from the assembled sequences. Bacteriophage diversity displayed distinct depth and regional patterns. There were clear differences in the bacteriophage populations associated with the EAC and Tasman Sea euphotic zones, at both the taxonomic and functional level. In contrast, bathypelagic phages were similar across the two oceanic regions. These data provide the first characterisation of viral diversity across a dynamic western boundary current, which is an emerging model for studying the response of microbial communities to climate change.

Highly heterogeneous temporal dynamics in the abundance and diversity of the emerging pathogens Arcobacter at an urban beach.

While the significance of Arcobacter in clinical settings grows, the ecological dynamics of potentially pathogenic Arcobacter in coastal marine environments remains unclear. In this study, we monitored the temporal dynamics of Arcobacter at an urban beach subject to significant stormwater input and wet weather sewer overflows (WWSO). Weekly monitoring of bacterial communities over 24 months using 16S rRNA amplicon sequencing revealed large, intermittent peaks in the relative abundance of Arcobacter. Quantitative PCR was subsequently employed to track absolute abundance of Arcobacter 23S rRNA gene copies, revealing peaks in abundance reaching up to 108 gene copies L-1, with these increases statistically correlated with stormwater and WWSO intrusion. Notably, peaks in Arcobacter abundance were poorly correlated with enterococci plate counts, and remained elevated for one week following heavy rainfall. Using oligotyping we discriminated single nucleotide variants (SNVs) within the Arcobacter population, revealing 10 distinct clusters of SNVs that we defined as Arcobacter "ecotypes", with each displaying distinct temporal dynamics. The most abundant ecotype during stormwater and modelled WWSO events displayed 16S rRNA sequence similarity to A. cryaerophilius, a species previously implicated in human illness. Our findings highlight the diverse environmental drivers of Arcobacter abundance within coastal settings and point to a potentially important, yet overlooked exposure risk of these potential pathogens to humans.

Harnessing long-read amplicon sequencing to uncover NRPS and Type I PKS gene sequence diversity in polar desert soils.

The severity of environmental conditions at Earth's frigid zones present attractive opportunities for microbial biomining due to their heightened potential as reservoirs for novel secondary metabolites. Arid soil microbiomes within the Antarctic and Arctic circles are remarkably rich in Actinobacteria and Proteobacteria, bacterial phyla known to be prolific producers of natural products. Yet the diversity of secondary metabolite genes within these cold, extreme environments remain largely unknown. Here, we employed amplicon sequencing using PacBio RS II, a third generation long-read platform, to survey over 200 soils spanning twelve east Antarctic and high Arctic sites for natural product-encoding genes, specifically targeting non-ribosomal peptides (NRPS) and Type I polyketides (PKS). NRPS-encoding genes were more widespread across the Antarctic, whereas PKS genes were only recoverable from a handful of sites. Many recovered sequences were deemed novel due to their low amino acid sequence similarity to known protein sequences, particularly throughout the east Antarctic sites. Phylogenetic analysis revealed that a high proportion were most similar to antifungal and biosurfactant-type clusters. Multivariate analysis showed that soil fertility factors of carbon, nitrogen and moisture displayed significant negative relationships with natural product gene richness. Our combined results suggest that secondary metabolite production is likely to play an important physiological component of survival for microorganisms inhabiting arid, nutrient-starved soils.

Wave 2 strains of atypical Vibrio cholerae El Tor caused the 2009-2011 cholera outbreak in Papua New Guinea.

Vibrio cholerae is the causative agent of cholera, a globally important human disease for at least 200 years. In 2009-2011, the first recorded cholera outbreak in Papua New Guinea (PNG) occurred. We conducted genetic and phenotypic characterization of 21 isolates of V. cholerae, with whole-genome sequencing conducted on 2 representative isolates. The PNG outbreak was caused by an atypical El Tor strain harbouring a tandem repeat of the CTX prophage on chromosome II. Whole-genome sequence data, prophage structural analysis and the absence of the SXT integrative conjugative element was indicative that the PNG isolates were most closely related to strains previously isolated in South-East and East Asia with affiliations to global wave 2 strains. This finding suggests that the cholera outbreak in PNG was caused by an exotic (non-endemic) strain of V. cholerae that originated in South-East Asia.

Oceanographic boundaries constrain microbial diversity gradients in the South Pacific Ocean.

Marine microbes along with microeukaryotes are key regulators of oceanic biogeochemical pathways. Here we present a high-resolution (every 0.5° of latitude) dataset describing microbial pro- and eukaryotic richness in the surface and just below the thermocline along a 7,000-km transect from 66°S at the Antarctic ice edge to the equator in the South Pacific Ocean. The transect, conducted in austral winter, covered key oceanographic features including crossing of the polar front (PF), the subtropical front (STF), and the equatorial upwelling region. Our data indicate that temperature does not determine patterns of marine microbial richness, complementing the global model data from Ladau et al. [Ladau J, et al. (2013) ISME J 7:1669-1677]. Rather, NH4+, nanophytoplankton, and primary productivity were the main drivers for archaeal and bacterial richness. Eukaryote richness was highest in the least-productive ocean region, the tropical oligotrophic province. We also observed a unique diversity pattern in the South Pacific Ocean: a regional increase in archaeal and bacterial diversity between 10°S and the equator. Rapoport's rule describes the tendency for the latitudinal ranges of species to increase with latitude. Our data showed that the mean latitudinal ranges of archaea and bacteria decreased with latitude. We show that permanent oceanographic features, such as the STF and the equatorial upwelling, can have a significant influence on both alpha-diversity and beta-diversity of pro- and eukaryotes.

Systematic, continental scale temporal monitoring of marine pelagic microbiota by the Australian Marine Microbial Biodiversity Initiative.

Sustained observations of microbial dynamics are rare, especially in southern hemisphere waters. The Australian Marine Microbial Biodiversity Initiative (AMMBI) provides methodologically standardized, continental scale, temporal phylogenetic amplicon sequencing data describing Bacteria, Archaea and microbial Eukarya assemblages. Sequence data is linked to extensive physical, biological and chemical oceanographic contextual information. Samples are collected monthly to seasonally from multiple depths at seven sites: Darwin Harbour (Northern Territory), Yongala (Queensland), North Stradbroke Island (Queensland), Port Hacking (New South Wales), Maria Island (Tasmania), Kangaroo Island (South Australia), Rottnest Island (Western Australia). These sites span ~30° of latitude and ~38° longitude, range from tropical to cold temperate zones, and are influenced by both local and globally significant oceanographic and climatic features. All sequence datasets are provided in both raw and processed fashion. Currently 952 samples are publically available for bacteria and archaea which include 88,951,761 bacterial (72,435 unique) and 70,463,079 archaeal (24,205 unique) 16 S rRNA v1-3 gene sequences, and 388 samples are available for eukaryotes which include 39,801,050 (78,463 unique) 18 S rRNA v4 gene sequences.

Unicellular Cyanobacteria Are Important Components of Phytoplankton Communities in Australia's Northern Oceanic Ecoregions.

The tropical marine environments of northern Australia encompasses a diverse range of geomorphological and oceanographic conditions and high levels of productivity and nitrogen fixation. However, efforts to characterize phytoplankton assemblages in these waters have been restricted to studies using microscopic and pigment analyses, leading to the current consensus that this region is dominated by large diatoms, dinoflagellates, and the marine cyanobacterium Trichodesmium. During an oceanographic transect from the Arafura Sea through the Torres Strait to the Coral Sea, we characterized prokaryotic and eukaryotic phytoplankton communities in surface waters using a combination of flow cytometry and Illumina based 16S and 18S ribosomal RNA amplicon sequencing. Similar to observations in other marine regions around Australian, phytoplankton assemblages throughout this entire region were rich in unicellular picocyanobacterial primary producers while picoeukaryotic phytoplankton formed a consistent, though smaller proportion of the photosynthetic biomass. Major taxonomic groups displayed distinct biogeographic patterns linked to oceanographic and nutrient conditions. Unicellular picocyanobacteria dominated in both flow cytometric abundance and carbon biomass, with members of the Synechococcus genus dominating in the shallower Arafura Sea and Torres Strait where chlorophyll a was relatively higher (averaging 0.4 ± 0.2 mg m-3), and Prochlorococcus dominating in the oligotrophic Coral Sea where chlorophyll a averaged 0.13 ± 0.07 mg m-3. Consistent with previous microscopic and pigment-based observations, we found from sequence analysis that a variety of diatoms (Bacillariophyceae) exhibited high relative abundance in the Arafura Sea and Torres Strait, while dinoflagellates (Dinophyceae) and prymnesiophytes (Prymnesiophyceae) were more abundant in the Coral Sea. Ordination analysis identified temperature, nutrient concentrations and water depth as key drivers of the region's assemblage composition. This is the first molecular and flow cytometric survey of the abundance and diversity of both prokaryotic and picoeukaryotic phytoplankton in this region, and points to the need to include the picocyanobacterial populations as an essential oceanic variable for sustained monitoring in order to better understand the health of these important coastal waters as global oceans change.

Diversity and Activity of Diazotrophs in Great Barrier Reef Surface Waters.

Discrepancies between bioavailable nitrogen (N) concentrations and phytoplankton growth rates in the oligotrophic waters of the Great Barrier Reef (GBR) suggest that undetermined N sources must play a significant role in supporting primary productivity. One such source could be biological dinitrogen (N2) fixation through the activity of "diazotrophic" bacterioplankton. Here, we investigated N2 fixation and diazotroph community composition over 10° S of latitude within GBR surface waters. Qualitative N2 fixation rates were found to be variable across the GBR but were relatively high in coastal, inner and outer GBR waters, reaching 68 nmol L-1 d-1. Diazotroph assemblages, identified by amplicon sequencing of the nifH gene, were dominated by the cyanobacterium Trichodesmium erythraeum, γ-proteobacteria from the Gamma A clade, and δ-proteobacterial phylotypes related to sulfate-reducing genera. However, diazotroph communities exhibited significant spatial heterogeneity, correlated with shifts in dissolved inorganic nutrient concentrations. Specifically, heterotrophic diazotrophs generally increased in relative abundance with increasing concentrations of phosphate and N, while Trichodesmium was proportionally more abundant when concentrations of these nutrients were low. This study provides the first in-depth characterization of diazotroph community composition and N2 fixation dynamics within the oligotrophic, N-limited surface waters of the GBR. Our observations highlight the need to re-evaluate N cycling dynamics within oligotrophic coral reef systems, to include diverse N2 fixing assemblages as a potentially significant source of dissolved N within the water column.

Circular linkages between soil biodiversity, fertility and plant productivity are limited to topsoil at the continental scale.

The current theoretical framework suggests that tripartite positive feedback relationships between soil biodiversity, fertility and plant productivity are universal. However, empirical evidence for these relationships at the continental scale and across different soil depths is lacking. We investigate the continental-scale relationships between the diversity of microbial and invertebrate-based soil food webs, fertility and above-ground plant productivity at 289 sites and two soil depths, that is 0-10 and 20-30 cm, across Australia. Soil biodiversity, fertility and plant productivity are strongly positively related in surface soils. Conversely, in the deeper soil layer, the relationships between soil biodiversity, fertility and plant productivity weaken considerably, probably as a result of a reduction in biodiversity and fertility with depth. Further modeling suggested that strong positive associations among soil biodiversity-fertility and fertility-plant productivity are limited to the upper soil layer (0-10 cm), after accounting for key factors, such as distance from the equator, altitude, climate and physicochemical soil properties. These findings highlight the importance of surface soil biodiversity for soil fertility, and suggest that any loss of surface soil could potentially break the links between soil biodiversity-fertility and/or fertility-plant productivity, which can negatively impact nutrient cycling and food production, upon which future generations depend.

Multiple stressors in sediments impact adjacent hard substrate habitats and across biological domains.

Coastal systems are increasingly impacted by human activities. While the direct effects of individual contaminants have been investigated, the potential for multiple contaminants to impact adjacent hard substrate habitats is poorly understood. Sediment-bound contaminants pose a risk to water column organisms through resuspension and the fluxing of dissolved nutrients and metals. This study experimentally manipulated contaminated coastal sediments in mesocosms with additions of a common fertiliser to investigate the impact on both bacterial biofilms and macrofouling communities on nearby hard substrates. Field mesocosms were deployed sub-tidally for two weeks in a fully crossed design with two levels of metal contamination (ambient or high) and three levels of organic enrichment (ambient, low and high). Developing biofilm and macrofaunal communities were collected on acetate settlement sheets above the mesocosm sediments and censused with a combination of high-throughput sequencing (biofilm) and microscopy (macrofauna). Organic enrichment of sediments induced compositional shifts in biofilm communities, reducing their diversity, evenness and richness. Furthermore, co-occurrence networks built from microbial assemblages exposed to contaminated sediments displayed reduced connectivity compared to controls, suggesting a more stochastic assembly dynamic, where microbial interactions are reduced. Macrofouling community composition shifted in response to increased enrichment with separate and interactive effects of metals also observed for individual taxa. Specifically, antagonistic stressor interactions were observed for colonial ascidians and arborescent bryozoans; metal contamination decreased abundances of these taxa, except under high enrichment conditions. Together these micro- and macrofaunal responses indicate selection for depauperate, but contaminant-tolerant, communities and a potential breakdown in biotic connectivity through multiple stressor impacts across habitat boundaries.