Sea pen ( Pennatulacea ) aggregations on the northern Spanish shelf : distribution and faunal assemblages

The sea pens Funiculina quadrangularis, Pennatula aculeata, Pennatula phosphorea, Pteroeides spinosum and Veretillum cynomorium are known to be characteristic species of the habitat ‘pennatulacean communities on circalittoral and upper bathyal soft bottoms’, described recently on the northern Atlantic Spanish shelf under the terms of the EU directives. Distribution, abundance and environmental data collected from 1995 to 2010 in the Northern Spanish Shelf Groundfish Surveys in the Cantabrian Sea and off Galicia (SPNGFS) were analysed. Although low densities of sea pens were commonly found over the years, some aggregations were shown. Multivariate analysis separated the aggregations of F. quadrangularis and of Pennatula spp. found on circalittoral soft bottoms. In addition, evidence of aggregations of F. quadrangularis on upper bathyal sites was also found. Depth and longitude, which were closely related to organic matter and sediment particle size, were key factors in determining the distribution and abundance of these aggregations. However, very little information on distribution and density of the species P. spinosum and V. cynomorium was obtained. Details of community composition of F. quadrangularis and Pennatula spp. were studied and compared with those of the OSPAR habitat “sea pen and burrowing megafauna communities” in other areas. The results improve our knowledge concerning sea pen communities in the southern Bay of Biscay and contribute to a global assessment of the status of these communities in the northeast Atlantic.


INTRODUCTION
Recent studies and environmental legislation are involving sea pens in habitat conservation (OSPAR 2010, NOAA 2014, Kenchington et al. 2014).Despite their ecological value, pennatulacean coral habitats have been overlooked compared with cold-water coral ecosystems and studies of the latter mainly focus on stony corals (Order Scleractinia), sea fans and soft corals (Order Alcyonacea).Examples of sea pens forming fields (Aguilar et al. 2009, Porporato et al. 2014, Kenchington et al. 2014) and acting as nurseries for fish larvae (Pirtle 2005, Baillon et al. 2012) or as indicators of quality of mud habitats (MacDonald et al. 1996) highlight their role in benthic ecosystems.Like other corals, sea pens are slow-growing species that are particularly vulnerable because recovery after disturbances can take a long time (Kaiser et al. 2006, de Moura Neves et al. 2015).Moreover, in contrast with other corals, they are able to anchor in soft sediment (Hughes 1998), offering a shelter in a bare environment and acting as an indicator of the degradation of these sediments, which commonly sustain intensive fishing pressure or other disturbances (Thrush and Dayton 2002).The OSPAR Convention assessed the threatened habitat "sea pen and burrowing megafauna communities" in the northeast Atlantic affected by Nephrops trawling (Hughes 1998, Greathead et al. 2007, OSPAR 2010).Moreover, sea pens are involved in the soft-bottom coral garden habitat that is also considered threatened and/or declining (OSPAR 2008).Outside the OSPAR region, in eastern Canada, the important role of deep-sea pennatulacean corals as biogenic habitats has been reported (Kenchington et al. 2010, Baillon et al. 2014), and in the USA, under the reauthorized Magnuson-Stevens fishery Conservation and Management Act, deep-sea coral areas, including sea pens, have been protected from fishing (NOAA 2014).The Northwest Atlantic Fisheries Organization (NAFO) has also protected sea pen fields on the Flemish Cap.
After OSPAR recommendations to increase knowledge on the "sea pen and burrowing megafauna habitat" in areas where there are insufficient data and intensive fishing, such as the Cantabrian Sea (Sánchez and Olaso 2004, OSPAR 2010, ICES 2012), and under the European Marine Strategy Framework Directive (EU 2008) considerable advances in habitat definition have been made in Spain in the last few years, describing 'the pennatulacean communities on circalittoral and bathyal soft bottoms' in the North Atlantic Marine Area (MAGRAMA 2012).Despite these efforts, there is still limited information on this habitat in the area and it receives little attention and protection.
In the Bay of Biscay and nearby areas, 21 species of sea pens have been reported, 13 of them occurring off Spain south of the 44°N parallel (Altuna 2010, Supplementary Material Table S1).During several surveys undertaken every autumn by the Spanish Institute of Oceanography (IEO) in the Cantabrian Sea (southern sector of the Bay of Biscay) and off Galicia, five of these species, Funiculina quadrangularis, Pennatula aculeata, Pennatula phosphorea, Pteroeides spinosum and Veretillum cynomorium, were found on circalittoral and upper bathyal soft bottoms, sometimes in abundance.Although some of these species are of conservation significance (Hughes 1998, OSPAR 2010), knowledge of them is uneven in this area, with little information on their abundance, distribution and habitat requirements.
F. quadrangularis has a worldwide range, being common throughout the North Atlantic and the Mediterranean Sea (Manuel 1981, Gili andPagès 1987), with a wide bathymetric range from 20 m to depths of over 3000 m (Tixier-Durivault and d'Hondt 1975, Manuel 1981, López-González et al. 2001).It is a sensitive species considered to be of conservation significance.This species is reported to be especially vulnerable to disturbances because of its inability to withdraw rapidly into the sediment (Hughes 1998).According to the European Nature Information System (EUNIS), it is a characterizing species of the circalittoral biotope 'sea pen and burrowing megafauna', as well as being common in the bathyal facies of soft muds associated with Aporrhais serresianus in the Mediterranean Sea.Furthermore, with the implementation of the MSFD in Spain, specifically in the North Atlantic Marine Area (MAGRAMA 2012), this species is also described as a characterizing species of 'the pennatulacean communities on circalittoral and bathyal soft bottoms' in the northern Atlantic Spanish area.
The genus Pennatula has a worldwide distribution (Williams 2011), and three out of the five species listed by Williams (2014) from the European Register of Marine Species area occur in the Cantabrian Sea off Spain: namely P. aculeata, P. grandis and P. phosphorea.Records of P. rubra from this area are misidentifications.The species P. aculeata and P. phosphorea are often uncertainly classified (Manuel 1981), and the latter has been confused with other species (Hughes 1998).Thus, tracing their distributions is not easy despite the fact that they are accepted species (Williams 2014).Both species may co-occur in some areas and at some depths in the southern sector of the Bay of Biscay and Galicia (Supplementary Material Table S1).P. phosphorea is widespread in the North Atlantic and in the Mediterranean Sea, with several records from the whole Bay of Biscay (Tixier-Durivault andd'Hondt 1975, Gili andPagès 1987), and is also involved in conservation issues (OSPAR 2010).According to EUNIS classification, it occurs together with Virgularia mirabilis in the biotope 'sea pen and burrowing megafauna' in the North and Irish Sea, but also in the facies of sticky muds with V. mirabilis in the Mediterranean Sea.It may also be found in the biotope 'Myrtea spinifera and polychaetes in offshore circalittoral sandy mud' in the North and Irish Sea.
On the other hand, P. aculeata has been mainly reported in the North Atlantic, including the Bay of Biscay, but there are some Mediterranean records as well (Tixier-Durivault and d'Hondt 1975, Langton et al. 1990, Ocaña et al. 2000).A review of the literature (Supplementary Material Table S1) indicates that it can reach deeper depths than P. phosphorea in the Bay of Biscay and nearby areas, and it has not been associated with a given habitat.
The two genera Pteroeides and Veretillum are reported in the eastern Atlantic and West Pacific and described to be bathymetrically restricted to shallow waters (Williams 2011).Specifically, the geographical distribution of P. spinosum and V. cynomorium ranges from the Bay of Biscay to the Atlantic shores of Africa, including the Mediterranean Sea (Gili and Pagès 1987, López-González et al. 2001, Altuna, 2010).They are involved in the habitat 'the pennatulacean communities of circalittoral and bathyal soft bottoms', according to the description of the biodiversity in the North Atlantic Marine Area (MAGRAMA 2012) with the implementation of the MSFD.
The present work considers the biology and ecology of the most abundant sea pens on the shelf and in the upper bathyal of the southern Bay of Biscay and northwestern Spain (Galicia), mapping their distribution, identifying significant concentrations and describing the benthic/demersal communities associated with them.

Study area
The study area is the northern Spanish shelf in the southern Bay of Biscay (Fig. 1).This area is included in the OSPAR Region IV in the northeast Atlantic and in the ICES Ecoregion Bay of Biscay and the Iberian Coast.It is situated at temperate latitudes and is strongly influenced by the inflow of oceanic water from the Atlantic Ocean.The bottom topography is highly variable, with noteworthy features such as seamounts, banks and submarine canyons, and a highly diverse coastline.The westernmost area (Galician shelf) in particular is fed by many rias and the large Miño River and has a wider shelf, whereas the eastern area (Cantabrian Sea) is cut by many canyons, the shelf is narrower and the slope steep.The former area falls within ICES Division IXa and the latter within Division VIIIc.
In addition, the Le Danois Bank, the Galicia Bank and the Avilés Canyon have been detected as vulnerable ecosystems off the shelf in the study area and proposed as a Marine Protected Area (MPA).

Data source
The study area has been sampled annually for 33 years by a Spanish groundfish survey to collect demersal and benthic megafauna data, oceanography and sediment characteristics, as described by Olaso (1990) and the standard International Bottom Trawl Survey (IBTS) methodology (ICES 2010).These data were used to analyse changes in the distribution and abundance of the pennatulaceans from the families Funiculinidae, Pennatulidae and Veretillidae from 1995 to 2010.
Following the IBTS methodology (ICES 2010), hauls of 30 minutes (ca.126 hauls per survey) were performed at a speed of 3 knots with an otter trawl sampler (baca 44/60) with a 20/mm mesh cod-end liner and a mean swept area of 52504 m 2 .The sampling design used was random stratified to the area, with five geographical sectors and three depth strata, two on the circalittoral (70-120 m, 121-200 m) and one on the upper bathyal (201-500 m) (Fig. 1).In addition, some extra hauls were carried out to cover shallower (<70 m) and deeper (>500 m) grounds.
The environmental variables depth, near-bottom temperature and salinity, location (longitude) and sediment characteristics were analysed to explain differences among the samples.Temperature and salinity were measured with a CTD Seabird 25.Sedimentary characteristics were determined using a sediment collector attached to the groundgear.Particle size of sediments was analysed by combining dry sieving and sedimentation techniques and the variables weight percentage of gravel and coarse sands (GCS: >500 µm), weight percentage of medium, fine and very fine sands (MFS, 63-500 µm) and weight percentage of silt (<63 µm).In addition, weight percentage of organic matter (OM) in the sediment was estimated as weight loss of dried (100°C, 24 h) samples after combustion (500°C, 24 h).

Data analysis
Density (number of colonies per haul) of F. quadrangularis, Pennatula spp., P. spinosum and V. cynomorium were obtained from the 2019 hauls (standard and extra hauls) performed between 1995 and 2010.In addition, records of Pennatula spp.collected on the same surveys but from 2008 and 2014 were used to shed some light on the bathymetric distribution of the two species P. aculeata and P. phosphorea, although they were not included in the density analysis.
The aggregations at or above 10.5 sea pens/haul were studied to analyse faunal assemblages, according to the threshold defined under the MSFD, of species forming habitat when abundances are at or above 200 individuals km -2 or biomass at or above 4 kg km -2 in the same area as our study area, the North Atlantic Marine Area (MAGRAMA 2012).We chose the threshold based on abundance because the biomass of sea pens is not an accurate measurement at sea for these lightweighing species.
An abundance matrix of demersal species was obtained.Highly abundant mobile invertebrates such as the decapod Polybius henslowii, swimming pelagic crustaceans and nektonic cephalopods, namely squids, were removed from the data set.
Density of sea pens per haul was plotted on a histogram and on a spatial context in the area to study the aggregations.Species accumulation curves were represented with the "vegan" package in R software (R Core Team 2014) to ensure that the aggregations were adequately sampled.
Environmental differences between aggregations were studied with constrained correspondence analysis (CCA) and implemented with R (Borcard et al. 2011).Unlike in unconstrained ordination, in canonical ordination the data matrix is forced or constrained by the environmental variables and performs well with species abundance data (Palmer 1993).Forward selection was used for variable reduction and partial CCA to study the effect of correlation and the power of the most significative variables when the other variables were controlled.The significance of the F-statistics associated with all variables was tested using permutation tests, and the most significant explanatory variable was selected (Blanchet et al. 2008).Significance of the adjusted model with the most significant variables was tested using an ANOVA-like permutation test for CCA.CCA results were presented graphically in a tridimensional ordination diagram generated by triplot scaling focusing on inter-species distances, in which samples were represented by symbols, species by points and environmental variables by vectors.In addition, a Mantel correlogram of the Bray-Curtis similarity matrix of the aggregations was constructed to evaluate the spatial correlation as a function of the geographic distance classes.
Differences of the species composition between aggregations were studied with PRIMER version 6.1.13(Clarke and Gorley 2006).Analysis of Bray-Curtis similarity (ANOSIM) was used to detect the statistical significances of these differences, and the similarity of percentages procedure (SIMPER) was used to identify higher-contributing species (at least 80%) to intragroup Bray-Curtis similarity and intergroup dissimilarities.In addition, the frequency and density of the species composition of the aggregations were obtained.Categories of the percentage of occurrence were established according to biotopes classification for Britain and Ireland and categories of density were based on a modified SACFOR scale (semiquantitative abundance scale, Connor et al. 2004), but reduced to four categories (Abundant, >1000/haul; Common, 100-999/haul; Frequent, 10-99/haul; and Occasional, 1-9/haul) due to the low abundances found in the study area.

Distribution of sea pens in the area
Approximately 21% of the hauls were found to contain sea pens every year.The summary of these hauls with sea pens is shown in Table 1.In general, the density was low: around 65% of the hauls with sea pens showed an abundance of one colony per haul.Despite this predominance of low densities, evidence of aggregations was found during the time series (Table 1 and Fig. 2) and throughout the study area (Fig. 3).Approximately 20% of the hauls with sea pens were found to contain densities of sea pens at or above the threshold of 10.5 sea pens per tow.In addition, the Mantel correlogram showed a positive spatial correlation among near samples and negative significance in distant samples, supporting the hypothesis that sea pens form aggregations across space in the area (Supplementary Material Fig. S1).Species accumulation curves for these aggregations (Supplementary Material Fig. S2) indicated a suitable sampling.The results of CCA (Fig. 4A and Table 2) showed that the first axis (Eigenvalue=0.822)splits the aggregations into two groups, depth being the most contributing variable, with higher scores associated with deeper aggregations.The second axis (Eigenvalue=0.467) also splits the sites into two groups, but in this case the most contributing variables were mainly OM and longitude, with negative scores associated with western stations.The two groups of the first axis corresponded to the circalittoral (negative scores) and upper bathyal (positive scores) aggregations, and the two groups of the second axis corresponded to the aggregations on the circalittoral of F. quadrangularis (positive scores) and Pennatula spp.(negative scores).Therefore, three groups were defined in the study area as the results of the CCA: the aggregations of F. quadrangularis on the circalittoral, the aggregations of F. quadrangularis on the upper bathyal and the aggregations of Pennatula spp. on the circalitttoral.Some aggregations on the upper bathyal were slightly deeper than 201 m, so they were at the border between the circalittoral and the upper bathyal.In fact, they were closer to the groups of aggregations on the circalittoral than to the group on the upper bathyal, as is shown in the CCA.
Forward selection showed a significant model (p>0.005) with depth, OM and salinity.Partial CCAs highlighted that depth, and to a lesser extent longitude, were the variables which best defined the groups of aggregations when the other variables were partialled out.
Most of the aggregations (80%) were on the circalittoral between 94 and 196 m, but some aggregations  were also found on the upper bathyal from 194 m to 300 m.On the circalittoral, approximately 73% of the aggregations were dominated by F. quadrangularis, 25% by Pennatula spp., 1% by P. spinosum and 1% by V. cynomorium.On the upper bathyal, all aggregations were dominated by F. quadrangularis (Fig. 3B).
F. quadrangularis was the most abundant sea pen in the study area.This species represented about 83% of the pennatulacean mean abundance caught with an average of 23 colonies per haul.F. quadrangularis was widespread in the sampling area from 84 to 793 m.
Pennatula spp.represented about 14% of the mean abundance of pennatulaceans caught.Despite their low percentages, Pennatula spp. was found every year, with an average of nine colonies per haul.It was mainly found in the Galician area from 91 to 793 m, with two species identified, namely P. aculeata and P. phosphorea.All specimens verified as P. aculeata in the study area occurred deeper than 700 m, while those of P. phosphorea were caught on the shelf (Supplementary Material Table S1).Among the few colonies of P. aculeata identified, four were found at 793 m on the surrounding slope of the Le Danois Bank MPA and, out of the study time series, three colonies were found at 811 m in the southwest of the Le Danois Bank MPA in 2013 and 2014 (unpublished data).
P. spinosum represented a small percentage of the mean abundance of pennatulaceans caught (3%), but it was found every year from a minimum of 1 colony to a maximum of 42 colonies per haul and with an average of 3 colonies per haul.P. spinosum occurred from 50 m to 296 m, although it was mostly found shallower than 150 m.Two aggregations were found, one in Galician waters at a depth of about 48 m and one in the western Cantabrian Sea at about 120 m (Fig. 3B).
V. cynomorium was scarce in the sampling area and it was not found every year.A maximum of 118 colonies per haul was found in 2009 at about 50 m depth in the Galician area, although this species mainly ranged between 90 and 170 m in the easternmost area of Cantabrian Sea and exceptionally occurred around 370 m.Two aggregations were found, one together with the aggregation of P. spinosum in the southern Galician waters and one in the easternmost area of Cantabrian Sea at about 98 m (Fig. 3B).
The density of these two last species may be underestimated in the present study as depths shallower than 70 m were under-sampled.This scarcity prevented us from assessing their ecological importance.

Aggregations of Funiculina quadrangularis on circalittoral and upper bathyal soft bottoms
The aggregations of F. quadrangularis were mainly found on the circalittoral of the central area of the Cantabrian Sea between 95 and 301 m (80%), although, as stated above, some aggregations were also observed on the upper bathyal in the study area (Fig. 3B).
Differences in species composition between aggregations on the circalittoral and the upper bathyal and ecological information of species were obtained by  combining the results of CCA (Fig. 4B), the density data (Supplementary Material Table S2) and SIMPER (Table 3).ANO-SIM identified significant dissimilarities between the samples of F. quadrangularis on the circalittoral and the samples on the upper bathyal (R=0.678;p=0.001), and the SIMPER showed an intergroup dissimilarity between these groups of 93.1%.
The SIMPER results showed that the species Munida sarsi, F. quadrangularis, Macropipus tuberculatus and Sepiola sp.mainly drove the intergroup dissimilarities.In the shallower (circalittoral) aggregations, the species F. quadrangularis and the crinoid Leptometra celtica showed the highest densities, although the highest-contributing species to the similarity were F. quadrangularis, the seastar Astropecten irregularis, the cephalopod Eledone cirrhosa, the hermit crab Pagurus prideaux and L. celtica.In the deeper aggregations, M. sarsi overshadowed the rest of the species due to its high density and contribution to the communities.If we do not cut off the low contributions, we can see that the cephalopod Sepiola sp., F. quadrangularis and the portunid crab M. tuberculatus contributed 4.82%, 3.15% and 1.66%, respectively.Regarding density, M. tuberculatus was the second most abundant species despite its lowest contribution.
In addition, the CCA (Fig. 4B) confirmed the segregation of samples into two groups together with a segregation of species.The eurybathic species were found in both aggregations and located near the CCA centroid, but the positive axis 1 of CCA also showed species with deeper habits and associated with larger sediment (GCS) on upper bathyal soft bottoms, and the negative axis 1 of CCA showed species with shallower habits and associated with fine sediment (silt and MFS) on circalittoral soft bottoms.
The results of the Mantel correlogram (Supplementary Material Fig. S1) also confirmed a variability of species composition along a spatial gradient.The species composition of the aggregations was significantly spatially similar at short distances, whereas when they were spatially separated they did not show significant similarity.

Aggregations of Pennatula spp. on circalittoral soft bottoms
The aggregations of Pennatula spp.were mainly found on the circalittoral in the southern Galician area between 95 and 179 m.No evidence of aggregations was found on upper bathyal (Fig. 3B).
Differences of species composition between the aggregations on the circalittoral of F. quadrangularis and Pennaluta spp.were found (CCA [Fig. 4B], density data [Supplementary Material Table S2] and SIMPER [Table 3]).ANOSIM identified significant dissimilarity between these samples on the circalittoral (R=0.629;p=0.001), and the SIMPER showed an intergroup dissimilarity between groups of 89.4%.The SIMPER showed that the species F. quadrangularis, the crinoid L. celtica, Pennatula spp., the hermit crab P. prideaux, the cephalopod E. cirrhosa, the seastar A. irregularis, the cephalopod Sepiola sp., the anomuran crab Munida intermedia, the portunid crab Liocarcinus depurator, the anomuran crab M. sarsi, and the cephalopods Sepia orbygnyana and Rondeletiola minor mainly drove the intergroup dissimilarities between the aggregations of F. quadrangularis and of Pennatula spp. on the circalittoral in descending order of contribution.In the aggregations of Pennatula spp., the species, E. cirrhosa, A. irregularis and Sepiola sp. were the highest-contributing species to the similarity.

DISCUSSION
Our results provided further evidence of sea pens in the southern Bay of Biscay.Pennatulaceans have been previously reported in the area, not only in the proposed MPAs on the Galicia Bank (IEO 2014a) and in the Avilés Canyon (IEO 2014b), but also in the inner basin of the Le Danois Bank MPA (BOE 2011) and in Galician waters (Aguilar et al. 2009).Moreover, under the MSFD, Spain has recently identified pennatulacean communities as the most diverse and richest habitats on circalittoral and bathyal soft bottoms of the northern Atlantic Spanish area (MAGRAMA 2012).Our results also provided further evidence of sea pens forming aggregations on the northern Atlantic Spanish shelf.Although sea pens were broadly found at low density in the study area, the aggregations with a noteworthy high density in the central part of the Cantabrian Sea and in the southern Galician area and the high spatial autocorrelation among near samples supports the hypothesis of sea pens forming aggregations in the area.The hotspots in the Cantabrian Sea were of the species F. quadrangularis, while in the Galician area they were of Pennatula spp.Though the sea pens V. cynomorium and P. spinosum were scarce, they showed two hotspots in the southern Galician area and in the easternmost area of the Cantabrian Sea, respectively.
In addition, differences in the distribution of sea pen aggregations were supported by the fact that longitude was among the variables contributing most to the differentiation of the aggregations.This factor is a result of environmental differences between the western (Galicia area) and eastern (Cantabrian Sea) zones of the study area, the Galician and Cantabrian Sea, respectively.In general, there is higher primary production in Galician waters (López-Jamar et al. 1992, OSPAR 2000) and, in particular, a heterogeneous grain size and a higher OM content have been reported to the south of Cape Finisterre (López-Jamar et al. 1992), where Pennatula spp. was highly abundant.Pennatula spp.are known to be much stouter and fleshier than F. quadrangularis (Manuel 1981), which probably helps them to settle in less homogeneous sediments.In fact, both species are known to occur in muddy substrata, but Manuel (1981) also reported P. phosphorea in sandy substrata and Greathead et al. (2015) reported the absence of F. quadrangularis in sediment with gravel content above 30%.
In addition, longitude also means differences in topography, which leads to an apparent absence of aggregations of sea pens on the narrow shelf and steep slope of the eastern part of the Cantabrian Sea, with the exception of the easternmost shelf area.Therein, mud deposits (Jouanneau et al. 2008) could be supporting the settlement of the few aggregations of F. quadrangularis.
On the other hand, depth was also an influential factor in the distribution of the sea pens in the study area.The few aggregations of P. spinosum and V. cynomorium were found at about 50 m depth, although some specimens were recorded at unusual depths in the area if compared with data in the literature (Supplementary Material Table S1).We found that 70% of the occurrences of V. cynomorium in our study area were detected above 100 m and three colonies were even found at an unusual depth of 372 m, despite the fact that these two species are known to be continental-shelf species with a distribution related to terrigenous sediments (Gili and Pagès 1987, Williams 2011, Sardá et al. 2012).López-González et al. (2001) reported, in the eastern Atlantic, a wider depth range (25-188 m) than Williams (1990, 13-91 m), and Gili and Pagès (1987) showed an even deeper limit in the Mediterranean Sea (around 260 m).However, neither of them mentioned such a deep distribution for the species.In addition to this interesting result, the highest aggregation of colonies of V. cynomorium (118 colonies), together with some colonies of P. spinosum (19 colonies), were found near the Arousa Ría (Galicia, Spain), where fields of V. cynomorium have been previously reported (Fernández et al. 2005, Aguilar et al. 2009).This result provides further evidence of sea pens forming aggregations, widely recognized as key benthic habitats (OSPAR 2008, DFO 2010, Kenchington et al. 2011).
Another result of the influence of depth on the distribution of sea pen species was found in the Galician area, where Pennatula spp. was abundant from 100 to 200 m, while the few occurrences of F. quadrangularis were mainly between 200 and 300 m.Hughes (1998) and Greathead et al. (2007) also suggest the preference of F. quadrangularis for deeper habitats where the conditions may be more stable.In addition, also relating depth to sediment particle size and organic content, finer and organic-rich sediments have been reported on the inner shelf off the Rías Bajas (Galicia, Spain), where Pennatula spp. was more abundant, while on the outer shelf the sediments are sandier and have low organic content (López-Jamar et al. 1992).Nevertheless, Greathead et al. (2007) suggest that depth is a more important factor than sediment type in the distribution of F. quadrangularis.
Tracing the distributions and bathymetric ranges of P. phosphorea and P. aculeata in the Bay of Biscay and nearby areas from the literature is not easy, as we highlighted in the Introduction.However, our results provide evidence of different bathymetric distributions between the two species.Hence, the fact that P. aculeata extends deeper than P. phosphorea suggests that the aggregations of Pennatula spp.found in the study area do not include specimens of P. aculeata but only of P. phosphorea, although more data are still sorely needed.P. aculeata has also been found at similar depths in the southern Bay of Biscay, as a result of the bathyal surveys carried out by the ECOMARG project (www.ecomarg.com).During these surveys, P. aculeata was found only at 745 and 761 m in the Le Danois Bank MPA, at 786 m on the Galicia Bank (northwestern Spain) and at 828 m in the Avilés Canyon (northern Spain) (A.Altuna, personal communication).All these records of upper bathyal presence of P. aculeata, and to a large extent, the lack of reliable data of P. phosphorea from these areas, highlight the need for more research on the bathymetric distribution and ecological needs of these two species.
With respect to the sampling of sea pen communities in the southern Bay of Biscay, the present study adds long-term information on the distribution of sea pens in an extensive area, despite the use of trawling sampling.Techniques such as beam trawl and demersal otter trawl are often used in benthic community studies, although monitoring is extensively recommended for the management of sea pens and burrowing megafauna (Hughes 1998, Greathead et al. 2007).Monitoring and assessment of vulnerable habitats have been carried out in the area with non-destructive monitoring techniques such as remotely-operated vehicles and photogrammetric sledges (Aguilar et al. 2009, www.ecomarg.com, Sánchez et al. 2009), but they were short-term or spatially isolated studies.In addition, sampling of some sea pen communities is not easy and very inaccurate, as some pennatulaceans can withdraw below the sediment surface and bias the results of monitoring.The density of some such sea pens in a given area can change by over 200% in a few hours (Birkeland 1974).Virgularia mirabilis is known to withdraw into the sediment and P. phosphorea and P. aculeata are also able to do so, although not as quickly, but F. quadrangularis, P. spinosum and V. cynomorium do not show this feature (Langton et al. 1990, Hughes 1998).
The species composition of the aggregations of sea pens in the study area agrees with previous studies in the area, which showed depth as a key factor structuring communities (Olaso 1990, Serrano et al. 2006).Many species that occur in the aggregations of F. quadrangularis and Pennatula spp.coincided with the eurybathic species described in the area by Olaso (1990), such as the portunid crab Liocarcinus depurator, the spider crab Macropodia longipes and the sea snail Neptunea contraria, or with the characteristic species of the shelf (between 140 and 300 m), such as the sea anemone Actinauge richardi, the sea snail Ranella olearium and the sea cucumber Parastichopus regalis.Specifically in the shallower samples of F. quadrangularis, the higher-contributing species Astropecten irregularis, Eledone cirrhosa, Pagurus prideaux and Leptometra celtica were reported by Olaso (1990) on the upper edge of the shelf (between 90 and 140 m), whereas the anomuran Munida sarsi, as the most abundant and the highest-contributing species in the deeper samples of F. quadrangularis, was reported as a dominant species on the shelf.This author compared the different distribution of M. sarsi with the other common anomuran in the area, Munida intermedia and described both species as dominant species on the shelf, but M. sarsi showed greater density at 275 m in the Cantabrian Sea, and M. intermedia showed greater density at 175 m in the Galician area.This explains the higher abundance of M. sarsi in the deeper aggregations of sea pens in the area and the greater occurrence of M. intermedia in the aggregations of Pennatula spp., which are mainly distributed in Galician waters.In addition, Olaso (1990) and Serrano et al. (2006) showed depth and sediment characteristics to be priority factors determining species composition.
Large differences were found between the aggregations described in the present study and the habitat "sea pen and burrowing megafauna communities" defined by OSPAR.This habitat was reported in the North Sea, the Irish Sea and the Bay of Biscay and may include populations of V. mirabilis, P. phosphorea or F. quadrangularis with burrowing crustaceans such as Nephrops norvegicus, in fine mud at depths of between 15 and 200 m.In the southern Bay of Biscay, we found aggregations of F. quadrangularis on circalittoral soft bottoms but also on the upper bathyal deeper than 200 m.In addition, we also found aggregations of Pennatula spp. in the Galician area, all of them on circalittoral bottoms.Finally, V. mirabilis was absent and N. norvegicus was found at low densities in the study area, despite its general distribution between 90 and 600 m in the study area (ICES 2012).Other species not mentioned by OSPAR, M. sarsi, A. irregularis, L. depurator, Ophiura ophiura and M. intermedia, were frequent in the aggregations found in the study area.Although the southerner area and the deeper and trawling sampling in the present study marked differences from the methodology used to define this habitat in Northern Europe, and the community composition and density of the aggregations of sea pens on the northern Spanish shelf might be biased by some other factors.
Some of these frequent species, L. depurator, O. ophiura and A. irregularis, are known to be opportunistic species in trawled areas (Kaiser et al. 1998, Vergnon and Blanchard 2006, de Juan et al. 2007) like our study area which support an intensive fishing pressure (Sánchez andOlaso 2004, ICES 2012).In contrast to the North Sea, the Irish Sea and northern Bay of Biscay, where Nephrops is one of the most important fisheries, in the southern Bay of Biscay N. norvegicus has stopped being a target species over time due to its low abundance.Nevertheless, it is a valued component of the mixed bottom trawl fishery together with hake, anglerfish, megrim, horse mackerel, mackerel and blue whiting (Fariña and Herraiz 2003, Punzón et al. 2010, ICES 2012).The Nephrops stock is declining in the study area (Fariña and Herraiz 2003) and recovery plans were adopted in 2005 (Council Regulation (EC) 2166/2005).This decline in Nephrops stock could have coincided with a decline in sea pens in the past, which could partially explain the low densities of sea pens found in the study.Now, this mixed bottom trawling could be a similar threat for sea pens in the Cantabrian and Galician waters to that of Nephrops trawling in the North Sea, the Irish Sea and the northern Bay of Biscay, but this possibility has not been assessed.Evidences like the intensive fishing pressure, the declining Nephrops stock, the presence of opportunistic species in the aggregations and the low densities of sea pens found in the present study suggest that trawling might be affecting sea pen abundance in the southern Bay of Biscay and highlight the need for further information to assess the status of these important communities in this region.
the Cantabrian Sea and off Galicia and to the crew of the R/V Cornide de Saavedra.We also thank the ECO-MARG and INDEMARES projects for the information on benthic communities in the Le Danois Bank MPA, the Galicia Bank and the Avilés Canyon.
The manuscript was improved thanks to the careful revision and constructive comments of Ellen Kenchington.

Fig
Fig. 3. -Distribution of sea pens (A) and of the aggregations (B) in the study area.

Fig. 4 .
Fig. 4. -A, CCA of the aggregations (F.quadrangularis [dots] and Pennatula spp.[squares]) based on the species composition.The samples on the circalittoral in black and the samples on the upper bathyal in grey.B, CCA of the species composition of these aggregations.Species abbreviations appear in Supplementary Material Table S2.The arrows indicate the direction and magnitude of the environmental variables: depth, longitude, temperature, salinity, median particle diameter (GCS, weight percentage of gravel and coarse sands [>500 µm]; MFS, weight percentage of medium, fine and very fine sands [63-500 µm]; Silt, percentage weight of silt [<63 µm]) and weight percentage of organic matter (OM).

Table 1 .
-Summary of the collection data from 1995 to 2010.N, number of hauls analysed; Nseapens, number of hauls with seapens; Naggreg.C, number of hauls with densities at or above the threshold 10.5 sea pens/haul on the circalittoral; Naggreg.U, number of hauls with densities at or above the threshold 10.5 sea pens/haul on upper bathyal.

Table 2 .
-Canonical coefficients for each environmental variable on each canonical axis (CCA) on the aggregations of F. quadrangularis and Pennatula spp.GCS, weight percentage of gravel and coarse sands (>500 µm); MFS, weight percentage of medium, fine and very fine sands (63-500 µm); Silt, percentage weight of silt (<63 µm); OM, weight percentage of organic matter.

Table 3 .
-Results of the SIMPER analysis of the aggregations of F. quadrangularis on the circalittoral (1) and upper bathyal (3) and of Pennatula spp. on the circalittoral (2).Av abund, average abundance; Contrib.% and Cum.%, contribution and cumulative percentage respectively of the higher-contributing species.