Introduction Methods Results Discussion Acknowledgements References About the Author(s) Clement J. Arendse Rondevlei Scientific Services, Garden Route National Park, South African National Parks (SANParks), Sedgefield, South Africa Ian A. Russell Rondevlei Scientific Services, Garden Route National Park, South African National Parks (SANParks), Sedgefield, South Africa Citation Arendse, C.J. & Russell, I.A., 2020, ‘Morphometric regressions for the endangered Knysna seahorse, Hippocampus capensis, in the Swartvlei Estuary from mass stranding events’, Koedoe 62(1), a1632. https://doi.org/10.4102/koedoe.v62i1.1632 Short Communication Morphometric regressions for the endangered Knysna seahorse, Hippocampus capensis, in the Swartvlei Estuary from mass stranding events Clement J. Arendse, Ian A. Russell Received: 08 Apr. 2020; Accepted: 22 July 2020; Published: 28 Sept. 2020 Copyright: © 2020. The Author(s). Licensee: AOSIS. This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Introduction The Knysna seahorse, Hippocampus capensis, is listed as endangered on the IUCN Red List of Threatened Species, with the only known, apparently isolated, populations occurring in the Knysna, Swartvlei and Keurbooms estuaries along the south coast of South Africa. The Swartvlei Estuary population is the most genetically distinct (Mkare et al. 2017). Despite all three known populations occurring within protected areas and there being no commercial or subsistence use of the species, population sizes are still thought to be declining because of habitat degradation across much of its range (Pollom 2017). Environmental management authorities have outlined strategies for the conservation of H. capensis, including listing it as a species of special concern in the Estuarine Lower Level Plan of the Garden Route National Park (South African National Parks 2020) and setting management objectives for the Keurbooms population in the Keurbooms–Bitou Estuarine Management Plan (Western Cape Government 2018). Assessments of population size and habitat use of H. capensis have focussed mainly on the relatively larger and more stable population in Knysna Estuary (Bell et al. 2003; Claassens 2016; Claassens & Hodgson 2018; Teske et al. 2007). The lesser studied Swartvlei Estuary population is subject to large fluctuations in abundance, with past densities ranging from 0.04 individuals per m2 to 0.26 individuals per m2 (Lockyear et al. 2006). The Swartvlei Estuary is the only intermittently open host estuary, with elevated water levels and prolonged flooding of land adjacent to the shores common during closed periods. Periodic mass strandings and associated die-offs of H. capensis have been recorded in Swartvlei Estuary when the estuary is breached (Russell 1994; Skelton 1987), mainly by artificial means to prevent inundation of properties on the estuarine floodplain. The breaching protocol and criteria are contained with the Garden Route National Park Management Plan (SANParks 2020). South African National Parks (SANParks) as the management authority of Swartvlei Estuary has since 1991 intermittently collected dead H. capensis from the banks of the Swartvlei Estuary after breaching events. This now regular sampling has evolved into a SANParks co-ordinated citizen science programme. The only known morphometric descriptions of H. capensis come from small samples of flood-stranded individuals from the Knysna Estuary (n = 57) and Swartvlei Estuary (n = 39) (Toefy 2000), despite the usefulness of such data in allowing comparison between studies and the conversion of data where different measures of length are used (Lourie 2003). This article aims to use the seahorse morphometric data collected following two mass stranding events (> 500 individuals) in 1991 and 2017 to provide more robust estimates of the relationship between key morphometric characteristics of H. capensis. In addition, data from mass strandings are compared to those collected from seine netting of live individuals in 2019 to determine the feasibility of using these periodic events as a means of assessing population attributes. Methods Sample sites Hippocampus capensis were sampled from the Swartvlei Estuary (34° 01’ 48.58” S, 22° 47’ 50.66” E), situated along the south coast of South Africa. Collection protocol Dead H. capensis were collected after mass stranding events in February 1991 and November 2017. Specimens collected in 1991 were preserved in 10% neutral freshwater-buffered formalin immediately after collection (Russell 1994), whereas in 2017 dead individuals were frozen immediately after collection. In 2017, samples were collected within 1 day of stranded individuals first being noticed, and for 1 week after estuary breaching. Sampling effort was spread throughout the estuary, although repeat sampling, where sites were searched daily over the sampling period, was undertaken mostly in areas with previous high abundance (Lockyear et al. 2006). Approximately 1000 individuals were recorded, of which half (52%) were dead. All live individuals were returned to the water. Retained dead individuals showed various levels of desiccation, likely a result of differences in exposure time because of the intermittent stranding of H. capensis over several days following breaching. Samples were pooled for all sites, and all samples were preserved without manipulating (i.e. straightening) body parts. Live H. capensis were collected in March 2019 using a 10 m × 1.5 m seine net with a stretched mesh size of 8 mm. Three effective hauls were undertaken in ≤ 1.3 m water depth in areas of previous high H. capensis density, with each haul being approximately 3 m wide, ranging from 12 m to 21 m long, and spaced at least 20 m apart. Data were pooled for all replicates. Live samples were compared to preserved samples to determine if size classes and sex ratios differed between sampling techniques. Morphometrics Preserved samples were photographed and measured using ImageJ software. Where necessary, individuals were flattened before being photographed to minimise bias in measurements because of body arching (discussed in Sotola et al. 2019). All individuals were sexed and then measured using the protocol set by Lourie (2003). Head length (HdL), curved trunk length (TrL), curved tail length (TaL) and height (Ht) were measured, with standard length (SL) calculated as the sum of the first three measurements. Straight-line SL (SLstraight) and straight-line Ht (Htstraight) were also determined for a subsample of individuals whose heads were positioned at a right angle to their trunk. The trunk measurements for the straight-line sizes were measured as straight-line distances from the respective anterior starting points ending at the last trunk ring. Regression equations of size measurements of H. capensis included measurement of Htstraight and Ht because Htstraight is often the size measurement used when sampling live specimens (Foster & Vincent 2004; Lourie 2003), likely to avoid handling stress, and Ht may be useful if the snout of a dead seahorse is damaged. Live individuals collected in 2019 were sexed and Htstraight was measured before the animals were released at the location of capture. The effect of preservation techniques on morphometric relationships was investigated by comparing the regression equations between SL and HdL, TrL, TaL and Htstraight of formalin preserved (collected 1991) and frozen (collected 2017) specimens. Preservation of samples in a fixative can cause shrinkage, and different preservation techniques may result in different shrinkage rates of seahorse body parts (Nadeau, Curtis & Lourie 2009). Also, further shrinkage is known to occur in fish when multiple preservation techniques, in this case freezing or preservation in formalin after desiccation, are applied (Sotola et al. 2019). Statistical analysis Sex ratio bias within the 1991, 2017 and 2019 datasets was tested using an exact binomial test (Crawley 2007; Zar 2010). Sexes were analysed separately because of known sexual dimorphism (Bell et al. 2003). Normality of SL (between preservation techniques) and Htstraight (between alive and preserved specimens) measurements for years and sexes were tested using a Shapiro–Wilk test. Homogeneity of variances was tested using a Fligner–Killeen test before differences in means were tested using either a one-tailed Welch’s two-sample t-test or a one-tailed Wilcoxon rank-sum test. Differences in mean Htstraight between years were tested using a Kruskal–Wallis rank-sum test. Regression equations of SL and Htstraight and their component measurements for years and sexes were produced separately using ordinary least squares regressions and compared to determine if sex and preservation techniques influenced correlation between variables. Regression equations of SL and Htstraight, separated by year and sex, were tested using an analysis of covariance (ANCOVA) to determine if regression equations obtained differed significantly (Crawley 2007; Zar 2010). Data analyses were undertaken using the programming language R (R Core Team 2019). Results Sex ratios Sex ratios were biased towards females in the 1991 sample (1.025:1; n = 170) and males in the 2017 and 2019 samples (0.974:1; n = 304 and 0.814:1; n = 78, respectively), although none differed significantly from the predicted 1:1 ratio (p > 0.05). Size class distribution Mean length of H. capensis collected in 2017 (62 mm – 100 mm SL; n = 274) was significantly larger than the mean length of those collected in 1991 (33 mm – 100 mm SL; n = 161) when sexes were grouped (W = 7800; p < 0.05) as well as for females (W = 1650; p < 0.05) and males (W = 2227; p < 0.05) separately. The majority of individuals collected in 2017 were above the size of sexual maturity (≥ 65 mm SL [Whitfield 1995]), whereas most individuals collected in 1991 were likely sexually immature (Figure 1a). Mean SL differed significantly between sexes in both 1991 (W = 2255; p < 0.05) and 2017 (t = −2.989; p < 0.05), with the mean size of males being larger than that of females in both years. FIGURE 1: Box and whisker plots of Hippocampus capensis separated by sex showing the minimum, first quartile, median, third quartile and maximum sizes of (a) SL for preserved (1991 and 2017) samples, and (b) Htstraight for preserved (1991 and 2017) and live samples (2019). Measurements of Htstraight from live H. capensis collected in 2019 were normally distributed for both females (W = 0.946; p > 0.05) and males (W = 0.984; p > 0.05). Mean Htstraight did not differ significantly between sexes (t = −0.432; p > 0.05). The 2019 size range (38 mm – 71 mm Htstraight; n = 78) was similar to those collected in both 1991 (35 mm – 76 mm Htstraight; n = 53) and 2017 (54 mm – 77 mm Htstraight; n = 97), though mean sizes differed significantly for both females (chi-squared = 51.875; p < 0.05) and males between years (chi-squared = 52.689; p < 0.05) (Figure 1b). Morphometric regressions The fit of regression models between SL and HdL, TrL and TaL was better for the 1991 samples compared to those collected in 2017, as evidenced by the high R2 values in the former (Tables 1 and 2). Test statistics for the regressions between SL and HdL, TrL, TaL and Htstraight are presented in Table 3. HdL, TrL and TaL were significantly correlated to SL for both sexes when years were analysed separately, whilst sex effected all regressions except HdL in 1991. Interaction did not significantly effect the fit of the models in either year (p > 0.05) and the slope of the models did not differ significantly between sexes (p > 0.05), with the exception of HdL in 2017. When differences between years were calculated separately for sexes, interaction effects were evident for the female datasets for TrL and TaL and male datasets for HdL and TrL (p > 0.05) but not the female dataset for HdL or the male dataset for TaL. Year significantly effected all regressions except the male dataset for TrL. Slopes differed between years for all variables for both sexes. Preservation technique, therefore, influenced TrL for both sexes, whereas HdL only differed between males and TaL between females. When years were analysed individually, Htstraight significantly effected SL in 1991 and 2017 (p < 0.05), but sex effected the regressions only in 2017 (p < 0.05). No interaction effects were evident in either year (p > 0.05). When sexes were analysed separately, Htstraight and year significantly influenced the regressions for both sexes, with no interaction effects evident for females or males (p > 0.05). The linear regression for males collected in 1991 (Table 1) was used to calculate Htstraight at sexual maturity (56 mm) (Figure 1b). TABLE 1: Correlation between selected morphometric characteristics of Hippocampus capensis collected in 1991. TABLE 2: Correlation between selected morphometric characteristics of Hippocampus capensis collected in 2017. TABLE 3: Test statistics of regression models of standard length with head length, curved trunk length and curved tail length, showing interaction between the measured variables and sexes and years, respectively. Discussion The 1991 and 2017 datasets provided reliable correlations between the morphometric variables measured for the different preservation techniques. Both preservation techniques used in this study are known to result in significant shrinkage of fish bodies (Berbel-Filho, Jacobina & Martinez 2013; Buchheister & Wilson 2005; Nadeau et al. 2009) with the preservation of samples in formalin and ethanol resulting in ongoing shrinkage of body parts (Martinez, Berbel-Filho & Jacobina 2013; Sotola et al. 2019). To compound the complexity of using preserved specimens in morphometric investigations, shrinkage in fish may not be uniform across body parts and differences can even be observed at an individual level (Sotola et al. 2019). This effect was noted in this study with the regressions of SL with HdL, TrL and TaL of H. capensis collected in 2017, showing greater variation than those collected in 1991. Although the individuals collected in 2017 were frozen before processing, they showed differing levels of desiccation when collected. This, coupled with an extended, staggered collection period for the stranded seahorses, and the relatively shorter period between preservation and processing of samples compared to the 1991 dataset, could have contributed to the larger variation in the body length ratios. The greater correlation between morphometric measurements of SL and HdL, TrL and TaL in 1991 was therefore likely because of more even shrinkage resulting from prolonged exposure of the samples to the fixative. This result, together with the differences in slopes between the 2 years, suggests that different regression equations are needed for each preservation type and that comparisons between samples preserved in different fixatives may require correction factors to be calculated by measuring samples before and at several periods after preservation. In addition, using regression equations obtained from preserved specimens to describe morphometric relationships of live specimens is not recommended. The mean length (Htstraight) of H. capensis differed significantly between 1991, 2017 and 2019 for both sexes. Differences could relate to the timing of the breaching events, with both the 1991 and 2019 samples collected towards the end of the breeding period for H. capensis (austral summer – Whitfield [1995]), whereas the 2017 samples were collected at the start of the breeding period, potentially explaining the lower proportion of smaller individuals in that year. The duration of preservation of specimens could also have contributed to differences in mean length, with specimens collected in 1991 being subjected to longer exposure to a fixative likely having shrunk to a greater degree. Although differences in the size data existed, the slopes of the SL versus Htstraight regressions for 1991 and 2017 did not differ significantly from one another, indicating that the relationship between SL and Htstraight remains relatively constant regardless of preservation technique used. This was, however, not the case for the relationship between SL and HdL, TrL and TaL, where preservation technique, as well as sex, influenced the slopes of the regressions. As various length measurements are often used for measuring seahorses, standardisation of seahorse measurements is difficult (Foster & Vincent 2004). The morphometric regressions presented in Tables 1 and 2 can assist in comparing the results of previous studies where different measurements and preservation methods were used. Equal sex ratios observed in 1991 and 2017 were similar to those found in in situ studies of natural populations (Lockyear et al. 2006) and in 2019, suggesting that stranding affected both sexes equally. The use of stranding data for H. capensis in the Swartvlei Estuary could be a novel way of post hoc determination of population attributes, as data collected from these events are comparable to those collected from the study of in situ populations. The effects of preservation on morphometric variables also need additional consideration, and shrinkage correction equations should be developed for different preservation methods and durations. Acknowledgements The authors would like to thank the citizen scientists, SANParks staff and honorary rangers for their assistance with data collection. They would also like to thank Dr Mohlamatsane Mokhatla for his help with the statistical programme R. Competing interests The authors have declared that no competing interest exists. Authors’ contributions C.J.A. and I.A.R. conceptualised the study, conducted fieldwork and contributed to writing the article. C.J.A. performed all laboratory and statistical analysis. Ethical consideration This study fully meets the ethical standards for handling of Hippocampus capensis as set out by the SANParks Animal Use and Care Committee (004/17). Funding information This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors. Data availability statement Data can be requested via the SANParks data repository (http://dataknp.sanparks.org/sanparks/). Data use limitations may apply. Disclaimer The views and opinions expressed in this article are those of the authors and do not necessarily reflect the official policy or position of any affiliated agency of the authors. References Bell, E.M., Lockyear, J.F., McPherson, J.M., Marsden, A.D. & Vincent, A.C.J., 2003, ‘First field studies of an endangered South African seahorse, Hippocampus capensis’, Environmental Biology of Fishes 67(1), 35–46. https://doi.org/10.1023/A:1024440717162 Berbel-Filho, W.M., Jacobina, U.P. & Martinez, P.A., 2013, ‘Preservation effects in geometric morphometric approaches: Freezing and alcohol in a freshwater fish’, Ichthyological Research 60(3), 268–271. https://doi.org/10.1007/s10228-013-0339-x Buchheister, A. & Wilson, M.T., 2005, ‘Shrinkage correction and length conversion equations for Theragra chalcogramma, Mallotus villosus and Thaleichthys pacificus’, Journal of Fish Biology 67(2), 541–548. https://doi.org/10.1111/j.0022-1112.2005.00741.x Claassens, L., 2016, ‘An artificial water body provides habitat for an endangered estuarine seahorse species’, Estuarine, Coastal and Shelf Science 180, 1–10. https://doi.org/10.1016/j.ecss.2016.06.011 Claassens, L. & Hodgson, A.N., 2018, ‘Gaining insights into in situ behaviour of an endangered seahorse using action cameras’, Journal of Zoology 304(2), 98–108. https://doi.org/10.1111/jzo.12509 Crawley, M.J., 2007, The R Book, John Wiley & Sons, Ltd, The Atrium, Southern Gate, Chichester, West Sussex, England. Foster, S.J. & Vincent, A.C.J., 2004, ‘Life history and ecology of seahorses: Implications for conservation and management’, Journal of Fish Biology 65(1), 1–61. https://doi.org/10.1111/j.0022-1112.2004.00429.x Lockyear, J.F., Hecht, T., Kaiser, H. & Teske, P.R., 2006, ‘The distribution and abundance of the endangered Knysna seahorse Hippocampus capensis (Pisces: Syngnathidae) in South African estuaries’, African Journal of Aquatic Science 31(2), 275–283. https://doi.org/10.2989/16085910609503897 Lourie, S., 2003, Measuring seahorses, Project seahorse technical report no. 4, version 1.0, p. 15, Project seahorse, Fisheries Centre, University of British Columbia. Martinez, P.A., Berbel-Filho, W.M. & Jacobina, U.P., 2013, ‘Is formalin fixation and ethanol preservation able to influence in geometric morphometric analysis? Fishes as a case study’, Zoomorphology 132(1), 87–93. https://doi.org/10.1007/s00435-012-0176-x Mkare, T.K., Jansen van Vuuren, B. & Teske, P.R., 2017, ‘Conservation implications of significant population differentiation in an endangered estuarine seahorse’, Biodiversity and Conservation 26(6), 1275–1293. https://doi.org/10.1007/s10531-017-1300-5 Nadeau, J.L., Curtis, J.M.R. & Lourie, S.A., 2009, ‘Preservation causes shrinkage in seahorses: Implications for biological studies and for managing sustainable trade with minimum size’, Aquatic Conservation: Marine and Freshwater Ecosystems 19(4), 428–438. https://doi.org/10.1002/aqc.1002 Pollom, R.A., 2017, ‘Hippocampus capensis’, The IUCN Red list of threatened species 2017, e.T10056A54903534. R Core Team, 2019, R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, viewed 07 October 2019, from https://www.R-project.org/. Russell, I.A., 1994, ‘Mass mortality of marine and estuarine fish in the Swartvlei and Wilderness lake systems, Southern Cape’, Southern African Journal of Aquatic Sciences 20(1–2), 93–96. https://doi.org/10.1080/10183469.1994.9631353 Skelton, P.H., 1987, South African Red Data Book – Fishes, Pretoria, South Africa. Sotola, V.A., Craig, C.A., Pfaff, P.J., Maikoetter, J.D., Martin, N.H. & Bonner, T.H., 2019, ‘Effect of preservation on fish morphology over time: Implications for morphological studies’, PLoS One 14(3), e0213915. https://doi.org/10.1371/journal.pone.0213915 South African National Parks, 2020, Garden Route National Park Management Plan, p. 222, Groenkloof, South Africa. Teske, P.R., Lockyear, J.F., Hecht, T. & Kaiser, H., 2007, ‘Does the endangered Knysna seahorse, Hippocampus capensis, have a preference for aquatic vegetation type, cover or height?’, African Zoology 42(1), 23–30. https://doi.org/10.3377/1562-7020(2007)42[23:DTEKSH]2.0.CO;2 Toefy, Z., 2000, ‘A preliminary examination of the morphological and genetic structure within populations of the Knysna seahorse Hippocampus capensis (Pisces: Syngnathidae) from two South African estuaries’, M.Sc dissertation, University of the Western Cape, Bellville, South Africa, p. 83. Western Cape Government, 2018, Keurbooms-Bitou Estuarine Management Plan, Cape Town, South Africa, p. 139. Whitfield, A.K., 1995, ‘Threatened fishes of the world: Hippocampus capensis Boulenger, 1900’, Environmental Biology of Fishes 44(2), 362. https://doi.org/10.1007/BF00008251 Zar, J.H., 2010, Biostatistical analysis, 5th edn., Prentice Hall, Inc, Upper Saddle River, New Jersey.