1. Introduction
Heavy metals are a major concern for terrestrial and aquatic ecosystems because of their negative effects even in very low concentrations, their nondegradable nature, their high capacity for bioaccumulating, and their long-term persistence making them one of the main environmental problems of the 21st century worldwide [
1,
2,
3]. Some heavy metals have been classified in terms of biological function as beneficial or essential for living organisms based on their biological functions while others, such as lead (Pb), cadmium (Cd), mercury (Hg), and arsenic (As), are considered to be non-threshold micropollutants able to induce toxic effects in living organisms [
4,
5,
6,
7].
Pb is a highly toxic and non-biodegradable heavy metal listed in the Priority List of Hazardous Substances released by the Agency for Toxic Substances and Disease Registry (ATSDR) [
8,
9]. Although Pb is naturally occurring in the environment, its unregulated use in a number of human activities (e.g., mining and agriculture activities, paint pigments production) has resulted in an increase in its levels in all environmental comparts, raising great concern for both human and wildlife [
10]. Lead pollution of the aquatic environment occurs through agricultural, domestic, and industrial wastewater discharges, exerting a wide range of toxic effects on aquatic biota, including fish [
11]. Fish play a prominent role in the functioning and balance of aquatic ecosystems [
12] and are particularly sensitive to environmental pollutants [
13,
14,
15], thus representing good indicators of environmental water quality.
Fish assimilate Pb by direct ingestion (i.e., food and water), ion exchange across lipophilic membranes (e.g., the gills), or adsorption across specific tissue and membrane surfaces [
3]. In fish, numerous detrimental effects induced by Pb have been reported, among which are genotoxicity [
16], oxidative stress induction [
17,
18,
19,
20] change in the activities of immune-related enzymes and genes [
8,
21,
22] and histological changes of some tissues and organs [
16,
23,
24,
25].
Given its role in metal detoxification and excretion, the liver is a target organ for Pb [
21,
26]. In freshwater species, waterborne Pb mainly enters through the gills reaching the liver via the circulatory system [
26]. Therefore, the liver has been claimed as the main target for Pb accumulation in fish [
23,
27,
28], and the noxious effect induced in this complex organ have been investigated in several species. A systematic review of the literature provides plenty of evidence of adverse effects induced in the liver, including a modification of liver enzyme activity [
29], the occurrence of metabolic disorders [
30], alteration of liver morphology [
31,
32,
33,
34]. However, all these studies refer to high Pb concentrations, and there is a lack of knowledge on the effects of Pb at environmentally relevant concentrations [
35,
36].
It must be emphasized that fish exposed to very low concentrations of Pb may not show obvious signs of pathology, but the subtle morphological and functional alterations induced by the metal can reduce the health of individuals with important and dramatic repercussions at the population level.
Since the 1980s, zebrafish (
Danio rerio) has been used in a broad spectrum of research fields due to its small body size, short reproductive cycle, easy husbandry, and high homology to the human genome. All these studies provided a powerful basis for using zebrafish as a model organism for aquatic ecotoxicology [
37,
38]. It is surprising that a very limited number of researches investigated the hepatotoxic effects of Pb in zebrafish liver, and so far, only three studies are available focusing on the induction of metabolic disorder and oxidative stress [
39,
40] and morphological alterations following chronic exposure to high Pb concentrations (60 mg/L) [
41].
To fill the knowledge gap on Pb hepatotoxicity in fish, here we evaluated, for the first time, the effect induced in zebrafish liver by two very low and environmentally relevant concentrations of Pb (2.5 and 5 μg/L) after 48, 96, and 192 hours of exposure. The tested doses have been selected based on Pb concentrations found in aquatic environments worldwide and, particularly, in the concentrations range of Pb in surface waters; therefore, our results would also support the implementation of risk assessment protocols.
Given the general paucity of information about morphological and functional injuries induced by Pb in the fish liver, we first assessed the histological alterations, which are widely recognised as the best tool for assessing the effects of chemical contaminants, including heavy metals [
42]. Moreover, to allow a more reliable and objective comparison among experimental groups, we applied a semi-quantitative method for evaluating the severity of histological changes.
Since histopathological lesions represent an integration of the effects of prior biochemical and physiological perturbations [
43], in a second step, we analyzed the modulations of some genes involved in i) metal detoxification (metallothionein,
mtf1) and ii) oxidative stress defense (i.e., superoxide dismutase-
sod1 and catalase-
cat) to better clarify the molecular mechanisms underlying Pb hepatotoxicity.
Metallothionein (MTs) are low-molecular-weight proteins responsible for metal binding, which are effective in non-essential metal detoxification and protection from oxidative stress [
44]. The biosynthesis of MTs in fish is induced by a variety of metals, including Pb, providing an excellent biomarker of exposure to metals [
45,
46].
Pb toxicity in the liver could be mediated through different mechanisms, but the most common response in both fish and mammalian models is the imbalance between reactive oxygen species (ROS) production and the removal of such molecules [
47]. The cell counteracts ROS overproduction by induction of antioxidant molecules, which may be either enzymatic (e.g., catalase and superoxide dismutase) or non-enzymatic (e.g., glutathione). In fish exposed to Pb, the antioxidant responses have been demonstrated in several species, thus supporting the role of such molecules as biomarkers of oxidative stress induced by heavy metals [
21,
48,
49,
50].
To the best of our knowledge, this is the first study documenting the morphological, morphometric, and functional alterations of low Pb concentrations in fish liver. Our results, providing new insights into lead-induced hepatotoxicity, will also contribute to a better understanding of the risk posed by heavy metals to wildlife species under a realistic exposure scenario.
2. Materials and Methods
2.1. Fish maintenance
A total of 70 individuals of both sexes (length, 3.5±0.5 cm, and weight, 0.43±0.06 g) were obtained from a local fish retailer. For two weeks, fish were acclimatized under controlled conditions in aquaria filled with dechlorinated tap water (temperature=26±0.5°C, pH=7.3, conductivity=300 μs/cm, dissolved oxygen=8±1 mg/L, hardness=100 mg/L CaCO3, 14:10 light regime). During the acclimatization period, half of the water was renewed daily, and fish were fed daily with commercial fish food.
2.2. Exposure Conditions
To obtain the two nominal concentrations of 2.5 µg/L (low concentration) and 5 µg/L (high concentration), a stock solution (1000 µg/L) of lead acetate was prepared in distilled water [Pb(CH3CO2)2 3H2O, Sigma-Aldrich Chemical Co., St. Louis, MO, USA]; then an adequate amount was diluted in dechlorinated water.
To determine Pb in water samples, an Elan DRC-e inductively Coupled Plasma-Mass Spectrometry (ICP-MS) (PerkinElmer SCIEX, Woodbridge, ON, Canada) was used. Samples were diluted in ultrapure nitric acid (500 µL) and then introduced into the instrument system using a PerkinElmer AS-93 plus autosampler and a cross-flow nebulizer with a Scott-type spray chamber. Quantitative analysis was performed by constructing the calibration curve for the lead on five different Plasma-Mass (calibration range of 0.1-50 g/L). The analytical verification of the actual concentrations was performed throughout the experiment (starting from time 0 every 24 hours) (
Table S1). No evident variation was recorded in agreement with previous literature data [25 and references therein].
The two selected concentrations corresponded to 0.00146% and 0.00292% of the median lethal concentration at 96 hours (LC50
96h), respectively for adult zebrafish [
39]. Moreover, both doses were selected considering pre-existing data on the worldwide concentration of Pb in surface water and can be considered very low and environmentally realistic [
3,
35].
Fish were exposed to the low or the high Pb dose for 48, 96, and 192 hours, resulting in 6 experimental groups (n=10). The control group (n=10) was maintained in aquaria filled with dechlorinated tap water.
During the experiment, temperature, pH, conductivity, dissolved oxygen, hardness, and photoperiod were monitored daily and kept constant, as described for the acclimatization period. Fish were fed on alternate days, and food waste and debris were removed daily using a fine mesh.
After 48, 96, and 192 hours, fish were immersed in an anesthetic water bath containing ethyl 3-aminobenzoate methanesulfonate (20 mg/L MS 222, Sandoz, Sigma-Aldrich, St. Louis, MO, USA), and the liver was rapidly dissected and processed for subsequent analyses as reported below. For each experimental unit, including the control, two replicates were conducted. The use of animals in this study was approved by the Institutional Animal Care and Use Committee at the National University of Entre Rios and the Italian University Institute of Rosario (Rosario, Argentina; protocol N°028/12).
2.3. Histology and Histopathological assessment
Excised liver samples (n=4) were immediately fixed in 4% glutaraldehyde (Electron Microscopy Sciences, Hatfield, PA, USA) in phosphate-buffered saline solution (PBS 0.1 M, pH 7.2, 4°C) and post-fixed in osmium tetroxide (1% in PBS) for 2 hours. Samples were dehydrated through a graded ethanol series, placed in propylene oxide, and embedded in Epon-Araldite (Araldite 502/Embed 812, Electron Microscopy Sciences). Longitudinal serial semithin sections of 1 μm, obtained using a Leica UltraCut UCT (Leica Microsystems, Wetzlar, Germany), were mounted on glass slides, stained with toluidine blue, and observed under an LM Leitz Dialux 20 E.B. (Leica Microsystems, Wetzlar, Germany) equipped with a digital camera.
The prevalence of each histological alteration was obtained by doing the ratio between the number of fish affected by a specific alteration and the total number of fish. We also determined the histological changes' severity using a semi-quantitative method according to previous literature data [
51,
52]. Briefly, the alterations were attributed to a specific reaction pattern (circulatory disturbances, regressive changes, progressive changes, and inflammation). Then an importance factor was assigned to each observed alteration following the relevance of the change and its pathological importance (from 1, minimal pathological importance, to 3, marked pathological importance). A score value was then assigned based on the degree and extent of each lesion, as follows: 0 (unchanged), 2 (mild occurrence), 4 (moderate occurrence), and 6 (severe occurrence). (
Table S2). The organ index (I
org) representing the degree of organ damage was calculated using the importance factor and the score value according to the following formula:
where: org=organ, rp=reaction pattern, alt=alteration; a=score value; w=importance factor.
2.4. Lipid droplets content
Lipid droplets analysis has been performed on semithin sections (toluidine blue-stained). Four liver sections were photographed (100×) for each animal of both the control and Pb exposed groups (n=4) and evaluated for the percentage of area occupied by lipid droplets. The lipid granules were isolated in each micrograph using the free and open-source ImageJ software (NIH, developed at the National Institutes of Health, a part of the U.S. Department of Health and Human Services), and the total area occupied by the granules was quantified.
The results, expressed as the percentage of area occupied by the lipid granules in each section, were statistically compared using the two-way ANOVA followed by Tukey's multiple comparisons tests (at a significance level of 0.05). Data were checked for normality (Shapiro-Wilk test) and presented as mean±standard deviation.
2.5. Quantitative Real-Time PCR
Excised liver samples of animals of both treated and control groups (n=6) were promptly stored at -80°C for subsequent real-time PCR analyses. Total RNA was extracted using the PureLink RNA Mini Kit and the PureLink™ DNase Set (Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer's protocol. The quantity and quality of RNA were verified using a NanoDrop One spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and 1.5% agarose gel electrophoresis, respectively. 2 μg of total RNA was utilized for first-strand cDNA synthesis using the high capacity RNA to cDNA kit (Applied Biosystems, Foster City, CA, USA); the resulting cDNA was kept at -20°C. cDNA was used as a template for quantitative reverse transcription-polymerase chain reaction (RT-qPCR) analysis to quantify the expression of the metal-regulatory transcription factor 1 (mtf1, NCBI Reference Sequence NM_152981.1), the superoxide dismutase 1 (sod1, NCBI Reference Sequence NM_131294.1) and the catalase (cat, NCBI Reference Sequence NM_130912.2). RT-qPCR was performed in triplicate in a Light Cycler (Applied Biosystems StepOne, Real-Time PCR System, Foster City, CA, USA) using the TaqMan Gene Expression Assays (Thermo Fisher Scientific, Waltham, MA, USA). Each reaction contained 2 μL of cDNA, 10 μL of master mix (TaqMan Universal Master Mix II, Applied Biosystems), 1 μL of assay mix (TaqMan Gene Expression Assay), and 7 μL of RNase- and DNase-free water and was run according to the manufacturer's instructions: one cycle at 50°C for 2 minutes, 95°C for 10 minutes, 40 cycles at 95°C for 15 seconds, and 60°C for 1 minute.
The glyceraldehyde-3-phosphate dehydrogenase (gapdh, NCBI Reference Sequence: NM_001115114.1) and the actin beta 1 (actb1, NCBI Reference Sequence: NM_131031.2) were used as internal reference genes. The relative copy number of each analyzed gene was calculated according to the 2−ΔΔCT comparative CT method.
Author Contributions
Methodology, R.M., V.C., E.P., L.O.; investigation, R.M., V.C., E.B.; formal analysis, A.I.M.A., F.T., L.O., M.M.; data curation, F.T., S.S.; software, A.I.M.A., F.T. S.S.; writing—original draft preparation, R.M., V.C., E.B.; writing—review and editing, E.B.; project administration, E.B. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Light micrographs of Danio rerio liver under basal condition. (a) General organization of the liver parenchyma; note the bile ducts enclosed by cuboidal epithelium. (b) High magnification showing the space of Disse between the hepatocytes and the sinusoid wall. Note erythrocytes and a few macrophages in the lumen of the veins. bd=bile duct, s=sinusoid, v=vein, n=nucleus, arrow=glycogen granules, arrowhead=lipid droplets, sd=space of Disse.
Figure 1.
Light micrographs of Danio rerio liver under basal condition. (a) General organization of the liver parenchyma; note the bile ducts enclosed by cuboidal epithelium. (b) High magnification showing the space of Disse between the hepatocytes and the sinusoid wall. Note erythrocytes and a few macrophages in the lumen of the veins. bd=bile duct, s=sinusoid, v=vein, n=nucleus, arrow=glycogen granules, arrowhead=lipid droplets, sd=space of Disse.
Figure 2.
Light micrographs of Danio rerio liver after exposure to 2.5 µg/L of Pb. (a-c) After 48 hours of exposure, bile duct degeneration and the congestion of blood vessels and sinusoids were observed. Note the increase in lipid droplet content. (d-g) After 96 hours of exposure, the cuboidal epithelium lining the bile ducts was modified, and the detachment of the duct epithelium was evident. Note the congestion of blood vessels and sinusoids, macrophage proliferation, cytoplasmic vacuolization, the increase in lipid droplet content, and the appearance of both apoptotic and necrotic hepatocytes. Also, lysed areas were frequently observed. (h-j) After 192 hours of exposure, cytoplasm vacuolization, the congestion of vessels and sinusoids, and numerous lipid droplets were frequently detected. Note the detachment of bile duct epithelium, macrophage proliferation, and both apoptotic and necrotic hepatocytes. bdd=bile duct degeneration, vc=vessel congestion, sc=sinusoids congestion, black arrow= bile duct epithelial detachment, black arrowhead=lipid droplets, white arrowhead=cytoplasmic vacuolization, white arrow=macrophages proliferation, white star=apoptotic cell, black star=necrotic cell, la=lysed area.
Figure 2.
Light micrographs of Danio rerio liver after exposure to 2.5 µg/L of Pb. (a-c) After 48 hours of exposure, bile duct degeneration and the congestion of blood vessels and sinusoids were observed. Note the increase in lipid droplet content. (d-g) After 96 hours of exposure, the cuboidal epithelium lining the bile ducts was modified, and the detachment of the duct epithelium was evident. Note the congestion of blood vessels and sinusoids, macrophage proliferation, cytoplasmic vacuolization, the increase in lipid droplet content, and the appearance of both apoptotic and necrotic hepatocytes. Also, lysed areas were frequently observed. (h-j) After 192 hours of exposure, cytoplasm vacuolization, the congestion of vessels and sinusoids, and numerous lipid droplets were frequently detected. Note the detachment of bile duct epithelium, macrophage proliferation, and both apoptotic and necrotic hepatocytes. bdd=bile duct degeneration, vc=vessel congestion, sc=sinusoids congestion, black arrow= bile duct epithelial detachment, black arrowhead=lipid droplets, white arrowhead=cytoplasmic vacuolization, white arrow=macrophages proliferation, white star=apoptotic cell, black star=necrotic cell, la=lysed area.
Figure 3.
Prevalence of histological alterations in Danio rerio liver after exposure to 2.5 and 5 µg/L of Pb for 48, 96, and 192 hours. a significance of low concentration group with respect to the control group, b significance of the high concentration group with respect to the control group, c significance of the high concentration group with respect to the low concentration group. *p≤0.05; **p≤0.005; ****p≤0.0001.
Figure 3.
Prevalence of histological alterations in Danio rerio liver after exposure to 2.5 and 5 µg/L of Pb for 48, 96, and 192 hours. a significance of low concentration group with respect to the control group, b significance of the high concentration group with respect to the control group, c significance of the high concentration group with respect to the low concentration group. *p≤0.05; **p≤0.005; ****p≤0.0001.
Figure 4.
Light micrographs of Danio rerio liver after exposure to 5 µg/L of Pb. (a-c) After 48 hours of exposure, bile duct degeneration, bile duct epithelial detachment, wide lysed areas and numerous lipid droplets could be detected. Moreover, the vacuolization of hepatocytes cytoplasm, the congestion of blood vessels, and both apoptotic and necrotic were frequently visible. Note the proliferation of macrophages. (d-f) After 96 hours of exposure, numerous degenerations such as vessels and sinusoid congestion, macrophage proliferation, and apoptotic hepatocytes, were visible in all samples. Also, the increase in lipid contents, cytoplasmic vacuolization, bile duct degeneration, detachment of the bile duct epithelium and lysed areas were frequently detected. (g-i) After 192 hours of exposure, all the considered alterations were detected in all samples. bdd=bile duct degeneration, vc=vessel congestion, sc=sinusoids congestion, black arrow= bile duct epithelial detachment, black arrowhead=lipid droplets, white arrowhead=cytoplasmic vacuolization, white arrow=macrophages proliferation, white star=apoptotic cell, black star=necrotic cell, la=lysed area.
Figure 4.
Light micrographs of Danio rerio liver after exposure to 5 µg/L of Pb. (a-c) After 48 hours of exposure, bile duct degeneration, bile duct epithelial detachment, wide lysed areas and numerous lipid droplets could be detected. Moreover, the vacuolization of hepatocytes cytoplasm, the congestion of blood vessels, and both apoptotic and necrotic were frequently visible. Note the proliferation of macrophages. (d-f) After 96 hours of exposure, numerous degenerations such as vessels and sinusoid congestion, macrophage proliferation, and apoptotic hepatocytes, were visible in all samples. Also, the increase in lipid contents, cytoplasmic vacuolization, bile duct degeneration, detachment of the bile duct epithelium and lysed areas were frequently detected. (g-i) After 192 hours of exposure, all the considered alterations were detected in all samples. bdd=bile duct degeneration, vc=vessel congestion, sc=sinusoids congestion, black arrow= bile duct epithelial detachment, black arrowhead=lipid droplets, white arrowhead=cytoplasmic vacuolization, white arrow=macrophages proliferation, white star=apoptotic cell, black star=necrotic cell, la=lysed area.
Figure 5.
Percentage of the area occupied by lipid droplets in Danio rerio liver after exposure to 2.5 and 5 µg/L of Pb for 48, 96, and 192 hours. Graphs indicate the mean ± S.D. Asterisks indicate significant differences between the treated and control groups. Hashtags indicate significant differences between the high Pb concentration group and the low Pb concentration group. ***p≤0.001; ****p≤0.0001; ##p≤0.005; ###p≤0.001.
Figure 5.
Percentage of the area occupied by lipid droplets in Danio rerio liver after exposure to 2.5 and 5 µg/L of Pb for 48, 96, and 192 hours. Graphs indicate the mean ± S.D. Asterisks indicate significant differences between the treated and control groups. Hashtags indicate significant differences between the high Pb concentration group and the low Pb concentration group. ***p≤0.001; ****p≤0.0001; ##p≤0.005; ###p≤0.001.
Figure 6.
Gene expression in Danio rerio liver after exposure to 2.5 and 5 µg/L of Pb for 48, 96, and 192 hours. Graphs indicate the mean±S.D. (a) Metallothioneins (mtf1), (b) Catalase (cat), (c) Superoxide dismutase (sod1). Asterisks indicate significant differences between treated and control groups. Hashtags indicate significant differences between the high concentration and the low concentration groups. *p≤0.05; **p≤0.005; ***p≤0.001; ****p≤0.0001; ####p≤0.0001.
Figure 6.
Gene expression in Danio rerio liver after exposure to 2.5 and 5 µg/L of Pb for 48, 96, and 192 hours. Graphs indicate the mean±S.D. (a) Metallothioneins (mtf1), (b) Catalase (cat), (c) Superoxide dismutase (sod1). Asterisks indicate significant differences between treated and control groups. Hashtags indicate significant differences between the high concentration and the low concentration groups. *p≤0.05; **p≤0.005; ***p≤0.001; ****p≤0.0001; ####p≤0.0001.
Table 1.
Comparison in the organ index (mean ± SD) between control and Pb exposed groups.
Table 1.
Comparison in the organ index (mean ± SD) between control and Pb exposed groups.
|
2.5 µg/L Pb |
5 µg/L Pb |
CTRL |
0.00±0.00 |
0.00±0.00 |
48 hours |
2.66±1.15 |
16.00±5.29(a***)(b**)
|
96 hours |
14.66±3.05(a**)(c*)
|
42.00±2.40(a****)(b****)(c****)
|
192 hours |
35.33±4.61(a****)(c****)
|
53.33±3.05(a****)(b****)(c*)
|