Preprint
Article

Condition, Reproductive Fitness, and Fluctuating Asymmetry in Brook Stickleback: Responses to Anthropogenic Runoff

Altmetrics

Downloads

61

Views

34

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

24 October 2023

Posted:

26 October 2023

You are already at the latest version

Alerts
Abstract
Multiple indicators have been used to assess the degree of exposure of fish to anthropogenic chemicals in their stream habitats. We hypothesized that brook stickleback (Culaea inconstans) in a headwater stream receiving urban and agricultural runoff (South Fork Whitewater River, SFWR) would exhibit poorer condition, reduced reproductive fitness, and greater left side-to-right side morphological asymmetry (i.e., fluctuating asymmetry or FA) than would fish from a nearby headwater stream with a forested drainage basin (Garvin Brook). Male and female fish were collected from both streams just prior to spawning in 2013-2015. In 2013 and 2014, fish were assessed for overall condition (Fulton’s K) and internal measures of condition (hepatosomatic index, HSI) and reproductive fitness (gonadosomatic index [GSI], total egg count, egg weight). In 2015, measurements of head length, jaw length, eye diameter, pectoral fin length, and pelvic fin length were made on both sides of each fish for use in assessing degree of FA. We observed declining condition with fish size, increased liver size, and reduced egg counts and egg size in female brook stickleback in SFWR relative to those from Garvin Brook. SFWR females had significantly higher FA than Garvin females for all structures assessed except pelvic fin length. FA also was slightly higher for all structures in SFWR males compared to Garvin males, but differences were not significant. A composite FA index combining all measurements from an individual fish into a single value displayed highly significant differences for female fish (SFWR FA >> Garvin FA), but not for male fish (SFWR FA = Garvin FA). Exposure of brook stickleback to reduced water quality in SFWR during early development appears to increase morphological asymmetry in female (but not male) fish, and continuing exposure to compromised water quality throughout life impacts both general condition and reproductive fitness of stickleback, especially older female fish, in SFWR.
Keywords: 
Subject: Biology and Life Sciences  -   Aquatic Science

1. Introduction

Poor water quality can have significant effects on the health of aquatic organisms. Severe toxic spills, leaks, and releases of many substances can quickly kill entire communities of organisms in lakes, streams, and rivers [1,2,3], whereas milder conditions or events may cause more subtle, less obvious, but longer-term changes in organism behavior or fitness [4,5,6,7]. A wide range of human activities (e.g., logging, agriculture, mining, manufacturing, construction, general urbanization) within watersheds have been shown to have negative impacts on many different types of aquatic organisms [4,7,8,9].
Like many regions worldwide, numerous lakes, streams, and rivers in Minnesota, USA, are affected negatively by human activities occurring within their respective drainage basins. These activities have resulted in widespread impacts to surface waters, with Minnesota currently listing 2,904 water bodies as impaired (with 6,168 total impairments) due to compromised water quality [10]. Impairments encompass a wide range of water quality standards, including suspended sediments that reduce water clarity, excessive nutrients that increase algal production, bacteria that make water unsafe for contact or consumption, chemicals (e.g., mercury, PFOS) that can contaminate fish tissue and be passed on to humans when consumed, and generally unhealthy conditions for fish and aquatic macroinvertebrates [10].
Impaired water quality in streams and rivers may not directly kill aquatic fauna or prevent them from inhabiting certain systems, but their overall health and fitness may be reduced significantly [5,7,9,11]. Poor health and fitness in stream fish may be evident externally as small size caused by slower growth, lower weight at length, increased prevalence of DELT (deformities, eroded fins, lesions, tumors) anomalies, or morphological asymmetry (i.e., fluctuating asymmetry, FA) due to altered early development [5,11,12]. Internally, poor health and fitness in fish may be expressed as enlarged livers and smaller reproductive organs [5,12]. Such indicators of stress often are useful in detecting the presence of many pollutants at non-lethal, but chronic, levels [4,5,7].
In this study, we examined fish in two headwater streams in southeastern Minnesota to assess how differing degrees of water quality impairment may impact their fitness and reproductive condition. Specifically, we measured internal and/or external characteristics in brook stickleback (Culaea inconstans) just prior to their spawning season in three different years (expecting fish to exhibit peak fitness during these periods) to see if such variables responded to varying levels of stream impairment, and if they may be used to assess reduced water quality in other systems. Brook stickleback are broadly distributed across the northern United States and southern Canada [13], and have been shown to be physiologically and morphologically responsive to various chemical pollutants [11,14,15]. We expected that stickleback from a stream with more confirmed impairments resulting from agriculture and urbanization within the watershed would exhibit reduced condition and reproductive fitness and greater morphological asymmetry than fish collected from a nearby, but less impaired, stream.

2. Study sites

Brook stickleback were examined in two headwater streams located 35 km apart in adjacent, separate watersheds in southeastern Minnesota, USA. The South Fork Whitewater River (SFWR) site was located within the City of Eyota, a small (population size = 2047) rural community surrounded by agricultural lands and livestock rangelands (Figure 1). The stream originates as two spring seeps that create wetlands bordered by agriculture fields 2 km upstream from the fish collection site within the city. Agricultural fields are drained by subsurface drain lines that connect to the stream, as well as by grassed waterways that conduct surface runoff toward the stream during heavy precipitation events or spring snowmelt. Within the city, the SFWR receives runoff via storm sewers from city streets, parking lots, and other impervious surfaces. In addition to brook stickleback, which comprised 90% of the fish community at this site during our study, creek chub (Semotilus atromaculatus) and western blacknose dace (Rhinichthys obtusus) also were present. During other years, central stoneroller (Campstoma anomalum) and southern redbelly dace (Chrosomus erythrogaster) also have been collected from the SFWR study site.
The Garvin Brook site was located within public forest lands approximately 4 km downstream from the headwater (limestone) springs (Figure 1). Agricultural fields, rangelands, and several rural residences also are located in the watershed upstream from our collection site, with all agriculture and rangelands located on bluff tops elevated 100 m or more higher than the headwater springs. A forested riparian buffer 100 to 200 m wide borders both sides of the stream at the spring source and along the intermittent channel that extends for approximately 1 km above the headwater springs. The stream corridor between the spring source and our study site included public and private forest lands and a county park (comprised of forested lands and mowed turfgrass). Along with brook stickleback, brown trout (Salmo trutta), slimy sculpin (Cottus cognatus), and brook trout (Salvelinus fontinalis) were present at the Garvin Brook location.
We did not directly assess water quality at either stream site during our study, but inferred potential impacts based on observations of surrounding lands, publicly available water quality records and described impairments (https://www.pcs.state.mn.us/air-water-land-climate/water-quality-trends-and-data), consultations with governmental and agency personnel with knowledge of the stream sites, and previous studies within these waters. During the study period, the stream reach containing the SFWR site had four active water quality impairments (i.e., not meeting regional water quality standards) in effect: high fecal coliform bacteria count, high turbidity, reduced benthic macroinvertebrate community composition/structure (ranked as fair on a scale ranging from excellent to very poor), and reduced fish community composition/structure (fair). By comparison, the stream reach containing the Garvin Brook study site had a single active water quality impairment: high fecal coliform bacteria count. Several recent studies conducted on these two streams [3,16,17,18,19,20,21] further suggest that the SFWR has much poorer water quality than Garvin Brook.

3. Methods

Brook stickleback were collected from each study site during late-May and early-June, 2013-2015. These time periods were immediately prior to the onset of stickleback spawning in these waters. Fish were collected by a combination of backpack electrofishing and overnight sets of standard cylindrical minnow traps “baited” with activated chemical light sticks. Sticklebacks were over-anesthetized with a 100 mg/L solution of tricaine methanesulfonate (MS-222), fixed for 48 hours in an 8% formalin solution, and preserved and stored in 70% ethanol until examined.
Sticklebacks collected in 2013 and 2014 were used to assess fish condition and reproductive fitness. Fish were removed from the ethanol preservative, rinsed in running tap water for 1 to 2 minutes, then towel dried prior to weighing (wet weight, nearest 0.001 g) and measuring (total length [TL], nearest mm). Weights and TLs were used to calculate a Fulton’s condition factor (K, where K = [weight/TL3] * 100,000) for each individual fish [22]. Condition factors were compared between sites separately for each sex.
To assess reproductive fitness [23] and exposure to stressful environmental conditions [24], individual sticklebacks were dissected to remove gonads (ovaries in females, testes in males) and the liver. These organs were weighed (nearest 0.1 mg) to calculate a gonadosomatic index (GSI, where GSI = [mg gonad/mg total fish weight] * 100) and a hepatosomatic index (HSI, where HSI = [mg liver/mg total fish weight] * 100) for each fish. For female fish, egg number was estimated by doubling the egg count from a single ovary. Individual egg size (mass) was estimated by dividing combined ovary weights by total egg counts. After combining data from 2013 and 2014, ovary weights, testes weights, GSIs, HSIs, total egg counts, and egg sizes were compared between stream sites. Due to expected increases in gonad and liver weights, egg counts, and egg sizes with increasing fish size, we used analysis of covariance (ANCOVA) to control for fish size differences, to compare each of the fish length-internal fitness variable relationships between stream sites.
Stickleback collected in 2015 were anesthetized, preserved, rinsed, and towel-dried as described above. Fish TL (nearest mm) was measured, and sex was determined based on external sex characteristics (e.g., dark coloration of males, lighter coloration and swollen abdomens [mature ovaries] of females). FA was assessed for individual fish by taking duplicate measurements (nearest 0.1 mm) of five anatomical structures (lower jaw length, eye diameter, head length, pectoral fin length, pelvic fin length; Figure 2) on both sides of the fish. Duplicate measurements were averaged prior to analysis. Differences in structure size between sides were calculated (left value minus right value) and compared to a value of 0 to test for “ideal FA” (i.e., mean values for each test group compared to 0 using t-tests, with separate tests for each structure/sex/stream). To compensate for variation in ideal FA among different structures (i.e., some displayed ideal FA, but others did not), FA for each structure was standardized (difference between left and right values/left value, expressed as %) and compared between sites separately for each sex. We also created a composite FA index for each fish, based on an average of the standardized values for all five structures assessed.

4. Results

Across the three years, we examined 263 brook stickleback: 126 from Garvin Brook and 137 from SFWR (Table 1). Fish ranged in size from 26 to 69 mm TL, likely representing age 1 and age 2 fish (Becker 1983). All fish in both streams were sexually mature, based on examination of gonads. For all years combined, lengths of male and female fish did not differ within a stream, but both male and female fish were significantly larger in Garvin Brook than they were in SFWR (Figure 3A).
Overall condition of stickleback, based on Fulton’s K, did not differ between streams or between sexes when data from 2013 and 2014 were examined (Figure 3B), averaging 1.26 across streams and sexes. Although there were no significant changes in Fulton’s K with increasing fish length in either male or female fish from Garvin Brook (Figure 4A,B), female stickleback in SFWR displayed significantly declining condition with increasing fish length and male fish displayed a similar, nearly statistically significant decline (Figure 4C,D).
ANCOVA tests comparing internal measures of reproductive fitness and condition in stickleback (ovary and testes weights, GSI, HSI, egg number, egg weight) between the two study streams revealed several differences during the years examined (Table 2). Neither testes weights (Figure 5A) nor male GSI values (Figure 6A) differed significantly (Table 2) between the two streams after accounting for fish size. However, ovary weights (Figure 5B) and female GSI values (Figure 6B) displayed variable patterns in the different streams. Slopes of the total length-ovary weight regressions differed significantly between sites (Table 2), invalidating the ANCOVA test but demonstrating that larger females in Garvin Brook had larger ovaries than females of similar lengths in SFWR (Figure 5B). In contrast, ANCOVA revealed that Garvin Brook females had significantly greater GSI values than those from SFWR (Table 2, Figure 6B). ANCOVA also revealed significantly higher egg numbers in Garvin versus SFWR females across the entire range of fish size (Table 2, Figure 7A). Slopes of the total length-egg weight regressions differed significantly between sites (Table 2), again invalidating the ANCOVA test but demonstrating that larger females in Garvin Brook had larger eggs than females of similar lengths in SFWR (Figure 7B). Finally, ANCOVA for both male and female HSI (Table 2) indicated significantly higher HSI values (i.e., enlarged livers that represented a greater proportion of total body weight) in SWFR fish compared to Garvin fish (Figure 8A,B).
FA in Garvin Brook and SFWR stickleback, as assessed by left side-right side differences in head, jaw, eye, pectoral fin, and pelvic fin measurements, exhibited high variability during 2015 (Figure 9). Nine of the 20 stream/sex/structure categories examined did not display ideal FA, i.e., mean values differed significantly (P < 0.05) from 0 (Figure 9). Seven of these nine categories skewed positive (left-side structures more often longer than right-side structures), and six were observed in male stickleback. Due to these variations in ideal-versus-not ideal FA among the structures examined, we standardized all side-to-side differences for each individual fish prior to assessing possible differences in FA between streams.
Mean FA measures for all five structures examined in this study were higher for both male and female brook stickleback from SFWR than from Garvin Brook (Figure 10). These differences were significant for four of the five female characters (jaw length, head length, eye diameter, pectoral fin length) plus a composite index (Figure 10A). However, none of the male characters nor the male composite index displayed significant differences between Garvin Brook and SFWR (Figure 10B). Consequently, SFWR female stickleback displayed significantly higher degrees of FA than female fish from Garvin Brook, whereas male fish from both systems exhibited similar degrees of FA.

5. Discussion

Our study of brook stickleback condition and reproductive fitness in two headwater streams in Minnesota provides evidence that various indicators of stickleback health likely can be used to demonstrate fish exposure to impaired water quality. We found that stickleback in the more impaired SFWR displayed smaller size, declining Fulton’s condition with increasing fish size or age, and enlarged livers compared to stickleback from the less impaired Garvin Brook site. Female fish from SFWR also displayed smaller ovaries and smaller eggs, plus greater FA than females in Garvin Brook.
Brook stickleback morphology appears to respond to impaired water quality in several predictable ways. Previous studies [11,12,14,15] have reported that various stickleback species are sensitive both physiologically and morphologically to environmental perturbations. Typically, a relatively small number of freshwater aquatic species (e.g., fathead minnow Pimephales promelas, rainbow trout Oncorhynchus mykiss, cladocerans such as Daphnia magna, midges Chironomus tentans, scuds Hyallela azteca, zebra fish Danio rerio) are favored as models or bioindicators during traditional laboratory and/or field investigations used to assess the potential effects of various chemicals on aquatic life [25]. But other species, including brook stickleback, have been presented as either more sensitive bioindicators or more widely distributed and consequently more applicable over a broader geographic area [11,15,26] than many of the typical test species, especially under field conditions. Consequently, the results of our study add support to the recommendation [11] for using brook stickleback as a model for detecting environmental impacts on freshwater aquatic life.
Smaller size, reduced condition in larger (older) fish, and enlarged livers all suggest that brook stickleback in the SFWR were exposed to reduced water quality that produced some degree of chronic stress. Fish exposed to many pollutants often respond similarly to what we observed in stickleback in the SFWR [4,5,7,14]. Pollution stress can slow or suppress growth, leading to smaller overall body sizes for fish living in polluted waters [27,28,29,30,31]. Long-term exposure to poor water quality also may lead to reduced fish condition (e.g., lower Fulton’s K) in older fish, as low-level but continuing stress may reduce foraging ability/success, especially for older fish with higher maintenance needs [22,32,33]. The same impairments stressing fish also may suppress prey resources (i.e., fewer potential prey) [18] due to poor water quality, making it more difficult for stickleback to find and consume sufficient prey to remain healthy.
Both male and female sticklebacks in the SFWR had significantly higher HSI values compared to fish in Garvin Brook. Enlarged livers (hepatomegaly) often are indicative of contaminant stress in fishes [11,24,34,35,36]. Exposure to waterborne pollutants can induce greatly increased metabolic activity within the liver as it works to detoxify contaminants taken in through the gills or food, producing larger liver mass relative to body mass [24]. Since livers of female fish also enlarge during egg development due to the organ’s role in vitellogenesis [37], male fish (which do not undergo vitellogenesis except under rare conditions) provide better models for using relative liver size to assess exposure to water pollutants [24,38]. In our study, SFWR males averaged 67% higher HSI than Garvin males, whereas SFWR females averaged 30% higher values than Garvin females. These data suggest that fish in the SFWR likely had been exposed to more impaired water quality than those in Garvin Brook, and even vitellogenesis occurring in female fish did not mask the difference in liver size between the two systems.
The smaller eggs and ovaries present in female stickleback in the SFWR point toward some form of stress experienced by these fish during oogenesis. Egg development and growth in fishes is sensitive to many different environmental conditions, ranging from water temperature and oxygen levels to the presence of many different chemical pollutants [4,5,31]. Female fish are especially susceptible to estrogen-like chemicals that may either stimulate or suppress various reproductive processes [4,5,11,31] or interfere with sex determination [27]. The reduced GSI values for female stickleback in the SFWR parallel the declines observed in female stickleback GSIs after exposure to estradiol in the laboratory [11]. However, in some fish species the testes of males may be more impacted by some pollutants than the ovaries of the females, due primarily to the pollutants’ ability to alter various spermatogenesis processes [5]. It could be argued that smaller female fish in general produce fewer and smaller eggs than do larger females [39,40], complicating interpretation of the effects of water quality on stickleback egg number and size. However, our results demonstrated that female stickleback in the more impaired waters of the SFWR produced significantly fewer and smaller eggs than females of the same size in Garvin Brook, especially when the largest females were examined. Consequently, various reproductive characteristics of female brook stickleback appear to be sensitive to water quality, allowing for their use as bioindicators of pollutant exposure.
Elevated FA in female (but not male) brook stickleback in the SFWR demonstrates that some factor(s) impacted the early development of female embryos, interfering with the normal processes responsible for producing perfect bilateral symmetry. Many studies have reported that fish and amphibians exposed to various waterborne chemicals develop increased asymmetry, demonstrating the sensitivity of embryonic development to environmental stressors [7,9,14,15,41]. Several of the same anatomical structures measured in our study have been used when assessing FA in other fish species, and these structures have been shown to be impacted by poor water quality [7,9,14]. Increased FA in fish and other organisms generally is considered as a response to environmental stress, caused when a developing organism is forced to allocate more energy to stress response and has less energy to expend on developmental control [6,15].
FA has been criticized by some researchers as an inconsistent approach for assessing the presence of environmental stressors [6,15]. Criticisms include lack of repeatability, susceptibility of the procedure to sampling and measurement error, alternative causes of asymmetry not related to organismal stress, sample size issues, and others. FA methodologies and data analyses have been modified over time to compensate for several of these concerns [6,9,15,42]. Our study of brook stickleback used several of these modified approaches, and we were successful at demonstrating a consistent pattern of greater asymmetry in most of the characters assessed in female fish in SFWR (but not in male fish) relative to those in Garvin Brook.
Female fish in general may be more sensitive to environmental stressors, especially those water-quality impairments that are chemically similar to estrogen [4,5,11,27,31]. There also appear to be correlations between asymmetry and both egg number and size in female stickleback, with fish that exhibit asymmetry producing fewer and smaller eggs than those produced by symmetric individuals [14]. These correlations likely are not cause-and-effect, but instead are separate, organismal (physiological, developmental) responses to a common environmental stressor. In our study, we observed both reduced egg number/size and greater FA in female fish from the more impaired SFWR, although these were not assessed simultaneously in the same individual fish. However, these observations in female stickleback, coupled with limited to no responses (i.e., lack of FA, no change in testes weights or GSI values) in male fish to the poorer water quality in the SFWR, supports the contention that females are the more sensitive gender and have greater potential use as a bioindicator.

6. Conclusions

Our study of brook stickleback in two streams differing in water quality indicated that exposure of brook stickleback to reduced water quality from agricultural and urban runoff during early development resulted in increased morphological asymmetry in female (but not male) fish. In addition, continuing exposure to compromised water quality throughout life impacted both general condition and reproductive fitness of stickleback, especially older female fish. These results corroborate previous studies [11,14,15] which suggested that brook stickleback can serve as a bioindicator of water quality impairment, and that indicate the species can be used as an effective environmental monitor within its natural habitats across its wide geographic range in North America [11,15].

Author Contributions

NM and KH developed the concept for the paper. NM made all field collections, and NM and KH conducted the laboratory measurements. NM and KH analyzed the data and NM led the writing of the paper. Both authors have edited the manuscript and approve of the submitted version.

Funding

Partial funding for this project was provided by the Trout Unlimited Driftless Area Restoration Effort.

Institutional Review Board Statement

Fish collections were made under special permits (Numbers 18903, 19721, 20320) from the Minnesota Department of Natural Resources, Division of Fish and Wildlife, Section of Fisheries, and were conducted with the approval of the Winona State University Animal Care and Use Committee (1317064-1, 1310072-2). This research complied with all ethical standards.

Acknowledgments

We thank the City of Eyota and the Minnesota Department of Natural Resources, Division of Forestry, for access to the stream sites. We thank Jeff Hastings, Trout Unlimited, for supporting research on non-game animals in trout streams.

Conflicts of Interest

The authors declare no competing interests.

References

  1. Thronson, A.; Quigg, A. Fifty-five years of fish kills in coastal Texas. Estuar. Coasts 2008, 31, 802–813. [Google Scholar] [CrossRef]
  2. Phelps, N.B.D.; Bueno, I.; Poo-Muñoz, D.A.; Knowles, S.J.; Massarani, S.; Rettkowski, R.; Shen, L.; Rantala, H.; Phelps, P.L.F.; Escobar, L.E. Retrospective and predictive investigation of fish kill events. J. Aquat. Anim. Health 2019, 31, 61–70. [Google Scholar] [CrossRef]
  3. Schnaser, A.C.; Mundahl, N.D. Recovery of a headwater stream population of brown trout after a fish kill in southeastern Minnesota, USA. Environ. Biol. Fish. 2022, 105, 1179–1192. [Google Scholar] [CrossRef]
  4. Colborn, T.; vom Saal, F.S.; Soto, A.M. Developmental effects of endocrine-disrupting chemicals in wildlife and humans. Environ. Health Perspect. 1993, 101, 378–384. [Google Scholar] [CrossRef]
  5. Rohr, J.R.; McCoy, K.A. A qualitative meta-analysis reveals consistent effects of atrazine on freshwater fish and amphibians. Environ. Health Perspect. 2010, 118, 20–32. [Google Scholar] [CrossRef]
  6. Beasley, D.E.; Bonisoli-Alquati, A.; Mousseau, T.A. The use of fluctuating asymmetry as a measure of environmentally induced developmental instability: a meta-analysis. Ecol. Indic. 2013, 30, 218–226. [Google Scholar] [CrossRef]
  7. Simon, T.P.; Burskey, J.L. Deformity, erosion, lesion, and tumor occurrence, fluctuating asymmetry, and population parameters for bluntnose minnow (Pimephales notatus) as indicators of recovering water quality in a Great Lakes area of concern, USA. Arch. Environ. Contam. Toxicol. 2016, 70, 181–191. [Google Scholar] [CrossRef]
  8. McCoy, K.A.; Bortnick, L.J.; Campbell, C.M.; Hamlin, H.J.; Guillette, L.J. Jr.; St. Mary, C. Agriculture alters gonadal form and function in the toad Bufo marinus. Environ. Health Perspect. 2008, 116, 1526–1532. [Google Scholar] [CrossRef]
  9. Michaelsen, S.; Schaefer, J.; Peterson, M.S. Fluctuating asymmetry in Menidia beryllina before and after the 2010 Deepwater Horizon oil spill. PLOS ONE 2015. [Google Scholar] [CrossRef]
  10. Minnesota Pollution Control Agency. Minnesota’s impaired waters list. 2023. https://www.pca.state.mn.us/air-water-land-climate/minnesotas-impaired-waters-list.
  11. Muldoon, B.M.; Hogan, N.S. Biomarker responses to estrogen and androgen exposure in the brook stickleback (Culaea inconstans): a new bioindicator species for endocrine disrupting compounds. Comp. Biochem. Physiol., C Toxicol. Pharmacol. 2016, 180, 1–10. [Google Scholar] [CrossRef]
  12. Handy, R.D.; Runnalls, T.; Russell, P.M. Histopathologic biomarkers in three-spined sticklebacks, Gasterosteus aculeatus, from several rivers in southern England that meet the Freshwater Fisheries Directive. Ecotoxicology 2002, 11, 467–479. [Google Scholar] [CrossRef]
  13. Becker, G.C. Fishes of Wisconsin; University of Wisconsin Press: Madison, Wisconsin, USA, 1983; pp. 777–781. Available online: https://search.library.wisc.edu/digital/AWFIWUVZSK4EFH9B.
  14. Hecter, R.P.; Moodie, P.F.; Moodie, G.E.E. Pectoral fin asymmetry, dimorphism and fecundity in the brook stickleback, Culaea inconstans. Behaviour 2000, 137, 999–1009. [Google Scholar] [CrossRef]
  15. Lajus, D.L.; Golovin, P.V.; Yurtseva, A.O.; Ivanova, T.S.; Dorgham, A.S.; Ivanov, M.V. Fluctuating asymmetry as an indicator of stress and fitness in stickleback: a review of the literature and examination of cranial structures. Evol. Ecol. Res. 2019, 20, 83–106. [Google Scholar]
  16. Mundahl, N.D.; Hunt, A.M. Recovery of stream invertebrates after catastrophic flooding in southeastern Minnesota, USA. J. Freshw. Ecol. 2011, 26, 445–457. [Google Scholar] [CrossRef]
  17. Mundahl, N.D. 2022. Analysis of a catch-and-harvest fishing contest for brown trout spanning 50+ years: long-term trends influenced by fisheries management actions and angler behavior. N. Am. J. Fish. Manag, 2022; 42, 1653–1664. [Google Scholar]
  18. Mundahl, N.D.; Mundahl, E.D. Aquatic community structure and stream habitat in a karst agricultural landscape. Ecol. Process. 2022, 11, 18. [Google Scholar] [CrossRef]
  19. Mundahl, N.D.; Schnaser, A.C. Abundance, cover use, and clustering of brown trout spawning redds during stream habitat rehabilitation. Ecol. Freshw. Fish 2023, 32, 902–915. [Google Scholar] [CrossRef]
  20. Mundahl, N.D.; Varela, W.L.; Weaver, C.; Mundahl, E.D.; Cochran-Biederman, J.L. Stream habitats and aquatic communities in an agricultural watershed: changes related to a mandatory riparian buffer law. Environ. Manage. 2023, 72, 945–958. [Google Scholar] [CrossRef]
  21. Varela, W.L.; Mundahl, N.D.; Bergen, S.; Staples, D.F.; Cochran-Biederman, J.; Weaver, C.R. Physical and biological stream health in an agricultural watershed after 30+ years of targeted conservation practices. Water 2023, 15, 3475. [Google Scholar] [CrossRef]
  22. Ney, J.J. Practical use of biological statistics, pp. 167-191. In Kohler, C.C. and Hubert, W.A. (eds.) Inland fisheries management in North America, 2nd edition. 1999, American Fisheries Society, Bethesda, Maryland USA.
  23. Crim, L.W.; Glebe, B.D. Reproduction, pp. 529-553. In Schreck, C.B. and Moyle P.B. (eds.) Methods for fish biology. 1990, American Fisheries Society, Bethesda, Maryland USA.
  24. Popović, N.T.; Čižmek, L.; Babic, S.; Strunjak-Perović, I.; Čož-Rakovac, R. Fish liver damage related to the wastewater treatment plant effluents. Environ. Sci. Pollut. Res. 2023, 30, 48739–48768. [Google Scholar] [CrossRef]
  25. Rand, G.M.; Petrocelli, S.R. Fundamentals of aquatic toxicology: methods and applications; FMC Corporation: Princeton, New Jersey, USA, 1985. [Google Scholar]
  26. Gray, M.A.; Curry, R.A.; Arciszewski, T.J.; Munkittrick, K.R.; Brasfield, S.M. The biology and ecology of slimy sculpin: a recipe for effective environmental monitoring. FACETS 2018, 3. [Google Scholar] [CrossRef]
  27. Brown, A.R.; Owen, S.F; Peters, J.; Zhang, Y.; Soffker, M.; Paull, G.C.; Hosken, D.J.; Wahab, M.A.; Tyler, C.R. Climate change and pollution speed declines in zebrafish populations. PNAS, 2015; E1237–E1246. Available online: www.pnas.org/cgi/doi/10.1073/pnas.1416269112.
  28. Affandi, F.A.; Ishak, M.Y. Impacts of suspended sediment and metal pollution from mining activities oon riverine fish population – a review. Environ. Sci. Pollut. Res. 2019. [Google Scholar] [CrossRef]
  29. Deane, E.E.; and Woo, N.Y.S. Modulation of fish growth hormone levels by salinity, temperature, pollutants and aquaculture related stress: a review. Rev. Fish Biol. Fisheries 2009, 19, 97–120. [Google Scholar] [CrossRef]
  30. Kruitwagen, G.; Hecht, T.; Pratap, H.B.; Wendelaar Bonga, S.E. Changes in morphology and growth of the mudskipper (Periophthalmus argentilineatus) associated with coastal pollution. Mar. Biol. 2006, 149, 201–211. [Google Scholar] [CrossRef]
  31. Authman, M.M.N; Zaki, M.S.; Khallaf, E.A.; Abbas, H.H. Use of fish as bio-indicator of the effects of heavy metals pollution. J. Aquac. Res. Development 2015, 6, 4. [Google Scholar] [CrossRef]
  32. Sprague, J.B. Aquatic toxicology, pp. 491-528. In Schreck, C.B and Moyle P.B. (eds.) Methods for fish biology. 1990, American Fisheries Society, Bethesda, Maryland USA.
  33. Diana, J.S. Biology and ecology of fishes, 2nd edition. 2004, Biological Sciences Press, Cooper Publishing Group, Traverse City, Michigan USA. 498 pp.
  34. Slooff, W.; Van Kreijl, C.F.; Baars, A.J. Relative liver weights and xenobiotic-metabolizing enzymes of fish from polluted surface waters in the Netherlands. Aquat. Toxicol., 1983, 4, 1–14. [Google Scholar] [CrossRef]
  35. Fabacher, D.L.; Baumann, P.C. Enlarged livers and hepatic microsomal mixed-function oxidase components in tumor-bearing brown bullheads from a chemically contaminated river. Environ. Toxicol. Chem. 1985; 4, 703–710. [Google Scholar]
  36. Al-Ghais, S.M. Acetylcholinesterase, glutathione and hepatosomatic index as potential biomarkers of sewage pollution and depuration in fish. Mar. Pollut. Bull., 2013, 74, 183–186. [Google Scholar] [CrossRef]
  37. Long, S.M.; Tull, D.L.; De Souza, D.P.; Kouremenos, K.A.; Dayalan, S.; McConville, M.J.; Hassel, K.; Pettigrove, K.J.; Gagnon, M.M. Metabolomics provide sensitive insights into the impacts of low level environmental contamination on fish health – a pilot study. Metabolites 2020, 10, 24. [Google Scholar] [CrossRef]
  38. Duarte, R.M.; Sadauskas-Henrique, H.; de Almeida-Val, V.M.F.; Val, A.L.; Nice, H.E.; Gagnon, M.M. Biomarker responses and PAH ratios in fish inhabiting an estuarine urban waterway. Environ. Toxicol. 2017, 32, 2305–2315. [Google Scholar] [CrossRef]
  39. Moyle, P.; Cech Jr., J. Fishes: an introduction to ichthyology, 5th edition. 2003, Pearson/Prentice Hall, Upper Saddle River, New Jersey USA.
  40. Burton, D; Burton, M. Essential fish biology: diversity, structure and function. 2018, Oxford University Press, Oxford, United Kingdom.
  41. Hermita, Z.M.; Gorospe, J.G.; Torres, M.A.J.; Lumasag, G.J.; Demayo, C.G. Fluctuating asymmetry in the body shape of the mottled spinefoot fish, Siganus fuscescens (Houttuyn, 1782) collected from different bays in Mindanao Island, Philippines. Sci Int. (Lahore), 2013, 25, 857–861. [Google Scholar]
  42. Klingenberg, C.P.; McIntyre, G.S. Geometric morphometrics of developmental instability: analyzing patterns of fluctuating asymmetry with Procrustes methods. Evolution 1998, 52, 1363–1375. [Google Scholar] [CrossRef]
Figure 1. Aerial views of study site locations (indicated by stars) on the South Fork Whitewater River (SFWR) and Garvin Brook. The inset highlights the location of the study area (indicated by the star) within the southeastern corner of Minnesota, USA. The SFWR site was located within the City of Eyota, a small community surrounded by agricultural land. The Garvin Brook site was located within state forest lands.
Figure 1. Aerial views of study site locations (indicated by stars) on the South Fork Whitewater River (SFWR) and Garvin Brook. The inset highlights the location of the study area (indicated by the star) within the southeastern corner of Minnesota, USA. The SFWR site was located within the City of Eyota, a small community surrounded by agricultural land. The Garvin Brook site was located within state forest lands.
Preprints 88631 g001
Figure 2. Side view of a brook stickleback depicting the five morphological measurements (jaw length, head length, eye diameter, pectoral fin length, pelvic fin length) made on both left and right sides of each fish to assess fluctuating asymmetry.
Figure 2. Side view of a brook stickleback depicting the five morphological measurements (jaw length, head length, eye diameter, pectoral fin length, pelvic fin length) made on both left and right sides of each fish to assess fluctuating asymmetry.
Preprints 88631 g002
Figure 3. Box and whisker plots (means, medians, interquartile ranges, standard deviations, outliers) of total lengths (A, 2013-2015) and Fulton’s condition factors (B, 2013-2014) for male and female brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR). Total length bars not sharing a common letter beneath them are significantly different (ANOVA and least significant difference tests).
Figure 3. Box and whisker plots (means, medians, interquartile ranges, standard deviations, outliers) of total lengths (A, 2013-2015) and Fulton’s condition factors (B, 2013-2014) for male and female brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR). Total length bars not sharing a common letter beneath them are significantly different (ANOVA and least significant difference tests).
Preprints 88631 g003
Figure 4. Relationships between male (A, C) and female (B, D) brook stickleback total lengths and Fulton’s condition factors (K) from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream/sex combination.
Figure 4. Relationships between male (A, C) and female (B, D) brook stickleback total lengths and Fulton’s condition factors (K) from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream/sex combination.
Preprints 88631 g004
Figure 5. Total length-gonad weight relationships for male (A) and female (B) brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream/sex combination. Refer to Table 2 for results of ANCOVA comparisons between streams.
Figure 5. Total length-gonad weight relationships for male (A) and female (B) brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream/sex combination. Refer to Table 2 for results of ANCOVA comparisons between streams.
Preprints 88631 g005
Figure 6. Total length-gonadosomatic index (GSI) relationships for male (A) and female (B) brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream/sex combination. Refer to Table 2 for results of ANCOVA comparisons between streams.
Figure 6. Total length-gonadosomatic index (GSI) relationships for male (A) and female (B) brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream/sex combination. Refer to Table 2 for results of ANCOVA comparisons between streams.
Preprints 88631 g006
Figure 7. Total length-egg number (A) and total length-egg weight (B) relationships for female brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream. Refer to Table 2 for results of ANCOVA comparisons between streams.
Figure 7. Total length-egg number (A) and total length-egg weight (B) relationships for female brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream. Refer to Table 2 for results of ANCOVA comparisons between streams.
Preprints 88631 g007
Figure 8. Total length-hepatosomatic index (HSI) relationships for male (A) and female (B) brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream/sex combination. Refer to Table 2 for results of ANCOVA comparisons between streams.
Figure 8. Total length-hepatosomatic index (HSI) relationships for male (A) and female (B) brook stickleback from Garvin Brook and South Fork Whitewater River (SFWR), 2013-2014. Simple linear regression statistics are displayed for each stream/sex combination. Refer to Table 2 for results of ANCOVA comparisons between streams.
Preprints 88631 g008
Figure 9. Left side versus right side differences in measurements of five morphological structures (head length, jaw length, eye diameter, pectoral fin length, pelvic fin length) in male and female brook stickleback from Garvin Brook (A) and South Fork Whitewater River (B), 2015. Bars are means and vertical lines represent ± 1 standard deviation. Means highlighted by stars are skewed significantly different from 0 (i.e., do not display ideal fluctuating asymmetry).
Figure 9. Left side versus right side differences in measurements of five morphological structures (head length, jaw length, eye diameter, pectoral fin length, pelvic fin length) in male and female brook stickleback from Garvin Brook (A) and South Fork Whitewater River (B), 2015. Bars are means and vertical lines represent ± 1 standard deviation. Means highlighted by stars are skewed significantly different from 0 (i.e., do not display ideal fluctuating asymmetry).
Preprints 88631 g009
Figure 10. Standardized fluctuating asymmetry (difference between left side versus right side measurements, expressed as % of left side measurement) for five morphological structures (head length, jaw length, eye diameter, pectoral fin length, pelvic fin length) plus a composite index (see description in Methods) in female (A) and male (B) brook stickleback from Garvin Brook and South Fork Whitewater River, 2015. Bars are means and vertical lines represent ± 1 standard deviation. Between-stream comparisons highlighted by stars are significantly different.
Figure 10. Standardized fluctuating asymmetry (difference between left side versus right side measurements, expressed as % of left side measurement) for five morphological structures (head length, jaw length, eye diameter, pectoral fin length, pelvic fin length) plus a composite index (see description in Methods) in female (A) and male (B) brook stickleback from Garvin Brook and South Fork Whitewater River, 2015. Bars are means and vertical lines represent ± 1 standard deviation. Between-stream comparisons highlighted by stars are significantly different.
Preprints 88631 g010
Table 1. Numbers (N) and total lengths (TL, mm) of male (M) and female (F) brook stickleback collected from Garvin Brook and the South Fork Whitewater River (SFWR) for assessments of general and reproductive fitness, 2013-2015.
Table 1. Numbers (N) and total lengths (TL, mm) of male (M) and female (F) brook stickleback collected from Garvin Brook and the South Fork Whitewater River (SFWR) for assessments of general and reproductive fitness, 2013-2015.
Preprints 88631 i001
Table 2. Analysis of covariance (ANCOVA) and homogeneity of slopes test results comparing internal reproductive fitness and general condition variables in brook stickleback between Garvin Brook and South Fork Whitewater River sites, 2013 and 2014. Fish total length was used as a covariate. Significant (P < 0.05) tests are highlighted in bold font. GSI = gonadosomatic index, HSI = hepatosomatic index.
Table 2. Analysis of covariance (ANCOVA) and homogeneity of slopes test results comparing internal reproductive fitness and general condition variables in brook stickleback between Garvin Brook and South Fork Whitewater River sites, 2013 and 2014. Fish total length was used as a covariate. Significant (P < 0.05) tests are highlighted in bold font. GSI = gonadosomatic index, HSI = hepatosomatic index.
Preprints 88631 i002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated