Next Article in Journal
Morphological and Molecular Characterization of Discolaimus haridwarensis sp. n. (Nematoda: Dorylaimida: Qudsianematidae) from India
Previous Article in Journal
A New Species of Cyprinid Genus Opsariichthys (Teleostei: Cyprinidae) from the Pearl River, Southern China
Previous Article in Special Issue
Decapods of the Southern Tip of South America and the Marine Protected Area Namuncurá–Burdwood Bank: A Nearshore–Offshore Comparison
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Biomarker Analysis Uncovers High Spatio-Temporal Stability of a Subarctic Rhodolith (Lithothamnion glaciale) Bed Food Web

by
Sean Hacker Teper
,
Christopher C. Parrish
and
Patrick Gagnon
*,†
Department of Ocean Sciences, Memorial University of Newfoundland, St. John’s, NL A1C 5S7, Canada
*
Author to whom correspondence should be addressed.
These authors have contributed equally to this work.
Diversity 2024, 16(10), 597; https://doi.org/10.3390/d16100597
Submission received: 7 April 2024 / Revised: 3 September 2024 / Accepted: 4 September 2024 / Published: 27 September 2024
(This article belongs to the Special Issue Marine Nearshore Biodiversity—2nd Edition)

Abstract

:
We used lipid, fatty acid, and stable isotope analyses to investigate variation, over nine months, in the trophodynamics of 10 dominant cryptofaunal, macroalgal/algal, and environmental components from two sites within a rhodolith (Lithothamnion glaciale) bed in southeastern Newfoundland (Canada). There was an overall shift from a diatom-based food web following the spring phytoplankton bloom to a kelp/algae-based food web during fall, accompanied by preferred use of EPA (20:5ω3) over DHA (22:6ω3) in most cryptofauna. The food web contained three trophic levels that encompassed: (1) direct feeding relationships from primary producers (e.g., rhodoliths, macroalgae) to second-order consumers (e.g., sea stars, polychaetes); (2) trophic subsidy from within and outside the rhodolith bed via settlement, resuspension, and consumption of macroalgal fragments and other detrital organic matter; and (3) strong pelagic/benthic coupling. Riverine input did not affect cryptofaunal diets, as shown by the lack of terrestrial biomarkers at the study site nearest to the riverine input, and there were minor differences in trophodynamics between both study sites. The present study’s findings, applicable to relatively broad spatial and temporal domains, as well as those of complementary studies of the same rhodolith bed, uncover high spatio-temporal stability of the rhodolith bed framework and of resident cryptofaunal abundance, diversity, and trophodynamics.

1. Introduction

Dietary linkages in marine benthic ecosystems remain largely unresolved because of the broad range of diets of many species and the sometimes complex benthic–pelagic relationships, which can further vary spatially and temporally [1,2,3]. Fluctuations in environmental conditions, including sea temperature [4], salinity [5,6,7], and chlorophyll (an indicator of food availability) [8], can trigger significant changes in marine invertebrate behavior, abundance, biomass, and feeding patterns [9,10,11,12]. Invertebrates can also adapt physiologically to environmental variability by modifying their diet [3], which ultimately affects their lipid class [13], fatty acid [14], and stable isotope composition [15]. Therefore, studying temporal variation in organismal lipids, fatty acids, and stable isotopes has the potential to more fully elucidate feeding relationships and their connection with environmental variability across benthic communities.
Rhodoliths (free-living, non-geniculate red coralline algae growing as balls, branched twigs, or rosettes) often form dense aggregations, known as rhodolith beds, at depths of up to 150 m in tropical to polar seas [16,17]. The relatively complex morphology of rhodoliths creates suitable habitats for attachment [18,19,20], reproduction [6,19,21], and feeding [5,6,7] of highly diverse algal and faunal assemblages. The important contribution of rhodolith beds to marine biodiversity [6,22,23] and global calcium carbonate (CaCO3) production [24,25,26] has contributed to the recent increase in studies on the factors and processes regulating their structure and function [27,28,29,30].
Knowledge about rhodolith ecosystem dynamics is largely inferred from a few studies on biodiversity patterns [31,32,33,34] and the pervasive effects of human disturbances on rhodoliths or rhodolith bed structure [22,23,34,35]. In general, studies broadly characterize changes in rhodolith bed macrofaunal assemblages without necessarily investigating the influences of environmental variability on such changes. However, because phytoplankton blooms, macroalgae, and bacterial components may contribute just as significantly to nutrition as the macroinvertebrates themselves [36], quantifying solely diversity and abundance may not offer a complete representation of the interconnections in rhodolith communities [3]. As environmental conditions shift seasonally in temperate seas, so may community dynamics, in part because of dietary shifts as resource availability changes [37]. Understanding the complexity of the trophodynamic system requires evaluating the effects of phytoplankton seasonality and seabed composition to support production for different trophic levels.
The present study characterizes spatial and temporal variation in trophic relationships of a large (>500 m2) rhodolith (Lithothamnion glaciale) bed in St. Philip’s, Newfoundland. It builds on previous studies that: (1) characterized rhodolith morphology and invertebrate diversity [6,38]; (2) addressed how the dominant invertebrate bioturbators prevent rhodolith burial [29]; (3) concluded that rhodoliths are resilient to changes in sea temperature [39,40]; and (4) identified bottom-up forcing and benthic–pelagic coupling as a mechanism and pathway for nutrient and energy flow [30] in the St. Phillip’s bed. The latter study was the first to focus on a robust trophodynamic analysis of a rhodolith bed, but it lacked spatial and temporal complexity. Specifically, we use chlorophyll levels from phytoplankton blooms to time three separate collections spanning nine months and analyze rhodoliths and their underlying substratum from two distant sites in the St. Philip’s rhodolith bed to characterize seasonal and spatial variation in lipid composition, nutrition, and trophodynamic relationships of dominant rhodolith cryptofauna. We test the hypotheses that: (1) seasonal fluctuations in sea temperature and food availability affect the lipid composition and diets of cryptofauna; and (2) riverine input affects cryptofaunal diets.

2. Materials and Methods

Unless specified otherwise in the following sections, we used the same field, analytical, and statistical approaches as Hacker Teper et al. [30]. Accordingly, we kept the methodological descriptions brief.

2.1. Study Site and Selection of Focal Species

Field work was completed during the spring, summer, and fall of 2017 in a rhodolith (Lithothamnion glaciale) bed which extends ~5 to 30 m in depth along the coast of St. Philip’s in southeastern Newfoundland, Canada. Consistent with our broader objective of characterizing spatial and temporal variability in rhodolith bed trophodynamics and differences in marine environmental conditions in this area, we chose to study two sites in the bed fringing Broad Cove, hereafter termed the “South” site (47°35′36.5″ N, 52°53′31.0″ W) and the “North” site (47°35′39.6″ N, 52° 53′ 24.6″ W). Both sites and their relative placement are shown in Hacker Teper et al. [30]’s Figure 1A. Adding the analysis of the “North” site to the present study allowed us to investigate trophic variability because of the expected greater environmental variability at this site resulting from its closer proximity to freshwater input from the marina.
We focused our food web analyses on the same six animal, two macroalgal/algal, and two environmental components of the bed: (1) common sea star, Asterias rubens; (2) wrinkled rock-borer, Hiatella arctica (bivalve); (3) juvenile Nereis spp. (polychaetes); (4) daisy brittle star, Ophiopholis aculeata; (5) green sea urchin, Strongylocentrotus droebachiensis; (6) red molted chiton, Tonicella marmorea, and Atlantic red chiton, T. rubra, which were too difficult to distinguish morphologically and, hence, were pooled into the genus Tonicella spp.; (7) pieces of Laminaria digitata [kelp] from nearby kelp beds; (8) Lithothamnion glaciale (rhodoliths); (9) seawater (containing seston) from a few meters above the rhodolith bed; and (10) sediment (containing infauna) underlying the rhodoliths. The present study addresses trophic interactions at both sites across three seasons (spring, summer, fall) and thereby expands knowledge of spatial and temporal dimensions in the Broad Cove rhodolith bed.

2.2. Timing of Sampling

To examine how benthic–pelagic coupling might change during different phytoplankton bloom events, we sampled the rhodolith bed: (1) during and at the end of the annual spring phytoplankton bloom in southeastern Newfoundland, when we anticipated peak diatom abundance in the water column [41,42]; and (2) at the end (early December 2017) of the annual fall phytoplankton bloom, when we expected peak dinoflagellate abundance [43]. The spring bloom began in the last few days of March and persisted through 23 April 2017, when we sampled the rhodolith community and collected rhodoliths for food web analyses. The bloom was declining when we sampled it on 10 July 2017 (see Appendix A for a detailed description of bloom characterization).

2.3. Rhodolith Community

To broadly characterize the rhodolith community, scuba divers hand-collected, on 23 April 2017, all the rhodoliths from one 30 × 30 cm quadrat placed every 5 m along a 30 m long transect at a depth of ~15 m at the South site and North site (for a total of seven quadrats sampled at each site). We chose the North site for its greater proximity to freshwater input and its higher sediment abundance around rhodoliths (personal observations) than the South site. Briefly, we transported the rhodoliths to the Ocean Sciences Centre (OSC) of Memorial University of Newfoundland (MUN) and placed them in large (320 L) holding tanks supplied with running seawater pumped in from the adjacent Logy Bay. We measured, using calipers, the lengths of the longest, intermediate, and short axes, as well as the gross weight of the 247 rhodoliths from the South site and 163 rhodoliths from the North site. We calculated each rhodolith’s sphericity according to Graham and Midgley [44] and Sneed and Folk [45] and as described by Gagnon et al. [6] and Bélanger and Gagnon [38]. To calculate rhodolith abundance (biomass) and cryptofaunal abundance, we weighed each rhodolith, then extracted, identified, weighed, and counted all cryptofauna from the external surface and inside each rhodolith [30].

2.4. Rhodolith Percent Cover and Epiphyte Coverage

To investigate rhodolith epiphyte coverage differences between sites, scuba divers took photographs (with a Canon PowerShot D30, St. John’s, Canada) of rhodoliths from one 30 × 30 cm quadrat placed every 5 m along each of four 15 m long transects at a depth of ~15 m at the South site and North site (for a total of 24 photo quadrats each). From each image, we estimated the epiphyte coverage on rhodoliths and the percent cover of rhodolith, sediment, and macroalgae (kelp and algae) using the point intercept method. A 7 × 7 grid was overlaid onto each image in Adobe Photoshop, and the element present underneath each grid intercept was recorded (i.e., rhodolith, sediment, or macroalgal). We marked rhodolith intercepts for the presence or absence of epiphytes. We measured epiphyte coverage only using intercepts overlaying rhodoliths. The 7 × 7 grid yielded 49 intercepts on each image, each representing ~2.1% cover. The outline of the 30 × 30 quadrat frame provided a reference scale. We used Microsoft Excel to compute various descriptive statistics (see Section 2.10.4).

2.5. Collection and Preparation of Samples for Food Web Analyses

To prevent contamination from our rhodolith and cryptofaunal abundance investigation (see Section 2.3), we collected additional rhodoliths to use with our food web analyses. At both the South and North sites on 23 April, 10 July, and 2 December 2017, divers hand-collected in sealed plastic bags ~150 live rhodoliths measuring 8 to 10 cm along the longest axis. They also collected seawater a few centimeters above the rhodolith bed with two 12 L Niskin bottles deployed gently to prevent resuspension of sediment from the bed, and collected three sediment samples scooped from the top (~10 cm) layer of muddy sediment underneath rhodoliths with 15 mL centrifuge tubes. Briefly, we transported all samples to the OSC within 4 h of collection, where we transferred rhodoliths in their sealed plastic bags to large (320 L) holding tanks supplied with running seawater and kelp and sediment in their sealed plastic bags or centrifuge tubes to freezers at −80 °C. Seawater was processed immediately upon arrival. After 24 h, we removed and processed rhodoliths in preparation for food web analyses. We used a hammer and screwdriver to break open rhodoliths, and tweezers and forceps to extract cryptofauna inside each rhodolith. We scrubbed kelp and rhodolith samples and used a mortar and pestle to grind rhodoliths into a powder. Seawater was processed through a mechanized filtration system, which suctioned water with an aspirator through a 47 mm diameter GF/C filter (Whatman; General Electric, WHA1822047, Oakville, Canada) at the bottom of a Büchner funnel. We obtained separate cryptofauna, kelp, sediment, and seawater specimens for our lipid and FA analyses versus our stable isotope analyses [30].

2.6. Extraction and Characterization of Lipid Classes

The extraction of lipids followed the protocols of Folch et al. [46], with modifications by Parrish [47]. Briefly, we separated samples into two aqueous layers though centrifugation and sonication and used a double pipetting technique to remove only the bottom, organic layer [30]. We then used thin-layer chromatography with flame ionization detection (TLC-FID) (Mitsubishi Kagaku Iatron, Inc., Mark VI Iatroscan, Tokyo, Japan) to characterize lipid classes, using a 3-step method to separate lipid classes [48].

2.7. Preparation and Characterization of Fatty Acid Methyl Esters (FAME)

Fatty acid methyl esters (FAME) of lipids were prepared directly by transesterification from aliquots of lipid extract following a modified procedure described by Christie [49] and Hamilton [50]. This transesterification process created an organic layer of FAME. We analyzed these FAME samples using gas chromatography and flame ionization detection (GC-FID) (HP 6890 gas chromatograph, Agilent Technologies Canada, Inc., Mississauga, Canada), identified fatty acids by retention times on the GC column, and integrated the fatty acid peaks using Varian Galaxie Chromatography Data System V1.9.3.2 software [30].

2.8. Stable Isotope Preparation and Analysis

Samples in their scintillation vials were prepared separately for carbon and nitrogen analyses, transferred to desiccators, and subsequently taken to The Earth Resources Research and Analysis (TERRA) facility at MUN for analysis. Briefly, we analyzed bulk stable isotopes and complementary C and N elemental proportions (%) with the following quantities of solidified samples from each vial: 1 to 1.5 mg for animals; 4 to 4.5 mg for rhodoliths and kelp; 14 to 15 mg for sediment; 5 to 7 mg of filter for seawater δ13C and carbon content; and 8 to 10 mg of filter for seawater δ15N and nitrogen content [30].

2.9. Trophic Magnification of Fatty Acids

Stable isotope ratios are expressed in the conventional (δ) notation as parts per thousand (‰) per the equation of Minagawa and Wada [51]:
δ13C or δ15N (‰) = [(Rsample/Rstandard) − 1] × 1000
where Rsample and Rstandard are the ratios of 13C/12C or 15N/14N of a given sample and the corresponding standard, respectively. Results are reported relative to atmospheric N2 for nitrogen stable isotopes, and Vienna Pee Dee Belemnite (VPDB) for carbon stable isotopes. We calculated species trophic position (TPconsumer) with the equation used by Gale et al. [52] and developed originally by Cabana and Rasmussen [53]:
TPconsumer = [δ15Nconsumer − δ15Nbase)/∆15N] + TPbase
where δ15Nconsumer is the mean stable N isotope ratio of each species and Δ15N is the fractionation factor, which, to be consistent with rhodolith food web studies, is 3.4‰ [54]. δ15Nbase and TPbase represent the nitrogen stable isotope composition and trophic positions from the base of the food web, respectively. We then calculated a trophic magnification factor (TMF) for fatty acids (FA) correlated with δ15N. This factor quantitatively represents the biomagnification of compounds along a food web [55,56]. TMF, the average change in a specific FA proportion or marker, was calculated from:
FA% = em×TP
or
loge FA% = (m × TP) + b
and, therefore,
TMF = em
where FA% is the fatty acid (as a % of total identified fatty acids) or fatty acid marker of interest; TP is the trophic position; and m and b are the slope and intercept of the linear equation, respectively.

2.10. Statistical Analysis

2.10.1. Lipid Classes

We used two separate two-way permutational MANOVAs (PERMANOVA) (Euclidean distance matrices with 9999 permutations) to test total lipid (N = 139) and lipid classes (N = 147) for temporal and spatial differences with the factors ‘month’ (April, July, and December) and ‘site’ (North and South). We also used separate one-way PERMANOVAs with the factor ‘component’ (each of the 10 components of the food web studied (six animal species, two macroalgal species, and two environmental components)) to examine differences in proportions of total lipids (and overall lipid classes among samples, accounting for accidental loss of samples during the analyses). To limit extraneous data variability while focusing on the most significant lipid classes, only lipid classes present in over 50% of the samples were included in the analysis. For consistency with the companion study [30], and for comparison purposes, we pooled the data into the following five functional groups reflecting the three dominant feeding strategies of the six animal species, two macroalgal species, and two environmental components: (1) suspension/filter feeders (two species); (2) grazers (two species); (3) predators (two species); (4) kelp and rhodolith; and (5) seawater/sediment (samples combined because of expected benthic–pelagic coupling) (see Section 2.1 for details). Using these five groups, we then ran a one-way PERMANOVA with the factor ‘functional group’. We ran one-way ANOVAs to test for relationships between individual components and total lipid or individual lipid classes. We examined relationships between total lipids and each of the major lipid classes with conventional Spearman rank-order correlation tests [57].

2.10.2. Fatty Acids

We used the same statistical approach (one two-way PERMANOVA with the factors ‘month’ and ‘site’ and two one-way PERMANOVAs, one with the factor ‘component’ followed by one with the factor ‘functional group’), with data exclusion to account for samples lost during analyses, to examine differences in the proportions of fatty acids and their proportional sums among samples (N = 142). We then used a one-way SIMPER analysis (run on untransformed data with a Bray–Curtis similarity matrix) with the factor ‘component’ (each of the 10 components of the food web studied (six animal species, two macroalgal species, and two environmental components)) to identify potential food sources and the main fatty acids contributing to the multivariate lipid structure of each component [3,32]. To limit extraneous data variability while focusing on the most significant fatty acids, only fatty acids contributing to over 70% of the similarities were included in the SIMPER analysis. We used a follow-up principal coordinates analysis (PCO; also run on untransformed data with a Bray–Curtis similarity matrix) with the factor ‘component’ (same as above), mainly for visualization of the feeding relationships among specific groups of organisms [58,59]. To increase clarity on the PCO, only samples with a Pearson coefficient of correlation >65% were displayed. We ran two-way ANOVAs to test if individual fatty acid composition among all components changed by month or by site. We also ran individual two-way ANOVAs to test if the proportions of the 8 fatty acids (16:0, 16:1ω7, 18:0, 18:1ω9, 18:1ω7, 20:1ω11, 20:4ω6 [ARA], and 20:5ω3 [EPA]) exhibiting at least 70% correlation changed significantly either seasonally or temporally according to specific components.

2.10.3. Stable Isotopes

We examined differences in carbon (δ13C) and nitrogen (δ15N) isotope ratios with two one-way PERMANOVAs (one for each type of ratio; both were types based on Euclidean distance matrices with 9999 permutations) with the factor ‘component’ and two two-way PERMANOVAs with the factors ‘month and site’. All taxa (six animal species, kelp, and rhodoliths) and both environmental components (seawater and sediment) were included in the δ13C isotope ratio analysis (N = 128, accounting for accidental loss of samples during the analyses). All components were also included in the δ15N isotope ratio analysis (N = 161), accounting for accidental loss of samples during the analyses. We also carried out a cluster analysis (with the “Group Average” clustering method) on δ13C and δ15N isotope ratios simultaneously, and a complementary SIMPROF test (Euclidian distance matrix with 9999 permutations) (N = 120 because of a few unmatched pairs of δ13C and δ15N ratios) to group and map, in the form of a dendrogram, statistically different components of the food web [32,54]. Seven main isotopic groups emerged from the SIMPROF test. We therefore ran a follow-up one-way PERMANOVA with both isotopic ratios combined and complementary one-way ANOVAs and Tukey HSD post hoc tests with the factor ‘group’ to examine differences among these seven main trophic groups.

2.10.4. General Aspects of Statistical Tests

In all PERMANOVAs, data were untransformed and computed on Bray–Curtis similarity or Euclidian distance matrices (9999 permutations). All fatty acid multivariate data were computed using Bray–Curtis similarity matrices, while lipid and stable isotope multivariate data were computed using Euclidean distances matrices due to their better ability to handle missing data. We used pairwise comparisons to further explore PERMANOVA tests with multiple measured variables. We used PERMDISP (9999 permutations) to inform our decision (p = 0.3236); we tested for homogeneity of multivariate variances and confirmed that all variances were homogenous. We used PCO (principal coordinate analysis) instead of PCA (principal component analysis) to more efficiently account for missing data [60]. We combined the North site and the South site in all figures because of the absence of statistical differences (t = 0.723, P (perm) = 0.695). Multivariate analyses were also conducted on fatty acid, isotope, lipid class, fluorescence, and temperature data together using Bray–Curtis similarity. When all the chemical and environmental variables were combined, lipids were the major drivers of significant differences (see Appendix B for additional justification of the general statistical approach). In all ANOVAs, we verified the homogeneity of variance and normality of residuals by examining the distribution of the residuals and the normal probability plot of the residuals, respectively [61]. We used a significance level of 0.05 in all analyses and reported all means with standard deviations (mean ± SD) unless stated otherwise. We used standard error where applicable to maintain consistency with corresponding literature [6,56,62]. We used PRIMER v7 with PERMANOVA+ for multivariate statistical analyses, Minitab 18 for univariate statistical analyses, and Microsoft Excel for descriptive statistics.

3. Results

3.1. Rhodolith Community

Rhodolith biomass was similar between the South (19.5 ± 0.1 (SE) kg m−2) and North (20.2 ± 0.2 kg m−2) sites. However, the significantly fewer rhodoliths (N = 163, p < 0.001) at the North site were significantly larger (114.1 ± 9.9 cm3, p < 0.001) than the 247 rhodoliths at the South site (70.0 ± 4.5 cm3). Rhodolith shape did not vary significantly between sites; together, they were predominantly spheroidal and compact (~62%), but otherwise, they were platy (~11%), bladed (~16%), or elongated (~6%) (Figure 1). Sediment (p = 0.264) and rhodolith (p = 0.074) percent cover did not vary significantly between sites, but epiphyte coverage per rhodolith was significantly higher (p < 0.001) at the North site (33.0 ± 3.0%) than the South site (14.1 ± 2.4%) (Figure 2). The South site’s macrofaunal biomass (34.5 ± 4.3 g kg−1 rhodoliths) nearly doubled that of the North site (18.8 ± 8.4 g kg−1 rhodoliths) (p = 0.024). In total, the 2918 animals extracted from the rhodoliths from the North (N = 1727) and South (N = 1191) sites did not differ significantly (p = 0.226) in terms of total abundance, but included at least 21 species representing six phyla. Echinoderms (528 ± 61 individuals kg−1 rhodoliths) and molluscs (516 ± 53 individuals kg−1 rhodoliths) numerically dominated both sites (Table 1). Species included in the biochemical analyses were particularly abundant at both sites, including Ophiopholis aculeata (359 ± 58 individuals kg−1 rhodoliths), Tonicella marmorea/T. rubra (226 ± 25 individuals kg−1 rhodoliths), and Hiatella arctica (161.9 ± 31.5 individuals kg−1 rhodoliths). A few species not included in the analyses were also relatively abundant, including the brittle star Ophiura robusta (133 ± 27 individuals kg−1 rhodoliths), the caridean shrimp Pandalus borealis (64 ± 12 individuals kg−1 rhodoliths), and the polychaete Potamilla reniformis (35.5 ± 8.8 individuals kg−1 rhodoliths).

3.2. Total Lipid Content and Lipid Classes

Of the ten food web components included in the lipid analysis, the brittle star O. aculeata contained the highest average concentrations of total lipids (9.8 ± 4.1 mg g−1 [ww]), ranging from 8.9 ± 1.6 mg g−1 at the North site in December to 15.0 ± 5.5 mg g−1 at the South site in July (Table S1). Both the macroalgal and environmental groups had significantly lower total lipid concentrations (0.4 ± 0.1 to 2.7 ± 3.0 mg g−1) than the animal group (t > 4.29, p < 0.01 in all cases). Total lipid varied significantly by component (PERMANOVA; Pseudo-F9,129 = 12.225, P (perm) < 0.01), but not between sites (Pseudo-F1,138 = 0.321, P (perm) = 0.55) nor among months (Pseudo-F2,138 = 0.813, P (perm) = 0.16). The interaction was also not significant.
The ten food web components contained nine lipid classes (PL, TAG, FFA, ST, AMPL, HC, SE, KET, and ALC, see Appendix C), with six (PL, TAG, FFA, ST, AMPL, and HC) present in >50% of all samples. PL was present and the dominant lipid class in every component, except for O. aculeata (North site, 35.9 ± 3.4%) and Tonicella spp. (South, 36.8 ± 4.5%; North, 33.4 ± 2.4%) in July, with an average proportional contribution to total lipid concentration of 44% in O. aculeata to 70% in seawater (Table S1). Proportions of PL of all combined components did not vary by site nor month; however, Tonicella spp. (ANOVA, F2.14 = 10.84, p < 0.01), L. digitata (F2.7 = 35.98, p < 0.01), and seawater (F2.16 = 22.80, p < 0.01) each significantly varied individually by month. High average proportions of TAG were most common in animal species, with 30% in O. aculeata, 29% in H. arctica, and 28.2% in Tonicella spp. TAG proportions varied significantly among some individual components; Tonicella spp. (F2,14 = 57.47, p < 0.01), Nerididae (F2,16 = 4.25, p = 0.038), and S. droebachiensis (F2,15 = 4.61, p = 0.033). ST was most prevalent in animal species, with average high proportions of 20% in S. droebachiensis and 18% in Nereis spp. ST proportions varied significantly among some individual components, such as Tonicella spp. (F2.14 = 8.26, p = 0.006) and seawater (F2,16 = 5.48, p = 0.019). PL to ST ratios in animal species significantly varied from April to December (pairwise PERMANOVA, t = 2.6363, P (perm) = 0.005) and nearly significantly from April to July (t = 1.8437, P (perm) = 0.055), but not between July and December (t = 0.9321, P (perm) = 0.373). (Table S1).
Overall lipid class composition varied significantly among components (PERMANOVA; Pseudo-F9,137 = 40.385, P (perm) < 0.01) and five functional groups (suspension/filter feeder, grazer, predator, macroalgal, and environmental components) studied (Pseudo-F4,142 = 19.363, P (perm) < 0.01), but did not vary monthly (Pseudo-F2,146 = 0.568, P (perm) = 0.68), nor between sites (Pseudo-F1,146 = 0.121, P (perm) = 0.94). The interaction was also not significant. Total lipid concentration significantly correlated with PL (r = 0.094, p = 0.036, N = 138), TAG (r = 0.409, p < 0.01, N = 138), FFA (r = −0.115, p < 0.01, N = 138), and ST (r = −0.095, p = 0.02, N = 138).

3.3. Fatty Acid Profiles

The ten food web components included in the fatty acid (FA) analysis contained 63 FAs, with 50 present in >50% of all samples. Each component exhibited a distinct set of dominant fatty acids and biomarkers (Table S2, Figure 3). With a proportional contribution to FA profiles ranging from 14% in H. arctica in December to 32% in A. rubens in April, EPA (eicosapentaenoic acid, 20:5ω3; a typical diatom biomarker) was typically the dominant FA within each of the six animal species sampled (Table S2). In December, ARA (arachidonic acid, 20:4ω6; a kelp and amphipod biomarker) and DHA (docosahexaenoic acid, 22:6ω3; a dinoflagellate biomarker) were, respectively, the most abundant FA in S. droebachiensis (~25%) and H. arctica (~20 to 24%). Across all months, ARA, vaccenic acid (18:1ω7; an aerobe, bacteria, and marine vegetation biomarker), palmitoleic acid (16:1ω7; a diatom biomarker), and palmitic acid (16:0; a common flagellate and marine vegetation biomarker [3,63,64,65]) were the next most prominent FA among the animals, with contributions peaking between 8 and 25%.
Macroalgal components were rich in palmitic acid (18 to 31%) and EPA (9 to 23%). Additionally, large proportions of oleic acid (18:1ω9, 8 to 18%, a common crustacean, detritus, dinoflagellate, and brown seaweed biomarker [3,66,67,68,69]) and linoleic acid (18:2ω6, 4 to 14%, a vegetation biomarker) were present in kelp, while 16:1ω7 (5 to 8%) was abundant in rhodoliths. Seston FAs were dominated by stearic acid (18:0, 18 to 33%; a common detritus biomarker [69,70]), 16:0 (15 to 34%), and 18:1ω9 (8 to 33%) (Table S2). Overall, FA composition differed significantly among the ten food web components (PERMANOVA; Pseudo-F9,132 = 53.596, P (perm) < 0.01) and five functional groups (Pseudo-F4,137 = 27.898, P (perm) < 0.01) studied.
FA composition changed significantly monthly (Pseudo-F2,141 = 2.09, P (perm) = 0.025), but a subsequent pairwise test revealed that it changed specifically only between April and December (PERMANOVA, t = 1.769, P (perm) = 0.0123). FAs at each site were not significantly different (PERMANOVA, Pseudo-F1,141 = 1.575, P (perm) = 0.14), nor was the interaction effect. In addition to the slight significant changes between months, nine individual FA profiles significantly differed seasonally (Table 2). Of the eight FA which exhibited at least 70% correlation in the ten components (Figure 3), seven (16:1ω7, 18:0, 18:1ω9, 18:1ω7, 20:1ω11, 20:5ω3 [EPA], and 20:4ω6 [ARA]) showed significant changes among individual components (Table 3).
Figure 3. (A) PCO plot (based on Bray–Curtis similarity matrices) of the 8 fatty acids exhibiting at least 70% correlation and DHA, which had 21% correlation in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites. (B) Typical fatty acid trophic biomarkers for fatty acids included in the analysis (adapted from Parrish [10] and Legeżyńska et al. [71]) (Dinos = Dinoflagellates).
Figure 3. (A) PCO plot (based on Bray–Curtis similarity matrices) of the 8 fatty acids exhibiting at least 70% correlation and DHA, which had 21% correlation in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites. (B) Typical fatty acid trophic biomarkers for fatty acids included in the analysis (adapted from Parrish [10] and Legeżyńska et al. [71]) (Dinos = Dinoflagellates).
Diversity 16 00597 g003
Of the three essential fatty acids (EPA, DHA, and ARA), EPA was the most abundant, being present in all food web components, but particularly among the six animal species (Table S2). ARA was present in all components except for trace amounts in seston, peaking in December in S. droebachiensis (25%) at higher proportions than EPA. All components had DHA, which was typically less abundant than EPA and ARA, except in H. arctica, which peaked at both sites in December with proportions higher than EPA (20–24%) (Table S2). Animal FA profiles were richest in polyunsaturated FA (PUFA), which ranged from 38% in Tonicella spp. in July to 67% in A. rubens in December (Table S3). Monounsaturated FAs (MUFAs) were lower, varying from 17% in S. droebachiensis in April to 35% in Tonicella spp. in July, but were nevertheless still generally more abundant than saturated FAs (SFAs), ranging from 8% in A. rubens in December to 27% in Nereis spp. in April.
Environmental components and L. glaciale presented lower levels of PUFA than the animals, ranging from 15% in seston in December to 39% in L. glaciale in July, and typically contained higher proportions of MUFA (35–43% in infauna) and SFA (17–65% in seston). The average DHA/EPA ratio was significantly higher in December than in April (pairwise PERMANOVA, t = 2.780, P (perm) = 0.01) and July (t = 2.785, P (perm) = 0.01) and was highest in H. arctica and seston (0.8 and 0.6, respectively) with both peaking in December at the South site. Conversely, the DHA/EPA ratio showed low values in sediment, L. digitata, Tonicella spp., A. rubens, S. droebachiensis, and L. glaciale (0.1 to 0.3), and were below detection (0) in Nereis spp. and O. aculeata (Table S3).
Proportional sums of FA varied significantly among the 10 components (PERMANOVA, Pseudo-F9,132 = 37.054, P (perm) < 0.01) and five functional groups studied (Pseudo-F4,137 = 45.993, P (perm) < 0.01), but not by site (Pseudo-F1,141 = 1.759, P (perm) = 0.18) nor month (Pseudo-F2,141 = 0.41, P (perm) = 0.77). The interaction was also not significant.

3.4. Stable Isotopes and Trophic Magnification

Stable carbon and nitrogen isotope ratios (δ13C and δ15N) showed no significant seasonal or temporal changes, except in kelp, which was more carbon-δ13-enriched in July and December (−15.2‰) than in April (−21.0‰) (δ13C, PERMANOVA, Pseudo-F2,6 = 332.97, P (perm) = 0.013) and increased its nitrogen-δ15 from 1 to ~1.98 during the same period (δ15N, Pseudo-F2,6 = 300.97, P (perm) = 0.005).
Stable carbon isotope ratios (δ13C) differed significantly among the 10 food web components included in the carbon isotope analysis (Pseudo-F9,118 = 50.929, P (perm) < 0.001), ranging from most depleted in seawater in December at the North site (−27.5‰) to most enriched in L. digitata in July at the South site (−15.1‰) (Table S4).
Stable nitrogen isotope ratios (δ15N), which were lowest in L. digitata in April (3.4‰) and highest in A. rubens in April and December (11.4‰) (Table S4), also differed significantly among the ten food web components included in the nitrogen isotope analysis (Pseudo-F9,151 = 135.15, P (perm) < 0.001), indicating distinct trophic levels (see below).
Hierarchical clustering analysis of δ13C and δ15N separated the ten components in seven distinct groups (Psuedo-F6,113 = 127.29, P (perm) < 0.001; Figure 4 and Table 4). Two of these groups each contained all samples of a single food web component, namely, L. digitata (kelp) and A. rubens (sea star), hereafter termed Group 1 and Group 6, respectively (Figure 4 and Table 4).
Group 2 also had only a single food web component, seawater (seston), but was separated into two sub-groups (Groups 2a and 2b).
Group 3 contained four subgroups, of which only one was monospecific, but each included at least one of Nereis (polychaete) (Groups 3a and 3c), seawater (seston) (Groups 3a, 3b, and 3c), or sediment (Groups 3b and 3d) (Figure 4 and Table 4).
Group 5 also included four subgroups, with each containing at least one of A. rubens (Groups 5a and 5c), Nereis spp. (Groups 5a and 5b), O. aculeata (Groups 5a and 5b), S. droebachiensis (Groups 5c and 5d), or Tonicella spp. (Groups 5a and 5d) (Figure 4 and Table 4).
Groups 4 and 7 each had two distinct subgroups of monospecific components: H. arctica (Group 4a), L. glaciale (Group 4b), Tonicella spp. (Group 7a), and L. digitata (Group 7b) (Figure 4 and Table 4).
Group 1 contained L. digitata, sampled in April, while Group 7b contained L. digitata, sampled in July and December. 20:5ω3 [EPA], and 22:4ω6 exhibited a TMF > 1, and hence were biomagnified through trophic levels, whereas 21 had a TMF < 1, indicating biodilution (Table 5).

4. Discussion

We used lipid, fatty acid, and stable isotope analyses to investigate the spatial and temporal variation in the trophodynamics of 10 animal and environmental components from two locations (sites) within a well-studied rhodolith (L. glaciale) bed in southeastern Newfoundland. In a previous study of the same rhodolith bed, we identified macroalgal-based detritus, but not kelp, as a key rhodolith community food source; identified three distinct trophic levels (producers, suspension/filter feeders and grazers, and predators); and suggested a direct link between an algae-based diet and carbon source in three cryptofaunal species (H. arctica, Tonicella spp., and O. aculeata) [30]. The present study’s broader spatial and temporal domains further strengthen the evidence for benthic–pelagic coupling, while unveiling generally high spatio-temporal stability in trophic linkages despite a species-specific ability to adjust diet based on food (including kelp) availability.

4.1. Rhodolith Community

Macroalgal biomass and species composition commonly shift seasonally in rhodolith beds because of changing environmental conditions such as water temperature and nutrient availability [17,22]. In the Gulf of California, foliose red algae peaked in the winter [22], but were seen only in the fall (September and October) in the present study. This difference could reflect greater light availability in the Newfoundland summer [72] and similar temperature conditions in both studies, where annual summer water temperature in the St. Philip’s bed reaches up to 17 °C [73] and the annual winter water temperature in the Gulf of California reaches up to 18 °C [74]. Interestingly, foliose red algae cover was one of the few components that differed between the North (33%) and South (14%) sites in the present study. Smith et al. [75] attributed the increased biomass of marine epiphytic algae in coastal marine systems to eutrophication. We could hypothesize that the potential influence of freshwater feeding directly into the North site is due to its greater proximity to riverine input. However, the low abundance of terrestrial fatty acid biomarkers (18:2ω6 and 18:3ω3) and the larger rhodoliths at the North (~114 cm3) compared to the South (~70 cm3) site suggests otherwise. Bélanger and Gagnon [72] showed that eutrophication inhibits rhodolith growth. Based on preliminary diving observations of sediment patches between rhodoliths, we assumed that rhodoliths at the North site were larger because of a slightly greater sedimentary load from the riverine source stabilizing the rhodoliths [6]. However, our statistical analysis showed no difference in sediment cover between the North (26%) and South (31%) sites. Moreover, the larger rhodoliths at the North site contained a lower macrofaunal biomass (~18 g kg−1 rhodoliths) than the smaller rhodoliths at the South site (~34 g kg−1 rhodoliths). It is possible, therefore, that the larger and hollow, or nucleated but structurally less complex, rhodoliths at the North site provided less space or substratum for cryptofaunal recruitment than the smaller rhodoliths at the South site. Further studies of epiphytic cover and cryptofaunal biomass closer to the freshwater input are needed to better explain this pattern.

4.2. Lipid Content and Class Composition

Lipid content fluctuates in organisms through seasonal variability associated with changing environmental conditions, food availability, metabolic states, or reproductive status [76,77,78]. In the present study, although individual lipid classes of some individual components varied seasonally (TAG in Nereis spp., S. droebachiensis, L. digitata; ST in O. aculeata; AMPL in S. droebachiensis; PL in L. digitata; and HC in seawater), a significant overall lipid content change occurred only in Tonicella spp. We detected significantly higher total lipid (~13 mg g−1) and lower PL levels (~35%) simultaneously with elevated levels of TAG (~45%) in July than in April or December. This change in lipid content suggests a seasonal diet shift to accommodate reproductive timing. A decrease in PL and an increase in short-term reserves (TAG) indicate that PL was expended and replaced with TAG to enable rapid turnover of the energy required for reproduction. Lee at al. [79] suggest that phospholipids have a storage function in cold-ocean invertebrates. Langer [80] observed peak gonad indices in April and May in Maine (USA), indicating reproductive timing of Tonicella marmorea, one of our sampled species. Based on the geographic difference between Maine and Newfoundland and lipid class changes in our study, the reproductive season of T. marmorea may begin later in Newfoundland.
In addition to organisms changing their lipid compositions in response to reproductive needs, we also expected to see lipid composition changes based on prevailing environmental conditions. In Hacker Teper et al. [30], we suggested that high proportions of phospholipids and unsaturated fatty acids, as well as low proportions of sterols, reflected the cold-water (~0.3 °C) conditions of the Northwest Atlantic in April [73]. The warmer waters of July (~6–8 °C) and December (~4 °C) could have resulted in the need for a less fluid lipid membrane through lower proportions of phospholipids and unsaturated fatty acids and higher proportions of sterols and saturated fatty acids [81]. We observed a decrease in phospholipid-to-sterol ratios from April (9.1) to July (6.3) to December (5.9), but not in polyunsaturated-to-saturated fatty acid ratios (2.6 vs. 2.2 vs. 2.7). This trend could suggest that the more rigid cell membrane lipid structure in July and December reflects the warming temperatures, whereas the unchanged fatty acid composition suggests strong benthic–pelagic coupling and a reliance on the polyunsaturated fatty acids available from major food sources (EPA from diatoms or DHA from dinoflagellates) [82].

4.3. Fatty Acids

As discussed in Hacker Teper et al. [30], our study sites presented high levels of EPA (~17%) and low levels of DHA (~3%) in April, probably because of the spring diatom phytoplankton bloom (~0.4 mg chlorophyll m−3). The present study showed similarly high levels of EPA (~18%) and low levels of DHA (~3%) at the end of the phytoplankton bloom in July (~0.07 mg m−3), likely because of the continuous sinking and accumulation of phytoplankton and organisms that utilize them on the seabed. In shallow waters, diatoms often reach the benthos intact following the spring bloom [81]. However, after the fall dinoflagellate phytoplankton bloom in December (nearly 0 mg m−3), we observed lower levels of EPA (~14%) and higher levels of DHA (~4%), which caused a significantly higher ratio of DHA compared to EPA (0.47 vs. 0.19 and 0.23). Combined with a lack of overall seasonal change between April and July, this change in diet suggests that organisms prioritize EPA and its availability over DHA. Although these trends were consistent across our entire sample set, a few species (H. arctica, O. aculeata, and A. rubens) exemplify this phytoplankton and fatty acid biomarker shift. For example, in H. arctica, the proportion of EPA decreased from April (~20%) and July (~22%) to December (~15%), whereas the proportion of DHA increased substantially from April (~10%) and July (~12%) to December (~22%). This change in a fatty acid biomarker provides clear evidence of a strong benthic–pelagic coupling link that resulted in filter- and suspension-feeding invertebrates utilizing the available phytoplankton; feeding is likely based on a “bloom” vs. “no bloom” scenario rather than on distinct seasons [83].
In Hacker Teper et al. [30], we suggested that A. rubens preys on H. arctica given their similar FA compositions and A. rubens’ feeding preference for mollusks [84]. Consistent levels of DHA in A. rubens from April through December compared to high levels of DHA observed in H. arctica support our hypothesis that H. arctica represents a major food source for A. rubens. Even though we cannot explain why there was little DHA in the diets of other invertebrates, the increase in the overall abundance of ARA from April (~6%) and July (~5%) to December (~8%) may explain their low levels of DHA. For example, in O. aculeata, the proportion of EPA decreased from April (~22%) and July (~26%) to December (~15%), whereas the proportion of ARA doubled from April (~1%) and July (~1%) to December (~2%). Similar increases in proportional ARA also occurred in Tonicella spp., Nereis spp., and S. droebachiensis. These increases in ARA rather than DHA suggest that, in the absence of abundant diatoms (EPA), organisms in the Newfoundland rhodolith community actively select an algal and kelp-derived (ARA) diet rather than feeding on abundant dinoflagellates (DHA). Kelp also becomes more abundant in the summer and early fall, providing an alternative food source for primary consumers [85].
Based on our findings and on Hacker Teper et al. [30], we hypothesize that a resource-partitioning relationship may exist where multiple species utilize different forms of the same food (i.e., algae/kelp and particulate algae/kelp) through various feeding strategies (i.e., suspension feeding and grazing), where the organisms benefit from the residual food from their co-habitants (Figure 5). Accordingly, because suspension-feeding O. aculeata and grazing Tonicella spp. likely cannot utilize large pieces of drift kelp, they could potentially utilize residual feeding matter in the form of particulate kelp after S. droebachiensis or other larger organisms have grazed on it and broken it down [86,87]. Simultaneously, the complex physical structure of rhodoliths might trap algal particulates that subsequently become more enriched through microbial activity, resulting in increased nutritive value through increased nitrogen content [88,89]. Nereis spp. may not require large amounts of DHA [90], and, given the similarity between its fatty acid profile and that of O. aculeata and Tonicella spp., it likely receives ARA by preying or scavenging on organisms in the rhodolith community rather than feeding directly on a kelp source. However, depending on the specific Nereis species and population, we propose at least three different foraging strategies. Nereis spp. could: (1) prey on rhodolith biota [91]; (2) consume detrital and microbenthic algae in the uppermost portion of the sediment layer on which the rhodoliths sit [92]; and (3) filter feed on phytoplankton [93]. Some populations may also switch from one mode to another [92]. In fact, N. diversicolor can change its feeding strategy from predatory to suspension feeding [93], whereas juvenile N. virens feed on detritus [94] and plant matter, which promotes growth [95]. Regardless, the diverse fatty acid assemblage of Nereis spp. suggests that small, likely juvenile, individuals feed on a combination of prey, detritus, and phytoplankton.
We expected (marked) differences in cryptofaunal diets between sites given the proximity of the North site to riverine input and our observed between-site differences in rhodolith size, cryptofaunal abundance, and rhodolith epiphytic foliose coverage. However, this was not the case, with only a few, and rather small, differences in trophodynamics between sites and a lack of terrestrial biomarkers at the North site. This may be the outcome of potentially low-quality terrestrial sources, i.e., refractory carbon [96], compared to a high availability of nutritious phytoplankton. The few terrestrial biomarkers in our seawater samples could also indicate limited mixing between riverine and marine systems or the destruction of terrestrial matter upon entering the ocean [97]. A more comprehensive sampling design that includes water upstream of, at, and downstream of the mouth of the river would be needed to better speak to the abundance, fate, and contribution of terrestrial matter to this rhodolith bed’s trophydynamics.

4.4. Stable Isotopes

We showed that even though consumers like A. rubens, H. arctica, and O. aculeata altered their diets following the spring phytoplankton bloom, their overall δ13C and δ15N isotopic values did not change spatially or temporally. The lack of seasonal change in δ13C suggests a two-tier dietary preference whereby the overall components of an organisms’ food source (i.e., phytoplankton or bacteria) remain consistent throughout the year, but the specific components of the food source (i.e., diatom vs. dinoflagellate) change. The temporal stability of the δ15N values supports the notion that the consumers sampled did not alter their diets extensively, but rather adjusted to the available resources.
Although we found no major temporal change in δ13C or δ15N values, there was still some variation at the species-individual level, likely reflecting food obtained during different stages of the recycling pathway outlined above. Because POM degrades as it settles on the sea floor, organisms feeding during later stages of the pathway (i.e., consuming resuspended material or predation) typically exhibit greater enrichment in δ13C and δ15N than organisms feeding on sinking POM [83]. Furthermore, overlap in δ15N values indicates reliance on similar food sources and extensive competition for food [83]. For example, in the present study, Nereis spp. split into four different bulk isotope subgroups. We believe this variation in isotopic value reflects opportunistic feeding, where Nereis spp. is not a major predator in this system, but rather an opportunist. Likewise, the isotopic values of A. rubens varied greatly. Most notably, two A. rubens samples had more enriched δ13C values (~14‰) than the majority of A. rubens samples (~19‰). Asterias rubens samples with ~14‰ δ13C values likely obtained their food from sources outside the recyclable macroalgae food-partitioning pathway (outlined in Section 4.3) and, therefore, probably do not fall within the benthic–pelagic coupling link we proposed for other focal community members. Instead, given their enriched values, we propose that these individuals preyed on organisms in a kelp- rather than phytoplankton-based food web [98]. However, L. digitata, a species with seasonal growth patterns [99], had the largest seasonal isotopic difference, with δ13C values jumping from ~−21‰ in April to ~−15‰ in July and December, and its δ15N increased from ~3‰ in April to ~7‰ in July and December. Low values in April compared to both our seasonal samples and typical averages (~−12‰ to ~−15‰ for δ13C and ~7‰ for δ15N) [100,101] could reflect a March sea ice event inhibiting L. digitata photosynthesis and growth. Indeed, prolonged sea ice limits light availability and thus may inhibit metabolic processes, resulting in stunted growth, reduced carbon uptake, less enriched δ13C values, and lower δ15N values [102].

5. Conclusions and Future Research Directions

Our results support our first hypothesis that seasonal fluctuations in sea temperature and food availability affect the lipid composition and diets of rhodolith cryptofauna. This was in part demonstrated by the overall shift from a diatom-based food web following the spring phytoplankton bloom to a kelp/algae-based food web during fall, and the observation that EPA is of higher value than DHA to the rhodolith bed community studied. Our fatty acid and stable isotope analyses strongly support Hacker Teper et al. [30]’s conclusion that first- and second-order consumers share a common resource (diatoms or kelp) in this rhodolith food web. The lack of terrestrial biomarkers at the study site nearest to riverine inputs, and minor differences in trophodynamics between both study sites, refuted our second hypothesis that riverine inputs affect cryptofaunal diets, although a more comprehensive water sampling design that includes the river/ocean transition zone would be required in order to fully elucidate this aspect.
Based on our findings, we propose a partial food web (Figure 5) encompassing: (1) direct feeding relationships from primary producers (e.g., rhodoliths and macroalgae) to second-order consumers (e.g., sea stars and polychaetes); (2) trophic subsidy from within and outside the rhodolith bed via settlement, resuspension, and consumption of macroalgal fragments and other detrital organic matter; and (3) benthic–pelagic coupling, as proposed in Hacker Teper et al. [30] and unequivocally demonstrated with the present study. This food web appears spatially and temporally stable in part because of the ability of key consumers to adjust their diets to available resources sharing similar nutritional characteristics (e.g., diatoms vs. dinoflagellates during or after phytoplankton blooms) [34,103]. The high and stable abundance of a few key consumers, namely, chitons and brittle stars ([6], present study), along with their dominant feeding modes [38], i.e., grazing and suspension feeding, respectively, may be key to this stability by greatly limiting epiphytic growth on rhodoliths. Epiphyte-free rhodolith surfaces should facilitate light capture and gas/nutrient uptake, contributing to photosynthesis, growth, and maintenance of rhodolith shape and structural integrity. Complementary studies of the same rhodolith bed showed that: (1) rhodolith morphology (shape, size, hollowness) impacts rhodolith cryptofaunal abundance and diversity [6,38]; (2) rhodolith growth is adversely affected by coastal eutrophication [72], yet is resilient to changes in sea temperature within the 2–10 °C range [39]; and (3) bioturbation by a few dominant invertebrates within and outside rhodoliths is sufficient to prevent sediment accumulation on, and burial of, rhodoliths [29]. These and the present study, therefore, highlight the high spatio-temporal stability of the rhodolith bed framework and resident cryptofaunal abundance, diversity, and trophodynamics. Given a rapidly changing ocean climate, further laboratory and field studies should investigate the impacts and relative importance of physical, chemical, and biological forcing on rhodolith bed homeostasis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d16100597/s1. Table S1: Sample size (N), mean wet weight, mean total lipid, and mean proportion (%) of the six dominant lipid classes (PL: phospholipid; TAG: triacylglycerol; FFA: free fatty acid; ST: sterol; AMPL: acetone mobile polar lipid; and HC: hydrocarbon) in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites (see Table 1 for the entire list of animal species found). Each component group (animal, macroalgal, environmental) variable’s lowest and highest values are bolded.; Table S2: Mean proportion (%) of each of the 50 dominant fatty acids found in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites. Each component’s highest FA proportion is bolded. FA are listed in ascending order of retention time from the 30 m long ZB WAXplus (Phenomenex) GC column in the Varian Galaxie Chromatography Data System (see Hacker Teper et al. [30] for details).; Table S3: Sample size (N), mean proportional sum (Σ) of saturated (SFA), monounsaturated (MUFA), polyunsaturated (PUFA), ω3 (omega-3), and ω6 (omega-6) fatty acids, and mean ratios of polyunsaturated:saturated (P/S) and DHA:EPA (DHA/EPA), in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites. Each component group (animal, macroalgal, environmental) variable’s lowest and highest values are bolded.; Table S4: Sample size (N), bulk stable isotope ratio (δ13C and δ15N; ‰), and relative trophic position (TP) in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites. Trophic position is based on an isotopic model with a Δ15N fractionation factor of 3.4‰ (see Hacker Teper et al. [30]). Each component group (animal, macroalgal, environmental) variable’s lowest and highest values are bolded.

Author Contributions

Conceptualization, S.H.T., C.C.P. and P.G.; methodology, S.H.T., C.C.P. and P.G.; validation, S.H.T.; formal analysis, S.H.T.; investigation, S.H.T., C.C.P. and P.G.; field work, S.H.T. and P.G.; data curation, S.H.T.; resources, C.C.P. and P.G.; writing—original draft preparation, S.H.T., C.C.P. and P.G.; writing—review and editing, C.C.P. and P.G.; visualization, S.H.T.; supervision, C.C.P. and P.G.; project administration, C.C.P. and P.G.; funding acquisition, C.C.P. and P.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Sciences and Engineering Research Council (NSERC Discovery Grant, grant number RGPIN-2015-05343), Canada Foundation for Innovation (CFI Leaders Opportunity Funds, grant number 16940), and Research and Development Corporation of Newfoundland and Labrador (Ignite R&D, grant number 5003.070.002) grants to P.G., and by a Natural Sciences and Engineering Research Council (NSERC Discovery Grant, grant number RGPIN-2017-04639) grant to C.C.P.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the Supplementary Materials.

Acknowledgments

We are grateful to J. Wells and A. Pye for assistance with laboratory work, and to D. Bélanger, A. P. St-Pierre, S. Trueman, and the Ocean Sciences Centre Field Services Unit for assistance with fieldwork. We thank Memorial University’s CREAIT Network for their support with laboratory work. P. Snelgrove generously provided access to CTD equipment. We also thank P. Snelgrove and S. Dufour for constructive comments that helped to improve the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. Physicochemical Environment and Timing of Phytoplankton Bloom at the Study Site

The timing of the phytoplankton bloom was confirmed by casting, at the two study sites, a CTD (conductivity, temperature, and depth) instrument (SBE-19 PLUS; Seabird, Bellevue, WA, USA) equipped with PAR (QSP2300; Biospherical Instruments Inc., San Diego, CA, USA) and fluorescence (FLRT; WETstar fluorometer) sensors every few days from 1 to 29 March, 23 to 26 April, 30 May to 10 July, 12 to 26 October, and 2 December, 2017. A sheet of thick ice covered the sea surface from 30 March to 22 April, preventing CTD casts during this period. Logistical constraints prevented data collection between 27 October and 1 December, and malfunctioning of fluorescence sensor yielded inaccurate data during the fall collection up to 2 December. For each cast we lowered the instrument at a speed of ~1 m s−1 from the ocean surface down to just above the surface of the rhodolith bed. CTD data were plotted with Ocean Data View V4.0 (https://odv.awi.de/, accessed on 4 September 2024) and SigmaPlot V11.0. Fluorescence data within the first 3 m above the rhodolith bed (i.e., at depths between 12 and 15 m) were used to monitor the progression, and confirm the occurrence of, the bloom. Below we present fluorescence (Figure A1 and Figure A2), temperature (Figure A3, Figure A4 and Figure A5), and salinity (Figure A6, Figure A7 and Figure A8) profiles. We also present mean fluorescence, temperature, and salinity from these plots (Table A1).
Figure A1. Fluorescence between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from March to April 2017. The horizontal dashed line indicates the depth (12 m) below which fluorescence was considered in determining the timing of the phytoplankton bloom.
Figure A1. Fluorescence between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from March to April 2017. The horizontal dashed line indicates the depth (12 m) below which fluorescence was considered in determining the timing of the phytoplankton bloom.
Diversity 16 00597 g0a1
Figure A2. Fluorescence between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from May to July 2017. The horizontal dashed line indicates the depth (12 m) below which fluorescence was considered in determining the timing of the phytoplankton bloom.
Figure A2. Fluorescence between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from May to July 2017. The horizontal dashed line indicates the depth (12 m) below which fluorescence was considered in determining the timing of the phytoplankton bloom.
Diversity 16 00597 g0a2
Figure A3. Temperature between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from March to April 2017.
Figure A3. Temperature between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from March to April 2017.
Diversity 16 00597 g0a3
Figure A4. Temperature between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from May to July 2017.
Figure A4. Temperature between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from May to July 2017.
Diversity 16 00597 g0a4
Figure A5. Temperature between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from October to December 2017.
Figure A5. Temperature between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from October to December 2017.
Diversity 16 00597 g0a5
Figure A6. Salinity between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from March to April 2017.
Figure A6. Salinity between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from March to April 2017.
Diversity 16 00597 g0a6
Figure A7. Salinity between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from March to July 2017.
Figure A7. Salinity between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from March to July 2017.
Diversity 16 00597 g0a7
Figure A8. Salinity between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from October to December 2017.
Figure A8. Salinity between the sea surface and rhodolith bed at 15 m depth at the (A) South and (B) North sites from October to December 2017.
Diversity 16 00597 g0a8
Table A1. Mean temperature, fluorescence, and salinity readings within the first 3 m above the rhodolith bed (as averaged from data shown in Figure A1, Figure A2, Figure A3, Figure A4, Figure A5, Figure A6, Figure A7 and Figure A8).
Table A1. Mean temperature, fluorescence, and salinity readings within the first 3 m above the rhodolith bed (as averaged from data shown in Figure A1, Figure A2, Figure A3, Figure A4, Figure A5, Figure A6, Figure A7 and Figure A8).
TemperatureFluorescenceSalinity
SiteCollection Month (2017)°C (±SE)(mg/m3) (±SE)psu (±SE)
SouthApril0.3 (0.1)0.4 (0.1)31.4 (0.1)
July7.9 (0.5)0.1(0.1)31.2 (0.1)
December4.4 (0.1)0 (0.1)31.5 (0.1)
NorthApril0.3 (0.1)0.4 (0.1)31.4 (0.1)
July6.2 (0.2)0.1 (0.1)31.4 (0.1)
December0 (0.1)0 (0)31.5 (0.1)
Only fluorescence data within 3 m above the rhodolith bed were used in determining the timing of the bloom. This layer of water was deemed sufficiently narrow to capture benthic-pelagic trophic interactions relevant to the present study.
Inspection of the fluorescence data from March and April indicated that the spring phytoplankton bloom began in the last few days of March and was still ongoing on 23 and 26 April (Figure A1). Fluorescence above 12 m was more variable than deeper, with epibenthic levels gradually increasing throughout March and stabilizing in late April. Phytoplankton concentration on 23 April (when we sampled the rhodolith community and collected rhodoliths for food web analyses) was three to four times that on 8 March. Sea temperature varied more within the first 12 m than deeper, although the largest difference between coldest (14 March) and warmest (26 April) water across the entire water column did not exceed ~1.3 °C (Figure A3). Salinity varied by less than ~0.4 PSU across the entire water column from 1 to 20 March, but exhibited marked changes of up to 2 PSU within the first 5 m in the week that preceded the formation, and day that followed the retreat, of surface ice (Figure A6). Salinity below 12 m remained comparatively much more stable, with no obvious influence of sea ice and a faster return to pre-ice salinities.
Inspection of the fluorescence data from May to July indicated that the spring phytoplankton bloom dwindled from late May to early July (Figure A2). Fluorescence above 12 m was lower than deeper, indicating low photosynthesis at the surface. Epibenthic levels gradually decreased throughout June. Phytoplankton concentration on 10 July (when we sampled the rhodolith community) was about three times less than on 10 June. Sea temperature increased steadily from ~3 °C on 30 May to ~11 °C on 10 July, but with steep thermoclines on 6 and 10 July (Figure A4). The bottom temperature on day of final collection was over two times higher than at the beginning of the cast series. Salinity varied by less than ~0.8 PSU across the entire water column from 30 May to 10 July. Salinity below 12 m was comparatively less stable than in shallower water (Figure A7).
Inspection of the fluorescence data acquired on 2 December (when we sampled the rhodolith community) indicated that the fall phytoplankton bloom had terminated with ~0 mg phytoplankton m−3 above and below 12 m. Sea temperature on 2 December (~3 °C) was four times lower than on 12 October (~12 °C) (Figure A5), whereas salinity was ~1 psu higher on 2 December than on 12 October (Figure A8).

Appendix B. Justification for General Statistical Approach

Given our robust dataset, we only chose to report those statistics which were significant drivers of the studied rhodolith ecosystem. We ran separate analyses for fatty acids (Bray-Curtis similarity matrix), isotopes (Euclidian matrix), and lipid classes (Euclidian matrix), lipid classes and isotopes together (Euclidian matrix), and then fatty acids, isotopes, lipid classes, fluorescence, and temperature together (Bray-Curtis similarity matrix). We only used samples where all data were present and factored by “component” (components of the food web studied [six animal species, two macroalgal species, and two environmental components]), “site”, and “month”.
As reported in our study, species differences were preserved with almost all being highly significant, except Strongylocentrotus droebachiensis (PERMANOVA; Pseudo-F8,74 = 38.419, P (perm) < 0.01) (Table A2). We think S. droebachiensis was not highly significant because of limited sample size (N = 1 for this analysis). Temperature was within the top-ten contributors to dissimilarity in eight of the SIMPER comparisons, but it never ranked above six as a contributor. δ13C and δ15N were within the top eleven contributors to dissimilarity in five separate SIMPER comparisons, but they never ranked above eight as a contributor (Table A3).
Site showed no significant differences. Month, however, did show a small significant difference between April and July: P(perm) = 0.0339. It was driven mainly by lipid classes and EPA. Again, temperature was indeed a contributor, but was only the 9th most important one (Table A3).
When all variables were combined, lipids were the major drivers of overall significant differences, thus our main statistical focus throughout our study. Additional supporting SIMPER data is provided below:
Table A2. Summary of pair-wise analyses examining the statistical differences among component groups (the 8 components of the food web studied with no missing data [six animal species and two environmental components]). To use the Bray-Curtis similarity matrix with precision, we eliminated missing data (Laminaria digitata and Lithothamnion glaciale).
Table A2. Summary of pair-wise analyses examining the statistical differences among component groups (the 8 components of the food web studied with no missing data [six animal species and two environmental components]). To use the Bray-Curtis similarity matrix with precision, we eliminated missing data (Laminaria digitata and Lithothamnion glaciale).
Component GroupP-Perm
Asterias rubens vs. Seawater<0.001
Asterias rubens vs. Sediment<0.001
Hiatella arctica vs. Asterias rubens<0.001
Hiatella arctica vs. Nereididae<0.001
Hiatella arctica vs. Ophiopholis aculeata<0.001
Hiatella arctica vs. Strongylocentrotus droebachiensis0.0692
Hiatella arctica vs. Seawater<0.001
Hiatella arctica vs. Sediment<0.001
Nerididae vs. Asterias rubens<0.001
Nerididae vs. Strongylocentrotus droebachiensis0.1037
Nerididae vs. Seawater<0.001
Nerididae vs. Sediment<0.001
Ophiopholis aculeata vs. Asterias rubens<0.001
Ophiopholis aculeata vs. Nerididae<0.001
Ophiopholis aculeata vs. Strongylocentrotus droebachiensis0.2208
Ophiopholis aculeata vs. Seawater<0.001
Ophiopholis aculeata vs. Sediment<0.001
Strongylocentrotus droebachiensis vs. Asterias rubens0.0813
Strongylocentrotus droebachiensis vs. Seawater0.0606
Strongylocentrotus droebachiensis vs. Sediment0.2037
Seawater vs. Sediment<0.001
Tonicella spp. vs. Asterias rubens<0.001
Tonicella spp. vs. Hiatella arctica<0.001
Tonicella spp. vs. Nerididae<0.001
Tonicella spp. vs. Ophiopholis aculeata<0.001
Tonicella spp. vs. Strongylocentrotus droebachiensis0.0739
Tonicella spp. vs. Seawater<0.001
Tonicella spp. vs. Sediment<0.001
Table A3. Relative contribution to observed dissimilarity and mean proportion of all measured variables among samples for the three sampled months in pairwise comparisons (April–July, April–December, July–December). Only variables contributing to at least 70% similarity were included in the dissimilarity percentage analysis (SIMPER).
Table A3. Relative contribution to observed dissimilarity and mean proportion of all measured variables among samples for the three sampled months in pairwise comparisons (April–July, April–December, July–December). Only variables contributing to at least 70% similarity were included in the dissimilarity percentage analysis (SIMPER).
April–July
VariableDissimilarity (%)Average Abundance
Phospholipid (%)14.745.7
Free fatty acid (%)12.715.2
Triacylglycerol (%)9.914.9
20:5ω3 (EPA)5.217.7
18:1ω94.47.3
18:04.08.3
16:03.412.9
Acetone mobile polar lipid (%)3.48.7
Temperature3.30.3
Sterol (%)3.210.7
20:4ω6 (ARA)3.15.4
Total lipid3.06.6
April–December
VariableDissimilarity (%)Average Abundance
Phospholipid (%)14.145.7
Free fatty acid (%)12.015.2
Triacylglycerol (%)7.314.9
20:5ω3 (EPA)6.017.7
18:1ω94.67.3
16:04.112.9
Sterol (%)3.710.7
16:1ω73.77.2
18:03.48.3
Acetone mobile polar lipid (%)3.48.7
22:6ω3 (DHA)3.33.1
20:4ω6 (ARA)3.25.4
Total lipid2.66.6
July–December
VariableDissimilarity (%)Average Abundance
Phospholipid (%)14.240.1
Free fatty acid (%)14.120.6
Triacylglycerol (%)9.922.2
20:5ω3 (EPA)5.519.2
18:1ω94.46.8
Sterol (%)3.710.9
18:03.58.2
22:6ω3 (DHA)3.43.4
16:03.411.4
16:1ω73.25.1
Total lipid3.18.2
Acetone mobile polar lipid (%)2.73.1

Appendix C. List of Abbreviations and Symbols

AMPLAcetone mobile polar lipid(s)
ALCAlcohol(s)
ARAArachidonic acid
°CDegrees Celsius
cmCentimetres
CTDConductivity, Temperature, and Depth Instrument
DHADocosahexaenoic acid
DHA/EPADocosahexaenoic acid:eicosapentaenoic acid ratio
EPAEicosapentaenoic acid
EAElemental analyzer
FAFatty acid(s) [notation: A:Bωn (e.g.,. 20:5ω3)]
FAMEFatty acid methyl ester(s)
A:BωnFatty acid notation
FIDFlame ionization detection
FFAFree fatty acid(s)
gGrams
GCGas chromatography
GC-FIDGas chromatography and flame ionization detection
hHours
HCHydrocarbon(s)
KETKetone(s)
kgKilograms
LLitres
MSMass Spectrometer
MUNMemorial University of Newfoundland
mMetres
μLMicrolitres
mLMillilitres
mmMillimetres
minMinutes
MUFAMonounsaturated fatty acid(s)
ω3Omega-3 fatty acid(s)
ω6Omega-6 fatty acid(s)
OSCOcean Sciences Centre
PLPhospholipid(s)
PARPhotosynthetically Active Radiation
PUFAPolyunsaturated fatty acid(s)
P/SPolyunsaturated:saturated fatty acid ratio
PCOPrincipal coordinates analysis
rpmRevolutions per minute
SFASaturated fatty acid(s)
δ13CStable carbon isotope ratio
δ15NStable nitrogen isotope ratio
SDStandard deviation
SEStandard error
SE (lipid class)Steryl ether
STSterol(s)
TERRAThe Earth Resources Research and Analysis facility
TCDThermal conductivity meter
TLC-FIDThin-layer chromatography with flame ionization detection
TAGTriacylglycerol(s)
TMFTrophic multiplication factor
TPTrophic position

References

  1. Kharlamenko, V.I.; Kiyashko, S.I.; Imbs, A.B.; Vyshkvartzev, D.I. Identification of food sources of invertebrates from the seagrass Zostera marina community using carbon and sulfur stable isotope ratio and fatty acid analyses. Mar. Ecol. Prog. Ser. 2001, 220, 103–117. [Google Scholar] [CrossRef]
  2. Pitt, K.A.; Connolly, R.M.; Meziane, T. Stable isotope and fatty acid tracers in energy and nutrient studies of jellyfish: A review. Hydrobiologia 2009, 616, 119–132. [Google Scholar] [CrossRef]
  3. Kelly, J.R.; Scheibling, R.E. Fatty acids as dietary tracers in benthic food webs. Mar. Ecol. Prog. Ser. 2012, 446, 1–22. [Google Scholar] [CrossRef]
  4. Schiel, D.R.; Steinbeck, J.R.; Foster, M.S. Ten years of induced ocean warming causes comprehensive changes in marine benthic communities. Ecology 2004, 85, 1833–1839. [Google Scholar] [CrossRef]
  5. Steneck, R.S. The ecology of coralline algal crusts: Convergent patterns and adaptative strategies. Annu. Rev. Ecol. Syst. 1986, 17, 273–303. [Google Scholar] [CrossRef]
  6. Gagnon, P.; Matheson, K.; Stapleton, M. Variation in rhodolith morphology and biogenic potential of newly discovered rhodolith beds in Newfoundland and Labrador (Canada). Bot. Mar. 2012, 55, 85–99. [Google Scholar] [CrossRef]
  7. Teichert, S.; Woelkerling, W.; Rüggeberg, A.; Wisshak, M.; Piepenburg, D.; Meyerhöfer, M.; Form, A.; Freiwald, A. Arctic rhodolith beds and their environmental controls (Spitsbergen, Norway). Facies 2014, 60, 15–37. [Google Scholar] [CrossRef]
  8. Górska, B.; Włodarska-Kowalczuk, M. Food and disturbance effects on Arctic benthic biomass and production size spectra. Prog. Oceanogr. 2017, 152, 50–61. [Google Scholar] [CrossRef]
  9. Kuklinski, P.; Berge, J.; McFadden, L.; Dmoch, K.; Zajaczkowski, M.; Nygård, H.; Piwosz, K.; Tatarek, A. Seasonality of occurrence and recruitment of Arctic marine benthic invertebrate larvae in relation to environmental variables. Polar Biol. 2013, 36, 549–560. [Google Scholar] [CrossRef]
  10. Parrish, C.C. Lipids in Marine Ecosystems. ISRN Oceanogr. 2013, 2013, 604045. [Google Scholar] [CrossRef]
  11. Czaja, R., Jr.; Beal, B.; Pepperman, K.; Espinosa, E.P.; Munroe, D.; Cerrato, R.; Busch, E.; Allam, B. Interactive roles of temperature and food availability in predicting habitat suitability for marine invertebrates. Estuar. Coast. Shelf Sci. 2023, 293, 108515. [Google Scholar] [CrossRef]
  12. Schertenleib, K.S.H. Quantifying the impacts of multiple stressors on the production of marine benthic resources. Ph.D. Thesis, Trinity College, Dublin, Irland, 2024. [Google Scholar]
  13. Honya, M.; Kinoshita, T.; Ishikawa, M.; Mori, H.; Nisizawa, K. Seasonal variation in the lipid content of cultured Laminaria japonica: Fatty acids, sterols, β-carotene and tocopherol. J. Appl. Phycol. 1994, 6, 25–29. [Google Scholar] [CrossRef]
  14. Nelson, M.M.; Phleger, C.F.; Nichols, P.D. Seasonal lipid composition in macroalgae of the northeastern Pacific Ocean. Bot. Mar. 2002, 45, 58–65. [Google Scholar] [CrossRef]
  15. Pakhomov, E.A.; McClelland, J.W.; Bernard, K.; Kaehler, S.; Montoya, J.P. Spatial and temporal shifts in stable isotope values of the bottom-dwelling shrimp Nauticaris marionis at the sub-Antarctic archipelago. Mar. Biol. 2004, 144, 317–325. [Google Scholar] [CrossRef]
  16. Foster, M.S. Rhodoliths: Between rocks and soft places. J. Phycol. 2001, 37, 659–667. [Google Scholar] [CrossRef]
  17. Foster, M.S.; McConnico, L.C.; Lundsten, L.; Wadsworth, T.; Kimball, T.; Brooks, L.B.; Medina-López, M.; Riosmena-Rodríguez, R.; Hernandez-Carmona, G.; Vásquez-Elizondo, R.M. Diversity and natural history of a Lithothamnion muelleri-Sargassum horridum community in the Gulf of California. Cienc. Mar. 2007, 33, 367–384. [Google Scholar] [CrossRef]
  18. Kamenos, N.A.; Moore, P.G.; Hall-Spencer, J.M. Attachment of the juvenile queen scallop (Aequipecten opercularis (L.)) to maërl in mesocosm conditions; juvenile habitat selection. J. Exp. Mar. Biol. Ecol. 2004, 306, 139–155. [Google Scholar] [CrossRef]
  19. Steller, D.L.; Cáceres-Martínez, C. Coralline algal rhodoliths enhance larval settlement and early growth of the pacific calico scallop Argopecten ventricosus. Mar. Ecol. Prog. Ser. 2009, 396, 49–60. [Google Scholar] [CrossRef]
  20. Riosmena-Rodriguez, R.; Medina-López, M.A. The role of rhodolith beds in the recruitment of invertebrate species from the southwestern Gulf of California, Mexico. Seaweeds Their Role Glob. Chang. Environ. 2010, 127–138. [Google Scholar]
  21. Kamenos, N.A.; Moore, P.G.; Hall-Spencer, J.M. Nursery-area function of maërl grounds for juvenile queen scallops Aequipecten opercularis and other invertebrates. Mar. Ecol. Prog. Ser. 2004, 274, 183–189. [Google Scholar] [CrossRef]
  22. Steller, D.L.; Riosmena-Rodriguez, R.; Foster, M.S.; Roberts, C.A. Rhodolith bed diversity in the Gulf of California: The importance of rhodolith structure and consequences of disturbance. Aquat. Conserv. Mar. Freshw. Ecosyst. 2003, 13, 5–20. [Google Scholar] [CrossRef]
  23. Riosmena-Rodríguez, R.; Nelson, W.; Aguirre, J. Rhodolith/maërl Beds: A Global Perspective; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  24. Amado-Filho, G.M.; Moura, R.L.; Bastos, A.C.; Salgado, L.T.; Sumida, P.Y.; Guth, A.Z.; Francini-Filho, R.B.; Pereira-Filho, G.H.; Abrantes, D.P.; Brasileiro, P.S. Rhodolith beds are major CaCO3 bio-factories in the tropical South West Atlantic. PLoS ONE 2012, 7, e35171. [Google Scholar] [CrossRef] [PubMed]
  25. Harvey, A.S.; Harvey, R.M.; Merton, E. The distribution, significance and vulnerability of Australian rhodolith beds: A review. Mar. Freshw. Res. 2017, 68, 411–428. [Google Scholar] [CrossRef]
  26. Teed, L.; Bélanger, D.; Gagnon, P.; Edinger, E. Calcium carbonate (CaCO3) production of a subpolar rhodolith bed: Methods of estimation, effect of bioturbators, and global comparisons. Estuar. Coast. Shelf Sci. 2020, 242, 106822. [Google Scholar] [CrossRef]
  27. Marrack, E.C. The relationship between water motion and living rhodolith beds in the southwestern Gulf of California, Mexico. Palaios 1999, 14, 159–171. [Google Scholar] [CrossRef]
  28. Hinojosa-Arango, G.; Maggs, C.A.; Johnson, M.P. Like a rolling stone: The mobility of maërl (Corallinaceae) and the neutrality of the associated assemblages. Ecology 2009, 90, 517–528. [Google Scholar] [CrossRef]
  29. Millar, K.R.; Gagnon, P. Mechanisms of stability of rhodolith beds: Sedimentological aspects. Mar. Ecol. Prog. Ser. 2018, 594, 65–83. [Google Scholar] [CrossRef]
  30. Hacker Teper, S.; Parrish, C.C.; Gagnon, P. Multiple trophic tracer analyses of subarctic Rhodolith (Lithothamnion glaciale) bed trophodynamics uncover bottom-up forcing and benthic-pelagic coupling. Front. Mar. Sci. 2022, 9, 899812. [Google Scholar] [CrossRef]
  31. Piazzi, L.; Pardi, G.; Cinelli, F. Structure and temporal dynamics of a macroalgal assemblage associated with a rhodolith bed of the Tuscan archipelago (Tyrrhenian Sea). Atti Soc. Toscana Sci. Nat. Pisa Mem. 2002, 109, 5–10. [Google Scholar]
  32. Gabara, S.S. Community Structure and Energy Flow within Rhodolith Habitats at Santa Catalina Island, Ca; San José State University: San José, CA, USA, 2014. [Google Scholar]
  33. McConnico, L.A.; Carmona, G.H.; Morales, J.S.M.; Rodríguez, R.R. Temporal variation in seaweed and invertebrate assemblages in shallow rhodolith beds of Baja California Sur, México. Aquat. Bot. 2017, 139, 37–47. [Google Scholar] [CrossRef]
  34. Bélanger, D. Growth Controls of Rhodoliths (Lithothamnion glaciale) and Relationships between Structural Complexity and Macrofaunal Diversity in Subarctic Rhodolith Beds; Memorial University of Newfoundland: St. John’s, NL, Canada, 2020. [Google Scholar]
  35. Gabara, S.S.; Hamilton, S.L.; Edwards, M.S.; Steller, D.L. Rhodolith structural loss decreases abundance, diversity, and stability of benthic communities at Santa Catalina Island, CA. Mar. Ecol. Prog. Ser. 2018, 595, 71–88. [Google Scholar] [CrossRef]
  36. Newell, R.I.E.; Marshall, N.; Sasekumar, A.; Chong, V.C. Relative importance of benthic microalgae, phytoplankton, and mangroves as sources of nutrition for penaeid prawns and other coastal invertebrates from Malaysia. Mar. Biol. 1995, 123, 595–606. [Google Scholar] [CrossRef]
  37. McMeans, B.C.; McCann, K.S.; Humphries, M.; Rooney, N.; Fisk, A.T. Food web structure in temporally-forced ecosystems. TREE 2015, 30, 662–672. [Google Scholar] [CrossRef] [PubMed]
  38. Bélanger, D.; Gagnon, P. Spatiotemporal variability in subarctic Lithothamnion glaciale rhodolith bed structural complexity and macrofaunal diversity. Diversity 2023, 15, 774–794. [Google Scholar] [CrossRef]
  39. Bélanger, D.; Gagnon, P. High growth resilience of subarctic rhodoliths (Lithothamnion glaciale) to ocean warming and chronic low irradiance. Mar. Ecol. Prog. Ser. 2021, 663, 77–97. [Google Scholar] [CrossRef]
  40. Arnold, C.L.; Bélanger, D.; Gagnon, P. Growth resilience of subarctic rhodoliths (Lithothamnion glaciale, Rhodophyta) to chronic low sea temperature and irradiance. J. Phycol. 2022, 58, 251–266. [Google Scholar] [CrossRef]
  41. Budge, S.M.; Parrish, C.C. Lipid biogeochemistry of plankton, settling matter and sediments in Trinity Bay, Newfoundland. II. Fatty acids. Org. Geochem. 1998, 29, 1547–1559. [Google Scholar] [CrossRef]
  42. Parrish, C.C.; Thompson, R.J.; Deibel, D. Lipid classes and fatty acids in plankton and settling matter during the spring bloom in a cold ocean coastal environment. Mar. Ecol. Prog. Ser. 2005, 286, 57–68. [Google Scholar] [CrossRef]
  43. Parrish, C.C.; McKenzie, C.H.; MacDonald, B.A.; Hatfield, E.A. Seasonal studies of seston lipids in relation to microplankton species composition and scallop growth in South Broad Cove, Newfoundland. Mar. Ecol. Prog. Ser. 1995, 129, 151–164. [Google Scholar] [CrossRef]
  44. Graham, D.J.; Midgley, N.G. Graphical representation of particle shape using triangular diagrams: An excel spreadsheet method. Earth Surf. Proc. Landf. 2000, 25, 1473–1477. [Google Scholar] [CrossRef]
  45. Sneed, E.D.; Folk, R.L. Pebbles in the Lower Colorado River, Texas: A study in particle morphogenesis. J. Geol. 1958, 66, 114–150. [Google Scholar] [CrossRef]
  46. Folch, J.; Lees, M.; Sloane Stanley, G.H. A simple method for the isolation and purification of total lipides from animal tissues. JBC 1957, 226, 497–509. [Google Scholar] [CrossRef]
  47. Parrish, C.C. Determination of Total Lipid, Lipid Classes, and Fatty Acids in Aquatic Samples. Lipids Freshw. Ecosyst. 1999, 4–20. [Google Scholar]
  48. Parrish, C.C. Separation of aquatic lipid classes by Chromarod thin-layer chromatography with measurement by Iatroscan flame ionization detection. Can. J. Fish. Aquat. Sci. 1987, 44, 722–731. [Google Scholar] [CrossRef]
  49. Christie, W.W. A simple procedure for rapid transmethylation of glycerolipids and cholesteryl esters. JLR 1982, 23, 1072–1075. [Google Scholar] [CrossRef]
  50. Hamilton, S. Extraction of lipids and derivative formation. Lipid Anal. Pract. Approach 1992, 13–64. [Google Scholar]
  51. Minagawa, M.; Wada, E. Stepwise enrichment of 15N along food chains: Further evidence and the relation between δ15N and animal age. Geochim. Cosmochim. Acta 1984, 48, 1135–1140. [Google Scholar] [CrossRef]
  52. Gale, K.S.P.; Hamel, J.-F.; Mercier, A. Trophic ecology of deep-sea Asteroidea (Echinodermata) from eastern Canada. Deep. Sea Res. Part 1 Oceanogr. Res. Pap. 2013, 80, 25–36. [Google Scholar] [CrossRef]
  53. Cabana, G.; Rasmussen, J.B. Comparison of aquatic food chains using nitrogen isotopes. Proc. Natl. Acad. Sci. USA 1996, 93, 10844–10847. [Google Scholar] [CrossRef]
  54. Grall, J.; Leloch, F.; Guyonnet, B.; Riera, P. Community structure and food web based on stable isotopes (δ15N and δ13C) analysis of a North Eastern Atlantic maërl bed. J. Exp. Mar. Bio. Ecol. 2006, 338, 1–15. [Google Scholar] [CrossRef]
  55. Borgå, K.; Kidd, K.A.; Muir, D.C.G.; Berglund, O.; Conder, J.M.; Gobas, F.A.P.C.; Kucklick, J.; Malm, O.; Powell, D.E. Trophic magnification factors: Considerations of ecology, ecosystems, and study design. IEAM 2012, 8, 64–84. [Google Scholar] [CrossRef] [PubMed]
  56. Connelly, T.L.; Deibel, D.; Parrish, C.C. Trophic interactions in the benthic boundary layer of the Beaufort Sea shelf, Arctic Ocean: Combining bulk stable isotope and fatty acid signatures. Prog. Oceanogr. 2014, 120, 79–92. [Google Scholar] [CrossRef]
  57. Zar, J.H. Biostatistical Analysis; Pearson Education India: Upper Saddle River, NJ, USA, 1999. [Google Scholar]
  58. Guest, M.A.; Nichols, P.D.; Frusher, S.D.; Hirst, A.J. Evidence of abalone (Haliotis rubra) diet from combined fatty acid and stable isotope analyses. Mar. Biol. 2008, 153, 579–588. [Google Scholar] [CrossRef]
  59. Drazen, J.C.; Phleger, C.F.; Guest, M.A.; Nichols, P.D. Lipid composition and diet inferences in abyssal macrourids of the eastern North Pacific. Mar. Ecol. Prog. Ser. 2009, 387, 1–14. [Google Scholar] [CrossRef]
  60. Rohlf, F.J. An empirical comparison of three ordination techniques in numerical taxonomy. Syst. Zool. 1972, 21, 271–280. [Google Scholar] [CrossRef]
  61. Snedecor, G.W.; Cochran, W.G. Statistical Methods. Affiliated East, 8th ed.; Iowa State University Press: Ames, IA, USA, 1994. [Google Scholar]
  62. Parzanini, C. An integrated approach to studying the trophic ecology of a deep-sea faunal assemblage from the Northwest Atlantic; Memorial University of Newfoundland: St. John’s, NL, USA, 2018. [Google Scholar]
  63. Sargent, J.R.; Parkes, R.J.; Mueller-Harvey, I.; Henderson, R.J. Lipid biomarkers in marine ecology. Microbes Sea 1987, 119–138. [Google Scholar]
  64. Graeve, M.; Kattner, G.; Wiencke, C.; Karsten, U. Fatty acid composition of Arctic and Antarctic macroalgae: Indicator of phylogenetic and trophic relationships. Mar. Ecol. Prog. Ser. 2002, 231, 67–74. [Google Scholar] [CrossRef]
  65. Wessels, H.; Karsten, U.; Wiencke, C.; Hagen, W. On the potential of fatty acids as trophic markers in Arctic grazers: Feeding experiments with sea urchins and amphipods fed nine diets of macroalgae. Pol. Biol. 2012, 35, 555–565. [Google Scholar] [CrossRef]
  66. Reuss, N.; Poulsen, L. Evaluation of fatty acids as biomarkers for a natural plankton community. A field study of a spring bloom and a post-bloom period off West Greenland. Mar. Biol. 2002, 141, 423–434. [Google Scholar]
  67. Dalsgaard, J.; St. John, M.; Kattner, G.; Müller-Navarra, D.; Hagen, W. Fatty acid trophic markers in the pelagic marine environment. Adv. Mar. Biol. 2003, 46, 225–340. [Google Scholar]
  68. Nyssen, F.; Brey, T.; Dauby, P.; Graeve, M. Trophic position of Antarctic amphipods—Enhanced analysis by a 2-dimensional biomarker assay. Mar. Ecol. Prog. Ser. 2005, 300, 135–145. [Google Scholar] [CrossRef]
  69. Søreide, J.E.; Falk-Petersen, S.; Hegseth, E.N.; Hop, H.; Carroll, M.L.; Hobson, K.A.; Blachowiak-Samolyk, K. Seasonal feeding strategies of Calanus in the high-Arctic Svalbard region. Deep Sea Res. Part 2 Top. Stud. Oceanogr. 2008, 55, 2225–2244. [Google Scholar] [CrossRef]
  70. Mayzaud, P.; Boutoute, M.; Noyon, M.; Narcy, F.; Gasparini, S. Lipid and fatty acid in naturally occurring particulate matter during spring and summer in a high arctic fjord (Kongsfjorden, Svalbard). Mar. Biol. 2013, 160, 383–398. [Google Scholar] [CrossRef]
  71. Legeżyńska, J.; Kędra, M.; Walkusz, W. Identifying trophic relationships within the high Arctic benthic community: How much can fatty acids tell? Mar. Biol. 2014, 161, 821–836. [Google Scholar] [CrossRef]
  72. Bélanger, D.; Gagnon, P. Low growth resilience of subarctic rhodoliths (Lithothamnion glaciale) to coastal eutrophication. Mar. Ecol. Prog. Ser. 2020, 642, 117–132. [Google Scholar] [CrossRef]
  73. de Young, B.; Sanderson, B. The circulation and hydrography of Conception Bay, Newfoundland. Atmos. Ocean 1995, 33, 135–162. [Google Scholar] [CrossRef]
  74. Robles-Tamayo, C.M.; Valdez-Holguín, J.E.; García-Morales, R.; Figueroa-Preciado, G.; Herrera-Cervantes, H.; López-Martínez, J.; Enríquez-Ocaña, L.F. Sea surface temperature (SST) variability of the eastern coastal zone of the gulf of California. Remote Sens. 2018, 10, 1434. [Google Scholar] [CrossRef]
  75. Smith, V.H.; Tilman, G.D.; Nekola, J.C. Eutrophication: Impacts of excess nutrient inputs on freshwater, marine, and terrestrial ecosystems. Environ. Pollut. 1999, 100, 179–196. [Google Scholar] [CrossRef]
  76. Lloret, J.; Planes, S. Condition, feeding and reproductive potential of white seabream Diplodus sargus as indicators of habitat quality and the effect of reserve protection in the northwestern Mediterranean. Mar. Ecol. Prog. Ser. 2003, 248, 197–208. [Google Scholar] [CrossRef]
  77. Meador, J.P. Bioaccumulation of PAHs in Marine Invertebrates. PAHs Ecotoxicological Perspect. 2003, 147–171. [Google Scholar]
  78. Parzanini, C.; Parrish, C.C.; Hamel, J.-F.; Mercier, A. Functional diversity and nutritional content in a deep-sea faunal assemblage through total lipid, lipid class, and fatty acid analyses. PLoS ONE 2018, 13, e0207395. [Google Scholar] [CrossRef] [PubMed]
  79. Lee, R.F.; Hagen, W.; Kattner, G. Lipid storage in marine zooplankton. Mar. Ecol. Prog. Ser. 2006, 307, 273–306. [Google Scholar] [CrossRef]
  80. Langer, P.D. Some aspects of the biology of three Northwestern Atlantic Chitons: Tonicella rubra, Tonicella marmorea, and Ischnochiton albus (Mollusca: Polyplacophora). Ph.D Thesis, University of New Hampshire, Durham, NH, USA, 1978. [Google Scholar]
  81. Graf, G. Benthic-pelagic coupling: A benthic view. OMBAR 1992, 30, 149–190. [Google Scholar]
  82. Hacker Teper, S. Assessing trophodynamics of a Newfoundland rhodolith (Lithothamnion glaciale) bed using lipid, fatty acid, and stable isotope analyses. Master’s Thesis, Memorial University of Newfoundland, St. John’s, NL, Canada, 2022. [Google Scholar]
  83. Iken, K.; Brey, T.; Wand, U.; Voigt, J.; Junghans, P. Food web structure of the benthic community at the Porcupine Abyssal Plain (NE Atlantic): A stable isotope analysis. Prog. Oceanogr. 2001, 50, 383–405. [Google Scholar] [CrossRef]
  84. Allen, P.L. Feeding behaviour of Asterias rubens (L.) on soft bottom bivalves: A study in selective predation. J. Exp. Mar. Bio. Ecol. 1983, 70, 79–90. [Google Scholar] [CrossRef]
  85. Ramshaw, B.C.; Pakhomov, E.A.; Markel, R.W.; Kaehler, S. Quantifying spatial and temporal variations in phytoplankton and kelp isotopic signatures to estimate the distribution of kelp-derived detritus off the west coast of Vancouver Island, Canada. L&O 2017, 62, 2133–2153. [Google Scholar]
  86. Krumhansl, K.A.; Scheibling, R.E. Production and fate of kelp detritus. Mar. Ecol. Prog. Ser. 2012, 467, 281–302. [Google Scholar] [CrossRef]
  87. Yorke, C.E.; Page, H.M.; Miller, R.J. Sea urchins mediate the availability of kelp detritus to benthic consumers. Proc. R. Soc. Lond. B 2019, 286, 20190846. [Google Scholar] [CrossRef]
  88. Tenore, K.R.; Hanson, R.B.; McClain, J.; Maccubbin, A.E.; Hodson, R.E. Changes in composition and nutritional value to a benthic deposit feeder of decomposing detritus pools. Bull. Mar. Sci. 1984, 35, 299–311. [Google Scholar]
  89. Duggins, D.O.; Eckman, J.E. Is kelp detritus a good food for suspension feeders? Effects of kelp species, age and secondary metabolites. Mar. Biol. 1997, 128, 489–495. [Google Scholar] [CrossRef]
  90. Narciso, L.; da Fonseca, L.C. Growth, survival and fatty acid profile of Nereis diversicolor (OF Müller, 1776) fed on six different diets. Bull. Mar. Sci. 2000, 67, 337–343. [Google Scholar]
  91. Copeland, M.; Wieman, H.L. The chemical sense and feeding behavior of Nereis virens. Sars. Biol. Bull. 1924, 47, 231–238. [Google Scholar] [CrossRef]
  92. Fauchald, K.; Jumars, P.A. The diet of worms: A study of polychaete feeding guilds. OMBAR 1979, 17, 193–284. [Google Scholar]
  93. Nielsen, A.M.; Eriksen, N.T.; Iversen, J.J.L.; Riisgård, H.U. Feeding, growth and respiration in the polychaetes Nereis diversicolor (facultative filter-feeder) and N. virens (omnivorous)—A comparative study. Mar. Ecol. Prog. Ser. 1995, 125, 149–158. [Google Scholar] [CrossRef]
  94. Olivier, M.; Desrosiers, G.; Retière, C.; Brêthes, J.-C. Variations spatio-temporelles de l’alimentation du polychète Nereis virens en zone intertidale (Estuaire maritime du saint-Laurent, Québec). Vie Mil. Life Environ. 1993, 43, 1–12. [Google Scholar]
  95. Olivier, M.; Desrosiers, G.; Caron, A.; Retière, C. Juvenile growth of the polychaete Nereis virens feeding on a range of marine vascular and macroalgal plant sources. Mar. Biol. 1996, 125, 693–699. [Google Scholar] [CrossRef]
  96. Schlacher, T.A.; Connolly, R.M.; Skillington, A.J.; Gaston, T.F. Can export of organic matter from estuaries support zooplankton in nearshore, marine plumes? Aquat. Ecol. 2009, 43, 383–393. [Google Scholar] [CrossRef]
  97. Hedges, J.I.; Keil, R.G.; Benner, R. What happens to terrestrial organic matter in the ocean? Org. Geochem. 1997, 27, 195–212. [Google Scholar] [CrossRef]
  98. Miller, R.J.; Page, H.M. Kelp as a trophic resource for marine suspension feeders: A review of isotope-based evidence. Mar. Biol. 2012, 159, 1391–1402. [Google Scholar] [CrossRef]
  99. Schaffelke, B.; Lüning, K. A circannual rhythm controls seasonal growth in the kelps Laminaria hyperborea and L. digitata from Helgoland (North Sea). Euro. J. Appl. Phycol. 1994, 29, 49–56. [Google Scholar] [CrossRef]
  100. Raven, J.A.; Johnston, A.M.; Kübler, J.E.; Korb, R.; McInroy, S.G.; Handley, L.L.; Scrimgeour, C.M.; Walker, D.I.; Beardall, J.; Vanderklift, M. Mechanistic interpretation of carbon isotope discrimination by marine macroalgae and seagrasses. Funct. Plant Biol. 2002, 29, 355–378. [Google Scholar] [CrossRef] [PubMed]
  101. Schaal, G.; Riera, P.; Leroux, C. Trophic ecology in a Northern Brittany (Batz Island, France) kelp (Laminaria digitata) forest, as investigated through stable isotopes and chemical assays. J. Sea Res. 2010, 63, 24–35. [Google Scholar] [CrossRef]
  102. Hepburn, C.D.; Holborow, J.D.; Wing, S.R.; Frew, R.D.; Hurd, C.L. Exposure to waves enhances the growth rate and nitrogen status of the giant kelp Macrocystis pyrifera. Mar. Ecol. Prog. Ser. 2007, 339, 99–108. [Google Scholar] [CrossRef]
  103. Griffiths, J.R.; Kadin, M.; Nascimento, F.J.A.; Tamelander, T.; Törnroos, A.; Bonaglia, S.; Bonsdorff, E.; Brüchert, V.; Gårdmark, A.; Järnström, M. The importance of benthic–pelagic coupling for marine ecosystem functioning in a changing world. Glob. Change Biol. 2017, 23, 2179–2196. [Google Scholar] [CrossRef]
Figure 1. (A) Ternary diagram of rhodolith (Lithothamnion glaciale) shape relative to purely spheroidal, discoidal, and ellipsoidal rhodoliths (North N = 156, South N = 247; one dot per rhodolith). Rhodoliths were collected in April 2017 at the South and North sites. The position of each rhodolith in the diagram is set by its sphericity, calculated from the length of its longest [L], intermediate [I], and shortest [S] axes. (B) Corresponding proportions of rhodoliths under each of 10 finer shape categories as defined by Sneed and Folk [45].
Figure 1. (A) Ternary diagram of rhodolith (Lithothamnion glaciale) shape relative to purely spheroidal, discoidal, and ellipsoidal rhodoliths (North N = 156, South N = 247; one dot per rhodolith). Rhodoliths were collected in April 2017 at the South and North sites. The position of each rhodolith in the diagram is set by its sphericity, calculated from the length of its longest [L], intermediate [I], and shortest [S] axes. (B) Corresponding proportions of rhodoliths under each of 10 finer shape categories as defined by Sneed and Folk [45].
Diversity 16 00597 g001
Figure 2. Cover of sediment, rhodolith epiphytes, and rhodolith (Lithothamnion glaciale) calculated from photographs of the bed at the South and North sites. The asterisk denotes a statistically significant difference between mean values for the “Rhodolith epiphytes” category.
Figure 2. Cover of sediment, rhodolith epiphytes, and rhodolith (Lithothamnion glaciale) calculated from photographs of the bed at the South and North sites. The asterisk denotes a statistically significant difference between mean values for the “Rhodolith epiphytes” category.
Diversity 16 00597 g002
Figure 4. Biplot of bulk carbon (δ13C) and nitrogen (δ15N) stable isotope ratios of the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites. Components are grouped (circled) based on agglomerative hierarchical cluster analysis, similarly to Figure S1 in [30].
Figure 4. Biplot of bulk carbon (δ13C) and nitrogen (δ15N) stable isotope ratios of the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites. Components are grouped (circled) based on agglomerative hierarchical cluster analysis, similarly to Figure S1 in [30].
Diversity 16 00597 g004
Figure 5. Partial rhodolith food web based on the present study’s analysis of fatty acid and bulk stable isotopes (δ13C and δ15N) of six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 (data pooled across the two study sites given noted high spatiotemporal stability). Energy flow follows the δ15N values on the isotope biplot (see Figure 4). Dotted lines represent settling organic matter (RM = resuspended matter for lines identified with a black dot).
Figure 5. Partial rhodolith food web based on the present study’s analysis of fatty acid and bulk stable isotopes (δ13C and δ15N) of six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 (data pooled across the two study sites given noted high spatiotemporal stability). Energy flow follows the δ15N values on the isotope biplot (see Figure 4). Dotted lines represent settling organic matter (RM = resuspended matter for lines identified with a black dot).
Diversity 16 00597 g005
Table 1. Taxonomical breakdown and abundance of invertebrate macrofauna associated with rhodoliths (Lithothamnion glaciale) collected in April 2017 at the South and North sites. Each phylum’s total abundance (bolded values) includes macrofauna that could not be identified to the genus level.
Table 1. Taxonomical breakdown and abundance of invertebrate macrofauna associated with rhodoliths (Lithothamnion glaciale) collected in April 2017 at the South and North sites. Each phylum’s total abundance (bolded values) includes macrofauna that could not be identified to the genus level.
Phylum/SpeciesMean (±SE) Density
(Individuals kg−1 Rhodoliths)
SouthNorth
Annelida75.9 (8.3)88.1 (21.5)
Myxicola spp.8.9 (4.7)6.7 (3.0)
Nerididae (including Nereis spp.)21.7 (4.7)17.8 (4.4)
Potamilla reniformis30.0 (10.3)41.0 (13.7)
Arthropoda72.8 (13.1)125.8 (19.2)
Amphipoda34.1 (12.4)29.7 (9.1)
Cancer irroratus1.8 (1.0)4.6 (2.2)
Pandalus borealis36.9 (7.5)91.5 (16.3)
Echinodermata452.7 (47.0)603.1 (97.4)
Asterias rubens19.7 (4.1)10.8 (4.2)
Ophiopholis aculeata336.7 (30.8)381.7 (68.9)
Ophiura robusta72.7 (13.8)193.5 (37.7)
Strongylocentrotus droebachiensis22.8 (3.7)17.2 (6.2)
Mollusca427.6 (39.6)603.7 (78.2)
Hiatella arctica152.8 (23.2)171.0 (26.6)
Lacuna vincta1.8 (1.0)2.6 (1.5)
Margarites costalis15.8 (3.6)33.9 (4.8)
Modiolus modiolus14.7 (3.5)22.1 (5.4)
Moelleria costulata8.6 (3.0)11.0 (3.6)
Puncturella noachina21.2 (5.4)51.3 (9.1)
Tonicella marmorea / T. rubra191.6 (27.2)260.9 (36.2)
Turbonilla spp.3.4 (2.3)5.6 (2.9)
Velutina velutina2.4 (1.5)3.5 (1.5)
Nemertea19.3 (5.5)19.5 (4.7)
Sipuncula5.8 (5.3)6.9 (3.6)
Table 2. Significant seasonal changes in all fatty acids in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites.
Table 2. Significant seasonal changes in all fatty acids in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites.
FAFactorF Valuep Value
14:1Month3.240.036
i15:0Month3.460.035
15:0Month6.84<0.010
ai17:0Month9.02<0.010
17:0Month3.900.023
17:1Month12.87<0.010
16:4ω3Month18.33<0.010
20:3ω6Month4.800.010
20:4ω3Month3.300.032
Table 3. Significant seasonal changes in specific components within the eight fatty acids exhibiting at least 70% correlation in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites.
Table 3. Significant seasonal changes in specific components within the eight fatty acids exhibiting at least 70% correlation in the six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites.
FAComponentFactorF Valuep Value
16:1ω7Tonicella spp.Month51.96<0.010
16:1ω7H. arcticaMonth50.22<0.010
16:1ω7S. droebachiensisMonth10.600.003
16:1ω7A. rubensMonth6.310.015
18:0Tonicella spp.Month8.60<0.010
18:0H. arcticaMonth49.64<0.010
18:0A. rubensMonth12.84<0.010
18:0SeawaterMonth12.82<0.010
18:1ω9S. droebachiensisMonth4.770.035
18:1ω7H. arcticaMonth19.62<0.010
18:1ω7A. rubensMonth9.75<0.010
18:1ω7SedimentMonth17.27<0.010
20:1ω11Tonicella spp.Month10.31<0.010
20:1ω11H. arcticaMonth12.26<0.010
20:1ω11Nereis spp.Month14.10<0.010
20:1ω11S. droebachiensisMonth35.5<0.010
20:1ω11SedimentMonth12.55<0.010
20:5ω3Tonicella spp.Month12.21<0.010
20:5ω3H. arcticaMonth16.80<0.010
20:5ω3Nereis spp.Month5.150.029
20:5ω3S. droebachiensisMonth12.99<0.010
20:5ω3A. rubensMonth18.68<0.010
20:5ω3SeawaterMonth35.86<0.010
20:5ω3SedimentMonth8.640.010
20:4ω6Tonicella spp.Month13.55<0.010
20:4ω6H. arcticaMonth12.41<0.010
20:4ω6Nereis spp.Month7.440.010
20:4ω6S. droebachiensisMonth29.14<0.010
20:4ω6L. digitataMonth13.310.017
Table 4. Tukey HSD test of bulk carbon (δ13C) and nitrogen (δ15N) stable isotope ratio groups based on agglomerative hierarchical cluster analysis. The six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites were included in the analysis.
Table 4. Tukey HSD test of bulk carbon (δ13C) and nitrogen (δ15N) stable isotope ratio groups based on agglomerative hierarchical cluster analysis. The six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites were included in the analysis.
ComponentGroupSub-Group13C)
Tukey HSD
15N)
Tukey HSD
L. digitata1-CD
Seawater2aAB
Seawaterb
Seawater, Nereis spp.3aBC
Seawater, Sedimentb
Seawater, Nereis spp.c
Sedimentd
H. arctica4aCB, C
L. glacialeb
A. rubens, Nereis spp., O. aculeata,
Tonicella spp.
5aCA
Nereis spp., O. aculeatab
A. rubens, S. droebachiensisc
Tonicella spp., S. droebachiensisd
A. rubens6-DA
Tonicella spp.7aDB, C
L. digitatab
Table 5. Trophic multiplication factor (TMF) of 33 fatty acids (FA) as calculated from the slope (m) of corresponding linear relationship between FA concentration and bulk nitrogen (δ15N) stable isotope ratio [30]. Six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites were included in the analysis. Only FAs with statistically significant correlation coefficients (r) are shown.
Table 5. Trophic multiplication factor (TMF) of 33 fatty acids (FA) as calculated from the slope (m) of corresponding linear relationship between FA concentration and bulk nitrogen (δ15N) stable isotope ratio [30]. Six animal species, two macroalgal species, and two environmental components sampled in April, July, and December 2017 at the South and North sites were included in the analysis. Only FAs with statistically significant correlation coefficients (r) are shown.
FATMFm (±SE)b (±SE)rp-Value
20:3ω32.761.0 (0.3)8.1 (0.2)0.3<0.001
20:1ω92.651.0 (0.1)7.3 (0.2)0.6<0.001
18:3ω62.580.9 (0.4)8.0 (0.2)0.30.015
20:2ω62.500.9 (0.2)7.3 (0.2)0.4<0.001
20:2a2.110.7 (0.1)7.7 (0.2)0.6<0.001
16:4ω31.410.3 (0)7.6 (0.2)0.5<0.001
17:11.330.3 (0.1)8.1 (0.2)0.1<0.001
20:1ω111.320.3 (0)7.4 (0.1)0.7<0.001
18:1ω111.240.2 (0.1)8.2 (0)0.40.023
22:4ω61.210.2 (0.1)8.2 (0.2)0.10.020
20:4ω61.170.2 (0)7.4 (0.2)0.6<0.001
20:5ω31.140.1 (0)6.2 (0.2)0.6<0.001
18:1ω90.95−0.1 (0)8.7 (0.2)−0.40.002
22:6ω30.94−0.1 (0)8.6 (0.2)0.00.044
16:00.87−0.1 (0)10.2 (0.2)−0.7<0.001
16:1ω70.87−0.1 (0)9.2 (0.2)−0.4<0.001
14:00.85−0.2 (0.1)8.9 (0.3)−0.20.032
i15:00.85−0.2 (0.2)9.1 (0.2)−0.5<0.001
18:4ω30.76−0.3 (0.1)8.9 (0.2)−0.3<0.001
21:5ω30.71−0.3 (0.1)8.6 (0.2)−0.20.011
16:1ω110.64−0.5 (0.1)8.6 (0.2)−0.40.001
18:4ω10.58−0.6 (0.2)8.6 (0.2)−0.20.028
18:3ω30.57−0.6 (0.1)8.8 (0.2)−0.3<0.001
ai15:00.51−0.7 (0.1)8.8 (0.2)−0.3<0.001
16:3ω30.50−0.7 (0.3)8.5 (0.2)−0.40.035
16:3ω40.48−0.7 (0.2)8.7 (0.2)−0.30.001
i16:00.34−1.1 (0.4)8.7 (0.2)−0.20.005
15:00.32−1.1 (0.4)9.0 (0.2)−0.10.004
16:2ω40.22−1.5 (0.4)9.1 (0.3)−0.30.001
16:1ω50.16−1.8 (0.3)9.1 (0.2)−0.4<0.001
16:1ω90.13−2.1 (0.5)8.9 (0.2)−0.4<0.001
14:10.07−2.6 (1.0)8.6 (0.2)−0.20.011
i17:00.07−2.7 (0.3)9.4 (0.2)−0.5<0.001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hacker Teper, S.; Parrish, C.C.; Gagnon, P. Multi-Biomarker Analysis Uncovers High Spatio-Temporal Stability of a Subarctic Rhodolith (Lithothamnion glaciale) Bed Food Web. Diversity 2024, 16, 597. https://doi.org/10.3390/d16100597

AMA Style

Hacker Teper S, Parrish CC, Gagnon P. Multi-Biomarker Analysis Uncovers High Spatio-Temporal Stability of a Subarctic Rhodolith (Lithothamnion glaciale) Bed Food Web. Diversity. 2024; 16(10):597. https://doi.org/10.3390/d16100597

Chicago/Turabian Style

Hacker Teper, Sean, Christopher C. Parrish, and Patrick Gagnon. 2024. "Multi-Biomarker Analysis Uncovers High Spatio-Temporal Stability of a Subarctic Rhodolith (Lithothamnion glaciale) Bed Food Web" Diversity 16, no. 10: 597. https://doi.org/10.3390/d16100597

APA Style

Hacker Teper, S., Parrish, C. C., & Gagnon, P. (2024). Multi-Biomarker Analysis Uncovers High Spatio-Temporal Stability of a Subarctic Rhodolith (Lithothamnion glaciale) Bed Food Web. Diversity, 16(10), 597. https://doi.org/10.3390/d16100597

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop