Abstract
Nephtheidae soft corals are a common component of benthic ecosystems in the Northwest (NW) Atlantic and in the Arctic, but little is known about their functional roles. Here we investigated their role as habitat for basket stars (Ophiuroidea). By examining over 4500 soft coral colonies of Drifa glomerata, Duva florida, Gersemia spp., and Pseudodrifa racemosa collected during trawl surveys in the regions of Newfoundland, Labrador and West Baffin Bay, we show that Nephtheidae soft corals are habitat for juvenile Gorgonocephalus sp. basket stars. Juvenile Gorgonocephalus sp. were found on all soft corals species, but Drifa glomerata was the soft coral species with the highest occurrence of associated juvenile basket stars (32%), with a maximum of 111 basket stars found on a single 4.8 cm tall (contracted state) colony. Individual basket star disc diameter did not change with colony height, but it varied with season for individuals found on D. glomerata. On the other hand, colony height was positively related with maximum disc diameter (D. glomerata), and presence and number of basket stars on a colony (D. glomerata and D. florida). One very young ophiuroid was found within one polyp of a D. florida colony collected in July. We also examined bottom trawl survey data from 2011 to 2017 from Newfoundland to Baffin Bay, and other locations in the Eastern Canadian Arctic and found that for our largest database, 45% of trawl sets that contained basket stars also contained corals, of which 67% were Nephtheidae soft corals. In this study we show that all four Nephtheidae soft coral taxa studied act as habitat for juvenile basket stars, which are conspicuous components of benthic communities in cold-water environments. Soft corals are commonly exposed and vulnerable to bottom fishing gear in the North Atlantic. We urge future studies evaluating the nature of this relationship, the potential role of soft corals as nurseries, and the potential degree of dependency/preference between basket stars and soft corals in relation to other organisms.
Introduction
The functional role of an organism can be defined as the functions provided to the ecosystem (Vaughn and Hakenkamp, 2008), such as ecosystem functions related to impacts on substrate, bentho-pelagic coupling, and associations with other organisms (e.g., provision of microhabitat) (). The function of habitat provision might also be accompanied by other functions, such as protection against predation, substrate for settlement, and camouflage (). Cold-water corals are well known for their habitat-provision potential, as they are often reported in association with other organisms. Gorgonians and black corals have drawn particular attention, as they can grow to large sizes and host a variety of organisms in symbiotic, parasitic, or simply epibiotic relationships (Wahl, 1999; Krieger and Wing, 2002; ,, ; Mosher and Watling, 2009; ; ). With the advent of imagery technologies, in situ evidence of conspicuous associations with cold-water corals have become more common. This is especially true for taxa that can attain large sizes like gorgonians, which are more noticeable from imagery than small taxa like Nephtheidae soft corals, which are relatively small (often < 15 cm tall), with flexible, soft bodies. In contrast to other octocorals, soft corals have even been considered unsuitable to host attached species (), and believed to not host symbionts, in comparison to their warm-water nephtheid counterparts (Watling et al., 2011).
Ophiuroids in the genus Gorgonocephalus are a good example of conspicuous organisms commonly reported in association with cold-water corals (Table 1), including the scleractinians Desmophyllum pertusum (published as Lophelia pertusa) (Jensen and Frederiksen, 1992; Rosenberg et al., 2005; Wheeler et al., 2005) and Madrepora oculata (), the stylasterid Errina antarctica (), the gorgonians Paragorgia arborea (; Metaxas and Davis, 2005; ; ), Paramuricea placomus (Rosenberg et al., 2005; ), Primnoa spp. (Jensen and Frederiksen, 1992; Krieger and Wing, 2002; Wild et al., 2008; ), and also sponges (Klitgaard, 1995). Adult Gorgonocephalus spp. can have disc diameters of 7–10 cm, and they can also be found as free-living organisms (Mortensen, 1927; Patent, 1970a). Juveniles, on the other hand, have been suggested to be strongly linked to cold-water Nephtheidae soft corals (hereinafter referred to as soft corals), as this association has been reported at several instances (e.g., ; Mortensen, 1927; Patent, 1970a; ; ; ).
TABLE 1
| Coral group | Specific coral taxa | Gorgonocephalus taxa | Location | Depth (m) | Notes from the publications | Source |
| Octocorallia | Nephtheidae sp. | G. lamarckii | Newfoundland continental slope | 200–1800 | Most specimens were entangled with a Nephtheidae coral in an apparently parasitic relationship “as noted by others.” | |
| Nephtheidae sp. (published as Eunephthya spp.) | G. eucnemis | North Atlantic | na | Lives parasitic on polyps of Eunephthya spp. When branching of arms begins they leave the coral and may seek adult specimens. | Mortensen, 1927 | |
| Drifa glomerata | Gorgonocephalus spp. | Russia, Norway | na | Observed juveniles of Gorgonocephalus spp. within the polyps of Drifa glomerata (Published as Gersemia glomerata). | ; | |
| Duva florida | Gorgonocephalus sp. | Continental shelf and slope, Atlantic Canada | na | Observed younger individuals amongst branches of soft coral Duva florida. | ||
| Gersemia sp. | G. eucnemis | San Juan Islands, Washington State, USA | 37–183 | Very small (<1 mm disc diameter) specimens found living inside polyps of Gersemia. | Patent, 1970a | |
| Gorgonians | Gorgonocephalus spp. | na | na | Often found clinging to gorgonians and can reach up to 10 cm (disc diameter) | Mortensen, 1927 | |
| Octocoral sp. | G. eucnemis | Guadalupe Island, Mexico | 287.6 | One specimen (100 mm) observed lying on remnants of dead octocoral colony. | ||
| Paragorgia aborea | G. caputmedusae | na | na | General information. | Storm, 1901; Mortensen, 1927; referenced in | |
| Paragorgia aborea | G. lamarckii | na | na | General information. | Storm, 1901; Mortensen, 1927; referenced in | |
| Paragorgia aborea | G. lamarckii | Continental shelf and slope, Atlantic Canada | 300–600 | G. lamarckii observed on outer branches of P. arborea. Observed on young parts of colony. Predominantly large individuals with disc diameter of 15–60 mm. Use of host: substrate. | ; | |
| Paragorgia arborea | G. arcticus | Northeast Channel, Northwest Atlantic | 380–498 | General information. | Metaxas and Davis, 2005 | |
| Paragorgia arborea | G. caputmedusae | na | na | General information. | ||
| Paragorgia arborea | G. eucnemis | na | na | General information. | ||
| Paragorgia arborea | Gorgonocephalus sp. | Norway | 200–400 | Young sections of Paragorgia sp. used as a feeding site by Gorgonocephalus sp. | ||
| Paragorgia arborea | Gorgonocephalus sp. | Continental margin, Northern Norway | 116–863 | Use corals as a substrate for attachment and shelter, for feeding, parasitism. | ||
| Paragorgia sp. | Gorgonocephalus sp. | West of Faroe Islands | 375 | Gorgonocephalus sp. found in same trawl catch with Lophelia and Paragorgia. | ||
| Paramuricea placomus | Gorgonocephalus sp. | Continental margin, Northern Norway | 116–863 | Use corals as a substrate for attachment and shelter, for feeding, parasitism. | ||
| Paramuricea placomus | G. caputmedusae | Oslofjord, Norway | 85–120 | 35 of 90 gorgonocephalids were attached to the gorgonian on one 400 m long transect. | Rosenberg et al., 2005 | |
| Primnoa resedaeformis | G. lamarckii | Continental shelf and slope, Atlantic Canada | 300–600 | Use of host: substrate. | ||
| Primnoa resedaeformis | Gorgonocephalus sp. | Continental margin, Northern Norway | 116–863 | Use corals as a substrate for attachment and shelter, for feeding, parasitism | ||
| Primnoa spp. | G. eucnemis | Gulf of Alaska | 161–365 | G. eucnemis either fed on Primnoa polyps, used the branches for suspension feeding, or sought protection. | Krieger and Wing, 2002 | |
| Primnoidae | Gorgonocephalus sp. | Argentina | 135 | Gorgonocephalus sp. Accompanied by different species of the octocoral family Primnoidae. | Portela et al., 2012 | |
| Ellisellidae | G. chilensis | South Brazil | 20 | One single specimen (28 mm disc diameter) collected associated with the octocoral. | ||
| Stylasteriidae | Errina antarctica | G. chilensis | Chile | 18 | Gorgonocephalus chilensis is commonly encountered on E. antarctica. | |
| Scleractinia | Lophelia pertusa, gorgonians | G. caputmedusae | Norway | na | Facultative, commensal. | Storm, 1901; Mortensen, 1927; ; ; referenced in |
| Lophelia banks | Gorgonocephalus sp. | Faroe Islands | na | The ophiuroid Gorgonocephalus sp. is known from trawl hauls on Lophelia banks in the Faroes. | Jensen and Frederiksen, 1992; Reference BIOFAR trawl surveys not conducted as part of this paper – see | |
| Lophelia pertusa | G. caputmedusae | Oslofjord, Norway | 85–120 | G. caputmedusae were found associated with live and dead L. pertusa and coral rubble. | Rosenberg et al., 2005 | |
| Lophelia pertusa | Gorgonocephalus sp. | Northwest Porcupine Bank, west Ireland | 664 | Photograph observation. | Wheeler et al., 2005 | |
| Lophelia sp. | Gorgonocephalus sp. | West of Faroe Islands | 375 | Gorgonocephalus sp. found in same trawl catch with Lophelia, and Paragorgia. | ||
| Madrepora oculata | Gorgonocephalus sp. | Norway | 200–400 | Submersible observation (photograph) of Gorgonocephalus sp. wrapped around M. oculata. | ||
| Oculina varicosa | Asteroporpa annulata | East-central Florida | 56–84 | Simple-armed gorgonocephalid, observed on live and dead coral mounds (109–207 mm disc diameter) |
Reported associations between Gorgonocephalidae basket stars and cold-water corals.
and Patent (1970a) revealed that the relationship between juvenile basket stars and soft corals begins very early in the basket stars’ ontogeny. observed juvenile Gorgonocephalus sp. (<0.5 mm disc diameter) protruding from polyps of live Drifa glomerata colonies from Kola Fjord (Barents Sea, published as Gersemia glomerata). Decades later, Patent (1970a) reported the presence of juvenile Gorgonocephalus eucnemis inside polyps of live Gersemia sp. colonies from Washington State (West coast of North America), in both August and September 1965. Patent (1970a) suggested that basket star embryos are passively ingested by the soft corals between spawning and 8–9 days of age, and develop within the polyps. Since Patent’s work, studies have focused on reporting the association: found Gorgonocephalus lamarckii individuals entangled with colonies of a Nephtheidae coral in Newfoundland, and young individuals of Gorgonocephalus sp. have been reported on D. florida () and D. glomerata (). Young basket stars have also been observed in situ on colonies of Gersemia sp. from coastal areas in Labrador (∼50–60 m, personal communication B. M. Neves and R. Rangeley, Oceana Canada).
Both soft corals and basket stars are common in the Northwest (NW) Atlantic and the Eastern Canadian Arctic, being found at wide latitudinal ranges and depths of <10–1287 m (Wareham and Edinger, 2007; this study) and 40–1324 m (Neves and Wareham Hayes, 2019; this study), respectively. Soft coral colonies can be found in high densities, being often the most frequent coral taxa in bottom trawl bycatch samples in this region (Wareham and Edinger, 2007; Wareham, 2009; Murillo et al., 2011). Similarly, basket stars can be found in high densities at some locations (e.g., Neves and Wareham Hayes, 2019).
Despite commonly reported in the literature, the basket star-soft coral association has never been quantified. Considering the potential importance of soft corals in the early stages of basket star development, the objectives of this study are: (1) to describe the relationship between the two taxa in terms of epibiont basket star size ranges, occurrence and position on a colony, and (2) to investigate the potential influence of colony height and environmental variables on epibiont basket star occurrence and abundance by soft coral taxa. We examined over 4500 physical colonies from four soft coral taxa (Drifa glomerata, Duva florida, Gersemia spp., and Pseudodrifa racemosa) from the regions of Newfoundland, Labrador, and W Baffin Bay (NW Atlantic). Using 7 years of bottom trawl bycatch data we also investigated the co-occurrence of soft corals and basket stars in the region.
Materials and Methods
Sampling
Most soft coral samples were collected during multispecies trawl surveys conducted by the Department of Fisheries and Oceans Canada (DFO) between the years 2003–2015 in the NW Atlantic and Eastern Canadian Arctic (Table 2 and Figure 1A). These surveys are conducted by DFO’s Newfoundland and Labrador (NL) and Ontario and Prairie (OP, previously Central and Arctic) Regions. Surveys have a stratified random design, with tow speeds of 3 knots and duration of 15 minutes, although OP tows can range between 15 and 30 min.
TABLE 2
| Source | Gear | Dur | Sets (N) | General region | Years | Depth (m) | Months | Ncol |
| Physical samples | ||||||||
| DFO-NL | Campelen (STD) | 15 min | 768 | Labrador, Newfoundland | 2003–2015 | 41–1452 | Jan, April–Jul, Oct–Dec | 2999 |
| DFO-OP | Campelen (MOD) | 15 min | 1 | Ungava Bay | 2011 | 171 | Sept | 1 |
| DFO-OP | Cosmos | 15 min | 114 | W Baffin Bay, Ungava Bay, West of Hudson Strait (WHS) | 2006–2014 | 104–908 | Aug–Nov | 323 |
| DFO-OP | Alfredo | 30 min | 110 | W Baffin Bay and Davis Strait | 2006–2014 | 435–1452 | Sept–Nov | 467 |
| DFO-OP | Gillnet | na | 1 | Baffin Bay | 2006 | 119 | Aug | 1 |
| NSRF | Campelen (MOD) | 15 min | 197 | W Baffin Bay and Davis Strait | 2005–2013 | 116–666 | April, Jul–Aug | 445 |
| Academic | Agassiz | 1.5–3 | 3 | Baffin Bay, Hudson Bay | 2018 | 91–309 | June–July | 33 |
| FOP | Variable | Variable | 113 | Baffin Bay, Labrador, Newfoundland | 2004–2013 | 46–1320 | Jan–Dec | 199 |
| CFERa | Longline, shrimp trap | Variable | 2 | Arctic Bay, Resolute | 2015 | 178, 717 | Aug–Sept | 2 |
| GEACb | Campelen (STD) | 15 min | 1 | S Gulf of St. Lawrence, NE Scotian Shelf, SW Grand Banks | 2014 | na | na | 2 |
| Unknownb | na | na | na | na | na | na | na | 55 |
| Total | 1322 | 4527 | ||||||
| Trawl database | ||||||||
| DFO-NL | Campelen (STD) | 15 min | 8096 | Labrador, Newfoundland | 2011–2017 | 36–1480 | Jan–Dec (except Feb) | Na |
| DFO-OP | Campelen (MOD) | 15 min | 1494 | Davis Strait, Ungava Bay, N Labrador | 2011–2017 | 104–745 | Jul–Sept | Na |
| DFO-OP | Cosmos | 15 min | 151 | W Baffin Bay, Davis Strait, Western Hudson Strait, Ungava Bay | 2011–2013 | 103–950 | Aug–Nov | Na |
| DFO-OP | Alfredo | 30 min | 854 | W Baffin Bay, Davis Strait | 2011–2017 | 298–1492 | Aug–Nov | Na |
| Academic | Agassiz | 1.5–3 min | 88 | Eastern Canadian Arctic (including W Baffin Bay) | 2007–2018 | 35–1453 | June–Oct | Na |
| Total | 10,683 |
Gear (mostly trawls), mean tow duration, number of sets/stations, regions surveyed, years, months, and number of colonies per survey (Ncol) for physical samples and trawl databases.
aExperimental fisheries (CFER, Marine Institute, Memorial University). bGroundfish Enterprise Allocation Council (GEAC) samples and source unknown represent samples whose specific origin data was lost. These samples were only included in the descriptive analyses that did not involve environmental data.
FIGURE 1
The NL Region conducts annual Multispecies Stock Assessment Surveys twice a year, between April-May and October-December. But in some cases, surveys took place during months outside of this range. In this region, surveys use a Campelen 1800 shrimp trawl equipped with rockhopper footgear, with tight rubber disks and spacers mounted along the footrope. The spring survey covers areas within NAFO Divisions 3PL, to a maximum depth of 750 m, while the fall survey covers deep and more northern areas (NAFO Divisions 2HJ, 3KLMNOP) to a maximum depth of 1500 m.
The OP Region conducts annual Groundfish Stock Assessment Surveys twice a year, in August and September, targeting Greenland halibut (Reinhardtius hippoglossoides) and shrimp (Pandalus spp.). These surveys use different trawl types depending on the area being surveyed. Cosmos and Alfredo trawls are used from Hatton Basin to Baffin Bay (NAFO Div. 0AB), to maximum depths of 1500 m. A Campelen trawl is used to survey SE Baffin Shelf - Hatton Basin area (NAFO Div. 0B).
Samples collected during Northern Shrimp Research Foundation (NSRF)/DFO Stock Assessment Surveys were also analyzed. These surveys are conducted from July to September (2006-present), and target Pandalus spp. in Hudson Strait and Ungava Bay, and Eastern Arctic Regions (NAFO Divisions 2G, 0B). This survey uses a modified Campelen 1800 shrimp trawl, sampling at depths < 750 m. A small number of samples from the Baffin Bay region were obtained from miscellaneous surveys, which include the Fisheries Observer Program (FOP), other fisheries surveys (e.g., GEAC), academic research surveys, and experimental fisheries (Table 2).
Over the years, coral samples collected during these surveys have been kept frozen at −20°C at the Northwest Atlantic Fisheries Center (St. John’s, Canada). For this study, soft coral samples were randomly selected from a number of bags containing a miscellaneous of specimens from different trips. Not all specimens from a trip were available, as not all samples caught are kept during the surveys. A total of 4527 soft coral colonies were examined for the presence of clinging epibiont ophiuroids. Samples include colonies of Duva florida (N = 1192), Drifa glomerata (N = 959), Gersemia spp. (N = 2154), and Pseudodrifa racemosa (N = 222) from Newfoundland, Labrador, and Baffin Bay (Figure 1A), from depths ranging between 41 and 1452 m, sampled between January and December. Only whole colonies were considered in this study. Due to challenges in the current taxonomy of Nephtheidae soft corals, we opted for pooling colonies of Gersemia, although we believe to have at least two taxa in this group (Gersemia rubiformis and Gersemia fruticosa), as per current classification (Utinomi, 1961).
Height was measured for each soft coral colony. Ophiuroid positions on the colony, their mass, and disc diameter were also determined. Ophiuroid position was standardized in relation to colony height due to variation in colony sizes. Position was measured as the distance from the center of the basket star’s disc to the colony distal extremity, with relative values of 100% representing individuals found at the very tip of the colony. Position of ophiuroid on a colony was not annotated for 11% of the specimens. This was merely circumstantial, as position was not a variable of interest in the early days of the project. Basket stars were photographed under a stereomicroscope and disc diameters were measured using the image analysis software Image-Pro Plus 7.0. Disc diameters were measured from the distal edge of the radial shields to the edge of the opposite interradial, following Stöhr et al. (2012). A total of 28 ophiuroids (1.4% of measured specimens) could not be measured, as they were deemed too damaged. These individuals were still counted. We provide the final N used in all analyses.
To investigate the presence of embryo or juvenile ophiuroids inside the corals, a subset of the studied colonies of D. glomerata (N = 178) and D. florida (N = 145) collected between April and December was also dissected and examined internally. A total of 10 polyps from three regions in each colony (i.e., top, middle, and base) were dissected and examined under a stereomicroscope (total 30 polyps per colony). Polyps of only a few Gersemia spp. (N = 13) and P. racemosa (N = 11) colonies were examined, as they were often retracted or in a poor physical condition due to their preservation method (i.e., freezing).
Due to the very small size of most ophiuroids, 19 randomly chosen specimens were submitted to the Natural Resources DNA Profiling and Forensic Centre (NRDPFC, Trent University, Canada) for DNA barcoding species identification. COI was the targeted gene, following Layton et al. (2016). Not all samples were amplified using Layton et al.’s primers, so additional primers produced by the NRDPFC (not published) and Kim and Shin (2015) were used.
To gain further insights into the relationship between soft coral and basket stars, we also performed δ13C and δ15N stable isotopes analyses and elemental C and N (%) on seven basket stars (disc diameters between 1.3 and 15 mm) and their soft coral hosts (D. florida and one P. racemosa). Sample preparation followed Sherwood et al. (2008) for both soft corals and basket stars. Briefly, randomly chosen individual soft coral polyps were removed to make a composite sample, while whole basket stars were processed due to their small size (one specimen had an arm sampled). Both dried and crushed basket stars and soft corals were acidified with 5% HCl prior to δ13C analyses to remove carbonates. Analysis of C/N ratios indicated that samples might not have been completely decalcified; therefore only results of the δ15N analyses are presented. Samples were analyzed at the Stable Isotopes Laboratory at Memorial University (St. John’s, Canada) and results are given in ‰.
Although the inspection of colonies for clinging basket stars was done carefully and diligently, we acknowledge that extremely small basket stars could have been missed during inspection, which could have led to false absences. Similarly, the loss of clinging juvenile basket stars during sampling (i.e., trawling) cannot be excluded. However, basket stars were found strongly intermingled with the corals, such that in some cases it was hard to remove them without breaking their arms. Ultimately, we believe that these potential false absences were not common and should not have largely influenced our analyses. In spite of that, we call for caution on the interpretation of absences.
Trawl Database
We have also examined DFO’s multispecies trawl survey and Agassiz trawl Arctic survey databases (Figure 1B) to investigate relationships between soft coral and basket star occurrence in trawl sets (Table 2). A total of 10595 trawl sets were analyzed for coral and basket star presence/absence using the DFO database (NL and OP), from mean depths ranging between 36 and 1492 m, between 2011 and 2017. For the Agassiz trawl Arctic database, data from 88 stations were analyzed, most of which have been resampled through the years during the months of June to October (2007–2018, with exception of 2012). Mean depth for these surveys ranged between 35 and 1453 m. For the trawl datasets we treated soft corals at the family level (i.e., Nephtheidae) due to difficulties with species identification at lower taxonomic levels in the field. Basket stars are not routinely measured in the field, so information on their ontogeny was not available for this study. Generally, the presence of juvenile basket stars on corals is also not accounted for as basket star presence in a trawl set, as soft corals are not inspected at sea.
Five different bottom trawls were used in the surveys included in the trawl database portion of this study: Alfredo, Campelen standard, Campelen modified, Cosmos, and Agassiz (Table 2). These trawls are routinely deployed in different regions, have different dimensions, and can be deployed for different periods of time. For instance, Campelen standard trawls are normally deployed for 15 min during DFO multispecies surveys in the Newfoundland Region, while the much smaller Agassiz trawls are usually deployed for less than 5 min in the Eastern Canadian Arctic as part of academic research surveys. However, our objective was to investigate the presence of soft corals and basket stars in a same set. In this sense, we refrained our analyses to descriptive statistics and considered that any emerging patterns should be relatively consistent across these trawls. Due to a more substantial difference in the gear dimensions, we present the data from the DFO trawls separately from the Agassiz trawl surveys.
Data Analysis
Data analyses were performed using RStudio 1.1.383 (RStudio Team, 2016). Pseudodrifa racemosa and Gersemia spp. were not included in the statistical analyses because of the small number of basket stars associated with them and/or very unbalanced number of observations per season (in some cases not enough to include co-variables). However, the raw data for these taxa is described and graphically presented. Analyses were performed separately for D. glomerata and D. florida, as descriptive statistics revealed that basket star occurrence rates between the two taxa were quite distinct.
Our main explanatory variables of interest were colony height and season. However, in the initial models we also included depth, latitude, year, and interactions (height:depth, height:latitude, season:latitude). We ran different models that included all or some of these variables, and performed model selection by comparing Akaike Information Criteria (AIC) values. Models with the lowest AIC values were selected, except for model 2b shown below, for which the diagnostic plots looked slightly improved than those for the model with the lowest AIC.
The relationship between individual basket star disc diameter (log-transformed) and colony height and season was investigated through Linear Mixed Effects Models (R package nlme, lme function in R, Pinheiro et al., 2020), with colony height and season as fixed effects, and colony id nested in trawl set (station) as random effects (Model 1). Basket star disc diameter was log-transformed to reduce its left-skewed impact on the analyses. We also examined the influence of colony height and season on basket star maximum disc diameter per colony. For D. glomerata the selected model also included depth and latitude (Model 2a), while for D. florida the simplest model (height and season) was selected (Model 2b). For this analysis, only trawl set was included as a random effect, as in this case there was only one disc diameter measurement per colony.
Basket star presence-absence (PA) on a colony was investigated through binomial Generalized Additive Mixed Models (GAMMs) with a clog-log link, with colony height and season as fixed effects, and trawl set (station) as a random effect (Model 3). The clog-log link function was chosen because there were considerably more zeros than ones (Zuur et al., 2009). We also used GAMMs with a Poisson family to investigate the relationship between number of basket stars per colony (in colonies with basket stars) and colony height and season as fixed effects, and trawl set (station) as a random effect (Model 4a). For D. glomerata the selected model also included depth and latitude (Model 4b). We did not include year or interactions in the initial GAMMs, as including these terms generated model convergence errors, likely due to the high complexity of the models for the available data. We used restricted maximum likelihood (REML) as the estimator for the mixed linear models (Zuur et al., 2009). The final models were:
Model 1: lme (logdiscdiam∼height + season, random = ∼1| station/colony, method = “REML”).
Model 2a (D. florida): lme (logmaxdd ∼ height + season, random = ∼1| station, method = “REML”).
Model 2b (D. glomerata): lme (logmaxdd ∼ height∗depth+latitude∗season+year, random = ∼1| station, method = “REML”).
Model 3: gamm(PA∼s(height)+season, random = list(station = ∼1), family = binomial (cloglog)).
Model 4a (D. florida): gamm(abundance∼s(height)+season, random = list (station = ∼1), family = poisson).
Model 4b (D. glomerata): gamm(abundance∼s(height)+season +depth+s(latitude), random = list (station = ∼1), family = poisson).
GAMMs were applied using the gamm function of the R package mgcv (Wood, 2011). When season was identified as a significant source of variation, we used the R function summary to obtain pair-wise comparisons (e.g., between seasons), with the function relevel used to change the order of the baseline, accordingly.
We tested for the presence of outliers and influential observations using Cook’s distance for linear mixed models (CookD function in predictmeans R package). Influential observations were removed as they can have too much influence in the models. For instance, for D. glomerata, analysis on the number of basket stars per colony excluded entries with > 100 basket stars (N = 2), and only considered basket stars < 15 mm in disc diameter (where disc diameter was included, both species). Assumptions of normality, independence, and homogeneity of residuals were assessed through the visualization of diagnostic plots for the linear mixed models (Kabacoff, 2015). For the GAMMs, validation tools include the examination of normalized residuals versus fitted values and explanatory variables (Zuur et al., 2009). The presence of spatial correlation was assessed by adding a spatial correlation structure to the models, and comparing these to those without a correlation structure using AIC values (Zuur et al., 2009). Significance level for all analyses was set at α < 0.05. We included season, depth, and latitude in two selected D. glomerata models (Models 2b and 4b); however since DFO trawl surveys per region depend on time of the year, effects of these variables might be confounded.
In terms of the trawl databases, in order to account for the unbalanced number of trawl sets per region and season, data on soft coral and basket star presence per trawl set is presented as relative proportions. For the Newfoundland and Labrador regions, the data is presented by Ecosystem Production Units (EPUs) (Pepin et al., 2014), and by general region where there is no defined EPU (e.g., Davis Strait, Ungava Bay, West of Hudson Bay, and Baffin Bay, Figure 1B). For the Agassiz trawl Arctic database, data was grouped in 18 main regions (Figure 1B).
Results
Ophiuroid DNA Barcoding
The molecular analysis of 19 basket star specimens yielded 99–100% GenBank matches to Gorgonocephalus arcticus. Two brittle stars were also identified: Ophiopholis aculeata and Ophiocten gracilis (99% matches). Since most specimens were not submitted for molecular analysis (N = 2003 specimens in total), we opted for treating the specimens in this study as Gorgonocephalus sp., but knowing that all specimens included in the molecular analysis yielded very close matches to G. arcticus. GenBank accession numbers for COI sequences of the analyzed specimens can be found in Supplementary Table S1.
General Description
A total of 2003 ophiuroids were found on 543 colonies (12% of all examined colonies), although of these only seven were brittle stars, which were found on only five colonies (in some cases both taxa on a same colony). For this reason, we focus the remaining of the study on basket stars only (N = 1996).
Basket stars were found clinging on colonies of all four studied soft coral taxa (e.g., Figure 2), collected from all studied regions (although limited in Baffin Bay), from depths ranging between 53 and 1317 m (Figure 3). Basket stars were present on colonies collected at every month of the year, except for February when 18 sets were deployed between 2005 and 2006 (one D. glomerata, 24 D. florida, and two Gersemia sp. colonies). No basket stars < 1 mm in disc diameter were found on colonies collected in February, March, or April. Basket stars were mainly found between the 2nd and 3rd quartile of colony height (upper half) in all four soft coral taxa (51–75%, Figure 4A).
FIGURE 2
FIGURE 3
FIGURE 4
Basket Star Disc Diameter
Disc diameters ranged between 0.42 and 19 mm, however, two basket stars 22.2 and 23.6 mm in disc diameter were found on broken D. florida and D. glomerata colonies, respectively, which were not included in the analyses. Disc diameters were classified and displayed in five size classes: ≤ 1, > 1–5, > 5–10, > 10–15, > 15 mm (Figure 4B). In most cases, disc diameters ranged between 1 and 5 mm for basket stars found on all soft coral species (Figure 4B), and for most colonies basket stars belonged to a single size class, although in D. florida and D. glomerata basket stars from up to three size classes on a same colony were identified (Figure 4C).
Individual basket star disc diameter did not significantly vary with colony height for both D. florida and D. glomerata, but it did for season in D. glomerata (Table 3). For this species, individual basket star disc diameters were significantly larger in the Spring than in the Summer and Fall (Figure 5A and Supplementary Table S2).
TABLE 3
| Response | Effect | Num df | Den df | F-value | p-value |
| Log(disc diameter) | |||||
| Drifa glomerata | Intercept | 1 | 1220 | 1086.43 | <0.0001 |
| Height | 1 | 107 | 0.49 | 0.49 | |
| Season | 3 | 188 | 5.09 | 0.002 | |
| Duva florida | Intercept | 1 | 140 | 322.90 | <0.0001 |
| Height | 1 | 42 | 0.50 | 0.48 | |
| Season | 3 | 80 | 1.62 | 0.19 | |
| Maximum disc diameter | |||||
| Drifa glomerata | Intercept | 1 | 183 | 1298.84 | <0.0001 |
| Height | 1 | 106 | 13.63 | <0.0001 | |
| Season | 3 | 183 | 1.01 | 0.39 | |
| Depth | 1 | 183 | 0.77 | 0.38 | |
| Latitude | 1 | 183 | 12.59 | <0.0001 | |
| Year | 1 | 183 | 0.58 | 0.45 | |
| Height:depth | 1 | 106 | 1.74 | 0.19 | |
| Latitude:season | 3 | 183 | 0.96 | 0.41 | |
| Duva florida | Intercept | 1 | 81 | 340.67 | <0.0001 |
| Height | 1 | 45 | 3.64 | 0.06 | |
| Season | 3 | 81 | 1.83 | 0.15 |
ANOVA of fixed effects for linear mixed models estimating colony height and season effects on juveniles basket star disc diameter associated with the Nephtheidae soft corals Drifa glomerata and Duva florida in the Northwest Atlantic.
F-values obtained using the anova function in R. Num df nd Den df = numerator and denominator degrees of freedom, respectively. Values in bold represent significant p-values.
FIGURE 5
On the other hand, maximum disc diameter per colony was significantly related to colony height, and latitude in D. glomerata (Table 3). For D. florida, maximum disc diameter was not significantly related to colony height or season (Table 3).
Basket Star Presence-Absence on Colonies
Basket stars were present on colonies ranging between 10.5 and 147 mm in height (Figure 7). They were more frequent in larger colonies for both D. glomerata, D. florida, and Gersemia spp. (Table 4 and Figure 7). The binomial GAMM indicated that for D. glomerata and D. florida the probability of having associated basket stars increases with colony height up to ∼70 mm for both taxa, when it starts to level off for the former and to slightly decrease for the latter (Figures 6A,B). Descriptive statistics corroborate these observations, as the inter-quartile range of colony height – for colonies where basket stars were present – was 38.8–65.6 mm for D. glomerata and 44.3–69.9 mm for D. florida.
TABLE 4
| Response | Effect | edf∗ | F-value | p-value |
| Basket star presence-absence | Smoother term | |||
| Drifa glomerata | s(colony height) | 2.57 | 24.78 | <0.0001 |
| Parametric terms | Df | F-value | p-value | |
| Season | 3 | 10.79 | <0.0001 | |
| Duva florida | Smoother term | edf | F-value | p-value |
| S(colony height) | 2.30 | 6.19 | 0.002 | |
| Parametric terms | df | F-value | p-value | |
| Season | 3 | 4.26 | 0.005 | |
| Basket star abundance | Smoother term | edf | F-value | p-value |
| Drifa glomerata | s(colony height) | 2.75 | 24.82 | <0.0001 |
| S(latitude) | 4.06 | 5.13 | <0.0001 | |
| Parametric terms | Df | F-value | p-value | |
| Season | 3 | 2.15 | 0.09 | |
| Depth | 1 | 1.32 | 0.25 | |
| Duva florida | Smoother term | edf | F-value | p-value |
| S(colony height) | 2.86 | 4.51 | 0.01 | |
| Parametric terms | df | F-value | p-value | |
| Season | 3 | 0.76 | 0.52 |
Statistical outputs of GAMM models on the effects of soft coral colony height (smoother) and season on juvenile basket star abundance and presence-absence on Nephtheidae soft corals in the Northwest Atlantic.
∗edf = effective degrees of freedom. Approximate significance of smooth terms. Values in bold represent significant p-values.
FIGURE 6
FIGURE 7
Colonies with basket stars were more frequent during the Spring for D. glomerata and D. florida, and during the Winter for Gersemia spp. and P. racemosa (Figure 5B). Basket star presence-absence on D. glomerata colonies significantly changed across seasons (Table 4), except for the Summer-Fall pair (Supplementary Table S3 and Figure 5B). For D. florida, basket star presence-absence on a colony was also significantly different across seasons (Table 4), except for the pairs Fall-Winter and Spring-Summer (Supplementary Table S3 and Figure 5B).
Examination of at least 10410 dissected polyps from 347 colonies yielded the presence of one juvenile ophiuroid (0.34 mm disc diameter) within a D. florida polyp (Figure 8), whose colony was collected in July in the Newfoundland Shelf region, at 306 m. This colony had no epibiont basket stars, and there was no record of a basket star collected in that trawl set. The closest set with a basket star reported was at ∼34 km from this site.
FIGURE 8
Basket Star Abundance
Most colonies had between 1 and 5 clinging basket stars, and only D. glomerata and D. florida had > 20 individuals on a single colony (Figure 4D). Drifa glomerata had both the highest ophiuroid occurrence (31.5%) and maximum number of ophiuroids per colony (N = 111), followed by D. florida and P. racemosa (11 and 10.5%, N = 28 and 6, respectively), and Gersemia spp. (3.9%, N = 8).
Number of basket stars per colony (in colonies with basket stars) significantly varied with colony height, but not with season for both D. glomerata and D. florida (Table 4). There was also an effect of latitude for basket star abundance on D. glomerata (Figure 6C), but not of depth (variables not tested for D. florida). The GAMM for both species indicated that, similarly, to the observed for presence-absence, the number of basket stars per colony increases with colony height up to ∼70 mm, when it levels off (Figures 6D,E).
Trawl Data
When examining the DFO-RV database pooled (2011–2017, Figure 9), 45% of trawl sets containing basket stars (N = 2624) also had corals (N = 1186), of which 67% were Nephtheidae soft corals only (N = 796), followed by a mix of coral taxa that could also include soft corals (16%). The total percentage of trawl sets with basket stars that contained soft corals was 35%. Overall, Labrador shelf was the region with the highest co-occurrence of basket stars and soft corals per trawl set (54% of sets, N = 166), followed by Newfoundland Shelf (48%, N = 206), Grand Banks (36%, N = 275), and southern Newfoundland (31%, N = 120, Figure 10A). Ungava Bay, Davis Strait, and W Baffin Bay had the lowest overall co-occurrence of basket stars and soft corals (18, 19, and 28%, respectively, Figure 10A). Maximum co-occurrence across years took place in the Labrador Shelf region in 2012 (84%, N = 50), and the lowest co-occurrence was identified in Baffin Bay, Davis Strait, and Ungava Bay, where in some years no soft corals were reported in sets with basket stars.
FIGURE 9

Trawl data on presence of basket stars, Nephtheidae soft corals, and both per trawl set (2011–2017). Dots represent DFO RV trawl surveys, while stars represent Agassiz trawl surveys. The maps were cropped to limit the view to the distribution of samples.
FIGURE 10

Relative proportion of coral presence-absence in DFO RV trawl sets containing basket stars by region (number of trawl sets with/without corals relative to total number of sets per region – including basket star absences). Number of sets per region is shown.
Since DFO trawl surveys in the different regions are directed by time of the year (e.g., Southern Newfoundland only during Spring), we could not assess the influence of season on basket star-soft coral co-occurrence per region. But in the Grand Banks region, for which there were sets deployed in every season, basket star presence at depths > 100 m was more frequent in the Summer, followed by Winter, Fall, and Spring (Figure 11). However, basket star-soft coral co-occurrence was more frequent in the Fall, followed by Winter, Summer, and Spring (Figure 11).
FIGURE 11

Relative proportion of DFO trawl sets with basket stars and Nephtheidae soft corals, and sets with basket stars per season in the Grand Banks region (pooled 2011–2017). Total number of sets with basket star and coral co-occurrence is shown on bottom of bars, and number of sets with basket stars is shown on top of bars.
Of the 88 stations sampled with the Agassiz trawl in the Eastern Canadian Arctic, 33% contained both basket stars and soft corals, which all together is comparable to the proportions found in Baffin Bay (28%) based on the larger trawls (see Figure 10). When examining this database by specific regions in the Arctic, stations (sometimes more than one set) with both basket stars and soft corals were identified in the Beaufort Sea (32%, N = 25), Amundsen Gulf (18%, N = 11), Peel Sound (N = 1), Queen Maud Gulf (37.5%, N = 8), Barrow Strait (100%, N = 3), Lancaster Sound (33%, N = 3), North Water Polynya (55.6%, N = 9), and Scott Inlet (50%, N = 2). Unfortunately, physical samples of soft coral colonies from these regions were not inspected for this study and could not be included in our analyzes. However, three physical samples from Hudson Bay/Strait area (not included in trawl database), where there was no co-occurrence between the two taxa based on the trawl database had juvenile basket stars associated to them.
Stable Isotope Analyses
Results of the stable isotope analysis indicate that absolute differences in δ15N content between basket star and soft coral host (i.e., D. florida and one P. racemosa) ranged between 0.28 and 1.1‰ (Supplementary Table S4). δ15N data show both cases of basket stars with higher and lower δ15N than the soft coral hosts (Supplementary Table S4).
Discussion
Basket Star and Soft Coral Co-occurrence
In this study we found juvenile basket stars associated with colonies of all four studied soft coral taxa, from Newfoundland to northern Baffin Bay, at occurrence rates ranging between 4 and 32%, depending on the host taxon. The presence of juvenile basket stars on soft corals is already evidence that the two taxa somehow spatially co-occur. Furthermore, examination of the pooled DFO-RV database showed that 45% of sets with basket stars also had corals, of which 67% were represented by soft corals only. In fact, 35% of the sets with basket stars contained soft corals, which is comparable to the highest proportion of basket star occurrence on physical samples (i.e., 32% in D. glomerata). Examination by region showed the highest relative co-occurrence in the Labrador Shelf Region, at 54%. Although this does not necessarily imply that basket star distribution depends on soft coral presence, our data does show a high degree of co-occurrence between the two taxa.
The northern regions (Ungava Bay, Davis Strait, Baffin Bay) had the lowest co-occurrence of soft coral and basket star in a trawl set (DFO trawls). Although we do not show distribution data by genus for the trawl datasets, examination of soft coral genera distribution based on the examined physical samples indicates that the taxa with the highest basket star occurrence on a colony (i.e., D. glomerata and D. florida) might have their distributions more limited in northern waters. The observation that basket star presence on a colony was limited in Baffin Bay (restricted to latitudes south of 68°N) might be related to the fact that most colonies available/examined from this region were Gersemia spp., which had the lowest basket star occurrence rate (i.e., 4%). In addition, for this taxon colonies with basket stars in the southern sites were more frequent during the Winter, when this region is rarely sampled. On the other hand, trawl sets with both basket stars and nephtheid soft corals were identified in Baffin Bay and elsewhere in the Eastern Canadian Arctic (although at lower rates), indicating that the two taxa do co-occur in this region (even northern of 68°N). The Agassiz trawl data pooled for all sampled Arctic stations indicate a comparable proportion of co-occurrence of the two taxa, in relation to the Baffin Bay data from the larger trawls.
Juvenile basket star occurrence rates on physical colonies were comparable to those found in a study that examined Gorgonocephalus spp. in situ occurrence (i.e., video data) with gorgonian corals in northern Norway (
The preference of juvenile basket stars for colonies of D. glomerata might be associated with different factors. First, our trawl database indicates that basket stars are mainly found along the continental shelf. Unpublished information on soft coral distribution (at the species level) in the region indicates that D. glomerata seems to follow a more similar distribution to that of basket stars, while D. florida is mainly found along the slope break (in preparation, V. Wareham Hayes, personal communication). Therefore, the higher occurrence of basket stars associated with D. glomerata could be related to a higher co-occurrence in the distribution of these two taxa.
Colonies of D. glomerata have the tissue heavily covered with club-shaped spinous CaCO3 sclerites (Utinomi, 1961), whose roughness might increase attachment capacity (
Basket stars were noticeably blended in color with colonies of D. glomerata. The similarity in color between host and basket star could be an advantage for the latter (e.g., camouflage against predators), and protective coloration was also suggested by Patent (1970a). However, colonies of Gersemia spp. and Primnoa resedaeformis gorgonians are also orange in tone, and had the lowest occurrence rate of Gorgonocephalus spp. among the studied colonies in this study and in the study by
Interestingly, young basket stars were commonly observed in situ on live colonies of Gersemia sp. from coastal areas in Labrador, where other soft coral species were not observed (∼50–60 m, personal communication B. M. Neves, and R. Rangeley – Oceana Canada). However, these colonies were observed in situ. It is possible that the difference between alive and dead in the case of this taxon is a factor influencing occurrence rates, since colonies of Gersemia spp. can contract considerably more than the other taxa. This might also corroborate our trawl data observations that co-occurrence rates vary regionally.
We found that most colonies had between 1 and 5 clinging basket stars per colony. Patent (1970a) found a maximum of five basket stars per colony on Gersemia sp., which is comparable to our numbers for this taxon (N = 8). On the other hand, we found a maximum of 111 juvenile basket stars on a single D. glomerata colony.
Soft Corals as a Platform for Recruitment and Nature of the Relationship
In this study we showed that most Gorgonocephalus sp. specimens found associated with soft corals had disc diameters < 5 mm (86%), characterizing them as juveniles. Our largest specimen was 23.6 mm in disc diameter, but these large individuals were not common, and likely still constitute juvenile basket stars, as Patent (1970a) indicated that most adults are usually between 700 and 850 mm in disc diameter. Our observations indicate that soft corals represent an important host for basket stars in the early stages of their lives.
The observation that for most colonies basket stars were in the same size class indicates that these individuals came from recruitment events that occurred close in time. For instance, in the D. glomerata colonies with 111 and 104 individuals, 84 and 76% of them, respectively, had disc diameters between 1 and 5 mm.
Our data show that basket stars find their coral hosts at a very young age, and we believe that soft coral passive ingestion and basket star development within polyps, as suggested in those early studies, is still a plausible hypothesis. It is possible that the fact that our samples were not alive when examined had an influence on our lack of success finding additional individuals within polyps. In fact, the trawl vibration, speed and overall stress could have had an influence on basket star expulsion/evacuation. Furthermore, very few colonies of Gersemia spp. or P. racemosa were inspected for endobiont basket stars. Gersemia sp. was the species studied by Patent (1970a). Soft corals are very small for adult Gorgonocephalus spp. to climb on, and it is hard to conceive that individuals as small as 0.4 mm in disc diameter with undeveloped arms would be able to actively swim or crawl from the seafloor to climb on soft corals. The embryology of Gorgonocephalus spp. is still not well understood, but Patent (1969) suggested that euryalids might not develop pelagic larvae.
Although
Basket Star Size, Abundance, and Position
Colony height had an effect on maximum disc diameter (D. glomerata) and number of basket stars per colony (both species). Although coral size was found not to be a limiting factor in the size of individual juvenile basket stars that it can host, our data indicate that there is a maximum basket star size usually supported by a colony. As mentioned before, larger specimens were found at low frequencies, and our largest specimen was ∼24 mm in dd. According to
Maximum disc diameter also significantly decreased with latitude. Despite the potential for confounding effects between latitude and season, in this case the effect of latitude on maximum disc diameter seems dissociated from season, since there was no significant interaction between latitude and season. Since we do not know whether our Gorgonocephalus spp. represent more than one species, we cannot infer on whether changes in disc diameter with latitude might be related to species-specific differences or to other variables.
Number of basket stars per colony was not influenced by seasonality, but rather by colony height (both species) and latitude (D. glomerata only). For both taxa, number of basket stars and the probability of basket star presence on a colony increased with colony height up to heights of ∼70 mm. Maximum D. glomerata height among the studied colonies was 147 mm, and the two colonies hosting > 100 basket stars were approximately 100 and 120 mm in height. These colonies were not included in the statistical analyses, as they were rare and had too much influence on the data. The GAMM plot indicated that number of basket stars per colony in D. glomerata reaches the lowest point at a latitude of ∼55°N (Labrador, Makkovik area). However, in this analysis, season was weakly not significant (p = 0.09). Considering that borderline p-values should be taken with extra caution in these types of models (Zuur et al., 2009), it is possible that the effect detected for latitude in this case could be associated to season, because of the survey design previously mentioned (location depends on time of the year).
The presence of more basket stars on larger colonies might be a response to space availability, but also increased distance from the seafloor, as most basket stars were found on the upper half of the colonies (as discussed below), indicating a preference for taller surfaces. However, colonies > 70 mm might not necessarily host more basket stars. Furthermore, no data is available on colony height variation between expanded and contracted states (studied species). A better understanding of these variations, which can be significant in soft corals (e.g.,
Individuals of Gorgonocephalus sp. were also found more at the top on the gorgonians P. placomus and P. reseadaeformis, although nearer the base on P. arborea (
Patent (1970a) suggested that basket stars might be feeding on detritus or food from the polyps, as well as on coral eggs, as
Our stable isotopes data showed no indication that increasing size in juvenile basket stars leads to an increase in trophic level. However, this might be the result of our limited sample size, and it might also just indicate that juvenile basket stars of different sizes feed on similar sources. Other studies have yielded enriched δ15N values for individuals of Gorgonocephalus spp. in Arctic waters, in comparison to our results (
Effects of Seasonality
When evaluating the potential effects of seasonality here, it is important to consider that the location of DFO trawl surveys is dependent on time of the year. For instance, Spring surveys take place in southern Newfoundland, while Fall surveys take place in Labrador. Therefore, we cannot exclude that variables affected by season might also be affected by latitude-related dynamics. However, phytoplankton and zooplankton data clearly show the importance of seasonality in this region (
Seasonality – with the constraint highlighted above – seems to play a role on the disc diameter of basket stars found on D. glomerata, as specimens were significantly larger in the Spring than in the Summer and Fall, but no difference was identified between Spring and Winter. Similarly, the relative number of D. glomerata colonies with basket stars peaked during the Spring, with the lowest number of occurrences during Winter. Since we have no information on early growth rates in these ophiuroids, we cannot infer on the timing when most basket stars (embryo/larvae) arrived to the colonies.
In the Grand Banks region, Fall was the season with the highest co-occurrence of basket stars-Nephtheidae soft corals in the trawl sets, followed by Winter, Summer, and Spring. Because soft corals are not mobile organisms, this variation on co-occurrence rates ought to be driven by the basket stars. Seasonality might have an influence on the distribution of adult basket stars, who might migrate seasonally, as do other echinoderms (e.g.,
In the NE Pacific, G. eucnemis (published as G. caryi) has a spawning season of 6 months, from June to November (Patent, 1969), but to our knowledge there is currently no published information on reproduction of Gorgonocephalus spp. in the NW Atlantic. Our observations that juvenile basket stars (<5 mm in dd) were found on coral colonies collected at every month of the year indicate that the spawning season for Gorgonocephalus in this region might be longer than in the Pacific. On the other hand, no basket stars < 1 mm dd were collected in February, March, or April. Reproduction of Gorgonocephalus spp. was not in the scope of this study, and therefore more dedicated studies are needed to shade light on the subject.
The lack of a statistically significant influence of seasonality on basket star metrics for D. florida (i.e., abundance, disc diameter) could be related to their lower basket star occurrence in comparison to D. glomerata. However, despite non-significance, the number of D. florida colonies with basket stars followed the same trend as observed for D. glomerata, with a peak during the Spring and a low in the Winter time. For Gersemia spp., the peak was during Winter time, with the other seasons showing similar trends to one another.
The divergence in co-occurrence per season between physical samples and the trawl sets might be related to the different ontogenetic stages of the basket stars. While physical coral samples contained only juveniles, individuals accounted for in trawl sets are usually large enough to be picked up by hand (Wareham Hayes, personal communication), with sizes not comparable to those found on soft corals. On the other hand, trawl data was considered at the family level, and specific differences by genus (e.g., season of highest frequency of co-occurrence) were not accounted for.
Conservation
Gorgonocephalus arcticus has been found in stomach contents of Greenland sharks (Somniosus microcephalus) from Cumberland Sound, Nunavut (McMeans et al., 2012), and Svalbard, Norway (Leclerc et al., 2012), as well as in stomachs of Northern Wolfish Anarhichas denticulatus (
Basket stars might be relatively long-lived, as some ophiuroids can live for decades (
In this study we showed that in the Northwest Atlantic and Eastern Canadian Arctic juvenile basket stars are consistently found on soft corals, and we believe that further studies are necessary to investigate this association; in particular the potential degree of dependency of basket stars in relation to soft corals and the potential role of the latter as nurseries. This information is paramount to assess the potential impacts that soft coral damage and/or removal (e.g., via bottom-contact fishing) might have on juvenile basket stars, their development, and the success of adult Gorgonocephalus populations. Cold-water Nephtheidae soft corals appear to play an essential role as habitat for juvenile basket stars.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material. Data not present in tables is available upon request to the corresponding author.
Author contributions
BMN and VW contributed to conception and design of the study. BMN performed laboratory work, organized the database, performed the statistical analysis, and wrote the manuscript. EH and VW performed laboratory work and literature review. KH, CG, and PA contributed the Arctic databases on trawl data. All authors contributed to manuscript revision, read and approved the submitted version.
Acknowledgments
We thank the Canadian Coast Guard captains and crew, and the Fisheries and Oceans Canada crew involved in sample collection and processing in the field. Samples collected aboard CCGS Amundsen (Agassiz trawl samples and data) were covered by ArcticNet funding. Laura Wheeland and Brynn Devine provided CFER samples. Special thanks to all involved with the Fisheries Observer Program, which provided invaluable material. We also thank Paula Hawkins, Sarah Olson, and Rebecca Raymonds (DFO-NL) for assistance in the lab, Alison Pye (Memorial University) for the stable isotope analyses, Matthew Harnden for the DNA barcoding analysis, and Hannah Munro for GenBank submissions. Mr. Klaus Hansen for translating the paper by
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.547896/full#supplementary-material
References
1
AfanasyevY. D.NezlinN. P.KostianoyA. G. (2001). Patterns of seasonal dynamics of remotely sensed chlorophyll and physical environment in the Newfoundland region.Remote Sens. Environ.76268–282. 10.1016/S0034-4257(00)00209-1
2
AlbikovskayaL. K. (1983). Feeding characteristics of wolffishes in the labrador-newfoundland region.NAFO Sci. Coun. Stud.635–38.
3
BarbozaC. A. M.MendesF. M.DalbenA.TommasiL. R. (2010). Echinodermata, Ophiuroidea, Gorgonocephalus Leach, 1815: first report of the genus for the Brazilian continental margin.Check List.6289–291. 10.15560/6.2.289
4
BellJ. J. (2008). The functional roles of marine sponges.Estuar. Coast. Shelf Sci.79341–353. 10.1016/j.ecss.2008.05.002
5
BruntseG.TendalO. S. (1999). Marine Biological Investigations and Assemblages of Benthic Invertebrates from the Faroe Islands.Faroe Islands: Kaldbak Marine Biological Laboratory.
6
Buhl-MortensenL.MortensenP. B. (2004a). Crustaceans associated with the deep-water gorgonian corals Paragorgia arborea (L., 1758) and Primnoa resedaeformis (Gunn., 1763).J. Nat. Hist.381233–1247. 10.1080/0022293031000155205
7
Buhl-MortensenL.MortensenP. B. (2004b). Symbiosis in deep-water corals.Symbiosis3733–61.
8
Buhl-MortensenL.MortensenP. B. (2005). “Distribution and diversity of species associated with deep-sea gorgonian corals off Atlantic Canada,” in Cold-Water Corals and Ecosystems, edsFreiwaldA.RobertsJ. M., (Berlin: Springer), 849–879. 10.1007/3-540-27673-4_44
9
Buhl-MortensenL.VanreuselA.GoodayA. J.LevinL. A.PriedeI. G.Buhl-MortensenP.et al (2010). Biological structures as a source of habitat heterogeneity and biodiversity on the deep ocean margins.Mar. Ecol.3121–50. 10.1111/j.1439-0485.2010.00359.x
10
Buhl-MortensenP.Buhl-MortensenL.PurserA. (2017). “Trophic ecology and habitat provision in cold-water coral ecosystems,” in Marine Animal Forests: The Ecology of Benthic Biodiversity Hotspots, edsRossiS.BramantiL.GoriA.CovadongaO., (Cham: Springer), 919–944. 10.1007/978-3-319-21012-4_20
11
Burdon-JonesC.Tambs-LycheH. (1960). Observations on the fauna of the North Brattholmen stone-coral reef near Bergen. Årbok for Universitetet i Bergen, Matematisknaturvitenskaplig Serie4, 1–24.
12
ClippeleL. H.Buhl-MortensenP.Buhl-MortensenL. (2015). Fauna associated with cold water gorgonians and sea pens.Cont. Shelf Res.10567–78. 10.1016/j.csr.2015.06.007
13
DodouD.BreedveldP.De WinterJ. C. F.DankelmanJ.Van LeeuwenJ. L. (2011). Mechanisms of temporary adhesion in benthic animals.Biol. Rev.8615–32. 10.1111/j.1469-185X.2010.00132.x
14
DonsC. (1944). Norges korallrev. Det Kongelige Norske Videnskabers Selskabs Forhandlinger16, 37–82.
15
EmsonR. H.MladenovP. V.BarrowK. (1991). The feeding mechanism of the basket star Gorgonocephalus arcticus.Can. J. Zool.69449–455. 10.1139/z91-070
16
FederH. M.IkenK.BlanchardA. L.JewettS. C.SchonbergS. (2011). Benthic food web structure in the southeastern chukchi sea: an assessment using δ13C and δ15N analyses.Polar Biol.34521–532. 10.1007/s00300-010-0906-9
17
FedotovD. M. (1924). Einige beobachtungen tiber die biologie und metamorphose von Gorgonocephalus.Zool. Anz.61303–311.
18
Food and Agriculture Organization of the United Nations [FAO] (2020). VME indicators, thresholds and encounter responses adopted by R(F)MOs in force during 2019. Available at: http://www.fao.org/in-action/vulnerable-marine-ecosystems/vme-indicators/en/(accessed March 31, 2020).
19
FörsterraG.HäussermannV.LaudienJ. (2017). “Animal forests in the chilean fjords: discoveries, perspectives, and threats in shallow and deep waters,” in Marine Animal Forests: The Ecology of Benthic Biodiversity Hotspots, edsRossiS.BramantiL.GoriA.CovadongaO., (Cham: Springer), 277–313. 10.1007/978-3-319-21012-4_3
20
FourgonD.JangouxM.EeckhautI. (2007). Biology of a “babysitting” symbiosis in brittle stars: analysis of the interactions between Ophiomastix venosa and Ophiocoma scolopendrina.Invertebr. Biol.126385–395. 10.1111/j.1744-7410.2007.00108.x
21
FreemanS. M.RichardsonC. A.SeedR. (2001). Seasonal abundance, spatial distribution, spawning and growth of Astropecten irregularis (Echinodermata?: Asteroidea).Estuar. Coast. Shelf Sci.5339–49. 10.1006/ecss.2000.0758
22
GageJ. D. (1990). Skeletal growth markers in the deep-sea brittle stars Ophiura ljungmani and Ophiomusium lymani.Mar. Biol.104427–435. 10.1007/BF01314346
23
GiliJ.-M.PetraitisP. S. (2009). Seasonal Dynamics.Berlin: Springer, 191–200.
24
GirardF.FuB.FisherC. R. (2016). Mutualistic symbiosis with ophiuroids limited the impact of the Deepwater Horizon oil spill on deep-sea octocorals.Mar. Ecol. Prog. Ser.54989–98. 10.3354/meps11697
25
GislénT. (1924). Echinoderm studies. Zool. Bidr. Upps. 9, 1–316.
26
HaedrichR.MaunderJ. (1984). “The echinoderm fauna of the Newfoundland continental slope,” in Proceedings of the 5th International Echinoderm Conference, Galway, 37–46. 10.1201/9781003079224-7
27
HäussermannV.FösterraG. (2007). Extraordinary abundance of hydrocorals (Cnidaria, Hydrozoa, Stylasteridae) in shallow water of the Patagonian fjord region.Polar Biol.30487–492. 10.1007/s00300-006-0207-5
28
HellströmM.BenzieJ. A. H. (2011). Robustness of size measurement in soft corals.Coral Reefs30:787. 10.1007/s00338-011-0760-4
29
HendlerG.GrygierM. J.MaldonadoE.DentonJ. (1999). Babysitting brittle stars: heterospecific symbiosis between ophiuroids (Echinodermata).Invertebr. Biol.118190–201. 10.2307/3227060
30
HendlerG.MillerJ. E. (1984). Feeding behavior of Asteroporpa annulata, a gorgonocephalid brittlestar with unbranched arms.Bull. Mar. Sci.34449–460.
31
Herrero-PérezrulM. A. D.Granja-FernándezR.Hoyos-PadillaM.Reyes-BonillaH. (2014). New record of the basket star Gorgonocephalus eucnemis (Ophiuroidea: Gorgonocephalidae) at the Pacific coast of Mexico.Mar. Biodivers. Rec.7:e33. 10.1017/S1755267214000323
32
HobsonK. A.FiskA.KarnovskyN.HolstM.GagnonJ.-M.FortierM. (2002). A stable isotope (δ13C, δ15N) model for the North Water food web: implications for evaluating trophodynamics and the flow of energy and contaminants.Deep Sea Res. Part II Top. Stud. Oceanogr.495131–5150. 10.1016/S0967-0645(02)00182-0
33
IkenK.BluhmB.DuntonK. (2010). Benthic food-web structure under differing water mass properties in the southern Chukchi Sea.Deep Sea Res. Part II Top. Stud. Oceanogr.5771–85. 10.1016/j.dsr2.2009.08.007
34
JensenA.FrederiksenR. (1992). The fauna associated with the bank-forming deepwater coral Lophelia pertusa (Scleractinaria) on the Faroe shelf.Sarsia7753–69. 10.1080/00364827.1992.10413492
35
JørgensenL. L.PlanqueB.ThangstadT. H.CertainG. (2015). Vulnerability of megabenthic species to trawling in the Barents Sea.ICES J. Mar. Sci. J. du Cons.73i84–i97. 10.1093/icesjms/fsv107
36
KabacoffR. (2015). R in Action: Data Analysis and Graphics with R, 2nd Edn. Shelter Island, NY: Manning.
37
KimD.ShinS. (2015). A newly recorded basket star of genus Gorgonocephalus (Ophiuroidea: Euryalida: Gorgonocephalidae) from the East Sea, Korea.Anim. Syst. Evol. Diversity31:311. 10.5635/ased.2015.31.4.311
38
KlitgaardA. B. (1995). The fauna associated with outer shelf and upper slope sponges (Porifera, Demospongiae) at the Faroe Islands, northeastern Atlantic.Sarsia801–22. 10.1080/00364827.1995.10413574
39
KriegerK. J.WingB. L. (2002). Megafauna associations with deepwater corals (Primnoa spp.) in the Gulf of Alaska.Hydrobiologia47183–90. 10.1023/A:1016597119297
40
LaytonK. K. S.CorstorphineE. A.HebertP. D. N. (2016). Exploring Canadian echinoderm diversity through DNA barcodes.PLoS One11:e0166118. 10.1371/journal.pone.0166118
41
LeclercL.-M. E.LydersenC.HaugT.BachmannL.FiskA. T.KovacsK. M. (2012). A missing piece in the Arctic food web puzzle? Stomach contents of Greenland sharks sampled in Svalbard, Norway.Polar Biol.351197–1208. 10.1007/s00300-012-1166-7
42
LongS.Sparrow-ScinoccaB.BlicherM. E.Hammeken ArboeN.FuhrmannM.KempK. M.et al (2020). Identification of a soft coral garden candidate Vulnerable Marine Ecosystem (VME) using video imagery, Davis Strait, West Greenland.Front. Mar. Sci.7:460. 10.3389/fmars.2020.00460
43
MaunaA. C.BottoF.FrancoB.SchwartzJ. M.AchaE. M.LastaM. L.et al (2011). Shifts in an epibenthic trophic web across a marine frontal area in the Southwestern Atlantic (Argentina).J. Sea Res.66248–255. 10.1016/j.seares.2011.08.005
44
McMeansB. C.ArtsM. T.FiskA. T. (2012). Similarity between predator and prey fatty acid profiles is tissue dependent in Greenland sharks (Somniosus microcephalus): implications for diet reconstruction.J. Exp. Mar. Bio. Ecol.42955–63. 10.1016/j.jembe.2012.06.017
45
MetaxasA.DavisJ. (2005). Megafauna associated with assemblages of deep-water gorgonian corals in Northeast Channel, off Nova Scotia, Canada.J. Mar. Biol. Assoc. U.K.851381–1390. 10.1017/S0025315405012567
46
MortensenT. H. (1927). Handbook of the Echinoderms of the British Isles.Humphrey Milford: Oxford University Press.
47
MosherC. V.WatlingL. (2009). Partners for life: a brittle star and its octocoral host.Mar. Ecol. Prog. Ser.39781–88. 10.3354/meps08113
48
MurilloF. J.DuránM. P.AltunaA.SerranoA. (2011). Distribution of deep-water corals of the Flemish Cap, Flemish Pass, and the Grand Banks of Newfoundland (Northwest Atlantic Ocean): interaction with fishing activities.ICES J. Mar. Sci. J. du Cons.68319–332. 10.1093/icesjms/fsq071
49
NadtochyV. A.KolpakovN. V.KorneichukI. A. (2017). The distribution of macrozoobenthos taxa, as potential indicators of Vulnerable Marine Ecosystems in the Western Bering Sea: 1. Anadyr Bay area.Russ. J. Mar. Biol.43555–567. 10.1134/S1063074017070033
50
NevesB. M.Wareham HayesV. (2019). “Distribution patterns of basket stars (Gorgonocephalidae) in the Northwest Atlantic,” in Proceedings of the Benthic Ecology Meeting, St. John’s.
51
PatentD. H. (1969). The reproductive cycle of Gorgonocephalus caryi (Echinodermata; Ophiuroidea).Biol. Bull.136241–252. 10.2307/1539817
52
PatentD. H. (1970a). Life history of the basket star, Gorgonocephalus eucnemis (Müller & Troschel) (Echinodermata; Ophiuroidea).Ophelia8145–159. 10.1080/00785326.1970.10429556
53
PatentD. H. (1970b). The early embryology of the basket star Gorgonocephalus caryi (Echinodermata, Ophiuroidea).Mar. Biol.6262–267. 10.1007/BF00347235
54
PepinP.ColbourneE.MailletG. (2011). Seasonal patterns in zooplankton community structure on the Newfoundland and Labrador Shelf.Prog. Oceanogr.91273–285. 10.1016/j.pocean.2011.01.003
55
PepinP.HigdonJ.Koen-AlonsoM.FogartyM.OllerheadN. (2014). Application of Ecoregion Analysis to The Identification of Ecosystem Production Units (EPUs) in the NAFO Convention Area. NAFO SCR Doc. 14/069. Dartmouth: NAFO.
56
PinheiroJ.BatesD.DebRoyS.SarkarD.R Core Team. (2020). nlme: Linear and Nonlinear Mixed Effects Models. R Package Version 3. 1–149. Available online at: https://CRAN.R-project.org/package=nlme
57
PortelaJ.AcostaJ.CristoboJ.MuozA.ParraS.IbarrolaT.et al (2012). “Management strategies to limit the impact of bottom trawling on VMEs in the High Seas of the SW Atlantic,” in Marine Ecosystems, ed.AcostaJ., (Rijeka: IntechOpen).
58
RosenbergR.DupontS.LundälvT.SköldH. N.NorkkoA.RothJ.et al (2005). Biology of the basket star Gorgonocephalus caputmedusae (L.).Mar. Biol.148:43. 10.1007/s00227-005-0032-3
59
RStudio Team (2016). RStudio: Integrated Development Environment for R. Available at: http://www.rstudio.com/
60
SherwoodO. A.JamiesonR. E.EdingerE. N.WarehamV. E. (2008). Stable C and N isotopic composition of cold-water corals from the Newfoundland and Labrador continental slope: examination of trophic, depth and spatial effects.Deep Sea Res. Part I Oceanogr. Res. Pap.551392–1402. 10.1016/j.dsr.2008.05.013
61
SimpsonM. R.SherwoodG. D.MelloL. G. S.MiriC. M.KulkaD. W. (2013). Feeding Habits and Trophic Niche Differentiation in Three Species of Wolffish (Anarhichas sp.) Inhabiting Newfoundland and Labrador Waters. DFO Can. Sci. Advis. Sec. Res. Doc. 2013/056. Canada: DFO.
62
SlatteryM.McClintockJ. B.HeineJ. N. (1995). Chemical defenses in Antarctic soft corals: evidence for antifouling compounds.J. Exp. Mar. Bio. Ecol.19061–77. 10.1016/0022-0981(95)00032-M
63
StöhrS.O’HaraT. D.ThuyB. (2012). Global diversity of brittle stars (Echinodermata: Ophiuroidea).PLoS One7:e0031940. 10.1371/journal.pone.0031940
64
StormV. (1901). “Oversigt over Throndheimsfjordens fauna (med et kort),” in Trondhjems Biologiske Station, Meddelelser fra stationsanleggets arbeidskomite, ed.Moe’s BogH. (Trondhjem: Accidentstrykkeri), 20.
65
TuronX.CodinaM.TarjueloI.UrizM. J. (2000). Mass recruitment of Ophiothrix fragilis (Ophiuroidea) on sponges: settlement patterns and post-settlement dynamics.Mar. Ecol. Prog. Ser.200201–212. 10.3354/meps200201
66
UtinomiH. (1961). A revision of the nomenclature of the family Nephtheidae (Octocorallia: Alcyonacea) II. The boreal genera Gersemia, Duva, Drifa and Pseudodrifa.Publ. Seto Mar. Biol. Lab.9229–246.
67
Vander ZandenM. J.RasmussenJ. B. (2001). Variation in δ15N and δ13C trophic fractionation: implications for aquatic food web studies.Limnol. Oceanogr.462061–2066. 10.4319/lo.2001.46.8.2061
68
VaughnC. C.HakenkampC. C. (2008). The functional role of burrowing bivalves in freshwater ecosystems.Freshw. Biol.461431–1446. 10.1046/j.1365-2427.2001.00771.x
69
WahlM. (1999). The predominantly facultative nature of epibiosis: experimental and observational evidence.Mar. Ecol. Prog. Ser.18759–66.
70
WarehamV. E. (2009). “Updates on deep-sea coral distributions in the Newfoundland and Labrador and Arctic Regions, Northwest Atlantic,” in The Ecology of Deep-Sea Corals of Newfoundland and Labrador Waters: Biogeography, Life History, Biogeochemistry, and Relation to Fishes, edsGilkinsonK.EdingerE., (London: Canadian Technical Report on Fisheries Aquatic Sciences).
71
WarehamV. E.EdingerE. N. (2007). Distribution of deep-sea corals in the Newfoundland and Labrador region, Northwest Atlantic Ocean.Bull. Mar. Sci.81289–313.
72
WatlingL.FranceS. C.PanteE.SimpsonA. (2011). “Biology of deep-water octocorals,” in Advances in Marine Biology Advances in Marine Biology, ed.SheppardC., (Cambridge, MA: Academic Press), 41–122.
73
WheelerA. J.BeckT.ThiedeJ.KlagesM.GrehanA.MonteysF. X. (2005). “Deep-water coral mounds on the Porcupine Bank, Irish Margin: preliminary results from the Polarstern ARK-XIX/3a ROV cruise,” in Cold-Water Corals and Ecosystems, edsFreiwaldA.RobertsJ. M., (Berlin: Springer), 393–402.
74
WildC.MayrC.WehrmannL.SchöttnerS.NaumannM.HoffmannF. (2008). Organic matter release by cold water corals and its implication for fauna–microbe interaction.Mar. Ecol. Prog. Ser.37267–75.
75
WilliamsD.AndersenR. J.Van DuyneG. D.ClardyJ. (1987). Cembrane and pseudopterane diterpenes from the soft coral Gersemia rubiformis.J. Org. Chem.52332–335. 10.1021/jo00379a002
76
WoodS. N. (2011). Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models.J. R. Stat. Soc. Ser. B Stat. Methodol.733–36. 10.1111/j.1467-9868.2010.00749.x
77
ZuurA.IenoE. N.WalkerN.SavelievA. A.SmithG. M. (2009). Mixed Effects Models and Extensions in Ecology with R.New York, NY: Springer.
Summary
Keywords
cold-water corals, Nephtheidae, functional role, habitat, basket stars, Gorgonocephalus, juveniles
Citation
Neves BM, Wareham Hayes V, Herder E, Hedges K, Grant C and Archambault P (2020) Cold-Water Soft Corals (Cnidaria: Nephtheidae) as Habitat for Juvenile Basket Stars (Echinodermata: Gorgonocephalidae). Front. Mar. Sci. 7:547896. doi: 10.3389/fmars.2020.547896
Received
01 April 2020
Accepted
25 August 2020
Published
09 October 2020
Volume
7 - 2020
Edited by
Santiago Herrera, Lehigh University, United States
Reviewed by
Steven Y. Litvin, Monterey Bay Aquarium Research Institute (MBARI), United States; Les Watling, University of Hawai‘i at Mānoa, United States; Fanny Girard, Institut Français de Recherche pour l’Exploitation de la Mer (IFREMER), France
Updates

Check for updates
Copyright
© 2020 Neves, Wareham Hayes, Herder, Hedges, Grant and Archambault.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bárbara de Moura Neves, barbara.neves@dfo-mpo.gc.ca
This article was submitted to Deep-Sea Environments and Ecology, a section of the journal Frontiers in Marine Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.