ORIGINAL RESEARCH article

Front. Ecol. Evol., 29 September 2022

Sec. Ecophysiology

Volume 10 - 2022 | https://doi.org/10.3389/fevo.2022.1014814

Differential sensitivity of offspring from four species of goodeine freshwater fish to acute exposure to nitrates

  • 1. Programa Institucional de Doctorado en Ciencias Biológicas, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Mexico

  • 2. CONACYT-CIAD, AC Mazatlán Unit for Aquaculture and Environmental Management, Mazatlan, Mexico

  • 3. Laboratorio de Biología Acuática, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Mexico

  • 4. Laboratorio de Ecofisiología Animal, Instituto de Investigaciones Sobre Recursos Naturales, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Mexico

Abstract

Nitrate-nitrogen (NO3-N) pollution related to anthropogenic activities is increasing in freshwater ecosystems. Knowledge about NO3-N sensitivity in freshwater wild fish is needed to understand the differential tolerance between species. Goodeinae is a subfamily of 41 endemic fishes that inhabit central Mexico, with 33 species in the IUCN red list and three extinct. Distributional patterns suggest tolerant and sensitive goodeines related to the conservation gradient of freshwater ecosystems. Four species with a differential distribution and tolerance were selected to evaluate their physiological responses to NO3-N. Fish were exposed to different NO3-N concentrations for 96 h and the median lethal concentration (LC50) was determined. Swimming disorders plus gill and liver histopathological indexes were estimated and incorporated into an Integrated Biomarker Response (IBR) for each species. Skiffia lermae (LC50 = 474.332 mg/L) and Xenotoca variata (LC50 = 520.273 mg/L) were more sensitive than Goodea atripinnis (LC50 = 953.049 mg/L) and Alloophorus robustus (LC50 = 1537.13 mg/L). The typical histological damage produced by NaNO3-N exposure was fusion of secondary lamellae in gills. This was present in all species and cellular degeneration was observed at the highest concentrations. Secondary lamellae aneurysms were only observed in G. atripinnis. Liver alterations included vascular dilation in hepatic sinusoids, hyperemia and nuclear hypertrophy; higher concentrations produced hepatocyte cytoplasmic vacuolation and reduced frequency of cell nuclei. Behavioral and histopathological alterations could explain the differential species sensitivity. The results suggest that species which preserve gill function and transfer the task of detoxification to the liver might have the best chance of surviving in polluted environments. Moreover, species previously considered as tolerant may be highly susceptible to NaNO3-N exposure. Therefore, it is necessary to closely monitor NaNO3-N concentrations in freshwater ecosystems and, if possible, reduce their levels to avoid the loss of wild populations.

Introduction

One of the most representative fish groups in central Mexico are the goodeines. The major diversification within this group occurs in the Río Lerma-Santiago basin, one of the most polluted basins in the world (, ; ). The Goodeinae subfamily is composed of 41 species from 19 genera (; ). This group presents unique reproductive characteristics such as internal fertilization, matrotrophy, and viviparity (), and is considered one of the most at-risk groups in the world (). During the last decades, their populations have decreased related to pollution and other anthropogenic activities (). According to the , 13 species are critically endangered, 14 endangered, six vulnerable, two extinct in the wild, and one extinct. Most Goodeinae species are endemic or micro-endemic to a single spring, portion of a river or drainage, whereas others are widely distributed in the region (). These distributional patterns suggest tolerant and sensitive goodeines related to the conservation gradient of freshwater ecosystems (; ). For example: Skiffia lermae and Alloophorus robustus have been classified as endangered, sensitive, and vulnerable species, respectively; since they only are found in spring areas, lake shorelines and small tributaries (; ). Xenotoca variata and Goodea atripinnis are considered as least concern and tolerant species, because some areas they inhabit are classified as highly polluted aquatic environments (; ; ).

Nitrates (NO3) constitute the most stable and abundant form of dissolved inorganic nitrogen in aquatic ecosystems (). Although nitrate-nitrogen (NO3-N) may be naturally present, they accumulate in the ecosystem due to human activities such as livestock and agricultural use of inorganic nitrogenous fertilizers. Municipal and industrial wastewater effluents also increase nitrogenous compounds (; ). Despite the abundance of NO3-N, their impact on fish welfare has been underestimated since they are less toxic than other nitrogenous compounds (; ; ). However, NO3-N negatively impact fish health (; ; ; , , ; ).

The standard sensitivity evaluation to acute toxicity is the median lethal concentration (LC50). This index is usually complemented with other responses to increase the understanding of physiological trade-offs under acute exposure or predict the long-term permanence of fish in polluted environments (; ). Exposure bioassays to NO3-N in fish revealed more sensitivity during early life stages (), likely due to interference with biological processes such as body growth and reproductive development or behavior (). The gills are the main route of NO3-N absorption (). Physiological responses to pollutants aim to maintain gill activity preceding the loss of function related to the degree of damage (). Inside the body, NO3-N causes physiological and behavioral alterations (; ; ). The liver is the detoxifying organ and its condition indicates environmental quality (; ; ). Toxic pollutants cause functional alterations of the gills and liver, such as gas exchange failure and metabolic damage. These changes are associated with altered swimming behavior, such as loss of balance, top or bottom position in the water column and lethargy (; ; ; ).

Given the conservation relevance of goodeines and based on the premise that differential persistence in freshwater bodies could be related to differential sensibility to environmental pollutants, the main objective of this study was to evaluate NO3-N sensitivity in offspring of S. lermae, X. variata, G. atripinnis, and A. robustus. Differential sensitivity was evaluated by LC50 values, histopathological indexes from the gills and liver, as well as swimming behavior alterations under acute exposure (96 h). Integrated Biomarker Response (IBR) analysis was performed to contrast their sensitivity.

Materials and methods

Ethical statements

Sampling and laboratory fish handling protocols were reviewed and approved (SEMARNATSGPA/DGVS/00012/19, PPF/DGOPA-014/20). The protocols also followed the . Toxicity was evaluated with the Fish Acute Toxicity Test ().

Field collection and fish maintenance

Gravid females of S. lermae (n = 25), X. variata (n = 20), G. atripinnis (n = 7), and A. robustus (n = 5) were collected from Zacapu Lake, central Mexico (19°49′26.88″N; 101°46′37.32″W), during years 2020 and 2021, using 2 mm aluminum mesh minnow traps (Gee-minnow-traps® G-40, USA). Gravid females were transported to the laboratory with water from the collection site. Females were acclimatized for 15 days in 120 L aquariums (loading capacity: 6 L per fish), with dechlorinated water and 7 mg/L of dissolved oxygen, at 22°C, under a natural light photoperiod. They were fed two times each day with Artemia sp. and commercial Wardley fish flakes (). On post-natal days 1–7, offspring were separated from their mothers, fed and housed under the same conditions previously described. Offspring remain in these conditions until post-natal 21, when they were subjected to the experiment (see below).

Water quality monitoring and solutions

Temperature, dissolved oxygen, pH, salinity, total dissolved solids and electric conductivity were recorded every 24 h using a multiparameter probe (YSI EXO2. Ohio, USA). The toxicity assay solution was prepared with sodium nitrate (NaNO3, 97%, Sigma-Aldrich) according to a stoichiometric calculation (). Concentrations of ammonium-nitrogen (NH4-N), nitrite-nitrogen (NO2-N), and NO3-N were measured following the methodologies described in the American Public Health Association ().

Toxicity test

Each fry was acclimatized for 7 days, weighted with an analytical precision balance (OHAUS™ Adventurer d = 0.0001 g, China), their length was measured using a digital Vernier caliper (Thermo Fisher Scientific™ S/N 1366162, USA) and then randomly assigned to a treatment. Three replicates were used per treatment, with 10 fish in 10 L tanks (loading capacity: 1 L per fish). A control group was included for each species.

During NO3-N exposure, survival was monitored at 0.5, 1, 2, 3, 4, 6, 9, 12, 24, 36, 48, 72, and 96 h. Mortality was confirmed by the absence of response to stimulation. Swimming behavior (loss of balance and position in the bottom, middle, or surface of the water column) was recorded during 10 min by counting the animals that exhibited these actions at the same intervals mentioned above. Changes in swimming behavior were categorized as “low effect,” indicating an alteration in 10–33% of individuals; “moderate effect” 34–63% and “high effect,” 64–100%. Changes to swimming patterns were recorded only if they occurred in 10% of the fish within each aquarium (). At the end of the experiment, surviving fish were anesthetized using clove oil (50 μl L-1) and processed for histology. Fish in the highest NO3-N concentrations were removed from the aquariums immediately after they became immobile.

Histological evaluation

Three individuals from three treatments (control, that closest to the LC50 and the most lethal concentration, LC100) per species were randomly selected and processed by routine histology. Gill and liver 5 μm transversal sections were obtained and stained with hematoxylin-eosin (; ). Slides were observed with a Leica microscope (DM3000, Germany), photo-documented using a camera (Leica DFC310 FX, Germany) and analyzed using NIH ImageJ software.

Gill and liver histopathological indexes

The histopathological condition index (HI) was calculated for each organism (; ). Briefly, according to distinctive histological features, gill and liver alterations were classified into circulatory, progressive and regressive patterns. A pathological importance factor (w) was assigned for each type of damage (Supplementary Table 1). Similarly, a score value (a) was assigned for the occurrence of alterations from 0 to 5 (where 0 is no change; 1, changes in 1–20% of the sections; 2, 21–40%; 3, 41–60%; 4, 61–80%; 5, 81–100%, respectively). The HI for each rp was determined using the formula: HIrp org = Σalt (w × a) where HIrp org is the reaction index of the organ for each reaction pattern (rp) and alt is the alteration ().

Liver damage tissue index

For the liver, tissue damage was assessed by nuclear density (nD = number of nuclei/mm2) and nuclear area (nA = area of hepatocyte nuclei in μm2; ). The average value of controls (nDc and nAc) was used as a reference to calculate a relative value for each feature. Since an increase in the nuclear area could indicate hypertrophy and a reduction could be related to pyknotic nuclei, the absolute value of the difference to 1 multiplied by 0.5 was used. Finally, the values were integrated into the following index: LDTI = (nD*1/nDc) + [| 1 – (nA*1/nAc)| × 0.5], where values near 1 indicate a normal or reference condition and values closest to 0 indicate major tissue damage.

Statistical analysis and integrated biomarkers response

All analyses were performed with R software (, version 4.1.3). Variables were tested for normality and homoscedasticity with Shapiro-Wilk and Levene’s tests, respectively. The survival rate for each species was calculated by the Kaplan Meier method and compared with the log-rank test using the “survival” and “survminer” packages. Acute toxicity to NaNO3-N was determined by the LC50 for each species (95% confidence intervals after 96 h) using the Spearman-Karber method and considered significantly different when not overlapping using the “tsk” package. Relationships between swimming disorders and concentrations at each exposure time were analyzed with Spearman correlations using the “stats” package; graphs were produced using the “ggpubr” package. Histomorphometric analyses were performed using the “car” and “pgirmess” packages. LDTI data distribution was analyzed by a Kruskal-Wallis test. Tukey multiple range post-hoc tests were used to discriminate means (). The values obtained for each parameter were expressed as mean ± standard error of the mean. Survival rates, LC50, swimming behavior, GHI, LHI, and LDTI were incorporated into an IBR to visualize and compare species-specific responses to acute NaNO3-N exposure using the “IBRtools” package ().

Results

NaNO3-N exposure and water quality

Concentrations for acute toxicity assays were determined for each species. Thus, each species was exposed to slightly different ranges; S. lermae: 100, 200, 400, 800, and 1,200 mg NO3-N/L; X. variata: 200, 400, 800, 1,200, and 1,600 mg NO3-N/L; G. atripinnis: 400, 800, 1,200, 1,600, 2,000, and 2,400 mg NO3-N/L; and A. robustus: 400, 800, 1,200, 1,600, 2,000, 2,400, and 2,800 mg NO3-N/L. The body mass, length and total number of offspring for each species was: 0.028 ± 0.0004 g; 12.265 ± 0.058 mm; n = 180 for S. lermae, 0.072 ± 0.001 g; 15.656 ± 0.103 m; n = 180 for X. variata, 0.184 ± 0.001 g; 20.122 ± 0.047 mm; n = 210 for G. atripinnis and 0.182 ± 0.001 g; 21.382 ± 0.055 mm; n = 240 for A. robustus.

Water quality conditions (temperature: 22.059 ± 0.003°C; dissolved oxygen: 7.069 ± 0.001 mg/L; pH: 8.930 ± 0.004 log units; ammonium: 0.003 ± 4.88104E-06 mg/L and nitrite: 0.042 ± 0.002 mg/L) in the experimental tanks were similar for all groups and maintained in a range similar to the natural variations observed in Zacapu lake (). Salinity (0.281 ± 0.001–12.041 ± 0.096 g/L), total dissolved solids (379.033 ± 1.913–13,017.266 ± 84.070 mg/L) and electric conductivity (551.290 ± 2.691–19,103.63 ± 141.191 μS/cm) increased as NaNO3-N concentrations increased.

Survival rate and lethal concentration 50%

Survival rates for the four Goodeinae species decreased as NaNO3-N concentrations increased and were related to exposure time. No mortality was observed at control concentrations (Figure 1). Survival rates for S. lermae and X. variata showed significant differences above 200 mg NaNO3-N/L vs. control [S. lermae: X2(5) = 306; p ≤ 0.001; X. variata: X2(5) = 149; p ≤ 0.001]; no fish survived after 96 h at 1,200 mg NaNO3-N/L for S. lermae and at 1,600 mg NaNO3-N/L for X. variata (Figures 1A,B). Survival rates for G. atripinnis showed significant differences above 800 mg NaNO3-N/L vs. control [X2(6) = 131; p ≤ 0.001] and no fish survived at 2,400 mg NaNO3-N/L (Figure 1C). For A. robustus significant differences were observed at 2,000 mg NaNO3-N/L vs. control [X2(7) = 940; p ≤ 0.001] and all fish died at 2,800 mg NaNO3-N/L (Figure 1D).

FIGURE 1

The calculated LC50 values and corresponding 95% confidence intervals for each species followed the survival trend (Figure 2). Skiffia lermae was the most sensitive, LC50 = 474.332 (373.850–601.820) mg NaNO3-N/L, followed by X. variata LC50 = 520.273 (380.390–711.594) mg NaNO3-N/L, G. atripinnis LC50 = 953.049 (722.835–1256.584) mg NaNO3-N/L, and A. robustus LC50 = 1537.130 (1203.223–1963.699) mg NaNO3-N/L.

FIGURE 2

Swimming behavior

Swimming disorders were related to exposure time and concentration (Tables 1, 2 and Supplementary Figure 1). The first alteration in S. lermae offspring was that they remained on the bottom of the water column at 24 h of exposure to the lowest concentration (100 mg NaNO3-N/L). The other three species showed surface permanence as the first alteration at 48 h of exposure to 200 mg NaNO3-N/L in X. variata; while for A. robustus and G. atripinnis this occurred at 36 and 84 h of exposure to 400 mg/L, respectively (Table 1). The prevailing behavior during acute exposure was permanence at the bottom without exploring the aquarium, but some fish also showed surface preference.

TABLE 1

SpeciesTime (h)Concentration (mg/L)
S. lermae01002004008001,200
3/+
4/+
5/+
6/+
12/+++
24++/++/++++++
36++/++++++++
48/++/+++++++++++
60/++/++++++++++++
72/+++++++++++++++
84+++++++++++++++
96++++++++++++x
X. variata02004008001,2001,600
12+
24+/+++++
36+++++++
48//+++++++++
60/++/+++++++++++
72/++/+++++++++++
84/++++++++++++x
96/++++++++++++x
G. atripinnis04008001,2001,6002,0002,400
12+
24/+/+/+++/++++++
36/+/++/+++/++++++
48/+/++++++++++++
60/+/++++++++++++
72/+/++++++++++++
84//++/++++++++++++
96/+/+++++++++++x
A. robustus04008001,2001,6002,0002,4002,800
5/
6/+
12//++
24/++++++++++
36///+/++++++++++
48//++/+++/++++++++++++
60/+/++/+++/++++++++++++
72/+++/+++++++++++++++x
84/++++++++++++++++++
96/++++++++++++++++++

Aquarium distribution.

Swimming behavior was recorded every hour during the six first hours; however, only the times when unusual behaviors were observed are shown, n = 30. Normal behavior in all fish (−); low affectation, unusual behavior was observed in 10–33% of the population (/top, + bottom); moderate affectation, 34–63% (//top, ++bottom); high affectation, 64–100% (///top, +++bottom); all individuals dead (x).

TABLE 2

SpeciesTime (h)Concentration (mg/L)
S. lermae01002004008001,200
24++++
36++++++
48++++++++
60+++++++++
72++++++++++
84++++++++++
96+++++++x
X. variata02004008001,2001,600
24+++++
36+++++++
48++++++++++
60++++++++++++
72++++++++++++
84++++++++++++++
96+++++++++++x
G. atripinnis04008001,2001,6002,0002,400
24+++++++++
36++++++++++++
48+++++++++++++
60++++++++++++++
72++++++++++++++
84++++++++++++++
96+++++++++++x
A. robustus04008001,2001,6002,0002,4002,800
24+++++
36++++++
48++++++++++++++
60+++++++++++++++
72++++++++++++++++x
84+++++++++++++++++
96++++++++++++++++++

Loss of equilibrium.

Swimming behavior was recorded every hour during the six first hours; however, only the times when unusual behaviors were observed are shown, n = 30. Normal swimming in all fish (−); low affectation, loss of equilibrium was observed in 10–33% of the population (+); moderate affectation, 34–63% of fishes showed loss of equilibrium and some swam laterally (++); high affectation, 64–100% of fishes swam laterally (+++); all individuals dead (x).

Loss of balance was related to exposure time and showed a similar trend, with S. lermae being more sensitive, followed by other species (Table 2 and Supplementary Figure 1). Goodea atripinnis did not show a loss of equilibrium at the lowest concentration.

Histological evaluation

Gill histopathological lesions integrated in the GHI showed a clear concentration-response pattern in all species (Table 3 and Figure 3). Circulatory index lesions were in a similar range; however, G. atripinnis was the only species that showed secondary lamellae aneurysms at the value closest to the LC50 (Figure 3). Lamellae aneurysms constitute the maximal circulatory damage associated with toxic exposure (w = 3 Supplementary Table 1). Secondary lamellae fusion was typical in S. lermae (Figure 3), where the progressive index reached the highest values (w = 3 Supplementary Table 1). Cellular degeneration was observed in the gills of all fish exposed to the lethal concentration. Thus, the regressive index showed the highest values at these concentrations for each species due to branchial degeneration resulting in death (w = 3 Supplementary Table 1).

TABLE 3

SpeciesNaNO3-N conc. (mg/L)Circulatory indexProgressive indexRegressive indexGHICirculatory indexProgressive indexRegressive indexLHINuclear area (μ m2)No. of cell nuclei/mm2LDTI
S. lermaeControl2.400 ± 0.3052.000 ± 0.7210.000 ± 0.0004.4000.333 ± 0.0660.000 ± 0.0000.000 ± 0.0000.33313.454 ± 0.24710319.783 ± 108.4771.092a
400 mg/L2.133 ± 0.29013.266 ± 1.9870.466 ± 0.13315.8663.133 ± 0.2902.00 ± 0.2302.866 ± 0.4808.00018.785 ± 0.2956834.688 ± 63.4880.864c
1,200 mg/L1.200 ± 0.1156.066 ± 1.57611.200 ± 1.10118.4663.666 ± 0.1330.333 ± 0.13316.600 ± 0.90120.6008.919 ± 0.230922.927 ± 22.6450.272f
X. variataControl3.333 ± 0.3330.733 ± 0.1760.000 ± 0.0004.0660.266 ± 0.0660.000 ± 0.0000.333 ± 0.0660.60013.179 ± 0.1309975.610 ± 104.5091.048ab
400 mg/L4.200 ± 0.4616.600 ± 0.8081.066 ± 0.40511.8661.400 ± 0.1150.733 ± 0.2907.400 ± 0.6119.53316.724 ± 0.2826325.203 ± 65.8540.788d
1,600 mg/L1.533 ± 0.2404.133 ± 1.63412.333 ± 0.94018.0002.333 ± 0.4800.333 ± 0.06614.266 ± 3.09016.93314.273 ± 0.3082128.211 ± 123.5680.316f
G. atripinnisControl2.666 ± 0.4371.200 ± 0.4000.000 ± 0.0003.8660.466 ± 0.0660.000 ± 0.0000.000 ± 0.000.46612.433 ± 0.17011761.518 ± 97.9381.067a
800 mg/L5.333 ± 1.5502.800 ± 0.8081.133 ± 0.3339.2661.866 ± 0.2402.266 ± 0.1762.333 ± 0.1766.46613.398 ± 0.18710883.469 ± 114.4830.999b
2,400 mg/L2.800 ± 0.72110.133 ± 0.4377.933 ± 1.8720.8662.600 ± 0.2300.333 ± 0.1338.866 ± 1.96411.80010.736 ± 0.2753989.807 ± 290.3730.458e
A. robustusControl2.600 ± 0.3050.533 ± 0.1330.000 ± 0.0003.1330.266 ± 0.0660.000 ± 0.0000.400 ± 0.1150.66611.816 ± 0.1598119.241 ± 117.6431.064a
1,600 mg/L2.533 ± 0.2902.466 ± 0.8961.533 ± 0.1766.5330.933 ± 0.1331.466 ± 0.1335.000 ± 1.5147.40016.470 ± 0.2045349.593 ± 86.6090.858c
2,800 mg/L1.866 ± 0.2403.933 ± 0.76811.466 ± 0.83517.2661.466 ± 0.4661.00 ± 0.2309.000 ± 0.98611.46612.501 ± 0.2053371.274 ± 36.2250.501e

Gill (GHI) liver histopathological (LHI) and tissue damage indexes (LTDI) of fish exposed to NaNO3-N (control, concentration closest to the LC50 and LC100), n = 3.

Superscript letters indicate statistical differences (Kruskal-Wallis-test with Tukey’s range post-hoc test, p ≤ 0.05).

FIGURE 3

In the liver, a concentration-response relationship for LHI also occurred (Table 3 and Figure 4). Vascular dilation in hepatic sinusoids occurred mainly in S. lermae (w = 1 Supplementary Table 1), while other species showed hyperemia as circulatory liver damage (Figure 4 and w = 3 Supplementary Table 1). Nuclear hypertrophy, as quantified in the progressive index (w = 1 Supplementary Table 1), was frequent at those concentrations close to the LC50 in all species. Hepatocyte cytoplasmic vacuolation and a relative reduction in nuclear frequency was observed in all species, especially in fish exposed to the lethal concentration, where the regressive index showed higher values than in other treatments (w = 3 Supplementary Table 1).

FIGURE 4

Quantitative hepatocyte density and nuclear area analysis were performed for the LTDI (values close to 0 indicate relevant tissue damage), which showed a concentration-response pattern (Table 3). The lowest hepatocyte density was observed at the lethal concentration, while hepatocyte nuclear area was higher at concentration closest to the LC50 in all species. However, the nuclear area in G. atripinnis showed similar values to their respective control group (Table 3).

Integrated biomarker response

For comprehensive evaluation, all traits analyzed were incorporated into an IBR for each species (Figure 5 and Supplementary Table 2). IBR values for each species were well weighted; control fish showed lower values and the lethal concentration had the highest score in all cases. Intercomparisons at doses closest to the LC50 indicated that the most sensitive species was S. lermae (IBR = 2.208), with the highest value, followed by X. variata (IBR = 2.124) at 400 mg NaNO3-N/L. A. robustus showed a higher IBR value than G. atripinnis; these were the most tolerant species to NaNO3. However, they differed by twice the concentration, in G. atripinnis it was 800 mg NaNO3-N/L (IBR = 0.977) and in A. robustus 1,600 mg NaNO3-N/L (IBR = 1.560). An important factor for the IBR in tolerant species was the liver (LHI), contrary to the most sensitive species where the important factor was gill affectation (GHI); swimming behavior was a promising biomarker in all species.

FIGURE 5

Discussion

There was a clear concentration-response relationship, also dependent on time, in all exposed fish. This was evidenced by the Kaplan-Meier survival curves and their differential sensitivity. Low concentrations produced mortality in S. lermae and X. variata, while G. atripinnis and A. robustus were more tolerant to higher NaNO3-N concentrations. The level of affectation, due to NO3-N exposure, between goodeine species suggests differential tolerance. Reference parameters for NaNO3-N acute toxicity in freshwater fish offspring are consistent with the results of G. atripinnis and A. robustus (Supplementary Table 3). However, the values for S. lermae and X. variata suggest they are the most sensitive species currently reported.

Sensitivity, ontogenetic stage, and body size

Fish responses to NO3-N depend on adsorption (accumulation time in plasma), distribution, metabolism, and excretion (). A relationship between sensitivity and ontogenetic stage has been found. suggested that offspring were more tolerant to elevated NO3-N than adults, contrary to other studies where earlier life stages were more sensitive (). Also, species with a bigger body size are less sensitive to NO3-N (). In fact, a recent study in relatively larger species found that low NO3-N concentrations did not alter growth and reproduction ().

In this study, the four Goodeinae species had the same age during NaNO3-N exposure but differed in body size and developmental stage. Skiffia lermae was the smallest species, followed by X. variata, while G. atripinnis and A. robustus were the largest species. Likewise, S. lermae have earlier maturation and begin their reproductive life approximately 3 months after birth (body length: 29.5 ± 5.7 mm in females and 25.2 ± 5.4 mm in males. ). In contrast, G. atripinnis start its reproductive stage at 8 months (), with a body length of 65.8 ± 17.6 mm in females and 65.6 ± 20.7 mm in males (). Thus, the smallest and most precocious species were the most sensitive, while more tolerant responses were found in larger species.

Swimming behavior reflects differential sensitivity to nitrate-nitrogen

Swimming behavior alterations have been considered as an attempt to face unfavorable environmental conditions, from the regular swimming pattern that occurs in the middle of the water column (). The behavioral patterns observed in response to acute NaNO3-N were: location on the surface or bottom of the water column, loss of balance and lethargy. Skiffia lermae showed swimming behavior alterations earlier and at lower NaNO3-N concentrations, followed by X. variata, A. robustus, and G. atripinnis. Permanence in the surface of the water column is a compensatory response to inefficient oxygen uptake caused by damage to the gills or by conversion of oxygen-carrying pigments (i.e., hemoglobin) to forms that are incapable of carrying oxygen (i.e., methemoglobin ). Loss of balance and position at the bottom of the water column have been associated with a decreased capacity for respiratory gas exchange, as well as with the depletion of liver lipid reserves (). Other studies have reported a loss of balance and lethargy in fish exposed to NO3-N (; ; ; ).

Histopathological alterations of gills and liver could reflect different sensitivity to nitrate-nitrogen in goodeine species

Remarkably, the LC50 value for A. robustus is similar to C. clupeaformis and O. mykiss, two euryhaline species having proficient ion exchange regulation to tolerate changes between freshwater and marine environments (; ). This is interesting given that the gills are one of the main organs involved in ion regulation and ammonia excretion, where specific glycoproteins are involved (). At the histological level, A. robustus showed preserved gill structure after NaNO3-N exposure. This suggests that gill function is maintained due to regressive, progressive and circulatory alterations restricted to the apical region. This was contrary to S. lermae which showed significant damage (highest GHI at minor concentrations), mainly lamellar fusion and cell hyperplasia; or X. variata and G. atripinnis which showed moderate gill alterations, like edema and aneurysms. Gill histopathological alterations in response to water conditions and pollutants, including NO2-N, have been described (). The gills play a critical physiological role in nitrogenous waste detoxification; however, the response to NO3-N has not been fully explained (). The branchial epithelium responds to non-optimal environmental changes through circulatory and progressive alterations oriented to preserve gill function. If ecological conditions remain unfavorable, degenerative mechanisms like lamellar fusion are triggered leading to loss of gill function (; ). Thus, lamellar fusion observed in S. lermae fry could lead to ion-exchange failure and consequent fish anoxia. The different degrees of gill tissue damage observed in X. variata, G. atripinnis, and A. robustus suggest different mechanisms to face NO3-N exposure.

Liver histopathology also confirmed differential tolerance due to variation in circulatory, progressive and regressive alterations. These constitute the LHI, which showed similar values between species. Xenotoca variata and A. robustus showed the highest regressive index associated with hepatocyte cytoplasmic vacuolation, pyknosis, and a relative decrease in nuclear density. Meanwhile, S. lermae and G. atripinnis showed circulatory and progressive alterations, respectively. The liver has an essential role in metabolic pathways and pollutant detoxification (). Exposure to environmental pollutants could result in erythrocyte accumulation in sinusoids (hyperemia), vascular dilation and hemorrhage, as well as hepatocyte hypertrophy. The degenerative indicators were hepatocyte cytoplasmic vacuolation, pyknosis, and decreased nuclear density. An increased size of hepatocyte nuclei (nuclear hypertrophy) reveals higher transcriptional and metabolic activity in response to pollutants. In contrast, a decreased cell density is related to liver degeneration (). Here the hepatocyte nuclear area and number of cell nuclei were quantified and integrated into the LDTI. All species showed an increased cell nuclear area and a decreased number of cell nuclei. However, X. variata presented the highest liver damage, followed by A. robustus and S. lermae, whereas G. atripinnis showed less liver damage.

Combined histopathological evaluation suggested that gill and liver damage were associated with alteration of critical metabolic processes, such as inhibition of protein synthesis, depletion of glycogen reserves and reduced ability to detoxify environmental pollutants (). In this regard, S. lermae showed the most severe gill damage and a relatively well-preserved liver structure. On the contrary, X. variata showed a high GHI and the highest LHI, whereas G. atripinnis and A. robustus showed lower GHI and LHI. Interestingly, G. atripinnis showed a higher GHI than A. robustus but a lower LHI. This supports the hypothesis that species that preserve gill structure and function, letting the liver carry out detoxification, likely have better survival opportunities. It has been demonstrated that lamellar fusion in response to acute NO3-N exposure limits the absorption of toxic substances and inhibits ion exchange ().

The inclusion of all evaluated responses into IBR radar graphs allowed identifying the GHI, LHI, and swimming behavior as the main factors that contribute to explain the differential sensitivity of Goodeinae species to NO3-N. Sensitive species (S. lermae and X. variata) showed similar IBR values at the concentration closest to the LC50, the first species showed the most severe gill damage. A certain tolerance level to NO3-N exposure was found in G. atripinnis and A. robustus species which could be related to the activation of physiological detoxifying strategies for survival. In both tolerant species, gill damage was minor, suggesting that fish species showing gill NO3-N tolerance have higher survival opportunities in polluted freshwater bodies.

The results show that NO3-N are toxic since they produce behavioral, branchial, and hepatic alterations at a lower concentration in goodeines. Thus, it is recommended that their concentrations be monitored in freshwater bodies where goodeines inhabit. Their toxic effect may be increased by alterations in other water quality parameters like pH, hardness, and temperature. In this study, most water quality variables were close to Zacapu lake parameters; except salinity, TDS and conductivity, which are directly associated with the sodium ions in solution. Due to the large difference in progeny between species (S. lermae and X. variata females produce a mean of 12 and 38 hatchlings, respectively; whereas G. atripinnis and A. robustus females produce larger progeny of up to 150 offspring. ), this study evaluated offspring from a different number of gravid females. These differences in fecundity between goodeine species, in addition to differential sensitivity to NaNO3-N, could increase the ecological risk for S. lermae and X. variata in wild freshwater ecosystems. However, population comparisons were outside the scope of this study.

Conclusion

Goodeinae species show a differential NaNO3-N sensitivity. Since the most sensitive species showed higher gill damage, differential sensitivity could be related to different physiological mechanisms to face NO3-N. Thus, species that preserve gill function and transfer the task of detoxification to the liver might have the best chance of surviving in polluted environments. Future studies should evaluate ecologically relevant NO3-N concentrations to confidently predict goodeines persistence in polluted freshwater bodies.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The animal study was reviewed and approved by the Mexican Ministry of Environmental and Natural Resources (SEMARNATSGPA/DGVS/00012/19) and the Secretariat of Agriculture and Rural Development (SAGARPA: PPF/DGOPA-014/20).

Author contributions

IV-V: investigation, methodology, formal analysis, and writing—original draft. BY-R: formal analysis and writing—review and editing. RR-J: conceptualization and methodology. MH-V and RH-M: methodology. OD-D: funding acquisition, conceptualization, and writing—review and editing. EM-H: funding acquisition, data curation, writing—original draft, and writing—review and editing. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Coordinación de la Investigación Científica-UMSNH to EM-H and OD-D and Chester Zoo Garden, Association Beauval Nature Pour la Conservation et la Recherche, Willema Zoo, Ostrava Zoo, The Mohammed Bin Zayed Species Conservation Found, Haus des Meeres—Aqua Terra Zoo, Poecilia Scandinavia, zoologische gesellschaft für arten- und populationsschutz, Goodeid Working Group, American Livebearers Association to OD-D. IV-V was Ph.D. fellow from CONACYT (grant no. 743476).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2022.1014814/full#supplementary-material

References

Summary

Keywords

ecotoxicology, LC50 96 h, histopathological changes, swimming behavior, gill, liver

Citation

Villa-Villaseñor IM, Yáñez-Rivera B, Rueda-Jasso RA, Herrera-Vargas MA, Hernández-Morales R, Meléndez-Herrera E and Domínguez-Domínguez O (2022) Differential sensitivity of offspring from four species of goodeine freshwater fish to acute exposure to nitrates. Front. Ecol. Evol. 10:1014814. doi: 10.3389/fevo.2022.1014814

Received

09 August 2022

Accepted

12 September 2022

Published

29 September 2022

Volume

10 - 2022

Edited by

Georgina Alexandra Rivera-Ingraham, Universidad Católica del Norte, Chile

Reviewed by

Wenyi Zhang, Nanjing Normal University, China; Ahasan Habib, University of Malaysia Terengganu, Malaysia

Updates

Copyright

*Correspondence: Esperanza Meléndez-Herrera,

This article was submitted to Ecophysiology, a section of the journal Frontiers in Ecology and Evolution

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics