RNA/DNA and derived condition indices for anchovy and hake larvae as relevant information for comprehensive fisheries management

INTRODUCTION

The Argentine anchovy, Engraulis anchoita (Hubbs and Marini, 1935), represents the pelagic resource with the highest biomass in the southwest Atlantic Ocean, with a total catch of approximately 9000 t reported in 2021 (Ministry of Agroindustry 2022Ministry of Agroindustry. 2022. https://www.magyp.gob.ar/sitio/areas/pesca_maritima/desembarques/ ). Its high densities and wide distribution make it the most ecologically important fish in the Argentine Sea. It plays a key role in the food web, sustaining several species of commercial value within this region (Angelescu 1982Angelescu V. 1982. “Ecología trófica de la anchoíta del Mar Argentino (Engraulidae, Engraulis anchoita). Parte II. Alimentación, comportamiento y relaciones tróficas en el ecosistema.” Contr. INIDEP. 409: 1-83., Hansen 2004Hansen J.E. 2004. Anchoíta (Engraulis anchoita). In: Sánchez R.P., Bezzi S.I., Boschi E.E. (eds), El mar argentino y sus recursos pesqueros. Publicaciones especiales INIDEP: 101-115.). Two populations of anchovy are separated by the 41°S parallel (Sánchez and Ciechomski 1995Sánchez R.P., Ciechomski J.D. 1995. Spawning and nursery grounds of pelagic fish species in the sea-shelf off Argentina and adjacent areas. Sci. Mar. 59: 455-478). The spawning and breeding areas of this species are related to upwelling regions, estuarine, tidal and shelf break fronts characteristic of the Argentine continental shelf. The reproductive activity of the Patagonian stock, distributed south of 41°S, begins in November and shows maximum spawning in December (Sánchez and Ciechomski 1995Sánchez R.P., Ciechomski J.D. 1995. Spawning and nursery grounds of pelagic fish species in the sea-shelf off Argentina and adjacent areas. Sci. Mar. 59: 455-478). The Argentine hake Merluccius hubbsi is one of the most important fishing resources for the Argentine bottom trawling fleet, with a total catch of approximately 287000 t reported in 2021 (Ministry of Agroindustry 2022Ministry of Agroindustry. 2022. https://www.magyp.gob.ar/sitio/areas/pesca_maritima/desembarques/ ). It inhabits the waters of the southwest Atlantic Ocean between Cabo Frio in Brazil (22°S) and southern Argentina (55°S), at depths between 50 and 500 m (Cousseau and Perrota 1998Cousseau M.B., Perrota R.G. 1998. Peces marinos de Argentina. Biología, distribución y pesca. INIDEP, Mar del Plata, 163 pp.). There are two main fishing stocks in Argentina, north, and south of 41°S. Owing to the increase in fishing pressure in recent decades, a decrease in spawning biomass has been observed in this species, as well as variations in the age structure of the parental stock and the location of spawning schools (Macchi et al. 2005Macchi G.J., Pájaro M., Madirolas A. 2005. Can a change in the spawning pattern of Argentine hake (Merluccius hubbsi) affect its recruitment? Fish. Bull. 103: 445-452., 2021Macchi G.J., Diaz M.V., Leonarduzzi E., et al. 2021. Temperature, maternal effects and density-dependent processes during early life stages of Argentine hake as relevant recruitment drivers. Fish. Res. 238: 105898. https://doi.org/10.1016/j.fishres.2021.105898 ). The reproductive activity of the southern or Patagonian stock, which is the one with the highest population abundance, occurs on the North Patagonian shelf mainly between November and April and peaks in January (Ehrlich 1998Ehrlich M.D. 1998. Los primeros estadios de vida de la merluza Merluccius hubbsi, Marini 1933, en el Mar Argentino como aporte al conocimiento de su reclutamiento y estructura poblacional. Doctoral thesis, Univ. Buenos Aires, 318 pp., Macchi et al. 2004Macchi G.J., Pájaro M., Ehrlich M. 2004. Seasonal egg production pattern of the Patagonian stock of Argentine hake (Merluccius hubbsi). Fish. Res. 67: 25-38. https://doi.org/10.1016/j.fishres.2003.08.006 , Macchi et al. 2007Macchi G.J., Pájaro M., Dato C. 2007. Spatial variations of the Argentine hake (Merluccius hubbsi) spawning shoals in the Patagonian area during a reproductive season. Rev. Biol. Mar. Oceanogr. (Chile). 42: 345-356. https://doi.org/10.4067/S0718-19572007000300013 ).

As in other fish species, variability in the recruitment of both anchovy and hake is affected by processes that operate on different spatial and temporal scales. This variability depends on physical and thermodynamic factors that determine survival during the early stages of life (Houde 2008Houde E.D. 2008. Emerging from Hjort’s Shadow. J. Northwest Atl. Fish. Sci. 41: 53-70. https://doi.org/10.2960/J.v41.m634 ). The study of larval nutritional condition allows us to evaluate the individual physiological state, which at the same time reflects the environmental context to which the larvae have been exposed (Chícharo and Chícharo 2008Chícharo M.A., Chícharo L. 2008. RNA: DNA ratio and other nucleic acid derived indices in marine ecology. Int. J. Mol. Sci. 9: 1453-1471. https://doi.org/10.3390/ijms9081453 ). It is also a useful instrument for determining favourable breeding areas (Diaz and Pájaro 2012Diaz M.V., Pájaro M. 2012. Protocolo para la determinación de las concentraciones de ácidos nucleicos en larvas de peces. Inf. Inv. INIDEP N°20, Mar del Plata, 9 pp.). Various nutritional condition indices have been used to estimate mortality due to starvation in marine fish larvae (Buckley 1984Buckley L. 1984. RNA/DNA ratio: an index of larval fish growth in the sea. Mar. Biol. 80: 291-298. https://doi.org/10.1007/BF00392824 , Clemmesen 1994Clemmesen C. 1994. The effect of food availability, age or size on the RNA/DNA ratio of individually measured herring larvae: laboratory calibration. Mar. Biol. 118: 377-382. https://doi.org/10.1007/BF00350294 ). The RNA/DNA ratio (RD) stands as one of the best indicators of the nutritional status of various marine organisms (Clemmesen 1994Clemmesen C. 1994. The effect of food availability, age or size on the RNA/DNA ratio of individually measured herring larvae: laboratory calibration. Mar. Biol. 118: 377-382. https://doi.org/10.1007/BF00350294 , Folkvord et al. 1996Folkvord A., Ystanes L., Moksness E. 1996. RNA:DNA ratios and growth of herring (Clupea harengus) larvae reared in mesocosms. Mar. Biol. 126: 591-602. https://doi.org/10.1007/BF00351326 ) and is currently the biochemical index most widely used as an indicator of the nutritional condition of fish larvae (Chícharo and Chícharo 2008Chícharo M.A., Chícharo L. 2008. RNA: DNA ratio and other nucleic acid derived indices in marine ecology. Int. J. Mol. Sci. 9: 1453-1471. https://doi.org/10.3390/ijms9081453 ). The RD ratio varies with age, developmental stage and size under different environmental conditions (Bulow 1970Bulow F.J. 1970. RNA-DNA ratios as indicators of recent growth rates of a fish. J. Fish. Res. Board Can. 27: 2343-2349. https://doi.org/10.1139/f70-262 ). It has also been shown to respond to changes in the concentration of available prey (McGurk et al. 1992McGurk M.D., Warburton H.D., Galbraith M., Kusser W.C. 1992. RNA-DNA ratio of herring and sand lance larvae from Port Moller, Alaska: Comparison with prey concentration and temperature. Fish. Oceanogr. 1: 193-207. https://doi.org/10.1111/j.1365-2419.1992.tb00038.x , Chícharo and Chícharo 1995Chícharo L., Chícharo M.A. 1995. The RNA/DNA ratio as a useful indicator of the nutritional condition in juveniles of Ruditapes decussatus. Sci. Mar. 59 (suppl. 1): 95-101.) among other factors. The monitoring of the larval state in situ over time could be a useful tool for determining favourable breeding areas for the species and for developing a time series to assess the effects of climate change in these areas. Detecting these favourable areas and periods for larval survival makes a valuable contribution to the comprehensive management of a population subjected to fishing exploitation (Viladrich et al. 2016Viladrich N., Rossi S., Lopez-Sanz A., Orejas C. 2016. Nutritional condition of two coastal rocky fishes and the potential role of a marine protected area. Mar. Ecol. 37: 46-63. https://doi.org/10.1111/maec.12247 ).

The northern Patagonia region is hydrographically characterized by the existence of a tidal front, the Northern Patagonian Frontal System or NPFS (Guerrero et al. 1997Guerrero R.A., Acha E.M., Framiñan M.B., Lasta C.A. 1997. Physical oceanography of the Río de la Plata Estuary, Argentina. Cont. Shelf Res. 17: 727-742. https://doi.org/10.1016/S0278-4343(96)00061-1 , Martos and Sánchez 1997Martos P., Sánchez R. 1997. Caracterización oceanográfica de regiones frontales en la plataforma patagónica en relación con áreas de desove y cría de la anchoíta (Engraulis anchoita). 10° Coloquio Argentino de Oceanografía, 4-5 Septiembre, IADO-CONICET, Bahía Blanca., Sabatini and Martos 2002Sabatini M.E., Martos P. 2002. Mesozooplankton features in a frontal area off northern Patagonia (Argentina) during spring 1995 and 1998. Sci. Mar. 66: 215-232. https://doi.org/10.3989/scimar.2002.66n3215 ). This system is characterized by the formation of a seasonal thermocline, particularly during the austral summer, which gives rise to a homogeneous coastal zone and an increasing stratification towards the offshore zone. Figure 1 depicts three typical vertical temperature profiles observed in this area during summer: coastal homogeneous stations (Fig. 1D-F), frontal stratified stations (Fig. 1G) and stratified offshore stations (Fig. 1H-I). The dynamic that characterizes the frontal system causes a high availability of nutrients that are mainly due to upwelling and concentration and favour primary and secondary productivity (Bakun and Parrish 1991Bakun A., Parrish R.H. 1991. Comparative studies of coastal pelagic fish reproductive habitats: the anchovy (Engraulis anchoita) of the Southwestern Atlantic. ICES J. Mar. Sci. 48: 343-361. https://doi.org/10.1093/icesjms/48.3.343 , Bakun 1997Bakun A. 1997. Patterns in the ocean: ocean processes and marine population dynamics. Oceanogr. Lit. Rev. 5: 530.), generating major phytoplankton blooms (Carreto and Benavídez 1990Carreto J.I., Benavídez H.R. 1990. Synopsis on the reproductive biology and early life of Engraulis anchoita, and related environmental conditions in Argentine waters. Phytoplankton. IOC. Worksh. Rep. 65. Annex V: 2-5.) and large aggregations of copepods (Derisio et al. 2014Derisio C., Alemany D., Acha E.M., Mianzan H.W. 2014. Influence of a tidal front on zooplankton abundance, assemblages and life histories in Península Valdés, Argentina. J. Mar. Syst. 139: 475-482. https://doi.org/10.1016/j.jmarsys.2014.08.019 , Temperoni et al. 2014Temperoni B., Viñas M.D., Martos P., Marrari M. 2014. Spatial patterns of copepod biodiversity in relation to a tidal front system in the main spawning and nursery area of the Argentine hake Merluccius hubbsi. J. Mar. Syst. 139: 433-445. https://doi.org/10.1016/j.jmarsys.2014.08.015 , Temperoni and Viñas 2015Temperoni B., Viñas M.D. 2015. Disponibilidad de presas zooplanctónicas para larvas de Merluccius hubbsi en el área de desove. Resultados de la campaña EH-01/14. Inf. Inv. INIDEP N°50, Mar del Plata, 15 pp.), in addition to a great diversity of other zooplanktonic organisms (Mianzan and Guerrero 2000Mianzan H.W., Guerrero R.A. 2000. Environmental patterns and biomass distribution of gelatinous macrozooplankton. Three study cases in the Southwestern Atlantic Ocean. Sci. Mar. 64: 215-224. https://doi.org/10.3989/scimar.2000.64s1215 , Schiariti 2008Schiariti A. 2008. Historia de vida y dinámica de poblaciones de Lychnorhyza lucerna (Scyphozoa) ¿un recurso pesquero alternativo? Doctoral thesis, Univ. Buenos Aires, 220 pp., Schiariti et al. 2015Schiariti A., Betti P., Dato C., et al. 2015. Medusas y ctenóforos de la región norpatagónica I: diversidad y patrones de distribución. Inf. Inv. INIDEP N°21, Mar del Plata, 16 pp.). The water circulation in this area also favours the retention of the first stages of life of both fish and invertebrates (Álvarez Colombo et al. 2011Álvarez-Colombo G., Dato C., Macchi G., et al. 2011. Distribution and behavior of Argentine hake larvae: Evidence of a biophysical mechanism for self-recruitment in northern Patagonian shelf waters. Cienc. Mar. 37, 633-657. https://doi.org/10.7773/cm.v37i4B.1777 ). The NPFS is therefore a propitious area for the reproduction and breeding of many species during spring and summer, including important fishing resources such as the Patagonian stocks of hake and anchovy (Hansen et al. 2001Hansen J.E., Martos P., Madirolas A. 2001. Relationship between spatial distribution of the Patagonian stock of Argentine anchovy, Engraulis anchoita, and sea temperatures during late spring to early summer. Fish. Oceanogr. 10: 193-206. https://doi.org/10.1046/j.1365-2419.2001.00166.x , Macchi et al. 2004Macchi G.J., Pájaro M., Ehrlich M. 2004. Seasonal egg production pattern of the Patagonian stock of Argentine hake (Merluccius hubbsi). Fish. Res. 67: 25-38. https://doi.org/10.1016/j.fishres.2003.08.006 , Pájaro et al. 2005Pájaro M., Macchi G.J., Martos P. 2005. Reproductive pattern of the Patagonian stock of Argentine hake (Merluccius hubbsi). Fish. Res. 72: 97-108. https://doi.org/10.1016/j.fishres.2004.09.006 ). Given that the survival during the first phases of life of these species is affected by the existence of the NPFS, the spatial and temporal coincidence of these organisms with this tidal front during the larval stage could be one of the most important factors that explain the observed variability in recruitment. Because of the importance of the NPFS for the conservation and management of major fishery resources and the aforementioned hydrographic characteristics, part of this area is currently under study for the implementation of a new marine protected area (MPA). For this reason, it is essential to carry out a comprehensive analysis of the NPFS, covering the environmental and biological aspects of this region.

In December 2018, a research survey was carried out to analyse the oceanographic and biological conditions in the NPFS region within the project “Strengthening the Management and Protection of Coastal Biodiversity Marine in Key Ecological Areas and the Application of the Ecosystem Approach to Fisheries” (GEF-FAO). Within the framework of this survey, our main objective was to evaluate the nutritional condition of E. anchoita and M. hubbsi larvae during the beginning of the austral summer in the northern Patagonia region and to analyse the spatial variation of this parameter. Because the sampling was performed early in the spawning period for hake, a small number of larvae of this species were collected. Thus, the study is mainly based on anchovy larvae, with some additional information on the nutritional condition of initial larval stages of hake. We hypothesized that the presence of the NPFS provides adequate features for larval growth and survival evidenced in a better nutritional condition of anchovy and hake larvae than that observed in the offshore area. To test this hypothesis, we assessed the RD index and a derived index of growth performance. The proportion of individuals in starvation and in optimal growth conditions was estimated. Larval condition was mapped in the study area to determine the existence of favourable areas for growth and larval survival in relation to the NPFS.

MATERIALS AND METHODS

Sample collection

 ⌅

The samples were taken during a research survey carried out on board the vessel Victor Angelescu of the National Institute for Fisheries and Development (INIDEP) in the northern Patagonia area between 4 and 16 December 2018 (Fig. 1). At each station, temperature and salinity measurements were made using the CTD Seabird Electronics system. The data collected by the staff of the Physical Oceanography Office of INIDEP are part of the BaRDO database. A stratified plankton sampling was performed with the HydroBios Model Midi Multinet (0.5×0.5 m), equipped with three opening and closing nets (300 µm pore size) and soft collectors. Oblique trawls were carried out at different levels in the water column. The towing speed during the ascent was maintained between 2 and 3 knots, with a duration that varied between 5 and 7 minutes per sampling strata. When the station was in the homogeneous area the nets were operated covering strata of equal width, while in the frontal and stratified regions the nets were operated above, on and below the thermocline. Five transects or oceanographic sections were made (Fig. 1: S1-S5). The stations of sections S2 and S3 were grouped for the analysis in relation to the position of the frontal system: the stations in the homogeneous or internal zone (Fig. 1: S2i, S3i) and the stations of the stratified offshore or external zone (Fig. 1: S2e, S3e). Larvae were not found in section S5.

Fig. 1.  Spatial distribution of the sampling stations. (A) Five transects or oceanographic sections (S1-S5) were made. The schematic position of the Northern Patagonian Frontal System (NPFS) during December 2018 is indicated in blue (calculated according to Simpson 1981Simpson J.H., Bowers D. 1981. Models of stratification and frontal movement in shelf seas. Deep-Sea Res. 28: 727-738. https://doi.org/10.1016/0198-0149(81)90132-1 ; by Martos and collaborators, INIDEP Physical Oceanography Office). For the analysis, the stations of sections S2 and S3 were grouped in relation to the position of the NPFS: the stations in the homogeneous zone of the front or internal zone (S2i, S3i) and the stations of the stratified offshore or external zone (S2e, S3e). (B) Horizontal isolines of temperature (ºC) at the surface and (C) at the bottom are indicated. The vertical gradients of temperature (ºC) with respect to depth (m) are indicated in the stations with an “x” in Section 2 (S2), (D-I) from coast to offshore stations. Data obtained from the BaRDO-INIDEP regional oceanographic database. Dotted lines indicate 50 and 100 m isobaths.

Once the plankton samples (N=126) were obtained, they were inspected on board to detect and separate anchovy and hake larvae. These were extracted from the sample and placed in labelled cryotubes and then stored in an ultrafreezer (-80°C) for studies of nutritional condition and growth. A representative larval sample was taken at each sampling station with a maximum of 100 larvae in those where anchovy larvae were very abundant, including all sizes present in the entire sample. In total, 1045 anchovy larvae were collected from 29 stations, and 16 hake larvae were obtained from four hauls. The rest of the plankton sample was fixed in 5% formalin in seawater to be later analysed under a binocular stereoscope in the INIDEP laboratories.

Data analysis

 ⌅

The RD values obtained in this study were standardized (RDs) according to Caldarone et al. (2006)Caldarone E.M., Clemmesen C.M., Berdalet E., et al. 2006. Intercalibration of four spectrofluorometric protocols for measuring RNA/DNA ratios in larval and juvenile fish. Limnol. Oceanogr-Meth. 4: 153-163. https://doi.org/10.4319/lom.2006.4.153 , using 2.4 as the reference value for the slope of the calibration curves of ultrapure DNA and RNA standards. This procedure allows direct comparison with other published RDs results, avoiding inter-laboratory differences caused by analytical protocols. The average value for the slope of the calibration curves obtained in the present study was 2.19 (±0.71 SD).

Growth rate (G) was estimated for each larva using the RDs-T-G model developed by Buckley et al. (2008)Buckley B.A., Caldarone E.M., Clemmesen C.M. 2008. Multi-species larval fish growth model based on temperature and fluorometrically derived RNA/DNA ratios: results from a meta-analysis. Mar. Ecol. Prog. Ser. 371: 221-232. https://doi.org/10.3354/meps07648 . For hake larvae, the model developed by Buckley for gadiform fish was used, and the general multi-specific model was used to determine the instantaneous growth of anchovy larvae according to the following expressions:

$\mathrm{G}=0.0254\mathrm{x}\mathrm{R}\mathrm{D}\mathrm{s}+0.0037\mathrm{x}\mathrm{T}\mathrm{x}\mathrm{R}\mathrm{D}\mathrm{s}-0.0873\mathrm{ }\left(\mathrm{h}\mathrm{a}\mathrm{k}\mathrm{e}\right)$  (1)

$\mathrm{G}=0.0145\mathrm{x}\mathrm{R}\mathrm{D}\mathrm{s}+0.0044\mathrm{x}\mathrm{T}\mathrm{x}\mathrm{R}\mathrm{D}\mathrm{s}-0.0780\mathrm{ }\left(\mathrm{a}\mathrm{n}\mathrm{c}\mathrm{h}\mathrm{o}\mathrm{v}\mathrm{y}\right)$  (2)

where G is the instantaneous growth rate and T is the temperature measured at the sampling site, which corresponds to the average of the stratum in which the larvae were collected.

A critical value of the RD index was determined assuming null larval growth (CRD), that is G=0. The percentage of larvae below the CRD was calculated, and it was assumed that these specimens were in starvation (% starvation, Fig. 2A).

In addition, the growth performance (Gpf) was determined as a derived condition index. This index represents the quotient of the observed G and a reference growth (Gref) rate achieved by a larva under optimal environmental and feeding conditions. Due to the lack of a Gref for the studied species, larval growth rates were compared with a Gref that was calculated according to Houde and Zastrow (1993)Houde E.D., Zastrow C.E. 1993. Ecosystem- and taxon-especific dynamic and energetics properties of larval fish assemblages. Bull. Mar. Sci. 53: 290-335., who established a multi-specific model based on 80 marine and estuarine species:

$\mathrm{G}\mathrm{r}\mathrm{e}\mathrm{f}=0.0106\mathrm{x}\mathrm{T}-\mathrm{ }0.0203$  (3)

Larvae with Gpf higher than or equal to 1 were assumed to be in optimal growth condition. The percentage of larvae with Gpf higher than 1 (optimal %, Fig. 2B) was calculated. A Student t-test was used to compare mean RD values of pre-flexion and flexion anchovy larvae. ANOVA was used to compare mean RD values of each size class of anchovy larvae. Analyses were performed to evaluate variability in the position of the frontal system (calculated according to Simpson (1981)Simpson J.H., Bowers D. 1981. Models of stratification and frontal movement in shelf seas. Deep-Sea Res. 28: 727-738. https://doi.org/10.1016/0198-0149(81)90132-1 , see Fig. 1), time of the day and depth where the larvae were collected. The mean values of RD and Gpf obtained for the larvae were compared by ANCOVA, using the standard length of the specimens as a co-variable. When significant differences were found, a Tukey test was performed. The results are then expressed according to oceanographic section (S1-S4), oceanographic section grouping the internal and external stations in relation to the position of the tidal front (S1, S2i, S2e, S3i, S3e, S4), time of day (M, morning; A, afternoon; N, night), and depth (B, bottom; T, thermocline; S, surface).

Fig. 2.  Theoretical distribution for (A) RNA/DNA index (RD) and (B) growth Pperformance (Gpf) values for fish larvae. (A) The curve corresponding to RD shows a critical value (CRD) used to detect starving larvae: above this value, the larvae are not in starvation (white colour) and below it they are in starvation (red colour). (B) The Gpf curve represents the growth performance being the quotient between the observed growth rate and the larval growth rate under optimal environmental and feeding conditions. Values greater than 1 (green colour) indicate more than optimal growth, and values less than 1 (white colour) indicate less than optimal growth. Modified from Alves et al. (2022)Alves N.M., Braverman M.S., Temperoni B. et al. 2022. Primeros estudios de condición nutricional en juveniles de Micropogonias furnieri en Bahía Samborombón durante dos temporadas del año (cálida y fría). Inf. Inv. INIDEP N°67, Mar del Plata, 18 pp..

RESULTS

A large number of the anchovy larvae obtained (N=959) were classified in pre-flexion (SL<8 mm) and in flexion (8-12.9 mm) and a small proportion in post-flexion (SL>12.9 mm). The RD index showed a significant increase throughout the anchovy larval ontogeny, with an average value of 1.84±1.39 (N=739) and 2.77±1.50 (N=220) in pre-flexion and flexion, respectively (t-test, T (N= 957) = 8.54; p<0.0001). A small number of Merluccius hubbsi larvae were collected (N=15) and classified in pre-flexion (SL<6.49 mm), and the RD obtained was 1.64±0.55.

The mean values ​​of the RD index for size class showed a positive trend in both species. For anchovy larvae, the mean RD by size class also showed an increase towards the post-flexion stage (Fig. 3A), but no significant differences were observed when the mean values ​​were compared by size class, except for the class 14 mm (ANOVA, F (11, 947) = 15.47; p<0.0001); however, only one larva of this size class was collected.

The scatterplot of the RD indices as a function of the standard length (SL, mm) of the anchovy larvae shows the proportion of individuals below the critical RD index (28%), that is, in starvation (Fig. 3B), and the proportion of individuals with a growth performance (Gpf) above one (26%), that is, in optimal growth (Fig. 3C).

Significant differences were observed in the mean RD values in the oceanographic sections for anchovy larvae. A co-variance analysis was performed using the SL of the specimens as a co-variable (Fig. 4, Table 1), and differences were observed when the stations were grouped according to the position of the NPFS. The larvae from the internal stations of the frontal system showed higher condition indices than the larvae from the external stations. (Fig. 4 and 5A-B, Table 2). A better nutritional condition was observed in the larvae collected in the upper stratum (surface), but no differences were observed at the different times of the day (Fig. 4, Table 3).

Finally, the percentage of larvae in starvation (RDs<CRDs) and the percentage of specimens in optimal growth (Gpf>1) were mapped in both species. A higher percentage of anchovy larvae in starvation was observed at the stations outside the position of the NPFS, and a higher percentage of larvae in optimal growth was detected at the internal stations (Fig. 5C-F, Table 4). It was observed that between 18% and 32% of anchovy larvae were in starvation and between 24% and 44% ​​above optimal growth inside the NPFS. At offshore stations between 59% and 66% of the anchovy larvae were in starvation and the percentage above optimal growth was low (0% to 6%). Of the hake larvae, it was observed that 25% were in starvation and 12.5% ​​above optimal growth.

DISCUSSION

Little information is available on the nutritional condition of the anchovy southern stock, so the results presented here are highly important. Somewhat lower RD values ​​ were observed than those previously reported in the study area (Diaz et al. 2016Diaz M.V., Do Souto M., Peralta M., et al. 2016. Comer o ser comido: factores que determinan la condición nutricional de larvas de Engraulis anchoita de la población patagónica de la especie. Ecología Austral 26:120-133. https://doi.org/10.25260/EA.16.26.2.0.71 , Do Souto et al. 2019Do Souto M., Brown D.R., Leonarduzzi E., et al. 2019. Nutritional condition and otolith growth of Engraulis anchoita larvae: the comparison of two life traits indexes. J. Mar. Syst. 193: 94-102. https://doi.org/10.1016/j.jmarsys.2019.01.008 ), and a large percentage of anchovy larvae were recorded below the critical value of the RD index in the zone external to the NPFS. The internal area of the NPFS showed a lower proportion of individuals under starvation and a higher incidence of individuals showing optimal growth.

Because the survey research was carried out during the beginning of the reproductive period of the Argentine hake, the number of hake larvae analysed was very low (Macchi et al. 2004Macchi G.J., Pájaro M., Ehrlich M. 2004. Seasonal egg production pattern of the Patagonian stock of Argentine hake (Merluccius hubbsi). Fish. Res. 67: 25-38. https://doi.org/10.1016/j.fishres.2003.08.006 ), and was also reflected in the small size of the larvae. However, the RD values ​​obtained for the nutritional condition of this species were similar to those recorded in previous studies, as were the percentages of hake larvae in starvation and in optimal growth conditions (Diaz et al. 2014Diaz M.V., Olivar M.P., Macchi G.J. 2014. Larval condition of Merluccius hubbsi (Marini, 1933) in the northern Patagonian spawning ground. Fish. Res. 160: 60-68. https://doi.org/10.1016/j.fishres.2013.11.009 , Cohen et al. 2021Cohen S., Díaz A.O., Diaz M.V. 2021. Morphological and biochemical approaches to assess the nutritional condition of the Argentine hake Merluccius hubbsi larvae from two different nursery areas. J. Fish Biol. 98: 132-141. https://doi.org/10.1111/jfb.14563 ).

Various studies have inferred the importance of the NPFS as a nursery area for the southern anchovy stock (Bakun and Parrish 1991Bakun A., Parrish R.H. 1991. Comparative studies of coastal pelagic fish reproductive habitats: the anchovy (Engraulis anchoita) of the Southwestern Atlantic. ICES J. Mar. Sci. 48: 343-361. https://doi.org/10.1093/icesjms/48.3.343 , Diaz et al. 2016Diaz M.V., Do Souto M., Peralta M., et al. 2016. Comer o ser comido: factores que determinan la condición nutricional de larvas de Engraulis anchoita de la población patagónica de la especie. Ecología Austral 26:120-133. https://doi.org/10.25260/EA.16.26.2.0.71 , Do Souto et al. 2018Do Souto M., Spinelli M., Brown D.R., et al. 2018. Benefits of frontal waters for the growth of Engraulis anchoita larvae: The influence of food availability. Fish. Res. 204: 181-188. https://doi.org/10.1016/j.fishres.2018.02.019 ), with better feeding conditions for anchovy and hake larvae in the area associated with the NPFS (Viñas and Ramírez 1996Viñas M.D., Ramírez F.C. 1996. Gut analysis of first-feeding anchovy larvae from Patagonian spawning area in relation to food availability. Arch. Fish. Mar. Res. 43: 231-256., Temperoni et al. 2014Temperoni B., Viñas M.D., Martos P., Marrari M. 2014. Spatial patterns of copepod biodiversity in relation to a tidal front system in the main spawning and nursery area of the Argentine hake Merluccius hubbsi. J. Mar. Syst. 139: 433-445. https://doi.org/10.1016/j.jmarsys.2014.08.015 ). Our results support this idea of better feeding conditions towards the NPFS.

Within the analysed area, the internal zone to the position of the tidal front seems to respond to the “Bakun triad” hypothesis (Bakun and Parrish 1991Bakun A., Parrish R.H. 1991. Comparative studies of coastal pelagic fish reproductive habitats: the anchovy (Engraulis anchoita) of the Southwestern Atlantic. ICES J. Mar. Sci. 48: 343-361. https://doi.org/10.1093/icesjms/48.3.343 ) in that the frontal structure would guarantee the stability of the water column as a result of vertical stratification, nutrient enrichment and retention of spawning products within a favourable habitat. It is evident that frontal systems play a fundamental role in the ecological processes of the ocean because they allow a high primary production, they are suitable areas for the reproduction and feeding of many nektonic species, they offer a suitable breeding environment for the feeding of the early stages of fish development, and they act as retention zones (Acha et al. 2004Acha E.M., Mianzan H.W., Guerrero R.A., et al. 2004. Marine fronts at the continental shelves of austral South America: Physical and ecological processes. J. Mar. Syst. 44: 83-105. https://doi.org/10.1016/j.jmarsys.2003.09.005 ). NPFS could also be an area with a high concentration of prey, predators and competitors of anchovy and hake eggs and larvae (Mianzan and Guerrero 2000Mianzan H.W., Guerrero R.A. 2000. Environmental patterns and biomass distribution of gelatinous macrozooplankton. Three study cases in the Southwestern Atlantic Ocean. Sci. Mar. 64: 215-224. https://doi.org/10.3989/scimar.2000.64s1215 , Álvarez-Colombo et al. 2003Álvarez-Colombo G., Mianzan H., Madirolas A. 2003. Acoustic characterization of gelatinous-plankton aggregations: four case studies from the Argentine Continental shelf. ICES J. Mar. Sci. 60: 650-657. https://doi.org/10.1016/S1054-3139(03)00051-1 , Diaz et al. 2016Diaz M.V., Do Souto M., Peralta M., et al. 2016. Comer o ser comido: factores que determinan la condición nutricional de larvas de Engraulis anchoita de la población patagónica de la especie. Ecología Austral 26:120-133. https://doi.org/10.25260/EA.16.26.2.0.71 , 2020Diaz M.V., Do Souto M., Betti P., et al. 2020. Evaluating the role of endogenous and exogenous features on larval hake nutritional condition. Fish Oceanogr. 29: 584-596. https://doi.org/10.1111/fog.12497 ).

Regarding daily variations in nucleic acid content, the available information is not consistent: some authors have found differences in the RD ratio throughout the day, and others have observed no detectable pattern (Buckley et al. 1999Buckley B.A., Caldarone E.M., Ong T.L. 1999. RNA-DNA ratio and other nucleic acid-based indicators for growth and condition of marine fishes. Hydrobiologia 401: 265-277. https://doi.org/10.1023/A:1003798613241 ). For example, Rooker and Holt (1996)Rooker J.K., Holt G.J. 1996. Application of RNA:DNA ratios to evaluate the condition and growth of larval and juvenile red drum (Sciaenops ocellatus). Mar. Freshw. Res. 47: 283-290. https://doi.org/10.1071/MF9960283 while studying croaker larvae, found a marked daily pattern in the RD values. These authors observed that the RD values were higher during the day and lower at night and suggested that this phenomenon was due to differences in metabolic rates throughout the day, food requirements or food digestion. In this study, we observed an increase during the afternoon hours, but no statistically significant differences were observed during the morning and night. Previous studies have reported that E. anchoita larvae feed during daylight hours (Viñas and Ramírez 1996Viñas M.D., Ramírez F.C. 1996. Gut analysis of first-feeding anchovy larvae from Patagonian spawning area in relation to food availability. Arch. Fish. Mar. Res. 43: 231-256.), which could explain the increase in the RD values during the afternoon hours.

There is also conflicting evidence in the literature about the variability observed in the nutritional condition of larvae in the water column. Grønkjær et al. (1997)Grønkjær P., Clemmesen C.M., St. John M. 1997. Nutritional condition and vertical distribution of Baltic cod larvae. J. Fish Biol. 51: 352-369. https://doi.org/10.1111/j.1095-8649.1997.tb06108.x found that protein growth rates of Gadus morhua larvae were significantly higher for all age groups of larvae in the upper layers. In contrast, Dänhardt et al. (2007)Dänhardt A., Peck M.A., Clemmesen C.M., Temming, A. 2007. Depth-dependent nutritional condition of sprat Sprattus sprattus larvae in the central Bornholm Basin, Baltic Sea. Mar. Ecol. Prog. Ser. 341: 217-228. https://doi.org/10.3354/meps341217 , studying different indicators of nutritional condition in Sprattus sprattus larvae, did not find a better nutritional condition in the superficial layers, although this was not the case in all the condition indicators used. Palomera (1991)Palomera I. 1991. Vertical distribution of eggs and larvae of Engraulis encrasicolus in stratified waters of the western Mediterranean. Mar. Biol. 111: 37-44. https://doi.org/10.1007/BF01986343 observed that E. encrasicolus larvae mainly occurred above the level of the thermocline, and the highest abundances were recorded in the first 10 m of the water column. The thermocline could function as an upper or lower barrier for larval distribution, favouring the permanence in superficial layers with greater availability of prey associated with the depth of the thermocline (Smith and Suthers 1999Smith K.A, Suthers I.M. 1999. Displacement of diverse ichthyoplankton assemblages by a coastal upwelling event on the Sydney shelf. Mar. Ecol. Prog. Ser. 176: 49-62. https://doi.org/10.3354/meps176049 ). In our work, we found differences in anchovy larvae condition at different levels with respect to the vertical stratification of the water column. This fact should be considered in the future in designing the collection of samples and making comparisons of larvae condition. However, it is not easy to make comparisons between different studies because these aspects are very dependent on variables such as the zone, period, prey availability and species.

Although the anchovy is currently a low-exploitation species, it has a key position in the food web, regulating the systems towards lower and higher levels (Do Souto et al. 2018Do Souto M., Spinelli M., Brown D.R., et al. 2018. Benefits of frontal waters for the growth of Engraulis anchoita larvae: The influence of food availability. Fish. Res. 204: 181-188. https://doi.org/10.1016/j.fishres.2018.02.019 ). Changes in the annual abundances and mean lengths of this species have recently been recorded, so studying its life traits during its early ontogeny is important for understanding the variability of its recruitment (Orlando et al. 2019Orlando P., Buratti C., Garciarena A.D. 2019. Diagnóstico de la población de anchoita bonaerense (Engraulis anchoita) y estimación de captura biológicamente aceptable durante el año 2019. Inf. Téc. INIDEP Nº 24, Mar del Plata, 29 pp.). Currently, with a continuous increase in fishing effort, many of the main resources are exploited to the limit of their possibilities. The impact of fishing on ecosystems leads to a decline in commercial and non-commercial species. There is a clear need for an integrated vision for the correct management of fishing resources, including both commercial and non-commercial species in the analysis (Coll and Palomera 2007Coll M., Palomera I. 2007. Hacia el estudio y la gestión pesquera basada en los ecosistemas. Ecología política 21: 87-89., Pauly 2009Pauly D. 2009. Beyond duplicity and ignorance in global fisheries. Sci. Mar. 73: 215-224. https://doi.org/10.3989/scimar.2009.73n2215 ). Therefore, the measures adopted must consider integral care of the ecosystem, including zonation of the oceans and the generation of new MPAs (Pauly 2009Pauly D. 2009. Beyond duplicity and ignorance in global fisheries. Sci. Mar. 73: 215-224. https://doi.org/10.3989/scimar.2009.73n2215 ).

Previous studies have shown that the application of MPAs has positively influenced the nutritional condition of fish species of commercial interest (Viladrich et al. 2016Viladrich N., Rossi S., Lopez-Sanz A., Orejas C. 2016. Nutritional condition of two coastal rocky fishes and the potential role of a marine protected area. Mar. Ecol. 37: 46-63. https://doi.org/10.1111/maec.12247 ). The RD index and its derivatives are highly sensitive, which has made it possible to map the larval condition and establish favourable areas in relation to the NPFS. Studying nutritional condition allows us to estimate the probabilities of survival of the organisms and their potential recruitment to fisheries. The results presented here are relevant as a baseline for future studies that consider the evolution of the larval condition of these two species and evaluate the potential of areas as MPAs. The RD index is a simple and useful tool that provides complementary information to that provided by indicators of diversity and abundance usually used in strategies aimed at ensuring the conservation and comprehensive management of these fishery resources in the Argentine Sea.

In conclusion, the results presented herein are of great importance because little information is available on the nutritional condition of E. anchoita larvae from the southern stock. In this species, lower RD values were observed​​ than those previously recorded in the study area. Nevertheless, as was hypothesized, the larvae collected in the areas influenced by the NPFS showed a good nutritional condition. The area inside the NPFS showed a lower proportion of individuals in starvation and a higher proportion of individuals in optimal growth, whereas the area outside the NPFS showed a high percentage of anchovy larvae below the critical value of the RD index. It was also observed that the hake RD values and the percentages of larvae in starvation at the beginning of the spawning period in the austral summer of 2018 were similar to those recorded in previous studies. Finally, the RD ratio and its derived index are highly sensitive, making it possible to “map” the larval condition and establish the most favourable areas for survival during the early phases of the life cycle.