Original Research ARTICLE
A Multidisciplinary Experimental Study on the Effects of Breeders Diet on Newborn Seahorses (Hippocampus guttulatus)
- 1Department of Ecology and Marine Resources, Instituto de Investigaciones Marinas (CSIC), Vigo, Spain
- 2Dipartimento di Scienze della Vita e dell’Ambiente, Università Politecnica delle Marche, Ancona, Italy
- 3Marine Products Chemistry, Instituto de Investigaciones Marinas (CSIC), Vigo, Spain
- 4Dipartimento di Scienze Veterinarie, Università di Messina, Messina, Italy
The global trade in seahorses is increasing and new rearing techniques are currently available for a few species. One of the main bottlenecks in seahorse production is reproductive success, which is dependent upon the administration of high quality diets to breeders. However, little is known about the most adequate diets, on how they should be administered and on how they might affect newborn performance. In this study, three homogeneous groups of adult seahorses Hippocampus guttulatus were maintained in captivity. Each group was fed on one of the following diets since the onset of the breeding period: unenriched adult Artemia (diet A), unenriched adult Artemia and mysidaceans (diet AM) and mysidaceans (diet M). As a positive control, a wild pregnant male was used. The mixed diet AM provided the highest overall breeding success and reasonable fatty acid profiles, and better resembled newborn from a wild male. The use of mysidaceans solely (Diet M) provided the worse results in spite of the high DHA content. Significant differences were observed in newborn characteristics and fatty acid profiles among dietary groups and along the breeding season, especially in total n-3 HUFA and DHA. Monospecific diets including Artemia or mysidaceans very likely suffered nutritional deficiencies, which were reflected in morphological alterations of the muscle tissue (diet A) or abnormal early mortalities in newborn (diet M). Three progressive stages were identified along the whole breeding period: initial mixed capital-income period (100–120 days since the onset of the breeding period), followed by an inflection period (short transitional stage – income sources), and finally a long period characterized by the use of income sources and the progressive exhaustion of body reserves. Special attention should be deserved to the progressive changes in fatty acid profiles along the breeding season, resulting in a decrease in newborn performance. Considering the time required for a diet to be reflected in newborn, we recommend Hippocampus guttulatus breeders to be fed on a high quality diet for at least 3 months prior the breeding period.
The global trade in seahorses is estimated to exceed 20 million animals per year to sustain the aquarium trade, the traditional Chinese medicine and as souvenirs (Vincent, 1996; Planas et al., 2017a). In addition, because of the destruction and degradation of their coastal habitats (seagrasses, coral reefs, and mangroves), wild populations are progressively declining (Vincent, 1996; Martin-Smith et al., 2004). At present, many are the concerns over their long-term viability in nature (Vincent, 1996; Vincent et al., 2011; Kumaravel et al., 2012) and all seahorses species are included in Appendix II list of endangered species by CITES (Convention for the International Trade in Endangered Species).
Aquaculture represents a valid alternative to wild captured seahorses but ex-situ production still needs to face several biological and technical challenges (Cohen et al., 2016; Planas et al., 2017a). Both, low spawning quality and juvenile survival rates are the most critical factors for optimizing commercial seahorse aquaculture (Adams et al., 2001; Chang and Southgate, 2001; Lin et al., 2007).
Reproduction represents a key aspect in breeding success, and broodstock feeding/nutrition is doubtless, one of the key factors to sustain reproduction in teleost (Migaud et al., 2013). The dietary effect on the output performance (egg and fry quality and quantity) is well known in fishes, and many studies have reported the importance of lipids and fatty acids (FAs) (Rainuzzo, 1993; Bromage, 1995; Izquierdo et al., 2001; Sargent et al., 2002; Tocher, 2003, 2010; Pavlov et al., 2004). Broodstock preparation is thus a critical step in the success of the spawning process of a certain species. On this regard, many studies have been directed toward egg quality in oviparous fishes and the underlying characteristics to produce viable fry (Kjørsvik et al., 1990; Brooks et al., 1997). The egg’s potential to produce viable fry is determined by several physical, genetic, and chemical factors, including the initial physiological processes occurring in the egg itself (Kjørsvik et al., 1990).
The ongoing process of oocytes maturation is one of the most metabolic demanding activities in fish and the quality of parental nutrition is thus of primary importance for the appropriate allocation of macromolecules into the oocytes and their maturation (Izquierdo et al., 2001; Volkoff and London, 2018). However, there is a scarcity of studies focused on breeding performance and egg or fry quality in fishes with male pregnancy such as syngnathids (pipefish and seahorses). Due to the uncommon nature of reproduction, syngnathids offer a unique opportunity to examine pathways of repercussion of breeder’s dietary resources on breeding performance and the resulting newborn quality. Reproduction in seahorses has received increasing attention in recent years, mostly on the effect of environment (e.g., photoperiod regime and temperature) (Lin et al., 2006, 2008b; Planas et al., 2010, 2013) zootechnical conditions (e.g., sex ratio, aquaria design, and seahorse size) (Woods, 2000; Dzyuba et al., 2006; Faleiro et al., 2008, 2016; Planas et al., 2008) and biological/physiological features (Boisseau, 1967; Carcupino et al., 2002; Poortenaar et al., 2004; Stölting and Wilson, 2007; Scobell and MacKenzie, 2011; Zhang et al., 2019; Wittington and Friesen, 2020). Only a few studies have focused on the effects of feeding and dietary sources on breeding success and newborn features in seahorses under ex-situ conditions (Lin et al., 2007; Faleiro and Narciso, 2010; Binh and Serrano, 2012; Otero-Ferrer et al., 2012, 2016, 2020; Palma et al., 2012, 2017; Saavedra et al., 2015).
Seahorses (Hippocampus spp.) are batch spawners with cyclical egg maturation and repeated mates within the breeding season. The eggs produced by females are transferred to the brood pouch of males, where embryogenesis will take place for 2–4 weeks, depending on the species and temperature (Foster and Vincent, 2004). The sealed brood pouch of the male works as a pseudo-placenta, not only providing protection and oxygen to the developing embryos, but also nutrients (Wilson et al., 2001; Carcupino et al., 2002; Stölting and Wilson, 2007). As a consequence both female and male nutrition play a role on the embryonic and fry quality and development, but unfortunately, this topic remains one of the less studied subjects in seahorses (Olivotto et al., 2011b; Cohen et al., 2016).
The European long-snouted seahorse Hippocampus guttulatus Cuvier has been identified as a potential candidate for the marine ornamental trade (Cohen et al., 2016). Rearing techniques are currently available for this species but the breeding success (especially mating events) are still not optimized (Olivotto et al., 2011b; Planas et al., 2017a). Male competition is rather frequent resulting in mating interruption, and juvenile number and size decrease over the time in captive conditions (Faleiro et al., 2008). A bottleneck in breeding performance relies on both type and nutritional quality of diets. In captivity, both frozen and alive feed are usually used to nourish the broodstock, including mysid shrimps, amphipods and Artemia brine shrimps (Woods and Valentino, 2003; Olivotto et al., 2008; Palma et al., 2008; Murugan et al., 2009; Planas et al., 2017a). As concerns mysid shrimps, they are generally considered an optimal diet for seahorses (Woods and Valentino, 2003; Otero-Ferrer et al., 2012). Nutritionally enriched Artemia has received great attention and is presently widely employed (Woods, 2001; Wong and Benzie, 2003; Woods and Valentino, 2003; Planas et al., 2008) especially because it is easily cultured up to the adult stage, provides satisfactory n-3 HUFA and docosahexaenoic acid (DHA) contents after a long-time enrichment (Planas et al., 2017b) and supports high growth rates in H. guttulatus (Planas et al., 2008). However, it has been reported that the breeding performance in H. hippocampus decreased when fed on Artemia compared to diets including mysidaceans (Otero-Ferrer et al., 2012).
In the present study, the following questions have been addressed: (1) how experimental diets are assimilated (stable isotope analysis) by breeders and reflected in general features and biochemical composition of H. guttulatus newborn? (2) do mixed diets perform better than monospecific ones? and (3) Does H. guttulatus follow a mixed income-capital breeding strategy? The main objectives of the study were (1) to propose a breeder’s diet enhancing reproduction performance in the selected species, and (2) to ascertain when and how long the diet should be administered to provide its beneficial effect. To our knowledge, this is the first multiapproach study on the effect of breeder’s diet on reproductive performance and newborn quality in H. guttulatus.
Materials and Methods
Microalgae and Live Prey
Microalgae (Phaeodactylum tricornutum, Rhodomonas lens, and Isochrysis galbana) were cultured at 22 ± 1°C in 80 L plastic bags containing sterilized seawater supplemented with F2P (100 g L–1) media (VarAqua). Additionally, silicates were added to P. tricornutum cultures, and 200 μL F2P media to the R. lens culture flasks.
Artemia cysts (AF, Inve, Spain) were hatched at 28°C for 20 h in 20 L units, and the freshly hatched nauplii (614 ± 140 μm in length) gently rinsed with tapwater, collected on a 125 μm mesh, rinsed, and offered to seahorse juveniles.
Adult non-enriched Artemia was produced to feed adult seahorses during the breeding period. For that, the nauplii were hatched from MC450 cysts (Ocean Nutrition, Spain) and grown in 100 L units, at 26–28°C with gentle aeration and constant light. Adult Artemia was long-time enriched (3–6 days) to feed breeders prior to the start of the experiment but not during the experimental period. The enrichment was carried out in Artemia from day 16 onward on a mixture consisting on live microalgae P. tricornutum and I. galbana (107 cells mL–1), red pepper (0.015 g L–1), and dried Spirulina (0.03 g L–1) (Planas et al., 2017a).
Copepods (Acartia tonsa; 614 ± 140 μm in length) were cultivated for the early feeding of seahorse juveniles in 700 L tanks at 26–27°C and 38 salinity, at an initial density of 1 copepod mL–1. Copepods were fed every two days on the microalgae R. lens (103 cells mL–1). Siphoning of the culture tanks and water renewals (10% of the total volume) were carried out three times per week.
Seahorse Breeding and Rearing
Animal capture, handling and sampling were conducted in compliance with all bioethics standards on animal experimentation of the Spanish Government (Real Decreto 1201/2005, 10th October 2005) and the Regional Government Xunta de Galicia (REGA ES360570202001/15/FUN/BIOL.AN/MPO01).
The adults of H. guttulatus seahorses used in the present study were reared in captivity or collected in Galicia (NW Spain) from August 2016 to January 2017 with permission of the Regional Government Xunta de Galicia. Prior to the start of the experiment, seahorses were maintained at the facilities of Instituto de Investigaciones Marinas (CSIC) in Vigo (Spain) in 320 L aquaria in a semi-closed system with a 10–14% daily renewal of seawater (Planas et al., 2008) and submitted to temperature and photoperiod natural-like regimes (Planas et al., 2010, 2013) fluctuating from 15°C (Winter) to 19°C (Summer) and from 10L:14D (Winter) to 16L:8D (Summer), respectively. Water quality was checked periodically for NO2, NO3, and NH4/NH3 content (0 mg L–1). Salinity and pH levels were 37 ± 2 and 8.0 ± 0.2, respectively, for both species. Wastes and uneaten food were removed daily by siphoning the bottom of aquaria.
Males and females were maintained separate until the start of the experiment in mid-March 2017, when all seahorses were randomly distributed (1:1 sex ratio; six pairs per aquaria) in three 320 L aquaria and maintained for a whole breeding season (until early October 2017). All aquaria were maintained under the same conditions except for the diet administered. Temperature and photoperiod regimes were as described above. The average size (curved standard length) and weight of seahorse breeders (17.4 ± 2.6 cm and 11.50 ± 5.0 g, respectively) did not differ across experimental groups (ANOVA, p = 0.702, 17.4 ± 2.6 cm; ANOVA, p = 0.470, 11.5 ± 5.0 g). Seahorse broodstocks were continuously monitored to check for newborn release from male’s pouch.
Prior to the breeding period, all breeders were fed for 2–5 months on a diet based on enriched adult Artemia and captured/frozen mysidaceans (Siriella armata and Leptomysis sp.). The length of the prey was 5–7 mm for Artemia and about 1–1.5 cm for mysidaceans. During the breeding period, each husbandry aquaria received two daily doses (ad libitum) of different diets (Figure 1):
Figure 1. General scheme of the experimental framework (breeders fed on diets A, AM, and M) and sampling (dorsal fin and newborn). Seahorse illustration by Sergio Picazo Ferro (https://sergiopicazo.com/).
– Diet A: Cultivated adult non-enriched Artemia.
– Diet M: Captured (Siriella armata and Leptomysis sp.) and frozen (Neomysis sp.; Ocean Nutrition, United States) mysidaceans (1:1).
– Diet AM: Mixture (1:1) of diets A and M.
Small portions of the breeders dorsal fin (partial fin-clipping; about 20 m2 surface) were sampled at the onset and at the end of the breeding period, and properly stored for further stable isotopic analysis (SIA) (Valladares and Planas, 2012).
Newly released juveniles were carefully collected by siphoning, counted and transferred (2–5 juveniles L–1) to 30 L pseudo-kreisel aquaria connected to a semi-opened recirculation system (Blanco et al., 2014). Each batch was cultivated in duplicate (whenever possible) until day 7 after male’s pouch release (DAR) at 19°C under a constant 14L:10D photoperiod. The aquaria were filled with seawater filtered by a series of filter-cartridges (20, 10, 5, and 1 μm) and UV treated (76 w; 16 L min–1) (JR1/50). The rearing system included a degasifying column and two 50 L chambers including mechanical (up to 20 μm) and biological filters (perforated plastic bio-balls) and aerators. From the biofilter unit, the seawater was pumped to 36 w UV units (AquaMedic®, Germany) and then to a 50 L reservoir aquarium, being finally routed by gravity toward the rearing aquaria (Planas et al., 2012).
Newborn seahorses were fed according to the following optimized schedule (Blanco and Planas, 2015):
– 0–5 DAR: Two daily doses of Acartia tonsa (0.67 copepods mL–1 dose–1).
– 6–7 DAR: One daily dose of Acartia tonsa (0.67 copepods mL–1 dose–1) and one daily dose of Artemia nauplii (1 Artemia mL–1 dose–1).
Dead seahorses were removed (8:00 am and 15:00 pm) and counted daily. Final survivals of juveniles were recorded at 7 DAR (M50 at 18–21°C occurs at 6.3–6.7 DAR) (Planas et al., 2012).
For each experimental group, samples of newborn juveniles (see further sections for details) were randomly collected before first feeding, euthanized with Tricaine MS-222 (0.1 mg L–1, Sigma Aldrich), rinsed with distilled water, and conserved at 80°C for further FAs and isotopic analyses.
For length and weight measurements, sampled seahorses (n = 15) were transferred to Petri dishes, photographed and weighed (pooled individuals) on a Sartorius microbalance MC210P (±0.01 mg). Curved standard lengths (SL = head + trunk + curved tail) were measured from digital photographs using an image processing software (NIS Elements, Nikon).
Total lipids from live prey and juveniles (pools of 10–20 mg dry weight per sample) (i.e., 25–40 newborn) were extracted according to Bligh and Dyer (1959). Aliquots of total lipid extracts with known lipid content were centrifuged, resuspended with 0.5 M ammonium formate solution, freeze-dried and stored at −80°C until further FAs analyses. Total lipid content was quantified gravimetrically (Herbes and Allen, 1983). FA composition of lipids was analyzed by gas-chromatography (GC) according to Christie (1982). Lipids were transmethylated (Lepage and Roy, 1986) and Fas analyzed by GC (Perkin Elmer, Clarus 500 gas chromatograph) as described in Planas et al. (2010). Samples were analyzed in duplicate.
Carbon (δ13C) and nitrogen (δ15N) isotope signatures and elemental composition (total C and N) were analyzed in both prey and seahorse samples (fin of breeders and bulk juveniles). Frozen juveniles (n > 5; pooled individuals) were dried for 48 h (60°C), weighted and manually homogenized. The analyses were performed on sub-samples of 0.2–1.0 mg dry weight biomass at Servizos de Apoio á Investigación (SAI) of the University of A Coruña (Spain) as reported by Valladares et al. (2016). As defatting of samples prior to isotopic analysis is recommended when lipid content exceeds 5% weight (C:N > 3.56) (Post, 2002) samples were not defatted but arithmetical corrected factors (Supplementary Table 1) were applied considering defatting (copepods, Artemia nauplii and metanauplii, and seahorse fins or juveniles) or defatting + acidification (10% HCl) (Mysidacea and adult Artemia). The samples were measured by continuous flow isotope ratio mass spectrometry using a FlashEA1112 elemental analyzer (Thermo Finnigan, Italy) coupled to a Delta Plus mass spectrometer (FinniganMat, Bremen, Germany) through a Conflo II interface. Stable isotope abundances are expressed as permil (‰) in conventional delta relative to VPDB (Vienna Pee Dee Belemnite) and Atmospheric Air. The precision (standard deviation) for the analysis of δ13C and δ15N of the laboratory standard (acetanilide) was ± 0.15‰ (1-sigma, n = 10). Standards were run every 10 biological samples.
Transmission electron microscopy (TEM) was applied to abdominal portions of seahorses (three newborns per dietary treatment). The samples were processed for Durcupan ACM (Fluka) resin embedding as follows: pieces fixed in 2.5% glutaraldehyde in 0.2 M phosphate-buffered saline (pH 7.4), for 2 h at 4°C, were washed repeatedly in the same buffer and postfixed in 1% osmium tetroxide in 0.2 M phosphate-buffered saline for 1 h at 4°C. Then, the tissues were dehydrated with increasing alcohol concentrations. The dehydrated pieces were embedded in Durcupan ACM (Squadrito et al., 2017). Finally, the sections were obtained with a Reichert Jung Ultracut E. Semi-thin sections (1 μm) were stained with toluidine blue and examined with a light microscope. Ultrathin sections (740 Å) were obtained from selected areas of the semi-thin sections, stained with uranyl acetate and lead citrate, and examined and photographed with a transmission electron microscope (JEOL JEM 1400 Flash) (Viña et al., 2014).
All statistical analyses were conducted in R v.3.6.1 (R Core Team, 2014). Variability across treatments was examined by ANOVA or non-parametric Kruskal–Wallis test for independent samples, depending on significant deviations from normality and homogeneity of variance, as examined by the Shapiro–Wilks and Levene tests, respectively (ggpubr v0.2.5 R package) (Kassambara, 2020a). When significant, A Kruskal–Wallis test followed by Wilcoxon range comparisons were applied (Pgirmess v1.6.9 package of R) (Giraudoux et al., 2018). Significance levels were set at p < 0.05.
The Jonckheere–Terpstra trend test (J–T test) was applied to individual FA data for an ordered alternative hypothesis within an independent sample (between-seahorse groups) design (clinfun v1.0.15 package in R) (Venkatraman, 2018). The test has more statistical power than the Kruskal–Wallis test, incorporating information about whether the order of the groups is meaningful (Siegel and Castellan, 1988). Hence, we applied the test to check whether the more proportions of mysidaceans in the diet (order: A – AM – M), the more the independent variable tested will go down or up.
Correlations between variables in the whole dataset and the resulting network plots of variables were explored with ggraph v2.0.2 (Pedersen, 2020) and corrr v0.4.2 (Kuhn et al., 2020) R packages. Hierarchical clustering of both variables and newborn batches was performed with Ward’s method (provided the highest agglomerative coefficients) using factoextra v1.0.7 (Kassambara, 2020b) and ComplexHeatmap v3.11 (Gu et al., 2016) packages in R. Additionally, principal component analyses (PCA) were performed to summarize and visualize the information of the datasets. For that, we used factoMineR v2.3 (Husson et al., 2020), factoextra v1.0.7 (Kassambara, 2020b), and corrplot v0.8.4 (Wei et al., 2017) packages in R. The data values were standardized (mean = 0; sd = 1) for clustering and PCA.
The protein content in breeder’s diets increased from diets A to M, whereas lipid content decreased due to the lower lipid content in mysidaceans compared to adult Artemia (% C and C/N in Table 1). Prey (and diets) offered to seahorse breeders also differed notably in FA profiles (Table 2 and Supplementary Figure 1). Compared to mysidaceans, adult non-enriched Artemia was richer in saturated FA (35.3%) (16:0, 17:0, and 18:0), n-6 FA (18:2n-6c and 18:3n-6c) (21.5%), and n-7 FA (18:1n-7) (18.8%). Captured and frozen mysidaceans were remarkably richer than Artemia in n-3 HUFA (42.6–45.2%, respectively) (EPA and especially DHA) and showed much higher DHA/EPA, DHA/ARA, and n-3/n-6 ratios.
Table 1. Average isotopic profiles (δ13C‰ and δ15N‰), total C and N, and C/N ratios in prey (A – Artemia, Mc – Captured mysidaceans, and Mf – Frozen mysidaceans) (n = 3) and in diets for breeders (diets A, AM, and M) of Hippocampus guttulatus.
Table 2. Total lipids (% dry weight) and fatty acids (relative percentage) in prey (n = 2), diets and newborn seahorses.
The incorporation of the experimental diets in breeder tissues was analyzed considering dietary isotopic profiles (Table 1) and the resulting isotopic changes in dorsal fins (initial vs final) (Table 2). Differences in isotopic signatures of diets corresponded approximately to one trophic level, with values ranging from −22.12 to −20.53‰ for δ13C (Δ = 1.59) and from 7.67 to 10.83‰ in δ15N (Δ = 3.16).
Initial isotope values for δ13C and δ15N in dorsal fin at the onset of the experimental breeding period did not differ significantly (ANOVA, P = 0.192 and 0.143, respectively) (Table 2). At the end of the experiment, changes in signatures for both isotopes in fins were inversely correlated (linear regression, R2 = 0.994) across treatments (Supplementary Figure 2) and strongly agreed with isotope values of the corresponding diets, with significantly lower δ15N signatures in the group fed uniquely on Artemia (K–W test, P = 0.044). Isotopic changes (Δδ) in fins during the breeding season were 2.97, −0.20, and −1.89 for δ13C, and −1.42, 0.88, and 2.58 for δ15N, respectively (Table 3). Isotopic changes between final isotopic values in dorsal fins and diets A, AM, and M were 8.13, 4.14, and 2.50 for δ13C, and 2.62, 4.25, and 3.91 for δ15N, respectively.
Table 3. Isotope profiles (δ13C‰ and δ15N‰) in dorsal fins (start and end of the experimental breeding period) of Hippocampus guttulatus breeders fed on diets A (Artemia), AM (Artemia + Mysidacea) and M (Mysidacea), and in newborn juveniles.
The total amount of newborn batches produced in treatments A, AM, and M were 4, 7, and 3, respectively, and the total production of newborns were 534, 1,416, and 248, respectively (Table 4). Average batch size in treatment AM (202 ± 82 newborns; max = 343) was higher than in treatments A (134 ± 83) and M (84 ± 72) but differences across treatments were not significant (ANOVA, P = 0.331) due to the large standard deviations of means. Newborn production started earlier and extended over a longer period for group AM breeders (from day 72 to 210 of the breeding period). The opposite occurred in group M (from days 152 to 179). In group A, newborn batches were released from day 121 to day 202.
Table 4. Batches produced, batch size, standard length (mm), dry weight (mg), and survival percentage (7 DAR) of newborn juveniles released by Hippocampus guttulatus breeders fed on diets A (Artemia), AM (Artemia + Mysidacea), and M (Mysidacea).
Newborn differed in length (K–W test, p = 0.018) but not in dry weight (0.60–0.74 mg) (K–W test, P = 0.440) (Table 4 and Supplementary Figure 3). Newborns from treatment A (12.9 ± 0.7 mm) were significantly smaller in length than those from treatments AM (14.7 ± 0.7 mm) and M (14.1 ± 0.2 mm).
Survivals at 7 DAR were high but not significantly different across treatments (K–W test, P = 0.066) (Table 4), ranging from 94.0 ± 5.3% in treatment A to 61.4 ± 19.2% in treatment M. First mortalities in treatment M occurred at 2 DAR, with a peak at 4–5 DAR (Figure 2). Similarly to newborn reared from a batch released by a wild male, first mortalities in treatments A and AM were delayed until 3–4 DAR.
Figure 2. Initial survival of newborn seahorses released by a wild male and by captive males fed on diets A, AM, and M.
Pooled samples of newborn differed in δ15N across treatments (K–W test, P = 0.044) but not in δ13C, ranging from 11.21 to 14.96‰ and from −17.31 to −13.86‰, respectively (Table 2 and Supplementary Figure 3). Mean δ15N value in group A was significantly lower than those from the other groups. However, mean isotopic values did not reflect the pattern of change across treatments during the breeding period. The low number of batches released in treatments A and M did not allow a precise analysis of isotopic changes but a trend could be inferred from treatment AM, in which a progressive change in δ15N toward those in the diet was noticed throughout the breeding season, with signals decreasing from 16.0 to 12.3‰.
Total lipids in newborn were similar in all experimental groups (12.17–12.99% DW) (K–W test, P = 0.580), but significant differences were noticed for some FA in total lipids (Table 2). Significant differences across treatments (K–W test, P < 0.05) were obtained for 15:0, 17:0, 18:1n-7, 20:1n-9, 20:2n-6, 20:4n-6, and 22:6n-3, total polyunsaturated FA and DHA/EPA ratio. J–T test revealed that increasing proportions of mysidaceans in the diet of breeders promoted significant increasing trends in newborn for 14:0, 15:0, 17:0, 18:1n-9t, 20:1n-9, 20:2n-6, 20:4n-6, 22:6n-3, and n-3 HUFA and DHA/EPA ratio, and declines in 16:1n-7, 18:1n-7, and 18:2n-6c (Table 2). The FA profile of newborn released by a wild male better resembled those from groups AM or M.
From the correlation network developed for treatment AM (Figure 3) it can be observed that time (days elapsed from the onset of the breeding season) negatively affected survival, the relative content in saturated FA (particularly 14:0) and δ15N values, but increased the content in total n-6 FA (mainly 18:2n-6c) and 18:3n-3. Length and dry weight of newborn were positively correlated but they did not correlate with time. Size and weight were negatively correlated with δ13C. In addition, weight was negatively and positively correlated with total monoenes and 18:0, respectively. The full results on correlation analyses are provided in Supplementary Figures 4, 5. Sample similarities were remarkably high (R2 > 0.7). However, early released newborn batches from group A were more similar to intermediate batches from group AM, and late released batches were similar, independently of the group considered. The only batch available by a wild male was similar to early batches released by males fed on diet AM.
Figure 3. Correlation network plots for newborn seahorses (Diet AM) showing the more highly correlated variables (Pearson’s r > 0.6), which appear closer together and joined by stronger paths. Paths are also colored by their sign (blue for positive and red for negative). The proximity of the points were determined using multidimensional clustering. Original variable names were modified by R software. Variables: (A) time (days elapsed since the onset of the breeding period), SL, W, Survival at 7 DAR (Surv), total C (C.), total N (N.), C/N, δ13C (C13), δ15N (N15), total lipids (Lip.DW), saturated FA (S), monounsaturated FA (M), polyunsaturated FA (P), n-3 HUFA, n-3 FA, n-6 FA, n-3/n-6 (n3.n6), DHA/EPA (DHA.EPA) and DHA/AA (DHA.AA); (B) time, SL, W, Survival at 7 DAR (Surv) and individual FA.
The results obtained with PCA differed depending on the dataset used (Figure 4 and Supplementary Figure 6). PCA performed with the datasets including general newborn characteristics (Figure 4A) and FA (Figure 4B) showed the proximity of newborn from the wild male to males from groups AM and, to a lesser extent, M. The first two components of PCA accounted for 58.1 and 50.4% of total variation in the dataset. The first component PC1 clearly separated samples from group A (left side) from the others (right side). The most significantly associated variables to PC1 (increasing order of significance) were δ15N, polyunsaturated FA, DHA/AA ratio, n-3 HUFA, DHA/EPA ratio and SL for positive associations, and monounsaturated FA, δ13C and time for negative associations (Supplementary Figures 4A,B). Main FA relevant to PC1 were 20:2n-6, 22:6n-3, 15:0, 20:1n-9, 17:0, 18:1n-9t, and SL (right side), and 16:1n-7, 20:5n-3 and survival at 7 DAR (left side). For PC2, the most relevant proximate variables were n-3/n-6, lipid content, C/N ratio (positive association), and n-6 FA and SL (negative association). For FA, main positive associations with PC2 were time, 20:0, 20:1n-9, 18:1n-9t, 20:3n-6, and 18:2n-6c (positive), and 24:0 and survival (negative).
Figure 4. Factor score plots for the principal component analysis (PCA) on newborn seahorses released by wild (W) or cultivated (A,AM, and M) adults. Only the variables with the highest contributions (cos2 > 0.6) on factors 1 and 2 maps are shown. Analyses performed on datasets (A) – all diets and general characteristics of newborn, (B) – all diets and all individual FA, (C) – all diets and variables selected according to K-M and J–T test, and (D) – diet AM (Same as for C). Original variable names were automatically modified by R software. Variables: days elapsed since the onset of the breeding period (Time), SL, W, Survival at 7 DAR (Surv), total C (C.), total N (N.), C/N (C.N), δ13C (C13), δ15N (N15), total lipids (Lip.DW), saturated FA (S), monounsaturated FA (M), polyunsaturated FA (P), n-3 HUFA, n-3 FA, n-6 FA, n-3/n-6 (n3.n6), DHA/EPA, DHA/AA, and individual FA. The ellipse corresponding to 95% confidence for group AM is shown in yellow shadow. Other ellipses were not computed due to insufficient data.
Variability explained by factors PC1 and PC2 increased up to 69.8% in the PCA performed on all diets and including exclusively significant variables in K-M and J–T tests (Figure 4C). In increasing order of importance, the most significantly associated variables (Supplementary Figure 4C) to PC1were DHA, 20:2n-6, total polyunsaturated FA, total N, δ15N, DHA/EPA ratio, 17:0, AA and 20:1n-9 for positive associations, and 18:1n-7, 16:1n-7 and survival for negative associations. Survival, 14:0 and saturated FA were positively associated to PC2, with negative association for time and monounsaturated 18:1n-9t and 20:1n-9.
Finally, the PCA performed for selected variables in AM dataset (other groups were not analyzed due to insufficient data) explained 69.2% of total variability with PC1 and PC2 (Figure 4D). PC1 was positively associated to DHA, 20:2n-6, total polyunsaturated FA, 15:0, δ15N, DHA/EPA ratio, 17:0, 20:4n-6 and 20:1n-9, and negatively associated to 16:1n-7, 18:1n-8 and survival (Supplementary Figure 4D). For PC2, the mostly positive contributions were due to survival, 14:0, and total saturated FA, with time, 20:1n-9 and 18:1n-9 as negatively associated factors.
Transmission Electron Microscopy
Ultrathin sections of abdominal portions of newborn seahorses were observed by TEM in order to analyze the ultrastructural features of muscular tissue. Generally, TEM micrographs of skeletal muscle sections revealed a regular morphology with a normal arrangement of sarcomeres, resulting in a regularly striated appearance with a normal alternation of dark anisotrope band A (myosin and actin) and light isotrope band I (actin filament) (Figures 5A,B,D). In the longitudinal sections of TEM micrographs of skeletal muscle sections, the myofibrils appeared tightly packed in the scarce amount of cytoplasm. The nuclei were located in the peripheral layer of the fiber and mitochondria beneath the sarcolemma and among the myofibrils. Besides, the regular repeating pattern of myofilaments in the myofibrils were detected displaying dense Z lines, forming a regular zigzag (Figures 5A,B,D) at the center of the light I bands. The A bands showed a central light H zone divided by a dark M line (Figure 5A). In addition, the organization of the sarcoplasmic reticulum and the sarcoplasmic triad, consisting of two terminal cisterns and the transverse tubule T (formed by the invagination of the sarcolemma) were also observed (Figures 5A,C). Scattered satellite cells surrounded by volumetrically modest cytoplasm were often seen in muscle tissue (Figures 5E,6B).
Figure 5. Transmission electron micrographs of newborn seahorse muscle (longitudinal sections) from groups M, AM and A. (A) Group AM: normal sarcomere with light isotrope band (I); Z line; dark anisotrope band (A), central H-zones and in the middle M lines (arrow); triad located at the Z line level. (B) Group A: myofibrils with a slight disarrangement. (C) Group AM: higher magnification. Note tubule T between two terminal cisternae close to each other just forming the longitudinal connection. (D) Group M: a regular zigzag configuration of the Z line of the sarcomere. (E) Group A: interruption of Z line (circled). (F) Higher magnification interruption of above circled Z line. SR, sarcoplasmatic reticule; G, glycogen; Tc, terminal cisternae; T, Tubule T; S, satellite cells.
In the cross sections, the disposition of the two sets of myofilaments in different bands of the sarcomere were analyzed. The myofibrils were well delineated by the close apposition of the sarcoplasmic elements (Figures 6A,B). The nuclei were located in the peripheral layer of the fiber and mitochondria beneath the sarcolemma and among the myofibrils.
Figure 6. Transmission electron micrographs of newborn seahorse muscle (cross sections) from groups AM, M, and A. (A) Group AM: normal structure. (B) Group M: normal structure. (C) Group A: disarrangement of myofibrils (asterisks); myofibrillary interstitial spaces (white arrow heads); discontinuity of the sarcolemma (black arrows); myelin-like bodies (black arrow heads). (D) Group AM: nerve ending. (E) Group M: nerve ending. (F) Group A: nerve ending. S, satellite cells; SR, sarcoplasmic reticule; Sarcolemma (white arrows); Mf, myofibrils; M, mitochondria; G, glycogen; V, vesicles.
The general morphology observed in newborn released by males from diet A appeared similar to the other groups (M and AM). However, some ultrastructural muscle changes were noticed. In particular, group A displayed regressive changes in some muscle fibers characterized by a disarrangement of myofibrils (Figure 5B), with evident myofibrillary interstitial spaces. Besides, a discontinuity of the sarcolemma and myelin-like bodies were present (Figure 6C). In addition, a disorganization of the sarcoplasmic elements and of the invagination of the sarcolemma (i.e., tubules T) was also noticed. An interruption of the Z line was also occasionally evidenced (Figures 5E,F) and neuromuscular junctions were observed between motor nerve endings and muscle fibers. Particularly, the terminal axons of the motor neuron were spherical or oval in shape and fit into smooth depressions in each muscle fiber (Figures 6D–F). At the presynaptic zone, several mitochondria and small clear synaptic vesicles were noticed. Other morphological alteration were not detected at the level of the above structures.
In this study, the isotopic results achieved in dorsal fins of the seahorse H. guttulatus confirmed the effective assimilation of experimental diets by breeders, with effects on average isotopic signatures and FA profiles in newborn seahorses across groups. The overall breeding performance (extension of the effective breeding period, batches produced, and newborn released) resulted enhanced in parents fed on the mixed diet AM (non-enriched Artemia and mysidaceans), suggesting a better nutritional condition than in groups A and M. Unexpectedly, the use of mysidaceans as only prey (diet M) resulted in the worse breeding success among the three experimental diets.
Brood size in some seahorse species (e.g., H. kuda and H. erectus) has been related to the length of brood pouches of males and trunk of females (Dzyuba et al., 2006; Lin et al., 2012). The newborn size might be correlated to brood size, gestation time and nutrient supply in the same males and females’ investment (Lin et al., 2008b, 2012). Besides this, larger H. erectus males (larger brood pouches) do not produce larger broods, and male size can be positively correlated with juvenile size (Lin et al., 2012). In H. guttulatus, it has been pointed out that the size of newborn seahorses can be, in part, paternally determined (limited carrying capacity of the male), as male size is negatively correlated with embryo density and positively correlated with newborn size (Faleiro et al., 2016). As breeders average size did not differ across our experimental groups, we assume that differences between treatments in our study were due to dietary dissimilarities. Differences in the capability of seahorses to digest and assimilate a given prey is another subject to consider when assessing diet quality. The monitoring of prey in feces of H. guttulatus by 18SrDNA amplification showed that gut passage times of crustaceans (Artemia, Palaemonetes, and mysidaceans Leptomysis sp., and Siriella sp.) are species-dependent (longer passages times in mysidaceans), probably due to differences in their digestibility (lower digestibility in Palaemonetes) (Corse et al., 2015).
Mysidaceans have been considered an excellent prey for the on-growing and breeding of some seahorse species (e.g., H. abdominalis and H. hippocampus). Otero-Ferrer et al. (2012) studied the effect of feeding H. hippocampus breeders exclusively on Artemia enriched with DHA or captured mysids (Leptomysis sp.). The authors reported better results on reproduction (spawning events and brood size) and newborn size in breeders fed on mysids, presumptively due to a higher nutritional quality of mysidaceans (basically in DHA content). Those results strongly agree with those with diet AM, but not with the poor results achieved with diet M. we also recommended (but not tested) the use of a mixed diet including Artemia and mysids.
In the present study, premature newborn release or egg clutch dropping were not observed. Such issues were observed in previous studies performed on H. hippocampus fed on Artemia (Otero-Ferrer et al., 2012) or commercial mysidaceans (i.e., much lower DHA and HUFA levels than in the present study) (Saavedra et al., 2015) and in H. guttulatus fed on suboptimal quality of Artemia (i.e., high lipid and low DHA content) (Planas et al., 2009b, 2010). In the later, the issue was solved by optimizing DHA and EPA content in Artemia (Planas et al., 2017b) and isolating pregnant males a few days prior to newborn release (Olivotto et al., 2011b).
Small initial lengths and reduced batch sizes in newborn from group A were accompanied by high survivals at 7 DAR. This result was unexpected as newborn survival commonly depends on the dietary quality provided to broodstock. Interestingly, it has been reported in H. reidi that reducing the nutritional quality of a male’s diet since one month before conception resulted in abnormally large and heavy newborns with extremely reduced survivals, independently of the dietary quality in females (Otero-Ferrer et al., 2016). These findings partially agree with those reported in the pipefish Syngnathus typhle, in which a higher nutritional condition in males brooding heavier embryos was correlated with a higher embryo survival (Sagebakken et al., 2017).
Even though the growth of adult H. guttulatus fed on Artemia may be satisfactory (Planas et al., 2008, 2009a; Segade et al., 2015), it seems that its high lipid content and suboptimal FA profile might result in an excessive accumulation of lipids in the liver and an abnormal (not necessarily pathological) displacement of hepatocytes core to the periphery of the cells (Segade et al., 2016). The presence of ultrastructural muscle alterations at the level of the contractile component and presence of myelin-like bodies (expression of the lysis of the membranes involved by the regressive phenomena) in newborns released by males fed on diet A can be related to some nutritional deficiencies in unenriched Artemia. Polyunsaturated fatty acids (PUFAs) are essential FAs in regulating normal growth and development in many fish species (Anderson et al., 1990; Olivotto et al., 2011a) and promoting muscle fiber development (Cong-cong et al., 2019). Specifically, DHA is a structural key component of cell membranes involved in their biosynthesis (Olivotto et al., 2011a; Randazzo et al., 2018). The low levels of DHA in diet A might thus be responsible for the alterations found in some muscle fibers. Furthermore, an interruption, observed only occasionally of the Z line, could be the expression of the structural disarrangement of myofibrils demonstrated in muscle fibers with regressive phenomena.
The supposed superior quality of several small species of marine crustaceans (e.g., amphipods and mysidaceans) over Artemia has been reported in seahorse culture (Palma et al., 2008; Woods, 2009; Otero-Ferrer et al., 2012; Buen-Ursue et al., 2015; Vargas-Abúndez et al., 2018). Those sources are among the mostly consumed preys by wild seahorses (Manning et al., 2019). However, the main dietary sources for H. guttulatus adults in our nearby region (NW Iberian Peninsula) are amphipoda (caprellids and gammarids), whereas mysidaceans are insignificant prey in the bulk diet composition (Valladares et al., 2016). We recently demonstrated that diets including mysidaceans (alone or mixed with Artemia) provide the best trophic enrichment factors (laboratory derived) to infer the composition of bulk diet in wild syngnathids by means of isotopic signatures (Planas et al., 2020). Hence, even though mysidaceans might be considered a high quality prey for some seahorse species, it is very likely not suitable for breeding H. guttulatus when administered as the sole prey. Unfortunately, nutritional requirements in seahorses are unknown, and provision of suitable food sources for ex-situ culture often represents a serious constraint (Chang and Southgate, 2001) as not all species seem to have the same nutritional needs, as suggested by inter-specific differences in FA profiles (Lin et al., 2008a). Consequently, mixed diets (present study) and/or dietary supplementation would potentially reduce some nutritional deficiencies of single prey (Palma et al., 2012, 2017).
Seahorses are synchronous batch spawners, as several oocyte populations are present during the reproductive season and clutches of oocytes develop, mature, and are released simultaneously (Fernández-Palacios et al., 2011). During the last part of the secondary growth of oocytes (vitellogenesis), lipids and yolk proteins are incorporated and stored in large amounts in the yolk. That process and further embryogenesis of fertilized eggs in males cyclically occurs along the breeding season, requiring a continuous expenditure of energy and nutrients. Seahorse females and males are time-synchronized for egg batch transfer from females to males and newborn release in males, and in H. guttulatus kept at 19–20°C both processes occur at 25–27 days time-intervals (Planas et al., 2010). On this regard, it is very important to ascertain the origin (capital, income or mixed capital-income reproduction strategy) of the energy and nutrients necessary for reproduction (Williams et al., 2017) and to evaluate whether and when egg/fry composition can be affected by the type of diet provided. Hence, the identification of both diet quality and effective feeding periods (the time period needed to change the biochemical composition of eggs/fry) for optimum broodstock performance is of pivotal importance in enhancing newborn quality (Fernández-Palacios et al., 2011; Griffen, 2018). The present study and the pattern of isotopic changes discussed elsewhere (Planas et al., 2020) suggest that the reproduction type followed by H. guttulatus corresponds to a mixed capital-income strategy. This strategy usually implies a progressive drop in newborn quality in successive batches (McBride et al., 2015; Griffen, 2018). In order to properly assess the breeding strategy followed by H. guttulatus, the effects of experimental diets must be analyzed over the entire breeding period. Given that the pattern of change in the biochemical components analyzed in newborn did not follow a mathematical model for fitting, the effective feeding periods were inferred from the smoothed trends in the most remarkable components along the breeding period (particularly for group AM) (Figure 7). In most cases, trend inflexions occurred within the period comprising 100 and 140 days, depending on the variable and seahorse group considered. We identified three stages along the whole breeding period:
Figure 7. Progressive changes in survival at 7 DAR, length (SL), fatty acids (area percentage), and FA ratios along the breeding period in newborn released by seahorses (Wild and groups A, AM and M). Smoothed trends are shown for groups AM (95% confidence intervals also shown) and A. Green shadowed areas show inflection periods indicating a remarkable change (transition from stage 1 to stage 2). Stages: S1 (Stage 1), IP (Inflection period), and S2 (Stage 2).
– Stage 1 (mixed capital-income sources): progressive changes occur since the dietary shift at the onset of the breeding season. Energy and nutritional resources would mostly originate from resources gained prior to the breeding period (capital breeding), accompanied by an increasing use of sources derived from the experimental diet (income breeding). This stage would extend until days 100–140 since the onset of the breeding period.
– Inflection period (income sources): this is a transitory period indicating the end of the capital-income breeding strategy and the use exclusively of resources from the experimental diet. In this period, differences among treatment were clearly perceptible.
– Stage 2 (income sources and exhaustion of body reserves): in most cases, differences among groups increased or decreased (depending on the variable considered) beyond the inflection period. Maximal divergences between treatments were reached during this period, notably by the end of the stage. The trend followed by certain variables (e.g., DHA and DHA/EPA) in late newborn batch releases suggests a more or less intense exhaustion of body resources, accompanied by decreases in survivals at 7 DAR and newborn size. This statement is also supported by the trends in newborn isotopic signatures (Planas et al., 2020).
Fry quality can be manipulated to a certain extend by modifying the nutritional composition of the diet during the breeding season of fishes (Fernández-Palacios et al., 2011) including seahorses (Otero-Ferrer et al., 2020). However, it is also important to assess the effect of preconception diet on reproduction success and the contribution of males and females, as reported by Otero-Ferrer et al. (2016) in H. hippocampus. According to our above hypothesis on the stages recognized along the breeding season, the diet provided preceding the reproduction season would contribute for a certain period to the initial composition and quality of H. guttulatus newborn. Hence, the initial quality of eggs/newborn could be improved by feeding the breeders on a high quality diet for at least 3 months in advance of the onset of the breeding season, which is in agreement with previous findings in the tropical species H. reidi (Planas and Chamorro, 2019). That finding is of primary importance to female diets as it is noteworthy that embryos are mostly dependent on maternal nutrition (yolk quality) until later stages of development, when the embryos are mostly dependent on paternal contribution (Otero-Ferrer et al., 2020).
Breeders diet must be optimized to ensure good early development (i.e., growth, survival, and health) (Izquierdo et al., 2001; Migaud et al., 2013). Currently, nutrients requirements are rather well known in many marine species with interest in ornamental trade, and optimal dietary formulations are available (Craig et al., 2017). However, there are two main issues in the feeding of seahorses in captivity: the lack of knowledge on their nutritional needs and the rejection of inert formulated diets (Olivotto et al., 2011b; Palma et al., 2012; Planas et al., 2017a). As the broodstock diet is strongly reflected in the composition of eggs/fry, particularly for lipids (Wiegand, 1996), the composition of eggs or fry can be used as a proxy to infer nutritional requirements in fishes to improve breeding success (Sargent et al., 1989, 1999; Rainuzzo, 1993; Fernández-Palacios et al., 1995; Rainuzzo et al., 1997). The content in lipids and the FA profiles of H. guttulatus eggs under different dietary regimes were reported by Planas et al. (2008, 2010) and Faleiro and Narciso (2010). Essentially, FA profiles in eggs and newborn do not differ significantly when the content is referred to relative percentage of total FA, but they markedly differ when the content is given as dry weight percentage (Faleiro and Narciso, 2010). In the present study, we focused exclusively on newborn, including a batch released by a wild pregnant male for comparative purposes.
The similarity between the wild batch (released early in the season) and the early and mid-term batches released from group AM was mainly due to DHA, total polyunsaturated FA, δ15N, survival at 7 DAR, DHA/EPA ratio, AA and, to a lesser extent, to minority 20:1n-9 and 20:2n-6 (<1%). The DHA/EPA ratio reflected the higher and more stable content in DHA thorough the breeding season, compared to treatment A. The essentiality of DHA, EPA, and AA in early developmental stages of marine fish has been previously discussed in depth (Rainuzzo, 1993; Bromage, 1995; Wiegand, 1996; Rainuzzo et al., 1997; Estévez et al., 1999; Tocher, 2010). There is no doubt about the nutritional essentiality of certain HUFAs in marine fish. These FA are important sources of energy in fish during early embryogenesis (Sargent et al., 1999) when DHA and AA are highly conserved in most teleost (Watanabe et al., 1989). However, Faleiro and Narciso (2010) analyzed the lipid and FA dynamics during early development of H. guttulatus (breeders fed on adult Artemia; Artemia diet not provided) and reported that EPA, DHA, and palmitic acid were not preferentially retained in embryos as they constitute the main energy sources, with EPA being catabolized at a high rate. Consequently, an exceptionally low lipid content (5.0% DW) and extremely low EPA levels (5.0%) were reported in newborn. Compared to that study, EPA content and total lipids in newborn from our study were two-fold higher, in agreement with previous results in freshly released eggs (Planas et al., 2008, 2010).
In spite of the extremely low levels of DHA in diet A, its average content in newborns was high (with a sharp decrease in late batches) although lower than in groups AM and M. This finding indicates the essentiality of DHA in seahorses. Overall, FA profile and DHA/EPA ratio in newborns released from the wild male showed a higher similarity with newborns released by males fed on diets including mysidaceans. Hence, the variable effect of breeder’s diet on the allocation of FA resources in newborn from different experimental assays is confirmed in H. guttulatus, deserving further attention, especially with regard to dietary DHA requirements and further effects on survival and size of newborn.
Higher dietary contents in total n-3 HUFA and DHA does not necessarily imply higher quality. Excessive levels of dietary n-3 HUFA can reduce eggs production and survival, and cause yolk sac hypertrophy in fish larvae (Fernández-Palacios et al., 1995). As the TEM did not show any morphological alteration in tissues of newborn from diet M, we suspect that n-3 HUFA and DHA levels in diet M were below detrimental levels, and that the low breeding success in that treatment was due in part to energetic and/or nutritional dietary deficiencies. It is known that egg viability may decrease in lower quality diets (Saavedra et al., 2015) which is in agreement with the reduced newborn production in diets A and M. As eggs dropping did not occur in our study and female-male egg transfers were not checked, viability in eggs transferred to males is not available. Hence, the reasons for that low performance are uncertain. In addition to potential negative effects of low quality diets on mating displays, a drop in reproductive success due to the strategy of cryptic choice should not be discarded. According to this strategy, males could increase rates of offspring abortion in pregnancies to absorb nutrients nutrients for future reproductive opportunities (Paczolt and Jones, 2010; Sagebakken et al., 2010).
In the experimental study presented here, isotopic signatures of experimental diets supplied during the breeding season were reflected in breeder’s dorsal fins (efficient assimilation of diets) and subsequently in newborn, with higher δ15N signals in breeders and newborn from the diets including mysidaceans. Monospecific diets including Artemia or mysidaceans very likely suffered from nutritional deficiencies, reflected in morphological alterations of the muscle tissue or abnormal early mortalities in newborn. Those deficiencies seemed to be overcome with the mixed diet AM (unenriched Artemia and mysidaceans), which provided the highest overall breeding success and reasonable FA profiles, and better resembled newborn from a wild male. However, special attention deserve the progressive changes in FA profiles along the breeding season, resulting in a decrease in newborn performance. Considering the time required for a diet to be reflected in newborn, we recommend that H. guttulatus breeders be fed on a high quality diet at least 3 months in advance of the onset of the breeding period.
In order to further optimize captive breeding and rearing of seahorses, further research should focus on wild populations with particular emphasis on broodstock feeding habits, on how feeding preferences possibly change during the breeding season and on potential shifts depending on the geographical distribution of the fish. Data collected from the wild will help scientist in the selection of the more suitable diet able to sustain the best egg and fry quality.
Data Availability Statement
The datasets analyzed for this study can be found in the Pangaea collection (https://doi.pangaea.de/10.1594/PANGAEA.919218). Data not listed in this collection is available from the authors upon reasonable request.
The animal study was reviewed and approved by Regional Government Xunta de Galicia (REGA ES360570202001/15/FUN/BIOL.AN/MPO01).
MP as leading author has done the conception of the work, breeding and rearing, sampling, lipid and stable isotope analysis, interpretation of data, and most part in drafting. IO and RL assisted in manuscript writing. MG was responsible for fatty acid analysis. RL and MZ were responsible for TEM and this part of the draft. All authors revised the manuscript. All authors contributed to the article and approved the submitted version.
This study was financially supported by the Spanish Government with Project Hippoeco (Ref. CGL2015-68110-R, Ministerio de Ciencia, Innovación y Universidades and Fondos FEDER).
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.
We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI). We are grateful to Alexandro Chamorro, Alex Paltrinieri, David Costas, Isabel Medina, and the staff of the Aquaria Service at IIM-CSIC for assistance in the maintenance of seahorses in captivity and technical support. We are also grateful to the Regional Government Xunta de Galicia for providing the permission for the capture of wild adult seahorses Hippocampus guttulatus. Finally, we are very grateful to Sergio Picazo (https://sergiopicazo.com/) for providing usage rights of seahorse illustration in Figure 1.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.00638/full#supplementary-material
Adams, M. B., Powell, M. D., and Purser, G. J. (2001). Effect of acute and chronic ammonia and nitrite exposure on oxygen consumption and growth of juvenile big bellied seahorse. J. Fish Biol. 58, 848–860. doi: 10.1111/j.1095-8649.2001.tb00535.x
Anderson, G. J., Connor, W. E., and Corliss, J. D. (1990). Docosahexaenoic acid is the preferred dietary n-3 fatty acid for the development of the brain and retina. Pediatr. Res. 27, 89–97. doi: 10.1203/00006450-199001000-00023
Binh, M. N., and Serrano, A. E. (2012). Optimum level of fish liver oil as enrichment for Artemia fed to the tiger tail seahorse Hippocampus comes for reproduction and juvenile survival. AACL Bioflux 5, 249–258.
Blanco, A., Clamorro, A., and Planas, M. (2014). Implications of physical key factors in the early rearing of the long-snouted seahorse Hippocampus guttulatus. Aquaculture 433, 214–222. doi: 10.1016/j.aquaculture.2014.06.019
Blanco, A., and Planas, M. (2015). Mouth growth and prey selection in juveniles of the European long-snouted seahorse, Hippocampus guttulatus. J. World Aquac. Soc. 46, 596–607. doi: 10.1111/jwas.12240
Bromage, N. R. (1995). “Broodstock management and seed quality – general considerations,” in Broodstock Management and Egg and Larval Quality, eds N. R. Bromage and R. J. Roberts (London: Blackwell), 1–24.
Buen-Ursue, S. M. A., Azuma, T., Arai, K., and Coloso, R. M. (2015). Improved reproductive performance of tiger tail seahorse, Hippocampus comes, by mysid shrimp fed singly or in combination with other natural food. Aquacult. Int. 23, 29–43. doi: 10.1007/s10499-014-9795-1
Carcupino, M., Baldacci, A., Mazzini, M., and Franzoi, P. (2002). Functional significance of the male brood pouch in the reproductive strategies of pipefishes and seahorses: a morphological and ultrastructural comparative study on three anatomically different pouches. J. Fish Biol. 61, 1465–1480. doi: 10.1111/j.1095-8649.2002.tb02490.x
Chang, M., and Southgate, P. C. (2001). Effects of varying dietary fatty acid composition on growth and survival of seahorse, Hippocampus sp., juveniles. Aquarium Sci. Conserv. 3, 205–214. doi: 10.1023/A:1011363807074
Cohen, F. P. A., Valenti, W. C., Planas, M., and Calado, R. (2016). Seahorse aquaculture, biology and conservation: knowledge gaps and research opportunities. Rev. Fish. Sci. Aquac. 25, 100–111. doi: 10.1080/23308249.2016.1237469
Cong-cong, W., Wen-bin, L., Yang-yang, H., Xi, W., Xiang-fei, L., Ding-dong, Z., et al. (2019). Dietary DHA affects muscle fiber development by activating AMPK/Sirt1 pathway in blunt snout bream (Megalobrama amblycephala). Aquaculture 518:734835.
Corse, E., Valladares, S., Planas, M., Chamorro, A., and Pintado, J. (2015). Analysis of the diet of the long-snouted seahorse Hippocampus guttulatus by 18SrDNA amplification of prey in faeces. Aquac. Nutr. 21, 528–540. doi: 10.1111/anu.12189
Craig, R. S., Gardner, T. R., and Carnevali, O. (2017). “Growout and broodstock nutrition,” in Marine Ornamental Species Aquaculture, eds R. Calado, I. Olivotto, M. Planas, and G. J. Holt (West Sussex: Wiley Blackwell), 139–158.
Dzyuba, B., Van Look, K. J. W., Cliffe, A., Koldewey, H. J., and Holt, W. V. (2006). Effect of parental age and associated size on fecundity, growth and survival in the yellow seahorse Hippocampus kuda. J. Exp. Biol. 209, 3055–3061. doi: 10.1242/jeb.02336
Estévez, A., McEvoy, L. A., Bell, J. G., and Sargent, J. R. (1999). Growth, survival, lipid composition and pigmentation of turbot (Scophthalmus maximus) larvae fed live-prey enriched in arachidonic and eicosapentaenoic acids. Aquaculture 180, 321–343. doi: 10.1016/S0044-8486(99)00209-4
Faleiro, F., Almeida, A. J., Ré, P., and Narciso, L. (2016). Size does matter: an assessment of reproductive potential in seahorses. Anim. Reprod. Sci. 170, 61–67. doi: 10.1016/j.anireprosci.2016.04.003
Faleiro, F., and Narciso, L. (2010). Lipid dynamics during early development of Hippocampus guttulatus seahorses: searching for clues on fatty acid requirements. Aquaculture 37, 56–64. doi: 10.1016/j.aquaculture.2010.07.005
Faleiro, F., Narciso, L., and Vicente, L. (2008). Seahorse behaviour and aquaculture: How to improve Hippocampus guttulatus husbandry and reproduction? Aquaculture 28, 33–40. doi: 10.1016/j.aquaculture.2008.05.038
Fernández-Palacios, H., Izquierdo, M. S., Robaina, L., Valencia, A., Salhi, M., and Vergara, J. M. (1995). Effect of n - 3 HUFA level in broodstock diets on egg quality of gilthead sea bream (Sparus aurata L.). Aquaculture 132, 325–337. doi: 10.1016/0044-8486(94)00345-O
Fernández-Palacios, H., Norberg, B. G., Izquierdo, M., and Hamre, K. (2011). “Effects of broodstock diets on eggs and larvae,” in Larval Fish Nutrition, ed. J. Holt (London: Wiley-Blackwell), 153–181. doi: 10.1002/9780470959862.ch5
Giraudoux, P., Antonietti, J. P., Beale, C., Pleydell, D., and Treglia, M. (2018). Spatial Analysis and Data Mining for Field Ecologists v. 1.6.9. Available online at: https://cran.r-project.org/web/packages/pgirmess/index.html (accessed January 16, 2018).
Gu, Z., Eils, R., and Schlesner, M. (2016). Complex Heatmaps Reveal Patterns and Correlations in Multidimensional Genomic Data. Bioinformatics. Available online at: http://bioconductor.org/packages/release/bioc/html/ComplexHeatmap.html (accessed May 11, 2020).
Husson, F., Josse, J., Le, S., and Mazet, J. (2020). FactoMineR: Multivariate Exploratory Data Analysis and Data Mining. Available online at: https://CRAN.R-project.org/package=FactoMineR (accessed March 20, 2020).
Kassambara, A. (2020a). ggpubr: ’ggplot2’ Based Publication Ready Plots. Available online at: https://CRAN.R-project.org/package=clinfun (accessed February 23, 2020).
Kassambara, A. (2020b). factoextra: Extract and Visualize the Results of Multivariate Data Analyses. Available online at: https://CRAN.R-project.org/package=factoextra (accessed February 23, 2020).
Kuhn, M., Jackson, S., and Cimentada, J. (2020). Correlations in R. Available online at: https://CRAN.R-project.org/package=corrr (accessed March 27, 2020).
Lin, Q., Gao, Y. L., Sheng, J. Q., Chen, Q. X., Zhang, B., and Lu, J. Y. (2007). The effect of food and the sum of effective temperature on the embryonic development of the seahorse, Hippocampus kuda Bleeker. Aquaculture 262, 481–492. doi: 10.1016/j.aquaculture.2006.11.011
Lin, Q., Li, G., Qin, G., Lin, J., Huang, L., Sun, H., et al. (2012). The dynamics of reproductive rate, offspring survivorship and growth in the lined seahorse, Hippocampus erectus Perry, 1810. Biol. Open 1, 391–396. doi: 10.1242/bio.2012398
Lin, Q., Lin, J., Lu, J., and Li, B. (2008a). Biochemical composition of six seahorse species, Hippocampus sp., from the Chinese, Coast. J. World Aquac. Soc. 29, 225–234. doi: 10.1111/j.1749-7345.2008.00159.x
Lin, Q., Lu, J., Gao, Y., Shen, L., Cai, J., and Luo, J. (2006). The effect of temperature on gonad, embryonic development and survival rate of juvenile seahorses, Hippocampus kuda Bleeker. Aquaculture 254, 701–713. doi: 10.1016/j.aquaculture.2005.11.005
Manning, C. G., Foster, S. J., and Vincent, A. C. J. (2019). A review of the diets and feeding behaviours of a family of biologically diverse marine fishes (Family Syngnathidae). Rev. Fish Biol. Fish. 29, 197–221. doi: 10.1007/s11160-019-09549-z
Martin-Smith, K. M., Samoilys, M. A., Meeuwig, J. J., and Vincent, A. C. J. (2004). Collaborative development of management options for an artisanal fishery for seahorses in central Philippines. Ocean Coast. Manag. 47, 165–193. doi: 10.1016/j.ocecoaman.2004.02.002
McBride, R. S., Somarakis, S., Fitzhugh, G. R., Albert, A., Yaragina, N. A., Wuenschel, M. J., et al. (2015). Energy acquisition and allocation to egg production in relation to fish reproductive strategies. Fish Fish. 16, 23–57. doi: 10.1111/faf.12043
Migaud, H., Bell, G., Cabrita, E., McAndrew, B., Davie, A., Bobe, J., et al. (2013). Gamete quality and broodstock management in temperate fish. Rev. Aquac. 5(Suppl. 1), S194–S223. doi: 10.1111/raq.12025
Murugan, A., Dhanya, S., Sreepada, R. A., Rajagopal, S., and Balasubramanian, T. (2009). Breeding and mass scale rearing of three spotted seahorse, Hippocampus trimaculatus Leach under captive conditions. Aquaculture 290, 87–96. doi: 10.1016/j.aquaculture.2009.01.033
Olivotto, I., Avella, M., Sampaolesi, G., Piccinetti, C. C., Ruiz, P., and Carnevali, O. (2008). Breeding and rearing the longsnout seahorse Hippocampus reidi: rearing and feeding studies. Aquaculture 283, 92–96.
Olivotto, I., Di Stefano, M., Rosetti, S., Cossignani, L., Pugnaloni, A., Giantomassi, F., et al. (2011a). Live prey enrichment, with particular emphasis on HUFAs, as limiting factor in false percula clownfish (Amphiprion ocellaris, Pomacentridae) larval development and metamorphosis: molecular and biochemical implications. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 159, 207–218. doi: 10.1016/j.cbpa.2011.02.004
Olivotto, I., Planas, M., Simões, N., Holt, G. J., Avella, M. A., and Calado, R. (2011b). Advances in breeding and rearing marine ornamentals. J. World Aquac. Soc. 42, 135–166. doi: 10.1111/j.1749-7345.2011.00453.x
Otero-Ferrer, F., Izquierdo, M., Fazeli, A., Band, F., and Holt, W. V. (2016). Embryonic developmental plasticity in the long-snouted seahorse (Hippocampus reidi, Ginsburg 1933) in relation to parental preconception diet. Reprod. Fertil. Dev. 28, 1020–1028. doi: 10.1071/rd14169
Otero-Ferrer, F., Lättekivi, F., Ord, J., Reimann, E., Kõks, S., Izquierdo, M., et al. (2020). Time-critical influences of gestational diet in a seahorse model of male pregnancy. J. Exp. Biol. 2020:jeb210302. doi: 10.1242/jeb.210302
Otero-Ferrer, F., Molina, L., Socorro, J., Fernández-Palacios, H., and Izquierdo, M. (2012). Effect of different live prey on spawning quality of short-snouted seahorse, Hippocampus hippocampus (Linnaeus, 1758). J. World Aquac. Soc. 43, 174–186. doi: 10.1111/j.1749-7345.2012.00550.x
Palma, J., Andrade, J. P., and Bureau, D. P. (2012). Growth, reproductive performances, and brood quality of long snout seahorse, Hippocampus guttulatus, fed enriched shrimp diets. J. World Aquac. Soc. 43, 802–813. doi: 10.1111/j.1749-7345.2012.00611.x
Palma, J., Andrade, J. P., and Bureau, D. P. (2017). The impact of dietary supplementation with astaxanthin on egg quality and growth of long snout seahorse (Hippocampus guttulatus) juveniles. Aquac. Nutr. 23, 304–312. doi: 10.1111/anu.12394
Palma, J., Stockdale, J., Correia, M., and Andrade, J. P. (2008). Growth and survival of adult long snout seahorse (Hippocampus guttulatus) using frozen diets. Aquaculture 278, 55–59. doi: 10.1016/j.aquaculture.2008.03.019
Pavlov, D., Kjørsvik, E., Refsti, T., and Anderson, O. (2004). “Brood stock and egg production,” in Culture of Cold-Water Marine Fish, eds E. Moksness, E. Kjorsvik, and Y. Olsen (Oxford: Blackwell), 129–203.
Pedersen, T. L. (2020). ggraph: An Implementation of Grammar of Graphics for Graphs and Networks. Available online at: https://CRAN.R-project.org/package=ggraph (accessed March 26, 2020).
Planas, M., Blanco, A., Chamorro, A., Valladares, S., and Pintado, J. (2012). Temperature-induced changes of growth and survival in the early development of the seahorse Hippocampus guttulatus. J. Exp. Mar. Biol. Ecol. 438, 154–162. doi: 10.1016/j.jembe.2012.10.003
Planas, M., Burhans, R., and Simões, N. (2017a). “Seahorses and pipefish,” in Marine Ornamental Species Aquaculture, eds R. Calado, I. Olivotto, M. Planas, and G. J. Holt (West Sussex: Wiley Blackwell), 299–326.
Planas, M., Silva, C., Quintas, P., Chamorro, A., and Piñero, S. (2017b). Ongrowing and enhancement of n-3 HUFA profile in adult Artemia: short- vs long-time enrichment. J. Appl. Phycol. 29, 1409–1420. doi: 10.1007/s10811-016-1016-z
Planas, M., and Chamorro, A. (2019). “Inheritance of stable isotopes in seahorses (Hippocampus reidi): effect of parental diet preceding breeding season,” in Proceedings of the 43rd Annual Larval Fish Conference, Palma. doi: 10.13140/RG.2.2.16922.08648
Planas, M., Chamorro, A., Paltrinieri, A., Campos, S., Jiménez, A., Nedelec, K., et al. (2020). Effect of diet on breeders and inheritance in syngnathids: application of isotopic experimentally derived data to field studies. Mar. Ecol. Prog. Ser. doi: 10.3354/meps13315
Planas, M., Chamorro, A., Quintas, P., and Vilar, A. (2008). Establishment and maintenance of threatened long-snouted seahorse, Hippocampus guttulatus, broodstock in captivity. Aquaculture 283, 19–28. doi: 10.1016/j.aquaculture.2008.06.023
Planas, M., Quintas, P., and Chamorro, A. (2013). Maturation of Hippocampus guttulatus and Hippocampus hippocampus females by manipulation of temperature and photoperiod regimes. Aquaculture 388-391, 147–152. doi: 10.1016/j.aquaculture.2013.01.030
Planas, M., Quintas, P., Chamorro, A., and Silva, C. (2010). Female maturation, egg characteristics and fatty acids profile in the seahorse Hippocampus guttulatus. Anim. Reprod. Sci. 122, 66–73. doi: 10.1016/j.anireprosci.2010.07.008
Planas, M., Quintas, P., and Chamorro, A. (2009a). “Growth of adult seahorses Hippocampus guttulatus fed exclusively on enriched adult Artemia (Project Hippocampus),” in Proceedings of the Book of Abstracts World Aquaculture Society World Aquaculture 2009, Veracruz, 669.
Planas, M., Quintas, P., and Silva, C. (2009b). “Snout abnormalities in young seahorses (Hippocampus guttulatus) (Project Hippocampus),” in Proceedings of the Book of Abstracts World Aquaculture Society Meeting, World Aquaculture 2009, Veracruz, 668.
Poortenaar, C. W., Woods, C. M. C., James, P. J., Giambartolomei, F. M., and Lokman, P. M. (2004). Reproductive biology of female big-bellied seahorses. J. Fish Biol. 64, 717–725. doi: 10.1111/j.1095-8649.2004.00341.x
Rainuzzo, J. R. (1993). “Fatty acid and lipid composition of fish eggs and larvae,” in Fish Farming Techonology, eds H. Reinertsen, L. A. Dahle, J. Jørgensen, and K. Tvinnereim (Rotterdam: Balkema), 43–48.
Randazzo, B., Rolla, L., Ofelio, C., Planas, M., Gioacchini, G., Vargas Abúndez, J. A., et al. (2018). The influence of diet on the early development of two seahorse species (H. guttulatus and H. reidi): traditional and innovative approaches. Aquaculture 490, 75–90. doi: 10.1016/j.aquaculture.2018.02.029
Saavedra, M., Batista, H., and Pousão-Ferreira, P. (2015). Dietary fatty acid enrichment during the spawning season increases egg viability and quality in Hippocampus hippocampus. Aquac. Nutr. 22, 343–351. doi: 10.1111/anu.12252
Sagebakken, G., Ahnesjö, I., Mobley, K. B., Braga Gonçalves, I., and Kvarnemo, C. (2010). Brooding fathers, not siblings, take up nutrients from embryos. Proc. R. Soc. B 277, 971–977. doi: 10.1098/rspb.2009.1767
Sargent, J., McEvoy, L., Estévez, A., Bell, G., Bell, M., Henderson, J., et al. (1999). Lipid nutrition of marine fish during early development: current status and future directions. Aquaculture 179, 217–229. doi: 10.1016/S0044-8486(99)00191-X
Sarkar, D., Andrews, F., Wright, K., Klepeis, N., and Murrell, P. (2020). lattice: Trellis Graphics for R. Available online at: https://cran.r-project.org/web/packages/lattice/index.html (accessed September 17, 2019).
Segade, A., Robaina, L., Novelli, B., Otero-Ferrer, F., and Molina Domínguez, L. (2016). Effect of the diet on lipid composition and liver histology of short snout seahorse Hippocampus hippocampus. Aquac. Nutr. 22, 1312–1319. doi: 10.1111/anu.12341
Segade, A., Robaina, L., Otero-Ferrer, F., García Romero, J., and Molina Domínguez, L. (2015). Effects of the diet on seahorse (Hippocampus hippocampus) growth, body colour and biochemical composition. Aquac. Res. 21, 807–813. doi: 10.1111/anu.12202
Squadrito, F., Micali, A., Rinaldi, M., Irrera, N., Marini, H., Puzzolo, D., et al. (2017). Polydeoxyribonucleotide, an adenosine-A2 A receptor agonist, preserves blood testis barrier from cadmium-induced injury. Front. Pharmacol. 10:537. doi: 10.3389/fphar.2016.00537
Venkatraman, E. S. (2018). Clinical Trial Design and Data Analysis Functions. Available online at: https://CRAN.R-project.org/package=clinfun (accessed February 4, 2020).
Viña, E., Parisi, V., Abbate, F., Cabo, R., Guerrera, M. C., Laurà, R., et al. (2014). Acid-sensing ion channel 2 (ASIC2) is selectively localized in the cilia of the non-sensory olfactory epithelium of adult zebrafish. Histochem. Cell Biol. 143, 59–68. doi: 10.1007/s00418-014-1264-4
Watanabe, T., Izquierdo, M. S., Takeuchi, T., Satoh, S., and Kitajima, C. (1989). Comparison between eicosapentaenoic and docosahexaenoic acids in terms of essential fatty acid efficacy in larval red seabream. Nippon Suisan Gakk. 55, 1635–1640. doi: 10.2331/suisan.55.1635
Wei, T., Simko, V., Levy, M., Xie, Y., Jin, Y., and Zemla, J. (2017). corrplot: Visualization of a Correlation Matrix. Available online at: https://CRAN.R-project.org/package=corrplot (accessed March 6, 2020).
Wickham, H., Chang, W., Henry, L., Pedersen, T. L., Takahashi, K., Wilke, C., et al. (2020). ggplot2: Create Elegant Data Visualisations Using the Grammar of Graphics. Available online at: https://CRAN.R-project.org/package=ggplot2 (accessed September 17, 2019).
Williams, C. T., Klaassen, M., Barnes, B. M., Buck, C. L., Arnold, W., Giroud, S., et al. (2017). Seasonal reproductive tactics: annual timing and the capital-to-income breeder continuum. Philos. Trans. R Soc. B 372:20160250. doi: 10.1098/rstb.2016.0250
Wilson, A. B., Vincent, A., Ahnesjö, I., and Meyer, A. (2001). Male pregnancy in seahorses and pipefishes (Family Syngnathidae): rapid diversification of paternal brood pouch morphology inferred from a molecular phylogeny. J. Heredity 92, 159–166. doi: 10.1093/jhered/92.2.159
Wong, J. M., and Benzie, J. A. H. (2003). The effects of temperature, Artemia enrichment, stocking density and light on the growth of juvenile seahorses, Hippocampus whitei (Bleeker 1855), from Australia. Aquaculture 228, 107–121. doi: 10.1016/S0044-8486(03)00320-X
Woods, C. M. C., and Valentino, F. (2003). Frozen mysids as an alternative to live Artemia in culturing seahorses Hippocampus abdominalis. Aquac. Res. 34, 757–763. doi: 10.1046/j.1365-2109.2003.00882.x
Zhang, H., Chen, L., Zhang, B., and Lin, Q. (2019). Molecular identification of GnIH and its potential role in reproductive physiology and male pregnancy of the lined seahorse (Hippocampus erectus). Gen. Comp. Endocr. 279, 196–202. doi: 10.1016/j.ygcen.2019.04.018
Keywords: syngnathidae, seahorse, Hippocampus, diet, breeding, newborn, quality
Citation: Planas M, Olivotto I, González MJ, Laurà R and Zarantoniello M (2020) A Multidisciplinary Experimental Study on the Effects of Breeders Diet on Newborn Seahorses (Hippocampus guttulatus). Front. Mar. Sci. 7:638. doi: 10.3389/fmars.2020.00638
Received: 09 June 2020; Accepted: 13 July 2020;
Published: 31 July 2020.
Edited by:Sami Souissi, Lille University of Science and Technology, France
Reviewed by:Khor Waiho, University of Malaysia Terengganu, Malaysia
Yun-Zhang Sun, Jimei University, China
Copyright © 2020 Planas, Olivotto, González, Laurà and Zarantoniello. 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: Miquel Planas, firstname.lastname@example.org