Seasonal study of the nutritional composition of unexploited and low commercial value fish species from the Portuguese coast

Abstract Target species diversification is essential for fisheries sustainability and fish market revitalization. Fish discards are a widely recognized problem resulting from fisheries worldwide, and are of major concern for all sector players, from administrations, to fishermen, and scientists. However, non‐target species are seldom studied, and information on nutritional profiles and seasonal changes in nutritional properties is generally lacking. This study assessed the seasonal nutritional composition of two unexploited (Serranus cabrilla, Capros aper) and three low commercial value fish species (Trachurus picturatus, Spondyliosoma cantharus, and Trigla lyra), captured on the Portuguese coast over 1 year. Significant seasonal variations were observed in the nutritional composition of all the species studied. Moisture and ash contents varied from 70% to 81% and from 5% to 13%, respectively. The maximum fat contents were 5% for C. aper and 4% for T. picturatus, allowing to classify all studied fishes as lean. The highest protein contents were recorded for C. aper (25%) and S. cantharus (20%). The unexploited and low commercial value fish species studied were shown to be good fat and protein sources, comparable to commonly consumed species, such as cod and salmon, having a great potential to become commonly consumed fish in Portugal.


| INTRODUC TI ON
The growing world population is increasing the nutritional demand, causing the overexploitation of food sources such as fish (FAO, 2020). Thus, discarded fish species could become new alternatives to the consumer's diet, bringing benefits to fishermen as well as to companies with the advantage of limiting the environmental impact. Fish discards may result in the rejection of nutritionally rich seafood, which may have similar contents of bioactive compounds, high-quality proteins, fats, and minerals, to those found in commonly consumed seafood (Blanco et al., 2018). Fish discards include low or no commercial value species, bycatches of fish under minimum conservation reference sizes, nonquota or exhausted quota fish species, and crushed or damaged fish (FAO, 1996). According to the last FAO's update, between 2010 and 2014, the annual discards from global marine capture fisheries were 9.1 million tons, of which 46% were from bottom trawls (Pérez Roda et al., 2019). Discarded fish have generally low survival rates leading to consequences in marine ecosystems (Guillen et al., 2018).
Therefore, the valorization of discarded fish is important not only for the fishery sector but also for marine systems sustainability. In order to avoid and/or reduce fish discards, a landing obligation regulation was implemented in 2013 for the European Fisheries (Regulation [EU] 1380[EU] , 2013, which includes all the species under TAC (total allowable catch) regulations.
Low or no commercial value species are commonly used for animal feeding purposes or as raw material for aquaculture feeding.
However, 12% of the global fish production fall in this category, demanding alternative solutions for such species (ICES, 2006). One of the approaches to promote the consumption of low or no commercial value species involves the assessment of their nutritional content, given the growing relevance of this topic to consumers and, consequently, to the industry. The main components of fish muscle are water, protein, and lipids, being important constituents when determining its nutritional value. Fish protein has high nutritional value, because it includes all the essential amino acids (e.g., methionine and lysine) and is also highly digestible when compared with terrestrial meat proteins (Tilami & Sampels, 2017). Fish is also an important source of lipids in the human diet, being rich in polyunsaturated fatty acids, such as omega-3, which contribute to the normal functioning of the cardiovascular, neural, and immunity systems (Tilami & Sampels, 2017). Moisture content is an important measure in food products because it is a quality and safety factor in food preservation (Nielsen, 2010).
The fish species that are more consumed in Europe are tuna (Thunnus sp.), Atlantic cod (Gadus morhua, Linnaeus, 1758), and Atlantic salmon (Salmo salar, Linnaeus, 1758) (European Commission, 2019). The most consumed fish species in Portugal are tuna, cod, and hake (Merluccius merluccius, Linnaeus, 1758) ( Barroso et al., 2020). However, the most caught fish species on the Portuguese coast are sardines, mackerel, tuna, and horse mackerel (in decreasing order) (INE, 2011(INE, -2021. This consumption preference is related to a lack of consumer awareness of other fish species, or due to the apprehensiveness to purchase and cook unfamiliar fish species. Thus, the valorization of unexplored species would be an advantage for the country's economy and sustainability. In this context, five fish species with low or no commercial value captured in the Portuguese coast were selected in order to study their nutritional value ( Figure 1). Among the species of low commercial value, blue jack mackerel (Trachurus picturatus, Bowdich, 1825), black seabream (Spondyliosoma cantharus, Linnaeus, 1758), and piper gurnard (Trigla lyra, Linnaeus, 1758) are particularly important for their landings and first auction prices (INE, 2011(INE, -2021. Regarding the species without commercial value, comber (Serranus cabrilla, Linnaeus, 1758) and boarfish (Capros aper, Linnaeus, 1758) are particularly abundant in the Portuguese coast (INE, 2011(INE, -2021. The rejection of these fish species by the consumer is mainly a cultural issue, as they are unknown by the majority of the population, despite being savory and nutritious. There are also other factors adding to the low popularity of some of these species, e.g., T. picturatus is often undervalued due to its darkened muscle coloration in comparison with Trachurus trachurus, and C. aper is a small fish.
Comber and boarfish are among the most important species in terms of abundance and discards, being mainly caught as by-catch.
S. cabrilla belongs to the Serranidae family and can be found in sandy and rocky bottoms at depths up to 500 m in eastern Atlantic from the Bay of Biscay to South Africa and the Mediterranean (Gordo et al., 2016). This hermaphroditic species is mainly caught as a by-catch in trawls and trammel nets and is discarded at sea or grouped in the "diverse" category at Portuguese auctions (Gordo et al., 2016;Sequeira et al., 2015). C. aper belongs to the Caproidae family and is the only species in its genus (Egerton et al., 2017). It is found from Norway to Senegal and also in superficial waters of the Mediterranean Sea, but it is more abundant on rocky bottoms at the edge of the continental slopes (Sequeira et al., 2015). In general, this species is discarded at sea whenever caught. T. picturatus belongs to the Carangidae family, is benthopelagic, and is found in the East Atlantic seas (including Azores Islands, Madeira, and the Canary Archipelagos). Landings on the mainland Portuguese coast, although highly variable, reached an average value of ca 3900 tons/ year in the last decade (INE, 2011(INE, -2021. In Madeira Archipelago, T. picturatus represents about 5% of the total landings, being part of the diet of local population and often used as bait for tuna fish catch (Vasconcelos et al., 2006). The Sparidae family, which includes S. cantharus, is the most targeted in Portuguese coastal fisheries, being caught by several different fishing gears, with mean annual landings of 220 tons in the last decade (INE, 2011(INE, -2021. S. cantharus is distributed in the Mediterranean Sea and eastern Atlantic Ocean (including the archipelagos of Madeira, Cape Verde, and Canaries), in rocky and sand substrates up to 300 m deep (Rodríguez et al., 2004;Neves et al., 2018). T. lyra belongs to the Triglidae family and is a demersal fish found in a wide variety of locations, including archipelagos of Azores and Madeira (Papaconstantinou, 1981), with mean annual landings of 370 tons in the last decade (INE, 2011(INE, -2021. The aim of this work was to show that the selected unexploited and low commercial value fish species have a great nutritional potential to be more regularly consumed. Hence, the nutritional composition of the five fish species was analyzed. The study was carried out throughout 1 year since the nutritional parameters may vary along seasons. In previous studies, the selected fish species were shown to be sensorially appealing , and good sources to be used in the development of marine-based food products with high acceptance from consumers (Silva et al., 2021).

| Sample collection and pretreatment
Fish samples, captured from the coast of Portugal (approximate capture zone between 39°N, 9.5°W and 40°N, 11°W), were acquired at Peniche fish market (center of Portugal) or purchased from local fishermen, every 2 months (Jan/Feb, Mar/Apr, May/Jun, Jul/Aug, Sep/ Oct, Nov/Dec), during 2019. Sample preparation was carried out according to Portuguese Standard NP 4424 (2004). Briefly, fishes were washed with deionized type III water, measured for the total length (from mouth to caudal fin, cm), gutted, and filleted, manually rejecting the skin. Then, fish were divided into fresh samples and freezedried samples. To obtain freeze-dried samples, fish were previously weighted and kept at −80°C for 12 h and then freeze-dried at −56°C for 60 h (Freeze-drier Labogene, CoolSafe 55-4). The freeze-dried biomass was then crushed using a kitchen robot (Bimby, Vorwerk, Thermomixer 31-1) and the resulting powder was stored in sampling bags in a fresh, dark, and dry room. Fresh samples were used to determine moisture and ash contents, and dried samples were homogenized to assess total fat and crude protein content. The number of individuals of each species used for the analyses was defined according to their length (NP 4424, 2004). The sample included 10 specimens for individuals up to 15 cm, 5 specimens for individuals up to 30 cm, and for those individuals over 30 cm, 3 cross-sections with 2.5 cm of each specimen were cut in the area after the pectoral fin, just before the caudal fin and the middle distance between both.
Thus, the study included 10 individuals for S. cantharus, T. lyra, and T. picturatus, 20 for C. aper, and 15 for S. cabrilla, for each period of analysis (for each 2 months).

| Moisture content
Moisture content assay was performed according to Portuguese Standard NP 282 (1991). Porcelain crucibles were incubated for 3 h at 105 °C, to remove moisture and other substances that could affect the effective crucible weight. After cooling to room temperature in a desiccator, 10 g of fresh sample were weighed into the pre-weighed crucibles (Analytical Balance Sartorius, ENTIRIS224I-1S). The crucibles with the samples were placed into a drying oven (Memmert, UF110) at the mentioned temperature and time. The samples were then removed and cooled down to room temperature in a desiccator (ca. 1 h). The samples were kept inside the desiccator and weighed daily until the values from successive weighings did not differ by more than 10 mg (3-4 days). Moisture content was expressed as g/100 g of fresh weight (FW) and calculated using the following equation: where W fresh is the fresh sample weigh and W dehydrated is the average weigh of the dehydrated sample (average of the three last weighings).

| Ash content
The ash content assay was performed according to Portuguese Standard NP 2032(1988. The crucibles incubation and cooling procedures were the same as described for the moisture assay. For this assay, 5 g of fresh sample were weighed, and the samples were placed in an incinerator (Nabertherm, B170) at 525°C for 4 h. The sample ashes were cooled to room temperature in a desiccator and their weight was recorded until the values from successive weighings did not differ by more than 1 mg. The ash content was expressed as g/100 g of FW and determined using the following equation: where W ash is the average ash weight and W fresh is the fresh sample weight.

| Total fat content
The total fat determination was adapted from Folch et al. (1957).
Briefly, 1 g of dried sample was mixed with 0.8 ml of distilled water and 5 ml of Folch reagent and mixed in a vortex (VWR, VV3) for 1 min. Folch reagent was prepared using CHCl 3 and MeOH (VWR, Fontenay-sous-Bois, France) in a proportion of 2:1 (v:v). Further 5 ml of Folch reagent were added and mixed for 5 min, and then for more 2 min after adding 1.2 ml of 0.8% NaCl solution (Biochem, Cosne sur loire, France). The samples were centrifuged (4347 × g, 10 min, 4°C; Centrifuge Eppendorf, 5810R) and the lower phase (organic) was filtered through hydrophilic cotton and anhydrous sodium sulfate (VWR, Leuven, Belgium) to round-bottomed glass flask. To avoid sample loss, more 5 ml of CHCl 3 were added, followed by a second centrifugation and filtration, under the same conditions. The solvent from the collected organic phases was removed through lowpressure evaporation using a rotary evaporator (Heidolph, Laborota 4000). After complete solvent evaporation, the samples were placed in an oven at 30°C for 4 h and then kept in a desiccator until a constant weight was registered (less than 10 mg between successive weighings). The total fat content was expressed as g/100 g of FW and determined using the following equation: where W fat is the average weight of the fat (weight of roundbottomed glass flask with the fat -weight of the empty roundbottomed glass flask), and W sample is the sample's initial weight. The crude protein content, represented by the sample's nitrogen content, is expressed as g/100 g of FW, and was calculated using the following equation: where 6.25 is the nitrogen-to-protein conversion factor (Jones, 1941), V a is the volume of HCl spent on sample titration, V b is the volume of HCl spent on blank assay, N is the HCl normality, and W is the sample weight.

| Average percentage variation over the year
The variation in the percentage of moisture, fat, and protein over the year was obtained as follows. The variation module between pairs of months was calculated according to the following equation: where M i is the first pair of months, M i+1 is the next couple of months, and max is the highest value between M i and M i+1 . Then, the annual variation was obtained with the average of the variations between pairs of months.

| Statistical analysis
One-way analysis of variance (ANOVA) was performed to identify significant differences in fish length and nutritional composition throughout the year (pairs of months/seasons) for each species.
All assumptions required to perform the ANOVA (that is, data normality and homogeneity of variance) were validated (by applying the Shapiro-Wilk and Levene tests, respectively), and the Tukey test was used in post hoc multi-comparisons (F value and p-value) (Zar, 2010). When assumptions were not met, the nonparametric Kruskal-Wallis test and the Games-Howell multi-comparison test were applied (Zar, 2010 (Bro and Smilde, 2014). The data were mean-centered by subtracting from every variable the corresponding mean level, and the variables were scaled by dividing them by the corresponding standard deviations. Thus, PCA proved to be adequate, since it allows to reduce the data dimensionality, maintaining relevant information contained therein (Bednárová, 2013). Although only the results concerning the first two components are present, the others were also analyzed. Analyses were performed using Canoco for Windows 4.5 software (Ter Braak and Smilauer, 2002).
All analyses were performed in triplicate. Whenever applicable, the results are presented by the mean ± standard deviation (SD).

| Characterization of the samples
Biotic, abiotic factors or even fishing pressure can influence the growth and consequently the fish length. Among the biotic factors, food availability and low competition for food are the most relevant.
Fishing pressure can lead to populations of smaller individuals since the larger ones are more susceptible to be caught by the different fishing gears (King, 2013). According to Webster et al. (2011), intraspecific morphological variations can be associated with environmental differences such as flow regime, water depth, water chemistry, and substrate type. Also, the water temperature can influence growth rate, body length, and reproductive schedules, since reduced temperatures result in lower initial growth rate and delayed maturation in larger body length (Trip, 2014).
The total length (cm) and number of specimens used for each species over the year are presented in Table 1. C. aper was the smallest species studied, with an annual average length of 12.7 ± 0.6 cm, while T. lyra and T. picturatus were the longer species, with annual average lengths of about 30 cm. T. lyra was also the species that presented the highest annual variation in total length. Trachurus picturatus also presented statistically significant differences in the length throughout the year, with the minimum length (25.6 ± 1.5 cm) in Jan/Feb and the maximum length (40.6 ± 1.2 cm) in May/Jun (Table 1). These values agree with the data reported for T. picturatus from Madeira and Azores archipelagos, with lengths ranging from 17 to 46 cm (Vasconcelos et al., 2006) and from 15 to 40 cm (Isidro, 1990), respectively.
The average length of C. aper through the year also presented statistically significant differences, presenting a constant length throughout the year (about 12 cm) and a maximum value in Sep/Oct ( 11.9 ± 0.9 a 18.9 ± 2.1 b 18.9 ± 4.0 BC 17 (4-34) Note: Values are presented as mean ± SD. Different lowercase letters in rows represent significant differences between pairs of months for each species and different uppercase letters in columns represent significant differences between species (p-value <.05). n is the number of specimens used for each pair of months. Annual variation (%) is presented as Av. (Min.-Max.), were Av. is the average of the variations between pairs of months, and Min and Max are the minimum and maximum variations observed between pairs of months, respectively.
[#] Data missing due to mechanical problems in the supply vessel.
this species captured off the Portuguese coast, with a total length ranging between 7.2 and 9.5 cm for immature specimens and 8.1 and 17.0 cm for mature ones . A maximum length of 23 cm was reported for C. aper from the Northeast Atlantic, but the length range for this species varies between 2.5 and 18 cm with 12 cm being a more commonly observed length .
Length at maturity was estimated to be at 9.7 cm in this region (Farrel et al., 2012), which is similar to the maximum length of an immature fish sampled in the Portuguese coast (9.5 cm).
As for S. cabrilla the length also presented statistically significant differences throughout the year, with the lowest length value (11.9 ± 0.9 cm) in Sep/Oct and a constant length of about 20 cm in the other periods (Table 1). These results are similar to those reported for this species captured off the Portuguese coast between July 2011 and September 2012, with a length variation between 12 and 26 cm (Gordo et al., 2016). 77.17 ± 0.38 g/100 g fresh wet mass and higher than female individuals (74.31 ± 1.78 g/100 g fresh wet mass).

| Ash content
The ash content refers to the inorganic residue remaining after incineration of fresh sample, representing the total mineral content in foods (Nielsen, 2010). Ash content by species and months is presented in Table 2. The variation in ash content throughout the year was statistically significant for S. cantharus, T. lyra, T. picturatus, and C. aper. The highest average variation in ash content throughout the year was observed for T. lyra (16%), while for the other species the variation was in average 8-11%. Ash content ranged between 5% and 13% for the species and period studied, which is higher than the ash content of 2.5% reported for raw fish on a wet weight basis (Nielsen, 2010). However, the values obtained are in agreement with what is reported in the literature, including ash content ranges from 0.7 to 5.3% (Bogard et al., 2015) and from 1.9 to 17.8% (Flowra et al., 2012). Note: Values are presented as mean ± SD, and expressed as g/100 g of fresh weight. Different lowercase letters represent significant differences between pairs of months for each species and different uppercase letters in columns represent significant differences between species for each parameter (p-value <.05). Annual variation (%) is the average of the variation between pairs of months.
[#] Data missing due to mechanical problems in the supply vessel.
The average ash content obtained for C. aper in the present study is higher than the value reported by Spitz et al. (2010) with 4.6% of wet total body mass, captured in the Bay of Biscay in autumn. The average ash content obtained for S. cabrilla was similar to that reported by Edirisinghe et al. (2013) for E. merra male and female specimens, captured in Sri Lanka, with 5.67 ± 0.88 and 5.36 ± 0.56 g/100 g dry mass, respectively.
The obtained range of ash content indicates that the species studied are a good source of minerals. However, the average ash content obtained for all the species is significantly higher than most of those reported in the literature. A possible explanation for this may have been the accidental inclusion of bones as edible parts in the manual filleting, which would have resulted on higher ash content in the samples.

| Fat content
Variations in the fat content among pairs of months were statistically significant for C. aper, S. cantharus, T. lyra, and T. picturatus, but not for S. cabrilla (Table 2) Aug. This could be related to an increase in food availability in the west coast of Portugal where strongest summer upwelling occurs (July to September) (Fiúza et al., 1982). Spitz et al. (2010) reported a fat content of 4.8% of wet total body mass for C. aper captured in the Bay of Biscay in autumn, which is higher than the value obtained for Nov/Dec in the present study (1.52%, Table 2), but similar to the maximum fat content of 5% obtained in May/Jun in this study. These variations may be due to the different environmental conditions and food availability in the two locations along the year. S. cabrilla had the lowest fat content among the fish species studied, with an average fat content of 0.68 ± 0.09 g/100 g and no statistical differences throughout the year (Table 2).
According to their fat content, fish may be classified as lean fish (less than 5% fat), such as cod, mid-fat fish (5-10% fat), such as bream, and fatty fish (10-25% fat), such as some salmon species, although the classification and the lipid ranges in fish may vary among publications (Taşbozan & Gökçe, 2017). The species studied herein presented average fat contents ranging between 0.7% and 3.7% and may be considered lean fishes according to the referred classification. Additionally, the fat contents obtained for the species in the present study are comparable to those reported for commonly consumed fish such as G. morhua, with a fat content of 0.1-0.9% (Murray & Burt, 1991), and M. merluccius, with a fat content of 0.27-1.24% (Küçükgülmez et al., 2008), and inferior to fatty fish, such as Thunnus sp., with 25% of fat (Murray & Burt, 1991 (Blanco et al., 2018). The authors reported a protein content of 17.270 ± 0.140%, which is similar to the average value obtained in the present study throughout the year (18%, Table 2). T. picturatus revealed a minimum protein content in Jan/Feb, and a maximum value in Jul/Aug. Similar values were reported for T. murphyi from Chile, with a protein content of 16.58 ± 0.88 g/100 g (Bastías et al., 2017), but a much higher protein content was reported for T. trachurus captured in Nigeria (65.65 ± 1.01%) (Adeyemi et al., 2013). The minimum protein content of C. aper was observed in Jan/Feb, reaching its maximum value in Jul/Aug. Lower values were reported for C. aper captured off the northwestern coast of Ireland with an average protein content of 6.1 ± 1.0% (Egerton et al., 2020). However, comparable protein contents (17.2%) were reported for C. aper captured in the Bay of Biscay in autumn (Spitz et al., 2010). The wide variations in protein content observed throughout the year for this species, with a maximum value in Jul/Aug, could be related to an increase in food availability in the western coast of Portugal in that period (Fiúza et al., 1982;Mittelstaedt, 1986), similar to what was observed for the variation of its fat content along the year. S. cabrilla protein content was minimum in Jan/Feb, reaching its maximum value in May/ Jun, which was not statistically different from the value obtained in Nov/Dec.
Fish protein content is usually between 15 and 20%, although some species may present protein contents out of this range (Murray & Burt, 1991). The protein contents obtained for the species in the present study (14-25%) are within the aforementioned protein contents. Moreover, the protein contents obtained are comparable to those reported for commonly consumed fish, such as G. morhua with 15% of protein (Murray & Burt, 1991), S. salar with 19-21% (Costa et al., 2015), M. merluccius with 18-21% (Küçükgülmez et al., 2008), and Thunnus sp. with 25% of protein (Murray & Burt, 1991).
Therefore, these unexploited and low commercial value fish species are a good protein alternative to the commonly consumed fish species.

| Seasonal variations in the nutritional composition
A PCA was conducted to determine possible variations in the nutritional composition of the sudied species (Figure 2). Since the analyses were performed every two months, the following association of months with seasons was considered: winter -Jan/Feb and Mar/  (Garcia et al., 2015;Vasconcelos et al., 2017). The minimum protein and fat levels are observed in winter, during the spawning season. In this period, both fat and protein are used to fulfill the higher energy demands. During breeding, the lipid fraction is the nutritional component with the widest variation, reaching its minimum value around that period. This occurs because the spawning requires higher energy intake, so the energy reserves are spent and most of the species do not ingest much food during spawning migration, enabling the intake of energy through feeding (Huss, 1995).
Higher levels of fat are also associated with C. aper in spring (May/Jun) and higher protein contents are associated with C. aper in summer (Jul/Aug). This species spawns all year round off the Portuguese coast, with a spawning peak from June to August (Sequeira et al., 2015). The spawning period for C. aper does not correspond to the period of the year when the lower fat and protein contents are observed, contrasting with what was observed for T. picturatus. The seasonal variations in protein and fat contents for this species could be related to an increase in food availability in the west coast of Portugal in that period (Fiúza et al., 1982;Mittelstaedt, 1986).
Regarding T. lyra, higher protein and fat levels are observed in autumn (Nov/Dec). This species has a short spawning season, from November to February with a peak in January (Sousa, 2019). Thus, the maximum fat and protein contents coincided with the beginning of the spawning season, while minimum levels of these components were observed in the summer.
As for S. cantharus, higher moisture levels were observed in the spring (May/Jun). In general, it was observed that muscle moisture content varied inversely with the amount of fat (Ali et al., 2020).
The increase in fat content and the simultaneous decrease in water content may be related with an increase in age (increasing size). The fluctuations in water, fat, protein, and ash contents with the age may also be related to migration and spawning. This was clearly observed for S. cantharus that presents maximum moisture levels in spring when the fat and protein contents are minimum.
This corresponded to the end of the spawning season of this species, which usually occurs between February and May (Neves et al., 2018).
Regarding S. cabrilla, this was the species that presented the smallest annual variation in the nutritional parameters studied (Table 2). As shown in Figure 2, no strong associations could be observed between any of the nutritional parameters and the months/seasons analyzed. Despite of this, the maximum protein and fat contents for this species were obtained in spring, which is coincident with the beginning of the spawning season of this species in the Portuguese coast that occurs from May to July (Sequeira et al., 2015).

| CON CLUS ION
This study assessed the seasonal nutritional composition of unexploited or low commercial value fish species captured on the Portuguese coast over one year. The species studied were mainly composed by water (69-82%), protein (14-25%), ash (5-13%), and fat (0.6-5%). Significant seasonal variations were observed in the nutritional composition of all the species studied. Fat was the parameter that presented the widest variations throughout the year, with notable differences in C. aper from different months. All the studied fish species were considered lean, according to their fat content. C. aper and S. cantharus were found to have high protein contents. The fish species studied were shown to have high fat and protein contents, comparable to those of commonly consumed species such as G. morhua and M. merluccius. Therefore, these species demonstrate a great potential to become more commonly consumed fish species, either for human or animal consumption.
The observed seasonal differences appeared to be related to the reproductive cycle of the species, but many other environmental factors, such as food availability, temperature, age, size, and sex, may also contribute to these variations. C. aper presented higher fat and protein contents in July/August, while for the other species these parameters were higher from September to December, except for