Hydrothermal activity lowers trophic diversity in Antarctic hydrothermal sediments

Abstract. Hydrothermal sediments are those in which hydrothermal fluid is discharged through sediments and are one of the least studied deep-sea ecosystems. We present a combination of microbial and biochemical data to assess trophodynamics between and within hydrothermal and background areas of the Bransfield Strait (1050–1647 m of depth). Microbial composition, biomass, and fatty acid signatures varied widely between and within hydrothermally active and background sites, providing evidence of diverse metabolic activity. Several species had different feeding strategies and trophic positions between hydrothermally active and inactive areas, and the stable isotope values of consumers were not consistent with feeding morphology. Niche area and the diversity of microbial fatty acids was lowest at the most hydrothermally active site, reflecting trends in species diversity. Faunal uptake of chemosynthetically produced organics was relatively limited but was detected at both hydrothermal and non-hydrothermal sites, potentially suggesting that hydrothermal activity can affect trophodynamics over a much wider area than previously thought.


Introduction
Hydrothermal sediments (also called sediment-hosted or sedimented hydrothermal vents), the product of subsurface mixing between hydrothermal fluid and ambient seawater within the sediment, are physically more similar to background deep-sea habitats than to high-temperature hard substratum vents (Bemis et al., 2012;Bernardino et al., 2012). This means that whilst they can host chemosynthetic obligate species, they can also be more easily colonised by nonspecialist fauna and potentially offer an important metabolic resource in the nutrient-limited deep sea (Levin et al., 2009;Dowell et al., 2016). Hydrothermal sediments have also been suggested to act as evolutionary bridges between hard substratum vents and methane seeps (Kiel, 2016). To utilise in situ production in hydrothermal sediments, fauna must overcome the environmental stress associated with hightemperature, acidic, and toxic conditions Gollner et al., 2015). The combination of elevated toxicity and in situ organic matter (OM) production results in a different complement of ecological niches between hydrothermal and background conditions that elicits compositional changes along a productivity-toxicity gradient (Bernardino et al., 2012;Gollner et al., 2015;Bell et al., 2016b). Hydrothermal sediments offer different relative abundances of chemosynthetic and photosynthetic organic matter, depending upon the supply of surface-derived primary productiv-ity and levels of hydrothermal activity (Tarasov et al., 2005). In shallow environments (< 200 m of depth) where the production of chemosynthetic and photosynthetic organic matter sources can co-occur, consumption may still favour photosynthetic OM over chemosynthetic OM as this does not require physiological adaptations to environmental toxicity (Kharlamenko et al., 1995;Tarasov et al., 2005;Sellanes et al., 2011). The limited data available concerning trophodynamics at deep-sea hydrothermal sediments in the Arctic indicate that diet composition can vary widely between species (Sweetman et al., 2013).
Hydrothermal sediments host diverse microbial communities (Teske et al., 2002;Kallmeyer and Boetius, 2004). Microbial communities are a vital intermediate between inorganic substrates and metazoan consumers, and thus their composition and isotopic signatures are of direct relevance to metazoan food webs. The heat flux associated with hydrothermal activity provides both benefits and constraints to microbial communities (Kallmeyer and Boetius, 2004;Teske et al., 2014) and accelerates the degradation of organic matter, giving rise to a wide variety of compounds including hydrocarbons and organic acids (Martens, 1990;Whiticar and Suess, 1990;Dowell et al., 2016). Microbial aggregations are commonly visible on the sediment surface in hydrothermal sediments (Levin et al., 2009;Sweetman et al., 2013;Dowell et al., 2016), but microbial activity also occurs throughout the underlying sediment, occupying a wide range of geochemical and thermal niches (reviewed by Teske et al., 2014). Sedimented chemosynthetic ecosystems may present several sources of organic matter to consumers (Bernardino et al., 2012;Sweetman et al., 2013;Yamanaka et al., 2015) and the diverse microbial assemblages can support a variety of reaction pathways, including methane oxidation, sulfide oxidation, sulfate reduction, and nitrogen fixation (Teske et al., 2002;Dekas et al., 2009;Jaeschke et al., 2014). Phospholipid fatty acid (PLFA) analysis can be used to describe recent microbial activity and δ 13 C signatures (Boschker and Middelburg, 2002;Yamanaka and Sakata, 2004;Colaço et al., 2007). Although it can be difficult to ascribe a PLFA to a specific microbial group or process, high relative abundances of certain PLFAs can be strongly indicative of chemoautotrophy (Yamanaka and Sakata, 2004;Colaço et al., 2007) and can support an understanding of microbial ecosystem function in hydrothermal sediments (e.g. in western Pacific vents; see Yamanaka and Sakata, 2004).
Stable isotope analysis (SIA) is a powerful tool to assess spatial and temporal patterns in faunal feeding behaviour and has been used to study trophodynamics and resource partitioning in other hydrothermal sediments, predominately in the Pacific (Fry et al., 1991;Levin et al., 2009;Portail et al., 2016). Stable isotopic analyses provide inferential measures of different synthesis pathways and can elucidate a wide range of autotrophic or feeding behaviours. Carbon and sulfur isotopes are used to delineate food sources and nitrogen to estimate trophic position. The signature of source isotope ratios (δ 13 C and δ 34 S) is influenced by the isotopic ratio of the chemical substrate and the fractionation associated with the metabolic process involved; thus, different fixation pathways can elicit different isotopic signatures, even when derived from a single source (e.g. DIC; Fry et al., 1991). Possible δ 13 C isotopic values of sources in the Bransfield Strait include ∼ −40 ‰ for thermogenic methane, ∼ −27 ‰ for suspended particulate matter, and ∼ −15 ‰ for ice algae (Whiticar and Suess, 1990;Mincks et al., 2008;Henley et al., 2012;Young et al., 2013). As an example, Siboglinum spp. can use a range of resources, including methane or dissolved organic matter (Southward et al., 1979;Schmaljohann et al., 1990;Thornhill et al., 2008;Rodrigues et al., 2013), making SIA an ideal way to examine resource utilisation in these settings (Levin et al., 2009;Soto, 2009). We also apply the concept of an isotopic niche (Layman et al., 2007) whereby species or community trophic activity is inferred from the distribution of stable isotopic data in two-or threedimensional isotope space.

Hypotheses
We used a combination of microbial diversity data based on sequencing and compound-specific isotopic analyses and bulk isotopic data from sediment, microbial, macrofaunal, and megafaunal samples to investigate resource utilisation, niche partitioning, and trophic structure at hydrothermal and background sites in the Bransfield Strait to test the following hypotheses: (1) chemosynthetic organic matter will be an important food source in hydrothermal sediments; (2) siboglinid species subsist upon chemosynthetically derived OM; (3) stable isotope signatures will reflect a priori functional designations defined by faunal morphology; and (4) fauna will have distinct niches between hydrothermal sites and background areas.
2 Materials and methods

Sites and sampling
Samples were collected during RRS James Cook cruise JC55 in the austral summer of 2011 (Tyler et al., 2011) from three raised edifices along the basin axis (Hook Ridge, the Three  Low activity (9 cm yr −1 ) Aquilina et al. (2013Aquilina et al. ( , 2014, Hook Ridge 2 (HR2) 1054 High activity (34 cm yr −1 ) Bell et al. (2016b) Sisters, and the Axe) and one off-axis site in the Bransfield Strait (1024-1311m of depth; Fig. 1; Table 1). We visited two sites of variable hydrothermal activity (Hook Ridge 1 and 2) and three sites where hydrothermal activity was not detected (Three Sisters, the Axe, and an off-axis site; Aquilina et al., 2013). Of the two hydrothermal sites, Hook Ridge 2 had higher fluid advection rates and pore fluid temperature but lower concentrations of sulfide and methane (Dählmann et al., 2001;Aquilina et al., 2013Aquilina et al., , 2014. Samples were collected with a Bowers and Connelly dampened megacorer (1024-1311 m of depth) and a single Agassiz trawl at Hook Ridge (1647 m of depth). With the exception of salps, all microbial and faunal samples presented here were from megacore deployments. For a detailed description of the megacore sampling programme and macrofaunal communities, see Bell et al. (2016b). Sampling consisted of 1-6 megacore deployments per site, with 2-5 cores pooled per deployment (Bell et al., 2016b). Cores were sliced into 0-5 and 5-10 cm partitions and macrofauna were retained on a 300 µm sieve. Residues were preserved in either 80 % ethanol or 10 % buffered formalin initially and then stored in 80 % ethanol after sorting (Bell et al., 2016b). Fauna were sorted according to species or morphospecies level (for annelid and bivalve taxa), family level (for peracarids), and higher levels for less abundant phyla (e.g. echiurans). Salps were collected using an Agassiz trawl and samples were immediately picked and frozen at −80 • C and subsequently freeze-dried.

Phospholipid fatty acids
Samples of 3-3.5 g of freeze-dried sediment from Hook Ridge 1 and 2, the off-vent site, and the Three Sisters were analysed at the James Hutton Institute (Aberdeen, UK) following the procedure detailed in Main et al. (2015), as summarised below. Samples were from the top 1 cm of sediment for all sites except Hook Ridge 2 where sediment was pooled from two core slices (0-2 cm) due to sample mass limitations. Lipids were extracted following a method adapted from Bligh (1959) using a single-phase mixture of chloroform : methanol : citrate buffer (1 : 2 : 0.8 v : v : v). Lipids were fractionated using 6 mL ISOLUTE SI SPE columns preconditioned with 5 mL of chloroform. Freeze-dried material was taken up in 400 µL of chloroform, vortex mixed twice, and allowed to pass through the column. Columns were washed in chloroform and acetone (eluates discarded) and finally 10 mL of methanol. Eluates were collected, allowed to evaporate under an N 2 atmosphere, and frozen (−20 • C). Fatty acids were derivatised with methanol and KOH to produce fatty acid methyl esters (FAMEs). Samples were taken up in 1 mL of a 1 : 1 (v : v) mixture of methanol and toluene; 1 mL of 0.2 M KOH (in methanol) was added with a known quantity of an internal standard (C19:0, nonadecanoic acid), vortex mixed, and incubated at 37 • C for 15 min. After cooling to room temperature, 2 mL of isohexane : chloroform (4 : 1 v : v), 0.3 mL of 1 M acetic acid, and 2 mL of deionised water was added to each vial. The solution was mixed and centrifuged, the organic phase was transferred to a new vial, and the remaining aqueous phase was mixed and centrifuged again to further extract the organic phase, which was combined with the previous. The organic phases were evaporated under an N 2 atmosphere and frozen at −20 • C.
Samples were taken up in isohexane to perform gas chromatography-combustion isotope ratio mass spectrometry (GC-C-IRMS). The quantity and δ 13 C values of individual FAMEs were determined using a GC Trace Ultra with a combustion column attached via a GC Combustion III to a Delta V Advantage isotope ratio mass spectrometer (Thermo Finnigan, Bremen). The δ 13 C VPDB values (‰) of each FAME were calculated with respect to a reference gas of CO 2 , traceable to IAEA reference material NBS 19 TS-Limestone. Measurement of the Indiana University reference material hexadecanoic acid methyl ester (certified δ 13 C VPDB −30.74 ± 0.01 ‰) gave a value of 30.91 ± 0.31 ‰ (mean ± SD, n = 51). Combined areas of all mass peaks (m/z 44, 45 and 46) following background correction were collected for each FAME. These areas relative to the internal C19:0 standard were used to quantify the 34 most abundant FAMEs and related to the FAs from which they are derived (Thornton et al., 2011).

Bulk stable isotopes
All bulk isotopic analyses were completed at the East Kilbride node of the Natural Environment Research Council Life Sciences Mass Spectrometry Facility. Specimens with carbonate structures (e.g. bivalves) were physically decarbonated and all specimens were rinsed in deionised water (e.g. to remove soluble precipitates such as sulfates) and cleaned of attached sediment before drying. Specimens were dried for at least 24 h at 50 • C and weighed (mg, correct to 3 dp) into tin capsules and stored in a desiccator whilst awaiting SIA. Samples were analysed with a continuous-flow isotope ratio mass spectrometer using a vario PYRO cube elemental analyser (Elementar) coupled with a Delta Plus XP isotope ratio mass spectrometer (Thermo Electron). Each of the runs of CN and CNS isotope analyses used laboratory standards (gelatin and two amino acid-gelatin mixtures) and the international standard USGS40 (glutamic acid). CNS measurements used the internal standards (MSAG2: (methanesulfonamide-gelatin and M1: methionine) and the international silver sulfide standards IAEA-S1, IAEA-S2, and IAEA-S3. All sample runs included samples of freezedried, powdered Antimora rostrata (ANR), an external reference material used in other studies of chemosynthetic ecosystems (Reid et al., 2013;Bell et al., 2016a) to monitor variation between runs and instruments (Supplement S1). Instrument precision (SD) for each isotope measured from ANR was 0.42, 0.33, and 0.54 ‰ for carbon, nitrogen, and sulfur, respectively. The reference samples were generally consistent except in one of the CNS runs, which showed unusual δ 15 N measurements (Supplement S1), so faunal δ 15 N measurements from this run were excluded as a precaution. Stable isotope ratios are all reported in delta (δ) per mil (‰) notation relative to international standards: VPDB (δ 13 C), air (δ 15 N), and VCDT (δ 34 S). Machine error relative to these standards ranged from 0.01-0.23 for δ 13 C, 0.01-0.13 for δ 15 N, and 0.13-3.04 for δ 34 S. One of the sulfur standards (Ag 2 S IAEA, Supplement Sect. S2) had a notable difference from the agreed measurements, suggesting either a compromised standard or poor instrument precision. This error was not observed in other standards or the reference material used, but given the uncertainty here only δ 34 S differences greater than 3 ‰ are considered as significant.
A combination of dual-isotope (δ 13 C and δ 15 N; 319 samples) and tri-isotope (δ 13 C, δ 15 N and δ 34 S; 83 samples) techniques was used to describe bulk isotopic signatures of 43 species of macrofauna (35 from non-hydrothermal sites, 19 from hydrothermal sites, and 11 from both) and 3 megafaunal taxa and sources of organic matter. Samples submitted for carbon and nitrogen (CN) analyses were pooled if necessary to achieve an optimal mass of 0.7 mg (±0.5 mg). Where possible, individual specimens were kept separate in order to preserve variance structure within populations, but in some cases low sample mass meant individuals had to be pooled (from individuals found in replicate deployments). Optimal mass for tri-isotope (CNS) measurements was 2.5 mg (±0.5 mg) and, as with CN analyses, specimens were preferentially submitted as individual samples or pooled if necessary. Samples of freeze-dried sediment from each site were also submitted for CNS analyses (untreated for NS and acidified with 6 M HCl for C). Acidification was carried out by repeated washing with acid and deionised water. Specimens were not acidified. A pilot study and the subsequent results presented here confirm that the δ 13 C measurements (0.1 and 1.0 M HCl) were within the untreated population range in both polychaetes and peracarids and that acidification did not notably or consistently reduce δ 13 C standard deviation (Table 2). In the absence of a large or consistent treatment effect, the low sample mass (particularly for CNS samples) was dedicated to increasing replication and preserving the integrity of δ 15 N and δ 34 S measurements instead of separating carbon and nitrogen-sulfur samples (Connolly and Schlacher, 2013).
Formalin and ethanol preservation effects can both influence the isotopic signature of a sample (Fanelli et al., 2010;Rennie et al., 2012). Taxa that had several samples of each preservation method from a single site (to minimise intraspecific differences) were examined to determine the extent of isotopic shifts associated with preservation effects. Carbon and nitrogen isotopic differences between ethanol-and formalin-preserved samples ranged from 0.1-1.4 ‰ and 0.4-2.0 ‰, respectively. Differences across all samples were not significant (paired t test, δ 13 C: t = 2.10, df = 3, p = 0.126; δ 15 N: t = 1.14, df = 3, p = 0.337). Given the unpredictable response of isotopic signatures to preservation effects (which also cannot be extricated from within-site, intraspecific variation), it was not possible to correct isotopic data (Bell et al., 2016a). This contributed an unavoidable but generally quite small source of error in these measurements.

Statistical analyses
All analyses were completed in the R statistical environment (R Core Team, 2013). CN stable isotopic measurements were divided into those from hydrothermal and non-hydrothermal sites, averaged by taxa, and used to construct a Euclidean distance matrix (Valls et al., 2014). A similarity profile routine (SIMPROF, 10 000 permutations, p = 0.05, Ward linkage) was applied to the distance matrix in the clustsig package (v1.0) (Clarke et al., 2008;Whitaker and Christmann, 2013) to detect significant structure. The resulting cluster assignations were compared to a priori feeding groups (Bell et al., 2016b) using a Spearman correlation test (with 9999 Monte Carlo resamplings) in the coin package (v1.0-24; Hothorn et al., 2015). Isotopic signatures of species sampled from both hydrothermal and non-hydrothermal sites were also compared with a one-way ANOVA with Tukey's HSD pairwise comparisons (following a Shapiro-Wilk normality test).
Mean faunal measurements of δ 13 C and δ 15 N were used to calculate Layman metrics for each site (Layman et al., 2007), sample-size-corrected standard elliptical area (SEAc), and Bayesian posterior draws (SEA.B, mean of 10 5 draws ±95 % credibility interval) in the SIAR package (v4.2) (Parnell et al., 2010;Jackson et al., 2011). Differences in SEA.B between sites were compared in mixSIAR. The value of p given is the proportion of ellipses from group A that were smaller in area than those from group B (e.g. if p = 0.02, then 2 % of posterior draws from group A were smaller than the group B mean) and is considered to be a semi-quantitative measure of difference in means .

Differences in microbial composition along a hydrothermal gradient
A total of 28 767, 35 490, and 47 870 sequences were obtained from the off-axis site and the hydrothermal sites Hook Ridge 1 and 2, respectively. Bacteria comprised almost the entirety of each sample, with Archaea being detected only in the Hook Ridge 2 sample (< 0.1 % of sequences; Fig. 2). Hook Ridge 1 was qualitatively more similar to the off-axis site than Hook Ridge 2. Both Hook Ridge 1 (hydrothermal) and the off-vent site were dominated by Proteobacteria (48 and 61 % of reads, respectively; Fig. 2), whereas Flavobacteriia dominated Hook Ridge 2 (43 %; 7-12 % elsewhere) Table 4. PLFA profiles from freeze-dried sediment (nM per g dry sediment). PLFA names relate to standard notation (i: iso; a: anti-iso; first number: number of carbon atoms in chain; ω =: double bond; Me: methyl group).
Bransfield off-vent Three Sisters with Proteobacteria accounting for a smaller percentage of sequences (36 %; Fig. 2). By sequence abundance, Flavobacteriia were the most clearly disparate group between Hook Ridge 2 and the other sites. Flavobacteriia were comprised of 73 genera at Hook Ridge 2, 60 genera at BOV, and 63 genera at HR1, of which 54 genera were shared between all sites. Hook Ridge 2 had 15 unique flavobacteriial genera, but these collectively accounted for just 0.9 % of reads, indicating that compositional differences were mainly driven by relative abundance rather than taxonomic richness. The most abundant genus from each site was Arenicella at BOV and HR1 (7.1 and 5.2 % of reads, respectively) and Aestuariicola at HR2 (6.9 % of reads; Table 3). The four most abundant genera at both BOV and HR1 were Arenicella (γ -Proteobacteria), Methylohalomonas (γ -Proteobacteria), Pasteuria (Bacilli), and Blastopirellula (Planctomycetacia), though not in the same order, and accounted for 17.2 and 16.0 % of reads, respectively. The four most abundant genera at HR2, accounting for 20.2 % of reads, were Aestuariicola, Lutimonas, Maritimimonas, and Winogradskyella (Flavobacteriia). The genera Arenicella and Pasteuria were the most relatively abundant across all sites (2.2-7.1 % and 1.7-5.0 % of reads, respectively; Table 3).  Table 4). PLFA profiles from each of the non-hydrothermal sites sampled (off-axis and the Three Sisters, 33 and 34 FAs, respectively) were quite similar (Table 4) and shared all but one compound (16:1ω11c, present only at the non-hydrothermal Three Sisters site). Fewer FAs were enumerated from Hook Ridge 1 and 2 (31 and 23, respectively), including 3 FAs not observed at the nonhydrothermal sites (br17:0, 10-Me-17:0, and 10-Me-18:0), which accounted for 0.5-1.2 % of the total at these sites. Polyunsaturated algal biomarkers (20:5ω3 and 22:6 ω3) were only detected at the non-hydrothermal site (0.83-1.57 % of total FA abundance). Hook Ridge 2 had the lowest number of FAs and the lowest total FA biomass of any site, though this was due in part to the fact that this sample had to be pooled from the top 2 cm of sediment (top 1 cm at other sites). Bacterial biomass was highest at Hook Ridge 1 and ranged from 85-535 mg C m −2 (Table 3).

Description of bulk isotopic signatures
Most faunal isotopic signatures were within a comparatively narrow range (δ 13 C: −30 to −20 ‰, δ 15 N: 5 to 15 ‰, and δ 34 S: 10 to 20 ‰) and more depleted isotopic signatures were usually attributable to siboglinid species ( Fig. 3; Bell et al., 2017b). The Siboglinum spp. (found at all nonhydrothermal sites) had mean δ 13 C and δ 15 N values of −41.4 and −8.9 ‰, respectively, and Sclerolinum contortum (predominately from Hook Ridge 1 but found at both hydrothermal sites) had values of −20.5 and −5.3 ‰, respectively. Some non-endosymbiont-bearing taxa (e.g. macrofaunal neotanaids from the off-axis site and megafaunal ophiuroids at Hook Ridge 2) also had notably depleted δ 15 N signatures (means −3.6 to 2.6 ‰, respectively; Fig. 3). Isotopic signatures of sediment organic matter were similar between hydrothermal and non-hydrothermal sites for δ 13 C and δ 15 N but δ 34 S was significantly greater at nonhydrothermal sites (p < 0.05, Table 5; Fig. 4). Variability was higher in hydrothermal sediments for all isotopic signatures. Faunal isotopic signatures for δ 13 C and δ 34 S ranged much more widely than sediment signatures and indicate that sediment organics were a mixture of two or more sources of organic matter. A few macrofaunal species had relatively heavy δ 13 C signatures that exceeded −20 ‰, suggesting either a heavy source of carbon or marine carbonate in residual exoskeletal tissue, particularly for peracarids (∼ 0 ‰). Samples of pelagic salps from Hook Ridge had mean values for δ 13 C of −27.4 ‰ (±0.9) and for δ 34 S of 21.5 ‰ (±0.8).

Comparing macrofaunal morphology and stable isotopic signatures
Isotopic data (mean of each species for δ 13 C, δ 15 N, and δ 34 S) were used to construct a Euclidean distance matrix and the resultant hierarchy was compared to classifications based upon morphology. Species were assigned to one of four clusters (SIMPROF, p = 0.05; Fig. S3). No significant correlation between a priori (based on morphology) and a posteriori cluster assignations (based on isotopic data) was detected (Spearman correlation test: Z = −1.34; N = 43; p = 0.18). Clusters were mainly discriminated based on δ 15 N values, and peracarids were the only taxa to be represented in all of the clusters, indicating relatively high trophic diversity. Several taxa found at both hydrothermal and nonhydrothermal sites were assigned to different clusters between sites. A total of 11 taxa were sampled from both hydrothermal and non-hydrothermal regions, of which 4 were assigned to different clusters at hydrothermal and nonhydrothermal sites. Neotanaids (Peracarida: Tanaidacea) had the greatest Euclidean distance between hydrothermal and non-hydrothermal samples (11.36), demonstrating clear differences in dietary composition (Fig. 5). All other species were separated by much smaller distances between regions (range: 0.24 to 2.69). Raw δ 13 C and δ 15 N values were also compared between hydrothermal and non-hydrothermal samples for each species (one-way ANOVA with Tukey HSD pairwise comparisons). Analysis of the raw data indicated that δ 13 C signatures were different for neotanaids only and δ 15 N signatures were different for neotanaids and an oligochaete species (Limnodriloides sp.; ANOVA, p < 0.01; Fig. 5).

Community-level trophic metrics
All site niches overlapped (mean 50, range 30-82 %) and the positions of ellipse centroids were broadly similar for all sites (Table 6; Fig. 6). Hydrothermal site ellipse areas were similar but significantly smaller than non-hydrothermal ellipses (SEA.B, n = 10 5 , p = < 0.05). There were no significant differences in ellipse area between any of the nonhydrothermal sites. Ranges in carbon sources (dCr) were higher for non-hydrothermal sites (Table 6), indicating a greater trophic diversity in background conditions. Nitrogen range (dNr; Table 6) was similar between hydrothermal and non-hydrothermal sites, suggesting a similar number of trophic levels within each assemblage. All site ellipses had broadly similar eccentricity (degree of extension along long axis) ranging from 0.85 to 0.97 (Table 6); however, theta (angle of long axis) differed between hydrothermal and non-hydrothermal sites (−1.43 to 1.55 at Hook Ridge, 0.67 to 0.86 at non-hydrothermal sites). The range in nitrogen sources was more influential at hydrothermal sites than Sclerolinum contortum, which had very low δ 15 N signatures but similar δ 13 C values when compared with non-endosymbiont- bearing taxa from the same sites. The strongly depleted δ 13 C measurements of Siboglinum spp. meant that ellipse theta was skewed more towards horizontal (closer to zero) for nonhydrothermal sites.

Microbial signatures of hydrothermal activity
Fatty acid profiles at the non-hydrothermal off-axis and Three Sisters sites indicated similar bacterial biomass. Bacterial biomass varied much more widely at Hook Ridge (Table 4). The Hook Ridge 2 sample is not directly comparable since it was sampled from sediment 0-2 cm b.s.f. (rather than 0-1 cm b.s.f. owing to sample mass availability). Organic carbon content, hydrogen sulfide flux, and taxonomic diversity were all lower at this site and may support the suggestion of a lower overall bacterial biomass (Aquilina et al., 2013;Bell et al., 2016b). The very high bacterial biomass at Hook Ridge 1 suggests a potentially very active bacterial community comparable to other hydrothermal sediments (Yamanaka and Sakata, 2004), but δ 13 C org was qualitatively similar to non-hydrothermal sites, implying that chemosynthetic activity was not the dominant source of organic carbon or that the isotopic signatures of the basal carbon source (e.g. DIC) and the fractionation associated with FA synthesis resulted in similar δ 13 C signatures.
A small number of the more abundant fatty acids had notable differences in relative abundance between hydrothermal and background sites (Table 4). For example, 16:1ω7 has been linked to sulfur-cycling pathways (Colaço et al., 2007) and comprised 14.0-15.2 % of abundance at nonhydrothermal sites and 20.0-23.5 % at hydrothermal sites. However, 18:1ω7, also a suggested PLFA linked to thiooxidation (McCaffrey et al., 1989;Colaço et al., 2007) occurred in lower abundance at hydrothermal sites (4.8-11.1 %) than non-hydrothermal sites (15.9-16.9 %) and was also abundant in deeper areas of the Antarctic shelf (Würzberg et al., 2011). Heavier carbon isotopic signatures (> −15 ‰) are generally associated with rTCA cycle carbon fixation (Hayes, 2001;Hugler and Sievert, 2011;Reid et al., 2013), suggesting that this pathway may have been active at the hydrothermal sites, albeit at probably quite low rates. Conversely, many of the lightest δ 13 C signatures (e.g. 19:1ω8, −56.6 ‰, off-axis site) were associated with the non-hydrothermal sites, although it should be noted that 19:1ω8 has not been definitively linked to a particular bacterial process (Koranda et al., 2013;Dong et al., 2015). Lower FA carbon isotope signatures with small ranges (e.g. Figure 5. Biplot of CN isotopic data from species sampled at both hydrothermal sites and non-hydrothermal background regions. Mean ± standard deviation, X-Y scales vary. −60 to −50 ‰) could also be indicative of methane cycling, but most FAs at all sites had δ 13 C of > −40 ‰. These results further suggest that chemosynthetic activity was relatively limited and support a rejection of hypothesis one; although there were differences between sites in PLFAs that are potentially indicative of chemosynthetic activity, these were not necessarily consistent between different PLFAs. The metabolic provenance of several of the more abundant PLFAs is also still uncertain. A number of fatty acids have been linked, though not exclusively, to chemoautotrophy, such as 10-Me-16:0 (Desulfobacter or Desulfocurvus, sulfate reducers) and 18:1ω7 (Yamanaka and Sakata, 2004;Colaço et al., 2007;Klouche et al., 2009;Boschker et al., 2014), and their presence may be consistent with the hydrothermal signature of the sediment microbial community. There were notable proportions of compounds nor-mally associated with sulfate-reducing bacteria (Kohring et al., 1994;Boschker et al., 2014). These included iC15:0, aiC15:0, 1C17:0, and aiC17:0, which together constituted ∼ 8-12 % of the FA suite. In addition, C16:1ω5c was relatively abundant (Fig. S1), and minor amounts of 10MeC16:0, C17:1ω8c, and cycloC17:0 were present. These have also been used as indicators of sulfate-reducing bacteria and sometimes of particular groups (e.g. Guezennec and Fiala-Medioni, 1996;Boschker et al., 2014). These compounds indicate the presence of sulfate-reducing bacteria, although perhaps not as the dominant group. Although the FA suite was indicative of active sulfur-cycling activity, it remains difficult to be conclusive about the origin of most FAs, even those which have been regularly observed in chemosynthetic contexts (e.g. 18:1ω7) and may still be abundant elsewhere (Würzberg et al., 2011). Together, C16:1ω7c and C18:1ω7 accounted for ∼ 25-35 % of the total FA suite and although they can be more generally associated with gram-negative eubacteria, they have frequently been linked to sulfur-oxidising bacteria in sediment samples (Pond et al., 1998;Yamanaka and Sakata, 2004;Boschker et al., 2014). Their dominance of the suite in the Bransfield Strait is similar to sediments from a vent in the Barbados Trench where C16:1ω7 and C18:1ω7 together contributed up to 50 % of FAs (Guezennec and Fiala-Medioni, 1996).
Long chain fatty acids (> C22) indicative of land plants (e.g. Yamanaka and Sakata, 2004) and typical indicators of marine phytoplankton production (e.g. C20:3ω5 and C22:6ω3) were very minor constituents, never accounting for more than 3 % of total PLFA mass and only detected at the non-hydrothermal sites: off-vent and Middle Sister. While their low abundance is at least partially accounted for by rapid degradation during sinking through the water column (Veuger et al., 2012), it also suggests that sedimentary FAs were predominantly of bacterial origin, whether due to bacterial reworking of photosynthetic organic matter or in situ production.
Chemotrophic bacterial sequences, such as Blastopirellula (Schlesner, 2015), and Rhodopirellula (Bondoso et al., 2014), were found at all sites in relatively high abundance, suggesting widespread and active chemosynthesis, though the lack of a particularly dominant bacterial group associated with chemosynthetic activity suggested that the supply of chemosynthetic OM was likely relatively limited. It re-mains difficult, however, to determine which FAs these bacterial lineages may be have been synthesising.
Some FAs also had marked differences in δ 13 C signatures, even where there was strong compositional similarity between sites (i.e. the non-hydrothermal sites). This suggested that either there were differences in the isotopic values of inorganic or organic matter sources or different bacterial metabolic pathways were active. Between the non-hydrothermal sites, these included PUFAs and MUFAs (polyunsaturated and monounsaturated fatty acids) such as 18:2ω6, 9 ( δ 13 C 24.4 ‰) and 19:1ω8 ( δ 13 C 19.1 ‰). Differences in PLFA δ 13 C between the hydrothermal sites also ranged widely, with the largest differences associated with PLFAs such as 16:1ω11t ( δ 13 C 17.2 ‰) and 10-Me-16:0 ( δ 13 C 11.0 ‰). However, it should be stressed that all PLFAs with larger δ 13 C differences between sites were comparatively rare and never individually exceeded 5 % of total abundance. Microbial signatures, whilst supporting the suggestion of chemosynthetic activity, are not indicative of chemosynthetic OM being the dominant source of organic matter to food webs at any site (hypothesis one). It is not possible to assess from PLFA data the relative importance of chemoautotrophic and photosynthetic OM sources, since PLFAs degrade quickly and therefore surface FA abundances are inevitably underestimated in deep water samples. Note: dSR reported only for Hook Ridge 1 and the off-vent site since δ 34 S values of siboglinids were only measured from these sites; hence dSr at other sites would be a considerable underestimate. As δ 34 S values were comparatively under-representative, these values were not used in the calculation of any other metric. Data rounded to 1 dp.

Siboglinids
Both species of infaunal siboglinid (Sclerolinum contortum from Hook Ridge and the Siboglinum sp. from the nonhydrothermal sites) appeared to subsist upon chemosynthetically derived organic matter, as evidenced by their morphology and also by their strongly 15 N-depleted isotopic signatures (see values with δ 15 N of < −2 ‰ in Fig. 3). Low δ 15 N signatures have also been observed in other siboglinids in a range of hydrothermal settings, such as Riftia pachyptila at the East Pacific Rise hard substratum vents (Rau, 1981). Diazotrophy has been detected previously in hydrothermal vents and cold seeps, typified by low δ 15 N values (e.g. Rau, 1981;Desai et al., 2013;Wu et al., 2014;Yamanaka et al., 2015). Diazotrophy in various reducing settings has been found associated with the anaerobic oxidation of methane (Dekas et al., 2009), methanotrophy (Mehta and Baross, 2006), and (in a non-marine cave) sulfate reduction (Desai et al., 2013). The latter is also consistent with the low δ 34 S signatures of both siboglinid species (Fig. 4), but gene expression analysis and/or isotopic tracing would be required to confirm this suggestion. Alternately, low δ 15 N signatures may be explained by the uptake of ammonium produced through dissimilatory nitrate reduction (Naraoka et al., 2008;Liao et al., 2014;Bennett et al., 2015) or strong isotopic fractionation during utilisation of ammonia (Naraoka et al., 2008;Liao et al., 2014;Bennett et al., 2015). Bulk faunal isotopic signatures are inadequate to determine which of these chemosynthesisrelated mechanisms is responsible for Siboglinum δ 15 N values, which would require analysis of the functional genes in the Siboglinum endosymbionts. The δ 15 N values for both siboglinids (δ 15 N Sclerolinum −5.3 ‰ ± 1.0, Siboglinum −8.9 ‰ ± 0.8) indicated reliance upon locally fixed N 2 (Rau, 1981;Dekas et al., 2009Dekas et al., , 2014Wu et al., 2014;Yamanaka et al., 2015) rather than utilisation of sediment organic nitrogen (δ 15 N = 5.7 ‰ ± 0.7). These values were also in contrast to the rest of the nonchemosynthetic obligate species, which generally had much heavier δ 15 N values. This supports hypothesis two that the siboglinid species were subsisting upon chemosynthetic OM, most likely supplied by their endosymbionts.
Carbon isotopic signatures in chemosynthetic primary production depend upon the mode of fixation and the initial 13 C of the inorganic substrate. Sclerolinum contortum δ 13 C J. B. Bell et al.: Hydrothermal activity lowers trophic diversity 5719 (−20.5 ‰ ± 1.0 ‰) was depleted in δ 13 C relative to Southern Ocean DIC by around 10 ‰ (Henley et al., 2012;Young et al., 2013), giving it a signal within the fractionation range of the reverse tricarboxylic acid cycle (Yorisue et al., 2012), but the concentration and isotopic composition of DIC can undergo considerable alteration in hydrothermal sediments (Walker et al., 2008). Therefore, without measurements of δ 13 C in pore fluid DIC, it was not possible to determine which fixation pathway(s) were being used by S. contortum endosymbionts.
Sulfur isotopic signatures in S. contortum were very low and quite variable (−26.7 ‰ ± 3.5 ‰). Sclerolinum endosymbionts may have been utilising sulfide either from hydrothermal fluid, microbial sulfate reduction, or redissolved from hydrothermal precipitates. Mineral sulfide was present at Hook Ridge that ranged between −28.1 and +5.1 ‰ (Petersen et al., 2004), which is consistent with the relatively high δ 34 S variability in S. contortum. Alternatively, sulfide supplied as a result of microbial sulfate reduction (Canfield, 2001) may have been the primary source of organic sulfur, similar to that of solemyid bivalves in reducing sediments (mean δ 34 S of −30 to −20 ‰; Vetter and Fry (1998) and in cold seep settings (Yamanaka et al., 2015). Sulfate reduction can also be associated with the anaerobic oxidation of methane (Whiticar and Suess, 1990;Canfield, 2001;Dowell et al., 2016), suggesting that methanotrophic pathways could also have been important at Hook Ridge. (e.g. abundance of Methylohalomonas, 2.1-4.3 % of sequences at all sites; Table 3). Although endosymbiont composition data were not available for the Southern Ocean population, Sclerolinum contortum is also known from hydrocarbon seeps in the Gulf of Mexico (Eichinger et al., 2013;Eichinger et al., 2014;Georgieva et al., 2015) and the Håkon Mosby mud volcano in the Arctic Ocean where S. contortum δ 13 C ranged between −48.3 and −34.9 ‰ (Gebruk et al., 2003), demonstrating that this species can occupy several reducing environments and use a range of chemosynthetic fixation pathways, including sulfide oxidation and methanotrophy (Eichinger et al., 2014;Georgieva et al., 2015).
The Siboglinum sp. δ 13 C values (mean −41.4, range −45.7 to −38.1 ‰, n = 8) corresponded very closely to published values of thermogenic methane (−43 to −38 ‰) from the Bransfield Strait (Whiticar and Suess, 1990), strongly suggesting that methanotrophy was the dominant carbon source for this species. Biogenic methane, although present in the Bransfield Strait, typically has much lower δ 13 C values (Whiticar, 1999;Yamanaka et al., 2015), indicating a hydrothermal and/or thermogenic source of methane in the Bransfield Strait (Whiticar and Suess, 1990). Sulfur isotopic signatures were also very low in the Siboglinum sp.: δ 34 S −22.9 ‰ in one sample from 15 pooled individuals from the off-axis site, the lowest measurement of δ 34 S reported for this genus (Schmaljohann and Flügel, 1987;Rodrigues et al., 2013). Methanotrophy in Siboglinum spp. has been previously documented at seeps in the NE Pacific (Bernardino and Smith, 2010), the Norwegian margin (δ 13 C −78.3 to −62.2 ‰) (Schmaljohann et al., 1990), and Atlantic mud volcanoes (δ 13 C range −49.8 to −33.0 ‰; Rodrigues et al., 2013). Rodrigues et al. (2013) also reported a greater range in δ 15 N than observed in the Bransfield siboglinids (δ 15 N −1.3 to 12.2 ‰ and −10.2 to −7.6 ‰, respectively). This suggests that in comparison to Siboglinum spp. in Atlantic mud volcanoes, which seemed to be using a mixture of organic matter sources (Rodrigues et al., 2013), the Bransfield specimens relied much more heavily upon a single OM source, suggesting considerable trophic plasticity in this genus worldwide.
Off-vent methanotrophy using thermogenic methane potentially illustrates an indirect dependence upon hydrothermalism (Whiticar and Suess, 1990). Sediment methane production is thought to be accelerated by the heat flux associated with the mixing of hydrothermal fluid in sediment (Whiticar and Suess, 1990), and sediment and Siboglinum isotopic data suggest that the footprint of hydrothermal influence may be much larger than previously recognised, giving rise to transitional environments (Bell et al., 2016a;Levin et al., 2016). A clear contribution of methane-derived carbon to consumer diets was limited predominately to neotanaids, which is consistent with the relatively small population sizes (64-159 ind. m 2 ) of the Siboglinum sp. observed in the Bransfield Strait (Bell et al., 2016b).

Organic matter sources
Pelagic salps, collected from an Agassiz trawl at Hook Ridge (1647 m), were presumed to most closely represent a diet of entirely surface-derived material and were more depleted in 13 C and more enriched in 34 S than were sediments (salp δ 13 C −27.4 ‰ and δ 34 S 20.1; Hook Ridge sediment δ 13 C −26.2 ‰ and δ 34 S 14.3). Salp carbon isotopic signatures were also lighter than the majority of macrofauna or sedimentary organic carbon, both at Hook Ridge and the nonhydrothermal sites (Fig. 3), and similar to other suspensionfeeding fauna in the Bransfield Strait (Elias-Piera et al., 2013).
Fauna with more depleted δ 34 S and/or more enriched δ 13 C values were likely to have derived at least a small amount of their diet from chemosynthetic sources (potentially indirectly through the non-selective consumption of detrital OM) both at hydrothermal and background regions (Bell et al., 2017a). Carbon and sulfur isotopic measurements indicated mixed sources for most consumers between chemosynthetic OM and surface-derived photosynthetic OM. The low content of algal biomarkers (particularly at the hydrothermal sites) suggests that phytodetritus was probably quite degraded and thus challenging to detect using short-lived fatty acids. However, the Bransfield Strait can be subject to substantial export production and it is probable that surface production contributes much more to sea floor OM than is evident from the fatty acid composition. Non-hydrothermal sediments were more enriched in 34 S than hydrothermal sediments, an offset that probably resulted from greater availability of lighter sulfur sources such as sulfide oxidation at Hook Ridge, even if surface-derived OM remained the dominant source of organic matter at the hydrothermal sites (Bell et al., 2017a).
Samples of bacterial mat could not be collected during JC55 (Tyler et al., 2011) and without these end-member measurements, it was not possible to quantitatively model resource partitioning in the Bransfield Strait using isotope mixing models (Phillips et al., 2014). Bacterial mats from hightemperature vents in the Southern Ocean had δ 34 S values of 0.8 ‰ (Reid et al., 2013), and at sedimented areas of the Loki's Castle hydrothermal vents in the Arctic Ocean they have δ 34 S values of −4.9 ‰ (bulk sediment; Jaeschke et al., 2014). Therefore it is probable that low faunal δ 34 S values represent a contribution of chemosynthetic OM (from either siboglinid tissue or free-living bacteria). Inorganic sulfur can also be a source to consumers when sulfide is utilised by free-living bacteria (δ 34 S ranged from −7.3 to 5.4 ‰; Erickson et al., 2009) and although we could not analyse the δ 34 S of fluid sulfide, sulfide crusts have been found at Hook Ridge and may provide a proxy for typical isotopic composition (δ 34 S −28.1 to 5.1 ‰; Petersen et al., 2004). There were several species (e.g. tubificid oligochaetes) that had moderately depleted δ 34 S signatures, such as Limnodriloides spp. (δ 34 S 7.6 ‰ at hydrothermal sites, −1.2 ‰ at nonhydrothermal sites; Fig. 4), further supporting the hypothesis of different trophic positions between hydrothermal and nonhydrothermal regions (hypothesis two). This provides evidence of the coupled anaerobic oxidation of methane-sulfate reduction, but overall the contribution of δ 34 S-depleted bacterial production did not seem widespread (further rejecting hypothesis four).
Without samples of all OM sources we cannot quantitatively assert that the faunal utilisation of chemosynthetic OM was low in the Bransfield Strait. Although isotopic data were consistent with several OM sources, it seemed unlikely that chemosynthetic OM was a dominant source of OM to the vast majority of taxa. The apparently limited consumption of chemosynthetic OM suggested that either it was not widely available (e.g. patchy or low density of endosymbiont-bearing fauna; Bell et al., 2016b) or that the ecological stress associated with feeding in areas of in situ production was a significant deterrent to many species (Bernardino et al., 2012;Levin et al., 2013).

A priori vs. a posteriori trophic groups
Classifications based upon morphology did not prove to be an accurate predictor of isotopic data, suggesting that faunal behaviour is potentially more important in determining dietary composition than morphology (e.g. having or lacking jaws). Peracarid species that possessed structures adapted to a motile, carnivorous lifestyle were assigned to a carnivore or scavenger guild (Bell et al., 2016b) and were distributed throughout the food web both at hydrothermal sites and background regions, indicating more diverse feeding strategies than expected. Taxa presumed to be deposit feeders (largely annelids) also had a large range of δ 15 N values. This may reflect the consumption of detritus from both "fresh" and more recycled or refractory OM sources as observed in other non-hydrothermal sedimented deep-sea habitats (Iken et al., 2001;Reid et al., 2012) or reflect variability in trophic discrimination related to diet quality (Adams and Sterner, 2000). A range of foraminifera have now been shown to utilise denitrification, which results in their having heavier δ 15 N values (Pina-Ochoa et al., 2010;Jeffreys et al., 2015). The result is high δ 15 N values in taxa without predatory morphology (e.g. oligochaetes). Tubificid oligochaetes had higher δ 15 N values at the hydrothermal sites, suggesting that they fed upon more recycled organic matter, possibly owing to greater microbial activity at hydrothermal sites.
Several taxa (e.g. ophiuroids at Hook Ridge) had low δ 15 N values relative to sediment OM, suggesting preferential consumption of chemosynthetic OM (Rau, 1981;Dekas et al., 2014). In these taxa, it is likely that the widespread but patchy bacterial mats or Sclerolinum populations at Hook Ridge (Aquilina et al., 2013) were an important source of organic matter. Fauna from the non-hydrothermal sites with low δ 15 N (e.g. neotanaids) were likely subsisting in part upon siboglinid tissue (Siboglinum spp.). There were no video transects over the off-axis site but footage of the Three Sisters, which was similar in macrofaunal composition (Bell et al., 2016b), did not reveal bacterial mats (Aquilina et al., 2013), and hence it is unlikely that these were an important resource at non-hydrothermal sites.
It is clear that some fauna can exhibit a degree of trophic plasticity depending upon habitat (supporting hypothesis four). This is consistent with other hydrothermal sediments in which several taxa (e.g. the Prionospio sp. Polychaeta: Spionidae) had different isotopic signatures depending upon their environment (Levin et al., 2009), demonstrating differential patterns in resource utilisation. Alternatively, there could have been different δ 15 N baselines between sites, though if these differences were significant, we argue that it is likely that more species would have had significant differences in tissue δ 15 N. Conversely, samples of Aurospio foodbancsia at both hydrothermal and non-hydrothermal sites had broadly similar δ 15 N values to those of the west Antarctic Peninsula (8 and 7.9 ‰, respectively), albeit with higher variability (Mincks et al., 2008). The δ 13 C values of Aurospio were also broadly similar, implying that this species occupied a detritivorous trophic niche regardless of environmental conditions.

Impact of hydrothermal activity on community trophodynamics
Standard ellipse area was lower at Hook Ridge than elsewhere (Table 6), which is analogous to trends in macrofaunal diversity and abundance in the Bransfield Strait (Bell et al., 2016b) and changes in SEA.B along a gradient of methane flux at vent and seep ecosystems in the Guaymas Basin (Portail et al., 2016). This demonstrates that at community level, ellipse area can be associated with other macrofaunal assemblage characteristics. A concurrent decline in niche area and alpha diversity is consistent with the concept that species have finely partitioned niches, and greater total niche area permits higher biodiversity (McClain and Schlacher, 2015). This relationship may also suggest that the influence of disturbance gradients created by hydrothermalism can result in an impoverished community (McClain and Schlacher, 2015;Bell et al., 2016b). Productivity-diversity relationships in which higher productivity sustains higher diversity have also been suggested for deep-sea ecosystems (McClain and Schlacher, 2015;Woolley et al., 2016), but this is not supported by the Bransfield Strait sites (Bell et al., 2017a). We suggest that in the Bransfield Strait the environmental toxicity in hydrothermal sediments (from differences in temperature and porewater chemistry) causes a concomitant decline in both trophic and species diversity (Bell et al., 2016b) in spite of the potential for increased localised production (Bell et al., 2017a). However, we acknowledge that, owing to the high small-scale habitat heterogeneity apparent from video imagery over the hydrothermally influenced area, it is likely that the contribution of chemosynthetic organic matter varies widely over tens of metres at Hook Ridge. Community-based trophic metrics (Layman et al., 2007) indicated that, although measures of dispersion within sites were relatively similar between hydrothermal sites and background areas (Table 6), trophic diversity, particularly in terms of the range of carbon sources (dCr) and total hull area (TA), were higher at background sites due to the more depleted carbon and nitrogen signatures of Siboglinum spp. It is still unclear whether the assemblage isotopic niche really corresponds to its actualised trophic niche, and although the niche space was smaller at the hydrothermal sites, the potential for different trophic strategies was still potentially greater (Bell et al., 2017a).

Conclusions
In this study, we demonstrate the influence of sedimenthosted hydrothermal activity upon trophodynamics and microbial populations. Low-activity hydrothermal microbiota were more similar to the non-hydrothermal site than to highactivity populations, illustrating the effect of ecological gradients upon deep-sea microbial diversity. Despite widespread bacterial mats and populations of vent-endemic macrofauna, the utilisation of chemosynthetic OM amongst non-specialist macrofauna and megafauna seemed relatively low, with a concomitant decline in trophic diversity with increasing hydrothermal activity. Morphology was also not indicative of trophic relationships, demonstrating the effects of differential resource availability and behaviour. We suggest that, be-cause these sedimented hydrothermal sites are insufficiently active to host large populations of vent-endemic megafauna, the transfer of chemosynthetic organic matter into the metazoan food web is likely to be more limited than in other similar environments.

Ethics statement
In accordance with the Antarctic Act (1994) and the Antarctic Regulations (1995), necessary permits (S5-4/2010) were acquired from the South Georgia and South Sandwich Islands Government.
The Supplement related to this article is available online at https://doi.org/10.5194/bg-14-5705-2017-supplement. Competing interests. The authors declare that they have no conflict of interest.