Holothurians are nutrient recyclers, bioturbators, and prey for other taxa in diverse marine ecosystems; thus, they directly impact the distribution of organisms with which they are associated. However, little is known about annual changes in holothurian recruitment. Islas Marietas harbors a coral community with a rich assemblage of holothurians. From March 2011 to February 2012, bimonthly samplings were performed in order to identify the annual recruitment patterns of species Pentamera chierchia and Neothyone gibber. Density (indcm−2) was assessed for each of the 2 species, as was the micro-substrate abutting them. Both species had highest densities in the May and June period (N. gibber=0.0903indcm−2; P. chierchia=0.0396indcm−2); densities<0.01indcm−2 were documented over the remainder of the year. This is the first record of P. chierchia in the study area. Furthermore, no correlation between species density and substrate type was found (s≤0.61, n=7, p>0.05). The study of seasonal recruitment fluctuations of reef-associated holothurians and their relationship with factors such as reproduction and availability of coralline substrate may provide information on their life cycles, in order to initially understand some ecological patterns of these species that contribute to the maintenance of regional biodiversity.
Los holoturoideos cumplen funciones como recicladores de nutrientes, bioturbadores y presas de otros taxones en diversos ecosistemas marinos; por lo que afectan directamente la distribución de los organismos asociados a ellos. Pese a su importancia funcional, poco se sabe acerca de sus patrones anuales de reclutamiento. Las Islas Marietas albergan una comunidad coralina con un rico ensamblaje de holoturoideos. Con el fin de identificar los patrones de reclutamiento anuales de las especies Pentamera chierchia y Neothyone gibber, de marzo 2011 a febrero 2012 se realizaron 6 muestreos bimensuales. Se evaluó la densidad (indcm−2) para cada especie y el micro-sustrato asociado. Ambas presentaron las mayores densidades durante mayo-junio (N. gibber=0.0903 ind cm−2; P. chierchia=0.0396 ind cm−2); durante el resto del año la densidad fue<0.01 ind cm−2. Este es el primer registro de P. chierchia en el área de estudio. No hubo correlación entre densidad y micro-sustrato (s ≤ 0.61, n=7, p>0.05). El estudio de las fluctuaciones estacionales del reclutamiento de holoturias y su relación con factores como reproducción y disponibilidad de sustrato puede proporcionar información sobre sus ciclos de vida, dando mayor entendimiento a los patrones ecológicos de estas especies que contribuyen al mantenimiento de la biodiversidad regional.
Recruitment is an important process that ultimately determines the population structure of most marine invertebrate species (Caley et al., 1996; Pawlik, 1992; Rodríguez, Ojeda, & Inestrosa, 1993). Patterns of marine invertebrate reproduction and recruitment are regulated by abiotic factors, including seawater temperature, nutrient concentrations, and pH, as well as by biotic factors such as coralline substrate availability, presence of predators, and population density (Abelson & Denny, 1997; Broitman et al., 2008; Hadfield, 1986).
Coral communities are productive ecosystems that harbor a high biodiversity of both vertebrates and invertebrates (Glynn & Enochs, 2011; Sheppard, Davy, & Pilling, 2010). Holothurians represent one highly diverse and conspicuous group of such coral reef inhabitants, and they rely on reefs for food, shelter, and reproduction (Borrero-Pérez, Benavides-Serrato, & Díaz-Sánchez, 2012; Miller & Pawson, 1984; Uthicke, 1999). The Class Holothuroidea is represented by 6 orders: Apodida, Elasipodida, Aspidochirotida, Molpadiida, Dendrochirotida, and Dactylochirotida (Kerr, 2000; Pawson, 2007); all of them can recruit and inhabit a variety of substrates, such as mud, sand, rocks, algal beds, coral rubble, dead coral, or live coral (Alvarado & Solís-Marín, 2013; Solís-Marín, Arriaga-Ochoa, Laguarda-Figueras, Frontana-Uribe, & Durán-González, 2009). The order Dendrochirotida is particularly well represented in coral reef ecosystems, although these sedentary, cryptic organisms can also be found in temperate areas (Bakus, 1973; Birkeland, 1989). Holothurians are an essential component of the diet of fishes, crustaceans, polychaetes, and seastars (Birkeland, 1989; Francour, 1997; Jennings & Hunt, 2010).
Islas Marietas National Park (IMNP) represents one of the most important coral communities in the Central Mexican Pacific (CMP) region (Cupul-Magaña, Aranda, Medina-Rosas, & Vizcaíno, 2000). The benthic structure is characterized by the presence of branching and massive corals, with the genus Pocillopora demonstrating the highest cover amongst scleractinians (Conanp, 2007; Reyes-Bonilla, Carriquiry, Leyte-Morales, & Cupul-Magaña, 2002). Until present, a total of 38 species of holothurians have been found in reef ecosystems of the Eastern Pacific, and 13 of these are dendrochirotids (Alvarado & Solís-Marín, 2013; Alvarado, Solís-Marín, & Ahearn, 2010; Bastida-Zavala et al., 2013; Deichmann,1941, 1958; Maluf, 1987; Ríos-Jara, Galván-Villa, & Solís-Marín, 2008; Ríos-Jara, Galván-Villa, Rodríguez-Zaragoza et al., 2013; Solís-Marín et al., 2009). At IMNP, 3 species of dendrochirotids have been reported: Pseudocnus californicus, Cucumaria flamma, and Neothyone gibber (Conanp, 2007; Nepote-González, 1998; Solís-Marín & Laguarda-Figueras, 1999).
Despite the importance and prevalence of sea cucumbers in the CMP area, most prior studies have only focused on taxonomic descriptions (e.g., Caso, 1962; Conanp, 2007; Chávez-Dagostino, Nepote-González, Medina- Rosas, & Solís-Marín, 2000; Deichmann, 1938, 1941; Nepote-González, 1998; Solís-Marín et al., 2009; Théel, 1886) or on the basic ecology of commercially valuable species (e.g., Herrero-Pérezrul & Reyes-Bonilla, 2008). In contrast, few data exist about the reproductive cycle (Herrero-Pérezrul, Bonilla, García-Domínguez, & Cintra-Buenrostro, 1999), and limited studies of their recruitment patterns have been undertaken (e.g., Aguillón-Negreros, 2011). This knowledge dearth can be attributed to the low economic value of most of the species, or to their small size, cryptic behavior, and low densities (Alvarado & Solís-Marín, 2013). Indeed, information on invertebrate recruitment at IMNP is limited (Cruz-García, Cupul-Magaña, Hendrickx, & Rodríguez-Troncoso, 2013; Medina-Rosas, Carriquiry, & Cupul-Magaña, 2005) and is absent in the case of holothurians; only a few reports of the early life history stages of these organisms in the Mexican Pacific can be found (Aguillón-Negreros, 2011; Alvarado & Solís-Marín, 2013).
The study of recruitment patterns in holothurians is important for the understanding of population structure at adult stages (Balch & Scheibling, 2000; Caley et al., 1996); besides, this process influences the population dynamics of other species that are ecologically linked to sea cucumbers (Bakus, 1973; Birkeland, 1989; Glynn & Enochs, 2011). Coral communities provide optimal biotic and abiotic conditions for the recruitment of many sea cucumber species (Alvarado & Solís-Marín, 2013). In this study, the recruitment of Pentamera chierchia and N. gibber was monitored over an annual timescale in IMNP, and we also determined the influence of factors such as seasonal changes in sea surface temperature and coverage by functional groups in the benthic substrate on the recruitment pattern of both species.
Materials and methodsIslas Marietas is located at the northeastern end of Banderas Bay, México (20°41′56″–20°41′53″N, 105°35′02″–105°34′59″ W; Fig. 1). It consists of 2 islands of volcanic origin; with an important coral community composed by patches of Pocillopora, Porites, and Pavona, as well as sandy and rocky bottom areas (Conanp, 2007). The region is characterized by a sub-tropical climate with predominantly southeasterly winds during the day (García, 1998). The average sea surface temperature (SST) is 26.4°C (annual minimum=23°C and maximum=30°C) and is influenced by the California Current (CC), in which cold, low salinity waters are brought in from November to April. Then, the Costa Rica Coastal Current (CCR) brings in warm, high salinity water from May to October. As such, Banderas Bay is an important transitional oceanographic area (Conanp, 2007; Espino-Barr, Cabral-Solís, García-Boa, & Puente-Gómez, 2004; Ulloa-Ramírez et al., 2008; Wyrtki, 1965).
Sampling was performed as reported by Cruz-García et al. (2013) and Sotelo-Casas, Cupul-Magaña, and Rodríguez-Troncoso (2014) from March 2011 to February 2012 at one reef patch located at Larga Island. Briefly, in order to generate a substrate for holothurian recruitment, 144 dead Pocillopora fragments ∼5cm lengths were used. To eliminate any organic matter, all fragments were submerged in 10% sodium hypochlorite (in filtered seawater) for 24h; afterwards, the fragments were washed with freshwater for 8h and dried in an oven (Precision Scientific) at 60°C for 24h. Before placing fragments, 2 semi spherical concrete structures with 6 steel rods each were placed on a sandy bottom, located at ∼5m depth next to the healthy Pocillopora coral colonies. Initially, a total of 24 fragments were attached to both structures; every 2 months the fragments were removed and replaced by new ones. Upon collection, fragments were placed individually in plastic bags in situ, transported to the laboratory, and preserved at −20°C until further processing. Each fragment was examined, and all the holothurians were separated, identified, and counted using a stereoscopic microscope (Zeiss Stemi DR 1040) or a compound microscope (LABO), following the taxonomic criteria described by Solís-Marín et al. (2009) and Prieto-Ríos (2010). The total length of each specimen was measured and individuals were dissected to locate their reproductive structures. Specimens and spicules were photo-documented using an Olympus camera and finally preserved in ethanol 70%. Additionally, the micro-substrate adhering to the coral fragments was separated and classified into 4 functional groups: algal turf, crustose algae, bryozoans, and macroalgae, dried individually in a Precision Scientific oven at 60°C during 24h, and weighed using an analytical balance (Explorer Ohaus). Finally, the area of each coral fragment was determined using the aluminum foil method of Marsh (1970). The abundance of each species was expressed as the number of individuals per area (indcm−2), and the average mass percentage of each type of micro-substrate was calculated as: (dry weight of each substrate/dry weight of total substrate)*100.
The sea surface temperature (SST) was recorded in situ every 15-min during the 12 months of experiment using a HOBO® Pendant data logger located next to the artificial structures. Temperature data was expressed as mean temperature±standard error (°C±SE). Statistical analyses were performed using Sigmaplot (v. 11.0) software for Windows. All data were tested for normality and homogeneity of variance. As all data were not normally distributed (Shapiro–Wilk, ¿=0.05), a Kruskal–Wallis (non-parametric, one-way ANOVA) test was calculated to determine differences between sampling times of each species; for pairwise comparisons of densities, Dunn's post hoc tests (¿=0.05) were performed. We assess the relation between the density of recruitment of each species with season, micro-substrate cover, and temperature using Spearman's correlations (¿=0.05).
ResultsA total of 163 sea cucumber recruits were identified: 84 N. gibber and 79 P. chierchia. All the organisms were classified as juveniles based on the lack of reproductive structures. Over the course of the year, the average densities of N. gibber and P. chierchia were 0.0156 and 0.0114indcm−2, respectively. For both species, higher densities were documented between May and June (N. gibber=0.0903indcm−2; P. chierchia=0.0396indcm−2). Mean annual SST was 26.1±0.02°C. The minimum average temperature was recorded during March–April (20.9±0.01°C), and the maxima were register during September–October with values of 30.5±0.01°C. The maximum variation between sampling bimonthly was from March–April to May–June with an increase of ∼6°C (Fig. 2).
The holothurians showed different densities of recruitment over the course of the year (Kruskal–Wallis [effect of time] for P. chierchia p<0.001, Dunn's test, p<0.05; Kruskal–Wallis [effect of time] for N. gibber p<0.001, Dunn's test, p<0.05). Recruitment was significantly higher in both species during May–June than the rest of the year, coinciding with an important rise of SST (Fig. 2). Except for the recruitment peak observed during May–June, both P. chierchia (from 0.00092 to 0.00839indcm−2) and N. gibber (from 0 to 0.0129indcm−2) present low recruitment densities throughout the year. There was a total absence of recruits of N. gibber between November and February, coincident with the time at which lowest SSTs were registered (Fig. 2). Despite this apparent pattern, Spearman's test revealed no significant correlation between temperature and the density of either species (s≤0.58, n=6, p>0.05).
A change in micro-substrata was observed over the year, with a negative correlation between the proportion of cover of algal turf and crustose algae (s=−0.96, n=6, p<0.001). There was a decrease in the cover on algal turf, and the consequent increase of crustose algae during the warm season, which persisted until the temperate season during January–February (Fig. 3). Macro-algae and bryozoans substrates were found from March to June; however, they declined in the July–August sampling period and decreased to zero from September to February. There was no significant correlation between the density of recruits and the substrate (s≤0.61, n=6, p>0.05).
DiscussionIMNP harbors a biodiverse and structure-rich coral community. The substrate coverage of this ecosystem is not only comprised of reef-building corals, but also of algae, rock, and sand; all such substrates could permit settlement by marine invertebrates (Birkeland, 1989; Cupul-Magaña et al., 2000; Chávez-Dagostino et al., 2000; Graham & Nash, 2013). Coral reef holothurians have the ability to colonize diverse micro-habitats in the reef structure, such as crevices and holes, soft bottoms, empty shells, and burrows of other organisms (Birkeland, 1989; Glynn & Enoch, 2011). However, different species of holothurians demonstrate unique preferences for recruitment locations (Mercier, Battaglene, & Hamel, 2000; Williams & García-Sais, 2010), and the factors that promote the recruitment processes have not been previously studied at IMNP. Previous reports in the Mexican Pacific found that P. chierchia and N. gibber tend to recruit to calcareous substrates, such as rocks or dead coral (Alvarado & Solís-Marín, 2013; Deichmann, 1941, 1958; Solís-Marín et al., 2009). As such, the study area provided widespread areas for recruitment of both species, the former of which had never been documented at IMNP.
Although, adults of both species were rarely seen during the majority of the year due to their cryptic nature during the day, the 2 species are considered to have a wide distribution and have been reported from several localities from Gulf of California to Peru (Deichmann, 1941; Prieto-Ríos, 2010; Solís-Marín et al., 2009). The high densities of P. chierchia and N. gibber recruits at IMNP could be an indication of the role of the island as a sea cucumber nursery. The importance to identify and conserve nursery sites has been emphasized in previous studies on the Eastern Pacific (Aguillón-Negreros, 2011; Toral-Granda, 2005); and therefore the Island becomes a priority area for biodiversity conservation in the CMP. In addition, the presence of P. chierchia in the study area fills a critical distribution gap for this species along the CMP and represents an important advance in the knowledge of the biodiversity in the region (Honey-Escandón, Solís-Marín, & Laguarda-Figueras, 2008; Salcedo-Martínez, Green, Gamboa-Contreras, & Gómez, 1988; Solís-Marín, Laguarda-Figueras, Durán-González, Ahearn, & Vega, 2005; Toral-Granda, 2005).
Annual temperature fluctuations regulate certain physiological activities of invertebrates, such as reproduction and recruitment, which are triggered and favored during the warm periods (Cowen & Sponaugle, 2009; Somero, 2005). The timing of such events is determined by the local oceanographic conditions (Broitman et al., 2008; Pineda, Reyns, & Starczak, 2009; Roughgarden, Gaines, & Possingham, 1988). Some species are characterized by a single spawning period during summer, which is also driven by SST changes (Herrero-Pérezrul et al., 1999). In the study area, the recruitment density of holothurians presented seasonal changes, with peaks during May–June, and low values for both species for the rest of the year. As the rate of recruitment depends on larval supply, the high densities of recruits observed in May–June may be due to a recent spawning event that occurred after the March–April sampling time; such an event may have been associated with the strong increase in SST between March–April and May–June. The lack of correlation between SST and recruitment has been previously reported for echinoids, asteroids, and holothuroids (Miller, 1995), and only associated as a determinant factor for spawning during the spring season and subsequent recruitment in summer in species such as Isostichopus fuscus and Holothuria (Selenkothuria) lubrica distributed at the Mexican Pacific (Herrero-Pérezrul et al., 1999; Skarbnik-López, Herrero-Pérezrul, Reyes-Bonilla, Domínguez, & Morales, 2010).
Although both temperature and substrate type are important for recruitment in most invertebrates, oceanic currents, chemical signals, depth, and bottom complexity all play important roles (Abelson & Denny, 1997; Allani & Melloni, 1999; Broitman et al., 2008). In addition, biotic factors, such as predation, substrate preference, and food supply and endogenous factors, such as larval development strategy, duration of the competent stage, and motility (Allani & Melloni, 1999; Hadfield, 1986; Pawlik, 1992; Poulin, Boletzky, & Féral, 2001), have important implications for recruitment and may have accounted for some of the patterns observed at IMNP. Here, it is important to consider that the holothurians in this study settled on coral skeletal structures, suggesting this as an adequate substrate for recruitment, independent of the micro-substrate growing on top. The affinity of the sea cucumbers to this substratum can be assumed as a positive response to the available refuge, positive chemical signals or structural complexity that may be beneficial for larval survival (Poulin et al., 2001). In addition, the data show that the presence of these dendrochirotids did not depend on a specific micro-substrate, possibly due to their filter-feeding habits (Birkeland, 1989; Solís-Marín et al., 2009).
After recruitment and the initial growth phase, most sea cucumbers migrate from the site of first settlement to more secure regions, such as soft substrate, burrows, or crevices in the rocks or other cryptic environments, where there is enough food and adequate shelter from predators (Alvarado & Solís-Marín, 2013; Balch & Scheibling, 2000; Young & Chia, 1982). This cryptic behavior may explain the low densities of juveniles found after to May–June and the absence of adult holothurians in our samplings.
The data on recruitment described for P. chierchia and N. gibber provides the first insight into the understanding of recruitment dynamics of subtropical holothurians and a non-random preference for certain types of habitats as observed on other holothurids (Aguillón-Negreros, 2011; Wolff, Schuhbauer, & Castrejón, 2012). The presence of recruits of holothurids associated to the coral community gives an indirect evidence that these organisms recognize appropriate environmental conditions to reproduce and settle within the site. Furthermore, the recruitment process maintains not only these 2 holothurid populations but the organisms that depend on them. The evaluation of the recruitment pattern of both species represents an increase of the knowledge of sea cucumbers in the study area and evidence of the potential of IMNP as a source of feeding, reproduction, refuge and breeding for marine invertebrates in the region.
RCSC was supported by Conacyt grant No. 291281, project PROMEP 220265 for APRT and PIFI P/PIFI-2010-14MSU0010Z-10 for ALCM. The authors thank the authorities of the IMNP (Conanp) for assistance and use of facilities during the sampling periods. In addition, the authors would like to thank students who contributed to the sample processing, as well as Dr. Anderson Mayfield, who carefully proofread the English.
Peer Review under the responsibility of Universidad Nacional Autónoma de México.