Can weight/length relationship predict size at first maturity? A case study with two species of Characidae

Can weight/length relationship predict size at first maturity? A case study with two species of Characidae

Nelson F. Fontoura; Anamélia S. Jesus; Gabriel G. Larre; Juliana R. Porto

Faculdade de Biociências, Pontifícia Universidade Católica do Rio Grande do Sul - PUCRS. Av. Ipiranga, 6681, P. O. Box 1429, 90619-900 Porto Alegre, RS, Brazil. [e-mail]

ABSTRACT

A weight/length relationship was established for Astyanax jacuhiensis (Cope, 1894) (n = 370) and Cheirodon ibicuhiensis Eigenmann, 1915 (n = 701), from samples taken monthly in Fortaleza Lagoon, Cidreira, Rio Grande do Sul, from December 1991 through November 1992. Both species showed a polyphasic allometric growth pattern, each stanza described by an independent power equation controlled by a switch function. For C. ibicuhiensis, this change in the growth pattern occurred at 2.948 cm standard length (SL), close to published sizes for the attainment of female maturity. The change in the growth pattern of A. jacuhiensis (SL = 3.481 cm) was below the predicted size at first maturity, and merits further investigation. Although not conclusive, our data suggest that a complex growth pattern is frequent in nature, and perhaps is not usually identified because trends are obscured by natural variability. Despite the increased complexity resulting from the application of a more-complex equation, the identification of a change in the growth pattern could indicate important aspects of fish biology, including the attainment of sexual maturity.

Key words: Allometric growth, Polyphasic growth, Growth stanza, Reproduction, Fish.

RESUMO

A relação peso/comprimento para Astyanax jacuhiensis (Cope, 1894) (n = 370) e Cheirodon ibicuhiensis Eigenmann, 1915 (n = 701) foi estimada a partir de amostras mensais efetuadas na lagoa Fortaleza, Cidreira, Rio Grande do Sul (dezembro de 1991 a novembro de 1992). As espécies analisadas apresentaram padrão de crescimento alométrico polifásico, sendo cada fase descrita por uma equação potência independente controlada por uma função interruptora. Para C. ibicuhiensis foi identificada uma modificação no padrão de crescimento com 2,948 cm de comprimento padrão (CP), valor próximo à dados publicados relativos ao início da maturação sexual. A alteração observada no padrão de crescimento de A. jacuhiensis (CP = 3,481 cm) encontra-se abaixo do tamanho previsto para a primeira maturação e sugere a necessidade de investigação adicional. Embora não conclusivos, nossos dados sugerem que padrões complexos de crescimento são frequentes na natureza e, talvez, não identificados porque tendências sutis acabam por ser obscurecidas pela variabilidade natural dos dados. Apesar do aumento de complexidade resultante da aplicação de uma equação complexa, a identificação de mudanças no padrão de crescimento pode apontar para aspectos relevantes da biologia de peixes, incluindo a maturação sexual.

Introduction

The weight/length relationship is an important tool to analyze fish populations. Its applications range from simple estimates of an individual's weight to indication of fish body condition factor or inferences regarding sexual development (Le Cren, 1951). Knowledge on this relationship also helps to identify energy investments for growth or reproduction as a natural cyclic phenomenon of natural populations (Bolger & Connolly, 1989).

The weight/length relationship is adjusted through a power equation described by Huxley (1924): W = a.L^b, where W is the individual weight, L is the length, a is the proportionality coefficient or condition factor, and b is the allometric coefficient. As a simple model, the weight/length relationship is widely applied for descriptive purposes. Bervian et al. (2006), however, revised the Huxley relative growth function and proposed a polyphasic model in which the growth pattern can be divided into different stanzas connected by switch equations. Through analyses of the relative growth of otoliths in Micropogonias furnieri (Desmarest, 1893) (Sciaenidae), Bervian et al. (2006) identified a change in the growth pattern related to the attainment of the size at first maturity.

Because weight/length data are available for many fish species, extracting additional information from regular weight/length relationships may be useful. The present study aimed to identify the existence of polyphasic growth by means of analysis of the weight/length relationship of fish, and its possible linkages to the attainment of sexual maturity, following the general model proposed by Bervian et al. (2006). As a study case, data from two species of Characidae were applied. Cheirodon ibicuhiensis Eigenmann, 1915 is found in the Uruguay River basin, the Patos Lagoon basin, and the Tramandaí River system, and attains a maximum standard length (SL) of 4.2 cm (Lima et al., 2003). Reported minimum sizes of mature females range from 2.7 to 3.2 cm SL (Braun et al., 2000; Oliveira et al., 2002). Astyanax jacuhiensis (Cope, 1894) is found in the Jacuí and Uruguay River basins and the coastal drainages of southern Brazil, and can reach 11 cm SL (Lima et al., 2003). No published data concerning size at first maturity are available for this species. Using the general model proposed by Binohlan & Froese (2009), the size at first maturity for A. jacuhiensis is estimated as 6.8 cm SL, based on a maximum length of 11 cm SL.

Although few data are available for comparison, Table 1 presents data on maximum length and estimates of size at first maturity (L₅₀) for species of small characids.

Material and Methods

Monthly samples were taken from December 1991 through November 1992 in the Fortaleza Lagoon (Tramandaí basin, Coastal Drainage, southern Brazil, 30º08'S 50º13'W) by means of a beach seine net 15 m long, 3 m high, and with 5 mm mesh size. The hauls were made from a sandy beach and near the stands of macrophytes. All captured fish were preserved in 4% formalin. In the laboratory, fish were identified, weighed (W, to the nearest 0.001 g) and measured (SL, to the nearest 0.01 cm). Because gender information is not available (the data analyzed were obtained from the laboratory data bank), the data set was analyzed as a whole (merged gender). Voucher specimens for A. jacuhiensis and C. ibicuhiensis are respectively deposited in the Fish Collections of the Universidade Federal do Rio Grande do Sul (UFRGS 2003) and Museu de Ciências e Tecnologia da PUCRS (MCP 13893).

For each species, the first procedure was to adjust a regular weight/length relationship: W=a.L^b; where W is the expected weight (g) for a specific length, a is the proportionality coefficient, L is the standard length (cm), and b is the allometric coefficient. The objective was to identify possible trends (residual not evenly distributed along the y-axis zero baseline) by plotting proportional residuals as a function of the observed standard length. Proportional residuals ([observed weight predicted weight] / predicted weight) were used instead of regular residuals (observed weight - predicted weight) to minimize bias from heteroscedasticity. Also, to avoid bias from log-transformation (Smith, 1980), the standard Huxley allometric equations were adjusted by using the Solver routine of the Microsoft Excel (loss function using least squares of proportional residues, Newton search algorithm).

If the weight/length data follow a regular power function, proportional residuals should distribute randomly along the y-axis zero baseline. Any residual pattern different from zerocentered horizontal distribution indicates the model's inability to describe the growth pattern, and suggests the presence of a growth anomaly such as polyphasic growth. In this case, a break in the residual pattern is an indication of a point of change in the growth cycle.

Following Bervian et al. (2006), polyphasic growth can be divided into different stanzas, each stanza described by a regular power function. According to these authors, the Huxley power function is appropriate to describe relative growth, but due to environmental pressures or development strategy, a change in the growth pattern could be expected in the life history of animals. Therefore, a polyphasic growth function is no more then a series of regular power functions (Y = a.X^b), each one controlled by a switch function that turns a growth segment on or off: S_w =1/(1+e^(Tx*(L-SCL)); where S_w is the interruption factor (ranging from one to zero),T_x is the rate of change, L is the standard length, and SCL is the mean length at which individuals switch from one stanza to the next.

Considering a growth pattern with two stanzas, the weigh/length relationship can be described as follows: W = (a₁.L^b1).(1/(1+e^(Tx*(L-SCL))) + (a₂.L^b2).(1/(1+e^(-Tx*(L-SCL)); where W is the expected weight for a specific length L, a₁ and a₂ are the proportionality coefficients for the first and second stanzas, b₁ and b₂ are the allometric coefficients for the first and second stanzas, T_x is the rate of change from the first to the second stanza, and SCL is the Stanza Changing Length.

As a first approximation, each growth stanza was visually identified as a clump of residuals presenting any particular trend, and the Stanza Changing Length (SCL) established as the middle-point between two successive growth stanzas. Then a new weight/length relationship was individually adjusted for each stanza, below and above the Stanza Changing Length. Next, the parameters a and b for each growth stanza, as well as the size at the Stanza Changing Length (SCL) were used as seed values to adjust the function as a whole. In the absence of a known biological value, the rate of change (T_x) was set at 1.0 as the initial value. The function adjustment was made by using the Microsoft Excel Solver routine using Least Squares (proportional residues) as the loss function and the Newton search algorithm. Considering the large sample size, Standard Errors for each parameter were estimated by adjusting the polyphasic models 20 times, each time with only half of the data set, using random resampling without replacement, i. e., emulating a bootstrap procedure. Parameter Standard Errors were estimated as the Standard Deviations of all the 20 different computations. The Coefficient of Determination, R², was calculated for regular and polyphasic models as follows: R² = 1 - (residual Sum of Squares /total Sum of Squares) (Zar, 1999).

Considering that reproduction can change dramatically the energy budget directed for body growth, and could be reflected in the growth pattern, we suppose that the presence of a Stanza Changing Length (SCL) could be indicative of the size at first maturity (L₅₀). So, for comparison purpose, estimates of the Stanza Changing Length were compared to published values of L₅₀ for A. jacuhiensis and C. ibicuhiensis, and to L₅₀ estimated values according to the general model proposed by Binohlan & Froese (2009): L_{50 =} 10-^{0.1189+0.9157*Log(Lmax)}, where L₅₀ is the size at first maturity and L_max is the maximum recorded size of a species.

Results

The captured specimens of A. jacuhiensis (n = 370) showed standard lengths ranging from 1.8 to 9.2 cm. Figure 1a shows the weight/length relationship for the species (a = 0.01264; b = 3.464; residual Sum of Squares = 4.954) according to Huxley's traditional model. It is apparent that the species has two independent spot clusters, which characterize at least two age groups, and no fish were captured with standard lengths around 4 cm SL. Figure 1b represents the proportional residuals resulting from the Huxley adjustment for A. jacuhiensis. The proportional residuals of animals smaller than 4 cm were predominantly negative, clearly demonstrating the inadequacy of the applied model. There was also a concentration of positive residues in animals with lengths between 4 and 6 cm SL, and a new concentration of negative residues in animals larger than 6 cm SL, following a complex diagonal pattern. With distinct trends but with no obvious indication of a breakdown in the residual pattern, the Stanza Changing Length (SCL) was initially adjusted using 4.0 as the seed value.

 

Figures 1c-d show, respectively, the weight/length relationship and distribution of proportional residuals according to the polyphasic model for A. jacuhiensis. The adjusted parameters with their confidence intervals (1.96 SE) are presented in Table 2. The Stanza Changing Length (SCL) was estimated as 3.841 cm. Below this size the species shows a positive allometric growth pattern (b₁ = 3.112), changing to negative allometry for larger animals (b₂ = 2.939). Although at least two age groups were clearly identified as two data clouds (Fig. 1d), each group was well distributed around the zero line, which represents the predicted function. Nevertheless, the first group was still unbalanced, perhaps because of asymmetrical data distribution, as there was no pattern related to size. For the second group, even considering a well-balanced residual distribution for animals with lengths between 4 and 8 cm SL, large animals (> 8 cm SL) are clearly underestimated concerning weight to length. R² for both models was close to 0.99, indicating that this statistical parameter is not robust to identify minor differences in the model's descriptive power. On the other hand, the Proportional Residuals Sum of Squares decreased from 4.954 to 2.537 in the polyphasic model, indicating a better adjustment.

The 701 collected specimens of C. ibicuhiensis ranged from 1.7 to 4.5 cm SL. The weight/length relationship according to the Huxley model is shown in Fig. 2a (a = 0.01600; b = 3.287; residual Sum of Squares = 3.925). Although the trend line appears to fit the data set, fish longer than 4 cm SL clearly showed an overestimated weight (Fig. 2a). Proportional residuals derived from Huxley's adjustment for C. ibicuhiensis are shown in Fig. 2b. Proportional residuals are distributed apparently in a random pattern, except for larger animals (> 4 cm SL), which had negative residuals. Besides this anomaly, a clear breakdown in the residual pattern was not identified for C. ibicuhiensis. With no obvious indication, the Stanza Changing Length (SCL) was initially adjusted using 3.0 as the seed value. This value was chosen because of the bottleneck in the residual distribution as identified through Fig. 2b, suggesting a concentration of data derived from at least two age groups.

 

The weight/length relationship and the distribution of proportional residuals according to the polyphasic growth model are presented in Figs. 2c-d respectively. The adjusted parameters are given in Table 2. The SCL was estimated as 2.948 cm. Below this size, the species shows a positive allometric growth pattern (b₁ = 3.319), changing to an unexpected very slight negative allometry for larger animals (b₂ = 1.224), following an almost straight line but correcting the identified anomaly of larger animals (Fig. 2d). Determination Coefficients (R²) were not informative, close to 0.99 for both models. The Sum of Squares of Proportional Residuals decreased slightly in the polyphasic model, but no more than 10%.

Figure 3

Discussion

A change in the growth pattern of A. jacuhiensis and C. ibicuhiensis were identified through the weight/length relationship. Is this change related to the attainment of sexual maturity? Concerning A. jacuhiensis, there are no published values for L₅₀. The general model proposed by Binohlan & Froese (2009) estimates L₅₀ for A. jacuhiensis as 5.8 cm SL if used the maximum species size of 9.6 cm SL of the present work. This value is well above the 3.8 cm estimate for the Stanza Changing Length (SCL) of the polyphasic growth model, and direct estimates of L₅₀ for the species are necessary.

For C. ibicuhiensis, mature females are found from 2.7 (Oliveira et al., 2002) to 3.2 cm SL Braun et al., 2000). These values are compatible with the Stanza Changing Length estimated by the polyphasic growth model (SCL = 2.948 cm). This suggests that, similarly to M. funieri (Bervian et al., 2006), a change in the growth pattern of C. ibicuhiensis can be related to a change in energy allocation related to the start of reproductive activities. It is also interesting the maturity? Concerning A. jacuhiensis, there are no published For C. ibicuhiensis, mature females are found from 2.7 values for L₅₀. The general model proposed by Binohlan & (Oliveira et al., 2002) to 3.2 cm SL Braun et al., 2000). These Froese (2009) estimates L₅₀ for A. jacuhiensis as 5.8 cm SL if values are compatible with the Stanza Changing Length used the maximum species size of 9.6 cm SL of the present estimated by the polyphasic growth model (SCL = 2.948 cm). work. This value is well above the 3.8 cm estimate for the This suggests that, similarly to M. funieri (Bervian et al., Stanza Changing Length (SCL) of the polyphasic growth 2006), a change in the growth pattern of C. ibicuhiensis can model, and direct estimates of L₅₀ for the species are be related to a change in energy allocation related to the necessary. start of reproductive activities. It is also interesting the negative allometry of the second stanza of C. ibicuhiensis (b₂ = 1.224). According to Oliveira et al. (2002), C. ibicuhiensis has a long reproductive period from September through February, although Braun et al. (2000) reported mature females throughout the year. The species has also partial spawning (Oliveira et al., 2002), suggesting high energy allocation for reproduction, and perhaps decreased energy for somatic growth.

Although the L₅₀ /SCL association is apparently clear for C. ibicuhiensis, other biological factors could be related to the observed change in the growth pattern of A. jacuhiensis. Unfortunately, there is no published information concerning the biology of A. jacuhiensis that could explain the identified change in the growth pattern. Even the two-stanza growth model does not perfectly describe the full data series, as positive residuals above 8 cm suggest the need for a third growth stanza. Because of the few data available for this size (only nine animals), we decided not to model a more complex function.

From Fig. 3, it is apparent that the predicted trend line (thin line) of Binohlan & Froese (2009) is well above the L₅₀ estimates for most small characids. In reality, this model is a correction of an older one proposed by these same authors (Froese & Binohlan, 2000), with inflated L₅₀ values especially for small species. Nevertheless, for small characids, the empirical linear relationship of L₅₀ as half of the maximum size is also simple and apparently fits better than the general proposition for all teleosts, as the Binohlan & Froese (2009) model, explaining 68% of L₅₀ variability. This "rule", however, requires corroboration from a larger data set. Applied to A. jacuhiensis, this pattern indicates an L₅₀ of 4.8 cm, only one centimeter above the estimated Stanza Changing Length (3.8). In this case, a relationship between the estimated change in the growth pattern (SCL) and the attainment of sexual maturity of A. jacuhiensis could only be matched if, in some way, the species anticipates the change in the growth pattern as a way to allocate energy for the demands of reproduction.

Although not conclusive, our data suggest that a complex growth pattern is frequent in the natural environment, and perhaps is not usually identified because small trends are obscured by natural variability. Nevertheless, there is a tradeoff when using sophisticated models: any additional parameter can improve the adjustment at a cost of increased complexity. In this respect, there are several methods to support model selection, as the adjusted R², the likelihood ratio test, the Akaike information criterion, or the Schwarz criterion (for revision in this subject, see Johnson & Omland, 2004). Despite the increased complexity, the key factor by choosing to apply the polyphasic model is not strictly statistical, but the possibility to identify a change in the growth pattern that could indicate relevant biological aspects, such as sexual maturity. In this respect, the increased complexity of the polyphasic growth model should not be viewed as a way to merely describe the weight/length relationship. The practical gain in terms of adjustment quality could be minimal. The key factor is that some important biological facts could be reflected in the growth trend. Besides the attainment of sexual maturity, a diet shift, a change of habitat occupation, or the start of a migratory process are all biological facts that may cause a change in the energy budget of a species, what could be reflected in the weight/length relationship. These biological facts, however, could not be identified if simple data adjustments are applied. Therefore, the application of a more complex function, as polyphasic growth model, could point out to unknown anomalies and inspire new insights, guiding research programs to understand fish biology in detail.

Concluding, can weight/length relationship predict size at first maturity? The present study supports our previous hypothesis that the energetic demands related to reproduction could be reflected in the growth pattern of fish species, as for C. ibicuhiensis and maybe for A. jacuhiensis. Yet, any change in the growth pattern could also be reflecting a neglected aspect of fish biology, and the application of the proposed methodology should be encouraged as an additional tool to depict the complex biology of fish.

Literature Cited

Accepted August 25, 2010

• Andrade, V. X. L. 2003. Ictiofauna da represa municipal de São José do Rio Preto, SP, Rio Preto, afluente do Rio Turvo, drenagem do Rio Grande: diversidade e reprodução. Unpublished M.Sc. Dissertation. Universidade Estadual Paulista, São Paulo, 53p.
• Bervian, G., N. F. Fontoura & M. Haimovici. 2006. Statistical model of variable allometric growth: otolith growth in Micropogonias funieri (Actinopterygii, Sciaenidae). Journal of Fish Biology, 68: 196-208.
• Binohlan, C. & R. Froese. 2009. Empirical equations for estimating maximum length from length at first maturity. Journal of Applied Ichthyology, 25(5): 611-613.
• Bolger, T. & P. L. Connolly. 1989. The selection of suitable indices for the measurement and analysis of fish condition. Journal Fish Biology, 34: 171-182.
• Braun, A. S., D. S. Lewis & N. F. Fontoura. 2000. Biologia reprodutiva de Cheirodon ibicuhiensis (Eigenmann, 1915) na lagoa Fortaleza, Cidreira, Rio Grande do Sul, Brasil (Teleostei: Characidae: Cheirodontinae). Comunicações do Museu de Ciências e Tecnologia PUCRS, Série Zoologia, 13(2): 159-166.
• Fontoura, N. F., A. S. Braun & P. C. C. Milani. 2009. Estimating size at first maturity (L₅₀) from Gonadossomatic Index (GSI) data. Neotropical Ichthyology, 7(2): 217-222.
• Froese, R. & C. Binohlan. 2000. Empirical relationships to estimate asymptotic length, length at first maturity and length at maximum yield per recruit in fishes, with a simple method to evaluate length frequency data. Journal of Fish Biology, 56: 758-773.
• García-Alzate, C. A. & C. Román-Valencia. 2008. Biologia alimentaria y reproductiva de Hyphessobrycon poeciloides (Pisces: Characidae) en la cuenca del río La Vieja, Alto Cauca, Colombia. Revista del Museo Argentino de Ciencias Naturales, 10: 17-27.
• Huxley, J. S. 1924. Constant differential growth-ratios and their significance. Nature, 14: 896-897.
• Johnson, J. B. & K. S. Omland. 2004. Model selection in ecology and evolution. Trends in Ecology and Evolution, 19(2): 101-108.
• Le Cren, E. D. 1951. The length-weight relationship and seasonal cycle in gonad and conditions in the perch Perca fluviatilis Journal of Animal Ecology, 20(2): 201-219.
• Lima, F. C. T., L. R. Malabarba, P. A. Buckup, J. F. Pezzi da Silva, R. P. Vari, A. Harold, R. Benine, O. T. Oyakawa, C. S. Pavanelli, N. A. Menezes, C. A. S. Lucena, M. C. S. L. Malabarba, Z. M. S. Lucena, R. E. Reis, F. Langeani, L. Casatti, V. A. Bertaco, C. Moreira & P. H. F. Lucinda. 2003. Genera Incertae sedis in Characidae. Pp. 106-169. In: Reis, R. E., S. O. Kullander & C. J. Ferraris Jr. (Eds.). Check List of the Freshwater Fishes of South and Central America. Porto Alegre, Edipucrs, 729p.
• Marques, C. S., A. S. Braun & N. F. Fontoura. 2007. Estimativa de tamanho de primeira maturação a partir de dados de IGS: Oligosarcus jenynsii, Oligosarcus robustus, Hoplias malabaricus, Cyphocharax voga, Astyanax fasciatus (Characiformes), Parapimelodus nigribarbis, Pimelodus maculatus, Trachelyopterus lucenai, Hoplosternum littorale, Loricariichthys anus (Siluriformes) e Pachyurus bonariensis (Perciformes) no lago Guaíba e laguna dos Patos, RS. Biociências, 15(2): 230-256.
• Mazzoni, R, R. S. Mendonça & E. P. Caramaschi. 2005. Reproductive biology of Astyanax janeiroensis (Osteichthyes, Characidae) from the Ubatiba river, Maricá, RJ, Brazil. Brazilian Journal of Biology, 65(4): 643-649.
• Oliveira, C. L. C., C. B. Fialho & L. R. Malabarba. 2002. Período reprodutivo, desova e fecundidade de Cheirodon ibicuhiensis Eigenmann, 1915 (Ostariophysi: Characidae) do arroio Ribeiro, Rio Grande do Sul, Brasil. Comunicações do Museu de Ciências e Tecnologia PUCRS, Série Zoologia, 15(1): 3-14.
• Román-Valencia, C. & R. I. Ruiz. 2005. Diet and reproduction aspects of Astyanax aurocaudatus (Teleostei: Characidae) from the upper part of the Río Cauca, Colombia. Dahlia, 8: 9-17.
• Smith, R. J. 1980. Rethinking allometry. Journal of Theoretical Biology, 87: 97-111.
• Vazzoler, A. E. A. de M. 1996. Biologia reprodutiva de peixes teleósteos: teoria e prática. Maringá, Eduem, 169p.
• Veregue, A. M. L & M. L. Orsi. 2003. Biologia reprodutiva de Astyanax scabripinis paranae (Eigenmann)(Ostheichthyes, Characidae), do ribeirão das Marrecas, bacia do rio Tibagi, Paraná. Revista Brasileira de Zoologia, 20(1): 97-105.

Publication Dates

• Publication in this collection
  15 Oct 2010
• Date of issue
  2010

History

• Accepted
  25 Aug 2010