Determination of biological and physicochemical parameters of Artemia franciscana strains in hypersaline environments for aquaculture in the Colombian Caribbean

Background Artemia (Crustacea, Anostraca), also known as brine shrimp, are typical inhabitants of extreme environments. These hypersaline environments vary considerably in their physicochemical composition, and even their climatic conditions and elevation. Several thalassohaline (marine) environments along the Colombian Caribbean coast were surveyed in order to contribute to the knowledge of brine shrimp biotopes in South America by determining some vital biological and physicochemical parameters for Artemia survival. Additionally, cyst quality tests, biometrical and essential fatty acids analysis were performed to evaluate the economic viability of some of these strains for the aquaculture industry. Results In addition to the three locations (Galerazamba, Manaure, and Pozos Colorados) reported in the literature three decades ago in the Colombian Caribbean, six new locations were registered (Salina Cero, Kangaru, Tayrona, Bahía Hondita, Warrego and Pusheo). All habitats sampled showed that chloride was the prevailing anion, as expected, because of their thalassohaline origin. There were significant differences in cyst diameter grouping strains in the following manner according to this parameter: 1) San Francisco Bay (SFB-Control, USA), 2) Galerazamba and Tayrona, 3) Kangarú, 4) Manaure, and 5) Salina Cero and Pozos Colorados. Chorion thickness values were smaller in Tayrona, followed by Salina Cero, Galerazamba, Manaure, SFB, Kangarú and Pozos Colorados. There were significant differences in naupliar size, grouping strains as follows (smallest to largest): 1) Galerazamba, 2) Manaure, 3) SFB, Kangarú, and Salina Cero, 4) Pozos Colorados, and 5) Tayrona. Overall, cyst quality analysis conducted on samples from Manaure, Galerazamba, and Salina Cero revealed that all sites exhibited a relatively high number of cysts.g-1. Essential fatty acids (EFA) analysis performed on nauplii from cyst samples from Manaure, Galerazamba, Salina Cero and Tayrona revealed that cysts from all sites exhibited high arachidonic acid:20:4(n-6) (ArA) and eicosapentaenoic acid: 20:5(n-3) (EPA) levels comparable to the control sample (SFB). In contrast, most cysts collected (including SFB) at different locations, and during different months, presented low docosahexaenoic acid: 22:6(n-3) (DHA) levels (Manaure was the only exception with high DHA levels). Some variations in EPA and ArA levels were observed in all sites, contrasting with the much lower DHA levels which remained constant for all locations, except for Manaure which exhibited variable DHA levels. DHA/EPA ratio was overall very low for all sites compared to SFB cysts. All strains had a low DHA/ArA, but a high EPA/ArA ratio, including the control. Conclusion The Colombian A. franciscana habitats analyzed were determined to be thalassohaline, and suitable for A. franciscana development. EFA profiles demonstrated that Tayrona, Galerazamba, Manaure and Salina Cero strains are suitable food for marine fish and crustacean culture because of their high EPA/ArA ratio, but might have to be fortified with DHA rich emulsions depending on the nutritional requirements of the species to be cultured, because of their overall low DHA content. The relatively small nauplii are appropriate for marine larvaeculture. In contrast, the strains from Tayrona, Kangarú, Salina Cero, and Pozos Colorados may be of use but limited to Artemia small biomass production quantities, because of the small surface area of their respective locations; Artemia could be exploited at these locations for local aquaculture applications. In general, cyst quality evaluation for Manaure, Salina Cero and Galerazamba cysts revealed that cysts from these three locations could improve their quality by concentrating efforts on cyst processing techniques. Finally, most locations had great A. franciscana production potential and require different degrees of water quality and/or infrastructure management.


Background
Members of the genus Artemia (Crustacea, Anostraca), also known as brine shrimp, are typical inhabitants of extreme environments that have low species diversity and simple trophic structures [1]. These hypersaline environments vary considerably in terms of ionic composition, climatic conditions and altitude. As a general rule, chloride-rich lakes are the most adequate for Artemia development [2]; however, some strains require carbonate-(Mono lake Artemia, USA) or sulfate-rich waters (Tso Kar Lake Artemia, Tibet) for survival [1, [3][4][5]. Conversely, some other ions may be deleterious to Artemia; potassium could be very toxic because of its occurrence with sodium [6,7]. Since water composition is important for Artemia survival, ecosystems where it occurs were classified in three categories based on their anionic composition: chloride-, sulfate-and carbonate-rich [2].
To characterize Artemia strains the aquaculture industry employ an array of evaluation tools [13,14]: i) Cyst and nauplii biometry: cyst biometry assists in the determination of number of cysts.g -1 ; generally, 1 g from strains that produce small cysts contains more cysts.g -1 , thus usually pro-duce more nauplii.g -1 . Similarly, naupliar biometry is also an essential tool for quality evaluation. Shorter naupliar length is important, particularly, to feed fish larvae which contrary to crustaceans, has to engulf prey in a single bite; ii) Cyst hatching characteristics: it can be affected by environmental factors, genotypical conditions and/or improper processing/storage. An acceptable cyst product should contain minimal quantities of impurities (i.e. sand, salt crystals, etc.) [13]. Hatching efficiency (HE) and hatching percentage (H%) vary greatly among commercial batches and account for most of the price difference [13]. However, HE may be a better criterion than H% since HE considers impurities content (i.e. empty cyst shells). Hatching values for a commercial product may be as low as 100,000 nauplii.g -1 , but might yield ideally near 300,000 nauplii.g -1 (with H% > 90). Hatching synchrony must be high (T s = 12-16 h), and the last nauplii should hatch within 8 h after T 90 [13]. When T s is low (T 0 -T 100 > 10 h), first hatched nauplii will have consumed most of their energy reserves by the time that the last nauplii hatched and harvesting has ended [13]; and iii) Essential fatty acid (EFA) profiles: particular attention has been given by marine larvaeculture production facilities to search for EFA rich Artemia strains, i.e. arachidonic acid:20:4(n-6) (ArA), eicosapentaenoic acid: 20:5(n-3) (EPA), and docosahexaenoic acid: 22:6(n-3) (DHA) [15]. Since ArA and DHA are vital in marine fish nutrition [16][17][18] great effort has been devoted to incorporate high levels of ArA, DHA, and high ratios of DHA/EPA/ArA in live food. ArA is generally conserved during periods of starvation in marine fishes [19], and serves as the preferred precursor for eicosanoid biosynthesis [20]. Artemia fed with enriched n-3 and n-6 highly unsaturated fatty acids (HUFA) in turn results in better larval growth and survival of several marine species fed with it [18,[21][22][23][24][25]. The DHA/ EPA ratio is very variable in non-enriched Artemia, with values often lower than 1. Through the addition of DHA rich emulsions the DHA/EPA ratio increases up to 7 [15]. Enrichment success is strain dependent (i.e. particular Chinese strains), and linked to variations in DHA catabolism. Further, EFA nutritional requirements may vary between species and developmental stages [18]. In white bass larvae, the optimal DHA/ArA and EPA/ArA dietary ratios have been established at 2:1 and 1:1, respectively [26]. This contrasts with flat fish larvae (i.e. turbot and Atlantic halibut) which require much higher ratios of over 10:1 [27,28]. However, high ArA levels have been implicated in the malpigmentation of various flatfish species [29]. An additional consideration is the possibility of DHA deficiency in neural tissues (i.e. vision) in larvae fed fish oil-based diet, as has been observed in Atlantic halibut larvae and juvenile herring eyes [30,31].
Our objective was to characterize Artemia franciscana biotopes in Colombia and to evaluate the viability of some strains for commercial exploitation. The data presented here may be of importance for the aquaculture industry to find new Artemia sources.

Results
Along the Colombian Caribbean coast, nine potential Artemia habitats were explored with variable quantities (from very few cysts in the sediment to several pounds dispersed along the pond edges) of cysts and/or biomass and variable surface area (2.5 Tayrona to 4000 ha Manaure).
Overall, Manaure: had a high number of cysts.g -1 (Table 2) and HE, but a low H% and T s . Galerazamba: had a low number of cysts.g -1 compared to other commercial cyst types, and a low H%, contrasting with a high HE and T s . Salina Cero: had a high number of cysts.g -1 , but a low H%, HE and T s . SFB (control): had a high number of cysts.g -1 , HE and H%, but a low T s .
From the three (Cl -, SO 4 2and CO 3 2-) characteristic anions used to classify hypersaline ecosystems [2], Clwas the most abundant anion (Table 3) in all locations evaluated. The physicochemical parameters monitored (Table 4) presented some tendencies inherent to each site. Salinity in Pozos Colorados and Salina Cero had a tendency to maintain low salinities (rarely crystallizing), contrasting with Manaure that presented salinities close to crystallization in the evaporation portion of the salt production circuit. Similarly, pH in Galerazamba, Salina Cero and Manaure was towards the low end pH for Artemia production; while for Tayrona and Pozos Colorados it was towards the ideal pH (8.0 to 8.5). Percent O 2 saturation was overall normal in most sites, with the exception of Tayrona which was rather low in some months. Water temperature was at the upper limit in most sites, and extremely high only in Pozos Colorados. Nitrite was overall low in all sites, contrasting with high nitrate concentration in all sites. Phosphate was also low, except in Pozos Colorados where it was too high. Primary production (chlorophyll a) was rather at the low end for hypersaline ecosystems. Precipitation was high in the southern sites explored (as expected) and low at northernmost locations (dessert-like sites).
Cyst samples from locations (Table 5) where enough cysts were collected to perform FAME analysis (FAME was actually done on freshly hatched nauplii from cysts), exhibited high EPA and ArA levels comparable to control sample (SFB-ARC1258). In contrast, most cysts collected (including SFB) at different locations, and during different months, presented low DHA levels (Manaure was the only exception with high DHA levels). Some variations in EPA and ArA levels were observed in all sites, contrasting with much lower DHA levels which remained constant for all locations, except for Manaure which exhibited variable DHA levels. DHA/EPA ratio was overall very low for all sites compared to SFB cysts. All strains had a low DHA/ ArA ratio, but a high EPA/ArA ratio, SFB included.

Discussion
Vanhaecke and Sorgeloos [32] reported cyst diameter as small as 224 µm for the San Francisco Bay strain (California, USA), while Abatzopoulos et al. [33] reported cyst diameters as large as 330 µm for the bisexual species A. tibetiana, surpassing even the well known large cyst diameters of the polyploid parthenogenetic strains with a typical diameter near 280 µm. Cysts from Great Salt Lake (GSL-Utah, USA) have a larger cyst diameter (244.2 -252.5 µm) compared with those from SFB (California, USA) (223.9 -228.7 µm) [32]; cyst diameters in the Colombian Caribbean are more similar to cysts from GSL than SFB (control). Thus, a cyst diameter grouping in the following order is possible for strains (smallest to largest): 1) San Francisco Bay (SFB-Control ARC1258, USA), 2) Galerazamba and Tayrona, 3) Kangarú, 4) Manaure, and 5) Salina Cero and Pozos Colorados. However, within the same species, strains present different cyst diameter as well as different chorion thickness. The chorion thickness for GSL ranges from 4.7 to 5.7 µm and for SFB it is 7.1 to 8.3 µm [32]. In the case of cyst samples examined in this study, Tayrona were the thinnest, followed by Salina Cero,  Note: Ionic concentrations represent one single sample per pond or salt concentration basin at the sites where Artemia was reported and must not be used for comparison purposes (between locations), since ionic concentration might vary periodically especially in managed (saltworks with several concentration levels as well as different viscosities in the system) compared to unmanaged ecosystems (single evaporation basin).
The biometric analysis grouped strains according to naupliar size (Instar I) as follows (smallest to largest): 1) Galerazamba, 2) Manaure, 3) SFB, Kangarú, and Salina Cero, 4) Pozos Colorados, and 5) Tayrona. However, it is remarkable that the naupliar length measured by Vanhaecke and Sorgeloos [32] for the Galerazamba strain collected in 1977 greatly differed from that length in our study (480 ± 31.1 vs. 390.3 ± 24.5 µm); maybe because of some physicochemical effects (i.e. salinity) on the strain manifested over time, and/or food conditions and/or cysts harvested from different pond sizes (widely separated harvest sites), which have been reported [32] to also affect cyst size and chorion thickness. The same authors reported that SFB cysts produced at 180 mg.l -1 in vitro are significantly smaller than cysts produced at lower salinities. Biometrical studies performed on several strains from different geographical origins concluded that Artemia biometrical parameters were mainly strain specific [32]. These authors revealed that the length of Instar-I nauplii (both bisexual and parthenogenetic) may vary between 430 and 520 µm. Moreover, the North American Instar I nauplii (bisexual) tend to be at the lower previously given range. Thus, naupliar length for SFB (California) is between 428 ± 28.8 to 431 ± 23.7 µm, and for GSL strain (Utah) is between 486 ± 30.6 to 489 ± 29.2 µm [32]. Naupliar size is non-critical for the feeding of crustacean larvae, which can capture and manipulate nutritional particles with their feeding appendages [18]. In contrast, prey size is very critical for fish larvae, which do not have feeding appendages and must engulf particles. The correlation between naupliar size and fish larvae mortality indicates that 20% of the larvae die of hunger when being offered nauplii greater than 480 µm in the first stages of feeding [34]. Thus, depending on the developmental stage of the cultured fish larvae, selecting an appropriate naupliar size as live feed is critical.
Overall, samples from Manaure, Galerazamba and Salina Cero, according to the cyst quality test, exhibited a relatively high number of cysts.g -1 . H% might have been affected (low in all samples evaluated) because the cyst processing method used (1/3 HP air blower with no heating element or temperature control) could not maintain a constant drying temperature/airflow, or some impurities were still present in samples. H% is dependent on degree of diapause termination, cysts energy content and amount of dead/non-viable/abortic embryos due to improper processing and/or storage [13]. Furthermore, HE reflects three factors: 1) H%, 2) presence of other components (i.e. empty shells, salt, sand, cysts water content), and 3) individual cyst weight. The low hatching synchrony of cysts from Manaure and Salina Cero could be attributed to environmental factors (i.e. raining after cysts were dehydrated, salinity, etc.), and/or as mentioned before improper processing. Significant interactions have been reported among some physicochemical-biotic factors (salinity, percent O 2 saturation and chlorophyll a) and Artemia cyst production [35], factors which might affect consequently cyst quality.
The ionic analysis of all locations registered a Clanion predominance (Table 3), as expected [36], because of their thalassohaline origin. All hypersaline environments analyzed are suitable habitats for A. franciscana development [2]. Colombian hypersaline ecosystems sampled are similar to other American Cldominant hypersaline biotopes such as Leslie saltworks (California -USA) [37], La Sal del Rey (Texas -USA) [38] and GSL (Utah -USA) [39]. In nature Artemia is found at salinity levels between 60 to 220 g.l -1 (depending on the strain and/or species) and in neutral to alkaline waters, at temperatures generally below 34°C, and at rather low O 2 levels. The low salinities measured for several months at Pozos Colorados and Salina Cero could hinder overtime Artemia production by favoring nauplii production (tolerant to low salinities and pHs) and affecting adults survival [35]. Similarly, the low end pH registered in Galerazamba, Salina Cero and Manaure could affect Artemia biomass, and even cyst production overtime. As pH decreases below 7.0 naupliar growth decreases and in adults the overall appearance deteriorates [40]. The same authors concluded that the optimum pH for Artemia growth was from 8.0 to 8.5. In the case of cyst, hatching efficiency is greatly compromised at pHs below 8.0 [41]. The low O 2 levels registered at the Tayrona site could be attributed to the high accumulation of organic matter in the pond from surrounding vegetation.
Water temperature was at the upper limit in most sites, with extreme temperatures at Pozos Colorados due to the lack of water circulation and small ponds size (< 0.5 ha). The nutrients (N:P, ideally 15:1) ratio was generally maintained within expected limits, except in Pozos Colorados where it was too high (1:1). The low primary production (chlorophyll a) determined in all sites was characteristic of hyperhaline environments. The most photosynthetically productive hypersaline environments [36] are the hypohalines and mesohalines. Furthermore, primary productivity in Salina Cero and Pozos Colorados was the highest among all sites, because of their low salinities. It is widely accepted that salinities higher than 50 g.l -1 hinder considerably primary productivity in hypersaline ecosystems perhaps because of an ionic complex formation of the dissolved macronutrients or because of a generic biologic phenomenon of a drastic specific reduction of microalgae, also occurring at higher salinities [42]. The high precipitation in the southern sites explored (Salina Cero, Galerazamba, Pozos Colorados and Tayrona) affected constantly Artemia cyst and biomass production by decreasing salinity and affecting light intensity.
The high EPA and ArA levels ( Table 5) for Colombian cysts, determined by EFA analysis, compared to the control (SFB) have great potential for the aquaculture industry. In contrast, the low DHA level content in cysts from the Colombian sites (except Manaure) might be of concern if fed to marine larvae without further DHA enrichment, since DHA deficiency affects neural tissues development [30,31], survival and growth [18], particularly at the larval stage. The very low DHA/EPA ratios (<0.5) for all Colombian sites were as expected with values often lower than 1, but this ratios could be increased up to 7 [15] through the addition of DHA rich emulsions. All strains had a low DHA/ArA ratio, contrasting with a high EPA/ArA ratio (close to dietary ratio 2:1 and 1:1 for marine fish-white bass larvae, respectively [26]). Artemia ArA storage/usage mechanism might be similar to that of other marine organisms (marine fishes) which conserve ArA even during periods of starvation [19]. The observed fluctuations in ArA and EPA levels according to FAME analysis between, and even within strains (collected in different months), may be due to year-round variations in the biochemical composition of the primary producers available to adult Artemia [15].

Conclusion
The Colombian A. franciscana habitats analyzed are of marine (thalassohaline) origin; thus, all locations were expected to be Clrich.
The relatively small nauplii are appropriate for marine larvaeculture. In contrast, the strains from Tayrona, Kangarú, Salina Cero, and Pozos Colorados may be of use but limited to Artemia small biomass production quantities, because of the small surface area of their respective loca- tions; Artemia could be exploited at these locations for local aquaculture applications.
In general, cyst quality evaluation for Manaure, Salina Cero and Galerazamba strains determined that cysts from these three locations could improve their quality by concentrating efforts on cyst processing techniques. Further, cyst quality might have been affected by interactions among some physicochemical-biotic factors and Artemia cyst production in their ecosystem which could be improved by managing some of these key physicochemical-biotic factors and/or infrastructure management (e.g. brine concentration in the different basins and nutrients).
EFA profiles demonstrated that Tayrona, Galerazamba, Manaure and Salina Cero strains are suitable for marine aquaculture because of their high EPA/ArA ratio, but might have to be fortified with DHA rich emulsions depending on the nutritional requirements of the species to be cultured, because of their overall low DHA content.
Kangarú (11° 59' 28''N, 74° 32' 21''W): is less than a 4 ha natural thalassohaline saltwork comprised of three small ponds located in the northern region of the Salamanca Island National Natural Park, Magdalena department (Fig. 1). It was explored in July 2000. Salt has been occasionally exploited for decades. This locality, is an important bird migration spot, however, it lost importance because of mangrove destruction as consequence of building a highway through the park.
Pozos Colorados (11° 09' 45''N, 74° 13' 34''W): is an approximately 65 ha very old artificial thalassohaline saltwork, currently abandoned. Few studies have been conducted by local researchers in the past using samples from this location. It is located near Santa Marta city, Magdalena department (Fig. 1), contiguous to the road connecting to Barranquilla to Santa Marta city [11]. This saltwork consists of only five irregularly shaped, shallow ponds with only 4 ha water surface.
Tayrona NNP encompasses a small number of saline non crystallizing ponds, with the exception of Chengue, where Artemia has been reported to occur. The salt pond is flooded during most of the year and serves as a saltwork during summer [57]. Chengue Inlet, is located in the middle of the Tayrona NNP, it presents a series of small bays and inlets extending from Santa Marta to Cañaverales to the east. Chengue salt exploitation existed long before the prehispanic period [58].
Manaure (11° 46' 32''N, 72° 29' 27''W): is located to the west, contiguous to the town of Manaure, in the center of La Guajira department, near Riohacha city (Fig. 1). Studies have been conducted by several authors in the past using samples from this location [5,[11][12][13]35,43,47,48,50,52,54,55,59]. This saltwork is a thalassohaline, shallow water body extending over 4,000 ha. Water movement through the saltwork system is achieved both by pumping and through gravity. There are six pumping stations that increase water volume to a predetermined water level, thereafter water will flow by gravity. This zone was originally a natural lagoon surrounded in some areas by mangroves. The deposits were constructed using the natural topography of the terrain with some modifications. The levees were built by compacting large amounts of clay material brought from the margins of the saltwork [47].
Warrego (12° 19'N, 71° 54'W): is an approximately 600 ha (2 miles long) thalassohaline saltern located in the northern tip of La Guajira department, near Puerto Nuevo village (Fig. 1). Occasionally, the Wayu Indians extract salt when the brine crystallizes. Since it was completely dried up when we visited it (January 18,2000), no water samples were collected from this location and found few Artemia cysts.

Preparation and sampling
Sampling was conducted monthly and cysts batches were collected irregularly (whenever available) in nine thalassohaline locations aforementioned in the northern region of the Colombian Caribbean, from July 1998 to June 2000. Cyst processing was done following these steps: (i) size separation with brine, (ii) density separation in brine, (iii) washing in freshwater, (iv) density separation in freshwater, (v) drying below 40°C, and (vi) vaccum packing and refrigerating cyst at 4 ± 2°C.
Cyst diameter and chorion thickness were recorded from sites where sufficient cysts were collected, using SFB (USA, ARC1258) cysts as reference material. Cysts were incubated for 3 hr in 10 g.l -1 artificial sea water (Instant Ocean ® ) at 25 ± 0.5°C and pH 8.3 [13]. One percent lugol's solution (5 %) was added to the sea water to stop embryos from hatching and cysts were in the dark overnight. Cyst diameter (µm) was measured in 200 cysts with a precalibrated microscope. Mean value and standard deviation were calculated using the predetermined conversion factor. Decapsulated cyst diameter (µm): a small sample of cysts was hydrated in tapwater for 2 h. Cysts were then decapsulated with a NaOH and NaOCl solution. Cysts were rinsed well and incubated in 10 g.l -1 artificial sea water (Instant Ocean ® ) with 1 % lugol for 1 hr, at 25 ± 0.5°C, and pH 8.3 and was incubated for 1 h more. Afterwards, 1 % lugol was added again to the incubating solution and cysts were left overnight in the dark.
Decapsulated cyst diameter was measured for 200 cysts with a precalibrated microscope. Mean value and standard deviation were calculated using the predetermined conversion factor. Chorion thickness was calculated using this formula: (cyst diameter -decapsulated cyst diameter)/ 2 Naupliar length was determined on Instar I nauplii, following this procedure [13]: cysts were incubated and hatched under controlled conditions (25 ± 0.5°C, pH 8.3 and illumination: 1000 lux) in artificial sea water (Instant Ocean ® ) at 35 g.l -1 salinity [32]. Nauplii were sampled at Instar I considering the protopodite of each antennae which bears two endites with a single long bristle attached to each and their brownish-orange color due to yolk presence (Instar II is translucent) as the traits defining this stage [60,61]. Nauplii were harvested when 90% of the total number of hatchable nauplii had been produced [22]. Two hundred nauplii were fixed in lugol's solution (5%) and the length determined using a microscope with a pre-calibrated projection system. Cyst quality studies [13] were performed only on major saltworks. The following parameters were used to evaluate cyst quality: i) Hatching percentage (H%): number of nauplii that can be produced under standard hatching conditions from 100 full cysts (with embryos). H% = (N × 100)(N+U+E) -1 Cysts (1.6 g) were incubated in 800 ml of 32 g.l -1 microfiltered (<1 µm) seawater (Instant Ocean ® ) under continuous illumination (2000 lux) at 28°C, pH = 8.3, in a cylindroconical vessel (test was run in triplicate per strain) with bottom aeration (>2 mg.l -1 ). Vessels were suspended in a water bath in a 100 gal aquarium with a water heather and a mixer to maintain a well distributed temperature (± 1°C). After 24 h incubation six 250 µl subsamples were taken from each cone with a micropipet. Each subsample was pipetted into a small vial and nauplii were fixated by adding a few drops of lugol's solution (5%). Nauplii (n i ) and umbrella (u i ) stages were counted in each subsample under a disection microscope. Mean values (N = nauplii and U = umbrella) were calculated each for these two stages. Unhatched cysts were decapsulated and empty cyst shells were dissolved with a drop of NaOH solution (40 g.100 ml -1 distilled water) and five drops NaOCl (5.25% NaOCl) added to each vial. Unhatched (orange colored) embryos (e i ) were counted per cone (i = 6) and mean value (E) was calculated for each cone. H% value was calculated per cone, and mean value and standard deviation was calculated for three cones (final H% value).
ii) Hatching efficiency (HE): number of nauplii/g dry cysts that can be produced under standard hatching conditions. HE = (N × 4 × 800 ml)(1.6 g) -1 HE value was calculated, for each strain evaluated, per cone, and mean value and standard deviation was calculated for three cones (final HE value). Hatching vessels were left for another 24 h, subsequently subsamples were again taken to calculate H% and HE for 48 h incubation period.
iii) Hatching rate (HR): period from incubation (cyst hydration) to nauplii release (hatching). The following HR time intervals are considered: T 0 = Incubation time untill appearance of first free swimming nauplii T 10 = Incubation time untill appearance of 10 % of total hatchable nauplii T 90 = Incubation time untill appearance of 90 % of total hatchable nauplii T S = T 90 -T 10 ; this value gives an indication of the hatching synchrony Six 250 µl samples, for each strain evaluated, were taken 12 h after incubation and HE was calculated every 3 h until HE mean value remained constant for three consecutive sampling periods. Mean values per period were then expressed as percentages of the maximal HE. A hatching curve was plotted for each strain, and T 10 and T 90 were extrapolated from the graph. iv) Number of cysts.g -1 : this parameter is dependent on cysts diameter. Cysts (4 g) were placed in an aluminum plate and weighted and dried in a drying oven set at 60°C for 24 h. Cysts were then cooled down to room temperature for 4 h in a tightly sealed glass drying chamber with fresh desiccant. One g of cyst sample (triplicates) was then weighted (0.1 mg accuracy) in an aluminum plate; to determine average cyst weight, for a single cyst, ten subsamples (10.0 ± 0.1 mg) were taken from each replicate and counted; to find how many cysts were in the 1 g sample average cyst number in 10 mg was then extrapo-lated to the 1 g cyst sample. This procedure was repeated for each strain evaluated.
At each location we measured: salinity (temperature compensated refractometer), percent O 2 saturation and temperature (Oxymeter WTW ® 330), pH (pH meter WTW ® 330), nitrates, nitrites and phosphates (Hatch ® DREL 2010 spectrophotometer), and chlorophyll a. For the latter, we used the Seston method and read (Hatch ® DREL 2010 spectrophotometer). The ionic composition (Table  3) was determined using a Unicam 939/959 atomic absorption spectrophotometer. All samples were diluted with deionized water because of the high ionic concentration. A sample of nauplii from Galerazamba, Manaure, Salina Cero and Tayrona was taken for FAME; for this analysis we followed Sorgeloos et al. [13]: cysts were incubated and hatched under controlled conditions (25 ± 0.5°C, pH 8.3 and illumination: 1000 lux) in artificial sea water (Instant Ocean ® ) at 32 g.l -1 salinity. FAME methodology for freshly hatched nauplii (0.25 g) was a modification of the direct esterification described by Lepage and Roy [62]. The latter implicates a direct acid catalized transesterification without prior extraction of total fat, on dry sample (triplicates) amounts ranging from 10 to 150 mg. An internal standard 20:2 (n-6) was added prior to the reaction. FAME were extracted with hexane. After solvent evaporation FAME were prepared for injection by redissolving them in iso-octane (2 mg/ml). Quantitative determination was done by a Chrompack CP9001 gas chromatograph equipped with an autosampler and a TPOCI (Temperature programmable on-column injector). Injections (0.5 µl) were performed on column into a polar 50 m capillary column, BPX70 (SGE Australia), with a diameter of 0.32 mm and a layer thickness of 0.25 µm, connected to a 2.5 m methyl deactivated precolumn. The carrier gas was H 2 , at a pressure of 100 kPa and the detection mode FID. The oven was programmed to rise from the initial temperature of 85 to 150°C at a rate of 30°C/ min, from 150 to 152°C at 0.1°C/min, from 152 to 172°C at 0.65°C/min, from 172 to 187°C at 25°C/min and to stay at 187°C for 7 min. The injector was heated from 85 to 190°C at 5°C/sec and stayed at 190°C for 30 min. Identification was based on standard reference mixtures (Nu-Chek-Prep, Inc., USA). Integration and calculations were done on computer with a software program Maestro (Chrompack).
Any experimental research on animals that was reported in this study was performed with the approval of an appropriate ethics committee regulating animal research.

Calculations and statistics
Standard deviations were calculated for all cyst diameter and naupliar length measurements. Data obtained were